Search

The accurate location of tectonic tremors helps improve understanding of their underlying physical processes. However, current location methods often do not statistically evaluate uncertainties to a satisfactory degree and do not account for potential biases due to subsurface structures not included in the model. To address these issues, we propose a novel three-step process for locating tectonic tremors. First, the measured time- and amplitude differences between station pairs are optimized to obtain station-specific relative time and amplitude measurements with uncertainty estimates. Secondly, the time– and amplitude–distance relationships in the optimized data are used to roughly estimate the propagation speed (i.e. shear wave velocity) and attenuation strength. Linear regression is applied to each event, and the resulting velocity and attenuation strength are used for quality control. Finally, the tremor location problem is formulated within a Bayesian framework where the model parameters include the source locations, local site delay/amplification factors, shear wave velocity and attenuation strength. The Markov chain Monte Carlo algorithm is used to sample the posterior probability and is augmented by a parallel tempering scheme for an efficient global search. We tested the proposed method on ocean-bottom data indicating an intense episode of tectonic tremors in Kumano-nada within the Nankai Trough subduction zone. The results show that the range of the 95 per cent confidence interval is typically <7 km horizontally and <10 km vertically. A series of experiments with different inversion settings reveals that adopting amplitude data and site correction factors help reduce random error and systematic bias, respectively. Probabilistic sampling allows us to spatially map the probability of a tremor occurring at a given location. The probability map is used to identify lineaments of tremor sources, which provides insights into structural factors that favour tremor activity. [ABSTRACT FROM AUTHOR]

The evolution history of the Sea of Japan back-arc basin remains under debate, involving the opening of sub-basins such as the Japan and Yamato Basins. Detailed knowledge of the lithospheric structure will provide the key to understanding tectonic history. This study identifies the lithosphere–asthenosphere boundary (LAB) beneath the Sea of Japan back-arc basin using S-receiver functions (S-RFs). The study area, including the Japan and Yamato Basins, has been instrumented with broadband ocean-bottom seismometers (OBSs). S-RFs from these OBSs show negative Sp phases preceding the direct S arrivals, suggesting the LAB. The S-RFs also show abnormally reduced amplitudes. For further qualitative interpretation of these findings, we conduct transdimensional Bayesian inversion for S-wave velocity models. This less-subjective Bayesian approach clarifies that the low-velocity seafloor sediments and damped deconvolution contribute to the amplitude reduction, illuminating the necessity of such considerations for similar receiver function works. Inverted velocity structures show a sharp velocity decrease at the mantle depths, which we consider the LAB. The obtained LAB depths vary among sites: ~ 45 km beneath the Japan and Yamato Basins and ~ 70 km beneath the Yamato Rise, a bathymetric high between the two basins. The thick lithosphere beneath the Yamato Rise most likely reflects its continental origin. However, the thickness is still thin compared to that of eastern Asia, suggesting lithosphere extension by rifting. Notably, the Japan and Yamato Basins show a comparable lithospheric thickness, although the crustal thickness beneath the Yamato Basin is known to be anomalously thick. This consistency in the lithospheric thickness implies that both basins undergo similar back-arc opening processes. [ABSTRACT FROM AUTHOR]

