Fault Mechanics Simulator: Stress, Friction & Earthquake Recurrence

simulator intermediate ~10 min
Loading simulation...
T_r = 2000 yr — 50 km fault with Mw 7.1 potential

A 50 km fault with μ=0.6 and σn=100 MPa has a shear strength of 60 MPa. At 30 kPa/yr loading rate, the recurrence interval is ~2000 years, capable of producing a Mw 7.1 earthquake with ~1.5 m average slip.

Formula

τ_f = μ × σ_n — Coulomb failure criterion
T_recurrence = τ_f / (dτ/dt) — recurrence interval
Mw = (log10(L) + 4.86) / 1.32 — Wells & Coppersmith scaling

The Seismic Cycle

Earthquakes are not random events but the culmination of a slow, relentless stress-accumulation process driven by plate tectonics. As tectonic plates converge, diverge, or slide past each other at rates of centimeters per year, shear stress builds on locked fault segments. When the accumulated stress exceeds the fault's frictional strength, sudden slip occurs — an earthquake — releasing the stored elastic energy as seismic waves and resetting the stress clock.

Coulomb Friction

The Coulomb failure criterion τ_f = μ × σ_n provides the simplest model for fault strength. Laboratory experiments by Byerlee (1978) showed that most rocks follow μ = 0.6–0.85 regardless of rock type — 'Byerlee's law'. However, real faults may be weaker due to elevated pore fluid pressure, clay minerals in fault gouge, or dynamic weakening during rupture. The effective normal stress σ_n = σ_total - P_pore means that pore pressure plays a critical role in fault mechanics.

Stress, Slip, and Magnitude

The stress drop during an earthquake — typically 1–10 MPa — is a small fraction of the absolute stress level. Fault slip D scales with the fault dimensions: larger faults produce more slip and higher magnitudes. Empirical scaling relations (Wells & Coppersmith, 1994) relate fault length, width, slip, and magnitude, providing the basis for seismic hazard assessment. A 50 km fault typically produces 1–2 m of slip in a Mw 7 earthquake.

Beyond Simple Models

Real fault behavior is far more complex than the simple elastic rebound model. Rate-and-state friction laws describe how fault friction evolves with slip velocity and contact time, explaining phenomena like earthquake nucleation, aftershock sequences, and aseismic creep. Some fault segments slip quietly (slow-slip events), while others remain locked for centuries before rupturing catastrophically. Understanding this spectrum of fault behavior is the central challenge of modern earthquake science.

FAQ

What is the Coulomb failure criterion for faults?

The Coulomb criterion states that a fault slips when the shear stress τ exceeds the frictional strength μ×σ_n, where μ is the friction coefficient and σ_n is the effective normal stress (total stress minus pore pressure). Laboratory experiments show μ ≈ 0.6–0.85 for most rock types (Byerlee's law).

What controls earthquake recurrence?

Recurrence is governed by the ratio of fault strength to tectonic loading rate. Stronger faults with slow loading have long intervals between large events. The 'seismic cycle' model assumes stress builds linearly until reaching the strength, then drops suddenly during rupture (Reid's elastic rebound theory).

How does fault length relate to earthquake magnitude?

Fault length sets the upper bound on earthquake magnitude through empirical scaling relations like Mw ≈ (log10(L) + 4.86)/1.32 (Wells & Coppersmith, 1994). A 10 km fault produces at most ~Mw 6.2, a 100 km fault ~Mw 7.5, and a 1000 km fault ~Mw 9.0.

What is the elastic rebound theory?

Proposed by H.F. Reid after the 1906 San Francisco earthquake, elastic rebound theory states that tectonic motion gradually deforms rocks on either side of a locked fault, accumulating elastic strain energy. When stress exceeds the fault's frictional strength, sudden slip releases the stored energy as seismic waves and the rock 'rebounds' to its unstrained state.

Sources

Embed

<iframe src="https://homo-deus.com/lab/seismology/fault-mechanics/embed" width="100%" height="400" frameborder="0"></iframe>
View source on GitHub