Blog

Upper-plate conduits linked to plate boundary that hosts slow earthquakes | Nature Communications

Thank you for visiting nature.com. You are using a browser version with limited support for CSS. To obtain the best experience, we recommend you use a more up to date browser (or turn off compatibility mode in Internet Explorer). In the meantime, to ensure continued support, we are displaying the site without styles and JavaScript.

Nature Communications volume  14, Article number: 5101 (2023 ) Cite this article 1hz 3c Geophone

Upper-plate conduits linked to plate boundary that hosts slow earthquakes | Nature Communications

In shallow subduction zones, fluid behavior impacts various geodynamic processes capable of regulating slip behaviors and forming mud volcanoes. However, evidence of structures that control the fluid transfer within an overriding plate is limited and the physical properties at the source faults of slow earthquakes are poorly understood. Here we present high-resolution seismic velocity models and reflection images of the Hyuga-nada area, Japan, where the Kyushu-Palau ridge subducts. We image distinct kilometer-wide columns in the upper plate with reduced velocities that extend vertically from the seafloor down to 10–13 km depth. We interpret the low-velocity columns as damaged zones caused by seamount subduction and suggest that they serve as conduits, facilitating vertical fluid migration from the plate boundary. The lateral variation in upper-plate velocity and seismic reflectivity along the plate boundary correlates with the distribution of slow earthquakes, indicating that the upper-plate drainage system controls the complex pattern of seismic slip at subduction faults.

Plate boundaries in subduction zones portray a wide spectrum of slip behavior, from day-to-year-long slow slip events to instantaneous faulting, governed by diverse physical factors. Pore fluid is one of the most influential parameters of seismic slips, and pore fluid behavior has been intensively studied in recent years1,2,3. High fluid pressure is a prerequisite condition for slow earthquakes to nucleate at subduction faults4. At shallow parts of subduction zones (depth < 15 km), the pore fluid pressure is primarily controlled by the compaction/dehydration reactions of sediments entering the subduction zone3. Once released from the source, the fluid migrates upward via permeable zones into the forearc wedge owing to its buoyancy5,6, forming mud volcanoes and cold seeps on the forearc basins and trench slopes7. Subducting seamounts (topographic highs at the plate boundary) impact the hydrological system in subduction zones by causing a heterogeneous distribution of permeability in the upper plate and forming a subduction channel filled with fluid-rich sediments8,9,10. However, the evidence on fluid distribution directly linked to subducting seamounts is limited. In particular, the along-strike extent of the impact of the subducting seamount is poorly constrained. In addition, our understanding of how the complex structure formed by subducting seamounts modulates the slip pattern at subduction faults remains unclear11.

The Hyuga-nada subduction zone in southwest Japan is an excellent area to study the relationship between subduction earthquakes and forearc hydrology. Adjacent to the Nankai Trough, where megathrust earthquakes of magnitudes greater than 8 occur frequently, the Hyuga-nada area historically hosted several interplate earthquakes of magnitudes ~712,13, the source areas of which are delineated by long-term slow slip events on the downdip side14 (Fig. 1a). On its updip side, slow earthquakes have been well-documented in recent seafloor observations15,16,17,18. For example, the activity of low-frequency tremors exhibiting rapid epicentral migration (with a speed of tens to hundreds of kilometers per day)15,16 implies that these events occur at the plate interface, driven by the spatiotemporal variation in the fluid pressure therein19. Another fluid-related phenomenon in the Hyuga-nada area is abundant mud volcanoes on the seafloor20,21,22. Although geochemical analyses of core samples suggested that these mud volcanoes are fed by the dehydration of clay minerals at shallow crustal depths21, the origin of the fluid is poorly understood. Moreover, the hydrological system responsible for the transfer of deep fluids is unknown.

a Regional bathymetric map of the study area. Purple, blue, and green shaded areas indicate the source areas of large megathrust earthquakes in 194647, 196812, and 199613, respectively Orange shades indicate the source areas of slow slip events14,48. White dots, yellow squares, and red triangles denote the locations of low-frequency tremors14, repeating earthquakes38, and mud volcanoes20,21,22, respectively. The convergence direction of the Philippine Sea Plate relative to the Amurian plate is shown by the white arrow28. The thin black line indicates the seismic profile of this study, and the thick black line indicates the range shown in Figs. 2, 4 and 5. The black box indicates the location enlarged in panel b. b Layout of the ocean bottom seismographs (OBSs) and multichannel seismic (MCS) reflection profiles superimposed on the reduced-to-pole (RTP) magnetic anomaly map49. The yellow circles indicate the locations of OBSs. The MCS reflection data were acquired along the black line. The velocity model and reflection images shown in Fig. 2 correspond to the 100-km-profile where the OBSs were deployed.

