Three-dimensional topographic relief of the oceanic crust may control the occurrence of shallow very-low-frequency earthquakes in the Nankai Trough off Kumano

To explore a local relationship between geological structures and the occurrence of very-low-frequency earthquakes (VLFEs), a particular class of slow earthquakes with characteristic periods of 10–100 s, we investigated three-dimensional (3D) structural features using reprocessed 3D seismic data from the Nankai Trough off Kumano, southwestern Japan. In this region, VLFEs have been observed along the subducting Philippine Sea Plate. Although the detailed source distribution of VLFEs was estimated by means of recent land-based and offshore seismic networks, the relation with geological features is not well understood. First, we reprocessed the 3D seismic data with advanced techniques and reinterpreted the fault distribution in the sediment layer of the accretionary prism and tracked two key horizons: a décollement and the oceanic crust surface. In the accretionary prism sediments, multiple continuous reflectors of basal detachments in the underthrust sequence and conjugate faults cutting the shallow imbricated thrust sequence were identified. In contrast to the gentle variation in the décollement surface, the topographic relief of the oceanic crust was prominent, with ridges and surface displacement due to faults in the oceanic crust. Then, we compared the structural features with the VLFE source locations. Most VLFEs were located deep in the underthrust sediments where the sediments may consist of underconsolidated muds. Furthermore, a high spatial correlation was observed between the VLFE distribution and the oceanic crust topographic relief. The maximum stress direction, which was inferred from the conjugate faults in the imbricated thrust zone, was consistent with the spatial relation between the VLFE localization and the oceanic crust central ridge. Oceanic crust ridges may cause strain accumulation in the underthrust sediments on the landward sides of the ridges, and low-angle slow thrust movements might be caused using weak slip planes in the underthrust muddy sediments. That is, the topographic relief of the oceanic crust may control the occurrence of shallow VLFEs in the Nankai Trough.


Introduction
Recent observations using land-based and offshore seismic and geodetic networks have revealed spatial variations in the occurrence of slow earthquakes, which are seismic activities with intermediate type of fault slip that is transitional between the fast rupture of regular earthquakes and stable sliding along a megathrust fault interface (Obara and Kato 2016). Slow earthquakes are typically classified as slow-slip events (SSEs) with durations of days to years, very-low-frequency events (VLFEs) with periods of 10-100 s, impulsive low-frequency earthquakes (LFEs) with frequencies around several hertz, and low-frequency tremor (LFT), and their occurrence times and locations are correlated (e.g., Ito et al. 2007). Many investigations have been conducted to explore the mechanism of slow earthquakes and their connection with large earthquakes. Various types of slow earthquakes have been observed in plate subduction zones around the world: Cascadia (Rogers and Dragert 2003), Mexico (Kostoglodov et al. 2003), New Zealand (Wallace and Beavan 2010; Wallace et al. 2016), and Japan (Ito et al. 2007;Yamashita et al. 2015;Nishikawa et al. 2019). In the Nankai Trough, southwestern Japan, Ito and Obara (2006) detected VLFEs using a land-based seismic network. In the same region, Obana and Kodaira (2009) and Sugioka et al. (2012) detected LFTs and VLFEs, respectively, using long-term ocean bottom seismic observations. Araki et al. (2017) detected shallow SSEs in the Nankai trough using subseafloor borehole observatories with pressure gauges and revealed that the SSEs occurred repeatedly and that some of them were triggered by other earthquakes. Nakano et al. (2018) located the sources of VLFEs with a dense ocean bottom seismic network and revealed the spatial localization and temporal migration of VLFE occurrences in a 1-month period in April 2016. Takemura et al. (2019) also located VLFE sources from observations using the land-based seismic network of southwestern Japan in a wide area along the Nankai Trough during a long interval from 2016 to 2018. These studies suggest spatial variations in the distribution and their mechanism along the subducting Philippine Sea Plate boundary (Obara and Kato 2016).
To understand the spatial variations in slow earthquake occurrence, the effect of geological structures is an important factor to investigate in detail. On the northern Hikurangi margin, which is one of the well-studied areas, seamount subduction has been suggested to influence the mechanism of SSEs, including the effects of the stress state and fluid migration in the deformed upper plate (Todd et al. 2018;Shaddox and Schwartz 2019;Sun et al. 2020). In the Nankai Trough, Ito and Obara (2006) and Obana and Kodaira (2009) discussed the possibility of dynamic deformation by reverse faults within the accretionary prism. Sugioka et al. (2012) suggested that the décollement of the imbricated thrust zone was a favorable slip plane for VLFEs. Toh et al. (2018) investigated the effect of large subducting ridges imaged by two-dimensional (2D) seismic data in the Tonankai region (Park et al. 2003) on VLFE occurrence. Tonegawa et al. (2018) reported that VLFEs were activated in areas with sporadic distributions of low S-wave velocity bodies beneath seismic stations along the Nankai Trough. Kitajima and Saffer (2012) and Tsuji et al. (2014) noted that the underlying sedimentary sequence of the imbricated thrust zone in the accretionary prism was a potential area for VLFE occurrence based on their quantitative estimation of possible high pore pressure. Seismic data can provide high-resolution structural features, which are useful to compare the local effects on shallow VLFE activities. The spatial characteristics of VLFE occurrence, however, have not been well studied based on 3D high-resolution geological structures.
To further understand the local VLFE occurrence in the Nankai Trough accretionary prism, high-resolution information on the geological structure should be extracted from 3D seismic data. In our study, we investigated the structural features from newly reprocessed 3D seismic images from an offshore area of Kumano where VLFEs have been repeatedly observed. The 3D seismic data were acquired in 2006 (Moore et al. 2007) and reprocessed in 2016 to better understand the seismogenic zones (Shiraishi et al. 2019a). Since many VLFEs in and around the seismic survey area are well documented by Nakano et al. (2018), this is one of the best areas to investigate a local relationship between geological structures and VLFE occurrences.
In this paper, we further reprocessed the 3D seismic data in the area of the VLFE occurrences and reinterpreted the fault distribution in the shallow sediments in the accretionary prism, and we tracked two key horizons: a décollement and the oceanic crust surface. We compared the local structure with VLFE source distributions. We discussed the possible local effects of the oceanic crust topographic relief and the related fault distribution on VLFE occurrence in this area.

