Fault Friction Research Articles

Implications of a Reverse Polarity Earthquake Pair on Fault Friction and Stress Heterogeneity Near Ridgecrest, California

AbstractWe apply the Matrix Profile algorithm to 100 days of continuous data starting 10 days before the 2019 M 6.4 and M 7.1 Ridgecrest earthquakes from borehole seismic station B921 near the Ridgecrest aftershock sequence. We identify many examples of reversely polarized waveforms, but focus on one particularly striking earthquake pair with strongly negatively correlated P and S waveforms at B921 and several other nearby stations. Waveform‐cross‐correlation‐based relocation of these events indicates they are at about 10 km depth and separated by only 115 m. Individual focal mechanisms are poorly resolved for these events because of the limited number of recording stations with unambiguous P polarities. However, relative P and S polarity and amplitude information can be used to constrain the likely difference in fault plane orientation between the two events to be 5–20°. We explore possible models to explain these observations, including low effective coefficients of fault friction and short‐wavelength stress heterogeneity caused by prior earthquakes. Although definitive conclusions are lacking, we favor local stress heterogeneity as being more consistent with other observations for the Ridgecrest region.

Does the direct effect of friction increase continuously with absolute temperature?

Constitutive models of fault friction form the basis of physics-based simulations of seismic activity. A generally accepted framework for the slip-rate and state dependence of friction involves a thermally activated process, whereby the probability of slip along microasperities adheres to an Arrhenius law. This model, which has become widely adopted among experimentalists and theoreticians, predicts a continuous increase of the direct effect with absolute temperature, but is it observed experimentally? Leveraging comprehensive laboratory data across diverse hydrothermal, barometric, and lithological conditions, we demonstrate that, contrary to the classical view, the direct effect for a given deformation mechanism remains largely temperature-independent. Instead, the incremental shifts in the direct effect often coincide with the brittle to semi-brittle transition, across which distinct deformation mechanisms operate. These considerations challenge the validity of the classical model. Realistic constitutive laws for rock failure within the lithosphere must incorporate the contributions of multiple deformation mechanisms within active fault zones.

Frictional Strength and Frictional Instability of Glaucophane Gouges at Blueschist Temperatures Support Diverse Modes of Fault Slip From Slow Slip Events to Moderate‐Sized Earthquakes

AbstractFluid overpressure from the water released by subducted sediments and oceanic crust is an important mechanism for generating earthquakes via brittle failure and frictional instability. If unstable, such fault materials may also host diverse fault reactivation mechanisms from slow slip events to moderate‐sized earthquakes in cold subduction zones. We examine this hypothesis for glaucophane gouge ‐ a key index minerals for blueschist facies ‐ at lower confining stresses where behavior is poorly understood. We measure friction and stability at temperatures of 100°C–500°C and effective normal stresses of 50–200 MPa, to explore the controls of temperature, stresses and excess pore fluid pressures on fault friction. The frictional coefficient of glaucophane at representative temperatures and stresses is ∼0.70, insensitive to temperature but with a slight increase at lower effective stresses. Elevating temperature promotes a transition from velocity‐strengthening to weak velocity‐weakening behavior, indicating the destabilizing effect of high temperature downdip in subduction zones. Reducing effective normal stress or concomitantly elevating pore fluid pressure further strengthens the velocity‐weakening response and would be manifest as moderate‐sized earthquakes. Our results support the potential for enhanced unstable sliding of glaucophane gouges at lower effective stresses and blueschist facies temperatures ‐ with weak to moderate velocity‐weakening responses upon fluid pressurization vital for understanding the abundance of slow slip to moderate‐sized earthquakes apparent in cold subduction zones.

Unstable Sliding of Plagioclase Gouge and Deformation Mechanisms Under Hydrothermal Conditions With Effective Normal Stresses of 100–300 MPa