Active‐source seismic surveys have resolved the fine‐scale P‐wave velocity (Vp) of the subsurface structure in subduction forearcs. In contrast, the S‐wave velocity (Vs) structure is poorly resolved despite its usefulness in understanding rock properties. This study estimates Vp and Vs structures of the Nankai Trough forearc, by applying transdimensional inversion to high‐frequency teleseismic waveforms. As a result, a thin (∼1 km) low‐velocity zone (LVZ) is evident at ∼6 km depth beneath the sea level, which is located ~3 km seaward from the outer ridge. Based on its high Vp/Vs ratio (∼2.5) and comparison to an existing seismic reflection profile, we conclude that this LVZ reflects a high pore pressure zone at the upper portion of the underthrust sediment. We infer that this overpressured underthrust sediment hosts slow earthquake activities and that accompanied strain release helps impede coseismic rupture propagation further updip. Plain Language Summary: Many geophysical surveys have investigated the subsurface structures of shallow subduction zones by estimating the propagation speed of compressional waves emitted from artificial explosive sources. Although shear and compressional wave speeds are necessary to understand rock properties (e.g., water pressure) of the subsurface, it has been difficult to constrain the shear wave speed with a high‐resolution. In this study, we estimate both compressional and shear wave speeds by applying an advanced technique to earthquake waveforms recorded at ocean‐bottom seismometers deployed at the Nankai subduction zone, Japan. The results show a sufficiently high spatial resolution to detect a thin (~1 km) layer, which we interpret as water‐rich subducted sediment. This water‐rich zone may promote slow slips on the megathrust fault and work as a barrier against rupture propagations during large earthquakes. Key Points: We applied transdimensional inversion of teleseismic waveforms to seafloor cabled stations at the Nankai subduction zoneThe resultant high‐resolution velocity structures allow a direct comparison to active source surveysA low‐velocity zone beneath the outer ridge is evident, which is interpreted as an overpressured portion of underthrust sediment [ABSTRACT FROM AUTHOR]

Reference 1 Yamaya L, Mochizuki K, Akuhara T, Takemura S, Shinohara M, Yamada T. CMT inversion for small-to-moderate earthquakes applying to dense short-period OBS array at off Ibaraki region. 3 Result of the CMT inversions for the earthquake (Mw 4.0) that occurred at 9:06, March 16, 2011 (JST). a Centroid location of the earthquake and the OBS array. [Extracted from the article]

We conducted centroid moment tensor (CMT) inversion for small-to-moderate aftershocks in the off Ibaraki region of the 2011 off the Pacific coast of Tohoku earthquake. In this study, we used high-frequency (0.4–1.0 Hz) seismograms from a dense array of short-period ocean bottom seismometers (OBSs) and a reliable three-dimensional (3-D) seismic velocity model. Higher frequency analysis and dense OBS arrays offer CMT solutions with high spatial resolutions. Since our OBS array observed aftershocks occurring immediately following the 2011 off the Pacific coast of Tohoku earthquake, we determined 536 CMT solutions for small-to-moderate aftershocks in the off Ibaraki region (JMA-scale magnitudes of 2.5–4.0). According to our CMT solutions, characteristics of the aftershock activities in the off Ibaraki region are classified into 5 groups: (1) thrust earthquakes, which are considered as interplate earthquakes, separated by the large slip area of the 2011 Ibaraki-oki earthquake (the largest M w 7.6 aftershock of the 2011 off the Pacific coast of Tohoku earthquake) and the tectonic tremors; (2) intraslab strike-slip earthquakes located at the north of the fault area of the 2011 Ibaraki-oki earthquake; (3) intraslab normal-fault earthquakes, which suggest a tensional stress field within the subducting Pacific Plate due to the plate bending by the overriding Philippine Sea Plate; (4) various earthquake focal mechanisms above a subducting seamount, which suggest the 3-D complex fractures; and (5) normal-fault earthquakes shallower than the interplate earthquakes, which were possibly caused by the heterogeneity of the subducting seamount. [ABSTRACT FROM AUTHOR]

