skip to main content


Title: Simulation of slip transients and earthquakes in finite thickness shear zones with a plastic formulation
We perform numerical experiments of damped quasi-dynamic fault slip that include a rate-and-state behavior at steady state to simulate earthquakes and a plastic rheology to model permanent strain. The model shear zone has a finite width which represents a natural fault zone. Here we reproduce fast and slow events that follow theoretical and observational scaling relationships for earthquakes and slow slip events (SSEs). We show that the transition between fast and slow slip occurs when the friction drop in the shear zone is equal to a critical value, Δμc. With lower friction drops, SSEs use nearly all of mechanical work to accumulate inelastic strain, while with higher friction drops fast slips use some of the mechanical work to slip frictionally. Our new formulation replaces the state evolution of rate and state by the stress evolution concurrent with accumulation of permanent damage in and around a fault zone.  more » « less
Award ID(s):
1722650
NSF-PAR ID:
10104321
Author(s) / Creator(s):
Date Published:
Journal Name:
Nature communications
Volume:
9
ISSN:
2041-1723
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    Tectonic faults fail through a spectrum of slip modes, ranging from slow aseismic creep to rapid slip during earthquakes. Understanding the seismic radiation emitted during these slip modes is key for advancing earthquake science and earthquake hazard assessment. In this work, we use laboratory friction experiments instrumented with ultrasonic sensors to document the seismic radiation properties of slow and fast laboratory earthquakes. Stick‐slip experiments were conducted at a constant loading rate of 8 μm/s and the normal stress was systematically increased from 7 to 15 MPa. We produced a full spectrum of slip modes by modulating the loading stiffness in tandem with the fault zone normal stress. Acoustic emission data were recorded continuously at 5 MHz. We demonstrate that the full continuum of slip modes radiate measurable high‐frequency energy between 100 and 500 kHz, including the slowest events that have peak fault slip rates <100 μm/s. The peak amplitude of the high‐frequency time‐domain signals scales systematically with fault slip velocity. Stable sliding experiments further support the connection between fault slip rate and high‐frequency radiation. Experiments demonstrate that the origin of the high‐frequency energy is fundamentally linked to changes in fault slip rate, shear strain, and breaking of contact junctions within the fault gouge. Our results suggest that having measurements close to the fault zone may be key for documenting seismic radiation properties and fully understanding the connection between different slip modes.

     
    more » « less
  2. Abstract

    Fluids influence fault zone strength and the occurrence of earthquakes, slow slip events, and aseismic slip. We introduce an earthquake sequence model with fault zone fluid transport, accounting for elastic, viscous, and plastic porosity evolution, with permeability having a power‐law dependence on porosity. Fluids, sourced at a constant rate below the seismogenic zone, ascend along the fault. While the modeling is done for a vertical strike‐slip fault with 2D antiplane shear deformation, the general behavior and processes are anticipated to apply also to subduction zones. The model produces large earthquakes in the seismogenic zone, whose recurrence interval is controlled in part by compaction‐driven pressurization and weakening. The model also produces a complex sequence of slow slip events (SSEs) beneath the seismogenic zone. The SSEs are initiated by compaction‐driven pressurization and weakening and stalled by dilatant suctions. Modeled SSE sequences include long‐term events lasting from a few months to years and very rapid short‐term events lasting for only a few days; slip is ∼1–10 cm. Despite ∼1–10 MPa pore pressure changes, porosity and permeability changes are small and hence fluid flux is relatively constant except in the immediate vicinity of slip fronts. This contrasts with alternative fault valving models that feature much larger changes in permeability from the evolution of pore connectivity. Our model demonstrates the important role that compaction and dilatancy have on fluid pressure and fault slip, with possible relevance to slow slip events in subduction zones and elsewhere.

     
    more » « less
  3. null (Ed.)
    Slow slip events (SSEs) at the northern Hikurangi subduction margin, New Zealand, are among the best-documented shallow SSEs on Earth. International Ocean Discovery Program Expeditions 372 and 375 were undertaken to investigate the processes and in situ conditions that underlie subduction zone SSEs at the northern Hikurangi Trough. We accomplished this goal by (1) coring and geophysical logging at four sites, including penetration of an active thrust fault (the Pāpaku fault) near the deformation front, the upper plate above the SSE source region, and the incoming sedimentary succession in the Hikurangi Trough and atop the Tūranganui Knoll seamount; and (2) installing borehole observatories in the Pāpaku fault and in the upper plate overlying the slow slip source region. Logging-while-drilling (LWD) data for this project were acquired as part of Expedition 372, and coring, wireline logging, and observatory installations were conducted during Expedition 375. Northern Hikurangi subduction margin SSEs recur every 1–2 y and thus provide an ideal opportunity to monitor deformation and associated changes in chemical and physical properties throughout the slow slip cycle. In situ measurements and sampling of material from the sedimentary section and oceanic basement of the subducting plate reveal the rock properties, composition, lithology, and structural character of material that is transported downdip into the SSE source region. A recent seafloor geodetic experiment raises the possibility that SSEs at northern Hikurangi may propagate to the trench, indicating that the shallow thrust fault (the Pāpaku fault) targeted during Expeditions 372 and 375 may also lie in the SSE rupture area and host a portion of the slip in these events. Hence, sampling and logging at this location provides insights into the composition, physical properties, and architecture of a shallow fault that may host slow slip. Expeditions 372 and 375 were designed to address three fundamental scientific objectives: 1. Characterize the state and composition of the incoming plate and shallow fault near the trench, which comprise the protolith and initial conditions for fault zone rock at greater depth and which may itself host shallow slow slip; 2. Characterize material properties, thermal regime, and stress conditions in the upper plate directly above the SSE source region; and 3. Install observatories in the Pāpaku fault near the deformation front and in the upper plate above the SSE source to measure temporal variations in deformation, temperature, and fluid flow. The observatories will monitor volumetric strain (via pore pressure as a proxy) and the evolution of physical, hydrological, and chemical properties throughout the SSE cycle. Together, the coring, logging, and observatory data will test a suite of hypotheses about the fundamental mechanics and behavior of SSEs and their relationship to great earthquakes along the subduction interface. 
    more » « less
  4. null (Ed.)
    Slow slip events (SSEs) at the northern Hikurangi subduction margin, New Zealand, are among the best-documented shallow SSEs on Earth. International Ocean Discovery Program Expedition 375 was undertaken to investigate the processes and in situ conditions that underlie subduction zone SSEs at the northern Hikurangi Trough by (1) coring at four sites, including an active fault near the deformation front, the upper plate above the high-slip SSE source region, and the incoming sedimentary succession in the Hikurangi Trough and atop the Tūranganui Knoll Seamount, and (2) installing borehole observatories in an active thrust near the deformation front and in the upper plate overlying the slow slip source region. Logging-while-drilling (LWD) data for this project were acquired as part of Expedition 372 (26 November 2017–4 January 2018; see the Expedition 372 Preliminary Report for further details on the LWD acquisition program). Northern Hikurangi subduction margin SSEs recur every 1–2 years and thus provide an ideal opportunity to monitor deformation and associated changes in chemical and physical properties throughout the slow slip cycle. Sampling of material from the sedimentary section and oceanic basement of the subducting plate reveals the rock properties, composition, lithology, and structural character of material that is transported downdip into the SSE source region. A recent seafloor geodetic experiment raises the possibility that SSEs at northern Hikurangi may propagate all the way to the trench, indicating that the shallow thrust fault zone targeted during Expedition 375 may also lie in the SSE rupture area. Hence, sampling at this location provides insights into the composition, physical properties, and architecture of a shallow fault that may host slow slip. Expedition 375 (together with the Hikurangi subduction LWD component of Expedition 372) was designed to address three fundamental scientific objectives: (1) characterize the state and composition of the incoming plate and shallow plate boundary fault near the trench, which comprise the protolith and initial conditions for fault zone rock at greater depth and which may itself host shallow slow slip; (2) characterize material properties, thermal regime, and stress conditions in the upper plate above the core of the SSE source region; and (3) install observatories at an active thrust near the deformation front and in the upper plate above the SSE source to measure temporal variations in deformation, temperature, and fluid flow. The observatories will monitor volumetric strain (via pore pressure as a proxy) and the evolution of physical, hydrological, and chemical properties throughout the SSE cycle. Together, the coring, logging, and observatory data will test a suite of hypotheses about the fundamental mechanics and behavior of SSEs and their relationship to great earthquakes along the subduction interface. 
    more » « less
  5. 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