More Like this

1. Constitutive Behavior of Rocks During the Seismic Cycle

   https://doi.org/10.1029/2023AV000972  

   Barbot, Sylvain (September 2023, AGU Advances)

   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
2. Influence of Shear Heating and Thermomechanical Coupling on Earthquake Sequences and the Brittle‐Ductile Transition

   https://doi.org/10.1029/2020JB021394  

   Allison, Kali L.; Dunham, Eric M. (June 2021, Journal of Geophysical Research: Solid Earth)

   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
3. Acoustic Energy Release During the Laboratory Seismic Cycle: Insights on Laboratory Earthquake Precursors and Prediction

   https://doi.org/10.1029/2019JB018975  

   Bolton, David C.; Shreedharan, Srisharan; Rivière, Jacques; Marone, Chris (August 2020, Journal of Geophysical Research: Solid Earth)

   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
4. Strong asperities nucleate earthquakes on laboratory faults

   https://doi.org/10.1130/g52853.1  

   Barbery, Monica; Hirth, Greg; Tullis, Terry (February 2025, Geology)

   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. Modeling of rock inhomogeneity and anisotropy by explicit and implicit representation of microcracks

   Abedi, R.; Clarke, P.L. (January 2018, Proceeding 52th U.S. Rock Mechanics/Geomechanics Symposium)

   Fracture 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