An official website of the United States government
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
Physics‐Based Model Reconciles Caldera Collapse Induced Static and Dynamic Ground Motion: Application to Kīlauea 2018
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
- PAR ID:
- 10368392
- Publisher / Repository:
- DOI PREFIX: 10.1029
- Date Published:
- Journal Name:
- Geophysical Research Letters
- Volume:
- 49
- Issue:
- 8
- ISSN:
- 0094-8276
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Abstract Axial Seamount is an active submarine volcano located at the intersection of the Cobb hot spot and the Juan de Fuca Ridge (45°57′N, 130°01′W). Bottom pressure recorders captured co‐eruption subsidence of 2.4–3.2 m in 1998, 2011, and 2015, and campaign‐style pressure surveys every 1–2 years have provided a long‐term time series of inter‐eruption re‐inflation. The 2015 eruption occurred shortly after the Ocean Observatories Initiative (OOI) Cabled Array came online providing real‐time seismic and deformation observations for the first time. Nooner and Chadwick (2016,https://doi.org/10.1126/science.aah4666) used the available vertical deformation data to model the 2015 eruption deformation source as a steeply dipping prolate‐spheroid, approximating a high‐melt zone or conduit beneath the eastern caldera wall. More recently, Levy et al. (2018,https://doi.org/10.1130/G39978.1) used OOI seismic data to estimate dip‐slip motion along a pair of outward‐dipping caldera ring faults. This fault motion complicates the deformation field by contributing up to several centimeters of vertical seafloor motion. In this study, fault‐induced surface deformation was calculated from the slip estimates of Levy et al. (2018,https://doi.org/10.1130/G39978.1) then removed from vertical deformation data prior to model inversions. Removing fault motion resulted in an improved model fit with a new best‐fitting deformation source located 2.11 km S64°W of the source of Nooner and Chadwick (2016,https://doi.org/10.1126/science.aah4666) with similar geometry. This result shows that ring fault motion can have a significant impact on surface deformation, and future modeling efforts need to consider the contribution of fault motion when estimating the location and geometry of subsurface magma movement at Axial Seamount.more » « less
-
Fault friction is central to understanding earthquakes, yet laboratory rock mechanics experiments are restricted to, at most, meter scale. Questions thus remain as to the applicability of measured frictional properties to faulting in situ. In particular, the slip-weakening distance d c strongly influences precursory slip during earthquake nucleation, but scales with fault roughness and is challenging to extrapolate to nature. The 2018 eruption of Kīlauea volcano, Hawaii, caused 62 repeatable collapse events in which the summit caldera dropped several meters, accompanied by M W 4.7 to 5.4 very long period (VLP) earthquakes. Collapses were exceptionally well recorded by global positioning system (GPS) and tilt instruments and represent unique natural kilometer-scale friction experiments. We model a piston collapsing into a magma reservoir. Pressure at the piston base and shear stress on its margin, governed by rate and state friction, balance its weight. Downward motion of the piston compresses the underlying magma, driving flow to the eruption. Monte Carlo estimation of unknowns validates laboratory friction parameters at the kilometer scale, including the magnitude of steady-state velocity weakening. The absence of accelerating precollapse deformation constrains d c to be ≤ 10 mm, potentially much less. These results support the use of laboratory friction laws and parameters for modeling earthquakes. We identify initial conditions and material and magma-system parameters that lead to episodic caldera collapse, revealing that small differences in eruptive vent elevation can lead to major differences in eruption volume and duration. Most historical basaltic caldera collapses were, at least partly, episodic, implying that the conditions for stick–slip derived here are commonly met in nature.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
-
Abstract In 2018 Kı̄lauea volcano erupted a decade's worth of basalt, given estimated magma supply rates, triggering caldera collapse. Yet, less than 2.5 years later Kı̄lauea re‐erupted. At the 2018 eruption onset, pressure within the summit reservoir was ∼20 MPa above magmastatic. By the onset of collapse this decreased by ∼17 MPa. Analysis of magma surges at the 2018 fissures, following collapse events, implies excess pressure at the eruption end of only ∼1 MPa. Given the new vent elevation, ∼11–12 MPa pressure increase was required to bring magma to the surface in December 2020. Analysis of Global Positioning System data between 8/2018 and 12/2020 shows there was a 73% probability that this condition was met at the onset of the 2020 eruption. Given a plausible range of possible vent elevations, there was a 40%–88% probability of sufficient pressure to bring magma to the surface 100 days before the eruption.more » « less