AbstractPlagioclase feldspar is a major mineral in mafic crustal rocks. To better understand the deformation mechanism of plagioclase feldspar during frictional faulting, we conducted shearing experiments on simulated plagioclase gouge in a wide range of effective normal stress of 100–300 MPa, pore‐water pressure of 30–100 MPa, and temperatures ranging from 100°C to 600°C. The coefficient of friction is found to range from 0.65 to 0.74 across the entire temperature range, showing no significant thermal weakening process. Except for a case at 200°C with an effective normal stress of 300 MPa, the frictional sliding is velocity weakening over the whole temperature range, showing a steady‐state rate dependence (a−b) ranging from −0.5 × 10−3 to −8.6 × 10−3. This property facilitates nucleation of unstable slips in frictional faulting. Above 200°C, the direct rate effect parameter (a) and the evolution effect parameter (b) of friction increase with temperature up to a threshold of 400°C or 500°C, depending on the effective normal stress. This thermal enhancement suggests thermally activated creep at contact junctions governed by intergranular pressure solution, as evidenced by microstructural signatures indicating the prevalence of very fine precipitates formed at the surfaces of gouge particles as a result of pressure solution. In frictional sliding of plagioclase, a low effective normal stress of 100 MPa corresponds to a higher degree of velocity weakening and tends to facilitate seismic slip rather than slow slips, whereas the high effective normal stress of 300 MPa corresponds to a minor velocity weakening which may cause slow‐slip events in faults of limited size.

Adjoint-based inversion for stress and frictional parameters in earthquake modeling

We present an adjoint-based optimization method to invert for stress and frictional parameters used in earthquake modeling. The forward problem is linear elastodynamics with nonlinear rate-and-state frictional faults. The misfit functional quantifies the difference between simulated and measured particle displacements or velocities at receiver locations. The misfit may include windowing or filtering operators. We derive the corresponding adjoint problem, which is linear elasticity with linearized rate-and-state friction and, for forward problems involving fault normal stress changes, nonzero fault opening, with time-dependent coefficients derived from the forward solution. The gradient of the misfit is efficiently computed by convolving forward and adjoint variables on the fault. The method thus extends the framework of full-waveform inversion to include frictional faults with rate-and-state friction. In addition, we present a space-time dual-consistent discretization of a dynamic rupture problem with a rough fault in antiplane shear, using high-order accurate summation-by-parts finite differences in combination with explicit Runge–Kutta time integration. The dual consistency of the discretization ensures that the discrete adjoint-based gradient is the exact gradient of the discrete misfit functional as well as a consistent approximation of the continuous gradient. Our theoretical results are corroborated by inversions with synthetic data. We anticipate that adjoint-based inversion of seismic and/or geodetic data will be a powerful tool for studying earthquake source processes; it can also be used to interpret laboratory friction experiments.

Emergent fault friction and supershear in a continuum model of geophysical rupture

Important physical observations in rupture dynamics such as static fault friction, short-slip, self-healing, and the supershear phenomenon in cracks are studied. A continuum model of rupture dynamics is developed using the field dislocation mechanics (FDM) theory. The energy density function in our model encodes accepted and simple physical facts related to rocks and granular materials under compression. We work within a 2-dimensional ansatz of FDM where the rupture front is allowed to move only in a horizontal fault layer sandwiched between elastic blocks. Damage via the degradation of elastic modulus is allowed to occur only in the fault layer, characterized by the amount of plastic slip. The theory dictates the evolution equation of the plastic shear strain to be a Hamilton–Jacobi (H-J) equation, resulting in the representation of a propagating rupture front. A Central-Upwind scheme is used to solve the H-J equation. The rupture propagation is fully coupled to elastodynamics in the whole domain, and our simulations recover static friction laws as emergent features of our continuum model, without putting in by hand any such discontinuous criteria in the formulation. Estimates of material parameters of cohesion and friction angle are deduced. Short-slip and slip-weakening (crack-like) behaviors are also reproduced as a function of the degree of damage behind the rupture front. The long-time behavior of a moving rupture front is probed, and it is deduced that equilibrium profiles under no shear stress are not traveling wave profiles under non-zero shear load in our model. However, it is shown that a traveling wave structure is likely attained in the limit of long times. Finally, a crack-like damage front is driven by an initial impact loading, and it is observed in our numerical simulations that an upper bound to the crack speed is the dilatational wave speed of the material unless the material is put under pre-stressed conditions, in which case supersonic motion can be obtained. Without pre-stress, intersonic (supershear) motion is recovered under appropriate conditions.

Relocation of the 2021 MW 7.4 Maduo, Qinghai, China earthquake sequence and implications for seismogenic structure

