Nucleation of Stick-Slip Instability Within a Large-Scale Experimental Fault: Effects of Stress Heterogeneities Due to Loading and Gouge Layer Compaction

Creep fronts and complexity in laboratory earthquake sequences illuminate delayed earthquake triggering

Earthquakes occur in clusters or sequences that arise from complex triggering mechanisms, but direct measurement of the slow subsurface slip responsible for delayed triggering is rarely possible. We investigate the origins of complexity and its relationship to heterogeneity using an experimental fault with two dominant seismic asperities. The fault is composed of quartz powder, a material common to natural faults, sandwiched between 760 mm long polymer blocks that deform the way 10 meters of rock would behave. We observe periodic repeating earthquakes that transition into aperiodic and complex sequences of fast and slow events. Neighboring earthquakes communicate via migrating slow slip, which resembles creep fronts observed in numerical simulations and on tectonic faults. Utilizing both local stress measurements and numerical simulations, we observe that the speed and strength of creep fronts are highly sensitive to fault stress levels left behind by previous earthquakes, and may serve as on-fault stress meters.

Laboratory earthquake experiments reproduce delayed earthquake triggering, similar to aftershocks, as a result of propagating slow slip fronts. The speed of the fronts can be highly sensitive to fault stress levels left behind by previous earthquakes.

Intermittent lab earthquakes in dynamically weakening fault gouge

Large and destructive earthquakes on mature faults in Earth's crust occur as slip in a layer of a fine granular material-fault gouge-produced by comminution during sliding^1,2. A range of insights into the frictional resistance of faults-one of the main factors controlling earthquake nucleation, dynamic propagation and arrest, and hence the destructive ground shaking of earthquakes^2,3-has been obtained in experiments with spatially uniform slip imposed in small samples^2,4-21. However, how various features of gouge friction combine to determine spontaneous progression of earthquakes is difficult to study in the lab owing to substantial challenges with sample sizes and adequate imaging²². Here, using lab experiments, we show that spontaneously propagating dynamic ruptures navigate a fault region with fine rock gouge through complex, intermittent slip processes with dramatic friction evolution. These include repeated arrest of rupture propagation caused by friction strengthening at lower slip rates and dynamic earthquake re-nucleation enabled by pronounced rapid friction weakening at higher slip rates consistent with flash heating^8,12,23. The spontaneous repeated weakening and strengthening of friction in fine rock gouge highlights the fundamental dependence of friction on slip rate and associated processes, such as shear heating, localization and delocalization of shear, and dilation and compaction of the shear layer^6,7,9-21. Our findings expand experimental support^9,11 of the concept that co-seismic weakening may enable earthquake rupture to break through stable fault regions^24,25, with substantial implications for seismic hazard.

Earthquake hazard and risk analysis for natural and induced seismicity: towards objective assessments in the face of uncertainty

The fundamental objective of earthquake engineering is to protect lives and livelihoods through the reduction of seismic risk. Directly or indirectly, this generally requires quantification of the risk, for which quantification of the seismic hazard is required as a basic input. Over the last several decades, the practice of seismic hazard analysis has evolved enormously, firstly with the introduction of a rational framework for handling the apparent randomness in earthquake processes, which also enabled risk assessments to consider both the severity and likelihood of earthquake effects. The next major evolutionary step was the identification of epistemic uncertainties related to incomplete knowledge, and the formulation of frameworks for both their quantification and their incorporation into hazard assessments. Despite these advances in the practice of seismic hazard analysis, it is not uncommon for the acceptance of seismic hazard estimates to be hindered by invalid comparisons, resistance to new information that challenges prevailing views, and attachment to previous estimates of the hazard. The challenge of achieving impartial acceptance of seismic hazard and risk estimates becomes even more acute in the case of earthquakes attributed to human activities. A more rational evaluation of seismic hazard and risk due to induced earthquakes may be facilitated by adopting, with appropriate adaptations, the advances in risk quantification and risk mitigation developed for natural seismicity. While such practices may provide an impartial starting point for decision making regarding risk mitigation measures, the most promising avenue to achieve broad societal acceptance of the risks associated with induced earthquakes is through effective regulation, which needs to be transparent, independent, and informed by risk considerations based on both sound seismological science and reliable earthquake engineering.

Two end-member earthquake preparations illuminated by foreshock activity on a meter-scale laboratory fault

The preparation process of natural earthquakes is still difficult to quantify and remains a subject of debate even with modern observational techniques. Here, we show that foreshock activity can shed light on understanding the earthquake preparation process based on results of meter-scale rock friction experiments. Experiments were conducted under two different fault surface conditions before each run: less heterogeneous fault without pre-existing gouge and more heterogeneous fault with pre-existing gouge. The results show that fewer foreshocks occurred along the less heterogeneous fault and were driven by preslip; in contrast, more foreshocks with a lower b value occurred along the more heterogeneous fault and showed features of cascade-up. We suggest that the fault surface condition and the stress redistribution caused by the ongoing fault slip mode control the earthquake preparation process, including the behavior of foreshock activity. Our findings imply that foreshock activity can be a key indicator for probing the fault conditions at present and in the future, and therefore useful for assessing earthquake hazard.