An official website of the United States government
Heterogeneity in geometry, stress, and material properties is widely invoked to explain the observed spectrum of slow earthquake phenomena. However, the effects of length scale of heterogeneity on macroscopic fault sliding behavior remain underexplored. We investigate this question for subduction megathrusts, via linear stability analysis and quasi-dynamic simulations of slip on a dipping fault characterized by rate-and-state friction (RSF). Frictional heterogeneity is imposed through alternating velocity-strengthening (VS) and velocity-weakening (VW) patches, over length scales spanning from those representative of basement relief (several km) to the entrainment of contrasting lithologies (100s of m). The resulting fault behavior is controlled by: (1) the average frictional properties of the fault, and (2) the size of VW blocks relative to a critical length scale. Reasonable ranges of these properties yield sliding behaviors spanning from stable sliding, to slow and seismic slip events that are confined within VW blocks or propagate along the entire fault.
more »
« less
Not All Heterogeneity Is Equal: Length Scale of Frictional Property Variation as a Control on Subduction Megathrust Sliding Behavior
Abstract Heterogeneity in geometry, stress, and material properties is widely invoked to explain the observed spectrum of slow earthquake phenomena. However, the effects of length scale of heterogeneity on macroscopic fault sliding behavior remain underexplored. We investigate this question for subduction megathrusts, via linear stability analysis and quasi‐dynamic simulations of slip on a dipping fault characterized by rate‐and‐state friction. Frictional heterogeneity is imposed through alternating velocity‐strengthening and velocity‐weakening (VW) patches, over length scales spanning from those representative of basement relief (several km) to the entrainment of contrasting lithologies (100s of m). The resulting fault behavior is controlled by: (a) the average frictional properties of the fault, and (b) the size of VW blocks relative to a critical length scale. Reasonable ranges of these properties yield sliding behaviors spanning from stable sliding, to slow and seismic slip events that are confined within VW blocks or propagate along the entire fault.
more »
« less
- PAR ID:
- 10608287
- Publisher / Repository:
- DOI PREFIX: 10.1029
- Date Published:
- Journal Name:
- Geophysical Research Letters
- Volume:
- 52
- Issue:
- 8
- ISSN:
- 0094-8276
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
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 Many natural faults are believed to consist of velocity weakening (VW) patches surrounded by velocity strengthening (VS) sections. Numerical studies routinely employ this framework to study earthquake sequences including repeating earthquakes. In this laboratory study, we made a VW asperity, of lengthL, from a bare Poly(methyl methacrylate) PMMA frictional interface and coated the surrounding interface with Teflon to make VS fault sections. Behavior of this isolated asperity was studied as a function ofL(ranging from 100 to 400 mm) and the critical nucleation length, , which is inversely proportional to the applied normal stress (2–16 MPa). Consistent with recent numerical simulations, we observed aseismic slip for < 2, periodic slip for 2 < < 6, and non‐periodic slip for 10 < . Furthermore, we compared the experiments whereLwas contained by VS material to standard stick‐slip events whereLwas bounded by free surfaces (i.e.,L = the total sample length). The free surface case produced ∼10 times larger slip during stick‐slip events compared to the contained fault ruptures, even with identical . This disparity highlights how standard, complete‐rupture stick‐slip events differ from contained events expected in nature, due to both the free surface conditions and the heterogeneous normal stress along the fault near the free ends, as confirmed by Digital Image Correlation analysis. This study not only introduces the Teflon coating experimental technique for containing laboratory earthquake ruptures, but also highlights the utility of as a predictive parameter for earthquake behavior.more » « less
-
Abstract The empirical rate‐ and state‐dependent friction law is widely used to explain the frictional resistance of rocks. However, the constitutive parameters vary with temperature and sliding velocity, preventing extrapolation of laboratory results to natural conditions. Here, we explain the frictional properties of natural gouge from the San Andreas Fault, Alpine Fault, and the Nankai Trough from room temperature to ∼300°C for a wide range of slip‐rates with constant constitutive parameters by invoking the competition between two healing mechanisms with different thermodynamic properties. A transition from velocity‐strengthening to velocity‐weakening at steady‐state can be attained either by decreasing the slip‐rate or by increasing temperature. Our study provides a framework to understand the physics underlying the slip‐rate and state dependence of friction and the dependence of frictional properties on ambient physical conditions.more » « less
-
Whether Earth materials exhibit frictional creep or catastrophic failure is a crucial but unresolved problem in predicting landslide and earthquake hazards. Here, we show that field-scale observations of sliding velocity and pore water pressure at two creeping landslides are explained by velocity-strengthening friction, in close agreement with laboratory measurements on similar materials. This suggests that the rate-strengthening friction commonly measured in clay-rich materials may govern episodic slow slip in landslides, in addition to tectonic faults. Further, our results show more generally that transient slow slip can arise in velocity-strengthening materials from modulation of effective normal stress through pore pressure fluctuations. This challenges the idea that episodic slow slip requires a narrow range of transitional frictional properties near the stability threshold, or pore pressure feedbacks operating on initially unstable frictional slip.more » « less
