Competition between preslip and deviatoric stress modulates precursors for laboratory earthquakes

Abstract 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

Laboratory observations of permeability enhancement by fluid pressure oscillation of in situ fractured rock

We report on laboratory experiments designed to investigate the influence of pore pressure oscillations on the effective permeability of fractured rock. Berea sandstone samples were fractured in situ under triaxial stresses of tens of megapascals, and deionized water was forced through the incipient fracture under conditions of steady and oscillating pore pressure. We find that short-term pore pressure oscillations induce long-term transient increases in effective permeability of the fractured samples. The magnitude of the effective permeability enhancements scales with the amplitude of pore pressure oscillations, and changes persist well after the stress perturbation. The maximum value of effective permeability enhancement is 5 × 10^(−16) m^2 with a background permeability of 1 × 10^(−15) m^2; that is, the maximum enhanced permeability is 1.5 × 10^(−15) m^2. We evaluate poroelastic effects and show that hydraulic storage release does not explain our observations. Effective permeability recovery following dynamic oscillations occurs as the inverse square root of time. The recovery indicates that a reversible mechanism, such as clogging/unclogging of fractures, as opposed to an irreversible one, like microfracturing, is responsible for the transient effective permeability increase. Our work suggests the feasibility of dynamically controlling the effective permeability of fractured systems. The result has consequences for models of earthquake triggering and permeability enhancement in fault zones due to dynamic shaking from near and distant earthquakes

Acoustically-induced slip in sheared granular layers: application to dynamic earthquake triggering

A fundamental mystery in earthquake physics is ``how can an earthquake be triggered by distant seismic sources?'' Here, we use discrete element method simulations of a granular layer, during stick-slip, that is subject to transient vibrational excitation to gain further insight into the physics of dynamic earthquake triggering. Using Coulomb friction law for grains interaction, we observe delayed triggering of slip in the granular gouge. We find that at a critical vibrational amplitude (strain) there is an abrupt transition from negligible time-advanced slip (clock advance) to full clock advance, {\it i.e.}, transient vibration and triggered slip are simultaneous. The critical strain is order of $10^{-6}$, similar to observations in the laboratory and in Earth. The transition is related to frictional weakening of the granular layer due to a dramatic decrease in coordination number and the weakening of the contact force network. Associated with this frictional weakening is a pronounced decrease in the elastic modulus of the layer. The study has important implications for mechanisms of triggered earthquakes and induced seismic events and points out the underlying processes in response of the fault gouge to dynamic transient stresses

On the evolution of elastic properties during laboratory stick-slip experiments spanning the transition from slow slip to dynamic rupture

The physical mechanisms governing slow earthquakes remain unknown, as does the relationship between slow and regular earthquakes. To investigate the mechanism(s) of slow earthquakes and related quasi-dynamic modes of fault slip we performed laboratory experiments on simulated fault gouge in the double direct shear configuration. We reproduced the full spectrum of slip behavior, from slow to fast stick slip, by altering the elastic stiffness of the loading apparatus (k) to match the critical rheologic stiffness of fault gouge (kc). Our experiments show an evolution from stable sliding, when k>kc, to quasi-dynamic transients when k ~ kc, to dynamic instabilities when k<kc. To evaluate the microphysical processes of fault weakening we monitored variations of elastic properties. We find systematic changes in P wave velocity (Vp) for laboratory seismic cycles. During the coseismic stress drop, seismic velocity drops abruptly, consistent with observations on natural faults. In the preparatory phase preceding failure, we find that accelerated fault creep causes a Vp reduction for the complete spectrum of slip behaviors. Our results suggest that the mechanics of slow and fast ruptures share key features and that they can occur on same faults, depending on frictional properties. In agreement with seismic surveys on tectonic faults our data show that their state of stress can be monitored by Vp changes during the seismic cycle. The observed reduction in Vp during the earthquake preparatory phase suggests that if similar mechanisms are confirmed in nature high-resolution monitoring of fault zone properties may be a promising avenue for reliable detection of earthquake precursors

Slip-rate-dependent friction as a universal mechanism for slow slip events

A growing body of observations worldwide has documented fault slip transients that radiate little or no seismic energy. The mechanisms that govern these slow slip events (SSEs) and their wide range of depths, slip rates, durations, stress drops and recurrence intervals remain poorly known. Here we show that slow slip can be explained by a transition from rate-weakening frictional sliding at low slip rates towards rate-neutral or rate-strengthening behaviour at higher slip rates, as has been observed experimentally. We use numerical simulations to illustrate that this rate-dependent transition quantitatively explains the experimental data for natural fault rocks representative of materials in the source regions of SSEs. With a standard constant-parameter rate-and-state friction law, SSEs arise only near the threshold for slip instability. The inclusion of velocity-dependent friction parameters substantially broadens the range of conditions for slow slip occurrence, and produces a wide range of event characteristics, which include stress drop, duration and recurrence, as observed in nature. Upscaled numerical simulations that incorporate parameters consistent with laboratory measurements can reproduce geodetic observations of repeating SSEs on tectonic faults. We conclude that slip-rate-dependent friction explains the ubiquitous occurrence of SSEs in a broad range of geological environments

Foreshock properties illuminate nucleation processes of slow and fast laboratory earthquakes

Laboratory experiments demonstrate that prior to fast laboratory earthquakes the fault begins to unlock and creep, causing foreshocks to coalesce in both space and time. This demonstrates that the evolution of foreshocks is closely connected to the fault slip velocity.Understanding the connection between seismic activity and the earthquake nucleation process is a fundamental goal in earthquake seismology with important implications for earthquake early warning systems and forecasting. We use high-resolution acoustic emission (AE) waveform measurements from laboratory stick-slip experiments that span a spectrum of slow to fast slip rates to probe spatiotemporal properties of laboratory foreshocks and nucleation processes. We measure waveform similarity and pairwise differential travel-times (DTT) between AEs throughout the seismic cycle. AEs broadcasted prior to slow labquakes have small DTT and high waveform similarity relative to fast labquakes. We show that during slow stick-slip, the fault never fully locks, and waveform similarity and pairwise differential travel times do not evolve throughout the seismic cycle. In contrast, fast laboratory earthquakes are preceded by a rapid increase in waveform similarity late in the seismic cycle and a reduction in differential travel times, indicating that AEs begin to coalesce as the fault slip velocity increases leading up to failure. These observations point to key differences in the nucleation process of slow and fast labquakes and suggest that the spatiotemporal evolution of laboratory foreshocks is linked to fault slip velocity

Experimental Studies of Electrical Fields on a Breaking Rock Sample

When a rock sample is pressed by a force, the pressure on the crystal lattice generates an electrical field around the quart grains due to the piezoelectric effect. If a rock is saturated by conductive fluid, the relative motion between the pore fluid and the matrix solid generates an electromagnetic field due to seismoelectric conversion, and the permeating of fluid into new microcracks made by the pressure changes the fluid distribution and the natural potential level. In this paper, we measure the electrical fields on dry and water-saturated Westerly granite cylinder samples during their breaking. Experimental results show that there are two kinds of mechanisms that generate two kinds of electrical fields during rock breaking: (1) Pressure, or rock breaking, generates an electrical potential on the dry rock surface due to piezoelectric effect; and (2) the potential on a dry sample due to a piezoelectric effect is small, and its polarization depends on the characteristic and orientation of quartz grains around the measurement point. Experiments with water-saturated granite samples record two electrical fields: An electromagnetic wave due to seismoelectric conversion, and the dc or low-frequency electrical potential due to the piezoelectric effect, which is an important indicator of rock breaking