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   
                    
                            
                            The Role of Quartz Cementation in the Seismic Cycle: A Critical Review
                        
                    
    
            Abstract Because quartz veins are common in fault zones exhumed from earthquake nucleation temperatures (150°C–350°C), quartz cementation may be an important mechanism of strength recovery between earthquakes. This interpretation requires that cementation occurs within a single interseismic period. We review slip‐related processes that have been argued to allow rapid quartz precipitation in faults, including: advection of silica‐saturated fluids, coseismic pore‐fluid pressure drops, frictional heating, dissolution‐precipitation creep, precipitation of amorphous phases, and variations in fluid and mineral‐surface chemistry. We assess the rate and magnitude of quartz growth that may result from each of the examined mechanisms. We find limitations to the kinetics and mass balance of silica precipitation that emphasize two end‐member regimes. First, the mechanisms we explore, given current kinetic constraints, cannot explain mesoscale fault‐fracture vein networks developing, even incrementally, on interseismic timescales. On the other hand, some mechanisms appear capable, isolated or in combination, of cementing micrometer‐to‐millimeter thick principal slip surfaces in days to years. This does not explain extensive vein networks in fault damage zones, but allows the involvement of quartz cements in fault healing. These end‐members lead us to hypothesize that high flux scenarios, although more important for voluminous hydrothermal mineralization, may be of subsidiary importance to local, diffusive mass transport in low fluid‐flux faults when discussing the mechanical implications of quartz cements. A renewed emphasis on the controls on quartz cementation rates in fault zones will, however, be integral to developing a more complete understanding of strength recovery following earthquake rupture. 
        more » 
        « less   
        
    
                            - Award ID(s):
- 1951985
- PAR ID:
- 10446094
- Publisher / Repository:
- DOI PREFIX: 10.1029
- Date Published:
- Journal Name:
- Reviews of Geophysics
- Volume:
- 60
- Issue:
- 1
- ISSN:
- 8755-1209
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
- 
            
- 
            Abstract Progressive cementation and sealing of fault-localized fractures impact crustal mass transport and the recovery of fault strength following slip events. We use discrete fracture network (DFN) models to examine how fracture sealing during end-member cementation mechanisms (i.e., reaction- versus transported-limited cementation) influences the partitioning of fluid flow through time. DfnWorks was used to generate randomized fracture networks parameterized with fracture orientation data compiled from field studies. Single-phase flow simulations were performed for each network over a series of timesteps, and network parameters were modified to reflect progressive fracture sealing consistent with either reaction- or transport-limited crystal growth. Results show that when fracture cementation is reaction-limited, fluid flow becomes progressively channelized into a smaller number of fractures with larger apertures. When fracture cementation is transport-limited, fluid flow experiences progressive dechannelization, becoming more homogeneously distributed throughout the fracture network. These behaviors are observed regardless of the DFN parameterization, suggesting that the effect is an intrinsic component of all fracture networks subjected to the end-member cementation mechanisms. These results have first-order implications for the spatial distribution of fluid flow in fractured rocks and recovery of permeability and strength during fault/fracture healing in the immediate aftermath of fault slip.more » « less
- 
            Abstract Megathrust geometric properties exhibit some of the strongest correlations with maximum earthquake magnitude in global surveys of large subduction zone earthquakes, but the mechanisms through which fault geometry influences subduction earthquake cycle dynamics remain unresolved. Here, we develop 39 models of sequences of earthquakes and aseismic slip (SEAS) on variably‐dipping planar and variably‐curved nonplanar megathrusts using the volumetric, high‐order accurate codetandemto account for fault curvature. We vary the dip, downdip curvature and width of the seismogenic zone to examine how slab geometry mechanically influences megathrust seismic cycles, including the size, variability, and interevent timing of earthquakes. Dip and curvature control characteristic slip styles primarily through their influence on seismogenic zone width: wider seismogenic zones allow shallowly‐dipping megathrusts to host larger earthquakes than steeply‐dipping ones. Under elevated pore pressure and less strongly velocity‐weakening friction, all modeled fault geometries host uniform periodic ruptures. In contrast, shallowly‐dipping and sharply‐curved megathrusts host multi‐period supercycles of slow‐to‐fast, small‐to‐large slip events under higher effective stresses and more strongly velocity‐weakening friction. We discuss how subduction zones' maximum earthquake magnitudes may be primarily controlled by the dip and dimensions of the seismogenic zone, while second‐order effects from structurally‐derived mechanical heterogeneity modulate the recurrence frequency and timing of these events. Our results suggest that enhanced co‐ and interseismic strength and stress variability along the megathrust, such as induced near areas of high or heterogeneous fault curvature, limits how frequently large ruptures occur and may explain curved faults' tendency to host more frequent, smaller earthquakes than flat 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.more » « less
- 
            Faults are usually surrounded by damage zones associated with localized deformation. Here we use fully dynamic earthquake cycle simulations to quantify the behaviors of earthquakes in fault damage zones. We show that fault damage zones can make a significant contribution to the spatial and temporal seismicity distribution. Fault stress heterogeneities generated by fault zone waves persist over multiple earthquake cycles that, in turn, produce small earthquakes that are absent in homogeneous simulations with the same friction conditions. Shallow fault zones can produce a bimodal depth distribution of earthquakes with clustering of seismicity at both shallower and deeper depths. Fault zone healing during the interseismic period also promotes the penetration of aseismic slip into the locked region and reduces the sizes of fault asperities that host earthquakes. Hence, small and moderate subsurface earthquakes with irregular recurrence intervals are commonly observed in immature fault zone simulations with interseismic healing. To link our simulation results to geological observations, we will use simulated fault slip at different depths to infer the timing and recurrence intervals of earthquakes and discuss how such measurements can affect our understanding of earthquake behaviors. We will also show that the maturity and material properties of fault damage zones have strong influence on whether long-term earthquake characteristics are represented by single events.more » « less
 An official website of the United States government
An official website of the United States government 
				
			 
					 
					
