Ultralow frictional healing explains recurring slow slip events.

Plate motion on shallow subduction megathrusts is accommodated by a spectrum of tectonic slip modes. However, the frictional properties and conditions that sustain these diverse slip behaviors remain enigmatic. Frictional healing is one such property, which describes the degree of fault restrengthening between earthquakes. We show that the frictional healing rate of materials entrained along the megathrust at the northern Hikurangi margin, which hosts well-characterized recurring shallow slow slip events (SSEs), is nearly zero (<0.0001 per decade). These low healing rates provide a mechanism for the low stress drops (<50 kilopascals) and short recurrence times (1 to 2 years) characteristic of shallow SSEs at Hikurangi and other subduction margins. We suggest that near-zero frictional healing rates, associated with weak phyllosilicates that are common in subduction zones, may promote frequent, small-stress-drop, slow ruptures near the trench.

Creep fronts and complexity in laboratory earthquake sequences illuminate delayed earthquake triggering.

Earthquakes occur in clusters or sequences that arise from complex triggering mechanisms, but direct measurement of the slow subsurface slip responsible for delayed triggering is rarely possible. We investigate the origins of complexity and its relationship to heterogeneity using an experimental fault with two dominant seismic asperities. The fault is composed of quartz powder, a material common to natural faults, sandwiched between 760 mm long polymer blocks that deform the way 10 meters of rock would behave. We observe periodic repeating earthquakes that transition into aperiodic and complex sequences of fast and slow events. Neighboring earthquakes communicate via migrating slow slip, which resembles creep fronts observed in numerical simulations and on tectonic faults. Utilizing both local stress measurements and numerical simulations, we observe that the speed and strength of creep fronts are highly sensitive to fault stress levels left behind by previous earthquakes, and may serve as on-fault stress meters.

The High-Frequency Signature of Slow and Fast Laboratory Earthquakes.

Tectonic faults fail through a spectrum of slip modes, ranging from slow aseismic creep to rapid slip during earthquakes. Understanding the seismic radiation emitted during these slip modes is key for advancing earthquake science and earthquake hazard assessment. In this work, we use laboratory friction experiments instrumented with ultrasonic sensors to document the seismic radiation properties of slow and fast laboratory earthquakes. Stick-slip experiments were conducted at a constant loading rate of 8 µm/s and the normal stress was systematically increased from 7 to 15 MPa. We produced a full spectrum of slip modes by modulating the loading stiffness in tandem with the fault zone normal stress. Acoustic emission data were recorded continuously at 5 MHz. We demonstrate that the full continuum of slip modes radiate measurable high-frequency energy between 100 and 500 kHz, including the slowest events that have peak fault slip rates <100 µm/s. The peak amplitude of the high-frequency time-domain signals scales systematically with fault slip velocity. Stable sliding experiments further support the connection between fault slip rate and high-frequency radiation. Experiments demonstrate that the origin of the high-frequency energy is fundamentally linked to changes in fault slip rate, shear strain, and breaking of contact junctions within the fault gouge. Our results suggest that having measurements close to the fault zone may be key for documenting seismic radiation properties and fully understanding the connection between different slip modes.

Frequency-Magnitude Statistics of Laboratory Foreshocks Vary With Shear Velocity, Fault Slip Rate, and Shear Stress.

Understanding the temporal evolution of foreshocks and their relation to earthquake nucleation is important for earthquake early warning systems, earthquake hazard assessment, and earthquake physics. Laboratory experiments on intact rock and rough fractures have demonstrated that the number and size of acoustic emission (AE) events increase and that the Gutenberg-Richter b-value decreases prior to coseismic failure. However, for lab fault zones of finite width, where shear occurs within gouge, the physical processes that dictate temporal variations in frequency-magnitude (F/M) statistics of lab foreshocks are unclear. Here, we report on a series of laboratory experiments to illuminate the physical processes that govern temporal variations in b-value and AE size. We record AE data continuously for hundreds of lab seismic cycles and report F/M statistics. Our foreshock catalogs include cases where F/M data are not exponentially distributed, but we retain the concept of b-value for comparison with other works. We find that b-value decreases as the fault approaches failure, consistent with previous works. We also find that b-value scales inversely with shear velocity and fault slip rate, suggesting that fault slip acceleration during earthquake nucleation could impact foreshock F/M statistics. We propose that fault zone dilation and grain mobilization have a strong influence on foreshock magnitude. Fault dilation at higher shearing rates increases porosity and results in larger foreshocks and smaller b-values. Our observations suggest that lab earthquakes are preceded by a preparatory nucleation phase with systematic variations in AE and fault zone properties.

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.

Acoustic Energy Release During the Laboratory Seismic Cycle: Insights on Laboratory Earthquake Precursors and Prediction.

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.

Slow slip source characterized by lithological and geometric heterogeneity.

Slow slip events (SSEs) accommodate a significant proportion of tectonic plate motion at subduction zones, yet little is known about the faults that actually host them. The shallow depth (<2 km) of well-documented SSEs at the Hikurangi subduction zone offshore New Zealand offers a unique opportunity to link geophysical imaging of the subduction zone with direct access to incoming material that represents the megathrust fault rocks hosting slow slip. Two recent International Ocean Discovery Program Expeditions sampled this incoming material before it is entrained immediately down-dip along the shallow plate interface. Drilling results, tied to regional seismic reflection images, reveal heterogeneous lithologies with highly variable physical properties entering the SSE source region. These observations suggest that SSEs and associated slow earthquake phenomena are promoted by lithological, mechanical, and frictional heterogeneity within the fault zone, enhanced by geometric complexity associated with subduction of rough crust.