An MW 7.4 earthquake struck Maduo County, Qinghai, China on 22 May 2021. To better understand the seismogenic structure of this region, we collect local earthquake arrival time data at the MAD station in the China Earthquake Networks Center observational bulletins during 1 June to 20 September 2021, and manually pick P and S wave arrival times from high-quality seismograms recorded at 34 recently deployed MaduoArray portable seismic stations. Using these arrival times, we relocate the Maduo earthquake sequence using earthquake association, absolute location and relative location methods. Our results show that the 2021 Maduo earthquake sequence occurred along the Kunlun Mountain Pass-Jiangcuo fault zone in the NWW-SEE direction and the aftershocks are located on both sides of the mainshock, showing characteristics of bilateral rupture. Vertical cross-sections of the aftershock distribution illustrate a nearly vertical shape of the seismogenic fault that tilts toward the northeast and southwest in different sections, reflecting a complex geometry of the fault plane. There is a horsetail bifurcation phenomenon at the eastern end of the fault zone. A sparse area of aftershocks appears at about 30 km east of the mainshock epicenter, which may be associated with a uniform fault friction and sufficient release of rupture energy caused by super-shear rupture of the mainshock. Taking into account many geophysical results including seismic tomography and magnetotelluric soundings, we speculate that the occurrence of the Maduo earthquake could be affected by crustal fluids in the fault zone. The fluids may ascend from the lower crustal flow beneath northeastern Tibet.

Stable Recovery of Coefficients in an Inverse Fault Friction Problem

We consider the inverse fault friction problem of determining the friction coefficient in the Tresca friction model, which can be formulated as an inverse problem for differential inequalities. We show that the measurements of elastic waves during a rupture uniquely determine the friction coefficient at the rupture surface with explicit stability estimates.

Frictional Fault Strength Analysis of Palu-Koro and Matano Faults, Sulawesi, Indonesia, from Earthquake Focal Mechanism Data

ABSTRACT Sulawesi’s complex tectonic setting generates active Palu-Koro and Matano faults. Understanding earthquake mechanic is important for mitigation purpose. We performed fault strength analysis by examining its friction, through joint stress inversion and stress magnitude estimation approach, using focal mechanism data from International Seismological Centre bulletin. Friction was analyzed in four segments of the Palu-Koro Fault and two segments of the Matano Fault. Bootstrap and Jack-knife are applied to assess its uncertainty. From north to south, the friction pattern increases at the Palu-Koro Fault segments: Donggala (0.35<μs <0.56) to Palu segment (0.59<μs <0.76); then decreases sequentially: at Saluki (0.49<μs <0.55), Moa segment (0.31< μs <0.35), and further at the Matano Fault segment (Pamsoa-Ballawai) (0.26<μs <0.3). Regarding the 2018 Mw 7.5 Palu earthquake rupture, high friction at the Saluki segment coincides with and indicates its rupture termination. Seismicity records from 1977 to 2011 showed that the largest event occurred at Matano Fault has Mw less than 7.5. Low friction observed at Matano Fault correlates with small magnitude events and low fault strength.

Impact of fault slip induced by water injection on casing deformation: a numerical study

Fluid injection can cause fault reactivation and slip due to weakened cementation strength. The resulting fault slippage poses significant challenges to wellbore and casing integrity. To address this, we developed a multi-scale 3D finite element model, integrating reservoir, and near-wellbore components and applied to an offshore well in China’s ZD Reservoir, an area that has experienced fault reactivations (and casing deformations of wells intersecting reactivated faults) due to routinely water injection in nearby wells. Based on the actual formation’s wellbore-fault angle, a full-size fault is inserted into the model. The traction–separation law (TSL) describes progressive fault damage induced by waterflooding. The model also simulates the impact of fault slippage on cement failure and casing deformation. Research findings indicate that fluid injection reduces friction and shear cementation strength between faults, leading to reactivated fault slippage. The slip process involves three stages: low-level, sharply rising, and high-level stages, each exhibiting reduced fault friction and shear cementation strength. These outcomes provide insights to prevent casing damage during water flooding, optimizing injection points, or using better cement recipes at minimal cost.

Earthquake Cycle Mechanics During Caldera Collapse: Simulating the 2018 Kı̄lauea Eruption

