The evolution of rock friction is more sensitive to slip than elapsed time, even at near-zero slip rates

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 .[Formula: see text]. to less than [Formula: see text], 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.

Source time functions of earthquakes based on a stochastic differential equation

Source time functions are essential observable quantities in seismology; they have been investigated via kinematic inversion analyses and compiled into databases. Given the numerous available results, some empirical laws on source time functions have been established, even though they are complicated and fluctuated time series. Theoretically, stochastic differential equations, including a random variable and white noise, are suitable for modeling complicated phenomena. In this study, we model source time functions as the convolution of two stochastic processes (known as Bessel processes). We mathematically and numerically demonstrate that this convolution satisfies some of the empirical laws of source time functions, including non-negativity, finite duration, unimodality, a growth rate proportional to [Formula: see text], [Formula: see text]-type spectra, and frequency distribution (i.e., the Gutenberg-Richter law). We interpret this convolution and speculate that the stress drop rate and fault impedance follow the same Bessel process.

Machine Learning Predicts the Timing and Shear Stress Evolution of Lab Earthquakes Using Active Seismic Monitoring of Fault Zone Processes

Machine learning (ML) techniques have become increasingly important in seismology and earthquake science. Lab-based studies have used acoustic emission data to predict time-to-failure and stress state, and in a few cases, the same approach has been used for field data. However, the underlying physical mechanisms that allow lab earthquake prediction and seismic forecasting remain poorly resolved. Here, we address this knowledge gap by coupling active-source seismic data, which probe asperity-scale processes, with ML methods. We show that elastic waves passing through the lab fault zone contain information that can predict the full spectrum of labquakes from slow slip instabilities to highly aperiodic events. The ML methods utilize systematic changes in P-wave amplitude and velocity to accurately predict the timing and shear stress during labquakes. The ML predictions improve in accuracy closer to fault failure, demonstrating that the predictive power of the ultrasonic signals improves as the fault approaches failure. Our results demonstrate that the relationship between the ultrasonic parameters and fault slip rate, and in turn, the systematically evolving real area of contact and asperity stiffness allow the gradient boosting algorithm to "learn" about the state of the fault and its proximity to failure. Broadly, our results demonstrate the utility of physics-informed ML in forecasting the imminence of fault slip at the laboratory scale, which may have important implications for earthquake mechanics in nature.

Unconventional singularities and energy balance in frictional rupture

A widespread framework for understanding frictional rupture, such as earthquakes along geological faults, invokes an analogy to ordinary cracks. A distinct feature of ordinary cracks is that their near edge fields are characterized by a square root singularity, which is intimately related to the existence of strict dissipation-related lengthscale separation and edge-localized energy balance. Yet, the interrelations between the singularity order, lengthscale separation and edge-localized energy balance in frictional rupture are not fully understood, even in physical situations in which the conventional square root singularity remains approximately valid. Here we develop a macroscopic theory that shows that the generic rate-dependent nature of friction leads to deviations from the conventional singularity, and that even if this deviation is small, significant non-edge-localized rupture-related dissipation emerges. The physical origin of the latter, which is predicted to vanish identically in the crack analogy, is the breakdown of scale separation that leads an accumulated spatially-extended dissipation, involving macroscopic scales. The non-edge-localized rupture-related dissipation is also predicted to be position dependent. The theoretical predictions are quantitatively supported by available numerical results, and their possible implications for earthquake physics are discussed.

Ordinary cracks in bulk materials feature square root singular deformation fields near their edge. Here, the authors show that rupture fronts propagating along frictional interfaces, while resembling ordinary cracks in some respects, feature edge sigularity that differs from the conventional square root one.

Shear Controls Frictional Aging by Erasing Memory

We simultaneously measure the static friction and the real area of contact between two solid bodies. These quantities are traditionally considered equivalent, and under static conditions both increase logarithmically in time, a phenomenon coined aging. Here we show that the frictional aging rate is determined by the combination of the aging rate of the real area of contact and two memory-erasure effects that occur when shear is changed (e.g., to measure static friction.) The application of a static shear load accelerates frictional aging while the aging rate of the real area of contact is unaffected. Moreover, a negative static shear-pulling instead of pushing-slows frictional aging, but similarly does not affect the aging of contacts. The origin of this shear effect on aging is geometrical. When shear load is increased, minute relative tilts between the two blocks prematurely erase interfacial memory prior to sliding, negating the effect of aging. Modifying the loading point of the interface eliminates these tilts and as a result frictional aging rate becomes insensitive to shear. We also identify a secondary memory-erasure effect that remains even when all tilts are eliminated and show that this effect can be leveraged to accelerate aging by cycling between two static shear loads.

Thermo-Poroelastic Analysis of Induced Seismicity at the Basel Enhanced Geothermal System