To better understand the geodynamic system in the Hyuga-nada subduction zone, we performed a dense seismic refraction experiment using ocean bottom seismographs (OBSs) and a multichannel seismic (MCS) reflection survey in the summer of 2020. For the refraction study, we laid out 50 OBSs at 2-km intervals on the seismic line (HYU01) (Fig. 1b). The seismic line, aligned parallel to the regional trend of the Nankai Trough in the SW-NE direction, crossed the source area of the low-frequency tremor distributed between the shallow and deep seamounts (Fig. 1). In this study, we construct a P-wave velocity (Vp) model by applying the advanced full waveform inversion (FWI) technique to the OBS data (Methods section). Our results provide important structural constraints of the crust from the seafloor down to the subducting plate at a significantly higher spatial resolution than previous studies23,24,25. We compare the Vp model with the reflection image on the coincident line for structural interpretations. Based on the seismic results, we suggest the origin of deep fluids and their pathways up to the seafloor mud volcanoes. We also compare them with the spatial distribution of slow and regular earthquakes to discuss the structural controls on their occurrence. The structural characterization on the subduction megathrust faults using our study method can provide important information to evaluate the spatial extent of future megathrust earthquakes in the region.

The FWI Vp modeling along the HYU01 profile reveals a remarkably heterogeneous structure in the upper plate (Fig. 2a). The most important finding is the lateral variation in the Vp, highlighted by horizontally alternating high- and low-velocity anomalies (Fig. 2b), which is supported by the varying apparent velocities of first-arrival refraction phases in the OBS record sections (Supplementary Figs. 1–8). We detected several low-velocity columns having Vp values < 3.5 km/s that extended near-vertically from the seafloor to depths of over 10 km. These upper-plate velocities are significantly lower than those in other areas of the Nankai subduction zone located at similar distances from the trench axis (~50 km)25,26, suggesting that the subduction of the Kyushu-Palau Ridge increases the porosity of the forearc wedge over a broad region (about four times broader than the ridge itself). Interestingly, some of the low-velocity columns dip slightly toward the summit of the ridge located in the middle of the study area. Based on their fan-shaped geometry radiating from the seamount, we interpret that the low-velocity features represent an intensive fracture network; previous studies utilizing sandbox experiments explain that such networks consisting of subvertical faults can develop above subducting seamounts27. The dense distribution of fracture zones in the vicinity of the Kyushu-Palau Ridge may reflect the stress concentration in the region caused by seamount subduction, and their lower seismic velocities in the western part (compared with those in the eastern part; Fig. 2d–h) may indicate higher porosity and be attributed to the regional stress regime governed by the northwestward convergence of the incoming Philippine Sea Plate28.

a Full waveform inversion (FWI) P-wave velocity (Vp) model highlighting the upper-plate lower-velocity columns (shown by black arrows). Areas with poor checkerboard recovery (Supplementary Fig. 10) are shaded. b Vp perturbation from the 1D depth average (depth from the seafloor). c Prestack depth-migrated seismic reflection image. d–h Close-up views of the low-velocity columns. i–m Close-up reflection images, with a focus on the domed structures atop the low-velocity columns.

In enlarged views of the seismic reflection image, the upper-plate fracture zones form conical zones with atypical reflection patterns characterized by inverted-V-shaped reflectors (Fig. 2c). The horizontal reflectors corresponding to the sedimentary layers are highly disturbed and portray upwardly domed geometry (Fig. 2i–m) on the top of the low-velocity fracture zones (Fig. 2d–h). These seismic reflection features are typical of mud volcanoes29 and suggest the occurrence of active upwelling flows inside the low-velocity columns. The seafloor bathymetry of the area shows several circular mounds (with steep slopes) close to the seismic profile, indicating another typical feature of mud volcanoes (Fig. 3). Additionally, the data suggest an extremely complex pattern of fault development, characterized by several linear fault scarps. Based on the consistency with experimental studies27, we suggest that these bathymetric features are associated with the ridge subduction. Importantly, the upper-plate low-velocity columns are situated where these linear scarps are developed (Fig. 3). These observational facts demonstrate that the upper plate accommodates well-developed fluid conduits that efficiently facilitate fluid flow to the seafloor and feed the mud volcanoes.

