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
Dynamic Rupture Simulations of Caldera Collapse Earthquakes: Effects of Wave Radiation, Magma Viscosity, and Evidence of Complex Nucleation at Kı̄lauea 2018
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
- PAR ID:
- 10504634
- Publisher / Repository:
- American Geophysical Union
- Date Published:
- Journal Name:
- Journal of Geophysical Research: Solid Earth
- Volume:
- 129
- Issue:
- 4
- ISSN:
- 2169-9313
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
null (Ed.)ABSTRACT Large earthquakes on strike-slip faults often rupture multiple fault segments by jumping over stepovers. Previous studies, based on field observations or numerical modeling with a homogeneous initial stress field, have suggested that stepovers more than ∼5 km wide would stop the propagation of rupture, but many exceptions have been observed in recent years. Here, we integrate a dynamic rupture model with a long-term fault stress model to explore the effects of background stress perturbation on rupture propagation across stepovers along strike-slip faults. Our long-term fault models simulate steady-state stress perturbation around stepovers. Considering such stress perturbation in dynamic rupture models leads to prediction of larger distance a dynamic rupture can jump over stepovers: over 15 km for a releasing stepover or 7 km for a restraining stepover, comparing with the 5 km limit in models with the same fault geometry and frictional property but assuming a homogeneous initial stress. The effect of steady-state stress perturbations is stronger in an overlapping stepover than in an underlapping stepover. The maximum jumping distance can reach 20 km in an overlapping releasing stepover with low-static frictional coefficients. These results are useful for estimating the maximum length of potential fault ruptures and assessing seismic hazard.more » « less
-
Predicting the onset and timing of fault failure is one of the ultimate goals of seismology. However, our current understanding of the earthquake preparation and nucleation process is limited. One direction towards understanding this process is looking at precursory signals preceding large earthquakes. Previous laboratory experiments have studied robust precursory signals, observed as temporal changes in pressure and shear wave velocities during the seismic cycle. The effects of such precursory velocity changes on the seismic cycle are not well understood. We use numerical models to simulate fully-dynamic earthquake cycles in 2D strike-slip fault systems with antiplane geometry, surrounded by a narrow fault-parallel damage zone. By imposing shear wave velocity changes inside fault damage zones, we investigate the effects of these precursors on multiple stages of the seismic cycle, including nucleation, coseismic, postseismic, and interseismic stages. Our modeling results show a wide spectrum of fault-slip behaviors including fast earthquakes, slow-slip events, and variable creep. One primary effect of the imposed velocity precursor is the facilitation of the otherwise slow-slip event to grow into a fully dynamic earthquake. Furthermore, the onset time of these precursors have significant effects on the nucleation phase of the earthquakes, and earlier onset of precursors causes the earthquakes to nucleate earlier with a smaller nucleation size. Our results highlight the importance of short and long-term monitoring of fault zone structures for better assessment of regional seismic hazard.more » « less
-
Abstract Earthquake nucleation is a fundamental problem in earthquake science and has practical implications for forecasting seismic hazards. Laboratory experiments performed on large, meter‐scale fault systems offer unique insights into the nucleation process because the migration and expansion of the nucleation zone can be precisely detected, measured, and characterized using arrays of local strain and slip measurements. We report on a series of laboratory experiments conducted on a 1‐m direct shear machine. We sheared layers of quartz gouge between roughened acrylic forcing blocks over a range of normal stresses between 3 and 12 MPa, generating a spectrum of slip modes, ranging from aseismic creep to fast‐dynamic rupture. Co‐seismic slip, peak slip velocity, and high‐frequency acoustic energy content of laboratory earthquakes increases systematically with both cumulative fault slip and normal stress. Slower and smaller laboratory earthquake sequences have larger nucleation zones, creep more during their inter‐seismic period, and are deficient in high‐frequency energy compared to larger and faster rupture sequences. We find that the critical nucleation length scale,H*, scales inversely with cumulative fault slip and normal stress. A reduction inH*and an increase in event size can be explained by a decrease in the critical slip distance,Dc, or an increase in the frictional rate parameterb–aand is likely driven by shear localization. Together, our results indicate that homogeneous, mature fault zones that have undergone more cumulative fault slip are expected to have smallerH*and can more easily host dynamic instabilities, relative to immature faults.more » « less
-
null (Ed.)The fault damage zone is a well-known structure of localized deformation around faults. Its material properties evolve over earthquake cycles due to coseismic damage accumulation and interseismic healing. We will present fully dynamic earthquake cycle simulations to show how the styles of earthquake nucleation and rupture propagation change as fault zone material properties vary temporally. First, we will focus on the influence of fault zone structural maturity quantified by near-fault seismic wave velocities in simulations. The simulations show that immature fault zones promote small and moderate subsurface earthquakes with irregular recurrence intervals, whereas mature fault zones host pulse-like earthquake rupture that can propagate to the surface, extend throughout the seismogenic zone, and occur at regular intervals. The interseismic healing in immature fault zones plays a key role in allowing the development of aseismic slip episodes including slow-slip events and creep, which can propagate into the seismogenic zone, and thus limit the sizes of subsequent earthquakes by releasing fault stress. In the second part, we will discuss how the precursory changes of seismic wave velocities of fault damage zones may affect earthquake nucleation process. Both laboratory experiments and seismic observations show that the abrupt earthquake failure can be preceded by accelerated fault deformation and the accompanying velocity reduction of near-fault rocks. We will use earthquake cycle simulations to systematically test the effects of timing and amplitudes of such precursory velocity changes. Our simulations will provide new insights into the interplay between fault zone structure and earthquake nucleation process, which can be used to guide future real-time monitoring of major fault zones.more » « less
An official website of the United States government