Geothermal energy has emerged as an alternative to ensure a green energy supply while tackling climate change. Geothermal systems extract the heat stored in the Earth’s crust by warming up water, but the low rock permeability at exploitation depths may require the hydraulic stimulation of the rock fracture network. Enhanced Geothermal Systems (EGS) employ techniques such as hydro-shearing and hydro-fracturing for that purpose, but their use promotes anthropogenic earthquakes induced by the injection or extraction of fluids. This work addresses this problem through developing a computational 3D model to explore fault reactivation and evaluating the potential for earthquake triggering at preexisting geological faults. These are included in the model as frictional contacts that allow the relative displacement between both of its sides, governed by rate-and-state friction laws and fully coupled with thermo-hydro-mechanical equations. We apply our methodology to the Basel project, employing the on-site parameters and conditions. Our results demonstrate that earthquakes which occurred in December 2006 in Basel (Switzerland) are compatible with the geomechanical and frictional consequences of the hydraulic stimulation of the rock mass. The application of our model also shows that it can be useful for predicting fault reactivation and engineering injection protocols for managing the safe and sustainable operation of EGS.

Spatiotemporal Dynamics of Frictional Systems: The Interplay of Interfacial Friction and Bulk Elasticity

Frictional interfaces are abundant in natural and engineering systems, and predicting their behavior still poses challenges of prime scientific and technological importance. At the heart of these challenges lies the inherent coupling between the interfacial constitutive relation—the macroscopic friction law—and the bulk elasticity of the bodies that form the frictional interface. In this feature paper, we discuss the generic properties of a minimal macroscopic friction law and the many ways in which its coupling to bulk elasticity gives rise to rich spatiotemporal frictional dynamics. We first present the widely used rate-and-state friction constitutive framework, discuss its power and limitations, and propose extensions that are supported by experimental data. We then discuss how bulk elasticity couples different parts of the interface, and how the range and nature of this interaction are affected by the system’s geometry. Finally, in light of the coupling between interfacial and bulk physics, we discuss basic phenomena in spatially extended frictional systems, including the stability of homogeneous sliding, the onset of sliding motion and a wide variety of propagating frictional modes (e.g., rupture fronts, healing fronts and slip pulses). Overall, the results presented and discussed in this feature paper highlight the inseparable roles played by interfacial and bulk physics in spatially extended frictional systems.

Evolution of real contact area under shear and the value of static friction of soft materials

The frictional properties of a rough contact interface are controlled by its area of real contact, the dynamical variations of which underlie our modern understanding of the ubiquitous rate-and-state friction law. In particular, the real contact area is proportional to the normal load, slowly increases at rest through aging, and drops at slip inception. Here, through direct measurements on various contacts involving elastomers or human fingertips, we show that the real contact area also decreases under shear, with reductions as large as 30[Formula: see text], starting well before macroscopic sliding. All data are captured by a single reduction law enabling excellent predictions of the static friction force. In elastomers, the area-reduction rate of individual contacts obeys a scaling law valid from micrometer-sized junctions in rough contacts to millimeter-sized smooth sphere/plane contacts. For the class of soft materials used here, our results should motivate first-order improvements of current contact mechanics models and prompt reinterpretation of the rate-and-state parameters.

Stable and unstable development of an interfacial sliding instability

Examining a nonlinear instability of sliding rate on a frictional interface of elastic bodies, we investigate whether laboratory-constrained frictional relations suggest universal scaling under even the simplest of configurations. We find blowup solutions by solving an equivalent, classical problem in fracture mechanics. The solutions are fixed points of a dynamical system and we show that their stability is lost by a cascade of Hopf bifurcations as a single problem parameter is increased, leading to chaotic dynamics.

Modeling fast and slow earthquakes at various scales

Earthquake sources represent dynamic rupture within rocky materials at depth and often can be modeled as propagating shear slip controlled by friction laws. These laws provide boundary conditions on fault planes embedded in elastic media. Recent developments in observation networks, laboratory experiments, and methods of data analysis have expanded our knowledge of the physics of earthquakes. Newly discovered slow earthquakes are qualitatively different phenomena from ordinary fast earthquakes and provide independent information on slow deformation at depth. Many numerical simulations have been carried out to model both fast and slow earthquakes, but problems remain, especially with scaling laws. Some mechanisms are required to explain the power-law nature of earthquake rupture and the lack of characteristic length. Conceptual models that include a hierarchical structure over a wide range of scales would be helpful for characterizing diverse behavior in different seismic regions and for improving probabilistic forecasts of earthquakes.

Instabilities at frictional interfaces: creep patches, nucleation, and rupture fronts

The strength and stability of frictional interfaces, ranging from tribological systems to earthquake faults, are intimately related to the underlying spatially extended dynamics. Here we provide a comprehensive theoretical account, both analytic and numeric, of spatiotemporal interfacial dynamics in a realistic rate-and-state friction model, featuring both velocity-weakening and velocity-strengthening behaviors. Slowly extending, loading-rate-dependent creep patches undergo a linear instability at a critical nucleation size, which is nearly independent of interfacial history, initial stress conditions, and velocity-strengthening friction. Nonlinear propagating rupture fronts-the outcome of instability-depend sensitively on the stress state and velocity-strengthening friction. Rupture fronts span a wide range of propagation velocities and are related to steady-state-front solutions.