Background colors represent seafloor depth (panel a) and slope gradient (panel b) highlighting circular mounds that possibly indicate mud volcanoes (blue arrows) and linear scarps (magenta arrows) that radiate from the seamounts (marked by white/yellow contour lines of the RTP magnetic data). The black/red thick lines indicate the location of the low-velocity columns in the upper plate (Fig. 2). To create this plot, we compiled the high-resolution bathymetry data obtained by multiple cruises using the Japan Agency for Marine-Earth Science and Technology (JAMSTEC) research vessels. The slope gradient of the seafloor was calculated using a generic mapping tool50.

The seismic reflection image depicts reflective zones that extend laterally at depths of 10–13 km in the middle and eastern part of the seismic profile (40–100 km distance in Fig. 2c), with higher reflectivity in the flanks of the subducting ridge (40–45 km distance and 55–72 km distance in Fig. 2c) and away from the ridge (90–100 km distance in Fig. 2c). The reflective zones mostly trace the iso-Vp contour of 3.5–4.5 km/s (Fig. 2a) and also correspond to the areas having a positive velocity gradient in the depth direction, underlying those having a negative velocity gradient (Fig. 4). These layer structures probably represent low-velocity underthrust sediments atop the upper crust of the subducted Kyushu-Palau Ridge. In general, spatial variations in seismic reflectivity at the plate boundary is closely correlated with the physical properties and plate coupling conditions, which are highly influenced by in-situ fluid pressure30,31. At shallow depths where the temperature is below 150 °C, the pore fluid pressure is primarily controlled by the compaction and dehydration reactions of clay-rich sediments3. Notably, Vp values of 3.5–4.5 km/s imply that the reflective zones may be comprised of volcaniclastic sediments, which can hold a substantial amount of fluid before subduction and expel it when subducted10,32,33. Faint seismic reflections are observed in the western portion of the plate boundary. The disappearance of the reflection may suggest that fluids are starved there. Correspondingly, the westernmost low-velocity column is limited to the upper 8 km and unlikely to reach the plate boundary. The portion between the upper-plate low-velocity zones and the plate boundary increases the Vp to over 4.5 km/s and exhibits a sporadic reflection pattern with slightly higher reflectivity (0–40 km distance and 10–13 km depth in Fig. 2a and c).

This figure highlights velocity reversals, or low-velocity zones (marked by red-colored areas), at different depths. BSR indicates bottom simulating reflectors.

Submarine mud volcanoes are a well-known feature in forearc basins and provide a window to explore the subsurface hydrological system in subduction zones34. One of the key issues regarding mud volcanoes is the diversity of the fluid origins. Even though previous studies state that the fluids originate from the dehydration of clay-rich sediments within the accretionary prisms, which typically occurs at depths of several kilometer by rapid sedimentation and compaction34, recent studies report other potential sources, such as subducted sediments along deep thrust faults35 and forearc mantle wedges36. The seismic data recorded in our study offer new evidence for the previously derived conclusions that the fluid source for the mud volcanoes in the Hyuga-nada subduction zone is along the plate boundary at depths >10 km.

At the shallow parts of subduction zones, fluid is likely to migrate via the decollement and branching faults, because they act as the most permeable pathways for fluid flow. A significant difference in the fluid feeding system in the Hyuga-nada area from this conventional concept is that the fluids may be ascending upward through the near-vertical conduits. We suggest that the fracture network in the Hyuga-nada subduction zone is primarily paved by the subducting seamount; the Vp model indicates that the network consists of discrete low-velocity columns that are densely distributed atop and around the subducting ridge (Fig. 2a and b). Our results also reveal some velocity reversals at different depths within the upper plate (Fig. 4). These results suggest that fluids from the plate boundary may be captured and accumulate at multiple depths within the upper plate and/or the upper plate itself also generates some fluids by dehydration of clay-rich sediments37. The latter possibility is also consistent with the observation that the upper-plate conduits are linked to not only the fluid-rich plate boundary (40–72 km and 80–100 km distances in Fig. 5a) but also the fluid-starved portions (30–40 km and 72–80 km distance in Fig. 5a).

a Interpretations of the upper-plate fluid conduits and seismic reflectivity at the plate boundary. The background color is the Vp model in Fig. 2a. Iso-velocity contour of 7.5 km/s is used as a proxy for the slab Moho. b Frequency plot of tremors15,16 and very low-frequency earthquakes (VLFE)17. The numbers of events within 5-km distance from the seismic profile are projected in the figure. Colors of each curve corresponds to different tremor/VLFE episodes shown in panel c. c Map-view distribution of each episode of tremors and VLFEs (colored dots). Solid lines indicate the location of the HYU01 seismic profile.

