Jiaxuan Li

Continuous geodetic measurements near volcanic systems can image magma transport dynamics, yet resolving dike intrusions with high spatiotemporal resolution remains challenging. We introduce fiber-optic geodesy, leveraging low-frequency distributed acoustic sensing (LFDAS) recordings along a telecommunication fiber-optic cable, to track dike intrusions near Grindavík, Iceland, on a minute timescale. LFDAS reveals distinct strain responses from nine intrusive events, six resulting in fissure eruptions. Geodetic inversion of LFDAS strain reveals detailed magmatic intrusions, with inferred dike volume rate peaking systematically 15 to 22 min before the onset of each eruption. Our results demonstrate DAS’s potential for a dense strainmeter array, enabling high-resolution, nearly real-time imaging of subsurface quasi-static deformations. In active volcanic regions, LFDAS recordings can offer critical insights into magmatic evolution, eruption forecasting, and hazard assessment.

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 provide critical in-situ insights about the subsurface faulting geometry and stress state. For frequent small earthquakes (magnitude\textless 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.

Analysis of deep subduction-zone earthquakes, those at depths greater than 60 km, reveals the physical and chemical properties of a descending oceanic lithosphere at mantle depths. Over the past five decades, it has been observed that a large fraction of deep earthquakes has non-double-couple (non-DC) seismic radiation patterns. In contrast, shallow earthquakes tend to have DC radiation patterns due to mechanisms of shear faulting. These observations have been used to argue that deep earthquakes rupture differently from shallow earthquakes. Here we show that the observed global distribution of non-DC deep earthquakes could be caused by shear faulting mechanisms, but in a highly anisotropic laminated rock fabric that surrounds the deep earthquakes within subducted slabs. For intermediate-depth earthquakes (~60–300 km), we found a large shear-wave anisotropy of ~25%, possibly caused by laminated fabric or aligned melt pockets oriented parallel to the slab interface, which provides new supporting evidence for the metamorphic dehydration reactions in slabs. However, at deep-focus depths (\textgreater300 km), the putative metastable phase-change mechanism alone cannot explain the seismic anisotropy. Instead, our results and those from recent experiments suggest materials such as magnesite, or perhaps carbonatite melt, may play a role in generating deep-focus earthquakes.