Geological setting
The Nankai Trough is located southwest of Japan ( Fig. 1a), where the Philippine Sea Plate is subducting to the northwest beneath the Amurian Plate at rates of 40-65 mm/year and azimuths of approximately 300°-315° clockwise from the north (Seno et al. 1993;Miyazaki and Heki 2001;Bird 2003). The Shikoku Basin is the northern margin of the subducting Philippine Sea Plate and is filled with thick layers of hemipelagic muds and turbidites. The large volume of sediments has caused the rapid evolution of the accretionary prism throughout the various accretion phases (e.g., Taira 2001;Kimura et al. 2014).
Based on many seismic surveys using 2D and 3D multichannel seismic (MCS) systems and 2D wide-angle seismic surveys using linear arrays of ocean bottom seismographs (e.g., Kodaira et al. 2000;Park et al. 2003;Moore et al. 2007;Nakanishi et al. 2008), the regional crustal structures and fault distributions in this region are generally well known. The accretionary prism can be divided into the outer prism (seaward region) and the inner prism (landward region) around the seaward edge  Nakano et al. (2018). The white squares denote the locations of ocean bottom seismic stations of the DONET. The magenta circle denotes the location of IODP Site C0006. b A vertical section along line IL25000 with vertical exaggeration (VE) of 1.5. The accretionary prism is divided into the inner prism and outer prism by a splay fault. Our study area enclosed by the red rectangle is characterized by imbricated thrusts and an underthrust succession above the oceanic crust of the Kumano Basin by an out-of-sequence thrust, or a splay fault, branching from the plate interface (Fig. 1b). The splay fault is considered one of the primary coseismic faults of tsunamigenic earthquakes (e.g., Moore et al. 2007). Sakaguchi et al. (2011) also showed evidence of coseismic slip at the updip end of the plate boundary fault and along the décollement at the toe of the accretionary prism, from core samples at the Integrated Ocean Drilling Program (IODP) Site C0006 (Fig. 1a, magenta dot). In addition to historical tsunamigenic earthquakes that have occurred repeatedly with recurrence intervals of 100-200 years (Ando 1975), slow earthquake activity has been identified by recent observations as mentioned above (e.g., Ito and Obara 2006;Obana and Kodaira 2009;Araki et al. 2017;Nakano et al. 2018;Takemura et al. 2019).
In the outer accretionary prism, where shallow VLFEs and LFTs are frequently observed, the spatial and temporal evolution of the splay fault have been documented using 3D seismic data and drilling data (Strasser et al. 2009;Kimura et al. 2011;Yamada et al. 2013). The outer accretionary prism is actively deforming with imbricated thrust faults (Fig. 1b). In contrast to the alternating succession in the overlying sequence with imbricated thrusts, the underlying sequence is less reflective, and it is considered to contain underthrust sediments of massive hemipelagic muds deposited in the Shikoku Basin (e.g., Park et al. 2010;Uunderwood et al. 2010). The possibilities of low velocity (Park et al. 2010;Kamei et al. 2012) and pore-fluid overpressure (Kitajima and Saffer 2012;Tsuji et al. 2014) in the underthrust sediments have been suggested. Within the underthrust sediments, a clear positive reflection surface marking a décollement is recognized, which continues from the deep plate boundary megathrust fault to the trench (Fig. 1b).