The variations in the Vp and reflectivity along the plate boundary agree well with the distribution of tremors and very low-frequency earthquakes (VLFEs) in the region. For example, while tremors and VLFEs actively occur around the subducting Kyushu-Palau Ridge, where the plate boundary is reflective and exhibits low seismic velocities (Vp of 3.5–4.5 km/s), their activities disappear in the western part of the profile (0–40 km distance in Fig. 5b), corresponding to the increase in the Vp value (> 4.5 km/s) and the waning seismic reflectivity around the plate boundary (0–40 km distance in Fig. 5a). Instead, the western part of the plate boundary is dominated by repeating regular earthquakes38 (Figs. 1 and 5b). This indicates a rapid transition in the slip pattern along the plate boundary that corresponds well with the variations in the physical properties around the plate boundary faults. As the Vp value measured in this study at the transitional area (~4.5 km/s) is consistent with that in previous studies on other subduction zones39,40, the value of 4.5 km/s is considered a common threshold to characterize the slip behavior along the megathrust faults.

Several episodes of tremors and VLFEs are observed within a localized area of an ~30-km range15,16,17 (“Isolated” tremor/VLFE episodes at 40–70 km distance in Fig. 5) and, similar to the western limit of the tremor zone, the eastern limits of these episodes coincide well with the nonreflective portion of the plate boundary (72–80 km distance in Fig. 2c). On the other hand, this nonreflective portion sometimes hosts tremor activities (“Pervasive” tremor episodes in Fig. 5). Although factors controlling this variation in tremor occurrence are unclear, one possibility inferred from our structural observation is the role of the upper-plate conduit. The Vp model indicates that one of the low-velocity anomalies exits in the upper plate just above the non-reflective portion (at 72–80 km distance in Figs. 2b and 5a) and thus may suggest that the fluid pressure at the plate boundary is modulated by the partial leakage of fluids into the upper-plate conduits. In general, the upper-plate drainage system is known to play such a role4. Although this hypothesis is proposed for the tremor activity in deep subduction zones (i.e., depths > 30 km), it may also be true for shallow subduction zones (i.e., depths < 15 km).

The seismic data were acquired through two research cruises conducted in August 2020. During the KM20-05 cruise by the R/V Kaimei of Japan Agency for Marine-Earth Science and Technology (JAMSTEC), we deployed 50 OBSs along the HYU01 profile at 2-km intervals. Each OBS contained three‐component geophone and hydrophone sensors (with a natural frequency of 4.5 Hz; the sampling rate for recording was 200 Hz). The OBSs recorded the acoustic signals from 596 airgun shots triggered every 200 m along the 120-km-long seismic line. All the OBSs were recovered by the R/V Kairei during the KR20-10 cruise, providing the seismic records of the airgun signals. We also collected the MCS reflection data on the same line as the OBS refraction survey, using a 5.6-km-long 444-channel streamer cable onboard the Kaimei; the shooting interval for the reflection survey was 50 m. For both the surveys, the Kaimei’s tuned airgun arrays (with a total air volume of 10,600 cubic inch) were used as the source.

After the recovery, the OBS data were preprocessed for further analyses. First, we located the OBS positions at the seafloor, using the same method applied in our previous study41. Then, we determined the first breaks in the refraction waves manually, using the vertical component of the geophone data. The picking error ranging from 30 ms to 60 ms were assigned to each pick, depending on the offset and the noise level of the data. A total of 27,902 first arrivals were used for a traveltime tomography analysis42 to derive a preferable Vp model. The initial Vp model for the first-arrival traveltime tomography (FAT) was created by referring to the seismic refraction studies close to the survey line considered in this study25 (Supplementary Fig. 9a). The final Vp model developed using FAT that explained the traveltime data within the error range (root mean square [RMS] traveltime residual of 56 ms) was obtained after 15 iterations of inversion (Supplementary Fig. 9b).

