Outer trench slope extension to frontal wedge compression in a subducting plate

Faulting in subducting plates is a critical process that changes the mechanical properties the subducting lithosphere and serves as a carrier of surface materials into mantle wedges. Two intraplate earthquake sequences located in the northern Manila subduction system were investigated in this study, which revealed distinct fault planes but a contrasting seismogeny over the northern Manila Trench. The seismic sequences analyzed in this study were of small-to-moderate events. The events were separately acquired by two ocean-bottom seismometer networks deployed on the frontal accretionary wedge in 2005 and the outer trench slope in 2006. The retrieved seismicity in the frontal wedge (in 2005) mainly included the overpressured sequence, whereas that in the approaching plate (in 2006) was aftershocks of an extensional faulting sequence. The obtained seismic velocity models and Vp/Vs ratios revealed that the overpressure was likely caused by dehydration within the shallow subduction zone. By using the near-field waveform inversion algorithm, we determined focal mechanism solutions for a few relatively large earthquakes. Data from global seismic observations were also used to conclude that stress transfer may be responsible for the seismic activity in the study area in 2005–2006. In late 2005, the plate interface in the frontal wedge area was unlocked by overpressure effect with the thrusting-dominant sequence. This event changed the stress regime across the Manila Trench and triggered the normal fault extension at the outer trench slope in mid-2006. However, the hybrid focal solution indicating reverse and strike-slip mechanisms provided in this study revealed that the plate interface had become locked again in late 2006.

A method for tracking blue whales (Balaenoptera musculus) with a widely spaced network of ocean bottom seismometers

Passive acoustic monitoring is an important tool for studying marine mammals. Ocean bottom seismometer networks provide data sets of opportunity for studying blue whales (Balaenoptera musculus) which vocalize extensively at seismic frequencies. We describe methods to localize calls and obtain tracks using the B call of northeast Pacific blue whale recorded by a large network of widely spaced ocean bottom seismometers off the coast of the Pacific Northwest. The first harmonic of the B call at ~15 Hz is detected using spectrogram cross-correlation. The seasonality of calls, inferred from a dataset of calls identified by an analyst, is used to estimate the probability that detections are true positives as a function of the strength of the detection. Because the spacing of seismometers reaches 70 km, faint detections with a significant probability of being false positives must be considered in multi-station localizations. Calls are located by maximizing a likelihood function which considers each strong detection in turn as the earliest arrival time and seeks to fit the times of detections that follow within a feasible time and distance window. An alternative procedure seeks solutions based on the detections that maximize their sum after weighting by detection strength and proximity. Both approaches lead to many spurious solutions that can mix detections from different B calls and include false detections including misidentified A calls. Tracks that are reliable can be obtained iteratively by assigning detections to localizations that are grouped in space and time, and requiring groups of at least 20 locations. Smooth paths are fit to tracks by including constraints that minimize changes in speed and direction while fitting the locations to their uncertainties or applying the double difference relocation method. The reliability of localizations for future experiments might be improved by increasing sampling rates and detecting harmonics of the B call.

Basin inversion: reactivated rift structures in the central Ligurian Sea revealed using ocean bottom seismometers

The northern margin of the Ligurian Basin shows notable seismicity at the Alpine front, including frequent magnitude 4 events. Seismicity decreases offshore towards the Basin centre and Corsica, revealing a diffuse distribution of low-magnitude earthquakes. We analyse data of the amphibious AlpArray seismic network with focus on the offshore component, the AlpArray ocean bottom seismometer (OBS) network, consisting of 24 broadband OBSs deployed for 8 months, to reveal the seismicity and depth distribution of micro-earthquakes beneath the Ligurian Sea. Two clusters occurred between ∼ 10 km to ∼ 16 km depth below the sea surface, within the lower crust and uppermost mantle. Thrust faulting focal mechanisms indicate compression and an inversion of the Ligurian Basin, which is an abandoned Oligocene–Miocene rift basin. The basin inversion is suggested to be related to the Africa–Europe plate convergence. The locations and focal mechanisms of seismicity suggest reactivation of pre-existing rift-related structures. Slightly different striking directions of presumed rift-related faults in the basin centre compared to faults further east and hence away from the rift basin may reflect the counter-clockwise rotation of the Corsica–Sardinia block. High mantle S-wave velocities and a low Vp/Vs ratio support the hypothesis of strengthening of crust and uppermost mantle during the Oligocene–Miocene rifting-related extension and thinning of continental crust.

Migration-Based Local Event-Location Workflow for Ocean-Bottom Seismometer (OBS) Records in Subduction Zones: A Practical Approach for Addressing a Large Number of Events

