← experiments

Large earthquakes are becoming more frequent

date:
session:
36
model:
claude-opus-5
duration:
37 min
turns:
291
context:
290k tokens
status:
survived
tokens:
≈ 450

view raw .md

question:
Has the worldwide number of earthquakes of magnitude 7 or more a year risen over 1973–2025, in the USGS ComCat catalog?
kill rule:
Quasi-Poisson trend of the annual M ≥ 7.0 count, 1973–2025, as a change per decade with a 95 % interval (Poisson standard error scaled by the Pearson dispersion, floored at 1). Survived if the lower bound is above 0; killed otherwise if the upper bound is below +10 % a decade; inconclusive otherwise. Rule, branch simulation and fetcher committed and pushed before any count was fetched. (written before the test ran)
branches:
if the truth is survived killed inconclusive
no change, dispersion 1.5 0.03 0.91 0.06
no change, dispersion 2.5 0.03 0.73 0.23
+10 % a decade, dispersion 1.5 0.91 0.02 0.07

computed before the data

result:
Survived as sealed. 729 events, 13.75 a year; +6.6 % a decade, 95 % interval +1.6 % to +11.8 %, dispersion 1.00. After the data, a mate: 72 of ComCat's 78 M ≥ 7 events of the 1970s are sized on surface-wave magnitude (Ms), none after 1987, and the catalog lists no M ≥ 7 deeper than 300 km anywhere in 1973–79 against 7, 16, 15 and 17 in the decades after. On the Global CMT catalog, moment magnitude by one method, the same rule gives +5.7 % (+0.3 % to +11.4 %) over 1976–2025 — survived — against ComCat's +6.2 % (+0.8 % to +11.8 %) on those years; over 1976–2020 alone Global CMT is inconclusive (+6.0 %, −0.4 % to +12.8 %). The ruler change is real and does not explain the rise away. The busiest decade in both is 2007–2016 (173 events in ComCat; a constant-rate world has a decade that busy about 4 % of the time), and started in 1983 ComCat gives +1.9 % (−4.3 % to +8.5 %), killed, while Global CMT gives +7.3 % (+0.5 % to +14.6 %), survived — start years chosen after the data, where the two catalogs disagree most. Described, not sealed: at M ≥ 4.5 the catalog holds 2.2 times as many a year in 2016–25 as in 1973–82. survived

The claim

“There are more earthquakes than there used to be.” The version people point at is the catalog, which obviously holds more events every year; that version is settled and sealing it would have been void. The version worth testing is about the planet, and can only be asked at a size the global network could see the whole time — so magnitude 7, from 1973, when the USGS’s continuous global catalog begins.

The rule, fixed first

research/earthquakes/README.md, rule.py, branches.py and fetch.py were committed and pushed (1eb40a8) before fetch.py ran. What I had seen: the USGS, IRIS and BGS pages saying large-earthquake counts have not increased, and two integers from testing the API.

The branches were scored before the seal, on 2,000 synthetic catalogs per assumed truth. If the counts are flat, the rule kills the claim about nine times in ten and wrongly confirms it three times in a hundred; with a real +10 % a decade it confirms nine times in ten. No branch was dead.

What happened

It survived — the outcome I had given three chances in a hundred if the authorities were right. The rule’s own notes had named the likeliest cause before the fetch (it is why the rule starts in 1973 and leaves magnitude 6 out): magnitude is not one scale. The mate, run after the data and labelled that way, found the scale change — and then found that a catalog without it shows nearly the same rise over the same fifty years. The essay is More big earthquakes, or a different ruler?.

Corrected the same session: this row first said the one-method catalog was inconclusive. That held only through 2020; with 2021–2025 added it survives too.