The Sanriku ocean‐bottom seismometer system uses an optical fiber cable to guarantee real‐time observations at the seafloor. A dark fiber connected to a Distributed Acoustic Sensing (DAS) interrogator converted the cable in an array of 19,000 seismic sensors. We use these measurements to constrain the velocity structure under a section of the cable. Our analysis relies on 24 hr of ambient seismic field recordings. We obtain a high‐resolution 2‐D shear‐wave velocity profile by inverting multimode dispersion curves extracted from frequency‐wave number analysis. We also produce a reflection image from autocorrelations of ambient seismic field, highlighting strong impedance contrasts at the interface between the sedimentary layers and the basement. In addition, earthquake wavefield analysis and modeling help to further constrain the sediment properties under the cable. Our results show for the first time that ocean‐bottom DAS can produce detailed images of the subsurface, opening new opportunities for cost‐effective ocean‐bottom imaging in the future. Plain Language Summary: Distributed Acoustic Sensing (DAS) is a relatively new measurement method that has the potential to convert existing fiber optic communication infrastructure into arrays of thousands of seismic sensors. In this research, we connected a DAS to a cable that was originally installed at the bottom of the ocean to sustain a seismic and tsunami observatory in the Sanriku Region. We show that this new type of measurement can provide reliable information to image and explore the shallow subsurface under this fiber cable. This is the first time such analysis is performed in an oceanic environment, and our methods could be readily exportable to other fiber‐optic cables that are the backbones of our modern telecommunication. Key Points: Ocean‐bottom Distributed Acoustic Sensing is used to image shallow VS structureRayleigh wave phase velocity dispersion curves are extracted from frequency‐wave number analysisReflection image is obtained from autocorrelations of ambient seismic field [ABSTRACT FROM AUTHOR]

Slow earthquakes occur at deep and shallow plate boundaries along the Nankai Trough. Deep slow earthquakes are continuously distributed along the 30–40 km depth contours of the upper surface of the subducted Philippine Sea Plate. In contrast, shallow slow earthquakes occur in limited regions: Hyuga-nada, off Cape Muroto, and southeast off the Kii Peninsula. This review provides an overview of the up-to-date seismological, geodetic, geological, and experimental results in the shallow Nankai area for a unified understanding of the spot-like occurrence of shallow slow earthquakes. Shallow slow earthquakes tend to be distributed in transitional regions between the frictionally locked and stably sliding zones on the plate boundary. Based on geological and experimental studies, the lithology of incoming sediments and their friction coefficients can be variable along the Nankai Trough. Laboratory friction experiments revealed that sediments under shallow plate boundary conditions often exhibit positive (a − b) values, while negative (a − b) is possible via several processes. Subducted seamounts create complex fracture networks and stress shadows in their surrounding areas; however, not all subducted seamounts are related to shallow slow earthquake activities. This incomplete correlation suggests that alternative factors are required to explain the spot-like distribution of shallow slow earthquakes in the Nankai subduction zone. High pore fluid pressure conditions around shallow slow earthquake zones were interpreted based on seismological structural studies. In addition, ambient noise monitoring revealed temporal changes in seismic velocity structures associated with shallow slow earthquake migrations. This result suggests a close link between pore fluid migration and shallow slow earthquake episodes. Because transient changes in pore fluid pressure can lead to various slip behaviors, the episodic migration of pore fluid around the plate boundary could promote shallow slow earthquake activity along the Nankai Trough. [ABSTRACT FROM AUTHOR]

To obtain an internal S -wave velocity structure, we conducted a passive seismic campaign with 21 1-Hz seismometers on and around the Showa-Shinzan lava dome, which emerged during the 1943–1945 eruption of Usu Volcano, Japan. Before the campaign, we calibrated seismometers and found slight phase-response differences between seismometers of less than 1–2°. After the campaign, we extracted seismic wavefield by taking cross-correlations of vertical-component ambient noise records between seismic sites. We developed a new method to measure phase velocities of the Rayleigh wave automatically by assuming layered structure and finally obtained 1-D S -wave velocity models in summit, roof and base regions. The obtained S -wave velocity right beneath the intruded lava dome is higher than that in surrounding areas by a few tens of per cent down to a few hundred metres below sea level, indicating narrow but deep existence of the root of the lava dome. The obtained S -wave velocity at depths shallower than ∼50 m inside the lava dome in the summit area was ∼1 km s−1, significantly lower than that predicted from the density of  2.3 × 103 kg m−3 estimated in previous muon-radiography studies and a conventional scaling, indicating the effect of cracking in the lava dome. [ABSTRACT FROM AUTHOR]

