Shear localisation controls the dynamics of earthquakes

Fabian Barras,Nicolas Brantut

2023-08-01

Abstract:Earthquakes are produced by the propagation of rapid slip along tectonic faults. The propagation dynamics is governed by a balance between elastic stored energy in the surrounding rock, and dissipated energy at the propagating tip of the slipping patch. Energy dissipation is dictated by the mechanical behaviour of the fault, which is itself the result of feedbacks between thermo-hydro-mechanical processes acting at the mm to sub-mm scale. Here, we numerically simulate shear ruptures using a dual scale approach, allowing us to couple a sub-mm description of inner fault processes and km-scale elastodynamics, and show that the sudden localisation of shear strain within a shear zone leads to the emergence of classical cracks driven by a constant fracture energy. The fracture energy associated to strain localisation is substantially smaller than that predicted assuming uniform shearing. We show the existence of a unique scaling law between the localised shearing width and the rupture speed. Our results indicate that earthquakes are likely to be systematically associated to extreme strain localisation.

Geophysics

What problem does this paper attempt to address?

The problem that this paper attempts to solve is how the shear localization phenomenon in earthquake rupture dynamics controls the dynamic process of earthquakes. Specifically, the paper explores the phenomenon that during fault - slip, shear strain suddenly concentrates within a narrow area (i.e., within the shear zone), and how this phenomenon affects the propagation speed of earthquake rupture and the energy dissipation mechanism. The author studied the relationship between shear localization and earthquake rupture dynamics through numerical simulation, especially how shear localization leads to the formation of classical fractures, and these fractures are driven by constant fracture energy. In addition, the paper also reveals that there is a unique scaling law between the width of shear localization and the rupture speed, indicating that earthquakes are usually associated with extreme shear strain localization.

### Key Point Summary:

1. **Definition of Shear Localization**:
- Shear localization means that during fault - slip, shear strain concentrates within a very narrow area instead of being evenly distributed on the fault plane.

2. **Influence of Shear Localization**:
- Shear localization causes a sharp drop in shear stress within the fault plane, thus producing ruptures similar to classical fractures.
- These fractures are driven by constant fracture energy, and the fracture energy is significantly reduced compared to the situation assuming uniform shear.

3. **Relationship between Shear Localization and Rupture Speed**:
- The paper discovered a unique scaling law that describes the relationship between the width of shear localization and the rupture speed.
- This relationship indicates that as the rupture speed increases, the degree of shear localization also increases, causing the area of deformation, heating, and fluid pressure increase to become narrower.

4. **Numerical Simulation Method**:
- The author used a two - scale method, combining the description of internal fault processes at the sub - millimeter scale and the elastodynamic simulation at the kilometer scale.
- Through this method, the author was able to study in detail the influence of shear localization on rupture dynamics.

5. **Theoretical and Experimental Verification**:
- The theoretical predictions in the paper are consistent with the numerical simulation results, further confirming the important role of shear localization in earthquake rupture dynamics.
- The author also discussed how these results can be applied to explain complex phenomena in natural earthquakes and how to infer the dynamic process of earthquake rupture through the localization features in geological records.

### Formula Summary:

- **Relationship between Shear Stress and Strain Rate**:
\[
\tau \sim f(\dot{\gamma}, \vartheta)
\]
where \(\tau\) is shear stress, \(\dot{\gamma}\) is shear strain rate, and \(\vartheta\) is a diffusive quantity (such as temperature or pore pressure).

- **Diffusion Equation**:
\[
\dot{\vartheta} = \beta \tau \dot{\gamma} + \alpha \nabla^2 \vartheta
\]
where \(\beta\) is similar to the Taylor - Quinney coefficient, \(\alpha\) is the diffusivity, and \(\nabla^2\) is the Laplace operator.

- **Localization Condition**:
\[
\frac{h'_0}{g'_0}<0
\]
where \(g'_0 = \frac{\partial f}{\partial \dot{\gamma}}\bigg|_{t = t_0}\), \(h'_0=\frac{\partial f}{\partial \vartheta}\bigg|_{t = t_0}\).

- **Critical Wavelength**:
\[
\lambda > \frac{2\pi s}{\sqrt{-\frac{\alpha}{\beta f_0 g'_0 h'_0}}} \equiv \lambda_c
\]
where \(\lambda_c\) is the critical wavelength and \(\lambda\) is the perturbation wavelength.

- **Localization Width**:
\[
W_{\text{loc}} \sim \frac{\lambda_c}{2}
\]

- **Shear Band Width**