Search

This study presents an approach to better characterize the P‐wave and S‐wave velocity structure of the seafloor sediment layer using ocean bottom seismometers. The presence of low‐velocity seafloor sediment layers influences the observed seismic record at the seafloor over a broad frequency range, such that detailed knowledge of this sediment structure is essential to predict its effect on teleseismic records. We use the radial component of teleseismic P waves and autocorrelation functions of the radial, vertical, and pressure components of teleseismic P and S waves to obtain sediment layer models using the Markov chain Monte Carlo approach with parallel tempering. Synthetic tests show that the body waves constrain the P‐ and S‐wave impedances and travel times and the P‐ to S‐wave velocity ratio of the sediment layers. The proposed method resolves thin layers at a high resolution, including the uppermost thin (∼50 m to a few hundred meters) low S‐wave velocity layer. Real data applications at sites across the Pacific Ocean that are coincident with previous in situ studies demonstrate the effectiveness of this method in characterizing the seafloor sediment unit. The sediment models characterized by this new approach will allow us to more accurately predict and correct the effects of sediment layers in generating P‐ and S‐wave reverberations. Plain Language Summary: Ocean bottom seismometers are powerful tools for investigating the detailed structure of the oceanic mantle. However, most ocean bottom seismometers are placed above the low‐velocity water‐saturated sediment layer, such that their seismic records are inevitably affected by this seafloor sediment layer. Detailed constraints on the structure of the seafloor sediment layer allow us to accurately predict its effects on the observed seismic records. The seafloor sediment layer may consist of multiple discontinuities and vary regionally. This study proposes a new approach for characterizing the seismic properties of the seafloor sediment layer using teleseismic body waves recorded at ocean bottom seismometers. Synthetic tests show that the method can resolve thin sublayers within the sediment layer, such as a very thin low shear‐wave velocity layer and sharp seismic velocity jumps. Real data applications show that the sediment models obtained using this approach are in good agreement with those obtained from direct measurements and previous body‐wave analyses. Our proposed method will enable more detailed prediction and correction of the influence of the sediment layer on ocean bottom seismometer records, thereby providing more accurate constraints on deeper Earth structures. Key Points: The seismic properties of the seafloor sediment layer are characterized using ocean bottom seismometer teleseismic body‐wave recordsThe P‐ and S‐wave impedances and travel times and the VP/VS values of each sublayer within the sediment layer are resolvedIn situ seafloor sediment structure constraints will enable accurate correction and prediction of shallow structure effects [ABSTRACT FROM AUTHOR]

The oldest oceanic basin (160–180 Ma) in the western Pacific is the birthplace of the Pacific Plate and is thus essential for understanding the formation and evolution of the oceanic plate. However, the upper mantle structure beneath the region has not been thoroughly investigated because of the remoteness and difficulties of long‐term in situ seismic measurements at the ocean bottom. From 2018 to 2019, the Oldest‐1 experiment on the oldest seafloor was conducted as part of the international Pacific Array initiative. We present the first three‐dimensional P‐wave velocity structure down to a depth of 350 km based on the relative travel time residuals of teleseismic earthquakes recorded by 11 broadband ocean‐bottom seismometers operated during the Oldest‐1 experiment. Our result shows a fast P‐wave velocity anomaly (VP perturbation of 2%–4% faster than average) at a depth of 95–185 km beneath the northeast of the study area. This structure is interpreted as evidence of dry, viscous, and rigid materials at depths below the lithosphere. Two slow anomalies (VP perturbation of 2%–4% slower than average) are seen beneath the southwestern and eastern (the oldest seafloor >170 Ma) parts of the array site. The low‐velocity zones are found at depths of 95–305 km. The observed velocity structures can be indicative of plume activities that affected the upper mantle as the Pacific Plate migrated over hotspots from the southeast. Alternatively, the observed velocity features may provide seismic evidence for small‐scale sublithospheric convection. Key Points: One‐year ocean‐bottom geophysical investigation on the oldest Pacific provides seismic mantle structure of the regionDetailed 3‐D mantle structure implies complex evolution process of Pacific PlateOur model implies thermochemical modification of the upper mantle by plume interaction or small‐scale convection [ABSTRACT FROM AUTHOR]

