An official website of the United States government
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.
more »
« less
A Simulator of Earthquakes and Aseismic Slip on a Heterogeneous Strike‐Slip Fault (HFQsim) With Static/Kinetic Friction and Temperature‐Dependent Creep
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
- Award ID(s):
- 2311206
- PAR ID:
- 10614871
- Publisher / Repository:
- DOI PREFIX: 10.1029
- Date Published:
- Journal Name:
- Journal of Geophysical Research: Solid Earth
- Volume:
- 130
- Issue:
- 6
- ISSN:
- 2169-9313
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
The potential relationship between surface creep and deeper geological processes is unclear, even on one of the world’s best-studied faults. From June to August 2021, a large creep event with surface slip of more than 16 mm was recorded on the Calaveras fault in California, part of the San Andreas fault system. This event initially appeared to be accompanied by along-fault migration of seismicity, suggesting it penetrated to depth. Other studies have suggested that surface creep events are likely a shallow feature, unrelated to deep seismicity. To provide more detail on the relationship between earthquakes, surface creep, and potential aseismic slip at seismogenic depth, we tripled the number of earthquakes in the Northern California Earthquake Catalog in the region of the creep event for all of 2021. This was accomplished by implementing earthquake detection techniques based on both template matching (EQCorrscan) and AI-based automatic earthquake phase picking (PhaseNet). After manual inspection, the detected earthquakes were first located using Hypoinverse and subsequently relocated via GrowClust. Our enhanced catalog indicates that the spatiotemporal pattern of earthquakes here is not strongly influenced by the creep event and is better explained by structural heterogeneity than transient stress changes, indicating a decoupling of seismicity rate and surficial creep on this major fault.more » « less
-
Recent seismic and geodetic observations indicate that interseismic creep rate varies in both time and space. The spatial extent of creep pinpoints locked asperities, while its temporary accelerations, known as slow-slip events, may trigger earthquakes. Although the conditions promoting fault creep are well-studied, the mechanisms for initiating episodic slow-slip events are enigmatic. Here we investigate surface deformation measured by radar interferometry along the central San Andreas Fault between 2003 and 2010 to constrain the temporal evolution of creep. We show that slow-slip events are ensembles of localized creep bursts that aseismically rupture isolated fault compartments. Using a rate-and-state friction model, we show that effective normal stress is temporally variable on the fault, and support this using seismic observations. We propose that compaction-driven elevated pore fluid pressure in the hydraulically isolated fault zone and subsequent frictional dilation cause the observed slow-slip episodes. We further suggest that the 2004 Mw 6 Parkfield earthquake might have been triggered by a slow-slip event, which increased the Coulomb failure stress by up to 0.45 bar per year. This implies that while creeping segments are suggested to act as seismic rupture barriers, slow-slip events on these zones might promote seismicity on adjacent locked segments.more » « less
-
Abstract The Calaveras Fault (CF) branches from the San Andreas Fault (SAF) near San Benito, extending sub‐parallel to the SAF for about 50 km with only 2–6 km separation and diverging northeastward. Both the SAF and CF are partially coupled, exhibit spatially variable aseismic creep and have hosted moderate to large earthquakes in recent decades. Understanding how slip partitions among the main fault strands of the SAF system and establishing their degree of coupling is crucial for seismic hazard evaluation. We perform a timeseries analysis using more than 5 years of Sentinel‐1 data covering the Bay Area (May 2015–October 2020), specifically targeting the spatiotemporal variations of creep rates around the SAF‐CF junction. We derive the surface creep rates from cross‐fault InSAR timeseries differences along the SAF and CF including adjacent Sargent and Quien Sabe Faults. We show that the variable creep rates (0–20 mm/yr) at the SAF‐CF junction are to first order controlled by the angle between the fault strike and the background stress orientation. We further examine the spatiotemporal variation of creep rates along the SAF and CF and find a multi‐annual coupling increase during 2016–2018 the subparallel sections of both faults, with the CF coupling change lagging behind the SAF by 3–6 months. Similar temporal variations are also observed in bothb‐values inferred from declustered seismicity and aseismic slip rates inferred from characteristic repeating earthquakes.more » « less
-
Abstract Mature strike‐slip faults are usually surrounded by a narrow zone of damaged rocks characterized by low seismic wave velocities. Observations of earthquakes along such faults indicate that seismicity is highly concentrated within this fault damage zone. However, the long‐term influence of the fault damage zone on complete earthquake cycles, that is, years to centuries, is not well understood. We simulate aseismic slip and dynamic earthquake rupture on a vertical strike‐slip fault surrounded by a fault damage zone for a thousand‐year timescale using fault zone material properties and geometries motivated by observations along major strike‐slip faults. The fault damage zone is approximated asan elastic layer with lower shear wave velocity than the surrounding rock. We find that dynamic wave reflections, whose characteristics are strongly dependent on the width and the rigidity contrast of the fault damage zone, have a prominent effect on the stressing history of the fault. The presence of elastic damage can partially explain the variability in the earthquake sizes and hypocenter locations along a single fault, which vary with fault damage zone depth, width and rigidity contrast from the host rock. The depth extent of the fault damage zone has a pronounced effect on the earthquake hypocenter locations, and shallower fault damage zones favor shallower hypocenters with a bimodal distribution of seismicity along depth. Our findings also suggest significant effects on the hypocenter distribution when the fault damage zone penetrates to the nucleation sites of earthquakes, likely being influenced by both lithological (material) and rheological (frictional) boundaries.more » « less
