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.
- Award ID(s):
- 1947107
- NSF-PAR ID:
- 10334531
- Date Published:
- Journal Name:
- Annual Review of Earth and Planetary Sciences
- Volume:
- 50
- Issue:
- 1
- ISSN:
- 0084-6597
- Page Range / eLocation ID:
- 323 to 343
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
Abstract -
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.
-
Abstract The icy shells of Enceladus and Europa consist from top down of cold (~100 K) ice at low (<0.1 MPa) pressure, cold ice at high pressure up to ~10 MPa, and warm ice near 273 K the base. The pressure ~10 MPa and temperature near 273 K of basal ice within Enceladus and Europa are similar to that within terrestrial glaciers that are known to have seismicity (icequakes). Warm ice easily melts during sliding so its icequakes are qualitatively explained and expected on these satellites if the macroscopic strain rates are comparable to those within terrestrial glaciers subjected to oceanic tides. However, cold ice at high pressures does not readily macroscopically melt during sliding. A dynamic weakening mechanism for crustal faults in rock may be applicable to Enceladus and Europa. Micron‐scale real contacts support ~0.4‐GPa shear tractions and normal tractions on rapidly sliding ice faults. At sliding velocities above ~0.1 m/s, the asperity tips of the contacts become hot and weak in ice. The macroscopic friction depends on the average strength of the asperity tips during the lifetimes of contact. The strength of the asperity tips self‐organizes so that frictional heating balances the heat lost from the asperity tip by conduction. The macroscopic coefficient of friction at coseismic sliding velocities decreases to a modest fraction of the low‐velocity coefficient of friction, but does not approach zero. This velocity‐weakening mechanism likely allows major icequakes within the cold interiors of Europa and Enceladus.
-
Modeling of rock inhomogeneity and anisotropy by explicit and implicit representation of microcracksFracture patterns experienced under a dynamic uniaxial compressive load are highly sensitive to rock microstructural defects due to its brittleness and the absence of macroscopic stress concentration points. We propose two different approaches for modeling rock microstructural defects and inhomogeneity. In the explicit realization approach, microcracks with certain statistics are incorporated in the computational domain. In the implicit realization approach, fracture strength values are sampled using a Weibull probability distribution. We use the Mohr-Coulomb failure criterion to define an effective stress in the context of an interfacial damage model. This model predicts crack propagation at angles ±ɸch = ±(45 − ɸ/2) relative to the direction of compressive load, where ɸ is the friction angle. By using appropriate models for fracture strength anisotropy, we demonstrate the interaction of rock weakest plane and ɸch. Numerical results demonstrate the greater effect of strength anisotropy on fracture pattern when an explicit approach is employed. In addition, the density of fractures increases as the angle of the weakest planes approaches ±ɸch. The fracture simulations are performed by an h-adaptive asynchronous spacetime discontinuous Galerkin (aSDG) method that can accommodate crack propagation in any directions.more » « less
-
Abstract 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.