AbstractIn multiple observed caldera‐forming eruptions, the rock overlying a draining magma reservoir dropped downward along ring faults in sequences of discrete collapse earthquakes. These sequences are analogous to tectonic earthquake cycles and provide opportunities to examine fault mechanics and collapse eruption dynamics over multiple events. Collapse earthquake cycles have been studied with zero‐dimensional slider‐block models, but these do not account for the complicated interplay between fluid and elastic dynamics or for factors such as the heterogeneous fault properties and non‐vertical ring fault geometries often inferred at volcanoes. We present two‐dimensional axisymmetric mafic piston‐like collapse earthquake cycle models that include rate‐and‐state friction, fully‐dynamic elasticity, and compressible viscous fluid magma flow. We demonstrate that collapse earthquake intervals and magnitudes are highly sensitive to inertial effects, evolving stress fields, fault geometry, and depth‐varying fault friction. Given the consistent earthquake cycles observed in most eruptions, this suggests that ring faults can quickly stabilize and often become nearly vertical at depth. We use the well‐monitored 2018 collapse sequence at Kı̄lauea as a case study. Our model can produce many features of Kı̄lauea seismic and geodetic observations, except for a significant amount of interseismic slip, which cannot be readily explained with simple rate‐and‐state friction parameterizations.

Post-mining related reactivation potential of faults hosted in tight reservoir rocks around flooded coal mines, eastern Ruhr Basin, Germany

The cessation of hard coal mining in the Ruhr Basin in 2018 marked the region's transition to the post-mining phase. Controlled mine water rebound induces changes in the subsurface stress conditions, as pore pressure increases locally. Presently, mine water rebound is observed in the eastern Ruhr Basin (water province “Haus Aden”) along with associated microseismicity. Furthermore, post-mining challenges might comprise the potential risk of fault reactivation, which is addressed in this study by conducting a fault slip assessment. Based on subsurface coal seam mapping data, a 3D structural model for the NE part of the “Haus Aden” water province has been constructed to serve as the basis for identifying the most vulnerable fault trends and types of the structural inventory. Slip tendency analysis, considering normal faulting conditions, revealed NW-SE to NNW-SSE trending normal faults to be most susceptible to reactivation. Probabilistic fault slip assessment, focused on NW-SE to NNW-SSE trending normal faults mapped within the “Heinrich-Robert” colliery, show no fault reactivation potential for a mine water rebound up to a level of 640 m below ground. Assuming hydrostatic conditions in the vicinity of the faults, friction coefficients are only partially exceeded for high differential stresses. In addition, a novel workflow is used to model the spatial variability of the frictional fault strength as input for a fault stability analysis, exemplified for a selected NNW-SSE trending normal fault. For considering hydrostatic pore pressure, results show that the fault consists mainly of stable, but also unstable, horizontally elongated patches. These findings question the conventional simplified approach of using a single constant friction coefficient for fault stability analysis.

Assessment of borehole breakouts from acoustic image log and its geomechanical implications — A case study from Triassic-Ordovician interval of the Berkaoui Field, southeastern Algeria

In this study, we interpret a cumulative 600 m acoustic image log across the Triassic to Cambro-Ordovician interval in the Berkaoui oil field, Algeria. We interpret 40 distinct breakout zones that have a combined length of 210 m. These breakouts are aligned in the north-northeast–south-southwest direction, indicating a mean maximum horizontal stress (SHmax) azimuth of 110°N. The observed breakouts are ranked as “A-Quality” following the World Stress Map ranking guidelines. The angular width of each breakout has been inferred from the image log analysis and the same has been used to infer the SHmax gradient by stress polygon approach following the frictional faulting mechanism. The stress polygon across all the breakout intervals provides a practical Shmax range between 24.7 MPa/km and 31.1 MPa/km, with an average gradient of approximately 27 MPa/km. Considering the Shmin range across the studied intervals, we infer a SHmax/Shmin ratio dominantly between 1.40 and 1.65, which is a much narrower and better-constrained range when compared with the previously published ranges from nearby fields with the same stratigraphy. The relative magnitudes of the in situ stresses indicate a strike-slip faulting regime in the Berkaoui Field. This study presents the utility of image log analysis and the integration of breakout interpretation to obtain a more robust geomechanical model with reduced SHmax uncertainty.