Seismic waves from earthquakes recorded on the seafloor are composed of complex multiple arrivals. Here, distributed acoustic sensing (DAS) observations along a cable located offshore the Sanriku Coast, Japan, show that the local earthquake wavefield is particularly rich in Scholte waves. We introduce a processing pipeline to extract these surface waves from DAS records. We then invert hundreds of dispersion curves along a section of the cable to form a shallow high‐resolution shear‐wave velocity model. Moreover, we focus on the possible generation mechanisms of Scholte waves through a series of 2D and 3D full‐wavefield numerical simulations. We show that water phase reverberations greatly contribute to the generation of Scholte waves on the ocean floor. This study demonstrates the potential of DAS to observe and better understand a poorly known marine wave phenomenon and image the offshore shallow seismic structure with an unprecedented spatial resolution. Plain Language Summary: Distributed acoustic sensing (DAS) is a measurement technique that has recently demonstrated its utility for marine geophysics. DAS offers the possibility to observe the seismic wave at a scale and an extent previously unattainable with traditional passive seismic surveys. Here, we use a linear DAS array located offshore the Sanriku Coast, Japan. We propose a processing pipeline to extract surface waves from local earthquake DAS records and obtain hundreds of measurements along the cable. These measurements are used to infer a high‐resolution model of the near‐shore subsurface. Supported by a series of numerical simulations, we suggest that acoustic reverberations in the water column greatly contribute to the generation of surface waves on the ocean floor. This study further demonstrates that DAS can be used to understand marine wave phenomena better and image offshore seismic structures. Key Points: Passive high‐resolution Vs imaging of shallow sedimentary layers with ocean‐bottom distributed acoustic sensing and low‐magnitude earthquakesGridded slant‐stack method to extract Scholte waves from local earthquake wavefieldsFull‐wavefield simulations suggest that water phase reverberations generate Scholte waves at the ocean floor [ABSTRACT FROM AUTHOR]

An earthquake with a magnitude of 6.7 occurred in the Japan Sea off Yamagata on June 18, 2019. The mainshock had a source mechanism of reverse-fault type with a compression axis of WNW–ESE direction. Since the source area is positioned in a marine area, seafloor seismic observation is indispensable for obtaining the precise distribution of the aftershocks. The source area has a water depth of less than 100 m, and fishing activity is high. It is difficult to perform aftershock observation using ordinary free-fall pop-up type ocean bottom seismometers (OBSs). We developed a simple anchored-buoy type OBS for shallow water depths and performed the seafloor observation using this. The seafloor seismic unit had three-component seismometers and a hydrophone. Two orthogonal tiltmeters and an azimuth meter monitored the attitude of the package. For seismic observation at shallow water depth, we concluded that an anchored-buoy system would have the advantage of avoiding accidents. Our anchored-buoy OBS was based on a system used in fisheries. We deployed three anchored-buoy OBSs in the source region where the water depth was approximately 80 m on July 5, 2019, and two of the OBSs were recovered on July 13, 2019. Temporary land seismic stations with a three-component seismometer were also installed. The arrival times of P- and S-waves were read from the records of the OBSs and land stations, and we located hypocenters with correction for travel time. A preliminary location was performed using absolute travel time and final hypocenters were obtained using the double-difference method. The aftershocks were distributed at a depth range of 2.5 km to 10 km and along a plane dipping to the southeast. The plane formed by the aftershocks is consistent with the focal mechanism of the mainshock. The activity region of the aftershocks was positioned in the upper part of the upper crust. Focal mechanisms were estimated using the polarity of the first arrivals. Although many aftershocks had a reverse-fault focal mechanism similar to the focal solution of the mainshock, normal-fault type and strike–slip fault type focal mechanisms were also estimated. [ABSTRACT FROM AUTHOR]