skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


This content will become publicly available on June 1, 2026

Title: Nonlinear stability analysis of slip in a single degree of freedom elastic system with frictional evolution laws spanning ageing to slip
We present a nonlinear stability analysis of quasi-static slip in a spring-block model. The sliding interface is governed by rate-and-state friction, with an intermediate state evolution law spanning ageing and slip laws. We examine the robustness of prior results to changes in the evolution law, including the unconditional stability of the ageing law for spring stiffnesses above a critical value. We investigate two scenarios: a spring-block model with stationary and non-stationary point loading rate. In the former scenario, deviations from the ageing law lead to only conditional stability for spring stiffnesses above a critical value: finite perturbations can trigger instability, consistent with prior results for the slip law. For a given supercritical stiffness, the perturbation size required to induce instability grows as the state evolution law approaches the ageing law. By contrast, for a stationary point loading rate, there exists a maximum critical stiffness above which instability can never develop, for any perturbation size. This critical stiffness vanishes as the slip law is approached. For a limited extension of the results found for a spring-slider to continuum faults, we derive relations for an effective spring stiffness as a function of elastic moduli and a characteristic fault dimension or perturbation wavelength.  more » « less
Award ID(s):
1834696
PAR ID:
10640953
Author(s) / Creator(s):
;
Publisher / Repository:
Royal Society Publishing
Date Published:
Journal Name:
Proceedings of the Royal Society A: Mathematical, Physical and Engineering Sciences
Volume:
481
Issue:
2316
ISSN:
1364-5021
Page Range / eLocation ID:
20240917
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract Geophysical and geological studies provide evidence for cyclic changes in fault‐zone pore fluid pressure that synchronize with or at least modulate slip events. A hypothesized explanation is fault valving arising from temporal changes in fault zone permeability. In our study, we investigate how the coupled dynamics of rate and state friction, along‐fault fluid flow, and permeability evolution can produce slow slip events. Permeability decreases with time, and increases with slip. Linear stability analysis shows that steady slip with constant fluid flow along the fault zone is unstable to perturbations, even for velocity‐strengthening friction with no state evolution, if the background flow is sufficiently high. We refer to this instability as the “fault valve instability.” The propagation speed of the fluid pressure and slip pulse, which scales with permeability enhancement, can be much higher than expected from linear pressure diffusion. Two‐dimensional simulations with spatially uniform properties show that the fault valve instability develops into slow slip events, in the form of aseismic slip pulses that propagate in the direction of fluid flow. We also perform earthquake sequence simulations on a megathrust fault, taking into account depth‐dependent frictional and hydrological properties. The simulations produce quasi‐periodic slow slip events from the fault valve instability below the seismogenic zone, in both velocity‐weakening and velocity‐strengthening regions, for a wide range of effective normal stresses. A separation of slow slip events from the seismogenic zone, which is observed in some subduction zones, is reproduced when assuming a fluid sink around the mantle wedge corner. 
    more » « less
  2. null (Ed.)
    Abstract Observations of glacier slip over till beds, across a range of spatial and temporal scales, show abundant seismicity ranging from Mw∼−2 microearthquakes and tremor (submeter asperities and millisecond duration) to Mw∼7 slow-slip events (∼50  km rupture lengths and ∼30  min durations). A complete understanding of the mechanisms capable of producing seismic signals in these environments represents a strong constraint on bed conditions. In particular, there is a lack of experimental confirmation of velocity-weakening behavior of ice slipping on till, where friction decreases with increasing velocity—a necessity for nucleating seismic slip. To measure the frictional strength and stability of ice sliding against till, we performed a series of double-direct-shear experiments at controlled temperatures slightly above and below the ice melting point. Our results confirm velocity-strengthening ice–till slip at melting temperatures, as has been found in the few previous studies. We provide best-fit rate-and-state friction parameters and their standard deviations from averaging 13 experiments at equivalent conditions. We find evidence of similar velocity-strengthening behavior with 50% by volume debris-laden ice slid against till under the same conditions. In contrast, velocity-weakening and linear time-dependent healing of ice–till slip is present at temperatures slightly below the melting point, providing an experimentally supported mechanism for subglacial seismicity on soft-beds. The stability parameter (a−b) decreases with slip velocity, and evolution occurs over large (mm scale) displacements, suggesting that shear heating and melt buildup is responsible for the weakening. These measurements provide insight into subglacial stiffness in which seismicity of this type might be expected. We discuss glaciological circumstances pointing to potential field targets in which to test this frozen seismic asperity hypothesis. 
    more » « less
  3. We perform numerical experiments of damped quasi-dynamic fault slip that include a rate-and-state behavior at steady state to simulate earthquakes and a plastic rheology to model permanent strain. The model shear zone has a finite width which represents a natural fault zone. Here we reproduce fast and slow events that follow theoretical and observational scaling relationships for earthquakes and slow slip events (SSEs). We show that the transition between fast and slow slip occurs when the friction drop in the shear zone is equal to a critical value, Δμc. With lower friction drops, SSEs use nearly all of mechanical work to accumulate inelastic strain, while with higher friction drops fast slips use some of the mechanical work to slip frictionally. Our new formulation replaces the state evolution of rate and state by the stress evolution concurrent with accumulation of permanent damage in and around a fault zone. 
    more » « less
  4. We combine a slip-spring model with an ‘entangled kink dynamics’ (EKD) model for strong uniaxial extensional flows (with Rouse Weissenberg number W i R ≫ 1 ) of long ( M w > 1   Mkg / mol for polystyrene) entangled polymers in solutions and melts. The slip-spring model captures the dynamics up to the formation of a ‘kinked’ or folded state, while the kink dynamics simulation tracks the dynamics from that point forward to complete extension. We show that a single-chain slip-spring model using affine motion of the slip-spring anchor points produces unrealistically high tension near the center of the chain once the Hencky strain exceeds around unity or so, exceeding the maximum tension that a chain entangled with a second chain is able to support. This unrealistic tension is alleviated by pairing the slip links on one chain with those on a second chain, and allowing some of the large tension on one of the two to be transferred to the second chain, producing non-affine motion of each. This explicit pairing of entanglements mimics the entanglement pairing also used in the EKD model, and allows the slip spring simulations to be carried out to strains high enough for the EKD model to become valid. We show that results nearly equivalent to those from paired chains are obtained in a single-chain slip-spring simulation by simply specifying that the tension in a slip spring cannot exceed the theoretical maximum value of ζ ′ ϵ ˙ L 2 / 8 where ζ ′ , ϵ ˙ and L are the friction per unit length, strain rate and contour length of the chain, respectively. The effects of constraint release (CR) and regeneration of entanglements is also studied and found to have little effect on the chain statistics up to the formation of the kinked state. The resulting hybrid model provides a fast, simple, simulation method to study the response of high molecular weight ( M w > 1   Mkg / mol ) polymers in fast flows ( W i R ≫ 1 ), where conventional simulation techniques are less applicable due to computational cost. 
    more » « less
  5. Nearly all frictional interfaces strengthen as the logarithm of time when sliding at ultra-low speeds. Observations of also logarithmic-in-time growth of interfacial contact area under such conditions have led to constitutive models that assume that this frictional strengthening results from purely time-dependent, and slip-insensitive, contact-area growth. The main laboratory support for such strengthening has traditionally been derived from increases in friction during “load-point hold” experiments, wherein a sliding interface is allowed to gradually self-relax down to subnanometric slip rates. In contrast, following step decreases in the shear loading rate, friction is widely reported to increase over a characteristic slip scale, independent of the magnitude of the slip-rate decrease—a signature of slip-dependent strengthening. To investigate this apparent contradiction, we subjected granite samples to a series of step decreases in shear rate of up to 3.5 orders of magnitude and load-point holds of up to 10,000 s, such that both protocols accessed the phenomenological regime traditionally inferred to demonstrate time-dependent frictional strengthening. When modeling the resultant data, which probe interfacial slip rates ranging from 3 . μ m · s − 1 . to less than 10 − 5 μ m · s − 1 , we found that constitutive models where low slip-rate friction evolution mimics log-time contact-area growth require parameters that differ by orders of magnitude across the different experiments. In contrast, an alternative constitutive model, in which friction evolves only with interfacial slip, fits most of the data well with nearly identical parameters. This leads to the surprising conclusion that frictional strengthening is dominantly slip-dependent, even at subnanometric slip rates. 
    more » « less