Earthquake swarms are ubiquitous in volcanic systems, being manifestations of underlying nontectonic processes such as magma intrusions or volatile fluid transport. The Long Valley caldera, California, is one such setting where episodic earthquake swarms and persistent uplift suggest the presence of active magmatism. We quantify the long-term spatial and temporal characteristics of seismicity in the region using cluster analysis on a 25-year high-resolution earthquake catalog derived using leading-edge deep-learning algorithms. Our results show that earthquake swarms beneath the caldera exhibit enlarged families with statistically significant tendency for upward migration patterns. The ascending swarms tend to nucleate at the base of the seismogenic zone with a spatial footprint that is laterally constrained by the southern rim of the caldera. We suggest that these swarms are driven by the transport of volatile-rich fluids released from deep volcanic processes. The observations highlight the potential for extreme spatial segmentation of earthquake triggering processes in magmatic systems.
more »
« less
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.
more »
« less
- Award ID(s):
- 1848166
- PAR ID:
- 10497451
- Publisher / Repository:
- Science
- Date Published:
- Journal Name:
- Science Advances
- Volume:
- 9
- Issue:
- 42
- ISSN:
- 2375-2548
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Abstract All instrumented basaltic caldera collapses have generated Mw > 5 very long period earthquakes. However, previous studies of source dynamics have been limited to lumped models treating the caldera block as rigid, leaving open questions related to how ruptures initiate and propagate around the ring fault, and the seismic expressions of those dynamics. We present the first 3D numerical model capturing the nucleation and propagation of ring fault rupture, the mechanical coupling to the underlying viscoelastic magma, and the associated seismic wavefield. We demonstrate that seismic radiation, neglected in previous models, acts as a damping mechanism reducing coseismic slip by up to half, with effects most pronounced for large magma chamber volume/ring fault radius or highly compliant crust/compressible magma. Viscosity of basaltic magma has negligible effect on collapse dynamics. In contrast, viscosity of silicic magma significantly reduces ring fault slip. We use the model to simulate the 2018 Kı̄lauea caldera collapse. Three stages of collapse, characterized by ring fault rupture initiation and propagation, deceleration of the downward‐moving caldera block and magma column, and post‐collapse resonant oscillations, in addition to chamber pressurization, are identified in simulated and observed (unfiltered) near‐field seismograms. A detailed comparison of simulated and observed displacement waveforms corresponding to collapse earthquakes with hypocenters at various azimuths of the ring fault reveals a complex nucleation phase for earthquakes initiated on the northwest. Our numerical simulation framework will enhance future efforts to reconcile seismic and geodetic observations of caldera collapse with conceptual models of ring fault and magma chamber dynamics.more » « less
-
Abstract Inflationary deformation and very long period (VLP) earthquakes frequently accompany basaltic caldera collapses, yet current interpretations do not reflect physically consistent mechanisms. We present a lumped parameter model accounting for caldera block/magma momentum change, magma chamber pressurization, and ring fault (assumed vertical) shear stress drop. Pressurization of the underlying magma chamber is represented by a tri‐axial expansion source, and the combined caldera block/magma momentum change by a vertical single force. The model is applied to Kīlauea 2018 caldera collapse events, accurately predicting near field static/dynamic ground motions. In addition to the tri‐axial expansion source, the single force contributes significantly to the VLP waveforms. For an average collapse event with fully developed ring fault, Bayesian inversion constrains ring fault stress drop to ∼0.4 MPa and the pressure increase to ∼1.9 MPa. That the predictions fit both geodetic and seismic observations confirms that the model captures the dominant caldera collapse mechanisms.more » « less
-
Abstract The 2022 eruption at Mauna Loa, Hawai'i, marked the first extrusive activity from the volcano after 38 years of quiescence. The eruption was preceded by several years of seismic unrest in the vicinity of the volcano's summit. Characterizing the structure and dynamics of seismogenic features within Mauna Loa during this pre‐eruptive interval may provide insights into how pre‐ and co‐eruptive processes manifest seismically at the volcano. In particular, the extent to which seismicity may be used to forecast the location and timing of future eruptions is unclear. To address these questions, we construct a catalog of relocated seismicity on Mauna Loa spanning 2011–2023. Our earthquake locations image complex, sub‐kilometer‐scale seismogenic structures in the caldera and southwest rift zone. We additionally identify a set of streaks of seismicity in the volcano's northwest flank that are radially oriented about the summit. Using a rate‐and‐state friction model for earthquake occurrences, we demonstrate that the seismicity rate in this region can be modeled as a function of the stressing history caused by magma accumulation beneath the summit. Finally, we observe a mid‐2019 step change in the seismicity rate in the Ka'oiki region that may have altered the stress state of the northeast rift zone in the three years before the eruption. Our observations provide a framework for interpreting future seismic unrest at Mauna Loa.more » « less
-
Abstract During the 3 month long eruption of Kı̄lauea volcano, Hawaii in 2018, the pre‐existing summit caldera collapsed in over 60 quasi‐periodic failure events. The last 40 of these events, which generated Mw > 5 very long period (VLP) earthquakes, had inter‐event times between 0.8 and 2.2 days. These failure events offer a unique data set for testing methods for predicting earthquake recurrence based on locally recorded GPS, tilt, and seismicity data. In this work, we train a deep learning graph neural network (GNN) to predict the time‐to‐failure of the caldera collapse events using only a fraction of the data recorded at the start of each cycle. We find that the GNN generalizes to unseen data and can predict the time‐to‐failure to within a few hours using only 0.5 days of data, substantially improving upon a null model based only on inter‐event statistics. Predictions improve with increasing input data length, and are most accurate when using high‐SNR tilt‐meter data. Applying the trained GNN to synthetic data with different magma‐chamber pressure decay times predicts failure at a nearly constant stress threshold, revealing that the GNN is sensing the underling physics of caldera collapse. These findings demonstrate the predictability of caldera collapse sequences under well monitored conditions, and highlight the potential of machine learning methods for forecasting real world catastrophic events with limited training data.more » « less