We analyze seismic records collected at the oldest (170–180 Ma) Pacific seafloor using broadband dispersion array analysis. Using ambient noise and teleseismic waveforms, we measure Rayleigh‐wave phase velocities in a period range of 5–200 s that are inverted for array‐average one‐dimensional isotropic and azimuthally anisotropic shear‐wave velocity depth profiles from the crust to a depth of 300 km. The high‐velocity Lid and the low‐velocity zone are well‐resolved with velocity difference of ∼4%, whose transition occurs at depths between 80 and 100 km. The profile is compared with that obtained at the 130‐ and 140‐Ma seafloor. Accounting for the cooling effect due to the plate age difference, the low‐velocity zone of the oldest Pacific seafloor is ∼1.3% slower (∼110°C warmer) than that beneath the 140‐Ma seafloor, suggesting the occurrence of some reheating process beneath the oldest lithosphere. The azimuthal anisotropy at shallow depths (<50 km) is significantly different between the western and eastern areas of the array where the peak‐to‐peak amplitudes are estimated to be ∼2.8% and ∼1.6%, respectively. The fast direction is nearly parallel to the past seafloor spreading direction (perpendicular to the magnetic lineation) in the west, while it largely deviates in the east. The observed difference in azimuthal anisotropy may represent complicated evolution dynamics of the infant Pacific plate that involved a ridge‐ridge‐ridge triple junction. Plain Language Summary: One of the classical enigmas of plate tectonics is the cause of seafloor flattening where the seafloor depth becomes nearly constant for old ages (>∼70 Ma) while a simple cooling model of plate evolution predicts the age dependence. To solve this enigma, it is essential to resolve the subsurface structure of the oceanic plate and the asthenosphere below. As a part of an international effort, Pacific Array, that aims to cover the Pacific seafloor by arrays of ocean bottom seismometers, a team of Japanese and Korean scientists deployed an array in the oldest (170–180 Ma) Pacific seafloor. We determine isotropic and anisotropic shear‐wave velocity structures as a first step toward resolving the enigma. The results indicate that the asthenosphere beneath the oldest Pacific lithosphere is significantly slower than the prediction of the simple cooling and that the occurrence of some reheating process may be needed. Furthermore, azimuthal anisotropy (the dependence of the wave velocity on the horizontal propagation direction) in the shallow lithosphere shows that the anisotropy patterns differ between the eastern and western areas of the array, and simulation results indicate that the evolution process of the infant Pacific plate is not merely a simple triangular plate spreading. Key Points: Broadband Rayleigh‐wave dispersion array analysis has been conducted in the oldest (∼170 Ma) Pacific seafloorThe low‐velocity zone is significantly slower to indicate the occurrence of some reheating process beneath the oldest Pacific lithosphereThe spatial variation of the shallow azimuthal anisotropy indicates a complicated evolution of the infant Pacific plate [ABSTRACT FROM AUTHOR]

This study presents reflectivity images of the northeast (NE) Japan subduction zone continuous across the ocean and land. As nearly half of its forearc region is under the ocean, data from ocean bottom seismometers (OBSs) must be utilized to fully image the region by passive seismic analysis. The use of OBS data has been a challenge due to inherent characters of the ocean bottom observations: high noise level and effects of seafloor sediment. Now, decent imaging is possible in NE Japan overcoming the high level noise due to the accumulated data set of the OBSs. The low‐velocity of seafloor sediment significantly delays and amplifies S waves passing through them, and thus complicates teleseismic waveforms. We identify and correct these effects to produce coherent receiver function images throughout the amphibious subduction zone. Our images provide a potential for discussing new structural features and will help better understanding of the dynamics of the NE Japan subduction zone. Plain Language Summary: In northeast Japan, the Pacific plate is subducting beneath the Japanese arc, and delineating how the actual subduction is taking place at depths has been the target of many seismological studies. This study presents our effort of continuous imaging of subsurface interfaces across the ocean and land employing an amphibious data set. On the land side, there are dense seismic networks that have been used to reveal the seismic interfaces using converted phases of distant earthquakes. On the ocean side, the accumulation of data from temporary ocean bottom seismometers in the last decade and the recently installed ocean bottom cable network has made it possible to use a similar approach as on the land. However, the low‐velocity sediment beneath the seafloor causes systematic slower arrivals and larger amplitudes of the seismic phases of interest that result in highly disturbed images if we analyze them in the same way as the land. Here, we identify and correct the effects of the seafloor sediment to analyze the seismic interfaces across the ocean and land. As a result, we succeed in imaging the interfaces of the subducting oceanic crust throughout the region, as well as previously unidentified features under the ocean that may deserve future investigation. Key Points: The northeast Japan subduction zone is imaged using receiver functions (RFs) of an amphibious data setWe identify significant amplitude and time delay effects on RFs from seafloor sedimentCorrections of the sediment effects significantly improve images across an amphibious region [ABSTRACT FROM AUTHOR]

Strong long‐period seismic signals at periods around 25 and 18 s appear in the ambient noise cross‐correlation functions recorded by an array of ocean bottom seismometers (OBSs) in the western‐central Pacific. The signal amplitude varies from time to time, and the apparent travel times of the signals are typically smaller than those expected for the Rayleigh waves propagating along the great circle connecting the station pairs. From the cross‐correlation functions, the signal sources are located in the Vanuatu Arc. Local data analysis suggests the signals originate from two different sources possibly located at depths of ~0–1 km below the sea level beneath the active cones of Ambrym volcano. Plain Language Summary: Spatially localized sources that persistently generate long‐period seismic waves at periods longer than 10 s are reported worldwide. For some sources, the excitation mechanisms have been well‐investigated. A volcano in Japan generates such waves that travel more than thousands of kilometers. For others, however, the excitation origins are still poorly understood. These include the one observed in the Vanuatu Arc. A temporal observation network of seismometers installed on the seafloor in the western‐central Pacific observed persistent waves peaking at around two different periods. From the data analysis, we find the signals originate from an active volcano, Ambrym volcano, in the Vanuatu Arc. Besides, we find that those seismic waves are most likely originated from two different sources possibly located at shallow depths beneath the volcano. Key Points: Persistent long‐period (25 and 18 s) seismic signals appear in ambient noise cross‐correlations of OBS data in the western‐central PacificThe signals are originated from Ambrym volcano in the Vanuatu ArcSource locations of the two signals may be different, but both are possibly located at shallow depths right beneath the active cones [ABSTRACT FROM AUTHOR]

