Abstract Establishing a constitutive law for fault friction is a crucial objective of earthquake science. However, the complex frictional behavior of natural and synthetic gouges in laboratory experiments eludes explanations. Here, we present a constitutive framework that elucidates the rate, state, and temperature dependence of fault friction under the relevant sliding velocities and temperatures of the brittle lithosphere during seismic cycles. The competition between healing mechanisms, such as viscoelastic collapse, pressure‐solution creep, and crack sealing, explains the low‐temperature stability transition from steady‐state velocity‐strengthening to velocity‐weakening as a function of slip‐rate and temperature. In addition, capturing the transition from cataclastic flow to semi‐brittle creep accounts for the stabilization of fault slip at elevated temperatures. We calibrate the model using extensive laboratory data on synthetic albite and granite gouge, and on natural samples from the Alpine Fault and the Mugi Mélange in the Shimanto accretionary complex in Japan. The constitutive model consistently explains the evolving frictional response of fault gouge from room temperature to 600°C for sliding velocities ranging from nanometers to millimeters per second. The frictional response of faults can be uniquely determined by the in situ lithology and the prevailing hydrothermal conditions.
more »
« less
Velocity and Temperature Dependence of Steady‐State Friction of Natural Gouge Controlled by Competing Healing Mechanisms
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
- Award ID(s):
- 1848192
- PAR ID:
- 10525041
- Publisher / Repository:
- American Geophysical Union
- Date Published:
- Journal Name:
- Geophysical Research Letters
- Volume:
- 51
- Issue:
- 11
- ISSN:
- 0094-8276
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
The frictional properties of faults control the initiation and propagation of earthquakes and the associated hazards. Although the ambient temperature and instantaneous slip velocity controls on friction in isobaric conditions are increasingly well understood, the role of normal stress on steady‐state and transient frictional behaviors remains elusive. The friction coefficient of rocks exhibits a strong dependence on normal stress at typical crustal depths. Furthermore, rapid changes in normal stress cause a direct effect on friction followed by an evolutionary response. Here, we derive a constitutive friction law that consistently explains the yield strength of rocks from atmospheric pressure to gigapascals while capturing the transient behavior following perturbations in normal stress. The model explains the frictional strength of a variety of sedimentary, metamorphic, and igneous rocks and the slip‐dependent response upon normal stress steps of Westerly granite bare contact and synthetic gouges made of quartz and a mixture of quartz and smectite. The nonlinear normal stress dependence of the frictional resistance may originate from the distribution of asperities that control the real area of contact. The direct and transient effects may be important for induced seismicity by hydraulic fracturing or for naturally occurring normal stress perturbations within fault zones in the brittle crust.more » « less
-
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
-
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
-
Avouac, J-P (Ed.)The Pelona–Orocopia–Rand (POR) schists were emplaced during the Farallon flat subduction in the early Cenozoic and now occupy the root of major strike-slip faults of the San Andreas Fault system. The POR schists are considered frictionally stable at lower temperatures than other basement rocks, limiting the maximum depth of seismicity in Southern California. However, experimental constraints on the composition and frictional properties of POR schists are still missing. Here, we study the frictional behavior of synthetic gouge derived from Pelona, Portal, and Rand Mountain schist wall rocks under hydrothermal, triaxial conditions. We conduct velocity-step experiments from 0.04 to 1 μm/s from room temperature to 500ºC under 200 MPa effective normal stress, including a 30 MPa porefluid pressure. The frictional stability of POR schists in the lower crust is caused by a thermally activated transition from slip-rate- and state-dependent friction to inherently stable, rate-dependent creep between 300ºC and 500ºC, depending on sample composition and slip-rate. The mineralogy of POR schists shows much variability caused by different protoliths and metamorphic grades, featuring various amounts of phyllosilicates, quartz, feldspar, and amphibole. Pelona and Portal schists exhibit a velocity-weakening regime enabling the nucleation and propagation of earthquakes when exhumed in the middle crust, as in the Mojave section of the San Andreas Fault. The contrasted frictional properties of POR schists exemplify the lithological control of seismic processes and associated hazards.more » « less
An official website of the United States government

