Seismic arrival-time picking on distributed acoustic sensing data using semi-supervised learning.

Distributed Acoustic Sensing (DAS) is an emerging technology for earthquake monitoring and subsurface imaging. However, its distinct characteristics, such as unknown ground coupling and high noise level, pose challenges to signal processing. Existing machine learning models optimized for conventional seismic data struggle with DAS data due to its ultra-dense spatial sampling and limited manual labels. We introduce a semi-supervised learning approach to address the phase-picking task of DAS data. We use the pre-trained PhaseNet model to generate noisy labels of P/S arrivals in DAS data and apply the Gaussian mixture model phase association (GaMMA) method to refine these noisy labels and build training datasets. We develop PhaseNet-DAS, a deep learning model designed to process 2D spatio-temporal DAS data to achieve accurate phase picking and efficient earthquake detection. Our study demonstrates a method to develop deep learning models for DAS data, unlocking the potential of integrating DAS in enhancing earthquake monitoring.

An upper-crust lid over the Long Valley magma chamber.

Geophysical characterization of calderas is fundamental in assessing their potential for future catastrophic volcanic eruptions. The mechanism behind the unrest of Long Valley Caldera in California remains highly debated, with recent periods of uplift and seismicity driven either by the release of aqueous fluids from the magma chamber or by the intrusion of magma into the upper crust. We use distributed acoustic sensing data recorded along a 100-kilometer fiber-optic cable traversing the caldera to image its subsurface structure. Our images highlight a definite separation between the shallow hydrothermal system and the large magma chamber located at ~12-kilometer depth. The combination of the geological evidence with our results shows how fluids exsolved through second boiling provide the source of the observed uplift and seismicity.

The break of earthquake asperities imaged by distributed acoustic sensing.

Rupture imaging of megathrust earthquakes with global seismic arrays revealed frequency-dependent rupture signatures1-4, but the role of high-frequency radiators remains unclear3-5. Similar observations of the more abundant crustal earthquakes could provide critical constraints but are rare without ultradense local arrays6,7. Here we use distributed acoustic sensing technology8,9 to image the high-frequency earthquake rupture radiators. By converting a 100-kilometre dark-fibre cable into a 10,000-channel seismic array, we image four high-frequency subevents for the 2021 Antelope Valley, California, moment-magnitude 6.0 earthquake. After comparing our results with long-period moment-release10,11 and dynamic rupture simulations, we suggest that the imaged subevents are due to the breaking of fault asperities-stronger spots or pins on the fault-that substantially modulate the overall rupture behaviour. An otherwise fading rupture propagation could be promoted by the breaking of fault asperities in a cascading sequence. This study highlights how we can use the extensive pre-existing fibre networks12 as high-frequency seismic antennas to systematically investigate the rupture process of regional moderate-sized earthquakes. Coupled with dynamic rupture modelling, it could improve our understanding of earthquake rupture dynamics.

Earthquake focal mechanisms with distributed acoustic sensing.

Earthquake focal mechanisms provide critical in-situ insights about the subsurface faulting geometry and stress state. For frequent small earthquakes (magnitude< 3.5), their focal mechanisms are routinely determined using first-arrival polarities picked on the vertical component of seismometers. Nevertheless, their quality is usually limited by the azimuthal coverage of the local seismic network. The emerging distributed acoustic sensing (DAS) technology, which can convert pre-existing telecommunication cables into arrays of strain/strain-rate meters, can potentially fill the azimuthal gap and enhance constraints on the nodal plane orientation through its long sensing range and dense spatial sampling. However, determining first-arrival polarities on DAS is challenging due to its single-component sensing and low signal-to-noise ratio for direct body waves. Here, we present a data-driven method that measures P-wave polarities on a DAS array based on cross-correlations between earthquake pairs. We validate the inferred polarities using the regional network catalog on two DAS arrays, deployed in California and each comprising ~ 5000 channels. We demonstrate that a joint focal mechanism inversion combining conventional and DAS polarity picks improves the accuracy and reduces the uncertainty in the focal plane orientation. Our results highlight the significant potential of integrating DAS with conventional networks for investigating high-resolution earthquake source mechanisms.

Adaptation of Range-Doppler Algorithm for Efficient Beamforming of Monostatic and Multistatic Ultrasound Signals.

Algorithmic changes that increase beamforming speed have become increasingly relevant to processing synthetic aperture (SA) ultrasound data. In particular, beamforming SA data in a spatio-temporal frequency domain using the F-k (Stolt) migration have been shown to reduce the beamforming time by up to two orders of magnitude compared with the conventional delay-and-sum (DAS) beamforming, and it has been used in applications where large amounts of raw data make real-time frame rates difficult to attain, such as multistatic SA imaging and plane-wave Doppler imaging with large ensemble lengths. However, beamforming signals in a spatio-temporal Fourier space can require loading large blocks of data at once, making it memory-intensive and less suited for parallel (i.e., multithreaded) processing. As an alternative, we propose beamforming in a range-Doppler (RD) frequency domain using the range-Doppler algorithm (RDA) that has originally been developed for SA radar (SAR) imaging. Through simulation and phantom experiments, we show that RDA achieves similar lateral resolution and contrast compared with DAS and F-k migration. At the same time, higher axial sidelobes in RDA images can be reduced via (temporal) frequency binning. Like the F-k migration, RDA significantly reduces the overall number of computations relative to DAS, and it achieves ten times lower processing time on a single CPU. Because RDA uses only a spatial Fourier transform (FT), it requires two times less memory than the F-k migration to process the simulated multistatic data and can be executed on as many as a thousand parallel threads (compared with eight parallel threads for the F-k migration), making it more suitable for implementation on modern graphics processing units (GPUs). While RDA is not as parallelizable as DAS, it is expected to hold a significant speed advantage on devices with moderate parallel processing capabilities (up to several thousand cores), such as point-of-care and low-cost ultrasound devices.