Paleomagnitude Estimation Simulator: Sizing Ancient Earthquakes

simulator intermediate ~10 min
Loading simulation...
M ≈ 7.1 — from 3 m displacement and 60 km rupture

A paleoearthquake with 3 m surface displacement and 60 km rupture length yields a magnitude estimate of approximately 7.1 using Wells & Coppersmith (1994) scaling relations — a major earthquake capable of producing severe damage over a wide area.

Formula

M = 6.69 + 0.74 × log₁₀(D) (Wells & Coppersmith, displacement)
M = 5.08 + 1.16 × log₁₀(L) (Wells & Coppersmith, length)
M₀ = μ × L × W × D (seismic moment, N·m)

Sizing the Invisible

When an earthquake struck before seismometers existed, the only magnitude evidence lies in the rocks. Surface rupture length, fault displacement, and deformation patterns preserved in the geological record can be measured centuries or millennia later. Empirical scaling relations, calibrated against hundreds of instrumentally recorded earthquakes, translate these measurements into magnitude estimates — connecting the geological record to the quantitative framework of modern seismology.

Scaling Relations

The Wells & Coppersmith (1994) study analyzed 244 earthquakes with known magnitudes and measured surface parameters, establishing the field's standard scaling relations. The key insight is that magnitude scales logarithmically with displacement and rupture dimensions: a tenfold increase in displacement corresponds to roughly 0.7 magnitude units. These relations have been refined for different tectonic environments but the original formulations remain widely used for their simplicity and robustness.

Multiple Estimators

Using both displacement and rupture length provides independent magnitude estimates that should agree within ~0.3 units if the data are reliable. Significant disagreement may indicate incomplete rupture mapping, non-characteristic behavior, or distributed deformation not captured at a single trench site. Seismic moment, computed from rigidity, area, and displacement, provides a physically grounded estimate independent of the empirical scaling.

Uncertainty and Consequences

Paleomagnitude uncertainty of ±0.3 units sounds small but corresponds to a factor of ~2 in seismic energy release. For hazard assessment, this uncertainty propagates into ground motion predictions and ultimately into structural design requirements. This simulation lets you explore how the measurable geological evidence constrains the magnitude of ancient earthquakes and how different fault types shift the scaling relationships.

FAQ

How do scientists estimate magnitudes of prehistoric earthquakes?

Paleomagnitudes are estimated using empirical scaling relations that link measured geological evidence (surface displacement, rupture length, rupture area) to magnitude. The most widely used relations are from Wells & Coppersmith (1994), which compiled data from 244 historical earthquakes to establish log-linear relationships between these parameters and moment magnitude.

What are the Wells & Coppersmith scaling relations?

The key relations are: M = 6.69 + 0.74×log(D) for displacement, M = 5.08 + 1.16×log(L) for rupture length, and M = 4.07 + 0.98×log(A) for rupture area. These have standard deviations of 0.3–0.4 magnitude units, meaning paleomagnitude estimates carry inherent uncertainty of roughly ±0.3 M.

Does fault type affect magnitude estimates?

Yes. Strike-slip, normal, and reverse faults have slightly different scaling relations because their rupture geometries differ. Reverse (thrust) faults tend to produce larger displacements for a given magnitude, while strike-slip faults tend to produce longer ruptures. The differences are typically 0.1–0.3 magnitude units.

What are the uncertainties in paleomagnitude estimation?

Uncertainties arise from: incomplete preservation of displacement (erosion), unmapped portions of the rupture, using maximum vs. average displacement, and the scatter in the empirical scaling relations themselves. Total uncertainty is typically ±0.3–0.5 magnitude units, which translates to a factor of 2–3 in seismic energy.

Sources

Embed

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