ABSTRACT An efficient event-location workflow is highly desired to analyze large numbers of local earthquakes recorded by ocean-bottom seismometers (OBSs) in subduction zones. The present study proposes a migration-based event-location approach for evaluating OBS records to examine local subduction-zone earthquakes. This approach can significantly reduce the amount of manual time picks compared with conventional methods. The event-location workflow was designed to detect arrival onsets of both P and S phases. Synthetic tests have shown that the proposed migration-based event-location method is robust against different types of noise, such as environmental noise and local spike noise. This workflow was then applied to real OBS data in the off-Ibaraki region at the southern end of the Japan trench. The results show that this approach is applicable to real data from subduction-zone events: It gives reasonable agreement with manual time picks for both P and S waves and reasonable error bars, and it demonstrates a clear down-dip trend of seismicity. The results also show fair agreement with event distributions from previous studies of the off-Ibaraki region. This proposed workflow can be used to examine the seismicity of local earthquakes around the subduction zone using OBSs. This approach is especially effective when the seismicity is high and/or in cases in which long-term OBS monitoring has recorded a large number of events.

The Pacific OBS Research into Convecting Asthenosphere (ORCA) Experiment

The Pacific ocean-bottom seismometer (OBS) Research into Convecting Asthenosphere (ORCA) experiment deployed two 30-station seismic arrays between 2018 and 2020—a US contribution to the international PacificArray project. The “Young ORCA” array deployed on ∼40 Ma central Pacific seafloor had a ∼68% data recovery rate, whereas the “Old ORCA” array deployed on ∼120 Ma southwest Pacific seafloor had a ∼80% recovery rate. We detail here the seismic data quality, spectral characteristics, and engineering challenges of this experiment. We provide information to assist users of this dataset, including OBS orientations and tables of daily data quality for all channels. Preliminary analysis illustrates the utility of these data for surface- and body-wave seismic imaging.

Preliminary Model of P-Wave Tomography Beneath Central Java using FMTOMO

The tectonic setting of Java Island is mainly controlled by the collision of Indo-Australian plate subducting the Eurasian plate. The high collision activity of Eurasian and Indo-Australian plates often causes megathrust earthquakes and the rise of arc magmatism that includes volcanic eruption. This study aims to determine the tectonic pattern beneath Central Java based on P-wave tomography inversion. We used the fast-marching method as ray tracing and subspace inversion to image subsurface velocity model to a depth of 150 km. The data used in this study are catalogue events data derived from a temporary seismometer network MERAMEX installed around central Java and DOMERAPI installed surround Mt. Merapi and Mt. Merbabu. We also include events collected from the International Seismological Centre. In total, we processed 563 earthquake events to illustrate velocity structures under central Java. The checker-board model shows that good resolutions can be identified at shallow depth, including offshore south Java contributed from Ocean Bottom Seismometer data. In vertical axis, good resolution models can be expected down to a depth 150 km following rich events from the Benioff zone. Current P wave model show a distinct low velocity zone under Mt Merapi that can be seen down to a depth of 40 km, suggesting a possible separated deep magma reservoir. To the south of Mt Merapi area also shows a low-velocity band that may be related with the southern mountain arc. Additionally, the northern part of Mt. Merapi displays a band of strong low-velocity anomaly to the East and West with the anomaly in the Eastern Part seems to have a deeper extension to a depth of ~50 km. We related this anomaly with Merapi Lawu Anomaly and Kendeng basin. Our results show a similar result with the previous tomography models in this region.

Application of the Stacked Refraction Convolution Section to 2D Ocean Bottom Seismometer Wide-angle Seismic Data Along the Tamayo Through Basin, Gulf of California

The stacked refraction convolution section can be used as an interpretation tool in wide-angle refraction seismic data generated by air gun shooting and recorded by Ocean Bottom Seismometers (OBS). The refraction convolution section is a full-wave extension of the Generalized Reciprocal Method (GRM), a method frequently used in shallow refraction seismic interpretation, but not applied to deep crustal-scale studies. The sum of the travel times of the waves refracted in the same interface and recorded in a pair of forward and reverse profiles, time-corrected by the reciprocal time, is an estimation close to the two-way travel times of the multichannel seismic reflection sections, but with seismic rays illuminating the interfaces upwards. The sum of seismic traces is obtained with the convolution section. Furthermore, several pairs of convolved forward-reverse refraction recordings of the same area can be stacked together to improve the signal to noise ratio. To show the applicability of the refraction convolution section in OBS deep data, we interpreted the basement structure of the Tamayo Through Basin in the southern Gulf of California, offshore Mexico. We compared the results with both, a multichannel seismic section recorded in the same profile, and the previous interpretations of the same wide-angle seismic data modeled with ray tracing and tomography methods. The basement imaged by the stacked refraction convolution section is similar in geometry to that obtained by seismic reflection processing. The stacked refraction convolution section identifies the full extent of the basement and confirms the location of a nearly constant thickness volcanic layer in the northwestern half of the basin. However, only a small area of volcanic deposits is found in the shallower parts of the southwestern margin. We also show that the convolution process can be used to estimate the occurrence of lateral variations of seismic velocities in the basement, as a further application of the GRM to deep refraction data.

