The constitutive behavior of faults is central to many interconnected aspects of earthquake science, from fault dynamics to induced seismicity, to seismic hazards characterization. Yet, a friction law applicable to the range of temperatures found in the brittle crust and upper mantle is still missing. In particular, rocks often exhibit a transition from steady‐state velocity‐strengthening at room temperature to velocity‐weakening in warmer conditions that is poorly understood. Here, we investigate the effect of competing healing mechanisms on the evolution of frictional resistance in a physical model of rate‐, state‐, and temperature‐dependent friction. The yield strength for fault slip depends on the real area of contact, which is modulated by the competition between the growth and erosion of interfacial micro‐asperities. Incorporating multiple healing mechanisms and rock‐forming minerals with different thermodynamic properties allows a transition of the velocity‐ and temperature‐dependence of friction at steady‐state with varying temperatures. We explain the mechanical data for granite, pyroxene, amphibole, shale, and natural fault gouges with activation energies and stress power exponent for weakening of 10–50 kJ/mol and 55–150, respectively, compatible with subcritical crack growth and inter‐granular flow in the active slip zone. Activation energies for the time‐dependent healing process in the range 90–130 kJ/mol in dry conditions and 20–65 kJ/mol in wet conditions indicate the prominence of viscoelastic collapse of microasperities in the absence of water and of pressure‐solution creep, crack healing, and cementation when assisted by pore fluids.
- Award ID(s):
- 1848192
- NSF-PAR ID:
- 10324724
- Date Published:
- Journal Name:
- Earth, Planets and Space
- Volume:
- 73
- Issue:
- 1
- ISSN:
- 1880-5981
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
Abstract -
Abstract We performed deformation and grain growth experiments on natural olivine aggregates with olivine water contents (COH = 600 ± 300 H/106 Si) similar to upper mantle olivine, at 1000–1200°C and 1,400 ± 100 MPa confining pressure. Our experiments differ from published grain growth studies in that most were (1) conducted on natural olivine cores rather than hot‐pressed aggregates and (2) dynamically recrystallized prior to or during grain growth. We combine our results with similar experiments performed at 1200–1300°C and fit the data to a grain growth relationship, yielding a growth exponent (
p ) of 3.2, activation energy ( ) 620 ± 145 kJ mol−1(570 ± 145 kJ mol−1when accounting for the role of temperature on water content), activation volume (E G ) ~5V G × 10−6 m3mol−1, and rate constant (k 0) 1.8× 103 mp s−1. Our is within uncertainty of that predicted for dislocation creep of wet olivine (E G E* = 480 ± 40 kJ mol−1). Grain size in strain rate‐stepping samples adjusted to the olivine piezometer within 1.3–7.9% strain. The active grain boundary migration processes during deformation and dynamic recrystallization affect the kinetics of postdeformation grain growth, as grain boundary migration driven by strain energy density (ρ GBM) may delay the onset of grain growth driven by interfacial energy (γGBM). We compared our postdeformation grain growth rates with data from previously published hydrostatic annealing experiments on synthetic olivine. At geologic timescales, the growth rates are much slower than predicted by the existing wet olivine grain growth law. -
Abstract The rheology of the crust and mantle and the interaction of viscoelastic flow with seismic/aseismic slip on faults control the state of stress in the lithosphere over multiple seismic cycles. The rheological behavior of rocks is well constrained in a laboratory setting, but the
in situ properties of the lithosphere and its lateral variations remain poorly known. Here, we access the lower‐crustal rheology in Southern California by exploiting 8 years of geodetic postseismic deformation following the 2010 El Mayor‐Cucapah earthquake. The data illuminate viscoelastic flow in the lower crust with lateral variations of effective viscosity correlated with the geological province. We show that a Burgers assembly with dashpots following a nonlinear constitutive law can approximate the temporal evolution of stress and strain rate, indicating the activation of nonlinear transient creep before steady‐state dislocation creep. The transient and background viscosities in the lower crust of the Salton Trough are on the order of ~1018and ~1019 Pa s, respectively, about an order of magnitude lower than those in the surrounding regions. We highlight the importance of transient creep, nonlinear flow laws, and lateral variations of rheological properties to capture the entire history of postseismic relaxation following the El Mayor‐Cucapah earthquake. -
Abstract We use Sentinel‐1 and ALOS‐2 Interferometric Synthetic Aperture Radar (InSAR), and Global Navigation Satellite System (GNSS) data to investigate the mechanisms of coseismic and postseismic deformation due to the 2021 M7.4 Maduo (China) earthquake. We present a refined coseismic slip model constrained by the rupture trace and precisely located aftershocks. The InSAR time series corrected for the atmospheric and decorrelation noise reveal postseismic line of sight displacements up to ∼0.1 m. The displacements are discontinuous along the fault trace, indicating shallow afterslip and velocity‐strengthening friction in the top 2–3 km of the upper crust. The magnitude of shallow afterslip is however insufficient to compensate for the coseismic slip deficit, implying substantial off‐fault yielding. The observed surface deformation does not exhibit obvious features that could be attributed to poroelastic effects. We developed a fully coupled model that accounts for both stress‐driven creep on a deep localized shear zone and viscoelastic relaxation in the bulk of the lower crust. The mid‐ to near‐field data can be reasonably well explained by deep afterslip and/or non‐Maxwellian visco‐elasticity. Our results suggest a power‐law stress exponent of ∼4–4.5 assuming a power‐law rheology, and transient and steady‐state viscosities of 1018and 1019 Pa s, respectively, assuming a bi‐viscous (Burgers) rheology. However, a good fit to the GNSS data cannot be achieved assuming the bulk viscoelastic relaxation alone, and requires a contribution of deep afterlip and/or a localized shear zone extending through much of the lower crust.
-
Abstract Following the reanalysis of individual experimental runs of some widely cited studies (Jain et al., 2018,
https://doi.org/10.1002/2017JB014847 ), we revisit the global data analysis of Korenaga and Karato (2008,https://doi.org/10.1029/2007JB005100 ) with a significantly improved version of their Markov chain Monte Carlo inversion. Their algorithm, previously corrected by Mullet et al. () to minimize potential parameter bias, is further modified here to estimate more efficiently interrun biases in global data sets. Using the refined Markov chain Monte Carlo inversion technique, we simultaneously analyze experimental data on the deformation of olivine aggregates compiled from different studies. Realistic composite rheological models, including both diffusion and dislocation creep, are adopted, and the role of dislocation‐accommodated grain boundary sliding is also investigated. Furthermore, the influence of interrun biases on inversion results is studied using experimental and synthetic data. Our analysis shows that existing data can tightly constrain the grain‐size exponent for diffusion creep at ∼2, which is different from the value commonly assumed ( p = 3). Different data sets and model assumptions, however, yield nonoverlapping estimates on other flow‐law parameters, and the flow‐law parameters for grain boundary sliding are poorly resolved in most cases. We thus provide a few plausible candidate flow‐law models for olivine rheology to facilitate future geodynamic modeling. The availability of more data that explore a wider range of experimental conditions, especially higher pressures, is essential to improve our understanding of upper mantle rheology.