Gutenberg-Richter Recurrence Calculator

The Gutenberg-Richter law states that the cumulative number of earthquakes of magnitude M or larger in a region follows log10(N) = a - bM, a power law in the underlying energy. The a-value fixes the overall productivity and depends on the area surveyed, the length of the catalogue and the tectonic setting, so it is meaningless without those specifications. The b-value is the slope and describes the proportion of small to large events: at b equal to 1.0, each unit increase in magnitude reduces the event count tenfold. Observed b-values cluster near 1.0 globally, with roughly 0.6 to 0.8 in stressed asperity-dominated fault segments and subduction interfaces, and 1.2 to 2.5 or higher in volcanic and geothermal swarms where many small events accompany fluid movement. The annual rate N converts to a mean return period of 1/N years, and if events are treated as a Poisson process the probability of at least one exceedance in t years is 1 - exp(-N t). That Poisson assumption is the main modelling caveat, since it implies memoryless timing and is violated by aftershock clustering and by characteristic earthquakes on individual faults, which recur more regularly than a Poisson process predicts. Catalogue issues matter just as much: the fit is only valid above the magnitude of completeness Mc, the smallest magnitude reliably detected by the network, and extrapolating the line beyond the largest magnitude a fault system can physically host will overstate the hazard, so seismic hazard models usually impose a maximum magnitude taper.

The frequency-magnitude relation was recognized independently on two continents in the same era. Mishio Ishimoto and Kumizi Iida of the Tokyo Earthquake Research Institute published a power-law distribution of maximum trace amplitudes in 1939, an equivalent formulation in terms of amplitude rather than magnitude. Beno Gutenberg and Charles Richter, working at the Caltech Seismological Laboratory with the magnitude scale Richter had defined in 1935, presented the logarithmic frequency-magnitude form in their 1944 Bulletin of the Seismological Society of America paper on the frequency of earthquakes in California, and expanded it in their 1949 book Seismicity of the Earth. Estimating b reliably from finite catalogues became the next problem. Keiiti Aki derived the maximum likelihood estimator in a 1965 paper, showing that b follows from the mean magnitude above the completeness threshold, and Tokuji Utsu refined the treatment for binned magnitudes, giving the Aki-Utsu formulation used in practice. From the late 1960s the relation was embedded in probabilistic seismic hazard analysis, the framework introduced by C. Allin Cornell in 1968, which combines Gutenberg-Richter recurrence with ground motion attenuation to produce the hazard maps behind modern building codes. Later work by researchers including Steven Wiemer and Max Wyss mapped b-value variations in space and time as a proxy for differential stress, and observed b-value drops preceding some large ruptures.