Fault lubrication during earthquakes

Seismic fault slip at depths simulated by high-velocity friction experiments under hydrothermal conditions

Seismic fault slip and rupture propagation often occur at crustal depths in the presence of hot and pressurized aqueous fluids (i.e., hydrothermal conditions). Previous experiments investigated fault frictional properties under hydrothermal conditions, but at imposed subseismic fault slip velocities (V ~μm/s). Here, using a rotary-shear apparatus equipped with a hydrothermal pressure vessel, we study friction at seismic slip velocities (V = 1.5 m/s) of gabbro- and marble-built faults under temperatures of 40 to 400 °C and pore water pressure of 30 MPa. We find that with increasing initial water temperature (Tamb), the dynamic friction during initial slip acceleration and subsequent high-velocity sliding decreases for both gabbro- and marble-built faults, while the slip-weakening distance decreases for gabbro but increases for marble. Then, during rapid deceleration at the end of sliding, frictional strength recovery decreases for gabbro with increasing Tamb and increases for marble independently of Tamb. As in previous experiments performed at room Tamb, the mechanical and microstructural data, plus numerical modeling, suggest that the seismic fault weakening mechanisms shift from flash heating to bulk melting for gabbro, and from flash heating to grain boundary sliding accommodated by diffusion creep for marble, with their activation processes depending on Tamb. Our results demonstrate the effects of ambient temperature on seismic fault friction, which contribute to changes in fault strength and dynamic weakening processes at crustal depths and should be considered in earthquake rupture modeling.

The 2023 Mw 7.8-7.7 Kahramanmaraş earthquakes were loosely slip-predictable

Understanding the behavior of large earthquakes over multiple seismic cycles is limited by short time spans of observations compared to recurrence intervals. Most of large instrumentally-recorded earthquakes have occurred on faults lacking well-documented histories of past events. The 2023 Mw 7.8-7.7 Kahramanmaraş earthquake doublet is exceptional as it ruptured multiple segments of the East Anatolian Fault (EAF) system, where historical records of devastating earthquakes span over two millennia. Here, we use historical earthquake records, measurements of interseismic deformation, and published slip models of the 2023 events to evaluate the recurrence patterns of large earthquakes. We compare slip deficit that accrued on each fault segment since the respective penultimate events to the average coseismic slip of the 2023 doublet. We find that the coseismic slip equaled to or exceeded the accumulated slip deficit, suggesting that the slip-predictable recurrence model applies as a lower bound on strain release during the Kahramanmaraş earthquakes.

India-Eurasia convergence speed-up by passive-margin sediment subduction

The fast increase of convergence rate between India and Eurasia around 65 million years ago (Ma)-from approximately 8 cm yr-1 to a peak rate of approximately 18 cm yr-1-remains a complex geological event to explain1-8, given the inherent uncertainty surrounding the tectonic history and the intricate interplay of forces influencing plate speed9-11. Here we use a combination of geochemical analysis and geodynamic modelling to propose that this rapid convergence can be explained by sediment subduction derived from the northern Indian passive margin. Through isotope and trace element analysis, we find an enhanced contribution of terrigenous sediment melt to the mantle source of the Gangdese magmatic rocks around 65 Ma, concurrent with the acceleration of India-Eurasia convergence. Numerical experiments suggest that subduction of sediments more than 1 km thick covering an approximately 1,000-km-wide ocean basin abutting the northern Indian passive margin starting from 65 Ma could have spurred the increased convergence rate and further led to significant crustal extension, consistent with empirical observations. Our study implies that the acceleration of India-Eurasia convergence marks the arrival of passive-margin-derived sediments, constraining the initial India-Eurasia collision to be around 60 Ma. It further suggests that temporary accelerations in subduction rates might be a common feature at the final stage of continental assembly.

Physics-based assessment of earthquake potential on the Anninghe-Zemuhe fault system in southwestern China

The seismic hazard of a fault system is controlled by the maximum possible earthquake magnitude it can host. However, existing methods to estimate maximum magnitudes can result in large uncertainties or ignore their temporal evolution. Here, we show how the maximum possible earthquake magnitude of a fault system can be assessed by combining high-resolution fault coupling maps with a physics-based model from three-dimensional dynamic fracture mechanics confirmed by dynamic rupture simulations. We demonstrate the method on the Anninghe-Zemuhe fault system in southwestern China, where dense near-fault geodetic data has been acquired. Our results show that this fault system currently has the potential to generate Mw7.0 earthquakes with maximum magnitudes increasing to Mw7.3 by 2200. These results are supported by the observed rupture extents and recurrence times of historical earthquakes and the b values of current seismicity. Our work provides a practical way to assess the earthquake potential of natural faults.

Fault size-dependent fracture energy explains multiscale seismicity and cascading earthquakes

Earthquakes vary in size over many orders of magnitude, often rupturing in complex multifault and multievent sequences. Despite the large number of observed earthquakes, the scaling of the earthquake energy budget remains enigmatic. We propose that fundamentally different fracture processes govern small and large earthquakes. We combined seismological observations with physics-based earthquake models, finding that both dynamic weakening and restrengthening effects are non-negligible in the energy budget of small earthquakes. We established a linear scaling relationship between fracture energy and fault size and a break in scaling with slip. We applied this scaling using supercomputing and unveiled large dynamic rupture earthquake cascades involving >700 multiscale fractures within a fault damage zone. We provide a simple explanation for seismicity across all scales with implications for comprehending earthquake genesis and multifault rupture cascades.

Earthquake energy dissipation in a fracture mechanics framework

Earthquakes are rupture-like processes that propagate along tectonic faults and cause seismic waves. The propagation speed and final area of the rupture, which determine an earthquake's potential impact, are directly related to the nature and quantity of the energy dissipation involved in the rupture process. Here, we present the challenges associated with defining and measuring the energy dissipation in laboratory and natural earthquakes across many scales. We discuss the importance and implications of distinguishing between energy dissipation that occurs close to and far behind the rupture tip, and we identify open scientific questions related to a consistent modeling framework for earthquake physics that extends beyond classical Linear Elastic Fracture Mechanics.

Physical state of water controls friction of gabbro-built faults

Earthquakes often occur along faults in the presence of hot, pressurized water. Here we exploit a new experimental device to study friction in gabbro faults with water in vapor, liquid and supercritical states (water temperature and pressure up to 400 °C and 30 MPa, respectively). The experimental faults are sheared over slip velocities from 1 μm/s to 100 mm/s and slip distances up to 3 m (seismic deformation conditions). Here, we show with water in the vapor state, fault friction decreases with increasing slip distance and velocity. However, when water is in the liquid or supercritical state, friction decreases with slip distance, regardless of slip velocity. We propose that the formation of weak minerals, the chemical bonding properties of water and (elasto)hydrodynamic lubrication may explain the weakening behavior of the experimental faults. In nature, the transition of water from liquid or supercritical to vapor state can cause an abrupt increase in fault friction that can stop or delay the nucleation phase of an earthquake.

Frictional melting mechanisms of rocks during earthquake fault slip

