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 (
Pore pressure in aquifers confined below a cryosphere will increase as Mars cools and the cryosphere thickens. Increased pore pressure decreases the effective stress and hence promotes seismicity. We calculate the rate of pore pressure change from cooling of Mars's interior and the modulation of pore pressure from solar and Phobos tides and barometric loading. Using the time‐varying pressure and tidal stresses, we compute Coulomb stress changes and the expected seismicity rate from a rate‐and‐state friction model. Seismicity rate will vary by several tens of percent to 2 orders of magnitude if the mean pore pressure is within 0.2 and 0.01 MPa of lithostatic, respectively. Seismic events promoted by high pore pressure may be tremor‐like. Documenting (or not) tidally modulated shallow seismicity would provide evidence for (or against) water‐filled confined aquifers, that pore pressure is high, and that the state of stress is close to failure—with implications for processes that can deliver water to the Martian surface.
more » « less- NSF-PAR ID:
- 10375352
- Publisher / Repository:
- DOI PREFIX: 10.1029
- Date Published:
- Journal Name:
- Geophysical Research Letters
- Volume:
- 46
- Issue:
- 12
- ISSN:
- 0094-8276
- Page Range / eLocation ID:
- p. 6333-6340
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
more » « less
Abstract $$k = 10^{-13} - 10^{-12}$$ $$\hbox {m}{^2}$$ -
more » « less
Abstract The constitutive behavior of faults is central to many interconnected aspects of earthquake science, from fault dynamics to induced seismicity, to seismic hazards characterization. Yet, a friction law applicable to the range of temperatures found in the brittle crust and upper mantle is still missing. In particular, rocks often exhibit a transition from steady‐state velocity‐strengthening at room temperature to velocity‐weakening in warmer conditions that is poorly understood. Here, we investigate the effect of competing healing mechanisms on the evolution of frictional resistance in a physical model of rate‐, state‐, and temperature‐dependent friction. The yield strength for fault slip depends on the real area of contact, which is modulated by the competition between the growth and erosion of interfacial micro‐asperities. Incorporating multiple healing mechanisms and rock‐forming minerals with different thermodynamic properties allows a transition of the velocity‐ and temperature‐dependence of friction at steady‐state with varying temperatures. We explain the mechanical data for granite, pyroxene, amphibole, shale, and natural fault gouges with activation energies and stress power exponent for weakening of 10–50 kJ/mol and 55–150, respectively, compatible with subcritical crack growth and inter‐granular flow in the active slip zone. Activation energies for the time‐dependent healing process in the range 90–130 kJ/mol in dry conditions and 20–65 kJ/mol in wet conditions indicate the prominence of viscoelastic collapse of microasperities in the absence of water and of pressure‐solution creep, crack healing, and cementation when assisted by pore fluids.
-
more » « less
Abstract We investigate the relative importance of injection and production on the spatial‐temporal distribution of induced seismicity at the Raft River geothermal field. We use time‐series of InSAR measurements to document surface deformation and calibrate a hydro‐mechanical model to estimate effective stress changes imparted by injection and production. Seismicity, located predominantly in the basement, is induced primarily by poroelastic stresses from cold water reinjection into a shallower reservoir. The poroelastic effect of production from a deeper reservoir is minimal and inconsistent with observed seismicity, as is pore‐pressure‐diffusion in the basement and along reactivated faults. We estimate an initial strength excess of ∼20 kPa in the basement and sedimentary cover, but the seismicity rate in the sedimentary cover is four times lower, reflecting lower density of seed‐points for earthquake nucleation. Our modeling workflow could be used to assess the impact of fluid extraction or injection on seismicity and help design or guide operations.
-
more » « less
Abstract To better quantify how injection, prior seismicity, and fault properties control rupture growth and propagation of induced earthquakes, we perform a finite‐fault slip inversion on a
M w4.0 earthquake that occurred in April 2015, the largest earthquake in an induced sequence near Guthrie, Oklahoma. The slip inversion reveals a rupture with slip patches that are anti‐correlated to the locations of prior seismicity. The prior seismicity driven by low pore pressure changes and static stress changes occurred on weaker portions of the fault, while theM w4.0 earthquake likely ruptured relatively stronger portions of the fault. To resolve if pore pressure changes or the initial underlying stress distribution and fault strength controlled the final slip distribution of the GuthrieM w4.0 earthquake, we compare strike‐slip events of similar magnitude from tectonically active regions and previously inactive regions. Earthquakes on reactivated faults exhibit different slip distributions than active regions, they have more prominent and well separated slip patches, a behavior often associated with faults of lower fault maturity. Pore pressure shows little effect on the distributions. These observations suggest that the initial underlying stress distribution and fault strength of reactivated faults in low deformation regions is the primary controlling factor of the slip distribution with pore pressure perturbations and earthquake interactions being secondary. Therefore, GuthrieM w4.0 earthquakes slip distribution was enhanced by pore‐pressure perturbations and earthquake interactions by creating an optimal stress state for its failure, but the slip distribution itself is controlled by its fault's initial stress and strength state. -
more » « less
Abstract Fault-zone fluids control effective normal stress and fault strength. While most earthquake models assume a fixed pore fluid pressure distribution, geologists have documented fault valving behavior, that is, cyclic changes in pressure and unsteady fluid migration along faults. Here we quantify fault valving through 2-D antiplane shear simulations of earthquake sequences on a strike-slip fault with rate-and-state friction, upward Darcy flow along a permeable fault zone, and permeability evolution. Fluid overpressure develops during the interseismic period, when healing/sealing reduces fault permeability, and is released after earthquakes enhance permeability. Coupling between fluid flow, permeability and pressure evolution, and slip produces fluid-driven aseismic slip near the base of the seismogenic zone and earthquake swarms within the seismogenic zone, as ascending fluids pressurize and weaken the fault. This model might explain observations of late interseismic fault unlocking, slow slip and creep transients, swarm seismicity, and rapid pressure/stress transmission in induced seismicity sequences.