To improve the resolution of the Vp model, we further applied the FWI analysis to the OBS records. For the analysis, we used TOY2DAC, an acoustic 2D and frequency-domain FWI code developed by SEISCOPE26,43. For the modeling, the Vp structure of 100 × 23 km was gridded every 40 m. The result of FAT was used in the starting model. We applied the pre-processing to the vertical geophone data, such as muting before first arrival and time-damping with different time windows at each step of inversion. To stabilize the calculation, the inversion was initiated from a low frequency band of 2.5–3.5 Hz at the first step and then sequentially increased up to a higher frequency band of 2.5–7.5 Hz at the final step. We sampled the frequencies every 0.25 Hz within each frequency band. Other inversion parameters are summarized in Supplementary Table 1. The OBS records contained significant attenuations in the refraction phases, with shadow zones at some locations. To better reproduce these features in the OBS records, we introduced Vp-dependent attenuation factors (Q values) in the calculation. We used the Q values of 10,000, 50, 90, 350, 500, and 850 for the regions having Vp values < 1.55, 1.55–3.5, 3.5–4.5, 4.5–6.5, 6.5–7.5, and > 7.5 km/s, respectively44,45. The initial density model was converted from the Vp model, using the Gardner relation46. Both the Q and density values remained unchanged during inversion, and only the Vp values were updated. The overall waveform fits (Supplementary Fig. 1–8), reduction in the cost function (Supplementary Table 1), and consistency with the MCS reflection data were carefully examined to evaluate the reliability of the final model. The checkerboard resolution tests demonstrated that the FWI analyses was capable for resolving the velocity anomalies with 1.5 km (horizontally) x 1 km (vertically) down to the depth of ~ 15 km (Supplementary Fig. 10).

The primary processing of the MCS reflection data was carried out by DownUnder GeoSolutions (DUG). The data were processed following a standard prestack flow that included trace editing, tidal static correction, and low-cut filtering (of 1–3 Hz). The processing also included swell-noise attenuation, deghosting, source signature deconvolution, 2D surface-related multiple elimination, least-squares adaptive subtraction, and high-resolution radon demultiple for noise reduction. Velocity analyses and prestack Kirchhoff depth migration were performed to produce the final seismic reflection image (Fig. 2c).

