Frictional state evolution during normal stress perturbations probed with ultrasonic waves

Fault normal stress changes dynamically during earthquake rupture; however, the impact of these changes on dynamic frictional strength is poorly understood. Here we report on a laboratory study to investigate the effect of normal stress perturbations on the friction of westerly granite surfaces sheared under normal stresses of 1-25 MPa. We measure changes in surface friction and elastic properties, using acoustic waves, for step changes in normal stress of 1–50% and shearing velocities of 1-100 μm/s. We demonstrate that transmitted elastic wave amplitude is a reliable proxy for the real contact area at the fault interface at steady state. For step increases in normal stress, wave amplitude increases immediately and then continues to increase during elastic shear loading to a peak value from which it decreases as fault slip rate increases. Friction changes in a similar fashion, showing an inelastic increase over a characteristic shear displacement that is independent of loading rate. Perturbations in normal stress during shear cause excursions in the frictional slip rate that must be accounted for in order to accurately predict the evolution of fault strength and elastic properties. Our work improves understanding of induced seismicity and triggered earthquakes with particular focus on simulating static triggering and stress transfer phenomena using rate-and-state frictional formulations in earthquake rupture models

Frictional and lithological controls on shallow slow slip at the Northern Hikurangi Margin

Slow slip events (SSEs) have been identified at subduction zones globally as an important link in the continuum between elastodynamic ruptures and stable creep. The northern Hikurangi margin is home to shallow SSEs which propagate to within 2 km of the seafloor and possibly to the trench, providing insights into the physical conditions conducive to SSE behavior. We report on a suite of friction experiments performed on protolith material entering the SSE source region at the Hikurangi margin, collected during the International Ocean Discovery Program Expedition 375. We performed velocity stepping and slide-hold-slide experiments over a range of fault slip rates, from plate rate (5 cm/yr or 1.6 × 10−9 m/s) to ∼1 mm/s (10−3 m/s) and quantified the frictional velocity dependence and healing rates for a range of lithologies atdifferent stresses. The frictional velocity dependence (a-b) and critical slip distance DC increase with fault slip rate in our experiments. We observe atransition from velocity weakening to strengthening at slip rates of ∼0.3 µm/s. This velocity dependence of DC could be due to a combination of dilatant strengthening and a widening of the active shear zone at higher slip rates. We document low healing rates in the clay-rich volcaniclastic conglomerates, which lie above the incoming plate basement at least locally, and relatively higher healing rates in the chalk lithology. Finally, our experimental constraints on healing rates in different input lithologies extrapolated to timescales of 1–10 years are consistent with the geodetically inferred low stress drops and healing rates characteristic of the Hikurangi SSEs

Deep learning can predict laboratory quakes from active source seismic data