Development of a Compact Broadband Ocean-Bottom Seismometer

Studies of very-low-frequency earthquakes and low-frequency tremors (slow earthquakes) in the shallow region of plate boundaries need seafloor broadband seismic observations. Because it is expected that seafloor spatially high-density monitoring requires numerous broadband sensors for slow earthquakes near trenches, we have developed a long-term compact broadband ocean-bottom seismometer (CBBOBS) by upgrading the long-term short-period ocean-bottom seismometer that has seismic sensors with a natural frequency of 1 Hz and is being mainly used for observation of microearthquakes. Because many long-term ocean-bottom seismometers with short-period sensors are available, we can increase the number of broadband seafloor sensors at a low cost. A short-period seismometer is exchanged for a compact broadband seismometer with a period of 20 or 120 s. Because the ocean-bottom seismometers are installed by free fall, we have no attitude control during an installation. Therefore, we have developed a new leveling system for compact broadband seismic sensors. This new leveling system keeps the same dimensions as the conventional leveling system for 1 Hz seismometers so that the broadband seismic sensor can be installed conveniently. Tolerance for leveling is less than 1°. A tilt of up to 20° is allowed for the leveling operation. A microprocessor controls the leveling procedure. Some of the newly developed ocean-bottom seismometers were deployed in the western Nankai trough, where slow earthquakes frequently occur. The data from the ocean-bottom seismometers on the seafloor were evaluated, and we confirmed that the long-term CBBOBS is suitable for observation of slow earthquakes. The developed ocean-bottom seismometer is also available for submarine volcanic observation and broadband seafloor observation to estimate deep seismic structures.

Ocean-bottom and surface seismometers reveal continuous glacial tremor and slip

Shearing along subduction zones, laboratory experiments on analogue faults, and sliding along glacier beds are all associated with aseismic and co-seismic slip. In this study, an ocean-bottom seismometer is deployed near the terminus of a Greenlandic tidewater glacier, effectively insulating the signal from the extremely noisy surface seismic wavefield. Continuous, tide-modulated tremor related to ice speed is recorded at the bed of the glacier. When noise interference (for example, due to strong winds) is low, the tremor is also confirmed via analysis of seismic waveforms from surface stations. The signal resembles the tectonic tremor commonly observed during slow-earthquake events in subduction zones. We propose that the glacier sliding velocity can be retrieved from the observed seismic noise. Our approach may open new opportunities for monitoring calving-front processes in one of the most difficult-to-access cryospheric environments.

Multi-parameter model building for the Qiuyue structure using four-component ocean-bottom seismometer data

Thin sand-mud-coal interbedded layers and multiples caused by shallow water pose great challenges to conventional 3D multi-channel seismic techniques used to detect the deeply buried reservoirs in the Qiuyue field. In 2017, a dense ocean-bottom seismometer (OBS) acquisition program acquired a four-component dataset in East China Sea. To delineate the deep reservoir structures in the Qiuyue field, we applied a full-waveform inversion (FWI) workflow to this dense four-component OBS dataset. After preprocessing, including receiver geometry correction, moveout correction, component rotation, and energy transformation from 3D to 2D, a preconditioned first-arrival traveltime tomography based on an improved scattering integral algorithm is applied to construct an initial P-wave velocity model. To eliminate the influence of the wavelet estimation process, a convolutional-wavefield-based objective function for the preprocessed hydrophone component is used during acoustic FWI. By inverting the waveforms associated with early arrivals, a relatively high-resolution underground P-wave velocity model is obtained, with updates at 2.0 km and 4.7 km depth. Initial S-wave velocity and density models are then constructed based on their prior relationships to the P-wave velocity, accompanied by a reciprocal source-independent elastic full-waveform inversion to refine both velocity models. Compared to a traditional workflow, guided by stacking velocity analysis or migration velocity analysis, and using only the pressure component or other single-component, the workflow presented in this study represents a good approach for inverting the four-component OBS dataset to characterize sub-seafloor velocity structures.