The magnetic anomaly data were downloaded from the Marine Trackline Geophysical Database of the National Centers for Environmental Information (https://www.ngdc.noaa.gov/mgg/gdas/). The data were gridded at 1arc-min (bin size), and then were transformed to the vertical component of the magnetic field produced by the same source magnetized in the vertical direction [reduced-to-pole (RTP) anomaly]. The paleo latitude of the study area was assumed to be 26° N, while considering the northward convergence direction of the Philippine Sea Plate.

The OBS and MCS data that support the conclusions are available through the JAMSTEC Seismic Survey Database (http://www.jamstec.go.jp/obsmcs_db/e/).

FWI code used in this study is available from SEISCOPE Consortium (https://seiscope2.osug.fr). GMT code is available at the website of University of Hawaii (http://www.soest.hawaii.edu/gmt/).

Arnulf, A. F. et al. Physical conditions and frictional properties in the source region of a slow-slip event. Nat. Geosci. 14, 334–340 (2021).

Article  CAS  ADS  Google Scholar 

Sun, T., Saffer, D. & Ellis, S. Mechanical and hydrological effects of seamount subduction on megathrust stress and slip. Nat. Geosci. 13, 249–255 (2020).

Article  CAS  ADS  Google Scholar 

Saffer, D. M. & Tobin, H. J. Hydrogeology and mechanics of subduction zone forearcs: Fluid flow and pore pressure. Annu. Rev. Earth Planet. Sci. 39, 157–186 (2011).

Article  CAS  ADS  Google Scholar 

Nakajima, J. & Hasegawa, A. Tremor activity inhibited by well-drained conditions above a megathrust. Nat. Commun. 7, 13863 (2016).

Article  CAS  PubMed  PubMed Central  ADS  Google Scholar 

Tonegawa, T., Takemura, S., Yabe, S. & Yomogida, K. Fluid migration before and during slow earthquakes in the shallow Nankai subduction zone. J. Geophys. Res.: Solid Earth 127, e2021JB023583 (2022).

Shaddox, H. R. & Schwartz, S. Y. Subducted seamount diverts shallow slow slip to the forearc of the northern Hikurangi subduction zone, New Zealand. Geology 47, 415–418 (2019).

Article  CAS  ADS  Google Scholar 

Kopf, A., Klaeschen, D. & Mascle, J. Extreme efficiency of mud volcanism in dewatering accretionary prisms. Earth Planet. Sci. Lett. 189, 295–313 (2001).

Article  CAS  ADS  Google Scholar 

Mochizuki, K., Yamada, T., Shinohara, M., Yamanaka, Y. & Kanazawa, T. Weak interplate coupling by seamounts and repeating M ~7 earthquakes. Science 321, 1194–1197 (2008).

Article  CAS  PubMed  ADS  Google Scholar 

Wang, K. & Bilek, S. L. Invited review paper: Fault creep caused by subduction of rough seafloor relief. Tectonophysics 610, 1–24 (2014).

Chesley, C., Naif, S., Key, K. & Bassett, D. Fluid-rich subducting topography generates anomalous forearc porosity. Nature 595, 255–260 (2021).

Article  CAS  PubMed  ADS  Google Scholar 

Watts, A. B., Koppers, A. A. P. & Robinson, D. P. Seamount subduction and earthquakes. Oceanography 23, 166–173 (2010).

Yagi, Y., Kikuchi, M., Yoshida, S. & Yamanaka, Y. Source process of the Hyuganada earthquake of April 1, 1968 (MJMA7.5), and its relationship to the subsequent seismicity. Zisin 51, 139–148 (1998).

Yagi, Y., Kikuchi, M. & Sagiya, T. Co-seismic slip, post-seismic slip, and aftershocks associated with two large earthquakes in 1996 in Hyuga-nada, Japan. Earth Planets Space 53, 793–803 (2001).

Takagi, R., Uchida, N. & Obara, K. Along‐strike variation and migration of long‐term slow slip events in the western Nankai subduction zone, Japan. J. Geophys. Res.: Solid Earth 124, 3853–3880 (2019).

Yamashita, Y. et al. Migrating tremor off southern Kyushu as evidence for slow slip of a shallow subduction interface. Science 348, 676–679 (2015).

Article  CAS  PubMed  ADS  Google Scholar 

Yamashita, Y., Shinohara, M. & Yamada, T. Shallow tectonic tremor activities in Hyuga‑nada, Nankai subduction zone, based on long‑term broadband ocean bottom seismic observations. Earth Planets Space 73, 196 (2021).

Tonegawa, T. et al. Spatial relationship between shallow very low frequency earthquakes and the subducted Kyushu-Palau Ridge in the Hyuga-nada region of the Nankai subduction zone. Geophys. J. Int. 222, 1542–1554 (2020).

Yokota, Y. & Ishikawa, T. Shallow slow slip events along the Nankai Trough detected by GNSS-A. Sci. Adv. 6, eaay5786 (2020).

Article  PubMed  PubMed Central  ADS  Google Scholar 

Cruz-Atienza, V. M., Villafuerte, C. & Bhat, H. S. Rapid tremor migration and pore-pressure waves in subduction zones. Nat. Commun. 9, 2900 (2018).

Article  PubMed  PubMed Central  ADS  Google Scholar 

Ujiie, Y. Mud diapirs observed in two piston cores from the landward slope of the northern Ryukyu Trench, northwestern Pacific Ocean. Mar. Geol. 163, 149–167 (2000).

Article  CAS  ADS  Google Scholar 

Nakayama, N., Ashi, J., Tsunogai, U., Gamo, T. & Takahashi, M. Sources of pore water in a Tanegashima mud volcano inferred from chemical and stable isotropic studies. Geochem. J. 44, 561–569 (2010).

Article  CAS  ADS  Google Scholar 

Hoshino, T. et al. Atribacteria from the subseafloor sedimentary biosphere disperse to the hydrosphere through submarine mud volcanoes. Front. Microbiol. 8, 1135 (2017).

Article  PubMed  PubMed Central  Google Scholar 

Nishizawa, A., Kaneda, K. & Oikawa, M. Seismic structure of the northern end of the Ryukyu Trench subduction zone, southeast of Kyushu, Japan. Earth Planets Space 61, 37–40 (2009).

Yamamoto, Y. et al. Imaging of the subducted Kyushu-Palau Ridge in the Hyuga-nada region, western Nankai Trough subduction zone. Tectonophysics 589, 90–102 (2013).

Nakanishi, A. et al Three-dimensional plate geometry and P-wave velocity models of the subduction zone in SW Japan: Implications for seismogenesis, in Geology and Tectonics of Subduction Zones: A Tribute to Gaku Kimura, 69–86, eds Byrne, T., Underwood, M. B., Fisher, D., McNeill, L., Saffer, D., Ujiie, K. & Yamaguchi, A. Geological Society of America Special Paper 534 (2018).

Górszczyk, A., Operto, S. & Malinowski, M. Toward a robust workflow for deep crustal imaging by FWI of OBS data: The eastern Nankai Trough revisited. J. Geophys. Res.: Solid Earth 122, 4601–4630 (2017).

Dominguez, S., Lallemand, S. E., Malavieille, J. & von Huene, R. Upper plate deformation associated with seamount subduction. Tectonophysics 293, 207–224 (1998).

Miyazaki, S. & Heki, K. Crustal velocity field of southwest Japan: Subduction and arc-arc collision. J. Geophys. Res. 106, 4305–4326 (2001).

Shiraishi, K., Yamada, Y. & Nibe, T. Thermogenic petroleum potential of the Nankai subduction zone offshore SW Japan. J. Pet. Geol. 42, 417–434 (2019).

Ma, B. et al. Megathrust reflectivity reveals the updip limit of the 2014 Iquique earthquake rupture. Nat. Commun. 13, 3969 (2022).

Article  CAS  PubMed  PubMed Central  ADS  Google Scholar 

Arai, R. et al. Structure of the tsunamigenic plate boundary and low-frequency earthquakes in the southern Ryukyu Trench. Nat. Commun. 7, 12255 (2016).

Article  CAS  PubMed  PubMed Central  ADS  Google Scholar 

Johnson, K. et al. Intra-oceanic submarine arc evolution recorded in an ~1-km-thick rear-arc succession of distal volcaniclastic lobe deposits. Geosphere 17, 957–980 (2021).

Barnes, P. M. et al. Slow slip source characterized by lithological and geometric heterogeneity. Sci. Adv. 6, eaay3314 (2020).

Article  CAS  PubMed  PubMed Central  ADS  Google Scholar 

Kopf, A. J. Significance of mud volcanism. Rev. Geophys. 40, 1005 (2002).

Hensen, C., Wallman, K., Schmidt, M., Ranero, C. R. & Suess, E. Fluid expulsion related to mud extrusion off Costa Rica—A window to the subducting slab. Geology 32, 201–204 (2004).

Article  CAS  ADS  Google Scholar 

Nishio, Y. et al. Origins of lithium in submarine mud volcano fluid in the Nankai accretionary wedge. Earth Planet. Sci. Lett. 414, 144–155 (2015).

Article  CAS  ADS  Google Scholar 

Mitsutome, Y. et al. Estimation of the depth of origin of fluids using noble gases in the surface sediments of submarine mud volcanoes off Tanegashima Island. Sci. Rep. 13, 5051 (2023).

Article  CAS  PubMed  PubMed Central  ADS  Google Scholar 

Igarashi, T. Catalog of small repeating earthquakes for the Japanese Islands. Earth Planets Space 72, 73 (2020).

Sallarès, V. & Ranero, C. R. Upper-plate rigidity determines depth-varying rupture behaviour of megathrust earthquakes. Nature 576, 96–101 (2019).

Article  PubMed  ADS  Google Scholar 

Bassett, D. et al. Crustal structure of the Hikurangi margin from SHIRE seismic data and the relationship between forearc structure and shallow megathrust slip behavior. Geophys. Res. Lett. 49, e2021GL096960 (2022).

Arai et al. Three‐dimensional P wave velocity structure of the northern Hikurangi margin from the NZ3D experiment: Evidence for fault‐bound anisotropy. J. Geophys. Res.: Solid Earth 125, e2020JB020433 (2020).

Fujie, G. et al. Systematic changes in the incoming plate structure at the Kuril trench. Geophys. Res. Lett. 40, 88–93 (2013).

Métivier, L. & Brossier, R. The SEISCOPE optimization toolbox: A large-scale nonlinear optimization library based on reverse communication. Geophysics 81, F1–F15 (2016).

Hino, R. et al. Qp structure of the accretionary wedge in the Kumano Basin, Nankai Trough, Japan, revealed by long-offset walk-away VSP. Earth Planets Space. 67, (2015).

Furumura, T., Hayakawa, T., Nakamura, M., Koketsu, K. & Baba, T. Development of long-period ground motions from the Nankai Trough, Japan, earthquakes: Observations and computer simulation of the 1944 Tonankai (Mw 8.1) and the 2004 SE off-Kii Peninsula (Mw 7.4) earthquakes. Pure Appl. Geophys. 165, 585–607 (2008).

Gardner, G. H. F., Gardner, L. W. & Gregory, A. R. Formation velocity and density—The diagnostic basics for stratigraphic traps. Geophysics 39, 770–780 (1974).

Sagiya, T. & Thatcher, W. Coseismic slip resolution along a plate boundary megathrust: The Nankai Trough, southwest Japan. J. Geophys. Res. 104, 1111–1129 (1999).

Nishimura, T. Short-term slow slip events along the Ryukyu Trench, southwestern Japan, observed by continuous GNSS. Prog. Earth Planet. Sci. 1, 22 (2014).

Okino, K. Magnetic anomalies in the Philippine Sea: Implications for regional tectonics. J. Geogr. 124, 729–747 (2015).

Wessel, P. & Smith, W. H. F. New improved version of the Generic Mapping Tools released. Eos. Trans. Am. Geophys. Union 79, 579 (1998).

We thank the captains, ship crews, and onboard technicians of the R/V Kaimei and Kairei for their dedicated efforts in data acquisition during the KM20-05 and KR20-10 cruises. Discussions with A. Gorszczyk, T. Tonegawa, K. Shiraishi, M. Otsubo, and M. Furuichi were helpful in improving the construction of the Vp model and interpreting its results. Additionally, we are grateful to R. Miura and K. Takizawa for their assistance in MCS data processing. This study was supported by the Japan Society for the Promotion of Science (JSPS) KAKENHI program (Grant numbers: 16H06475 (S. M.), 19H04629, 21H05202, and 22K03789 (R.A.)).

Research Institute for Marine Geodynamics, Japan Agency for Marine-Earth Science and Technology, 3173-25 Showa-machi, Kanazawa-ku, Yokohama, Kanagawa, 236-0001, Japan

Ryuta Arai, Seiichi Miura, Yasuyuki Nakamura, Gou Fujie, Shuichi Kodaira & Yuka Kaiho

Earthquake Research Institute, The University of Tokyo, 1-1-1 Yayoi Bunkyo-ku, Tokyo, 113-0032, Japan

Kimihiro Mochizuki, Rie Nakata & Masataka Kinoshita

Lawrence Berkeley National Laboratory, Berkeley, USA

Faculty of Science and Technology, Kochi University, Akebonocho 2‐5‐1, Kochi, 780‐8520, Japan

Kochi Institute for Core Sample Research, Japan Agency for Marine-Earth Science and Technology, 200 Monobe Otsu Nankoku, Kochi, 783-8502, Japan

Atmosphere and Ocean Research Institute, The University of Tokyo, 5-1-5 Kashiwanoha Kashiwa, Chiba, 277-8564, Japan

You can also search for this author in PubMed  Google Scholar

You can also search for this author in PubMed  Google Scholar

You can also search for this author in PubMed  Google Scholar

You can also search for this author in PubMed  Google Scholar

You can also search for this author in PubMed  Google Scholar

You can also search for this author in PubMed  Google Scholar

You can also search for this author in PubMed  Google Scholar

You can also search for this author in PubMed  Google Scholar

You can also search for this author in PubMed  Google Scholar

You can also search for this author in PubMed  Google Scholar

You can also search for this author in PubMed  Google Scholar

You can also search for this author in PubMed  Google Scholar

R.A. analyzed the OBS and MCS reflection data, interpreted the results, and drafted the paper through the discussion with all the co-authors. S.M. and Y.N. organized the research cruises of KM20-05 and KR20-10 and served as chief scientists. G.F., S.K., and K.M. supervised the data acquisition and processing. Y.K. participated in the data-acquisition cruise. R.N., M.K., Y. Has., and Y. Ham. led the acquisition and processing of the MCS reflection data. K.O. analyzed the seafloor magnetic data. All the authors contributed to discussions of the study results.

The authors declare no competing interests.

Nature Communications thanks the anonymous, reviewers for their contribution to the peer review of this work. A peer review file is available.

Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.

Arai, R., Miura, S., Nakamura, Y. et al. Upper-plate conduits linked to plate boundary that hosts slow earthquakes. Nat Commun 14, 5101 (2023). https://doi.org/10.1038/s41467-023-40762-4

DOI: https://doi.org/10.1038/s41467-023-40762-4

Anyone you share the following link with will be able to read this content:

Sorry, a shareable link is not currently available for this article.

Provided by the Springer Nature SharedIt content-sharing initiative

By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.

Nature Communications (Nat Commun) ISSN 2041-1723 (online)

Upper-plate conduits linked to plate boundary that hosts slow earthquakes | Nature Communications

Geophone 0.75hz Sign up for the Nature Briefing newsletter — what matters in science, free to your inbox daily.