Rapid slip, at rates in the order of 1 m/s or more, may induce frictional melting in rocks during earthquakes. The short-lived melting has been thought to be a disequilibrium process, for decades. We conducted frictional melting experiments on acidic, basic, and ultrabasic silicate rocks at a slip rate of 1.3 m/s. The experiments and microstructural observations reveal that all minerals in the rocks are melted at temperatures below their known melting temperatures (Tm); e.g., quartz is melted at ~ 1000-1200 °C, not ~ 1720 °C, while olivine at ~ 1300 °C, rather than ~ 1700 °C. The low-temperature melting is incompatible with the conventional disequilibrium melting, and may be caused predominantly by grain size reduction and phase boundary reactions during the early and later stages of slip, respectively. The newly estimated Tm and the melting mechanisms should be considered for understanding the mechanics of earthquakes, landslides, and caldera collapses.

Dynamics, interactions and delays of the 2019 Ridgecrest rupture sequence

The observational difficulties and the complexity of earthquake physics have rendered seismic hazard assessment largely empirical. Despite increasingly high-quality geodetic, seismic and field observations, data-driven earthquake imaging yields stark differences and physics-based models explaining all observed dynamic complexities are elusive. Here we present data-assimilated three-dimensional dynamic rupture models of California's biggest earthquakes in more than 20 years: the moment magnitude (M_w) 6.4 Searles Valley and M_w 7.1 Ridgecrest sequence, which ruptured multiple segments of a non-vertical quasi-orthogonal conjugate fault system¹. Our models use supercomputing to find the link between the two earthquakes. We explain strong-motion, teleseismic, field mapping, high-rate global positioning system and space geodetic datasets with earthquake physics. We find that regional structure, ambient long- and short-term stress, and dynamic and static fault system interactions driven by overpressurized fluids and low dynamic friction are conjointly crucial to understand the dynamics and delays of the sequence. We demonstrate that a joint physics-based and data-driven approach can be used to determine the mechanics of complex fault systems and earthquake sequences when reconciling dense earthquake recordings, three-dimensional regional structure and stress models. We foresee that physics-based interpretation of big observational datasets will have a transformative impact on future geohazard mitigation.

Dueling dynamics of low-angle normal fault rupture with splay faulting and off-fault damage

Despite a lack of modern large earthquakes on shallowly dipping normal faults, Holocene M_w> 7 low-angle normal fault (LANF; dip<30°) ruptures are preserved paleoseismically and inferred from historical earthquake and tsunami accounts. Even in well-recorded megathrust earthquakes, the effects of non-linear off-fault plasticity and dynamically reactivated splay faults on shallow deformation and surface displacements, and thus hazard, remain elusive. We develop data-constrained 3D dynamic rupture models of the active Mai'iu LANF that highlight how multiple dynamic shallow deformation mechanisms compete during large LANF earthquakes. We show that shallowly-dipping synthetic splays host more coseismic slip and limit shallow LANF rupture more than steeper antithetic splays. Inelastic hanging-wall yielding localizes into subplanar shear bands indicative of newly initiated splay faults, most prominently above LANFs with thick sedimentary basins. Dynamic splay faulting and sediment failure limit shallow LANF rupture, modulating coseismic subsidence patterns, near-shore slip velocities, and the seismic and tsunami hazards posed by LANF earthquakes.

Dynamics between earthquakes, volcanic eruptions, and geothermal energy exploitation in Japan

Intruding magma brings high temperatures close to the surface, thus offering possibilities for harnessing large amounts of heat for geothermal exploitation. Mount Aso in southern Japan showed frequent volcanic activity during 2016, accompanied by significant earthquake activities with tens of thousands of aftershocks (Kumamoto sequence). Here we investigate the influence of earthquake/volcanic activity on the future productivity of nearby geothermal power plants to determine whether the activity is detrimental or beneficial to energy exploitation. Model results show an increase in [Formula: see text] pressure and temperature with a spatio-temporal correlation between modeled earthquake locations and aftershock decay rates along the entire sequence, showing that seismic activity opened pre-existing vertical cracks providing pathways for the ascending magma. Interestingly, the minor but still significant eruption of Mount Aso in October 2021 may have enhanced future geothermal power generation, indicating a vigorous and active system, possibly increasing the future geothermal power production.

Coseismic fault sealing and fluid pressurization during earthquakes

Earthquakes occur because faults weaken with increasing slip and slip rate. Thermal pressurization (TP) of trapped pore fluids is deemed to be a widespread coseismic fault weakening mechanism. Yet, due to technical challenges, experimental evidence of TP is limited. Here, by exploiting a novel experimental configuration, we simulate seismic slip pulses (slip rate 2.0 m/s) on dolerite-built faults under pore fluid pressures up to 25 MPa. We measure transient sharp weakening, down to almost zero friction and concurrent with a spike in pore fluid pressure, which interrupts the exponential-decay slip weakening. The interpretation of mechanical and microstructural data plus numerical modeling suggests that wear and local melting processes in experimental faults generate ultra-fine materials to seal the pressurized pore water, causing transient TP spikes. Our work suggests that, with wear-induced sealing, TP may also occur in relatively permeable faults and could be quite common in nature.

Integrated rupture mechanics for slow slip events and earthquakes

Slow slip events occur worldwide and could trigger devastating earthquakes, yet it is still debated whether their moment-duration scaling is linear or cubic and a fundamental model unifying slow and fast earthquakes is still lacking. Here, we show that the rupture propagation of simulated slow and fast earthquakes can be predicted by a newly-developed three-dimensional theory of dynamic fracture mechanics accounting for finite rupture width, an essential ingredient missing in previous theories. The complete spectrum of rupture speeds is controlled by the ratio of fracture energy to energy release rate. Shear stress heterogeneity can produce a cubic scaling on a single fault while effective normal stress variability produces a linear scaling on a population of faults, which reconciles the debated scaling relations. This model provides a new framework to explain how slow slip might lead to earthquakes and opens new avenues for seismic hazard assessment integrating seismological, laboratory and theoretical developments.

Heterogeneous Subgreenschist Deformation in an Exhumed Sediment-Poor Mélange

Many described subduction complexes (or mélanges) exhumed from seismogenic depths comprise thick, turbidite-dominated sequences with deformed zones containing clasts or boudins of more competent sandstone and/or basalt. In contrast, many active subduction zones have a relatively small thickness of sedimentary inputs (<2 km), turbidite sequences are commonly accreted rather than subducted, and the role of pelagic sediments and basalt (lavas and hyaloclastites) in the deforming zone near the plate interface at <20 km depth is poorly understood. Field investigation of Neoproterozoic oceanic sequences accreted in the Gwna Complex, Anglesey, UK, reveals repeated lenticular slices of variably sampled ocean plate stratigraphy (OPS) bounded by thin mélange-bearing shear zones. Mélange matrix material is derived from adjacent OPS lithologies and is either dominantly illitic, likely derived from altered siliciclastic sediment, or chloritic, likely derived from altered volcanics. In the illitic mélange, mutually cross-cutting phyllosilicate foliation and variably deformed chlorite-quartz-calcite veins suggest ductile creep was cyclically punctuated by transient, localized fluid pulses. Chlorite thermometry indicates the veins formed at 260 ± 10°C. In the chloritic mélange, recrystallized through-going calcite veins are deformed to shear strains of 4-5 within a foliated chlorite matrix, suggesting calcite veins in subducting volcanics may localize deformation in the seismogenic zone. Shear stress-strain rate curves constructed using existing empirical relationships in a simplified shear zone geometry predict that slip velocities varied depending on pore fluid pressure; models predict slow slip velocities preferentially by frictional sliding in chlorite, at pore fluid pressures greater than hydrostatic but less than lithostatic.

