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.
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- Award ID(s):
- 1947448
- NSF-PAR ID:
- 10193949
- Publisher / Repository:
- Nature Publishing Group
- Date Published:
- Journal Name:
- Nature Communications
- Volume:
- 11
- Issue:
- 1
- ISSN:
- 2041-1723
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
more » « less
Abstract -
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
-
null (Ed.)International Ocean Discovery Program (IODP) Expedition 372 combines two research topics, slow slip events (SSEs) on subduction faults (IODP Proposal 781A-Full) and actively deforming gas hydrate–bearing landslides (Proposal 841-APL). Our study area on the Hikurangi margin east of New Zealand provides unique locations for addressing both research topics. Gas hydrates have long been suspected of being involved in seafloor failure; not much evidence, however, has been found to date for gas hydrate–related submarine landslides. Solid, icelike gas hydrate in sediment pores is generally thought to increase seafloor strength, as confirmed by a number of laboratory measurements. Dissociation of gas hydrate to water and overpressured gas, on the other hand, may destabilize the seafloor, potentially causing submarine landslides. The Tuaheni Landslide Complex on the Hikurangi margin shows evidence for active, creeping deformation. Intriguingly, the landward edge of creeping coincides with the pinchout of the base of gas hydrate stability (BGHS) on the seafloor. We therefore hypothesize that gas hydrate may be linked to creeping by (1) repeated small-scale sliding at the BGHS, in a variation of the conventional model linking gas hydrates and seafloor failure; (2) overpressure at the BGHS due to a permeability reduction linked to gas hydrates, which may lead to hydrofracturing, weakening the seafloor and allowing transmission of pressure into the gas hydrate stability zone; or (3) icelike viscous deformation of gas hydrates in sediment pores, similar to onshore rock glaciers. The latter two processes imply that gas hydrate itself is involved in creeping, constituting a paradigm shift in relating gas hydrates to submarine slope failure. Alternatively, creeping may not be related to gas hydrates but instead be caused by repeated pressure pulses or linked to earthquake-related liquefaction. We have devised a coring and logging program to test our hypotheses. SSEs at subduction zones are an enigmatic form of creeping fault behavior. At the northern Hikurangi subduction margin (HSM), they are among the best-documented and shallowest on Earth. They recur about every 2 y and may extend close to the trench, where clastic and pelagic sediments about 1.0–1.5 km thick overlie the subducting, seamount-studded Hikurangi Plateau. The northern HSM thus provides an excellent setting to use IODP capabilities to discern the mechanisms behind slow slip fault behavior, as proposed in IODP Proposal 781A-Full. The objectives of Proposal 781A-Full will be implemented across two related IODP expeditions, 372 and 375. Expedition 372 will undertake logging while drilling (LWD) at three sites targeting the upper plate (midslope basin, proposed Site HSM-01A), the frontal thrust (proposed Site HSM-18A), and the subducting section in the trench (proposed Site HSM-05A). Expedition 375 will undertake coring at the same sites, as well as an additional seamount site on the subducting plate, and implement the borehole observatory objectives. The data from each expedition will be shared between both scientific parties. Collectively, the LWD and coring data will be used to (1) characterize the compositional, structural, thermal, and diagenetic state of the incoming plate and the shallow plate boundary fault near the trench, which comprise the protolith and initial conditions for fault zone rock associated with SSEs at greater depth, and (2) characterize the material properties, thermal regime, and stress conditions in the upper plate above the SSE source region. These data will be used during Expedition 375 to guide the installation of CORK observatories at the frontal thrust and in the upper plate above the SSE source to monitor temporal variations in deformation, fluid flow, seismicity, and physical and chemical properties throughout the SSE cycle (Saffer et al., 2017). Together, these 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
-
more » « less
Abstract Localized frictional sliding on faults in the continental crust transitions at depth to distributed deformation in viscous shear zones. This brittle‐ductile transition (BDT), and/or the transition from velocity‐weakening (VW) to velocity‐strengthening (VS) friction, are controlled by the lithospheric thermal structure and composition. Here, we investigate these transitions, and their effect on the depth extent of earthquakes, using 2D antiplane shear simulations of a strike‐slip fault with rate‐and‐state friction. The off‐fault material is viscoelastic, with temperature‐dependent dislocation creep. We solve the heat equation for temperature, accounting for frictional and viscous shear heating that creates a thermal anomaly relative to the ambient geotherm which reduces viscosity and facilitates viscous flow. We explore several geotherms and effective normal stress distributions (by changing pore pressure), quantifying the thermal anomaly, seismic and aseismic slip, and the transition from frictional sliding to viscous flow. The thermal anomaly can reach several hundred degrees below the seismogenic zone in models with hydrostatic pressure but is smaller for higher pressure (and these high‐pressure models are most consistent with San Andreas Fault heat flow constraints). Shear heating raises the BDT, sometimes to where it limits rupture depth rather than the frictional VW‐to‐VS transition. Our thermomechanical modeling framework can be used to evaluate lithospheric rheology and thermal models through predictions of earthquake ruptures, postseismic and interseismic crustal deformation, heat flow, and the geological structures that reflect the complex deformation beneath faults.
-
more » « less
Abstract We study the mechanical response of two‐dimensional vertical strike‐slip fault to coseismic damage evolution and interseismic healing of fault damage zones by simulating fully dynamic earthquake cycles. Our models show that fault zone structure evolution during the seismic cycle can have pronounced effects on mechanical behavior of locked and creeping fault segments. Immature fault damage zone models exhibit small and moderate subsurface earthquakes with irregular recurrence intervals and abundance of slow‐slip events during the interseismic period. In contrast, mature fault damage zone models host pulse‐like earthquake ruptures that can propagate to the surface and extend throughout the seismogenic zone, resulting in large stress drop, characteristic rupture extents, and regular recurrence intervals. Our results suggest that interseismic healing and coseismic damage accumulation in fault zones can explain the observed differences of earthquake behaviors between mature and immature fault zones and indicate a link between regional seismic hazard and fault structural maturity.
