Earthquake swarms attributed to subsurface fluid injection are usually assumed to occur on faults destabilized by increased pore-fluid pressures. However, fluid injection could also activate aseismic slip, which might outpace pore-fluid migration and transmit earthquake-triggering stress changes beyond the fluid-pressurized region. We tested this theoretical prediction against data derived from fluid-injection experiments that activated and measured slow, aseismic slip on preexisting, shallow faults. We found that the pore pressure and slip history imply a fault whose strength is the product of a slip-weakening friction coefficient and the local effective normal stress. Using a coupled shear-rupture model, we derived constraints on the hydromechanical parameters of the actively deforming fault. The inferred aseismic rupture front propagates faster and to larger distances than the diffusion of pressurized pore fluid.
more »
« less
The Role of Background Stress State in Fluid‐Induced Aseismic Slip and Dynamic Rupture on a 3‐m Laboratory Fault
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- Award ID(s):
- 1847139
- NSF-PAR ID:
- 10420773
- Publisher / Repository:
- DOI PREFIX: 10.1029
- Date Published:
- Journal Name:
- Journal of Geophysical Research: Solid Earth
- Volume:
- 127
- Issue:
- 8
- ISSN:
- 2169-9313
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Fluids are known to trigger a broad range of slip events, from slow, creeping transients to dynamic earthquake ruptures. Yet, the detailed mechanics underlying these processes and the conditions leading to different rupture behaviors are not well understood. Here, we use a laboratory earthquake setup, capable of injecting pressurized fluids, to compare the rupture behavior for different rates of fluid injection, slow (megapascals per hour) versus fast (megapascals per second). We find that for the fast injection rates, dynamic ruptures are triggered at lower pressure levels and over spatial scales much smaller than the quasistatic theoretical estimates of nucleation sizes, suggesting that such fast injection rates constitute dynamic loading. In contrast, the relatively slow injection rates result in gradual nucleation processes, with the fluid spreading along the interface and causing stress changes consistent with gradually accelerating slow slip. The resulting dynamic ruptures propagating over wetted interfaces exhibit dynamic stress drops almost twice as large as those over the dry interfaces. These results suggest the need to take into account the rate of the pore-pressure increase when considering nucleation processes and motivate further investigation on how friction properties depend on the presence of fluids.more » « less
-
more » « less
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
Abstract There is a growing recognition that subsurface fluid injection can produce not only earthquakes, but also aseismic slip on faults. A major challenge in understanding interactions between injection-related aseismic and seismic slip on faults is identifying aseismic slip on the field scale, given that most monitored fields are only equipped with seismic arrays. We present a modeling workflow for evaluating the possibility of aseismic slip, given observational constraints on the spatial-temporal distribution of microseismicity, injection rate, and wellhead pressure. Our numerical model simultaneously simulates discrete off-fault microseismic events and aseismic slip on a main fault during fluid injection. We apply the workflow to the 2012 Enhanced Geothermal System injection episode at Cooper Basin, Australia, which aimed to stimulate a water-saturated granitic reservoir containing a highly permeable (
$$k = 10^{-13} - 10^{-12}$$ $$\hbox {m}{^2}$$ -
more » « less
Oceanic transform faults display a unique combination of seismic and aseismic slip behavior, including a large globally averaged seismic deficit, and the local occurrence of repeating magnitude (M)
earthquakes with abundant foreshocks and seismic swarms, as on the Gofar transform of the East Pacific Rise and the Blanco Ridge in the northeast Pacific Ocean. However, the underlying mechanisms that govern the partitioning between seismic and aseismic slip and their interaction remain unclear. Here we present a numerical modeling study of earthquake sequences and aseismic transient slip on oceanic transform faults. In the model, strong dilatancy strengthening, supported by seismic imaging that indicates enhanced fluid-filled porosity and possible hydrothermal circulation down to the brittle–ductile transition, effectively stabilizes along-strike seismic rupture propagation and results in rupture barriers where aseismic transients arise episodically. The modeled slow slip migrates along the barrier zones at speeds ∼10 to 600 m/h, spatiotemporally correlated with the observed migration of seismic swarms on the Gofar transform. Our model thus suggests the possible prevalence of episodic aseismic transients in M rupture barrier zones that host active swarms on oceanic transform faults and provides candidates for future seafloor geodesy experiments to verify the relation between aseismic fault slip, earthquake swarms, and fault zone hydromechanical properties.
