An official website of the United States government
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.
more »
« less
This content will become publicly available on February 1, 2026
Constraining Dike Opening Models With Seismic Velocity Changes Associated With the 2023–2024 Eruption Sequence on the Reykjanes Peninsula
Abstract The stress field perturbation caused by magmatic intrusions within volcanic systems induces strain in the surrounding region. This effect results in the opening and closing of microcracks in the vicinity of the intrusion, which can affect regional seismic velocities. In late November 2023, we deployed a distributed acoustic sensing interrogator to convert an existing 100‐km telecommunication fiber‐optic cable along the coast of Iceland's Reykjanes peninsula into a dense seismic array, which has run continuously. Measuring changes in surface wave moveout with ambient noise cross‐correlation, we observe up to 2% changes in Rayleigh wave phase velocity following eruptions in the peninsula's 2023–2024 sequence that are likely associated with magmatic intrusions into the eruption‐feeding dike. We apply a Bayesian inversion to compute the posterior distribution of potential dike opening models for each eruption by considering measurements for varying channel pairs and frequency bands, and assuming this velocity change is tied to volumetric strain associated with dike‐opening. Our results are in agreement with those based on geodetic measurement and provide independent constraints on the depth of the dike, demonstrating the viability of this novel inversion and new volcano monitoring directions through fiber sensing.
more »
« less
- Award ID(s):
- 1848166
- PAR ID:
- 10601299
- Publisher / Repository:
- AGU
- Date Published:
- Journal Name:
- AGU Advances
- Volume:
- 6
- Issue:
- 1
- ISSN:
- 2576-604X
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
ABSTRACT Distributed acoustic sensing (DAS) technology is an emerging field of seismic sensing that enables recording ambient noise seismic data along the entire length of a fiber-optic cable at meter-scale resolution. Such a dense spatial resolution of recordings over long distances has not been possible using traditional methods because of limited hardware resources and logistical concerns in an urban environment. The low spatial resolution of traditional passive seismic acquisition techniques has limited the accuracy of the previously generated velocity profiles in many important urban regions, including the Reno-area basin, to the top 100 m of the underlying subsurface. Applying the method of seismic interferometry to ambient noise strain rate data obtained from a dark-fiber cable allows for generating noise cross correlations, which can be used to infer shallow and deep subsurface properties and basin geometry. We gathered DAS ambient noise seismic data for this study using a 12 km portion of a dark-fiber line in Reno, Nevada. We used gathered data to generate and invert dispersion curves to estimate the near-surface shear-wave velocity structure. Comparing the generated velocity profiles with previous regional studies shows good agreement in determining the average depth to bedrock and velocity variations in the analyzed domain. A synthetic experiment is also performed to verify the proposed framework further and better understand the effect of the infrastructural cover along the cable. The results obtained from this research provide insight into the application of DAS using dark-fiber lines in subsurface characterization in urban environments. It also discusses the potential effects of the conduit that covers such permanent fiber installations on the produced inversion results.more » « less
-
null (Ed.)SUMMARY Interfaces are important part of Earth’s layering structure. Here, we developed a new model parametrization and iterative linearized inversion method that determines 1-D crustal velocity structure using surface wave dispersion, teleseismic P-wave receiver functions and Ps and PmP traveltimes. Unlike previous joint inversion methods, the new model parametrization includes interface depths and layer Vp/Vs ratios so that smoothness constraint can be conveniently applied to velocities of individual layers without affecting the velocity discontinuity across the interfaces. It also allows adding interface-related observation such as traveltimes of Ps and PmP in the joint inversion to eliminate the trade-off between interface depth and Vp/Vs ratio and therefore to reduce the uncertainties of results. Numerical tests show that the method is computationally efficient and the inversion results are robust and independent of the initial model. Application of the method to a dense linear array across the Wabash Valley Seismic Zone (WVSZ) produced a high-resolution crustal image in this seismically active region. The results show a 51–55-km-thick crust with a mid-crustal interface at 14–17 km. The crustal Vp/Vs ratio varies from 1.69 to 1.90. There are three pillow-like, ∼100 km apart high-velocity bodies sitting at the base of the crust and directly above each of them are a low-velocity anomaly in the middle crust and a high-velocity anomaly in the upper crust. They are interpreted to be produced by mantle magmatic intrusions and remelting during rifting events in the end of the Precambrian. The current diffuse seismicity in the WVSZ might be rooted in this ancient distributed rifting structure.more » « less
-
Abstract The Hawaiian‐Emperor seamount chain in the Pacific Ocean has provided fundamental insights into hotspot generated intraplate volcanism and the long‐term strength of oceanic lithosphere. However, only a few seismic experiments to determine crustal and upper mantle structure have been carried out on the Hawaiian Ridge, and no deep imaging has ever been carried out along the Emperor seamounts. Here, we present the results of an active source seismic experiment using 29 Ocean‐Bottom Seismometers (OBS) carried out along a strike profile of the seamounts in the region of Jimmu and Suiko guyots. Joint reflection and refraction tomographic inversion of the OBS data show the upper crust is highly heterogeneous withPwave velocities <4–5 km s−1, which are attributed to extrusive lavas and clastics. In contrast, the lower crust is remarkably homogeneous with velocities of 6.5–7.2 km s−1, which we attribute to oceanic crust and mafic intrusions. Moho is identified by a strongPmParrival at offsets of 20–80 km, yielding depths of 13–16 km. The underlying mantle is generally homogeneous with velocities in the range 7.9–8.0 km s−1. The crust and mantle velocity structure has been verified by gravity modeling. While top of oceanic crust prior to volcano loading is not recognized as a seismic or gravity discontinuity, flexural modeling reveals a ∼5.0–5.5 km thick preexisting oceanic crust that is overlain by a ∼8 km thick volcanic edifice. Unlike at the Hawaiian Ridge, we find no evidence of magmatic underplating.more » « less
-
Abstract Liu et al. (2022,https://doi.org/10.1029/2021GL093691) used Rayleigh waves extracted from the cross‐correlation of ambient noise recorded by two stations to monitor the seismic velocity variations associated with the 2018 Kı̄lauea eruption. However, their study ignored the fact that the tremors on the Island of Hawai'i were dominated by a source at the Kı̄lauea summit before the eruption. Close inspection of the waveforms of the station pair PAUD‐STCD shows a simple, mistakenly identified wave traveling direction in Liu et al. (2022,https://doi.org/10.1029/2021GL093691). A correct wave traveling direction agrees with the noise source model, where the dominant tremor source should be at the Kı̄lauea summit. Because of the drastic change in the tremor source after the eruption, the cross‐correlation of the tremor records may reflect predominantly changes in the source rather than in the medium properties between the two stations.more » « less