Reprocessed 3D MCS PSDM images
To analyze the 3D structures, we used 3D seismic data that were acquired in a southeastern offshore area off the Kii Peninsula in 2006 ( Fig. 1, black polygon) (Moore et al. 2007). Two airgun arrays with a 75 m spacing, each totaling 50.64 L (3090 cubic inches) with 13.8 MPa (2000 psi), were alternately fired at 37.5 m shot intervals. Four 4.5-km hydrophone streamer cables with 360 receiver groups at a 12.5 m spacing were towed 150 m apart. The resulting 3D dataset covered an area of approximately 62.5 km × 12 km with irregular boundaries around the trench because of acquisition difficulties due to the high-speed Kuroshio current. The initial preprocessing and a full 3D prestack depth migration (PSDM) were performed in 2006 ). The 3D bins for the final migrated data were located every 12.5 m in the inline direction (NW-SE) and every 18.75 m in the crossline direction (SW-NE).
Based on the improved results from the reprocessing of the middle part from XL 4750 to XL 6750 in 2016 (Shiraishi et al. 2019a), we additionally applied the same methods to the southern part from XL 3000 to XL 4750 in 2018. We seamlessly merged the two reprocessed volumes of the neighboring parts to one volume. The preprocessing from the original data included an optimal combination of several methods to attenuate noise and multiple reflections, broadband processing with deghosting to recover desirable frequency responses, and data regularization to mitigate nonuniform fold distribution. Using the preprocessed data, a P-wave velocity model for PSDM and the PSDM volumes of both Kirchhoff migration (e.g., Sheriff 2002) and beam migrations (Sherwood et al. 2009) were produced (see details in Shiraishi et al. 2019a). The image quality of the new PSDM results was better than that of the legacy PSDM result (Fig. 2). The images from the Kirchhoff PSDM (Fig. 2b, e) were useful for interpreting the detailed deformation structures in the accretionary prism (Shiraishi et al. 2019a, b). The beam PSDM images (Fig. 2c, f ) with better continuity and less noise than in the Kirchhoff PSDM images were useful for interpreting deeper parts (Shiraishi et al. 2019a). In this study, we focused on the southern part XL 3200-5800 of the 3D survey box where VLFEs occurred ( Fig. 1a, b, red rectangles) and the southern end of the study area was limited because of poor images in the southernmost area due to irregular boundary of data acquisition. The extracted vertical sections of the beam PSDM in the NW-SE direction in our study area are shown in Fig. 3.

Structural interpretation by horizon picking and seismic attributes
To investigate the structural features, we interpreted horizons (i.e., stratigraphic surfaces) with the reprocessed 3D PSDM images. We tracked two horizons, i.e., a décollement and the top of the oceanic crust, by manually picking their reflections on the seismic sections at every 5 or 10 lines in the NW-SE direction and smoothly interpolating them for each continuous surface. The surfaces divided by displacement due to faults were treated separately, as shown in Fig. 4.
To analyze the fault distribution in shallow sediments, we generated a semblance-based coherence cube (Marfurt et al. 1998) from the Kirchhoff PSDM volume. The coherence is a measure of similarity among surrounding traces (e.g., Sheriff 2002) such that higher values indicate more continuous reflection events and lower values indicate more discontinuities, such as layer boundaries or faults. The common depth slices extracted from the