A proxy implementation of thermal pressurization for earthquake cycle modelling on rate-and-state faults

SUMMARY The reduction of effective normal stress during earthquake slip due to thermal pressurization of fault zone pore fluids is a significant fault weakening mechanism. Explicit incorporation of this process into frictional fault models involves solving the diffusion equations for fluid pressure and temperature outside the fault at each time step, which significantly increases the computational complexity. Here, we propose a proxy for thermal pressurization implemented through a modification of the rate-and-state friction law. This approach is designed to emulate the fault weakening and the relationship between breakdown energy and slip resulting from thermal pressurization and is appropriate for fully dynamic simulations of multiple earthquake cycles. It preserves the computational efficiency of conventional rate-and-state friction models, which in turn can enable systematic studies to advance our understanding of the effects of fault weakening on earthquake mechanics. In 2.5-D simulations of pulse-like ruptures on faults with finite seismogenic width, based on our thermal pressurization proxy, we find that the spatial distribution of slip velocity near the rupture front is consistent with the conventional square-root singularity, despite continued slip-weakening within the pulse, once the rupture has propagated a distance larger than the rupture width. An unconventional singularity appears only at shorter rupture distances. We further derive and verify numerically a theoretical estimate of the breakdown energy dissipated by our implementation of thermal pressurization. These results support the use of fracture mechanics theory to understand the propagation and arrest of very large earthquakes.

Numerical simulations of earthquake rupture induced by pressure perturbation

The subsurface fluid injection can cause pressure increase within faults, leading to earthquake occurrences. However, the factors controlling earthquake rupture due to pressure perturbation remain poorly understood. To resolve this problem, we simulate the physical processes of earthquake nucleation and rupture on strike-slip faults perturbated by pressure migration based on the slip-weakening law. Multiple kinds of factors, including background stress, fluid injection rates, the area of the pressurized region, fault geometry, and fault friction coefficients, are considered in our simulations. Our simulation results reveal that the ratio of shear stress to normal stress rather than their absolute values controls the rupture behavior. With the large stress ratios, high injection rates, and large pressurized areas, earthquakes are prone to propagate as runaway ruptures. Additionally, faults with large aspect ratios of length to width are also favorable for causing runaway ruptures. In contrast, the factors of fault strike, dip angles and friction coefficients have minor influence on rupture behavior.

Constraining maximum horizontal stress using wellbore breakouts — A case study from the Ordovician tight reservoir of the northeastern Oued Mya Basin, Algeria

In this study, we interpret the maximum horizontal stress (SHmax) azimuth from the breakout positions of the wellbore and attempt to constrain the SHmax gradient based on the interpreted breakout width. A cumulative of 110 m of breakouts are deciphered within the Ordovician Hamra Quartzite interval of the Oued Mya Basin from a 138 m acoustic image log. These breakouts are ranked as “A-Quality” following the World Stress Map ranking guidelines. We infer a mean SHmax orientation of N28°E ± 8°. Following the frictional faulting mechanism and stress polygon approach, measurement of the minimum horizontal stress (Shmin) from minifrac tests and observations of the compressive failures from the acoustic image log provide strong constraints on the SHmax magnitude in the reservoir interval in the absence of core-measured rock strength. Interpreted breakout widths exhibit a range between 32.6° and 90.81°, which indicates a SHmax range of 24.4–34.7 MPa/km. The average breakout width of 62.58° translates to a narrower SHmax gradient range, varying between 27.2 and 31.2 MPa/km. The relative magnitudes of the principal stresses indicate a strong strike-slip tectonic stress state. Considering all the uncertainties, we infer a SHmax/Shmin ratio of 1.41–1.81 within the Ordovician interval.

Rupture Jumping and Seismic Complexity in Models of Earthquake Cycles for Fault Stepovers with Off-Fault Plasticity

