skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: Influence of Shear Heating and Thermomechanical Coupling on Earthquake Sequences and the Brittle‐Ductile Transition
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
Award ID(s):
1947448
PAR ID:
10367023
Author(s) / Creator(s):
 ;  
Publisher / Repository:
DOI PREFIX: 10.1029
Date Published:
Journal Name:
Journal of Geophysical Research: Solid Earth
Volume:
126
Issue:
6
ISSN:
2169-9313
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Deformation experiments on hematite characterize its slip‐rate dependent frictional properties and deformation mechanisms. These data inform interpretations of slip behavior from exhumed hematite‐coated faults and present‐day deformation at depth. We used a rotary‐shear apparatus to conduct single‐velocity and velocity‐step experiments on polycrystalline specular hematite rock (∼17 μm average plate thickness) at slip rates of 0.85 μm/s to 320 mm/s, displacements of primarily 1–3 cm and up to 45 cm, and normal stresses of 5 and 8.5 MPa. The average coefficient of friction is 0.70; velocity‐step experiments indicate velocity‐strengthening to velocity‐neutral behavior at rates <1 mm/s. Scanning electron microscopy showed experimentally generated faults develop in a semi‐continuous, thin layer of red hematite gouge. Angular gouge particles have an average diameter of ∼0.7 μm, and grain size reduction during slip yields a factor of 10–100 increase in surface area. Hematite is amenable to (U‐Th)/He thermochronometry, which can quantify fault‐related thermal and mechanical processes. Comparison of hematite (U‐Th)/He dates from the undeformed material and experimentally produced gouge indicates He loss occurs during comminution at slow deformation rates without an associated temperature rise required for diffusive loss. Our results imply that, in natural fault rocks, deformation localizes within coarse‐grained hematite by stable sliding, and that hematite (U‐Th)/He dates acquired from ultracataclasite or highly comminuted gouge reflect minor He loss unrelated to thermal processes. Consequently, the magnitude of temperature rise and associated thermal resetting in hematite‐bearing fault rocks based on (U‐Th)/He thermochronometry may be overestimated if only diffusive loss of He is considered. 
    more » « less
  2. 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
  3. Water ice Ih exhibits brittle behavior when rapidly loaded. Under tension, it fails via crack nucleation and propagation. Compressive failure is more complicated. Under low confinement, cracks slide and interact to form a frictional (Coulombic) fault. Under high confinement, frictional sliding is suppressed and adiabatic heating through crystallographic slip leads to the formation of a plastic fault. The coefficient of static friction increases with time under load, owing to creep of asperities in contact. The coefficient of kinetic (dynamic) friction, set by the ratio of asperity shear strength to hardness, increases with velocity at lower speeds and decreases at higher speeds as contacts melt through frictional heating. Microcracks, upon reaching a critical number density (which near the ductile-to-brittle transition is nearly constant above a certain strain rate), form a pathway for percolation. Additional work is needed on the effects of porosity and crack healing. ▪ An understanding of brittle failure is essential to better predict the integrity of the Arctic and Antarctic sea ice covers and the tectonic evolution of the icy crusts of Enceladus, Europa, and other extraterrestrial satellites. ▪ Fundamental to the brittle failure of ice is the initiation and propagation of microcracks, frictional sliding across crack faces, and localization of strain through both crack interaction and adiabatic heating. 
    more » « less
  4. Abundant heterogeneity has been documented on faults in nature across a wide range of length scales, including structural, mineralogical, and roughness variations. The role of complex heterogeneity on fault mechanics and frictional stability is not well established, and experiments investigating heterogeneity have typically incorporated a single source of heterogeneity. Here, we conduct rock friction experiments on rough, bimaterial faults that are creeping, or steadily sliding, to explore the role of lithological heterogeneity on fault mechanics and stability. When strong asperities juxtapose weak gouge, stable sliding occurs with a low friction coefficient, µ. Encounters of strong diabase asperities on talc gouge lined faults initiate dramatic increases in µ and transitions to unstable sliding characterized by frequent stick-slip events (StSE). Seismic moments and stress drops of StSE decrease with increasing asperity abundance. Stress is concentrated at asperities during encounters, increasing with decreasing asperity abundance and leading to extensive mechanical damage. Interactions between strong, velocity weakening asperities provide a model to explain the nucleation of seismic and aseismic slip events on nominally stable, creeping faults. 
    more » « less
  5. SUMMARY This paper examines the linear stability of sliding on faults embedded in a 2-D elastic medium that obey rate and state friction and have a finite length and/or are near a traction-free surface. Results are obtained using a numerical technique that allows for analysis of systems with geometrical complexity and heterogeneous material properties; however only systems with homogeneous frictional and material properties are examined. Some analytical results are also obtained for the special case of a fault that is parallel to a traction-free surface. For velocity-weakening faults with finite length, there is a critical fault length $$L^{*}$$ for unstable sliding that is analogous to the critical wavelength $$h^{*}$$ that is usually derived from infinite fault systems. Faults longer than $$L^{*}$$ are linearly unstable to perturbations of any length. On vertical strike-slip faults or faults in a full-space $$L^{*} \approx h^{*}/e$$, where e is Euler’s number. For dip-slip faults near a traction-free surface $$L^{*} \le h^{*}/e$$ and is a function of dip angle $$\beta$$, burial depth d of the fault’s up-dip edge and friction coefficient. In particular, $$L^{*}$$ is at least an order of magnitude smaller than $$h^{*}$$ on shallowly dipping ($$\beta < 10^\circ$$) faults that intersect the traction-free surface. Additionally, $$L^{*} \approx h^{*}/e$$ on dip-slip faults with burial depths $$d \ge h^{*}$$. For sliding systems that can be treated as a thin layer, such as landslides, glaciers or ice streams, $$L^{*} = h^{*}/2$$. Finally, conditions are established for unstable sliding on infinitely-long, velocity-strengthening faults that are parallel to a traction-free surface. 
    more » « less