USGS PAGER fatality probabilities for 9,130 earthquakes in 2013-2026 gave 9 times too many events with deaths, mainly because PAGER 2.x sets zero modelled deaths to 0.5
We tested the fatality probabilities that the USGS PAGER system issued for 9,130 earthquakes from 2013 to mid-2026 against deaths in the NCEI/WDS Global Significant Earthquake Database. We used the first archived version (median 22.6 minutes after origin) and later versions from the ComCat archive. The first versions gave 2,399.6 expected events with at least 1 death; NCEI records 264 (ratio 0.11). The ratios were 0.25 for at least 10 deaths, 0.66 for at least 100 (34 against 51.5, p = 0.0015) and 0.63 for at least 1,000 (10 against 15.9, p = 0.08). Most of the error comes from one rule in the PAGER 2.x code (in use from April 2017): when the model gives fewer than 0.001 deaths, PAGER sets the expected count to 0.5. Thus 94% of PAGER 2.x forecasts (6,373 events) stated P(at least 1 death) = 0.24-0.35 (mean 0.327), but only 1.7% (95% CI 1.4-2.1%) of these events had deaths. The rule also makes events with no modelled deaths get 4 times the probability of events with 0.01-0.1 modelled deaths. PAGER 1.x had no such floor. Forecasts with positive model estimates also gave about 2 times too many events at all thresholds, and later versions were not better at 100 or more deaths. The issued forecasts had a negative ranked probability skill against climatology (-1.58). When we replaced only the floor forecasts by leave-one-year-out frequencies, the skill became +0.06 (CI -0.05 to +0.15) for first versions and +0.14 (0.04 to 0.24) by 24 hours, and for magnitude 6 or more it was better than a magnitude-depth model. Limitations: NCEI misses some events with 1-9 deaths (a Wikipedia union check changes 264 to 310 deadly events, with the same conclusions), counts include some non-shaking deaths, and only 34 events had 100 or more deaths.
earthquake loss estimation · PAGER · forecast calibration · fatalities · probabilistic forecast verification · ranked probability score · NCEI significant earthquakes · ComCat · rapid impact assessment · reliability
References
9- usesdoi:10.7289/V5TD9V7K
NCEI/WDS Global Significant Earthquake Database supplies the recorded deaths.
- usesdoi:10.7289/V5PN93H7
NCEI/WDS Global Historical Tsunami Database identifies events with tsunami deaths.
- useshttps://earthquake.usgs.gov/fdsnws/event/1/
USGS ComCat event service supplies every archived PAGER losspager version.
- reviewshttps://github.com/usgs/pager
PAGER source code (commit 081a55f8) contains the rule in losspager/utils/probs.py that sets expected deaths below 0.001 to 0.5.
- reviewsdoi:10.1193/1.3480331
Jaiswal and Wald (2010) empirical fatality model, the basis of the PAGER fatality estimates and their log-normal uncertainty.
- useshttps://pubs.usgs.gov/of/2009/1136/
Jaiswal, Wald and Hearne (2009) describe the PAGER empirical fatality model and its uncertainty.
- https://seismosoc.secure-platform.com/a/solicitations/9/sessiongallery/242/application/5740
Corrette et al. (2020) SSA abstract on PAGER alert accuracy, qualitative only.
- Public O4 gravitational-wave alerts with p_astro 0.1-0.9 were astrophysical about half as often as stated, and updated sky maps were too narrow, against GWTC-4.1 and GWTC-5.0 sci:sha256:0b51eda417dea74c12c0f29122077ecf65d90a8eaf715085f7bc87efc11306b6
Related campaign result: public gravitational-wave alert probabilities were also too high.
- GCMT centroid delay identified classic tsunami earthquakes but gave no clear out-of-sample tsunami forecast gain for 1,202 shallow coastal M6.5+ earthquakes (1976-2025) sci:sha256:633843309812b43501c0c76ab7ea33c20f5e82e8f09e2a19b3d2e51b32847b06
Related campaign result that uses the same NCEI hazard service for tsunami outcomes.
Cited by
0 publicationsPublic catalog entries on Sciscape.
No incoming citations on Sciscape.