ABSTRACT Fault stepovers are prime examples of geometric complexity in natural fault zones that may affect seismic hazard by determining whether an earthquake rupture continues propagating or abruptly stops. However, the long-term pattern of seismicity near-fault stepovers and underlying mechanisms of rupture jumping in the context of earthquake cycles are rarely studied. Leveraging a hybrid numerical scheme combining the finite element and the spectral boundary integral methods, FEBE, we carry out fully dynamic simulations of sequences of earthquakes and aseismic slip for both compressive and tensile stepovers with off-fault plasticity. We consider a rate-and-state friction law for the fault friction and pressure-sensitive Drucker–Prager plasticity for the off-fault bulk response. We observe that the accumulation of plastic deformation, an indication of off-fault damage, is significantly different in the two cases, with more plastic deformation projected in the overlapping region for the tensile stepover. The seismic pattern for a tensile stepover is more complex than for a compressive stepover, and incorporating plasticity also increases complexity, relative to the elastic case. A tensile stepover with off-fault plasticity shows rupture segmentation, temporal clustering, and frequent rupture jumping from one fault to another. These results shed light on possible mechanisms of rupture jumping in fault stepovers as well as the long-term evolution of the fault zone.

Frictional Aging of Single‐Asperity Nanoindentation Contacts in Quartz and Calcite

AbstractThe evolution of fault friction during the interseismic period affects the mechanics of a future earthquake on the same fault patch. Frictional aging has been previously tied to time‐dependent contact area growth through observations made on rock analogs. However, our understanding of the processes that control frictional aging is limited and is dependent on experiments that explore only numerous mechanisms. We conduct slide‐hold‐slide experiments with a dual‐axis nanoindenter on single‐crystal surfaces of quartz and calcite. Our results show that frictional aging in diamond‐quartz contacts is independent of time and contact area, in stark contradiction to past experiments done on quartz‐quartz contacts in rocks. Diamond‐calcite contacts show modest frictional aging, but still well below previous reported values from calcite‐calcite contacts. These results suggest that frictional aging of like‐on‐like minerals may be of chemical origin, as suggested in recent studies with atomic force microscopy and molecular dynamics simulations.

Gravity-induced seismicity modulation on planetary bodies and their natural satellites

Ground-based monitoring of seismicity and modulation by external forces in the field of planetary seismology remains equivocal due to the lack of natural observations. Constrained by the natural observations (including Earthquakes, Moonquakes, and Marsquakes) and theoretical models, we present the variation in gravitational acceleration “g” of different solar system objects, combined with external harmonic forcings that are responsible for seismicity modulation on the planetary bodies and their natural satellites. From the global diversity in seismicity modulation, it has been observed that the plate-boundary regions on the Earth exhibit both short and long-period seismicity modulation. In contrast, the stable plate interior regions appear to be more sensitive to long-period seismicity modulation, however, lacking in short-period modulation. The deep Moonquakes are susceptible for both the lunar tidal period (13.6 days and 27 days) and long-period pole wobble modulation (206 days), whereas shallow emergent type moonquakes show a seismic periodicity at the lunation period (29.5 days). Further, the seasonal variation with an annual seismicity burst and seismic periodicity at polar wobble periods for high-frequency Marsquakes captured by InSight lander indicate a natural origin. Whereas diurnal and semi-diurnal periodicity along with Phobos’ tidal period, indicate possible artifacts due to different detection probabilities and non-seismic noise in the Martian environment. We argue that, in the context of rate-state-dependent fault friction, the gravity-induced resonance destabilization model appears to be better agreement with the contrast and relative diversity in seismicity modulation linked to the Earth, Moon, and Mars.

Dynamic-Projection-Integrated Particle-Filtering-Based Identification of Friction Characteristic Curve for Train Wheelset on Slipping Fault Condition

This paper proposes a dynamic-projection-integrated particle-filtering-based identification strategy for the friction characteristic curve of a train wheelset under the slipping fault condition. This strategy aims to achieve the identification of the fault friction characteristic curve (FFCC) in the early slipping fault stage. First, a multi-dimensional integrated particle-filtering (MDIPF)-based parameters correction method is proposed. The MDIPF constructs an error particle state transition model encompassing multi-dimensional parameters, which integrates inter-particle correlation to facilitate error fusion during the state transition process. Then, a dynamic projection domain (DPD)-based particle refinement method is proposed. The DPD constructed the contraction factors to dynamically fine-tune the particle projection domain. Finally, a multi-level evaluation-based identification method for the FFCC is proposed. And the dynamic-projection-integrated particle-filtering-based identification strategy is validated, which can actualize the rapid and accurate identification of the FFCC.