Earthquake magnitude-frequency relation
Gutenberg-Richter Law
In many regions and time windows, the cumulative number of earthquakes decreases approximately tenfold for each unit increase in magnitude.
log10 N(M or greater) = a - bM
N is the number or rate of earthquakes at least magnitude M. Parameter a controls overall seismic productivity; b controls the relative proportion of large to small events and is often near 1, but must be estimated.
Illustrative parameters are a = 6 and b = 1, giving N = 10^(6-M). This is not a forecast for any real region; actual a, b, completeness, area, and time window must be estimated from a catalog.
(events/year)
- CHANGE
- Minimum earthquake magnitude
- WATCH
- event frequency
- MEANING
- Illustrative parameters are a = 6 and b = 1, giving N = 10^(6-M). This is not a forecast for any real region; actual a, b, completeness, area, and time window must be estimated from a catalog.
Large earthquakes are rare on a logarithmic ladder.
Equal magnitude steps correspond to multiplicative frequency changes. The straight line on log-count axes becomes a steep cascade when translated back into event counts.
What it actually says
The Gutenberg-Richter relation summarizes the size distribution of earthquakes rather than their timing. On a plot of cumulative log10 event count against magnitude, a catalog often forms an approximately straight line above its completeness threshold.
A b-value near 1 means about one-tenth as many earthquakes at magnitude M+1 or greater as at M or greater. Because magnitude itself is logarithmic and seismic moment grows faster than event count falls, rare large earthquakes can dominate total released moment.
"A useful law compresses a pattern. It does not erase the conditions that make the pattern true."
How the idea developed
The modern form emerged through observation, argument, and later refinement. The timeline separates the first insight from the version now used in textbooks and practice.[1]
Richter and Gutenberg define the local magnitude scale for southern California earthquakes.
Gutenberg and Richter publish a systematic magnitude-frequency relation for California.
Global catalogs and improved magnitude scales extend statistical seismology.
Magnitude-frequency models inform seismic hazard, aftershock analysis, induced seismicity, and catalog quality control.
How the pattern works
The relation becomes useful only when its mechanism, measurement process, and operating range are visible.
Fault systems produce many small ruptures and progressively fewer large ruptures across a range of scales.
N(M or greater) stabilizes noisy tail counts and leads to the conventional log-linear expression.
Parameter a varies with area, observation time, tectonic rate, and catalog definition.
Parameter b changes the large-to-small event ratio and can vary across regions, sequences, stress states, and methods.
N is the number or rate of earthquakes at least magnitude M. Parameter a controls overall seismic productivity; b controls the relative proportion of large to small events and is often near 1, but must be estimated.
Where it earns its keep
Applications are strongest when the law changes a decision, measurement, model, or experiment rather than merely providing an analogy.
Estimate recurrence rates by magnitude
ApplicationCatalog fits contribute to probabilistic forecasts of how often damaging events may occur.
Extrapolation beyond observed magnitudes needs fault physics and maximum-magnitude constraints.
Estimate completeness threshold
ApplicationThe roll-off of small-event counts can reveal where detection becomes incomplete.
Fitting below completeness biases b downward and creates false structure.
Track changing event populations
ApplicationRate and b-value estimates can summarize sequences associated with injection or extraction.
Short windows produce unstable estimates and do not identify causation by themselves.
Where it stops working
Catalogs miss small events, mix magnitude scales, change instrumentation, and contain aftershock clustering. The fitted range, declustering method, spatial boundary, time window, and magnitude uncertainty can materially alter a and b.
Individual faults or the largest events may deviate from a simple unbounded law. Characteristic-earthquake models, tapered distributions, finite fault dimensions, and maximum-magnitude constraints address behavior the straight line cannot.
"The law predicts when the next M7 will occur"
Better: It models frequency, not a deterministic event clock."b always equals exactly 1"
Better: It is estimated and varies with catalog and physical setting."A straight line proves one universal mechanism"
Better: Different rupture processes and mixtures can produce similar aggregate slopes."Tiny earthquakes release most seismic energy"
Better: Rare large events can dominate moment release despite lower counts.Sources and further reading
Original publications and serious secondary scholarship are prioritized over summaries.
- USGS - Calculating California Seismicity RatesOfficial statement and use of log N = a - bM in seismicity-rate estimation.https://www.usgs.gov/publications/calculating-california-seismicity-rates
- USGS National Seismic Hazard Model - Magnitude Frequency DistributionsUSGS implementation reference for Gutenberg-Richter rates.https://ghsc.code-pages.usgs.gov/nshmp/nshmp-lib/gov/usgs/earthquake/nshmp/mfd/Mfds.html
- USGS - Earthquake Magnitude, Energy Release, and Shaking IntensityOfficial explanation of magnitude scales and logarithmic earthquake size.https://www.usgs.gov/programs/earthquake-hazards/earthquake-magnitude-energy-release-and-shaking-intensity
- Parsons et al. - Characteristic Magnitude-Frequency Distributions on FaultsEvidence and cautions about departures from regional Gutenberg-Richter behavior.https://www.usgs.gov/publications/characteristic-earthquake-magnitude-frequency-distributions-faults-calculated