VLFE source distributions
To compare the obtained structural features with VLFE occurrence, we used an existing VLFE source catalog obtained by Nakano et al. (2018) (Fig. 1, yellow dots) in which the centroid moment tensors (CMTs) of VLFEs that occurred in April 2016 were obtained by waveform inversion using data from a permanent ocean bottom seismic network, the Dense Oceanfloor Network system for Earthquakes and Tsunamis (DONET; Kaneda et al. 2015;Kawaguchi et al. 2015). Figure 6a shows the epicenter distribution. The VLFEs were distributed in the outer prism parallel to the trough axis. Figure 6b shows the standard errors of the source locations in two horizontal directions and the depth direction estimated by bootstrapping . The average location errors were 0.0424° (approximately 3.95 km) and 0.0383° (approximately 4.26 km) in longitude and latitude, respectively, and 1.62 km in depth.
In our study, we used only the VLFE sources located with higher accuracy and having standard errors smaller than the average values in each direction described above. The resultant sources can be separated into several clusters (Fig. 6c). We successfully removed misaligned events located outside the observation network or small events for which few waveforms were available. Figure 6d shows a comparison of the source depth distributions before and after the screening, which ensures that the locations of deeper events that are less accurate have been removed. The selected events are mostly located between depths of 5.5 km and 7.5 km. We note that the source locations of Nakano et al. (2018) are based on the 1D velocity structure of P-and S-waves in this region, and systematic errors of a few km deep or landward may exist because of the heterogeneous structure in the horizontal directions in the accretionary prism.

Results
A complex thrust fault distribution is observed in the accretionary prism (Fig. 3). In contrast to the alternating sediments in the overlying imbricated thrust sequence, the thick underlying sequence is less reflective above the oceanic crust. A clear reflector within the low-reflectivity sequence has been considered a décollement (e.g., Moore et al. 2007). In the reprocessed images, several continuous reflectors are clearly identified in the thick low-reflectivity sequence above the décollement. Some of these reflectors can be traced as detachments that are continuous from shallow thrust faults to the décollement (Fig. 3b, d).
The two interpreted horizons, the décollement and the top of the oceanic crust, generally deepen from southeast to northwest (Fig. 4). The slope of the décollement surface is mostly gentle but partly steep in the southwestern area and approximately 6 km landward from the trench (Fig. 4a). Regarding the top of the oceanic crust (Fig. 4b), topographic relief is prominent with several ridges and faults. The topographic relief is caused by faults in the oceanic crust, which are part of the widely distributed thrust and strike-slip faults in this region and have grown in the oceanic crust before and during plate subduction in regional tectonics . A remarkable ridge in the oceanic crust surface with a topographic a b c d Fig. 6 a Distribution of VLFE sources . b Standard errors of the VLFE source locations in the longitude, latitude, and depth directions (black crosses). The red circles indicate VLFE sources whose errors are lower than the mean value of the standard deviation (STD) for each direction. c The distribution of VLFEs with location errors smaller than the mean value of the STD; the red dashed circles indicate VLFE clusters. d Histograms of the depth distribution of the original (gray) and selected (red) VLFEs height of approximately 1 km and an area of approximately 3 km × 8 km is present in the central part near IL2400-2500 and XL 4400 (Fig. 4b). The topographic change in the décollement above the central ridge in the oceanic crust is small, while the topographic high associated with the steep slope of the décollement in the southwestern area corresponds to another ridge in the oceanic crust.
In the horizontal slice of the coherence cube at a depth of 4 km from mean sea level (Fig. 5a, b), the typical features of imbricated fold-and-thrust structures and conjugate strike-slip faults cutting the in-sequence thrusts are identified (Fig. 5a). The vertical displacement of the conjugate faults and related faults can be observed in the SW-NE vertical sections (Fig. 5c, d). The conjugate faults on the horizontal slice imply the principal stress direction, which can be inferred as approximately 333°. The inferred principal stress direction is approximately perpendicular to the anticline axes and oblique to the direction of Philippine Sea Plate subduction at 300°-315° (Fig. 5b).
Projecting the selected VLFE sources on a vertical section extracted from the 3D PSDM volume, most of the VLFEs are located in the low-reflectivity sequence between the shallow imbricated thrust sequence and the oceanic crust (Fig. 7a). On the plane view of the oceanic crust surface with the VLFEs, a localized distribution of VLFE sources can be recognized in the area approximately 4 km northwest of the central ridge of the oceanic crust (Fig. 7b). The conjugate faults are located in the area approximately 5 km northwest of the steep slope of the décollement (Fig. 7c), while the fault distribution does not correspond to the VLFE sources and oceanic crust topography. The spatial relation between the central ridges and VLFEs is consistent with the direction of the maximum horizontal stress inferred from the conjugate faults rather than the direction of plate subduction (Fig. 7b). The slip directions of low-angle reverse fault movements from the CMT solutions of VLFEs are also consistent with the direction of the present maximum horizontal stress (Fig. 7a, b).

