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.
Inversions of InSAR ground deformation in the Delaware Basin have revealed an aseismic slip on semi‐optimally oriented normal faults located close to disposal wells. The slip, occurring over 3–5 years, extends approximately 1 km down‐dip, over 10 km along strike, and reaches 25 cm. We develop and calibrate 2D and pseudo‐3D coupled pore pressure diffusion and rate‐state models with velocity‐strengthening friction tailored to this unique height‐bounded fault geometry. Pressure diffusion is limited to a high‐permeability fault damage zone, and the net influx of fluid is adjusted to match the observed slip. A 1–2 MPa pressure increase initiates slip, with ∼5 MPa additional pressure increase required to produce ∼20 cm slip. Most slip occurs at approximately constant friction. Fault zone permeability must exceed ∼10−13 m2to match the along‐strike extent of slip. Models of the type developed here can be used to operationally manage injection‐induced aseismic slip.
more » « less- Award ID(s):
- 1947448
- NSF-PAR ID:
- 10444636
- Publisher / Repository:
- DOI PREFIX: 10.1029
- Date Published:
- Journal Name:
- Geophysical Research Letters
- Volume:
- 49
- Issue:
- 11
- ISSN:
- 0094-8276
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
more » « less
Abstract -
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.
-
null (Ed.)The Nankai Trough Seismogenic Zone Experiment (NanTroSEIZE) is a coordinated, multiexpedition International Ocean Discovery Program (IODP) drilling project designed to investigate fault mechanics and seismogenesis along subduction megathrusts through direct sampling, in situ measurements, and long-term monitoring in conjunction with allied laboratory and numerical modeling studies. The fundamental scientific objectives of the NanTroSEIZE drilling project include characterizing the nature of fault slip and strain accumulation, fault and wall rock composition, fault architecture, and state variables throughout the active plate boundary system. IODP Expedition 365 is part of NanTroSEIZE Stage 3, with the following primary objectives: 1. Retrieval of a temporary observatory at Site C0010 that began monitoring temperature and pore pressure within the major splay thrust fault (termed the “megasplay”) at 400 meters below seafloor in November 2010. 2. Deployment of a complex long-term borehole monitoring system (LTBMS) designed to be connected to the Dense Oceanfloor Network System for Earthquakes and Tsunamis (DONET) seafloor cabled observatory network postexpedition. The LTBMS incorporates multilevel pore pressure sensing, a volumetric strainmeter, tiltmeter, geophone, broadband seismometer, accelerometer, and thermistor string. Together with an existing observatory at Integrated Ocean Drilling Program Site C0002 and a planned future installation near the trench, the Site C0010 observatory allows monitoring within and above regions of contrasting behavior of the megasplay fault and the plate boundary as a whole. These include a site above the updip edge of the locked zone (Site C0002), a shallow site in the megasplay fault zone and its footwall (Site C0010), and a site at the tip of the accretionary prism (possible future installation at Integrated Ocean Drilling Program Site C0006). Together, this suite of observatories has the potential to capture deformation spanning a wide range of timescales (e.g., seismic and microseismic activity, slow slip, and interseismic strain accumulation) across a transect from near-trench to the seismogenic zone. Site C0010 is located 3.5 km along strike to the southwest of Integrated Ocean Drilling Program Site C0004. The site was drilled and cased during Integrated Ocean Drilling Program Expedition 319, with casing screens spanning a ~20 m interval that includes the megasplay fault, and suspended with a temporary instrument package (a “SmartPlug”), which included pressure and temperature sensors. During Integrated Ocean Drilling Program Expedition 332 in late 2010, the instrument package was replaced with an upgraded sensor package (the “GeniusPlug”), which included a set of geochemical and biological experiments in addition to pressure and temperature sensors. Expedition 365 achieved its primary scientific and operational objectives, including recovery of the GeniusPlug with a >5 y record of pressure and temperature conditions within the shallow megasplay fault zone, geochemical samples, and its in situ microbial colonization experiment; and installation of the LTBMS. The pressure records from the GeniusPlug include high-quality records of formation and seafloor responses to multiple fault slip events, including the 11 March 2011 Tohoku M9 and 1 April 2016 Mie-ken Nanto-oki M6 earthquakes. The geochemical sampling coils yielded in situ pore fluids from the splay fault zone, and microorganisms were successfully cultivated from the colonization unit. The complex sensor array, in combination with the multilevel hole completion, is one of the most ambitious and sophisticated observatory installations in scientific ocean drilling (similar to that in Hole C0002G, deployed in 2010). Overall, the installation went smoothly, efficiently, and ahead of schedule. The extra time afforded by the efficient observatory deployment was used for coring in Holes C0010B–C0010E. Despite challenging hole conditions, the depth interval corresponding to the screened casing across the megasplay fault was successfully sampled in Hole C0010C, and the footwall of the megasplay was sampled in Hole C0010E, with >50% recovery for both zones. In the hanging wall of the megasplay fault (Holes C0010C and C0010D), we recovered indurated silty clay with occasional ash layers and sedimentary breccias. Mudstones show different degrees of deformation spanning from occasional fractures to intervals of densely fractured scaly claystones of up to >10 cm thickness. Sparse faulting with low displacement (usually <2 cm) is seen in core and exhibits primarily normal and, rarely, reversed sense of slip. When present, ash was entrained along fractures and faults. In Hole C0010E, the footwall to the megasplay fault was recovered. Sediments are horizontally to gently dipping and mainly comprise silt of olive-gray color. The hanging wall sediments recovered in Holes C0010C–C0010D range in age from 3.79 to 5.59 Ma and have been thrust over the younger footwall sediments in Hole C0010E, ranging in age from 1.56 to 1.67 Ma. The deposits of the underthrust sediment prism are less indurated than the hanging wall mudstones and show lamination on a centimeter scale. The material is less intensely deformed than the mudstones, and apart from occasional fracturation (some of it being drilling disturbance), evidence of structural features is absent.more » « less
-
more » « less
Abstract It is widely recognized that fluid injection can trigger aseismic fault slip. However, the processes by which the fluid‐rock interactions facilitate or inhibit slip are poorly understood and some are oversimplified in most models of injection‐induced slip. In this study, we perform a 2D anti‐plane shear investigation of aseismic slip that occurs in response to fluid injection into a permeable fault governed by rate‐and‐state friction. We account for porosity and permeability changes that accompany slip, including dilatancy, and quantify how these processes affect pore pressure diffusion, which couples to aseismic slip. Fault response to injection has two phases. In the first phase, slip is negligible and pore pressure closely follows the standard linear diffusion model. Pressurization eventually triggers aseismic slip close to the injection site. In the second phase, aseismic slip front expands outward and dilatancy causes pore pressure to depart from the linear diffusion model. We quantify how prestress, injection rate, permeability and other fluid transport properties affect the slip front migration rate, finding rates ranging from 10 to 1,000 m/day for typical parameters. The migration rate is strongly influenced by the fault's closeness to failure and injection rate. The total slip on the fault, on the other hand, is primarily determined by the injected volume, with minimal sensitivity to injection rate. Additionally, we show that when dilatancy is neglected, slip front migration rate and total slip can be several times higher. Our modeling demonstrates that porosity and permeability evolution, especially dilatancy, fundamentally alters how faults respond to fluid injection.
-
null (Ed.)The Nankai Trough Seismogenic Zone Experiment (NanTroSEIZE) is a coordinated, multiexpedition International Ocean Discovery Program (IODP) drilling project designed to investigate fault mechanics and seismogenesis along subduction megathrusts through direct sampling, in situ measurements, and long-term monitoring in conjunction with allied laboratory and numerical modeling studies. The fundamental scientific objectives of the NanTroSEIZE drilling project include characterizing the nature of fault slip and strain accumulation, fault and wall rock composition, fault architecture, and state variables throughout the active plate boundary system. IODP Expedition 365 is part of NanTroSEIZE Stage 3, with the following primary objectives: (1) retrieval of a temporary observatory at Site C0010 that has been monitoring temperature and pore pressure within the major splay thrust fault (termed the “megasplay”) at 400 meters below seafloor since November 2010 and (2) deployment of a complex long-term borehole monitoring system (LTBMS) that will be connected to the Dense Oceanfloor Network System for Earthquakes and Tsunamis (DONET) seafloor cabled observatory network postexpedition (anticipated June 2016). The LTBMS incorporates multilevel pore pressure sensing, a volumetric strainmeter, tiltmeter, geophone, broadband seismometer, accelerometer, and thermistor string. Together with an existing observatory at Integrated Ocean Drilling Program Site C0002 and a possible future installation near the trench, the Site C0010 observatory will allow monitoring within and above regions of contrasting behavior of the megasplay fault and the plate boundary as a whole. These include a site above the updip edge of the locked zone (Site C0002), a shallow site in the megasplay fault zone and its footwall (Site C0010), and a site at the tip of the accretionary prism (Integrated Ocean Drilling Program Site C0006). Together, this suite of observatories has the potential to capture deformation spanning a wide range of timescales (e.g., seismic and microseismic activity, slow slip, and interseismic strain accumulation) across a transect from near-trench to the seismogenic zone. Site C0010 is located 3.5 km along strike to the southwest of Integrated Ocean Drilling Program Site C0004. The site was drilled and cased during Integrated Ocean Drilling Program Expedition 319, with casing screens spanning a ~20 m interval that includes the megasplay fault, and suspended with a temporary instrument package (a “SmartPlug”). During Integrated Ocean Drilling Program Expedition 332 in late 2010, the instrument package was replaced with an upgraded sensor package (the “GeniusPlug”), which included pressure and temperature sensors and a set of geochemical and biological experiments. Expedition 365 achieved its primary scientific and operational objectives, including recovery of the GeniusPlug with a >5 y record of pressure and temperature conditions within the shallow megasplay fault zone, geochemical samples, and its in situ microbial colonization experiment; and installation of the LTBMS. The pressure records from the GeniusPlug include high-quality records of formation and seafloor responses to multiple fault slip events, including the 11 March 2011 Tohoku M9 and 1 April 2016 Mie-ken Nanto-oki M6 earthquakes. The geochemical sampling coils yielded in situ pore fluids from the splay fault zone, and microbes were successfully cultivated from the colonization unit. The complex sensor array, in combination with the multilevel hole completion, is one of the most ambitious and sophisticated observatory installations in scientific ocean drilling (similar to that in Hole C0002G, deployed in 2010). Overall, the installation went smoothly, efficiently, and ahead of schedule. The extra time afforded by the efficient observatory deployment was used for coring in Holes C0010B–C0010E. Despite challenging hole conditions, the depth interval corresponding to the screened casing across the megasplay fault was successfully sampled in Hole C0010C, and the footwall of the megasplay was sampled in Hole C0010E, with >50% recovery for both zones. In the hanging wall of the megasplay fault (Holes C0010C and C0010D), we recovered indurated silty clay with occasional ash layers and sedimentary breccias. Some of the deposits show burrows and zones of diagenetic alteration/colored patches. Mudstones show different degrees of deformation spanning from occasional fractures to intervals of densely fractured scaly claystones of up to >10 cm thickness. Sparse faulting with low displacement (usually <2 cm) is seen in core and exhibits primarily normal and, rarely, reversed sense of slip. When present, ash was entrained along fractures and faults. On one occasion, a ~10 cm thick ash layer was found, which showed a fining-downward gradation into a mottled zone with clasts of the underlying silty claystones. In Hole C0010E, the footwall to the megasplay fault was recovered. Sediments are horizontally to gently dipping and mainly comprise silt of olive-gray color. The deposits of the underthrust sediment prism are less indurated than the hanging wall mudstones and show lamination on a centimeter scale. The material is less intensely deformed than the mudstones, and apart from occasional fracturation (some of it being drilling disturbance), evidence of structural features is absent.more » « less
