An official website of the United States government
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
Slide‐Hold‐Slide Protocols and Frictional Healing in Discrete Element Method (DEM) Simulations of Granular Fault Gouge
Abstract The empirical constitutive modeling framework of rate‐ and state‐dependent friction (RSF) is commonly used to describe the time‐dependent frictional response of fault gouge to perturbations from steady sliding. In a previous study (Ferdowsi & Rubin, 2020), we found that a granular‐physics‐based model of a fault shear zone, with time‐independent properties at the contact scale, reproduces the phenomenology of laboratory rock and gouge friction experiments in velocity‐step and slide‐hold (SH) protocols. A few slide‐hold‐slide (SHS) simulations further suggested that the granular model might outperform current empirical RSF laws in describing laboratory data. Here, we explore the behavior of the same Discrete Element Method (DEM) model in SH and SHS protocols over a wide range of sliding velocities, hold durations, and system stiffnesses, and provide additional support for this view. We find that, similar to laboratory data, the rate of stress decay during SH simulations is in general agreement with the “Slip law” version of the RSF equations, using parameter values determined independently from velocity step tests. During reslides following long hold times, the model, similar to lab data, produces a nearly constant rate of frictional healing with log hold time, with that rate being in the range of ∼0.5 to 1 times the RSF “state evolution” parameterb. We also find that, as in laboratory experiments, the granular layer undergoes log‐time compaction during holds. This is consistent with the traditional understanding of state evolution under the Aging law, even though the associated stress decay is similar to that predicted by the Slip and not the Aging law.
more »
« less
- Award ID(s):
- 1946434
- PAR ID:
- 10366924
- Publisher / Repository:
- DOI PREFIX: 10.1029
- Date Published:
- Journal Name:
- Journal of Geophysical Research: Solid Earth
- Volume:
- 126
- Issue:
- 12
- ISSN:
- 2169-9313
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Abstract Machine learning can predict the timing and magnitude of laboratory earthquakes using statistics of acoustic emissions. The evolution of acoustic energy is critical for lab earthquake prediction; however, the connections between acoustic energy and fault zone processes leading to failure are poorly understood. Here, we document in detail the temporal evolution of acoustic energy during the laboratory seismic cycle. We report on friction experiments for a range of shearing velocities, normal stresses, and granular particle sizes. Acoustic emission data are recorded continuously throughout shear using broadband piezo‐ceramic sensors. The coseismic acoustic energy release scales directly with stress drop and is consistent with concepts of frictional contact mechanics and time‐dependent fault healing. Experiments conducted with larger grains (10.5 μm) show that the temporal evolution of acoustic energy scales directly with fault slip rate. In particular, the acoustic energy is low when the fault is locked and increases to a maximum during coseismic failure. Data from traditional slide‐hold‐slide friction tests confirm that acoustic energy release is closely linked to fault slip rate. Furthermore, variations in the true contact area of fault zone particles play a key role in the generation of acoustic energy. Our data show that acoustic radiation is related primarily to breaking/sliding of frictional contact junctions, which suggests that machine learning‐based laboratory earthquake prediction derives from frictional weakening processes that begin very early in the seismic cycle and well before macroscopic failure.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
-
Nearly all frictional interfaces strengthen as the logarithm of time when sliding at ultra-low speeds. Observations of also logarithmic-in-time growth of interfacial contact area under such conditions have led to constitutive models that assume that this frictional strengthening results from purely time-dependent, and slip-insensitive, contact-area growth. The main laboratory support for such strengthening has traditionally been derived from increases in friction during “load-point hold” experiments, wherein a sliding interface is allowed to gradually self-relax down to subnanometric slip rates. In contrast, following step decreases in the shear loading rate, friction is widely reported to increase over a characteristic slip scale, independent of the magnitude of the slip-rate decrease—a signature of slip-dependent strengthening. To investigate this apparent contradiction, we subjected granite samples to a series of step decreases in shear rate of up to 3.5 orders of magnitude and load-point holds of up to 10,000 s, such that both protocols accessed the phenomenological regime traditionally inferred to demonstrate time-dependent frictional strengthening. When modeling the resultant data, which probe interfacial slip rates ranging from 3 . μ m · s − 1 . to less than 10 − 5 μ m · s − 1 , we found that constitutive models where low slip-rate friction evolution mimics log-time contact-area growth require parameters that differ by orders of magnitude across the different experiments. In contrast, an alternative constitutive model, in which friction evolves only with interfacial slip, fits most of the data well with nearly identical parameters. This leads to the surprising conclusion that frictional strengthening is dominantly slip-dependent, even at subnanometric slip rates.more » « less
-
Abstract The southern San Andreas fault is in its interseismic period, occasionally releasing some stored elastic strain during triggered slow slip events (SSEs) at <2.5 km depth. A distinct, shallowly exhumed gouge defines the fault and is present at SSE depths. To evaluate if this material can host SSEs, we characterize its mineralogy, microstructures, and frictional behavior with water‐saturated deformation experiments near‐in situ conditions, and we compare laboratory healing rates to natural SSEs. Our results show that slip localizes along clay surfaces in both laboratory and natural settings. The gouge is weak (coefficient of friction of ∼0.29), exhibits low healing rates (<0.001/decade), and transitions from unstable to stable behavior at slip rates above ∼1 μm/s. Healing rate and friction drop data from laboratory instabilities are comparable to geodetically‐constrained values for SSEs. Collective observations indicate this gouge could host shallow SSEs and/or localize slip facilitating dynamic rupture propagation to the surface.more » « less
