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
Repeating caldera collapse events constrain fault friction at the kilometer scale
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
- Award ID(s):
- 2040425
- PAR ID:
- 10316778
- Date Published:
- Journal Name:
- Proceedings of the National Academy of Sciences
- Volume:
- 118
- Issue:
- 30
- ISSN:
- 0027-8424
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
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
-
The science of volcanology advances disproportionately during exceptionally large or well-observed eruptions. The 2018 eruption of Kīlauea Volcano (Hawai‘i) was its most impactful in centuries, involving an outpouring of more than one cubic kilometer of basalt, a magnitude 7 flank earthquake, and the volcano's largest summit collapse since at least the nineteenth century. Eruptive activity was documented in detail, yielding new insights into large caldera-rift eruptions; the geometry of a shallow magma storage-transport system and its interaction with rift zone tectonics; mechanisms of basaltic tephra-producing explosions; caldera collapse mechanics; and the dynamics of fissure eruptions and high-volume lava flows. Insights are broadly applicable to a range of volcanic systems and should reduce risk from future eruptions. Multidisciplinary collaboration will be required to fully leverage the diversity of monitoring data to address many of the most important outstanding questions. ▪ Unprecedented observations of a caldera collapse and coupled rift zone eruption yield new opportunities for advancing volcano science. ▪ Magma flow to a low-elevation rift zone vent triggered quasi-periodic step-like collapse of a summit caldera, which pressurized the magma system and sustained the eruption. ▪ Kīlauea's magmatic-tectonic system is tightly interconnected over tens of kilometers, with complex feedback mechanisms and interrelated hazards over widely varying time scales. ▪ The eruption revealed magma stored in diverse locations, volumes, and compositions, not only beneath the summit but also within the volcano's most active rift zone. Expected final online publication date for the Annual Review of Earth and Planetary Sciences, Volume 52 is May 2024. Please see http://www.annualreviews.org/page/journal/pubdates for revised estimates.more » « less
-
Abstract The empirical rate‐ and state‐dependent friction law is widely used to explain the frictional resistance of rocks. However, the constitutive parameters vary with temperature and sliding velocity, preventing extrapolation of laboratory results to natural conditions. Here, we explain the frictional properties of natural gouge from the San Andreas Fault, Alpine Fault, and the Nankai Trough from room temperature to ∼300°C for a wide range of slip‐rates with constant constitutive parameters by invoking the competition between two healing mechanisms with different thermodynamic properties. A transition from velocity‐strengthening to velocity‐weakening at steady‐state can be attained either by decreasing the slip‐rate or by increasing temperature. Our study provides a framework to understand the physics underlying the slip‐rate and state dependence of friction and the dependence of frictional properties on ambient physical conditions.more » « less
-
Abstract Establishing a constitutive law for fault friction is a crucial objective of earthquake science. However, the complex frictional behavior of natural and synthetic gouges in laboratory experiments eludes explanations. Here, we present a constitutive framework that elucidates the rate, state, and temperature dependence of fault friction under the relevant sliding velocities and temperatures of the brittle lithosphere during seismic cycles. The competition between healing mechanisms, such as viscoelastic collapse, pressure‐solution creep, and crack sealing, explains the low‐temperature stability transition from steady‐state velocity‐strengthening to velocity‐weakening as a function of slip‐rate and temperature. In addition, capturing the transition from cataclastic flow to semi‐brittle creep accounts for the stabilization of fault slip at elevated temperatures. We calibrate the model using extensive laboratory data on synthetic albite and granite gouge, and on natural samples from the Alpine Fault and the Mugi Mélange in the Shimanto accretionary complex in Japan. The constitutive model consistently explains the evolving frictional response of fault gouge from room temperature to 600°C for sliding velocities ranging from nanometers to millimeters per second. The frictional response of faults can be uniquely determined by the in situ lithology and the prevailing hydrothermal conditions.more » « less
An official website of the United States government