Intermittent lab earthquakes in dynamically weakening fault gouge

Large and destructive earthquakes on mature faults in Earth's crust occur as slip in a layer of a fine granular material-fault gouge-produced by comminution during sliding^1,2. A range of insights into the frictional resistance of faults-one of the main factors controlling earthquake nucleation, dynamic propagation and arrest, and hence the destructive ground shaking of earthquakes^2,3-has been obtained in experiments with spatially uniform slip imposed in small samples^2,4-21. However, how various features of gouge friction combine to determine spontaneous progression of earthquakes is difficult to study in the lab owing to substantial challenges with sample sizes and adequate imaging²². Here, using lab experiments, we show that spontaneously propagating dynamic ruptures navigate a fault region with fine rock gouge through complex, intermittent slip processes with dramatic friction evolution. These include repeated arrest of rupture propagation caused by friction strengthening at lower slip rates and dynamic earthquake re-nucleation enabled by pronounced rapid friction weakening at higher slip rates consistent with flash heating^8,12,23. The spontaneous repeated weakening and strengthening of friction in fine rock gouge highlights the fundamental dependence of friction on slip rate and associated processes, such as shear heating, localization and delocalization of shear, and dilation and compaction of the shear layer^6,7,9-21. Our findings expand experimental support^9,11 of the concept that co-seismic weakening may enable earthquake rupture to break through stable fault regions^24,25, with substantial implications for seismic hazard.

The impact of melt versus mechanical wear on the formation of pseudotachylyte veins in accretionary complexes

Whether seismic rupture propagates over large distances to generate mega-earthquakes or is rapidly aborted mainly depends on the slip processes within the fault core, including particularly frictional melting or intense grain-size reduction and amorphization. The record of seismic slip in exhumed fault zones consists in many instances in Black Faults Rocks, dark and glass-like-filled aphanitic veins that have been interpreted as resulting from the quenching of frictional melts, i.e. pseudotachylytes. Such interpretation has nevertheless been questioned as similar macro to nano-microstructures have been observed either on intensely comminuted natural fault rocks or on slow creep experiments conducted on crustal rocks, where melting is absent. Here, we report a new dataset of Raman Spectroscopy of Carbonaceous Material analyses, aimed at discriminating the slip weakening processes operating in the fault core during slip. Using high spatial resolution profiles on natural Black Fault Rocks from exhumed accretionary complexes and an experimentally calibrated modelling of Raman intensity ratio evolution with temperature, we assessed different scenarios of temperature evolution during fault slip. None of them is able to account for the distribution of Raman signal, so that in the three studied Black Fault Rocks interpreted so far as natural pseudotachylytes, Raman Spectroscopy of Carbonaceous Material rather reflects the effect of intense and localized strain during fault slip. Furthermore, the absence of thermal imprint on Raman signal puts upper bounds on the temperature reached within the fault zone. If one cannot rule out the occurrence of high and short-lived temperature increase due to friction, the latter was not high enough as to melt the large quartz fraction of the fault zone rocks.

Fault rock heterogeneity can produce fault weakness and reduce fault stability

Geological heterogeneity is abundant in crustal fault zones; however, its role in controlling the mechanical behaviour of faults is poorly constrained. Here, we present laboratory friction experiments on laterally heterogeneous faults, with patches of strong, rate-weakening quartz gouge and weak, rate-strengthening clay gouge. The experiments show that the heterogeneity leads to a significant reduction in strength and frictional stability in comparison to compositionally identical faults with homogeneously mixed gouges. We identify a combination of weakening effects, including smearing of the weak clay; differential compaction of the two gouges redistributing normal stress; and shear localization producing stress concentrations in the strong quartz patches. The results demonstrate that geological heterogeneity and its evolution can have pronounced effects on fault strength and stability and, by extension, on the occurrence of slow-slip transients versus earthquake ruptures and the characteristics of the resulting events, and should be further studied in lab experiments and earthquake source modelling.

Nanometric flow and earthquake instability

Fault zones accommodate relative motion between tectonic blocks and control earthquake nucleation. Nanocrystalline fault rocks are ubiquitous in "principal slip zones" indicating that these materials are determining fault stability. However, the rheology of nanocrystalline fault rocks remains poorly constrained. Here, we show that such fault rocks are an order of magnitude weaker than their microcrystalline counterparts when deformed at identical experimental conditions. Weakening of the fault rocks is hence intrinsic, it occurs once nanocrystalline layers form. However, it is difficult to produce "rate weakening" behavior due to the low measured stress exponent, n, of 1.3 ± 0.4 and the low activation energy, Q, of 16,000 ± 14,000 J/mol implying that the material will be strongly "rate strengthening" with a weak temperature sensitivity. Failure of the fault zone nevertheless occurs once these weak layers coalesce in a kinematically favored network. This type of instability is distinct from the frictional instability used to describe crustal earthquakes.

Repeating caldera collapse events constrain fault friction at the kilometer scale

Fault friction is central to understanding earthquakes, yet laboratory rock mechanics experiments are restricted to, at most, meter scale. Questions thus remain as to the applicability of measured frictional properties to faulting in situ. In particular, the slip-weakening distance [Formula: see text] strongly influences precursory slip during earthquake nucleation, but scales with fault roughness and is challenging to extrapolate to nature. The 2018 eruption of K̄ılauea volcano, Hawaii, caused 62 repeatable collapse events in which the summit caldera dropped several meters, accompanied by [Formula: see text] 4.7 to 5.4 very long period (VLP) earthquakes. Collapses were exceptionally well recorded by global positioning system (GPS) and tilt instruments and represent unique natural kilometer-scale friction experiments. We model a piston collapsing into a magma reservoir. Pressure at the piston base and shear stress on its margin, governed by rate and state friction, balance its weight. Downward motion of the piston compresses the underlying magma, driving flow to the eruption. Monte Carlo estimation of unknowns validates laboratory friction parameters at the kilometer scale, including the magnitude of steady-state velocity weakening. The absence of accelerating precollapse deformation constrains [Formula: see text] to be [Formula: see text] mm, potentially much less. These results support the use of laboratory friction laws and parameters for modeling earthquakes. We identify initial conditions and material and magma-system parameters that lead to episodic caldera collapse, revealing that small differences in eruptive vent elevation can lead to major differences in eruption volume and duration. Most historical basaltic caldera collapses were, at least partly, episodic, implying that the conditions for stick-slip derived here are commonly met in nature.

Seismic Slip-Pulse Experiments Simulate Induced Earthquake Rupture in the Groningen Gas Field

