A supershear transition mechanism for cracks

How frictional slip evolves

Earthquake-like ruptures break the contacts that form the frictional interface separating contacting bodies and mediate the onset of frictional motion (stick-slip). The slip (motion) of the interface immediately resulting from the rupture that initiates each stick-slip event is generally much smaller than the total slip logged over the duration of the event. Slip after the onset of friction is generally attributed to continuous motion globally attributed to 'dynamic friction'. Here we show, by means of direct measurements of real contact area and slip at the frictional interface, that sequences of myriad hitherto invisible, secondary ruptures are triggered immediately in the wake of each initial rupture. Each secondary rupture generates incremental slip that, when not resolved, may appear as steady sliding of the interface. Each slip increment is linked, via fracture mechanics, to corresponding variations of contact area and local strain. Only by accounting for the contributions of these secondary ruptures can the accumulated interface slip be described. These results have important ramifications both to our fundamental understanding of frictional motion as well as to the essential role of aftershocks within natural faults in generating earthquake-mediated slip.

Transonic dislocation propagation in diamond

The motion of line defects (dislocations) has been studied for more than 60 years, but the maximum speed at which they can move is unresolved. Recent models and atomistic simulations predict the existence of a limiting velocity of dislocation motion between the transonic and subsonic ranges at which the self-energy of dislocation diverges, though they do not deny the possibility of the transonic dislocations. We used femtosecond x-ray radiography to track ultrafast dislocation motion in shock-compressed single-crystal diamond. By visualizing stacking faults extending faster than the slowest sound wave speed of diamond, we show the evidence of partial dislocations at their leading edge moving transonically. Understanding the upper limit of dislocation mobility in crystals is essential to accurately model, predict, and control the mechanical properties of materials under extreme conditions.

The break of earthquake asperities imaged by distributed acoustic sensing

Rupture imaging of megathrust earthquakes with global seismic arrays revealed frequency-dependent rupture signatures1-4, but the role of high-frequency radiators remains unclear3-5. Similar observations of the more abundant crustal earthquakes could provide critical constraints but are rare without ultradense local arrays6,7. Here we use distributed acoustic sensing technology8,9 to image the high-frequency earthquake rupture radiators. By converting a 100-kilometre dark-fibre cable into a 10,000-channel seismic array, we image four high-frequency subevents for the 2021 Antelope Valley, California, moment-magnitude 6.0 earthquake. After comparing our results with long-period moment-release10,11 and dynamic rupture simulations, we suggest that the imaged subevents are due to the breaking of fault asperities-stronger spots or pins on the fault-that substantially modulate the overall rupture behaviour. An otherwise fading rupture propagation could be promoted by the breaking of fault asperities in a cascading sequence. This study highlights how we can use the extensive pre-existing fibre networks12 as high-frequency seismic antennas to systematically investigate the rupture process of regional moderate-sized earthquakes. Coupled with dynamic rupture modelling, it could improve our understanding of earthquake rupture dynamics.

Laboratory earthquakes decipher control and stability of rupture speeds

Earthquakes are destructive natural hazards with damage capacity dictated by rupture speeds. Traditional dynamic rupture models predict that earthquake ruptures gradually accelerate to the Rayleigh wave speed with some of them further jumping to stable supershear speeds above the Eshelby speed (~[Formula: see text] times S wave speed). However, the 2018 M_w 7.5 Palu earthquake, among several others, significantly challenges such a viewpoint. Here we generate spontaneous shear ruptures on laboratory faults to confirm that ruptures can indeed attain steady subRayleigh or supershear propagation speeds immediately following nucleation. A self-similar analysis of dynamic rupture confirms our observation, leading to a simple model where the rupture speed is uniquely dependent on a driving load. Our results reproduce and explain a number of enigmatic field observations on earthquake speeds, including the existence of stable subEshelby supershear ruptures, early onset of supershear ruptures, and the correlation between the rupture speed and the driving load.

The role of pore fluids in supershear earthquake ruptures

The intensity and damage potential of earthquakes are linked to the speed at which rupture propagates along sliding crustal faults. Most earthquakes are sub-Rayleigh, with ruptures that are slower than the surface Rayleigh waves. In supershear earthquakes, ruptures are faster than the shear waves, leading to sharp pressure concentrations and larger intensities compared with the more common sub-Rayleigh ones. Despite significant theoretical and experimental advances over the past two decades, the geological and geomechanical controls on rupture speed transitions remain poorly understood. Here we propose that pore fluids play an important role in explaining earthquake rupture speed: the pore pressure may increase sharply at the compressional front during rupture propagation, promoting shear failure ahead of the rupture front and accelerating its propagation into the supershear range. We characterize the transition from sub-Rayleigh to supershear rupture in fluid-saturated rock, and show that the proposed poroelastic weakening mechanism may be a controlling factor for intersonic earthquake ruptures.

Signature of transition to supershear rupture speed in the coseismic off-fault damage zone