We performed waveform inversion of the P waveforms recorded by our BroadBand Ocean Bottom Seismometers (BBOBSs) deployed in the Northwestern Pacific. Consequently, the depth profile of the P velocity of the oceanic upper mantle, which has not been well resolved by previous surface wave or receiver function analyses, was revealed. We considered the azimuthal anisotropy in the lithosphere, which significantly improved the variance reduction from 34% to 44%. The resulting P model exhibited higher and lower velocities in the lithosphere and asthenosphere, respectively. The velocity contrast was found to be second/third of that observed in the previous S models; however the obtained model appeared to have some trade‐off with the VS structures in the vicinity of the source. We compared our P model with the previous S model obtained using our BBOBSs and obtained the VP/VS model. The resulting VP/VS model has two notable features. First, the lithosphere is characterized by a rapid increase in the VP/VS values with depth, which implies chemical stratification. Second, the VP/VS values in the vicinity of the lithosphere‐asthenosphere boundary (LAB) are larger than the synthetic values of any major upper mantle mineral predicted by considering the anharmonic effects, which suggests the effects of anelasticity or melt. Key Points: We inverted P waveforms recorded by our BBOBSs to obtain the depth profile of the P velocity beneath the older Pacific plateCombining our P model with a previous S model revealed two notable features in Vp/VS structuresThe obtained Vp/Vs model suggests chemical stratification in the lithosphere and the effects of anelasticity or melt at the LAB [ABSTRACT FROM AUTHOR]

We conducted broadband dispersion survey by deploying two arrays of broadband ocean bottom seismometers in the northwestern Pacific Ocean at seafloor ages of 130 and 140 Ma. By combining ambient noise and teleseismic surface wave analyses, dispersion curves of Rayleigh waves were obtained at a period range of 5–100 s and then used to invert for one‐dimensional isotropic and azimuthally anisotropic βV (VSV) profiles beneath each array. The obtained profiles show ~2% difference in isotropic βV in the low‐velocity zone (LVZ) at a depth range of 80–150 km in spite of the small difference in seafloor ages and the horizontal distance of ~1,000 km. Forward dispersion‐curve calculation for thermal models indicates that simple cooling models cannot explain the observed difference and an additional mechanism, such as sublithospheric small‐scale convection, is required. In addition, the fastest azimuths of azimuthal anisotropy in the LVZ significantly deviate from the current plate motion direction. We infer that these observations are consistent with the presence of small‐scale convection beneath the study area. As for azimuthal anisotropy in the Lid, the peak‐to‐peak intensity is 3–4% at the depth from Moho to ~40 km. The fastest direction is almost perpendicular to magnetic lineation in area A at 130 Ma and oblique to magnetic lineations in area B at 140 Ma, suggesting complex mantle flow beneath the infant Pacific Plate surrounded by three ridge axes. The intensity of azimuthal anisotropy in the LVZ is ~2%, indicating that radial anisotropy is stronger than azimuthal anisotropy therein. Key Points: Broadband ocean bottom seismic observation is conducted in the NW Pacific Ocean at two sites with seafloor ages of 130 and 140 MaSignificant seismic velocity variation in the LVZ of two areas is revealed that may indicate the presence of small‐scale convectionAzimuthal anisotropy depth profiles indicate that the fast directions do not necessarily coincide with the ancient spreading directions [ABSTRACT FROM AUTHOR]

Abstract: We have determined the three-dimensional shear wave speed structure of the upper mantle in and around the Philippine Sea region using seismograms recorded by dense land-based and long-term broadband ocean bottom seismographic stations. We used a surface wave tomography technique in which multimode phase speeds are measured and inverted for a 3-D shear wave speed structure by incorporating the effects of a finite frequency and ray bending. The new ocean bottom data provided us with improved spatial resolution (~300 km) in the Philippine Sea region. In the upper 120 km, the shear wave speed structure is well correlated with seafloor age. At depths greater than 160 km, fast anomalies of the subducting Pacific Plate are clearly defined. We also found slow speed anomalies beneath the Kyushu–Palau ridge at depths greater than 120 km. Along the Izu–Bonin(Ogasawara)–Mariana arc, we have detected three separate slow anomalies in the mantle wedge at depths shallower than 100 km beneath the back arc. Each anomaly has a width of ~500 km. Moreover, these three anomalies have a close relationship with the three groups of frontal and rear arc volcanoes having distinct Sr, Nd, and Pb isotope ratios. We suggest that each of the anomalies is a site of large-scale flow of upper mantle into the mantle wedge, and that each already contains a component from the adjacent subducting slab. [Copyright &y& Elsevier]