Rock materials show dramatic dynamic weakening in large-displacement (m), high-velocity (∼1 m/s) friction experiments, providing a mechanism for the generation of large, natural earthquakes. However, whether such weakening occurs during induced M3-4 earthquakes (dm displacements) is unknown. We performed rotary-shear experiments on simulated fault gouges prepared from the source-, reservoir- and caprock formations present in the seismogenic Groningen gas field (Netherlands). Water-saturated gouges were subjected to a slip pulse reaching a peak circumferential velocity of 1.2-1.7 m/s and total displacements of 13-20 cm, at 2.5-20 MPa normal stress. The results show 22%-81% dynamic weakening within 5-12 cm of slip, depending on normal stress and gouge composition. At 20 MPa normal stress, dynamic weakening from peak friction coefficients of 0.4-0.9 to 0.19-0.27 was observed, probably through thermal pressurization. We infer that similar effects play a key role during induced seismic slip on faults in the Groningen and other reservoir systems.

Fast and Localized Temperature Measurements During Simulated Earthquakes in Carbonate Rocks

The understanding of earthquake physics is hindered by the poor knowledge of fault strength and temperature evolution during seismic slip. Experiments reproducing seismic velocity (∼1 m/s) allow us to measure both the evolution of fault strength and the associated temperature increase due to frictional heating. However, temperature measurements were performed with techniques having insufficient spatial and temporal resolution. Here we conduct high velocity friction experiments on Carrara marble rock samples sheared at 20 MPa normal stress, velocity of 0.3 and 6 m/s, and 20 m of total displacement. We measured the temperature evolution of the fault surface at the acquisition rate of 1 kHz and over a spatial resolution of ∼40 µm with an optical fiber conveying the infrared radiation to a two-color pyrometer. Temperatures up to 1,250°C and low coseismic fault shear strength are compatible with the activation of grain size dependent viscous creep.

Resolution and uncertainties in estimates of earthquake stress drop and energy release

Our models and understanding of the dynamics of earthquake rupture are based largely on estimates of earthquake source parameters, such as stress drop and radiated seismic energy. Unfortunately, the measurements, especially those of small and moderate-sized earthquakes (magnitude less than about 5 or 6), are not well resolved, containing significant random and potentially systematic uncertainties. The aim of this review is to provide a context in which to understand the challenges involved in estimating these measurements, and to assess the quality and reliability of reported measurements of earthquake source parameters. I also discuss some of the ways progress is being made towards more reliable parameter measurements. At present, whether the earthquake source is entirely self-similar, or not, and which factors and processes control the physics of the rupture remains, at least in the author's opinion, largely unconstrained. Detailed analysis of the best recorded earthquakes, using the increasing quantity and quality of data available, and methods less dependent on simplistic source models is one approach that may help provide better constraints. This article is part of the theme issue 'Fracture dynamics of solid materials: from particles to the globe'.

Fracture mechanics of rate-and-state faults and fluid injection induced slip

Propagation of a slip transient on a fault with rate- and state-dependent friction resembles a fracture whose near tip region is characterized by large departure of the slip velocity and fault strength from the steady-state sliding. We develop a near tip solution to describe this unsteady dynamics, and obtain the fracture energy G_c, dissipated in overcoming strength-excursion away from steady state, as a function of the rupture velocity v_r. This opens a possibility to model slip transients on rate-and-state faults as singular cracks characterized by approximately steady-state frictional resistance in the fracture bulk, and by a stress singularity with the intensity defined in terms of G_c(v_r) at the crack tip. In pursuing this route, we develop and use an analytical equation of motion to study 1-D slip driven by a combination of uniform background stress and a localized perturbation of the fault strength with the net Coulomb force ΔT. In the context of fluid injection, ΔT is a proxy for the injection volume V_inj. We then show that, for ongoing fluid injection, the propagation speed of a transient induced on a frictionally stable fault is bounded by a large-time limiting value proportional to the injection rate dV_inj/dt, while, for stopped injection, the maximum slip run-out distance is proportional to [Formula: see text]. This article is part of the theme issue 'Fracture dynamics of solid materials: from particles to the globe'.

Fluid pressurisation and earthquake propagation in the Hikurangi subduction zone

In subduction zones, seismic slip at shallow crustal depths can lead to the generation of tsunamis. Large slip displacements during tsunamogenic earthquakes are attributed to the low coseismic shear strength of the fluid-saturated and non-lithified clay-rich fault rocks. However, because of experimental challenges in confining these materials, the physical processes responsible for the coseismic reduction in fault shear strength are poorly understood. Using a novel experimental setup, we measured pore fluid pressure during simulated seismic slip in clay-rich materials sampled from the deep oceanic drilling of the Pāpaku thrust (Hikurangi subduction zone, New Zealand). Here, we show that at seismic velocity, shear-induced dilatancy is followed by pressurisation of fluids. The thermal and mechanical pressurisation of fluids, enhanced by the low permeability of the fault, reduces the energy required to propagate earthquake rupture. We suggest that fluid-saturated clay-rich sediments, occurring at shallow depth in subduction zones, can promote earthquake rupture propagation and slip because of their low permeability and tendency to pressurise when sheared at seismic slip velocities.

The earthquake cycle in the dry lower continental crust: insights from two deeply exhumed terranes (Musgrave Ranges, Australia and Lofoten, Norway)

This paper discusses the results of field-based geological investigations of exhumed rocks exposed in the Musgrave Ranges (Central Australia) and in Nusfjord (Lofoten, Norway) that preserve evidence for lower continental crustal earthquakes with focal depths of approximately 25-40 km. These studies have established that deformation of the dry lower continental crust is characterized by a cyclic interplay between viscous creep (mylonitization) and brittle, seismic slip associated with the formation of pseudotachylytes (a solidified melt produced during seismic slip along a fault in silicate rocks). Seismic slip triggers rheological weakening and a transition to viscous creep, which may be already active during the immediate post-seismic deformation along faults initially characterized by frictional melting and wall-rock damage. The cyclical interplay between seismic slip and viscous creep implies transient oscillations in stress and strain rate, which are preserved in the shear zone microstructure. In both localities, the spatial distribution of pseudotachylytes is consistent with a local (deep) source for the transient high stresses required to generate earthquakes in the lower crust. This deep source is the result of localized stress amplification in dry and strong materials generated at the contacts with ductile shear zones, producing multiple generations of pseudotachylyte over geological time. This implies that both the short- and the long-term rheological evolution of the dry lower crust typical of continental interiors is controlled by earthquake cycle deformation. This article is part of a discussion meeting issue 'Understanding earthquakes using the geological record'.

Microphysical Modeling of Carbonate Fault Friction at Slip Rates Spanning the Full Seismic Cycle