Most earthquake ruptures propagate at speeds below the shear wave velocity within the crust, but in some rare cases, ruptures reach supershear speeds. The physics underlying the transition of natural subshear earthquakes to supershear ones is currently not fully understood. Most observational studies of supershear earthquakes have focused on determining which fault segments sustain fully grown supershear ruptures. Experimentally cross-validated numerical models have identified some of the key ingredients required to trigger a transition to supershear speed. However, the conditions for such a transition in nature are still unclear, including the precise location of this transition. In this work, we provide theoretical and numerical insights to identify the precise location of such a transition in nature. We use fracture mechanics arguments with multiple numerical models to identify the signature of supershear transition in coseismic off-fault damage. We then cross-validate this signature with high-resolution observations of fault zone width and early aftershock distributions. We confirm that the location of the transition from subshear to supershear speed is characterized by a decrease in the width of the coseismic off-fault damage zone. We thus help refine the precise location of such a transition for natural supershear earthquakes.

Pressure shock fronts formed by ultra-fast shear cracks in viscoelastic materials

Spontaneously propagating cracks in solids emit both pressure and shear waves. When a shear crack propagates faster than the shear wave speed of the material, the coalescence of the shear wavelets emitted by the near-crack-tip region forms a shock front that significantly concentrates particle motion. Such a shock front should not be possible for pressure waves, because cracks should not be able to exceed the pressure wave speed in isotropic linear-elastic solids. In this study, we present full-field experimental measurements of dynamic shear cracks in viscoelastic polymers that result in the formation of a pressure shock front, in addition to the shear one. The apparent violation of classic theories is explained by the strain-rate-dependent material behavior of polymers, where the crack speed remains below the highest pressure wave speed prevailing locally around the crack tip. These findings have important implications for the physics and dynamics of shear cracks such as earthquakes.

Propagating shear cracks in solids emit both shear and pressure waves, but it is usually thought that only shear waves coalesce to form shock fronts when the crack exceeds the shear wave speed. Here, the authors show that local material stiffening can further increase rupture speed and produce pressure shock fronts that hint at supersonic propagation.

The equation of motion for supershear frictional rupture fronts

The rupture fronts that mediate the onset of frictional sliding may propagate at speeds below the Rayleigh wave speed or may surpass the shear wave speed and approach the longitudinal wave speed. While the conditions for the transition from sub-Rayleigh to supershear propagation have been studied extensively, little is known about what dictates supershear rupture speeds and how the interplay between the stresses that drive propagation and interface properties that resist motion affects them. By combining laboratory experiments and numerical simulations that reflect natural earthquakes, we find that supershear rupture propagation speeds can be predicted and described by a fracture mechanics-based equation of motion. This equation of motion quantitatively predicts rupture speeds, with the velocity selection dictated by the interface properties and stress. Our results reveal a critical rupture length, analogous to Griffith's length for sub-Rayleigh cracks, below which supershear propagation is impossible. Above this critical length, supershear ruptures can exist, once excited, even for extremely low preexisting stress levels. These results significantly improve our fundamental understanding of what governs the speed of supershear earthquakes, with direct and important implications for interpreting their unique supershear seismic radiation patterns.

Interplay between Process Zone and Material Heterogeneities for Dynamic Cracks

Using an elastodynamic boundary integral formulation coupled with a cohesive model, we study the problem of a dynamic rupture front propagating along an heterogeneous plane. We show that small-scale heterogeneities facilitate the supershear transition of a mode-II crack. The elastic pulses radiated during front accelerations explain how microscopic variations of fracture toughness change the macroscopic rupture dynamics. Perturbations of dynamic fronts are then systematically studied with different microstructures and loading conditions. The process zone size is the intrinsic length scale controlling heterogeneous dynamic rupture. The ratio of this length scale to asperity size is proposed as an indicator to transition from quasihomogeneous to heterogeneous fracture. Moreover, we discuss how the shortening of the process zone size with increasing crack speed brings the front to interact with smaller details of the microstructure. This study shines new light on recent experiments reporting perturbations of dynamic rupture fronts, which intensify with crack propagation speed.

Twin motion faster than the speed of sound

Twin growth is commonly thought to be bounded by the velocity of shear waves C(T) at which the information about this mechanical process travels in the material. Here, we report on experimental evidence of twin growth faster than the material's speed of sound. Driven by an electric field, needle twins in a ferroelectric crystal grew at intersonic speed, with an estimated average velocity close to square root(2) C(T). These results strengthen recent theoretical indications of intersonic dislocation motion, and contribute to the understanding of several twin motion-related processes.

Laboratory Earthquakes: The Sub-Rayleigh-to-Supershear Rupture Transition

We report on the experimental observation of spontaneously nucleated supershear rupture and on the visualization of sub-Rayleigh-to-supershear rupture transitions in frictionally held interfaces. The laboratory experiments mimic natural earthquakes. The results suggest that under certain conditions supershear rupture propagation can be facilitated during large earthquake events.

Observation of long supershear rupture during the magnitude 8.1 Kunlunshan earthquake

The 2001 Kunlunshan earthquake was an extraordinary event that produced a 400-km-long surface rupture. Regional broadband recordings of this event provide an opportunity to accurately observe the speed at which a fault ruptures during an earthquake, which has important implications for seismic risk and for understanding earthquake physics. We determined that rupture propagated on the 400-km-long fault at an average speed of 3.7 to 3.9 km/s, which exceeds the shear velocity of the brittle part of the crust. Rupture started at sub-Rayleigh wave velocity and became supershear, probably approaching 5 km/s, after about 100 km of propagation.