# Stability of pulse-like earthquake ruptures

## TL;DR

This paper investigates the stability of pulse-like earthquake ruptures driven by thermal pressurisation, revealing their inherent instability and implications for fault strength and rupture complexity through simulations and theoretical modeling.

## Contribution

It introduces a dynamic elastodynamic model showing that steady-state slip pulses driven by thermal pressurisation are unstable, providing insights into rupture evolution and fault strength.

## Key findings

- Steady-state pulse-like ruptures are unstable under thermal pressurisation.
- Positive stress perturbations lead to rupture growth and transition to other modes.
- Unstable pulses influence minimum stress conditions for sustainable ruptures.

## Abstract

Pulse-like ruptures arise spontaneously in many elastodynamic rupture simulations and seem to be the dominant rupture mode along crustal faults. Pulse-like ruptures propagating under steady-state conditions can be efficiently analysed theoretically, but it remains unclear how they can arise and how they evolve if perturbed. Using thermal pressurisation as a representative constitutive law, we conduct elastodynamic simulations of pulse-like ruptures and determine the spatio-temporal evolution of slip, slip rate and pulse width perturbations induced by infinitesimal perturbations in background stress. These simulations indicate that steady-state pulses driven by thermal pressurisation are unstable. If the initial stress perturbation is negative, ruptures stop; conversely, if the perturbation is positive, ruptures grow and transition to either self-similar pulses (at low background stress) or expanding cracks (at elevated background stress). Based on a dynamic dislocation model, we develop an elastodynamic equation of motion for slip pulses, and demonstrate that steady-state slip pulses are unstable if their accrued slip $b$ is a decreasing function of the uniform background stress $\tau_\mathrm{b}$. This condition is satisfied by slip pulses driven by thermal pressurisation. The equation of motion also predicts quantitatively the growth rate of perturbations, and provides a generic tool to analyse the propagation of slip pulses. The unstable character of steady-state slip pulses implies that this rupture mode is a key one determining the minimum stress conditions for sustainable ruptures along faults, i.e., their ``strength''. Furthermore, slip pulse instabilities can produce a remarkable complexity of rupture dynamics, even under uniform background stress conditions and material properties.