Laboratory studies suggest that seismogenic rupture on faults in carbonate terrains can be explained by a transition from high friction, at low sliding velocities (V), to low friction due to rapid dynamic weakening as seismic slip velocities are approached. However, consensus on the controlling physical processes is lacking. We previously proposed a microphysically based model (the "Chen-Niemeijer-Spiers" [CNS] model) that accounts for the (rate-and-state) frictional behavior of carbonate fault gouges seen at low velocities characteristic of rupture nucleation. In the present study, we extend the CNS model to high velocities (1 mm/s ≤ V ≤ 10 m/s) by introducing multiple grain-scale deformation mechanisms activated by frictional heating. As velocity and hence temperature increase, the model predicts a continuous transition in dominant deformation mechanisms, from frictional granular flow with partial accommodation by plasticity at low velocities and temperatures, to grain boundary sliding with increasing accommodation by solid-state diffusion at high velocities and temperatures. Assuming that slip occurs in a localized shear band, within which grain size decreases with increasing velocity, the model results capture the main mechanical trends seen in high-velocity friction experiments on room-dry calcite-rich rocks, including steady-state and transient aspects, with reasonable quantitative agreement and without the need to invoke thermal decomposition or fluid pressurization effects. The extended CNS model covers the full spectrum of slip velocities from earthquake nucleation to seismic slip rates. Since it is based on realistic fault structure, measurable microstructural state variables, and established deformation mechanisms, it may offer an improved basis for extrapolating lab-derived friction data to natural fault conditions.

Propagation of large earthquakes as self-healing pulses or mild cracks

Observations suggest that mature faults host large earthquakes at much lower levels of stress than their expected static strength^1-11. Potential explanations are that the faults are quasi-statically strong but experience considerable weakening during earthquakes, or that the faults are persistently weak, for example, because of fluid overpressure. Here we use numerical modelling to examine these competing theories for simulated earthquake ruptures that satisfy the well known observations of 1-10 megapascal stress drops and limited heat production. In that regime, quasi-statically strong but dynamically weak faults mainly host relatively sharp, self-healing pulse-like ruptures, with only a small portion of the fault slipping at a given time, whereas persistently weak faults host milder ruptures with more spread-out slip, which are called crack-like ruptures. We find that the sharper self-healing pulses, which exhibit larger dynamic stress changes compared to their static stress changes, result in much larger radiated energy than that inferred teleseismically for megathrust events¹². By contrast, milder crack-like ruptures on persistently weak faults, which produce comparable static and dynamic stress changes, are consistent with the seismological observations. The larger radiated energy of self-healing pulses is similar to the limited regional inferences available for crustal strike-slip faults. Our findings suggest that either large earthquakes rarely propagate as self-healing pulses, with potential differences between tectonic settings, or their radiated energy is substantially underestimated, raising questions about earthquake physics and the expected shaking from large earthquakes.

Molecular origin of sliding friction and flash heating in rock and heterogeneous materials

It is generally believed that earthquakes occur when faults weaken with increasing slip rates. An important factor contributing to this phenomenon is the faults’ dynamic friction, which may be reduced during earthquakes with high slip rates, a process known as slip-rate weakening. It has been hypothesized that the weakening phenomenon during fault slip may be activated by thermal pressurization of pores’ fluid and flash heating, a microscopic phenomenon in which heat is generated at asperity contacts due to high shear slip rates. Due to low thermal conductivity of rock, the heat generated at the contact points or surfaces cannot diffuse fast enough, thus concentrating at the contacts, increasing the local contact temperature, and reducing its frictional shear strength. We report the results of what we believe to be the first molecular scale study of the decay of the interfacial friction force in rock, observed in experiemntal studies and attributed to flash heating. The magnitude of the reduction in the shear stress and the local friction coefficients have been computed over a wide range of shear velocities V. The molecular simulations indicate that as the interfacial temperature increases, bonds between the atoms begin to break, giving rise to molecular-scale fracture that eventually produces the flash heating effect. The frequency of flash heating events increases with increasing sliding velocity, leaving increasingly shorter times for the material to relax, hence contributing to the increased interfacial temperature. If the material is thin, the heat quickly diffuses away from the interface, resulting in sharp decrease in the temperature immediately after flash heating. The rate of heat transfer is reduced significantly with increasing thickness, keeping most of the heat close to the interface and producing weakened material. The weakening behavior is demonstrated by computing the stress–strain diagram. For small strain rates there the frictional stress is essentially independent of the materials’ thickness. As the strain rate increases, however, the dependence becomes stronger. Specifically, the stress–strain diagrams at lower velocities V manifest a pronounced strength decrease over small distances, whereas they exhibit progressive increase in the shear stress at higher V, which is reminiscent of a transition from ductile behavior at high velocities to brittle response at low velocities.

Flow-to-Friction Transition in Simulated Calcite Gouge: Experiments and Microphysical Modeling

A (micro)physical understanding of the transition from frictional sliding to plastic or viscous flow has long been a challenge for earthquake cycle modeling. We have conducted ring-shear deformation experiments on layers of simulated calcite fault gouge under conditions close to the frictional-to-viscous transition previously established in this material. Constant velocity (v) and v-stepping tests were performed, at 550°C, employing slip rates covering almost 6 orders of magnitude (0.001-300 μm/s). Steady-state sliding transitioned from (strong) v-strengthening, flow-like behavior to v-weakening, frictional behavior, at an apparent "critical" velocity (v cr ) of ~0.1 μm/s. Velocity-stepping tests using v < v cr showed "semi-brittle" flow behavior, characterized by high stress sensitivity ("n-value") and a transient response resembling classical frictional deformation. For v ≥ v cr , gouge deformation is localized in a boundary shear band, while for v < v cr , the gouge is well-compacted, displaying a progressively homogeneous structure as the slip rate decreases. Using mechanical data and post-mortem microstructural observations as a basis, we deduced the controlling shear deformation mechanisms and quantitatively reproduced the steady-state shear strength-velocity profile using an existing micromechanical model. The same model also reproduces the observed transient responses to v-steps within both the flow-like and frictional deformation regimes. We suggest that the flow-to-friction transition strongly relies on fault (micro)structure and constitutes a net opening of transient microporosity with increasing shear strain rate at v < v cr , under normal stress-dependent or "semi-brittle" flow conditions. Our findings shed new insights into the microphysics of earthquake rupture nucleation and dynamic propagation in the brittle-to-ductile transition zone.

Earthquake slip surfaces identified by biomarker thermal maturity within the 2011 Tohoku-Oki earthquake fault zone

Extreme slip at shallow depths on subduction zone faults is a primary contributor to tsunami generation by earthquakes. Improving earthquake and tsunami risk assessment requires understanding the material and structural conditions that favor earthquake propagation to the trench. We use new biomarker thermal maturity indicators to identify seismic faults in drill core recovered from the Japan Trench subduction zone, which hosted 50 m of shallow slip during the M_w9.1 2011 Tohoku-Oki earthquake. Our results show that multiple faults have hosted earthquakes with displacement ≥ 10 m, and each could have hosted many great earthquakes, illustrating an extensive history of great earthquake seismicity that caused large shallow slip. We find that lithologic contrasts in frictional properties do not necessarily determine the likelihood of large shallow slip or seismic hazard.

Weak phases production and heat generation control fault friction during seismic slip