Detachments in the underthrust sequence
We identified some novel continuous reflectors in the thick low-reflectivity sequence that may be part of the detachment of the thrust faults likely converging to the décollement. In terms of the low-reflectivity sequence, Park et al. (2010) and Kamei et al. (2012) noted that the seismic velocity in the underlying low-reflectivity sequence was lower than that in the overlying imbricated thrust sequence by analyzing the 3D MCS data and 2D wide-angle seismic data acquired with ocean bottom seismographs. The thick low-reflectivity sequence was inferred to represent underthrust sediments with currently forming antiformal stacking (Park et al. 2010;Tsuji et al. 2014). The underthrust sediments consist of hemipelagic muds accompanied by turbidites and volcanic ash that were deposited in the Shikoku Basin before subduction (e.g., Park et al. 2010;Uunderwood et al. 2010). The underthrust mud-dominant sediments are considered underconsolidated Park et al. 2010) and have high pore pressure (Kitajima and Saffer 2012;Tsuji et al. 2014).
The clear continuous reflectors in the underthrust sediments may be detachments of the in-sequence thrust faults that might have developed using thin weak layers within the massive underconsolidated muddy sediments. Tsuji et al. (2014) noted that discontinuous reflections above the décollement imply structures deformed by shear zones within the underthrust sequence. These continuous and discontinuous detachments might work as slip planes in addition to the major décollement in the underthrust sediments.

Stress field inferred from the conjugate faults cutting in-sequence thrusts
We identified the conjugate faults in the shallow imbricated thrust sequence and inferred that the maximum horizontal stress axis is oriented 333°. The inferred stress direction is almost consistent with the maximum stress direction of 330°, which was estimated from borehole breakout analyses at Site C0006 (Fig. 1, magenta dot) (Lin et al. 2016). This consistency in the maximum stress direction, which is oblique to the direction of plate subduction, implies that conjugate faults were generated in the present stress field within the accretionary prism. Furthermore, the conjugate faults with vertical displacements subsequently cutting the in-sequence thrust faults (Fig. 5) imply a possible change in the stress field from a reverse fault regime to a strike-slip fault regime with upward movement.
The conjugate faults are located in the area northwest of the topographic high of the décollement and their positional relationship is in the same direction as the maximum horizontal stress axis (Figs. 5d and 7c). This relationship suggests that a possible change in the stress regime was caused by the steeper slope of the décollement around the southwest area. Bilek (2011, 2016) reviewed the mechanism of seamount subduction and explained that geometric incompatibility must be overcome by the deformation of the surrounding rocks and the seamount itself with a complex network of fractures. A 3D numerical simulation study by Ruh et al. (2016) suggested that uplift movement with tectonic underpressure may generate extension in the shallow Fig. 7 Comparisons of the 3D seismic images, structural features, and VLFE source distributions. a CMT solutions within 4 km, approximately the maximum spatial error, from line IL2500 are projected on the seismic section. The vertical exaggeration (VE) of the seismic section is 1.5. Most of the selected VLFEs are located in low-reflectivity underthrust sediments on the landward side of the central ridge in the oceanic crust and are thought to consist of underconsolidated hemipelagic muds. b The oceanic crust topography and VLFE distribution with the CMT solutions. A cluster of VLFEs is located on the landward side of the central ridge in the oceanic crust. c The décollement topography and fault distribution in the shallow imbricated thrust sequence. The conjugate faults are also located on the landward side of the topographic high of the décollement. These comparisons imply that the maximum stress direction, which is inferred from the conjugate faults in the imbricated thrust zone, is consistent with the slip direction of VLFEs and the spatial relation between the VLFE distribution and the oceanic crust ridge. The consistencies suggest that the activity and distribution of VLFEs are controlled by the relief of the oceanic crust affecting the present stress field in the overlying sequence part of the rear side of an underthrusting seamount. Accordingly, conjugate faults with vertical displacements on the landward side of the steeper décollement may be generated in a manner similar to that of the stress field on the rear side of an underthrusting seamount.