Small changes in seismic wave properties foretell frictional failure in laboratory experiments and in some cases on seismic faults. Such precursors include systematic changes in wave velocity and amplitude throughout the seismic cycle. However, the relationships between wave features and shear stress are complex. Here, we use data from lab friction experiments that include continuous measurement of elastic waves traversing the fault and build data-driven models to learn these complex relations. We demonstrate that deep learning models accurately predict the timing and size of laboratory earthquakes based on wave features. Additionally, the transportability of models is explored by using data from different experiments. Our deep learning models transfer well to unseen datasets providing high-fidelity models with much less training. These prediction methods can be potentially applied in the field for earthquake early warning in conjunction with long-term time-lapse seismic monitoring of crustal faults, CO2 storage sites and unconventional energy reservoirs

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

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

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 (&lt;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

Expedition 372B/375 summary

Slow slip events (SSEs) at the northern Hikurangi subduction margin, New Zealand, are among the best-documented shallow SSEs on Earth. International Ocean Discovery Program Expeditions 372 and 375 were undertaken to investigate the processes and in situ conditions that underlie subduction zone SSEs at the northern Hikurangi Trough. We accomplished this goal by (1) coring and geophysical logging at four sites, including penetration of an active thrust fault (the Pāpaku fault) near the deformation front, the upper plate above the SSE source region, and the incoming sedimentary succession in the Hikurangi Trough and atop the Tūranganui Knoll seamount; and (2) installing borehole observatories in the Pāpaku fault and in the upper plate overlying the slow slip source region. Logging-while-drilling (LWD) data for this project were acquired as part of Expedition 372, and coring, wireline logging, and observatory installations were conducted during Expedition 375. Northern Hikurangi subduction margin SSEs recur every 1–2 y and thus provide an ideal opportunity to monitor deformation and associated changes in chemical and physical properties throughout the slow slip cycle. In situ measurements and sampling of material from the sedimentary section and oceanic basement of the subducting plate reveal the rock properties, composition, lithology, and structural character of material that is transported downdip into the SSE source region. A recent seafloor geodetic experiment raises the possibility that SSEs at northern Hikurangi may propagate to the trench, indicating that the shallow thrust fault (the Pāpaku fault) targeted during Expeditions 372 and 375 may also lie in the SSE rupture area and host a portion of the slip in these events. Hence, sampling and logging at this location provides insights into the composition, physical properties, and architecture of a shallow fault that may host slow slip. Expeditions 372 and 375 were designed to address three fundamental scientific objectives: Characterize the state and composition of the incoming plate and shallow fault near the trench, which comprise the protolith and initial conditions for fault zone rock at greater depth and which may itself host shallow slow slip; Characterize material properties, thermal regime, and stress conditions in the upper plate directly above the SSE source region; and Install observatories in the Pāpaku fault near the deformation front and in the upper plate above the SSE source to measure temporal variations in deformation, temperature, and fluid flow. The observatories will monitor volumetric strain (via pore pressure as a proxy) and the evolution of physical, hydrological, and chemical properties throughout the SSE cycle. Together, the coring, logging, and observatory data will test a suite of hypotheses about the fundamental mechanics and behavior of SSEs and their relationship to great earthquakes along the subduction interface

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

Preseismic fault creep and elastic wave amplitude precursors scale with lab earthquake magnitude for the continuum of tectonic failure modes

Tectonic faults fail in a continuum of modes from slow earthquakes to elastodynamic rupture. Precursory variations in elastic wavespeed and amplitude, interpreted as indicators of imminent failure, have been observed in limited natural settings and lab experiments where they are thought to arise from contact rejuvenation and microcracking within and around the fault zone. However, the physical mechanisms and connections to fault creep are poorly understood. Here we vary loading stiffness during frictional shear to generate a range of slip modes and measure fault zone properties using transmitted elastic waves. We find that elastic wave amplitudes show clear changes before fault failure. The temporal onset of amplitude reduction scales with lab earthquake magnitude and the magnitude of this reduction varies with fault slip. Our data provide clear evidence of precursors to lab earthquakes and suggest that continuous seismic monitoring could be useful for assessing fault state and seismic hazard potential

Competition between preslip and deviatoric stress modulates precursors for laboratory earthquakes

Variations in elastic wave velocity and amplitude prior to failure have been documented in laboratory experiments as well as in a limited number of crustal earthquakes. These variations have generally been attributed to fault zone healing, changes in crack density, or pore fluid effects modulated dilatation or fault slip. However, the relationships between amplitude and velocity variations during the seismic cycle, and the underlying mechanisms of precursors to failure remain poorly understood. Here, we perform frictional shear experiments and measure the evolution of elastic wave velocity and amplitude throughout the laboratory seismic cycle. We find that elastic amplitudes and velocities undergo clear preseismic variations prior to fault failure. While preseismic amplitude reduction occurs early in the interseismic period, wave speed reduces later, just prior to failure. We perform a complementary set of stress oscillation experiments to quantify the response of seismic amplitudes and velocities to variations in the stress tensor. Taken together, our results indicate that preseismic amplitude variations are primarily controlled by fault slip rate and acceleration. On the other hand, elastic velocity responds to a combination of fault preslip which reduces seismic wavespeed and increasing stress in the wallrock, which increases wavespeed. Our data show that precursory changes in seismic wave speed may be more common than previously thought because they are masked by changes in wallrock stress. These results underscore the importance of continuous and long-term time-lapse monitoring of crustal faults for seismic hazard assessment and potential precursors to failure