The triggering and magnitude of earthquakes is determined by the friction evolution along faults. Experimental results have revealed a drastic decrease of the friction coefficient for velocities close to the maximum seismic one, independently of the material studied. Due to the extreme loading conditions during seismic slip, many competing physical phenomena occur (like mineral decomposition, nanoparticle lubrication, melting among others) that are typically thermal in origin and are changing the nature of the material. Here we show that a large set of experimental data for different rocks can be described by such thermally-activated mechanisms, combined with the production of weak phases. By taking into account the energy balance of all processes during fault movement, we present a framework that reconciles the data, and is capable of explaining the frictional behavior of faults, across the full range of slip velocities (10^-9 to 10 m/s).

Across-Fault Velocity Gradients and Slip Behavior of the San Andreas Fault Near Parkfield

A long-lasting question in earthquake physics is why slip on faults occurs as creep or dynamic rupture. We compute passive measurements of the seismic P wave velocity gradient across the San Andreas Fault near Parkfield, where this transition of slip mode occurs at a scale of a few kilometers. Unbiased measurements are obtained through the application of a new Bayesian local earthquake tomographic code that avoids the imposition of any user-defined, initial velocity-contrast across the fault, or any damping scheme that may cause biased amplitude in retrieved seismic velocities. We observe that across-fault velocity gradients correlate with the slip behavior of the fault. The P wave velocity contrast decays from 20% in the fault section that experience dynamic rupture to 4% in the creeping section, suggesting that rapid change of material properties and attitude to sustain supra-hydrostatic fluid pressure are conditions for development of dynamic rupture. Low Vp and high Vp/Vs suggest that fault rheology at shallow depth is conversely controlled by low frictional strength material.

Breakdown of Scaling and Friction Weakening in Intermittent Granular Flow

Many materials are produced, processed and stored as grains, while granularity of matter can be crucial in triggering potentially catastrophic geological events like landslides, avalanches and earthquakes. The response of grain assemblies to shear stress is therefore of utmost relevance to both human and natural environment. At low shear rate a granular system flows intermittently by distinct avalanches. In such state the avalanche velocity in time is expected to follow a symmetrical and universal average behavior, whose dependence on the slip size reduces to a scale factor. Analyzing data from long lasting experiments, we observe a breakdown of this scaling: While in short slips velocity shows indeed a self-similar and symmetric profile, it does not in long slips. The investigation of frictional response in these different regimes evidences that this breakdown can be traced back to the onset of a friction weakening, which is of dynamical origin and can amplify instabilities exactly in this critical state, the most frequent state for natural hazards.

Energy-flux control of the steady-state, creep, and dynamic slip modes of faults

Faults exhibit a gamut of slip styles from stable sliding and creep events to earthquakes. These slip styles are affected by the fault properties, e.g., weakening or strengthening, and the properties of the loading system. Here, we investigate the poorly understood effect of energy-flux to the fault that should equal or exceed the energy-dissipation-rate along the slipping fault. We explore the relationship between energy-flux and slip style in shear experiments along granite and diorite laboratory faults, during which the faults were subjected to controlled energy-flux, and responded spontaneously to it. The monitored evolution of slip-velocity, shear stress, and slip-distance revealed three slip styles that depend on the applied energy-flux: (1) steady-state slip; (2) spontaneous creep events of small displacement with negligible weakening; and (3) spontaneous, unstable events with slip-velocities up to 0.8 m/s, slip-distances up to 0.5 m, and stress-drops up to 1 MPa, which are comparable to observed values of moderate earthquakes. These slip styles are similar in character to those observed along natural faults. We further propose that the rate of energy flow from crustal blocks can control the slip velocity during earthquakes.

Nanocrystalline Principal Slip Zones and Their Role in Controlling Crustal Fault Rheology

Principal slip zones (PSZs) are narrow (<10 cm) bands of localized shear deformation that occur in the cores of upper-crustal fault zones where they accommodate the bulk of fault displacement. Natural and experimentally-formed PSZs consistently show the presence of nanocrystallites in the <100 nm size range. Despite the presumed importance of such nanocrystalline (NC) fault rock in controlling fault mechanical behavior, their prevalence and potential role in controlling natural earthquake cycles remains insufficiently investigated. In this contribution, we summarize the physical properties of NC materials that may have a profound effect on fault rheology, and we review the structural characteristics of NC PSZs observed in natural faults and in experiments. Numerous literature reports show that such zones form in a wide range of faulted rock types, under a wide range of conditions pertaining to seismic and a-seismic upper-crustal fault slip, and frequently show an internal crystallographic preferred orientation (CPO) and partial amorphization, as well as forming glossy or “mirror-like” slip surfaces. Given the widespread occurrence of NC PSZs in upper-crustal faults, we suggest that they are of general significance. Specifically, the generally high rates of (diffusion) creep in NC fault rock may play a key role in controlling the depth limits to the seismogenic zone.

Diagenetic and shear-induced transitions of frictional strength of carbon-bearing faults and their implications for earthquake rupture dynamics in subduction zones

Subduction-related diagenetic reactions affect fault strength and are thus important for understanding earthquake rupture dynamics in subduction zones. Carbonaceous material (CM) is found worldwide in active plate-boundary and intracontinental faults, yet the effect of its transformation on frictional strength and rupture dynamics remains unknown. We conducted high-velocity friction experiments together with organochemical analyses on CM in the form of lignite, bituminous coal, anthracite and graphite. Results clearly show that an increase in CM maturity and crystallinity leads to a decrease in the peak friction coefficient (from 0.5 to 0.2). We also infer that friction applied to low-grade CM increases its maturity, but friction applied to high-grade CM reduces its maturity. These findings suggest that both diagenetic and shear-induced transformations of CM strongly affect the frictional strength of CM-bearing faults, potentially affecting the depth-dependences of frictional strength and rupture dynamics on plate-subduction faults.

Mechanical behaviour of fluid-lubricated faults

Fluids are pervasive in fault zones cutting the Earth's crust; however, the effect of fluid viscosity on fault mechanics is mainly conjectured by theoretical models. We present friction experiments performed on both dry and fluid-permeated silicate and carbonate bearing-rocks, at normal effective stresses up to 20 MPa, with a slip-rate ranging between 10 μm/s and 1 m/s. Four different fluid viscosities were tested. We show that both static and dynamic friction coefficients decrease with viscosity and that dynamic friction depends on the dimensionless Sommerfeld number (S) as predicted by the elastohydrodynamic-lubrication theory (EHD).Under favourable conditions (depending on the fluid viscosity (η), co-seismic slip-rate (V), fault geometry (L/H₀²) and earthquake nucleation depth (∝σ_eff)), EHD might be an effective weakening mechanism during natural and induced earthquakes. However, at seismic slip-rate, the slip weakening distance (D_c) increases markedly for a range of fluid viscosities expected in the Earth, potentially favouring slow-slip rather than rupture propagation for small to moderate earthquakes.

The effect of fluid viscosity on fault mechanics is mainly conjectured by theoretical models. Here, the authors present experimental data from rock friction experiments, showing both static and dynamic friction coefficients to decrease with viscosity and dynamic friction to depend on the Sommerfeld number.

Dynamic viability of the 2016 Mw 7.8 Kaikōura earthquake cascade on weak crustal faults