3D structural features and VLFE occurrence
A spatial comparison between the structural features and the VLFE source (Fig. 7) implies a high correlation of the VLFE occurrence with the central ridge in the oceanic crust topographic relief rather than with the décollement surface. Most of the VLFEs are located in the thick underthrust sequence at depths of 5.5-7.5 km (Fig. 7a), where the major décollement and other possible detachments may work as potential slip planes in the underthrust sediments. The VLFE sources in a cluster (Fig. 6c) are localized 3-4 km northwest of the central ridge in the oceanic crust (around XL4400-4800 in Fig. 7a). Both the slip directions of VLFEs and the spatial relation between the ridges and VLFEs are consistent with the direction of the maximum horizontal stress rather than the direction of plate subduction (Fig. 7a, b). These consistencies suggest that the activity and distribution of VLFEs are governed by the subducting ridges in the oceanic crust and the present stress state in the overlying sequence. In the area above the relatively flat surface of the oceanic crust on the landward side of the central ridge, seismic activity potentially occurs due to local compressive stress accumulation caused by the ridge with a topographic height of approximately 1 km (Fig. 7a, b). The 3D numerical simulation by Ruh et al. (2016) showed that tectonic overpressure (stress accumulation) was generated on the rear side (landward side) along the flank of a subducting seamount. Indeed, earthquakes with possible stress accumulation have been observed on the landward sides of seamounts (e.g., Mochizuki et al. 2008;Geersen et al. 2015). The mechanism of SSEs and tremors in the northern Hikurangi margin is influenced by subducted seamounts (e.g., Wallace et al. 2016;Todd et al. 2018;Shaddox and Schwartz 2019). In a similar manner, the VLFEs may occur due to local stress accumulation in the underthrust sediments caused by underthrusting ridges in the oceanic crust in our study area. The topographic relief with ridges and faults in the subducting oceanic crust may affect the local VLFE distribution in addition to subducting seamounts that affect the regional seismic activity.
Consequently, we propose that the topographic relief of the oceanic crust is an important factor controlling the occurrence of VLFEs. Most of the VLFEs, which were observed as low-angle thrust movements in April 2016   (Fig. 7a), may occur along the décollement or neighboring detachments in the underthrust muddy sediments under the condition of overpressure around the landward side of the oceanic crust ridges. We note that the source locations estimated with the 1D velocity structure may contain possible systematic errors in depth or landward position because of the heterogeneous structure in the accretionary prism. To reduce errors in estimating VLFE source locations, realistic 3D velocity models of both P-and S-waves and efficient modeling methods dealing with heterogeneous models are necessary in addition to improved dense and wide observations by offshore seismic networks. To further understand the mechanism of VLFEs and other seismic activity with respect to geological structures in the wide area of the Nankai Trough, detailed investigations with dense seismic surveys are required to characterize the oceanic crust topography and overlying sediments.

Conclusions
We showed a local relationship between 3D structures and shallow VLFEs in the Nankai Trough off Kumano. We identified new structural features based on the reinterpretation of the 3D seismic data and compared them with the reliable source locations of VLFEs.
1. VLFEs with low-angle thrust movements may occur along detachments in the underthrust sediments. The topographic relief of the oceanic crust and slip planes in the overpressured sediments on the landward sides of the oceanic crust ridges may control the occurrences of the shallow VLFEs in the Nankai Trough. 2. Multiple detachments of imbricated thrusts, which may provide potential slip planes in the low-reflectivity underthrust sediments, were identified above the décollement. In contrast to the gentle topography of the décollement, the topographic relief of the oceanic crust is prominent. Newly identified conjugate faults cutting the shallow in-sequence thrusts with vertical displacements imply that the maximum horizontal stress direction is oblique to the direction of plate subduction. 3. Most of the VLFEs are located in the deep underthrust sediments, which consist of massive underconsolidated muds, possibly with high pore pressure. The localized VLFE sources are found northwest of the oceanic crust ridges, and the spatial relationship between them is consistent with the direction of the maximum stress. The underthrusting ridges in the oceanic crust may affect the stress accumulation in the underthrust sediments.