An official website of the United States government
Abstract Fluid injection stimulates seismicity far from active tectonic regions. However, the details of how fluids modify on‐fault stresses and initiate seismic events remain poorly understood. We conducted laboratory experiments using a biaxial loading apparatus with a 3 m saw‐cut granite fault and compared events induced at different levels of background shear stress. Water was injected at 10 mL/min and normal stress was constant at 4 MPa. In all experiments, aseismic slip initiated on the fault near the location of fluid injection and dynamic rupture eventually initiated from within the aseismic slipping patch. When the fault was near critically stressed, seismic slip initiated only seconds after MPa‐level injection pressures were reached and the dynamic rupture propagated beyond the fluid pressure perturbed region. At lower stress levels, dynamic rupture initiated hundreds of seconds later and was limited to regions where aseismic slip had significantly redistributed stress from within the pressurized region to neighboring locked patches. We found that the initiation of slow slip was broadly consistent with a Coulomb failure stress, but that initiation of dynamic rupture required additional criteria to be met. Even high background stress levels required aseismic slip to modify on‐fault stress to meet initiation criteria. We also observed slow slip events prior to dynamic rupture. Overall, our experiments suggest that initial fault stress, relative to fault strength, is a critical factor in determining whether a fluid‐induced rupture will “runaway” or whether a fluid‐induced rupture will remain localized to the fluid pressurized region.
more »
« less
Microseismicity on Patches of Higher Compression During Larger‐Scale Earthquake Nucleation in a Rate‐and‐State Fault Model
Abstract While many large earthquakes are preceded by observable foreshocks, the mechanism responsible for the occurrence of these smaller‐scale seismic events remains uncertain. One physical explanation of foreshocks with growing support is that they are produced by the interaction of slow slip with fault heterogeneity. Inspired by the suggestion from laboratory experiments that foreshocks occur on fault asperities (bumps), we explore rate‐and‐state fault models with patches of higher normal stress embedded in a larger seismogenic region by conducting 3‐D numerical simulations of their behavior over long‐term sequences of aseismic and seismic slips. The models do produce smaller‐scale seismicity during the aseismic nucleation of larger‐scale seismic events. These smaller‐scale events have reasonable stress drops, despite the highly elevated compression assigned to the source patches. We find that the two main factors contributing to the reasonable stress drops are the significant extent of the rupture into the region surrounding the patches and the aseismic stress release just prior to the seismic events. The smaller‐scale seismicity can only occur if a sufficient separation in nucleation scales between the foreshock‐like events and mainshocks is achieved. Our modeling provides insight into the conditions conducive for generating foreshocks on both natural and laboratory faults.
more »
« less
- PAR ID:
- 10455930
- Publisher / Repository:
- DOI PREFIX: 10.1029
- Date Published:
- Journal Name:
- Journal of Geophysical Research: Solid Earth
- Volume:
- 124
- Issue:
- 2
- ISSN:
- 2169-9313
- Page Range / eLocation ID:
- p. 1962-1990
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
SUMMARY Earthquakes come in clusters formed of mostly aftershock sequences, swarms and occasional foreshock sequences. This clustering is thought to result either from stress transfer among faults, a process referred to as cascading, or from transient loading by aseismic slip (pre-slip, afterslip or slow slip events). The ETAS statistical model is often used to quantify the fraction of clustering due to stress transfer and to assess the eventual need for aseismic slip to explain foreshocks or swarms. Another popular model of clustering relies on the earthquake nucleation model derived from experimental rate-and-state friction. According to this model, earthquakes cluster because they are time-advanced by the stress change imparted by the mainshock. This model ignores stress interactions among aftershocks and cannot explain foreshocks or swarms in the absence of transient loading. Here, we analyse foreshock, swarm and aftershock sequences resulting from cascades in a Discrete Fault Network model governed by rate-and-state friction. We show that the model produces realistic swarms, foreshocks and aftershocks. The Omori law, characterizing the temporal decay of aftershocks, emerges in all simulations independently of the assumed initial condition. In our simulations, the Omori law results from the earthquake nucleation process due to rate and state friction and from the heterogeneous stress changes due to the coseismic stress transfers. By contrast, the inverse Omori law, which characterizes the accelerating rate of foreshocks, emerges only in the simulations with a dense enough fault system. A high-density complex fault zone favours fault interactions and the emergence of an accelerating sequence of foreshocks. Seismicity catalogues generated with our discrete fault network model can generally be fitted with the ETAS model but with some material differences. In the discrete fault network simulations, fault interactions are weaker in aftershock sequences because they occur in a broader zone of lower fault density and because of the depletion of critically stressed faults. The productivity of the cascading process is, therefore, significantly higher in foreshocks than in aftershocks if fault zone complexity is high. This effect is not captured by the ETAS model of fault interactions. It follows that a foreshock acceleration stronger than expected from ETAS statistics does not necessarily require aseismic slip preceding the mainshock (pre-slip). It can be a manifestation of a cascading process enhanced by the topological properties of the fault network. Similarly, earthquake swarms might not always imply transient loading by aseismic slip, as they can emerge from stress interactions.more » « less
-
Abstract Foreshocks can provide valuable information about possible nucleation process of a mainshock. However, their physical mechanisms are still under debate. In this study, we present a comprehensive analysis of the earthquake sequence preceding the 2010 Mw7.2 El Mayor‐Cucapah mainshock, including waveform detection of missing smaller events, relative relocation, and source parameter analysis. Based on a template matching method, we find a tenfold increase in the number of earthquakes than reported in the Southern California Seismic Network catalog. The entire sequence exhibits nearly continuous episodes of foreshocks that can be loosely separated into two active clusters. Relocated foreshocks show several seismicity streaks at depth, with a consistently active cluster at depths between 14 and 16 km where the mainshock was nucleated. Stress drop measurements from a spectral ratio approach based on empirical Green's functions show a range between 3.8 and 41.7 MPa with a median of 13.0 MPa and no clear temporal variations. The relocation results, together with the source patches estimated from earthquake corner frequencies, revealed a migration front toward the mainshock hypocenter within last 8 hr and a chain of active burst immediately 6 min prior to the mainshock. Our results support combined effects of aseismic slip and cascading failure on the evolution of foreshocks.more » « less
-
The temporal variation of elastic property of the bulk material surrounding the fault is considered an important contribution to the observed co-seismic velocity reduction and interseismic healing. Paglialunga et al. [2021] found that as fault normal stress increases, co-seismic velocity reduction becomes larger because more cracks reopen with higher stress drops. Larger normal stress can lead to smaller nucleation size and contribute to larger co-seismic slip. By contrast, with larger co-seismic velocity reduction and interseismic healing, more slow slip events can propagate in the seismogenic zone [Thakur and Huang, 2021], because the temporal velocity change related to fault zone damage modulates earthquake nucleation. Hence, fault normal stress and temporal damage zone structure evolution have opposite influences on the spatial distribution and recurrence intervals of earthquakes. We conducted 2-D anti-plane fully-dynamic seismic cycle simulations and explored the effects of fault normal stress on seismic cycle when there is coseismic damage and interseismic healing in the fault damage zone. The normal stress is in a range of 40-70 MPa and the co-seismic rigidity reduction is in a range of 5-8%. We find larger normal stress results in larger co-seismic slip and fewer slow slip events, while more co-seismic velocity reduction and interseismic healing leads to more partial ruptures as well as slow slip events. With the increase of both normal stress and seismic velocity change, more regular earthquakes occur and slow slip events gradually disappear. For the selected parameter space, the influence of seismic velocity change is not as significant as the effect of normal stress. However, fault zone maturity or the initial rigidity of fault damage zones should also affect the competitive relationship between normal stress and seismic velocity change, and we will characterize earthquakes and slow-slip events in immature and mature fault damage zones when both on-fault normal stress and off-fault seismic velocity vary over earthquake cycles.more » « less
-
Abstract We develop an earthquake simulator to study the partitioning of seismic/aseismic slip and dynamics of Earthquakes on a Heterogeneous strike‐slip Fault (HFQsim) using a generalized model of a discrete fault governed by static/dynamic friction and creep in an elastic half‐space. Previous versions of the simulator were shown to produce various realistic seismicity patterns (e.g., frequency‐magnitude event statistics, hypocenter and slip distributions, temporal occurrence) using friction levels and creep properties that vary in space but are fixed in time. The new simulator incorporates frictional heat generation by earthquake slip leading to temperature rises, subsequent diffusion cooling into the half space, and time‐dependent creep on the fault. The model assumes a power law dependence of creep velocity on the local shear stress, with temperature‐dependent coefficients based on the Arrhenius equation. Temperature rises due to seismic slip produce increased aseismic slip, which can lead to further stress concentrations, aftershocks, and heat generation in a feedback loop. The partitioning of seismic/aseismic slip and space‐time evolution of seismicity are strongly affected by the temperature changes on the fault. The results are also affected significantly by the difference between the static and kinetic friction levels. The model produces realistic spatio‐temporal distribution of seismicity, transient aseismic slip patterns, mainshock‐aftershock sequences, and a bimodal distribution of earthquakes with background and clustered events similar to observations. The HFQsim may be used to clarify relations between fault properties and different features of seismicity and aseismic slip, and to improve the understanding of failure patterns preceding large earthquakes.more » « less