We present a dynamic rupture model of the 2016 M_w 7.8 Kaikōura earthquake to unravel the event's riddles in a physics-based manner and provide insight on the mechanical viability of competing hypotheses proposed to explain them. Our model reproduces key characteristics of the event and constraints puzzling features inferred from high-quality observations including a large gap separating surface rupture traces, the possibility of significant slip on the subduction interface, the non-rupture of the Hope fault, and slow apparent rupture speed. We show that the observed rupture cascade is dynamically consistent with regional stress estimates and a crustal fault network geometry inferred from seismic and geodetic data. We propose that the complex fault system operates at low apparent friction thanks to the combined effects of overpressurized fluids, low dynamic friction and stress concentrations induced by deep fault creep.

Emergence of self-affine surfaces during adhesive wear

Friction and wear depend critically on surface roughness and its evolution with time. An accurate control of roughness is essential to the performance and durability of virtually all engineering applications. At geological scales, roughness along tectonic faults is intimately linked to stick-slip behaviour as experienced during earthquakes. While numerous experiments on natural, fractured, and frictional sliding surfaces have shown that roughness has self-affine fractal properties, much less is known about the mechanisms controlling the origins and the evolution of roughness. Here, by performing long-timescale molecular dynamics simulations and tracking the roughness evolution in time, we reveal that the emergence of self-affine surfaces is governed by the interplay between the ductile and brittle mechanisms of adhesive wear in three-body contact, and is independent of the initial state.

Surface roughness evolution with time is key for tribological applications. Here, the authors demonstrate by numerical simulations the evolution of sliding surfaces into self-affine morphologies during adhesive wear due to the formation of a third body trapped at the interface.

Empirical investigation of friction weakening of terrestrial and Martian landslides using discrete element models

Understanding what controls the travelling distance of large landslides has been the topic of considerable debate. By combining observation and experimental data with depth-averaged continuum modelling of landslides and generated seismic waves, it was empirically observed that lower effective friction had to be taken into account in the models to reproduce the dynamics and runout distance of larger volume landslides. Moreover, such simulation and observation results are compatible with a friction weakening with velocity as observed in earthquake mechanics. We investigate here as to whether similar empirical reduced friction should be put into discrete element models (DEM) to reproduce observed runout of large landslides on Earth and on Mars. First we show that, in the investigated parameter range and for a given volume, the runout distance simulated by 3D DEM is not much affected by the number (i.e. size) of grains once this number attains ~ 8000. We then calibrate the model on laboratory experiments and simulate other experiments of granular flows on inclined planes, making it possible for the first time to reproduce the observed effect of initial volume and aspect ratio on runout distances. In particular, the normalised runout distance starts to depend on the volume involved only above a critical slope angle > 16-19°, as observed experimentally. Finally, based on field data (volume, topography, deposit), we simulate a series of landslides on simplified inclined topography. The empirical friction coefficient, calibrated to reproduce the observed runout for each landslide, is shown to decrease with increasing landslide volume (or velocity), going down to values as low as 0.1-0.2. No distinguishable difference is observed between the behaviour of terrestrial and Martian landslides.

Earthquake lubrication and healing explained by amorphous nanosilica

During earthquake propagation, geologic faults lose their strength, then strengthen as slip slows and stops. Many slip-weakening mechanisms are active in the upper-mid crust, but healing is not always well-explained. Here we show that the distinct structure and rate-dependent properties of amorphous nanopowder (not silica gel) formed by grinding of quartz can cause extreme strength loss at high slip rates. We propose a weakening and related strengthening mechanism that may act throughout the quartz-bearing continental crust. The action of two slip rate-dependent mechanisms offers a plausible explanation for the observed weakening: thermally-enhanced plasticity, and particulate flow aided by hydrodynamic lubrication. Rapid cooling of the particles causes rapid strengthening, and inter-particle bonds form at longer timescales. The timescales of these two processes correspond to the timescales of post-seismic healing observed in earthquakes. In natural faults, this nanopowder crystallizes to quartz over 10s–100s years, leaving veins which may be indistinguishable from common quartz veins.

Tectonic faults weaken during slip in order to accelerate and produce earthquakes. Here the authors show a mechanism for weakening faults through the transformation of quartz to amorphous nanoparticulate wear powders that lubricate friction experiments, and transform back to quartz under geologic conditions.

Bimodal seismicity in the Himalaya controlled by fault friction and geometry

There is increasing evidence that the Himalayan seismicity can be bimodal: blind earthquakes (up to Mw ~ 7.8) tend to cluster in the downdip part of the seismogenic zone, whereas infrequent great earthquakes (Mw 8+) propagate up to the Himalayan frontal thrust. To explore the causes of this bimodal seismicity, we developed a two-dimensional, seismic cycle model of the Nepal Himalaya. Our visco-elasto-plastic simulations reproduce important features of the earthquake cycle, including interseismic strain and a bimodal seismicity pattern. Bimodal seismicity emerges as a result of relatively higher friction and a non-planar geometry of the Main Himalayan Thrust fault. This introduces a region of large strength excess that can only be activated once enough stress is transferred upwards by blind earthquakes. This supports the view that most segments of the Himalaya might produce complete ruptures significantly larger than the 2015 Mw 7.8 Gorkha earthquake, which should be accounted for in future seismic hazard assessments.

There is increasing evidence that the seismicity of large Himalayan earthquakes can be bimodal, but the underlying mechanisms are unclear. Here, the authors present a model and show that the bimodal seismicity results from a relatively higher friction and a non-planar geometry of the Himalayan megathrust.

A record of changes in the Gran Sasso groundwater before, during and after the 2016 Amatrice earthquake, central Italy

We performed continuous recordings (May 2015 – January 2017) of hydraulic pressure and electrical conductivity of groundwater in the 190 m-long horizontal S13 borehole drilled next to the deep underground laboratories of Gran Sasso (LNGS-INFN), located in the core of the Gran Sasso carbonate aquifer (central Italy) at a distance of about 39 km south-eastward from the 24 August 2016 Amatrice earthquake (6.0 M_w) epicenter. Using a 3-channel, 24-bit ADC we achieved a sampling rate of groundwater physical properties up to 50 Hz for each channel. We focused on the analysis of data recorded before, during and after the Amatrice earthquake, describing and discussing in detail the evidence for significant hydraulic pressure and electrical conductivity anomalies recorded before the main shock. We identified unambiguous signals in the hydraulic pressure data starting on 19 August, i.e. five days before the 24 August mainshock. A more careful analysis allowed us to detect the inception of a weak change up to 40 days before the Amatrice earthquake and a significant variation in the electrical conductivity data about 60 days before. The data revealed highly dynamic aquifer behaviour associated with the uprising of geogas probably related to the preparation stage of the Amatrice earthquake.

Friction of sea ice

