Geometrical dependence on the stress and slip tendency acting on the subduction megathrust of the Nankai seismogenic zone off Kumano
Masataka Kinoshita ORCID: orcid.org/0000-0003-1144-36511,
Kazuya Shiraishi2,
Evi Demetriou3,
Yoshitaka Hashimoto4 &
Weiren Lin5
We assess the dependence of megathrust geometry near the updip edge of the Nankai seismogenic zone according to the slip tendency of the megathrust based on a reprocessed 3D PSDM seismic volume image and International Ocean Discovery Program (IODP) NanTroSEIZE drilling data off the coast of the Kii Peninsula, central Japan. The plate boundary fault surface is manually picked from the 3D seismic volume and is divided into three groups: low dip (10°–15°) trending NW to N in the NE region, intermediate dip (~ 25°) trending north in the western region, and high dip (30°–40°) trending NW in the SE corner. We then calculate the overburden (Sv) by converting the modeled 3D velocity to bulk density. Here, Sv ranges from 100 MPa near the SW edge to 160 MPa on the NE corner. In order to derive horizontal principal stresses (SH and Sh) and the pore fluid pressure (Pp), we assign the ratio of horizontal-to-vertical principal stresses, r (= SH/Sv), and the ratio of pore fluid pressure-to-vertical stress, λ (= Pp/Sv) based on IODP drilling data. The directions of SH and the slip on the fault are set parallel to the plate convergence vector (N55W). Assuming a triaxial condition (SH > Sh = Sv), the slip tendency (Ts) can be calculated from the dip angle and dip azimuth of the fault surface, r, and λ. We found that Ts is sensitive to the variation in fault geometry. For r = 1.2 and λ = 0.7 (normalized pore pressure ratio λ* = 0.25), Ts is low (~ 0.1) in the low-angle dip region and higher (> 0.2) in the high-angle dip region. This suggests that the high-angle fault is optimally oriented under this condition. Low Ts in the low-angle dip region would correspond to a weaker region due to a low frictional strength, assuming that the fault surface should slip simultaneously. Using the high pore pressure ratio (λ* ~ 0.85) beneath the fault zone, the fault beneath IODP Site C0002 is very close to failure if r > 1.2.
According to the Coulomb failure theory, an earthquake or a slip along a fault can occur if the ratio of shear stress to effective normal stress exceeds the friction coefficient on the fault surface. This ratio is called the slip tendency (Ts). In order to estimate Ts along a pre-existing fault, we need to know the fault geometry (dip angle and dip azimuth), because the shear and normal stresses vary with dip.
Subduction zone earthquakes occur along a plate boundary megathrust and have geometry that varies either within a main fault trace or by fault branching. Such fluctuation can lead to variation in Ts within a megathrust system.
In the present paper, we consider the case for the Nankai Trough seismogenic zone off Kumano, central Japan, which is one of the most well-studied subduction zones. From the reprocessed 3D PSDM volume, we obtained the fine-scale geometry (dip angle and azimuth) of the updip region of plate boundary fault. Integrating this fault geometry with deep drilling data obtained within the 3D survey region, we calculate the slip tendency of the fault and assess the degree to which the fault geometry affects the slip tendency value. Then, we discuss the implication of geometrical variation, e.g., whether shallow regions of the megathrust exhibit different behaviors at different times or locations, or whether pore fluid pressures and friction coefficients vary in different regions of the fault surface.
Tectonic setting and previous research
The Nankai Trough is a plate convergent boundary where the Philippine Sea Plate subducts beneath western Honshu, Japan. Geodetic observations (Miyazaki and Heki 2001) show that the plate convergence occurs at a rate of 63–68 mm/year in the N55W direction. DeMets et al. (2010) reported a similar value of 58.4 ± 1.2 mm/year toward WNW.
In 1944, a 100-km-wide locked region of a subducting plate interface sheared, causing the devastating M8 Tonankai earthquake. The slip area of this event was estimated in several studies based on seismic inversion (Kikuchi et al. 2003; Ichinose et al. 2003) and tsunami inversion (Baba et al. 2006). The area with this coseismic slip (contours in Fig. 1) coincides with the present strong interseismic coupling, detected as a large slip deficit area inferred from seafloor geodetic measurements (Yokota et al. 2016). This indicates that the coupling/locking condition is related to seismogenic conditions (i.e., rupture nucleation, propagation, and arrest).
Index map showing the reprocessed area of the 3D seismic volume (yellow rectangle) and IODP NanTroSEIZE drill sites (white circles). Contours indicate the fault slip (contour increment: 0.5 m) of the 1944 Tonankai earthquake (modified from Kikuchi et al., 2003)
Intensive seismic monitoring using broadband seismometers or pressure gauges in the forearc region of the Nankai Trough off Kumano documented a wide spectrum of seismic activities, such as very-low-frequency events (VLFEs), low-frequency tremor (LFT) swarms, or slow slip events. Using the onshore broadband network (Hi-net), Obara and Ito (2005) detected the occurrence of VLFEs in the Nankai forearc region (including those off Kumano in 2004) with their focus located along the décollement or splay faults and thrust-type mechanisms. Obana and Kodaira (2009) observed low-frequency tremors near the shallowest region of a splay fault in the outer wedge of the accretionary prism off Kumano. Using four broad-band ocean bottom seismometers, Sugioka et al. (2012) detected VLFEs at the shallowest region of the plate boundary of the Nankai Trough in 2009, generated by slips along very low-angle thrusts. On April 1, 2016, the Mw 6.0 Mie-ken Nanto-Oki event occurred in the middle of the locked zone (e.g., Nakano et al. 2018). Although the focal depths differ among researchers, the main shock is supposed to have occurred on the plate interface with its mechanism being low-angle thrust-type dipping at either 14° or 76° (Wallace et al. 2016).
In order to understand the state and properties of the Nankai Trough seismogenic zone, a series of ocean scientific drillings has been conducted since 2007 as part of the Integrated Ocean Drilling Program (IODP) Nankai Trough Seismogenic Zone Experiments (NanTroSEIZE) (Fig. 1; Tobin and Kinoshita 2006; Kinoshita et al. 2009; Expedition 326 Scientists 2011; Henry et al. 2012; Kopf et al. 2011; Saffer et al. 2010; Saito et al. 2010; Screaton et al. 2009; Strasser et al. 2014; Tobin et al. 2009; Tobin et al. 2014; Kinoshita et al. 2014; Kopf et al. 2016). The ultimate goal of NanTroSEIZE is to core, log, and monitor the locked region of the plate boundary megathrust approximately 5 km below the seafloor at Site C0002.
One of the important components of NanTroSEIZE is to install borehole observatories for long-term monitoring of in situ pressure, temperature, strain, deformation, and seismicity. In 2010, the first permanent observatory was installed at approximately 900 mbsf at Site C0002 above the updip limit of the locked zone during IODP Expedition 332 (Kopf et al. 2011). In 2016, the second observatory was installed at the shallow region of the megasplay fault zone at Site C0010 (~ 400 mbsf) during IODP Expedition 365 (Kopf et al. 2011). Using data from both observatories, Araki et al. (2017) reported recurring slow-slip events (SSEs) on the plate interface immediately seaward of 1944 M8 earthquake rupture area, with their duration ranging days to weeks. They also reported swarms of low-frequency tremors following these SSEs.
A wide spectrum of seismicity suggests that the plate interface around the updip edge of the locked zone is sensitive to external stress because of the weakness of faults (Obara and Kato 2016). This is likely supported by the low-velocity zone extending along the décollement (Kamei et al. 2012; Park et al. 2010; Bangs et al. 2009).
Drilled core samples and in situ data (logging and downhole experiments) provide direct lines of information of the megathrust fault, but this has limited value unless this point-source information is properly extrapolated around the locked zone. Precise structural imaging at the plate boundary fault system and in accreted sediments is essential.
In 2006, 3D seismic reflection data acquisition and processing were carried out (Moore et al. 2007). However, the obtained volume could not necessarily resolve deep structures due to relatively short (4.5 km) streamers and the strong Kuroshio current. In 2016, this seismic volume was reprocessed using state-of-the-art advanced technology (e.g., better attenuation of noise and multiple reflections, broadband processing with deghosting, and 4D regularization to mitigate non-uniform fold distribution caused by the Kuroshio current) in order to obtain clearer images and a more reliable P-wave velocity (Vp) model (Shiraishi et al. 2016). Improved images in the shallow accretionary wedge revealed deformation features (e.g., branching of splay faults and thrusting of the lower Shikoku Basin formation). However, due to short streamer length, the spatial resolution in deeper regions was not improved significantly after reprocessing in 2016. This gives rise to an inherent uncertainty in the geometry of the plate boundary fault. Despite such limitations, we used the reprocessed 3D data, which are the best data available thus far around the drill site.
Figure 2 shows an outline of the method used in the present study and is explained in the following.
Outline of the method used in the present study. See the "Methods/Experimental" section for details
Slip tendency and slip likelihood
Slip tendency was first proposed by Morris et al. (1996) for the assessment of potential fault activity among multiple pre-existing faults under a uniform stress state:
$$ Ts=\tau /\left({\sigma}_n-{P}_p\right), $$
where τ is the shear stress, σn is the normal stress acting on a plane, and Pp is the pore fluid pressure (see Table 1 for definitions of parameters).
Table 1 List of parameters
Tong and Yin (2011) revised this model to include the cohesive shear strength and defined a new index, the reactivation tendency, as τ/(τ0 + μσn_e), where τ0 is the cohesion strength, μ is the intrinsic frictional coefficient, and σn_e is the effective normal stress.
Although this definition describes the "likelihood" of fault slip, in the present study, we use the definition of Morris et al. (1996) in Eq. (1) because, as of yet, we have no measurements on the intrinsic frictional coefficient on the fault, and because the cohesive strength is negligibly small in a reverse faulting regime in subduction zone megathrusts. A similar approach was taken by Miyakawa and Otsubo (2015) for the assessment of fault activity in NE Japan.
Alternatively, we use the slip likelihood (SL), which is defined as
$$ \mathrm{SL}\equiv \frac{1}{\mu}\bullet \frac{\tau }{\sigma_n-{P}_p}=\frac{1}{\mu}\bullet {T}_s. $$
According the Coulomb failure criterion, the fault will slip if SL = 1.
Templeton and Rice (2008) proposed an alternative definition of slip likelihood, closeness to failure (CF), which is defined as the distance between the Mohr circle and failure criteria. Since SL and CF are essentially identical, we herein use Ts and SL, mainly for reasons of simplicity.
Picking a fault reflector from the 3D seismic volume
The plate boundary fault was extracted from the reprocessed 3D seismic volume using Petrel® software in the area shown in Figs. 1 and 3 shows an example of a fault surface, which is manually picked by eye along the negative-polarity peaks (shown as blue) beneath the positive and generally strongest peaks. Reflection images for the reprocessed volume are clearer than the original images but still suffer from insufficient spatial resolution and migration noise because of the great depth of the volume. As shown in Fig. 3, the wavelength of reflected waves from the plate boundary fault is approximately 200 m. Since the seismic resolution is usually considered to be approximately 1/4 the dominant wavelength, the spatial resolution is approximately 50 m. We found some discontinuities in the reflectors, in which case, we attempted to follow as smooth a trend as possible. The picked surface is then smoothed in Petrel, but artificial fluctuations remain due to picking error or systematic offset among inlines. Overall, we assess that the relative uncertainty in the depth of the fault can be as large as ± 50 m.
Example of picked reflection from the fault surface. (Left) Shallow section (seafloor to the bottom simulating reflector). (Right) Deep section (plate boundary fault and top of oceanic crust). By comparing with those from seafloor (positive), BSR (negative), and top of oceanic crust (positive), the picked fault has a negative polarity
Vertical stress on the fault
The vertical stress applied to the fault surface is estimated from the 3D P-wave velocity model following the procedure in Fig. 2.
Erickson and Jarrard (1998) developed a global empirical relationship between P-wave velocity and porosity for consolidated and non-consolidated siliciclastic sediments, covering a velocity range of up to approximately 5 km/s and a full range of porosity (0–1). Kitajima et al. (2017) used this relation to determine the in situ porosity from Vp logs at IODP Site C0002. Above the picked fault surface, the Vp from the reprocessed 3D model ranges from 1.5 to 5.5 km/s (Shiraishi et al. 2016), and the estimated bulk density ranges from 1.5 to 2.7 g/cm3. In order to convert porosity to bulk density, the grain density is estimated for core samples at Site C0002. The grain density increases slightly from 2.68 g/cm3 at 950.5 mbsf to 2.72 g/cm3 at 3058.5 mbsf (Tobin et al. 2015). The vertical stress (Sv) is calculated by integrating the bulk density from the sea surface to the fault surface.
The effective vertical stress is calculated by subtracting the pore fluid pressure (Pp) from Sv. Tsuji et al. (2014) estimated the pore fluid pressure from the velocity anomaly obtained through full-wave inversion of the wide-angle seismic refraction data in this region. According to their estimation, the normalized pore pressure ratio λ*, defined as
$$ {\lambda}^{\ast}\kern0.5em =\kern0.5em \frac{P_p-{P}_{hyd}}{S_v-{P}_{hyd}}\kern0.75em {P}_{hyd}:\mathrm{hydrostatic}\ \mathrm{pressure}. $$
is 0 to 0.5 in the hanging wall accretionary prism, with a higher value (~ 0.7) around Site C0002. On the other hand, Kitajima et al. (2017) estimated λ* to be less than 0.2 (Fig. 4). In the present study, we follow the estimates of Kitajima et al., where λ* = 0.25, or λ = Pp/Sv ~ 0.7. This likely simulates the state on the hanging wall side of the fault. Since the fault zone itself and the footwall side can have a higher λ* value (~ 0.85, as predicted by Tsuji et al. 2014), we consider these two cases for Ts estimation.
(Left) Principal stress and pore fluid pressure estimated by Kitajima et al. (2017). (Right) Estimated pore pressure ratio, normalized pore pressure ratio, and stress ratio. SH: maximum (compressional) horizontal stress, Sh: minimum horizontal stress, Sv: vertical stress, Pp: pore fluid pressure, and Phyd: hydrostatic pore pressure. λ = Pp/Sv, λ* = (Pp-Phyd)/(Sv-Phyd). The dotted curves are the stress and pore fluid pressure models used to calculate the slip tendency, Ts
The effective vertical stress σv is shown as
$$ {\sigma}_v={S}_v-{P}_p=\left(1-\lambda \right){S}_v=\left(1-{\lambda}^{\ast}\right)\left(1-{\lambda}_h\right){S}_v=0.4\left(1-{\lambda}^{\ast}\right){S}_v, $$
where λh = Phyd/Sv ~ 0.6 (Fig. 4).
For convenience, the relation between λ and λ* in the present study is
$$ {\lambda}^{\ast }=\left(\lambda -{\lambda}_h\right)/\left(1-{\lambda}_h\right)\sim \left(\lambda -0.6\right)/0.4. $$
Horizontal principal stresses
Since we have the 3D geometry of the plate boundary fault, the stress acting on the fault must also be considered to be 3D. Assuming the vertical stress (Sv), the maximum and minimum horizontal stresses (SH and Sh) are the principal stresses. In the following, SH and Sh are given for a thrust-faulting regime.
Magnitude of S H
The magnitude of horizontal stresses in the crust is a long-standing issue in geomechanics. For a non-tectonic sedimentary basin, the state of stress is lithostatic (McGarr 1988), where all three principal stresses are identical (e.g., SH = Sh = Sv). Zoback (2007) indicated that the crust is in a state of incipient frictional failure. In such a case, we can predict the principal stress field if the faulting regime (normal, strike-slip, or thrust) and an appropriate value of the frictional coefficient are given.
Brown and Hoek (1978) reviewed and compiled the results of in situ stress measurements and concluded that the ratio of average horizontal-to-vertical stressed varies in a wide range (0.5 to > 3.5), although this ratio generally decreases with increasing depth. Rummel et al. (1986) summarized that the ratio of maximum compressional stress to the minimum compression stress was slightly above 1 below a depth of 1 km in the continental crust.
Based on geomechanical experiments on core samples drilled in this area, Kitajima et al. (2017) estimated the horizontal maximum principal stress (SH) and pore fluid pressure (Pp) along Site C0002 down to 3000 mbsf. Figure 4 indicates that the ratio of maximum-to-minimum compression stresses, or SH/Sv (= r) in this case, is 1.0 to 1.4 at Site C0002. Here, we assume r = 1.2 on the plate boundary fault. For effective stress σH, we make a similar assumption regarding Sv:
$$ {\sigma}_H={S}_H-{P}_p=\left(r-\lambda \right){S}_v=\left(r-0.4{\lambda}^{\ast }-0.6\right){S}_v. $$
Azimuth of S H and slip
We simply assume that SH is oriented parallel to the plate convergence vector. The plate convergence vector between the AM and PH plates is estimated as N55W at a rate of 63 to 68 mm/year (Miyazaki and Heki 2001) based on geodetic data. DeMets et al. (2010) reported a rate of 58.4 ± 1.2 mm/year in the WNW direction. In the present study, we adopt N55W as the azimuth of SH.
We also need to determine the slip direction in order to estimate the shear stress. A simple but reasonable assumption is that the slip direction is parallel to the maximum horizontal stress (SH), which is parallel to the plate convergence (N55W; Fig. 5).
Schematic diagram showing the geometry and direction of stress and slip. The plate slip direction (assumed to be equal to the SH direction) is set parallel to the plate conversion vector (N55W; Miyazaki and Heki 2001)
Horizontal minimum stress S h
Ito et al. (2009) estimated the stress ratio from the stress tensor inversion analysis of VLFEs and concluded that the intermediate stress (Sh for this case) is larger than the minimum principal stress (Sv) in the off-Kumano region.
On the other hand, the minimum horizontal stress (Sh) was estimated through downhole measurements (Leak-off test; Strasser et al. 2014; Tobin et al. 2015), which indicates that Sh is slightly (several MPa) smaller than Sv at depths of 2000 mbsf, suggesting a normal faulting or strike-slip regime.
Although the stress regime is different across the plate boundary fault, horizontal stresses should be larger than the vertical stress for a reverse fault. Since there are no constraints on the principal stresses near the plate interface at present, we assume a triaxial state, in which the minimum horizontal stress (Sh) is equal to Sv, which yields
$$ {\sigma}_H>{\sigma}_h={\sigma}_v, $$
where σH and σh are the effective maximum and minimum horizontal stresses, respectively.
Normal/shear stress and slip tendency
We now have all the necessary information to calculate the effective normal stress and shear stress: the effective principal stresses (magnitude, direction, and pore fluid pressure), geometry (dip angle and dip azimuth) of the plate boundary fault, and the potential slip direction (N55W). The effective normal stress and shear stress acting at each location on the plate boundary fault are
$$ {\sigma}_{n\_e}=\left[\ \left({\left(\sin \alpha \sin \delta \right)}^2+{\left(\cos \delta \right)}^2\right)\left(1-\lambda \right)+{\left(\cos \alpha \sin \delta \right)}^2\left(r-\lambda \right)\ \right]\ {S}_v. $$
$$ \tau =A\sin \kern0ex \delta \cos \kern0ex \delta \kern0.5em \left(r-1\right){S}_v. $$
where \( A={\left(\frac{\cos^2\delta }{\cos^2\alpha }+{\sin}^2\delta \right)}^{\raisebox{1ex}{$-1$}\!\left/ \!\raisebox{-1ex}{$2$}\right.} \)
(see Table 1 for notation, and Appendix for derivation of these equations).
From Eq. (1), the slip tendency is calculated as
$$ {T}_s=\frac{\sin \delta \cos \delta \kern0.5em \left(r-1\right)}{\sqrt{{\left(\frac{\cos \delta }{\cos \alpha}\right)}^2+{\sin}^2\delta}\kern0.5em \left[\left({\cos}^2\delta +{\sin}^2\alpha \kern0.3em {\sin}^2\delta \right)\left(1-\lambda \right)+{\left(\cos \alpha \kern0.3em \sin \delta \right)}^2\left(r-\lambda \right)\right]}. $$
Geometry of the plate boundary fault
Figure 6 shows the geometry of the picked plate boundary fault. The fault is shallowest (~ 3500 mbsf) near the SW edge and deepest (6500 mbsf) near the NE edge. On the NE to NW side, the dip angle is 10° to 15°, and the azimuth runs approximately NW, whereas on the NW side, the dip angle is approximately 25°, and the azimuth runs approximately N. On the NW side, the fault is generally shallower, which may be related to the shallow seafloor. Moreover, note that the strike of the outer ridge on the seafloor is ENE on the eastern side and changes to NE on the western side (Kimura et al. 2011).
a Dip angle and b dip azimuth along the plate boundary fault. Contours indicate the depth of fault. The white circle indicates the location of IODP Site C0002
On the seaward side, the fault geometry has a planer shape, where the dip angle is 30° to 40°, which is significantly larger than on the landward side, and the dip azimuth is NNW. This higher-angle structure appears to branch out from the low-angle reflector (as indicated above).
Principal stresses
Figure 7 shows the vertical stress (Sv) on the fault. Here, Sv is lowest around 105 MPa near the SW edge (3500 mbsf) and is higher than 160 MPa in the NE corner. The variation in Sv is mostly due to the variation of the fault depth. As mentioned in the "Introduction" section, Vp in the deeper region should have large uncertainties. The density anomaly, corresponding to a possible high Vp anomaly in the NE region above the fault (~ 1 km/s higher than the surrounding area; Shiraishi et al., this volume), is + 0.2 g/cm3. The density anomaly could produce an additional overburden of approximately 2 MPa (assuming that the Vp anomaly is 1 km thick).
Vertical stress, Sv, scaled between 90 MPa (purple) and 170 MPa (red). Contours indicate the depth of the fault (meters below the sea surface). Half-transparent colors on the inline PSDM image indicate the P-wave velocity scaled between ~ 3.5 km/s (blue) and 5 km/s (red). There is no vertical exaggeration. The yellow star at the drill site (Site C0002) indicates the current total drilled depth, and the white broken line indicates the planned drill path in 2018
The horizontal maximum stress (SH) is calculated using Eq. (6) and ranges from 100 to 200 MPa, which is attributed to the variation in Sv (note that r is kept constant in our model). As stated in the "Methods/Experimental" section, Sh = Sv (and thus σh = σv).
Slip tendency
Figure 8a through d indicates the effective normal stress, the shear stress, and the slip tendency of the fault. As described in the "Methods/Experimental" section, we used r = SH/Sv = 1.2 and λ* = 0.25 (λ = 0.7). For later discussion, the Ts value with λ* = 0.85 (λ = 0.94) is also plotted in Fig. 8.
a Effective normal stress image draped over the plate boundary fault. The image is scaled between 25 and 60 MPa. b Shear stress image scaled between 0 and 20 MPa. c Slip tendency image for λ* = 0.25 (λ = 0.7), scaled between 0 and 0.25. d Slip tendency image for λ* = 0.85 (λ = 0.94), scaled between 0 and 0.8. Note the difference in scale used for (c) and (d). The yellow stars on the drill site (Site C0002) in (a) and (b) indicate the current total drilled depths, and the white broken lines indicate the planned drill path in 2018
The effective normal stress ranges from 25 to 60 MPa. As recognized from Eq. (8), small values in the SW corner are due not only to the shallow depth (or small Sv) but also to the geometry (dip angle and azimuth). Note that we assume a near-hydrostatic pore fluid pressure based on Kitajima et al. (2017). For an overpressurized fault, the effective normal stress should decrease. The shear stress ranges from 0 to 20 MPa and depends on Sv, but is more sensitive to the variation in the dip angle than to Sv.
For a case with λ* = 0.25 (λ = 0.7) (Fig. 8), the slip tendency Ts varies from 0.05 to 0.1 on the N to NE side and from 0.15 to 0.2 in the center to the SW side. Moreover, Ts is larger than 0.2 on the SE side, mainly due to larger dip angle. For the case with λ* = 0.85 (λ = 0.94), Ts in the shallow dip region (N to NW side) varies from 0.2 to 0.6 and is larger than 0.6 on the SE side.
The slip tendency can also be shown on a lower-hemisphere equal-angle stereographic chart (Fig. 9), where the dip angle δ is shown along the radial axis, and the dip azimuth is shown along the azimuthal axis. Colored contours correspond to the amplitude of Ts with r = 1.2 and λ* = 0.25 (λ = 0.7). Black dots indicate the dip angle and azimuth distribution for the picked fault.
a Slip tendency Ts shown as colored contours on a lower-hemisphere equal-angle stereographic chart, where the dip angle δ is taken along the radial axis, and the dip azimuth is taken along the azimuthal axis. Color corresponds to the amplitude of Ts, with r = 1.2 and λ* = 0.25 (λ = 0.7). The black dots indicate the dip angle and azimuth distribution for the picked fault. The triangle corresponds to the fault geometry below Site C0002. The plate slip direction is equal to the SH direction (N55W). The slip direction of the 2016 Mie-ken Nanto-Oki event was obtained from the Japan Meteorological Agency (JMA) and U.S. Geological Survey (USGS). b The critical pore-pressure ratio λ*cr is the value at which the slip likelihood (SL) becomes 1, assuming that the intrinsic friction coefficient is 0.6 and Phyd/Sv = 0.6
The dots are distributed radially with most dip angles distributed at approximately 10°, where Ts is generally low (~ 0.1), and are distributed around NNW, where Ts is higher than 0.2. The triangle corresponds to the fault geometry below Site C0002, where dip angle is 6°, and thus Ts is smaller than 0.1. The color image in Fig. 9 shows the critical pore pressure ratio λ*cr at which the value of the slip likelihood (SL) becomes 1, assuming an intrinsic friction coefficient of 0.6 (e.g., Byerlee 1978). Here, λ*cr is larger than 0.85, where the dip angle is approximately 10° (suggesting that a higher pore fluid pressure is required for slip), whereas λ*cr is smaller, approximately 0.7 to 0.75, in the higher-angle region (most likely on the SE side).
Uncertainty in Ts
As shown in Eq. (10), the value of Ts depends on the fault dip angle and azimuth, SH/Sv (= r), as well as the pore pressure ratio (λ or λ*). Here, we estimate the uncertainties of these parameters and eventually assess the uncertainty in Ts.
The uncertainty due to picking the reflectors (± 50 m; see the "Methods/Experimental" section) leads to uncertainty in the dip angle, depending on the value of the dip angle. For a region with a shallow dip angle (approximately 10°), the uncertainty is ± 7°, resulting in Ts = 0.1 ± 0.05. For a steep dip region (approximately 30° in the SE corner), the uncertainty is ± 5°, resulting in Ts = 0.23 ± 0.05.
The uncertainty in the dip azimuth is difficult to assess, but based on its spatial distribution in Fig. 6, is probably approximately 20°. We recognize that Ts is not very sensitive to the variation in the dip azimuth (Fig. 9). We assess that the uncertainty in Ts due to that in the dip azimuth is less than 0.05.
The possible additional overburden of approximately 2 MPa due to high Vp above the fault would decrease r to 1.18 if SH remained unchanged. This might reduce Ts by 0.01 to 0.03 (Fig. 10). As stated in the previous section (Horizontal principal stresses), the stress ratio r should probably vary from 1.0 to 1.4. Figure 10 shows that the variation in Ts due to variation in r ranges from 0.01 to 0.16 for the nearly hydrostatic case (λ* = 0.25 or λ = 0.7), whereas Ts varies from 0.04 (r = 1) to 0.79 (r = 1.4) for the case of an overpressure (λ* = 0.85 or λ = 0.94).
Slip tendency Ts with respect to r = SH/Sv and λ = Pp/Sv for the fault geometry fixed at a dip angle of δ = 10° and α = 40°, where α is the angle between SH and the slip. The friction coefficient μ is set to be 0.6. The colder and warmer colors correspond to lower and higher values of Ts, respectively. Red indicates μ ≥ 0.6 (at which slip will occur)
Fault geometry is not optimally oriented
As described in the "Results" section, the peak population in the dip angle and the azimuth in Fig. 9 (dots in both panels) is significantly offset from the peak in Ts (color image). This indicates that most of the fault region is not optimally oriented. In other words, most regions of the fault plane must have some mechanism to "transform" their optimal dip (predicted from Ts) into the actual dip. In the following, we first discuss the offset in the dip angle, followed by the offset in the dip azimuth.
Fault dip angle is shallower than at the Ts peak
The main population of the fault dip angle is approximately 10° (except for ~ 30° in the SE corner), whereas the Ts peak [shown by the contours in Fig. 9a] corresponds to a dip angle of 30° to 40°. This can be reconciled if the frictional strength is weaker along the pre-existing, lower-angle (~ 10°) fault surface than the surrounding rock. Figure 11 shows the case in which the pre-existing fault (with dip angle δ2) is incipiently weaker than the host rock. For this case, if the pore fluid pressure increases (i.e., the Mohr circle moves leftward in Fig. 11) with time, then this pre-existing, low-angle thrust will be reactivated (leftmost circle).
Temporal change in stress and pore fluid pressure. μbed: incipient frictional coefficient of bedrock; μfault: incipient frictional coefficient of pre-existing fault with its dip angle δ2
In the present study, we assume that the principal stresses are oriented horizontally and vertically. This inference may be supported by the results of Ito et al. (2009), who estimates that σ1 (maximum principal stress) is oriented horizontally for the VLFEs that occurred near the trench axis in the forearc slope region off Kumano. However, in deeper fault zones, the σ1 axis can be tilted seaward. For the 2011 Tohoku earthquake case, Hasegawa et al. (2011) estimated the distributions of the P- and T-axes of the earthquake focal mechanisms before and after 2011 Tohoku earthquake. The distributions clearly indicate that most of the events have P-axes that dip seaward. If σ1 is tilted seaward by approximately 20° in the present study, the fault with a dip angle of approximately 10° will coincide with the Ts peak.
According to the critical taper model (e.g., Dahlen 1990), the stress state in accretionary wedges is controlled by tectonic/non-tectonic forces, such as a basal drag due to plate convergence, the gravitational force, and the compressional traction acting on the sidewall. The angle between σ1 and the slope on the seafloor can be derived theoretically. This angle is very sensitive to the effective friction along the basal décollement, as shown in Fig. 8 of Wang and Hu (2006). In the present study, the slope of the seafloor above the picked fault is less than 4°, and we assume that the pore fluid pressure ratio λ* is 0.85 (Tsuji et al. 2014), the intrinsic friction μ is 0.6 (Byerlee 1978), and the basal effective friction coefficient μb is 0.09. Figure 8 of Wang and Hu (2006) indicates that the corresponding angle of σ1 with respect to the seafloor is approximately 5°, or approximately horizontal. This may support our assumption that the principal stress is horizontal (σ1 = SH).
The shallow dip angle can be attributed to the weakness in the fault zone due to excess pore pressure. As shown in Fig. 8, if we use λ* = 0.85 (λ = 0.94), then the Ts value in the shallow dip region (N to NW side) increases from 0.05–0.2 to 0.2–0.6.
Fault dip azimuth is rotated from the plate convergence vector
The fault dip azimuth is N to NW, which is rotated clockwise from the plate convergence vector (N55W to WNW) or from the slip direction of the M6 earthquake (N70W; Wallace et al. 2016). This implies that the plate geometry is affected by factors other than the current stress state.
The shape of the subducting plate boundary for various scales has been revealed through seismic refraction (e.g., Nakanishi et al. 2008), seismic inversion (Kikuchi et al. 2003; Ichinose et al. 2003), tsunami inversion (Baba et al. 2002), and seismic tomography (Nakajima and Hasegawa 2007; Hirose et al. 2008; Yamamoto et al. 2017). Based on seismic analyses, Kodaira et al. (2006) reported a high-velocity doomed body to the south of Kii Peninsula, formed by breaking through the old accretionary prism sediment. A larger-scale geometry was compiled by F. Hirose (http://www.mri-jma.go.jp/Dep/sv/2ken/fhirose/ja/PHS.html) based on research by Baba et al. (2002), Nakajima and Hasegawa (2007), and Hirose et al. (2008). This geometry indicates that the subducting Philippine Sea Plate bends downward toward the Kii Peninsula. Other factors, such as the change in the plate convergence direction at approximately 3 Ma (Takahashi 2006), the hiatus in subduction at 13 to 6 Ma followed by a resumption in subduction (Tsuji et al. 2015), or the subduction of seamount (Kodaira et al. 2000; Kimura et al. 2014), should also be considered elsewhere.
Implications on the spatial variability of Ts
As defined in the "Methods/Experimental" section, we can assess the slip likelihood based on Ts and the intrinsic friction coefficient. Without knowledge of the friction coefficient, Ts may have two meanings, as follows.
Constant friction model
Assuming a constant friction everywhere on the fault plane, Ts exhibits a relative seismic likelihood. As indicated previously, Ts is highest on the high-angle fault in the SE corner, which is most likely to slip first, whereas other regions remain unslipped.
Simultaneous activation model
Assuming that slip is activated simultaneously over the fault plane, the variation of Ts should be compensated by the friction coefficient so that the slip likelihood value (SL) would be identical everywhere on the fault plane. As such, Ts is an indication of a relative strength, and the SE corner with highest Ts is likely to be strongest, resisting slip until the lower dip region of the fault reaches the critical state.
Tsuji et al. (2014) suggested a zone of very high pore pressure ratio beneath the plate boundary fault, based on the low P-wave velocity estimated from full-wave seismic inversion analysis (Kamei et al. 2012). In the Japan Trench, Ujiie et al. (2013) demonstrated the low friction of the fault material, either due to low intrinsic friction or to increased pore pressure. These inferences suggest that the low Ts on the main fault (excluding large-dip SE corner) corresponds to a reduction in effective friction strength either due to excess pore fluid pressure [larger λ value in Eq. (10)] or a low intrinsic friction coefficient. On the other hand, as shown in Fig. 7, the larger-dip SE region, which penetrates the high P-wave layer (Shiraishi et al., this volume), would probably have no (or a lower) excess pore fluid pressure, leading to a higher effective strength. Thus, we prefer model B to model A. However, model B must be tested through deep drilling at Site C0002.
Time-dependent Ts
Here, we consider how the value of Ts can evolve with time at a scale of less than 100 years. Among the factors affecting Ts (Sv, SH, and Pp), the overburden (Sv) remains approximately constant for this time scale. As such, we should consider an increase in SH due to tectonic loading, and an increase in Pp with time (e.g., Sibson 1992). These increases are translated as increases in r and λ (or λ*) in Eq. (10). The tectonic loading during the interseismic period would be a few to 10 MPa for subduction zone megathrusts, as inferred from the coseismic stress drop (e.g., Seno 2009) or a numerical simulation (Kinoshita et al. 2013). Accordingly, the value of r (SH/Sv) will increase from 1.2 (set for the present study) to approximately 1.3 for an Sv of approximately 100 MPa, in this case. The pore fluid pressure ratio can also increase above hydrostatic if any fluid flows into, or is produced in, the fault zone and remains undrained.
These temporal changes are illustrated in Fig. 11 in terms of Mohr's circle. Here, we assume that the incipient frictional coefficient of bedrock, μbed, is significantly larger than the incipient frictional coefficient of the pre-existing fault, μfault. The dip angle of the pre-existing fault is δ2. In stage A (post-/inter-seismic), the stress field is below the critical state (both for bedrock and a fault). As r and/or λ* increases with time, the diameter of Mohr's circle increases, and/or the circle moves to the left. In stage B, the circle touches the Coulomb failure line for the pre-existing fault, but the angle δ1 does not equal that for the fault (δ2). Thus, the fault does not reactivate yet. In stage C, in which the circle intersects the fault failure line with dip angle δ2, the pre-existing fault finally slips. However, the strength of the bedrock is much stronger than that of the fault. As such, the rock surrounding the fault will not break. Even if the strength of the bedrock is not very different from that for the fault, the bedrock will not break if its pore fluid pressure is not anomalous [in which case Mohr's circle remains at position (A)].
Based on this model, we investigate the temporal change in Ts with respect to r and λ* in Fig. 10. Here, the geometry is fixed at a dip of 10°, and the angle between SH (= slip direction) and the dip azimuth is fixed at 40°, because the geometry will not change with time. The color image indicates the magnitude of Ts. We recognize that Ts increases with increasing r and increasing λ*. If r increases from 1.0 to 1.4, then Ts will increase from 0.01 to 0.16 for a nearly hydrostatic case (λ* = 0.25 or λ = 0.7), whereas Ts increases from 0.04 to 0.79 for an overpressure case (λ* = 0.85 or λ = 0.94).
If we accept λ* = 0.85 on the fault (Tsuji et al. 2014) (λ ~ 0.94), the fault near Site C0002 is very close to failure for r > 1.2 (note that, here, the friction coefficient is assumed to be 0.6). In any case, drilling into the fault zone at Site C0002 will reveal its strength and stress. This will provide ground-truth evidence at this site, which can in turn be extrapolated to the 3D fault zone, and we can take an important step toward a better understanding of the slip likelihood of the Nankai seismogenic zone megathrust.
From the reprocessed 3D PSDM volume, we obtained a fine-scale geometry (dip angle and azimuth) of the plate boundary megathrust near the updip edge on the landward side of the Nankai seismogenic zone. This has enabled us to assess the degree to which the fault geometry affects the slip tendency.
The depth of the fault surface ranges from 3500 to 6500 mbsf. The geometry of the fault surface is divided into three groups: low dip (10°–15°) trending NW to N in the NE region, intermediate dip (approximately 25°) trending north in the western region, and high dip (30°–40°) trending NW in the SE corner. This higher-angle structure branches out from the low-angle reflector.
Shear and normal stresses and the slip tendency are estimated from fault dip and the azimuth for a geomechanical model with SH/Sv = 1.2, Sh = Sv, and Pp/Sv = 0.7 (the near-hydrostatic condition in which the normalized pore pressure ratio, λ*, is 0.25). The slip tendency is sensitive to the variation of fault geometry and is low in the low-angle dip region and high in the high-angle dip region. A low slip tendency region would correspond to a weaker region due to excess pore pressure or a low intrinsic friction coefficient, assuming that the entire fault is in a Coulomb-failure condition simultaneously.
Using the high λ* (~ 0.85) at the fault zone, the fault beneath IODP Site C0002 is very close to failure, assuming an intrinsic frictional coefficient of 0.6.
IODP:
mbsf:
Meters below sea floor
NanTroSEIZE:
Nankai Trough Seismogenic Zone Experiments
Ts :
Araki E, Saffer DM, Kopf AJ, Wallace LM, Toshinori Kimura T, Yuya Machida Y, Satoshi Ide S, Davis EE, IODP Expedition 365 shipboard scientists (2017) Recurring and triggered slow-slip events near the trench at the Nankai Trough subduction megathrust. Science 356: 1157–1160. doi:https://doi.org/10.1126/science.aan3120.
Baba T, Cummins PR, Hori T, Kaneda Y (2006) High precision slip distribution of the 1944 Tonankai earthquake inferred from tsunami waveforms: possible slip on a splay fault. Tectonophysics 426:119–134. https://doi.org/10.1016/j.tecto.2006.02.015.
Baba T, Tanioka Y, Cummins PR, Uhira K (2002) The slip distribution of the 1946 Nankai earthquake estimated from tsunami inversion using a new plate model. Phys Earth Planet Inter 132:59–73.
Brown ET, Hoek E (1978) Trends in relationship between measured in situ stresses and depth. Int J Rock Mech Min Sci Geomech Abst 15:211–215.
Byerlee J (1978) Friction of rocks. Pure Appl Geophys 116:615–626.
Dahlen FA (1990) Critical taper model of fold-and-thrust belts and accretionary wedges. Annual Rev Earth Planet Sci 18:55–99.
DeMets C, Gordon RG, Argus DF (2010) Geologically current plate motions. Geophysical J Int 181:1–80. https://doi.org/10.1111/j.1365-246X.2009.04491.x.
Erickson SN, Jarrard RD (1998) Velocity-porosity relationships for water-saturated siliciclastic sediments. J Geophys Res 103(B12):30,385–30,406. https://doi.org/10.1029/98JB02128.
Expedition 326 Scientists (2011) NanTroSEIZE stage 3: plate boundary deep riser: top hole engineering. IODP Preliminary Report 326. doi:https://doi.org/10.2204/iodp.pr.326.2011.
Hasegawa A, Yoshida K, Okada T (2011) Nearly complete stress drop in the 2011 Mw 9.0 off the Pacific coast of Tohoku earthquake. Earth Planets Space 63:703–707. https://doi.org/10.5047/eps.2011.06.007.
Henry P, Kanamatsu T, Moe KT (2012) the Expedition 333 Scientists Proc IODP (Vol. 333) Tokyo. Integrated Ocean Drilling Program Management International Inc. https://doi.org/10.2204/iodp.proc.333.2012.
Hirose F, Nakajima J, Hasegawa A (2008) Three-dimensional seismic velocity structure and configuration of the Philippine Sea slab in southwestern Japan estimated by double-difference tomography. J Geophys Res 113:B09315. https://doi.org/10.1029/2007JB005274.
Ichinose GA, Thio HK, Somerville PG, Sato T, Ishii T (2003) Rupture process of the 1944 Tonankai earthquake (Ms 8.1) from the inversion of teleseismic and regional seismograms. J Geophys Res 108 (B10):2497. doi:https://doi.org/10.1029/2003JB002393, 2003.
Ito Y, Asano Y, Obara K (2009) Very low-frequency earthquakes indicate a transpressional stress regime in the Nankai accretionary prism. Geophys Res Lett 36:L20309. https://doi.org/10.1029/2009GL039332.
Kamei R, Pratt RG, Tsuji T (2012) Waveform tomography imaging of a mega-splay fault system in the seismogenic Nankai subduction zone. Earth Planet Sci Lett 317–318:343–353.
Kikuchi M, Nakamura M, Yoshikawa K (2003) Source rupture processes of the 1944 Tonankai earthquake and the 1945 Mikawa earthquake derived from low-gain seismograms. Earth Planets and Space 55:159–172.
Kimura G, Hashimoto Y, Kitamura Y, Yamaguchi A, Koge H (2014) Middle Miocene swift migration of the TTT triple junction and rapid crustal growth in Southwest Japan: a review. Tectonics 33:1219–1238. https://doi.org/10.1002/2014TC003531.
Kimura G, Moore GF, Strasser M, Screaton E, Curewitz D, Streiff C, Tobin H (2011) Spatial and temporal evolution of the megasplay fault in the Nankai Trough. Geochem Geophys Geosyst 12:Q0A008. https://doi.org/10.1029/2010GC003335.
Kinoshita M, Kimura G, Saito S (2014) Chapter 4.4.2 Seismogenic processes revealed through the Nankai Trough seismogenic zone experiments: Core, log, geophysics and observatory measurements. In: Stein R, Blackman DK, Inagaki F, Larsen H-C (eds) Developments in Marine Geology Vol 7: Earth and Life Processes Discovered from Subseafloor Environment. Elsevier, pp 641–670. https://doi.org/10.1016/B978-0-444-62617-2.00021-9.
Kinoshita M, Tobin H, Ashi J, Kimura G, Lallement S, Screaton EJ, Curewitz D, Masago H, Moe KT (2009) The Expedition 314/315/316 Scientists Proc. IODP, 314/315/316. (Integrated Ocean Drilling Program Management International, Inc.), Washington, DC. https://doi.org/10.2204/iodp.proc.314315316.2009.
Kinoshita, M., Tobin, H.J. (2013) Interseismic stress accumulation at the locked zone of Nankai Trough seismogenic fault off Kii Peninsula, Tectonophysics, 600C, 153–164. https://doi.org/10.1016/j.tecto.2013.03.015
Kitajima H, Saffer D, Sone H, Tobin H, Hirose T (2017) In situ stress and pore pressure in the deep interior of the Nankai accretionary prism, integrated ocean drilling program site C0002. Geophy Res Lett 44:9644–9652. https://doi.org/10.1002/2017GL075127.
Kodaira S, Hori T, Ito A, Miura S, Fujie G, Park JO, Baba T, Sakaguchi H, Kaneda Y (2006) A cause of rupture segmentation and synchronization in the Nankai trough revealed by seismic imaging and numerical simulation. J Geophys Res 111:B09301. https://doi.org/10.1029/2005JB004030.
Kodaira S, Takahashi N, Nakanishi A, Miura S, Kaneda Y (2000) Subducted seamount imaged in the rupture zone of the 1946 Nankaido earthquake. Science 289:104–106.
Kopf A, Araki E, Toczko S, (2011) Expedition 332 Scientists NanTroSEIZE Stage 2: riserless observatory. International Ocean Discovery Program Prel Rept 332. doi:https://doi.org/10.2204/iodp.pr.332.2011.
Kopf A, Saffer D, Toczko S, (2016) Expedition 365 Scientists Expedition 365 Preliminary Report: NanTroSEIZE Stage 3: Shallow Megasplay Long-Term Borehole Monitoring System (LTBMS). International Ocean Discovery Program Prel Rept 332 doi:https://doi.org/10.14379/iodp.pr.365.2016.
Miyakawa A, Otsubo M (2015) Applicability of slip tendency for understanding long-term fault activity: a case study of active faults in northeastern Japan. J JSCE 3:105–114. https://doi.org/10.2208/journalofjsce.3.1_105.
Miyazaki S, Heki K (2001) Crustal velocity field of Southwest Japan: subduction and arc-arc collision. J Geophys Res 106:4305–4326. https://doi.org/10.1029/2000JB900312.
Moore GF, Bangs NL, Taira A, Kuramoto S, Pangborn E, Tobin HJ (2007) Three-dimensional splay fault geometry and implications for tsunami generation. Science 318:1128–1131. https://doi.org/10.1126/science.1147195.
Morris A, Ferrill DA, Henderson DB (1996) Slip-tendency analysis and fault reactivation. Geology 24:275–278. https://doi.org/10.1130/0091-7613(1996)024<0275:STAAFR>2.3.CO;2.
Nakajima J, Hasegawa A (2007) Subduction of the Philippine Sea plate beneath southwestern Japan: slab geometry and its relationship to arc magmatism. J Geophys Res 112:B08306. https://doi.org/10.1029/2006JB004770.
Nakanishi A, Kodaira S, Miura S, Ito A, Sato T, Park JO, Kido Y, Kaneda Y (2008) Detailed structural image around splay-fault branching in the Nankai sub-duction seismogenic zone: results from a high-density ocean bottom seismic survey. J Geophys Res 113:B03105. https://doi.org/10.1029/2007JB004974.
Nakano M, Hyodo M, Nakanishi A, Yamashita M, Hori T, Kamiya S, Suzuki K, Tonegawa T, Kodaira S, Takahashi N, Kaneda Y (2018) The 2016 off Southeast Mie prefecture earthquake (mw=5.9) as an indicator of preparatory processes of the next Nankai trough megathrust earthquake. Progress in earth and planetary science. Accepted in July 2018.
Obana K, Kodaira S (2009) Low-frequency tremors associated with reverse faults in a shallow accretionary prism. Earth Planet Sci Lett 287:168–174. https://doi.org/10.1016/j.epsl.2009.08.005.
Obara K, Ito Y (2005) Very low frequency earthquakes excited by the 2004 off the Kii peninsula earthquakes: a dynamic deformation process in the large accretionary prism. Earth Planets Space 57:321–326. https://doi.org/10.1186/BF03352570.
Obara K, Kato A (2016) Connecting slow earthquakes to huge earthquakes. Science 353:253–257. https://doi.org/10.1126/science.aaf1512.
Park JO, Fujie G, Wijerathne L, Hori T, Kodaira S, Fukao F, Moore GF, Bangs NL, Kuramoto S, Taira A (2010) A low velocity zone with weak reflectivity along the Nankai subduction zone. Geology 38:283–286. https://doi.org/10.1130/G30205.1.
Rummel F, Möhring-Erdmann G, Bäumgartner J (1986) Stress constraints and Hydrofracturing stress data for the continental crust. PAGEOPH 124:875–894.
Saffer D, McNeill L, Byrne T, Araki E, Toczko S, Eguchi N et al (2010) Proc IODP (Vol. 319) Tokyo: Integrated Ocean Drilling Program management international Inc. https://doi.org/10.2204/iodp.proc.319.104.2010.
Saito S, Underwood MB, Kubo Y (2010) The expedition 322 scientists proc IODP (Vol. 322) Tokyo. Integrated Ocean drilling program management international, Inc. https://doi.org/10.2204/iodp.proc.322.2010.
Screaton EJ, Kimura G, Curewitz D, (2009) et al. The expedition 316 scientists expedition 316 summary. In M. Kinoshita, H. Tobin, J. Ashi, G. Kimura, S. Lallemant, E. J. Screaton. (Eds.), Proc. IODP (Vol. 314/315/316). Washington, DC: Integrated Ocean Drilling Program Management International Inc doi:https://doi.org/10.2204/iodp.proc.314315316.131.2009.
Seno T (2009) Determination of the pore fluid pressure ratio at seismogenic megathrusts in subduction zones: implications for strength of asperities and Andean-type mountain building. J Geophys Res 114:B05405. https://doi.org/10.1029/2008JB005889.
Shiraishi K, Moore GF, Yamada Y, Kinoshita M, Sanada Y, Kimura G (2016) Improved 3D seismic images of dynamic deformation in the Nankai trough off Kumano, abstract T42C-06 presented at 2016 fall meeting, AGU, San Francisco, 12-16, Dec. Seismogenic zone structures revealed by improved 3D seismic image in the Nankai Trough off Kumano. in review SPEPS..
Sibson RH (1992) Implications of fault-valve behavior for rupture nucleation and recurrence. Tectonophysics 211:283–293. https://doi.org/10.1016/0040-1951(92)90065-E.
Strasser M, Dugan B, Kanagawa K, Moore GF, Toczko S, Maeda L et al (2014) Proc IODP (Vol. 338) Yokohama: Integrated Ocean Drilling Program 12: Q0AD17. https://doi.org/10.2204/iodp.proc.338.2014.
Sugioka H, Okamoto T, Nakamura T, Ishihara Y, Ito A, Obana K, Kinoshita M, Nakahigashi K, Shinohara M, Fukao Y (2012) Tsunamigenic potential of the shallow subduction plate boundary inferred from slow seismic slip. Nat Geosci 5:414–418. https://doi.org/10.1038/NGEO1466.
Takahashi M (2006) Tectonic development of Japanese Islands controlled by Philippine Sea plate motion (in Japanese with English abstract). J Geogr 115:116–123.
Templeton EL, Rice JR (2008) Off-fault plasticity and earthquake rupture dynamics: 1. Dry materials or neglect of fluid pressure changes. J Geophys Res 113:B09306. https://doi.org/10.1029/2007JB005529.
Tobin H, Henry P, Vannucchi P, Screaton E (2014) Chapter 4.4.1 Subduction zones: structure and deformation history. In: Stein R, Blackman DK, Inagaki F, Larsen H-C (eds) Developments in marine geology Vol 7: earth and life processes discovered from subseafloor environment. Elsevier, pp 599–640. https://doi.org/10.1016/B978-0-444-62617-2.00020-7.
Tobin H, Hirose T, Saffer D, Toczko S, Maeda L, Kubo Y, Boston B, Broderick A, Brown K, Crespo-Blanc A, Even E, Fuchida S, Fukuchi R, Hammerschmidt S, Henry P, Josh M, Jurado MJ, Kitajima H, Kitamura M, Maia A, Otsubo M, Sample J, Schleicher A, Sone H, Song C, Valdez R, Yamamoto Y, Yang K, Sanada Y, Kido Y, Hamada Y (2015) Site C0002. In: Tobin H, Hirose T, Saffer D, Toczko S, Maeda L, Kubo Y (eds) Expedition 348 scientists, proc IODP 348. TX (Integrated Ocean Drilling Program), College Station. https://doi.org/10.2204/iodp.proc.348.103.2015.
Tobin, H., Kinoshita, M., Ashi, J., Lallemant, S., Kimura, G., Screaton, E., (2009) et al.n.d.. NanTroSEIZE stage 1 expeditions 314, 315, and 316: first drilling program of the Nankai trough seismogenic zone experiment. Sci Drill 8: 4–17. https://doi.org/10.2204/iodp.sd.8.01.2009.
Tobin HJ, Kinoshita M (2006) NanTroSEIZE: the IODP Nankai trough seismogenic zone experiment. Sci Drill 2:23–27. https://doi.org/10.2204/iodp.sd.2.06.2006.
Tong H, Yin A (2011) Reactivation tendency analysis: a theory for predicting the temporal evolution of preexisting weakness under uniform stress state. Tectonophysics 503:195–200. https://doi.org/10.1016/j.tecto.2011.02.012.
Tsuji T, Ashi J, Strasser M, Kimura G (2015) Identification of the static backstop and its influence on the evolution of the accretionary prism in the Nankai trough. Earth Planet Sci Lett 431:15–25. https://doi.org/10.1016/j.epsl.2015.09.011.
Tsuji T, Kamei R, Pratt RG (2014) Pore pressure distribution of a mega-splay fault system in the Nankai trough subduction zone: insight into up-dip extent of the seismogenic zone. Earth Planet Sci Lett 396:165–178. https://doi.org/10.1016/j.epsl.2014.04.011.
Ujiie K, Tanaka H, Saito T, Tsutsumi A, Mori JJ, Kameda J, Brodsky EE, Chester FM, Eguchi N, Toczko S, Expedition 343 and 343T Scientists (2013) Low Coseismic shear stress on the Tohoku-Oki megathrust determined from laboratory experiments. Science 342:1211–1214. https://doi.org/10.1126/science.1243485.
Wallace LM, Araki E, Saffer D, Wang X, Roesner XA, Kopf A, Nakanishi A, Power W, Kobayashi R, Kinoshita C, Toczko S, Kimura T, Machida Y, Carr S (2016) Near-field observations of an offshore Mw 6.0 earthquake from an integrated seafloor and subseafloor monitoring network at the Nankai Trough, Southwest Japan. J Geophys Res Solid Earth 121:8338–8351. https://doi.org/10.1002/2016JB013417.
Wang K, Hu Y (2006) Accretionary prisms in subduction earthquake cycles: the theory of dynamic coulomb wedge. J Geophys Res 111:B06410. https://doi.org/10.1029/2005JB004094.
Yamamoto Y, Takahashi T, Kaiho Y, Obana K, Nakanishi A, Kodaira K, Kaneda Y (2017) Seismic structure off the Kii peninsula, Japan, deduced from passive-and active-source seismographic data. Earth Planet Sci Lett 461:163–175. https://doi.org/10.1016/j.epsl.2017.01.003.
Yokota Y, Ishikawa T, Watanabe S, Tashiro T, Asada A (2016) Seafloor geodetic constraints on interplate coupling of the Nankai trough megathrust zone. Nature 534:374–377. https://doi.org/10.1038/nature17632.
Zoback M (2007) Reservoir Geomechanics. Cambridge University Press. https://doi.org/10.1017/CBO9780511586477.
The authors would like to thank Gaku Kimura, Yasuhiro Yamada, Eiichiro Araki, and Kyuichi Kanagawa for their fruitful discussion. MK gratefully acknowledges Kotaro Fujita for his assistance in finalizing the figures.
The present study was supported by JSPS KAKENHI Grant Number JP15H05717, and by the Ministry of Education, Culture, Sports, Science and Technology (MEXT) through management expenses grants at the Research and Development Center for Ocean Drilling Science (ODS) and the Center for Deep Earth Exploration (CDEX) of Japan Agency for Marine-Earth Science and Technology (JAMSTEC).
The dataset supporting the conclusions of the present study is included in the manuscript and associated Additional file 1. For those who are interested in the original 3D seismic volume, please contact Masataka Kinoshita (masa@eri.u-tokoyo.ac.jp).
Earthquake Research Institute, The University of Tokyo, 1-1-1 Yayoi, Bunkyo-ku, Tokyo, 113-0032, Japan
Masataka Kinoshita
Research and Development Center for Ocean Drilling Science (ODS), Japan Agency for Marine-Earth Science and Technology (JAMSTEC), 3173-25 Showa-machi, Kanazawa-ku, Yokohama, Kanagawa, 236-0001, Japan
Kazuya Shiraishi
Kingston University, Kingston Upon Thames, UK
Evi Demetriou
Faculty of Science and Technology, Kochi University, Kochi, 780-8072, Japan
Yoshitaka Hashimoto
Graduate School of Engineering, Kyoto University, Kyoto, 606-8501, Japan
Weiren Lin
MK proposed the method and directed the overall project. KS processed the 3D seismic volume data and contributed to seismic interpretation. ED actually picked the fault surface under the supervision of MK. YH contributed to the geological/geotechnical interpretation of Ts contrast. WL helped in setting the regional principal stresses. All of the authors read and approved the final manuscript.
Correspondence to Masataka Kinoshita.
MK has been committed to the IODP NanTroSEIZE project, which was started in 2007, as a key member of the project coordination team.
Dataset supporting the conclusion. (ZIP 1803 kb)
Calculation of normal/shear stress and slip tendency
First, we define the y-axis as being parallel to σ1 = σH, the x-axis as being parallel to σh, and the z-axis as being parallel to σv. Considering the effective principal stresses, the stress tensor ∑ is expressed as
$$ \Sigma =\left[\begin{array}{ccc}{\sigma}_h& 0& 0\\ {}0& {\sigma}_H& 0\\ {}0& 0& {\sigma}_v\end{array}\right]=\left[\begin{array}{ccc} Sh- Pp& 0& 0\\ {}0& SH- Pp& 0\\ {}0& 0& {S}_v- Pp\end{array}\right]. $$
We set SH oriented to N55W, so that the angle (α) between SH (= slip) and the dip azimuth (ψ) at each point on the fault surface is defined as α = ψ – (− 55). Based on these definitions, the normal vector at each point on the fault is obtained as
$$ \mathbf{n}=\left[\begin{array}{c}\sin \alpha \sin \delta \\ {}\cos \alpha \sin \delta \\ {}\cos \delta \end{array}\right]. $$
The stress vector (p) acting at a point on the fault surface by the stress tensor ∑ is derived from Cauchy's formula:
$$ \mathbf{p}=\mathbf{\sum}\cdotp \mathbf{n}. $$
Thus, the normal stress (σn) acting at any point on the fault surface can be derived by taking the inner product of p and n.
For the shear stress, we need to define the direction of slip, which we assume to be parallel to the SH direction (N55 W), i.e., the unit vector parallel to the shear stress (τ) is
$$ \boldsymbol{t}=\left[\begin{array}{c}0\\ {}-A\cos \delta \kern0.3em {\cos}^{-1}\alpha \\ {}A\sin \delta \end{array}\right], $$
where \( A={\left(\frac{\cos^2\delta }{\cos^2\alpha }+{\sin}^2\delta \right)}^{\raisebox{1ex}{$-1$}\!\left/ \!\raisebox{-1ex}{$2$}\right.} \).
From (A1) and (A2), the effective normal stress is expressed as
$$ {\boldsymbol{\sigma}}_{\boldsymbol{n}\_\boldsymbol{e}}={\boldsymbol{n}}^T\bullet \left(\boldsymbol{\varSigma} \cdotp \boldsymbol{n}\right)={\left(\sin \alpha \sin \delta \right)}^2{\sigma}_h+{\left(\cos \alpha \sin \delta \right)}^2{\sigma}_H+{\left(\cos \delta \right)}^2{\sigma}_v. $$
Since we assume a triaxial condition (σh = σv), we have
$$ {\boldsymbol{\sigma}}_{\boldsymbol{n}\_\boldsymbol{e}}=\left({\left(\sin \alpha \sin \delta \right)}^2+{\left(\cos \delta \right)}^2\right){\sigma}_v+{\left(\cos \alpha \sin \delta \right)}^2{\sigma}_H. $$
The shear stress τ is
$$ \uptau ={\boldsymbol{t}}^T\bullet \left(\boldsymbol{\varSigma} \cdotp \boldsymbol{n}\right)=A\ \sin \delta \cos \delta \left({\sigma}_H-{\sigma}_v\right). $$
The shear stress is oriented in part seaward due to the seaward dipping region. This is likely due to artificial noise that could not be removed from the migration noise reduction process. In any case, we exclude such cases and set τ = 0 if τ < 0.
Slip tendency Ts is calculated by dividing (A4) by (A5):
$$ {T}_s=\frac{\sin \delta \cos \delta \left({\sigma}_H-{\sigma}_v\right)}{\sqrt{{\left(\frac{\cos \delta }{\cos \alpha}\right)}^2+{\sin}^2\delta}\kern0.5em \left[\left({\cos}^2\delta +{\sin}^2\alpha \kern0.3em {\sin}^2\delta \right){\sigma}_v+{\left(\cos \alpha \kern0.3em \sin \delta \right)}^2{\sigma}_v\right]}. $$
Here, we rewrite (A4) through (A6) by introducing r (= SH/Sv) and λ (= Pp/Sv):
$$ {\sigma}_{n\_e}=\left[\ \left({\left(\sin \alpha \sin \delta \right)}^2+{\left(\cos \delta \right)}^2\right)\left(1-\lambda \right)+{\left(\cos \alpha \sin \delta \right)}^2\left(r-\lambda \right)\ \right]\ {S}_v, $$
$$ \tau =A\ \sin \delta \cos \delta \left(r-1\right)\ {S}_v, $$
$$ {T}_s=\frac{\sin \delta \cos \delta \left(r-1\right)}{\sqrt{{\left(\frac{\cos \delta }{\cos \alpha}\right)}^2+{\sin}^2\delta}\kern0.5em \left[\left({\cos}^2\delta +{\sin}^2\alpha \kern0.3em {\sin}^2\delta \right)\left(1-\lambda \right)+{\left(\cos \alpha \kern0.3em \sin \delta \right)}^2\left(r-\lambda \right)\right]}. $$
Kinoshita, M., Shiraishi, K., Demetriou, E. et al. Geometrical dependence on the stress and slip tendency acting on the subduction megathrust of the Nankai seismogenic zone off Kumano. Prog Earth Planet Sci 6, 7 (2019). https://doi.org/10.1186/s40645-018-0253-y
Nankai Trough seismogenic zone
3D seismic volume
Stress tensor
4. Solid earth sciences
Subduction-zone megathrust earthquakes: New perspectives from insitu data & laboratory analyses