Static and kinetic friction play a fundamental role in sea-ice mechanics. The coefficient of static friction increases with hold time under normal load and is modelled in terms of creep and fracture of asperities in contact. The coefficient of kinetic friction exhibits velocity strengthening at lower speeds and velocity weakening at intermediate speeds. Strengthening is modelled in terms of asperity creep and hardness; weakening is modelled in terms of a progressive increase in the true area of contact wetted by meltwater produced through frictional heating. The concept is introduced of contact size distribution in which the smallest contacts melt first, leading to the onset of weakening; the largest melt last, leading to a third regime of kinetic friction and again to strengthening where hydrodynamics governs. Neither the static nor the kinetic coefficient is significantly affected by the presence of sea water. The paper closes with a few implications for sea-ice mechanics. The paper is based largely upon a critical review of the literature, but includes a more quantitative, physics-based analysis of velocity strengthening and a new analysis of velocity weakening that incorporates parameters that describe the (proposed) fractal character of the sliding interface.This article is part of the theme issue 'Modelling of sea-ice phenomena'.

Strain localization in dry sheared granular materials: A compactivity-based approach

Shear banding is widely observed in natural fault zones as well as in laboratory experiments on granular materials. Understanding the dynamics of strain localization under different loading conditions is essential for quantifying strength evolution of fault gouge and energy partitioning during earthquakes and characterizing rheological transitions and fault zone structure changes. To that end, we develop a physics-based continuum model for strain localization in sheared granular materials. The grain-scale dynamics is described by the shear transformation zone (STZ) theory, a nonequilibrium statistical thermodynamic framework for viscoplastic deformation in amorphous materials. Using a finite strain computational framework, we investigate the initiation and growth of complex shear bands under a variety of loading conditions and identify implications for strength evolution and the ductile to brittle transition. Our numerical results show similar localization patterns to field and laboratory observations and suggest that shear zones show more ductile response at higher confining pressures, lower dilatancy, and loose initial conditions. Lower pressures, higher dilatancy, and dense initial conditions favor a brittle response and larger strength drops. These findings shed light on a range of mechanisms for strength evolution in dry sheared granular materials and provide a critical input to physics-based multiscale models of fault zone instabilities.

Rapid migration of CO2-rich micro-fluids in calcite matrices

The transport of supercritical fluids is a determining factor for several geological processes and fundamental in predicting natural resource accumulation and distribution. Calcite, ubiquitous in most geological environments, may contain supercritical CO2 trapped under the form of fluid inclusions that may move through grain boundaries affecting the rock physical properties. However, despite macroscopic evidence for this process, until recent it was not possible to characterize this process at the nano-scale due to the difficulty of such observations. In this study, we report nanometer-scale observations on calcite crystal surfaces and demonstrate that stress with absence of visible deformation produces fluid leakage from fluid inclusions. Atomic Force Microscopy scanning experiments on freshly cleaved calcite crystals containing visible fluid inclusions revealed the spontaneous formation of nanometer-scale hillocks on flat crystal terraces in only a few minutes, without evidence of surface dissolution. The fact the hillocks formed on flat surface in a short time was unexpected and suggests deposition of material from the inner crystal to the surface through small-scale fluid migration. We estimated the rate of this fluid mobility is by several orders of magnitude higher than the diffusion rate through vacancies estimated in calcite crystals showing that CO2-rich fluids through micro-pore and nano-pore spaces is in reality much higher than previously assumed using current predictive models.

Carbonaceous Materials in the Fault Zone of the Longmenshan Fault Belt: 2. Characterization of Fault Gouge from Deep Drilling and Implications for Fault Maturity

In recent works on the determination of graphitization of carbonaceous materials (CM) within the principal slip zone (PSZ) of the Longmenshan fault (China), we demonstrated that the formation of graphite, resulted from strain and frictional heating, could be evidence of past seismic slip. Here we utilize Raman Spectroscopy of CM (RSCM) on the CM-bearing gouges in the fault zone of the Longmenshan fault belt, at the borehole depth of 760 m (FZ760) from the Wenchuan earthquake Fault Scientific Drilling project-1 (WFSD-1), to quantitatively characterize CM and further retrieve ancient fault deformation information in the active fault. RSCM shows that graphitization of CM is intense in the fault core with respect to the damage zone, with the graphitized carbon resembling those observed on experimentally formed graphite that was frictionally generated. Importantly, compared to the recognized active fault zone of the Longmenshan fault, the RSCM of measured CM-rich gouge shows a higher degree of graphitization, likely derived from high-temperature-perturbation faulting events. It implies that FZ760 accommodated numerous single-event displacement and/or at higher normal stresses and/or in the absence of pore fluid and/or along a more localized slip surface(s). Because graphite is a well-known lubricant, we surmise that the presence of the higher degree graphitized CM within FZ760 will reduce the fault strength and inefficiently accumulate tectonic stress during the seismic cycle at the current depth, and further infer a plausible mechanism for fault propagation at the borehole depth of 590 m during the Mw 7.9 Wenchuan earthquake.

Carbonaceous Materials in the Fault Zone of the Longmenshan Fault Belt: 1. Signatures within the Deep Wenchuan Earthquake Fault Zone and Their Implications

Graphitization of carbonaceous materials (CM) has been experimentally demonstrated as potential evidence of seismic slip within a fault gouge. The southern segment of the Longmenshan fault, a CM-rich-gouge fault, accommodated coseismic slip during the 2008 Mw 7.9 Wenchuan earthquake and potentially preserves a record of processes that occurred on the fault during the slip event. Here, we present a multi-technique characterization of CM within the active fault zone of the Longmenshan fault from the Wenchuan earthquake Fault Scientific Drilling-1. By contrast with field observations, graphite is pervasively and only distributed in the gouge zone, while heterogeneously crystallized CM are present in the surrounding breccia. The composite dataset that is presented, which includes the localized graphite layer along the 2008 Wenchuan earthquake principal slip zone, demonstrates that graphite is widely distributed within the active fault zone. The widespread occurrence of graphite, a seismic slip indicator, reveals that surface rupturing events commonly occur along the Longmenshan fault and are characteristic of this tectonically active region.

Dynamic earthquake rupture preserved in a creeping serpentinite shear zone

Laboratory experiments on serpentinite suggest that extreme dynamic weakening at earthquake slip rates is accompanied by amorphisation, dehydration and possible melting. However, hypotheses arising from experiments remain untested in nature, because earthquake ruptures have not previously been recognised in serpentinite shear zones. Here we document the progressive formation of high-temperature reaction products that formed by coseismic amorphisation and dehydration in a plate boundary-scale serpentinite shear zone. The highest-temperature products are aggregates of nanocrystalline olivine and enstatite, indicating minimum peak coseismic temperatures of ca. 925 ± 60 °C. Modelling suggests that frictional heating during earthquakes of magnitude 2.7–4 can satisfy the petrological constraints on the coseismic temperature profile, assuming that coseismic fluid storage capacity and permeability are increased by the development of reaction-enhanced porosity. Our results indicate that earthquake ruptures can propagate through serpentinite shear zones, and that the signatures of transient frictional heating can be preserved in the fault rock record.

Creeping serpentinite shear zones may host large earthquakes, but direct evidence of frictional heating and rupture have been missing. Here, the authors demonstrate via laboratory experiments that earthquake ruptures can propagate through serpentinite shear zones shown by high-temperature reaction products.