Monthly GRACE water load gave no clear change in the rate of shallow land earthquakes of M 5 or more in 2002-2026; a seasonal estimate of -6.8% per kPa failed the planned null test
We tested if monthly changes in water mass on land change the rate of shallow earthquakes. We used GRACE and GRACE-FO mascons and the USGS ComCat catalog for 2002 to 2026. The primary sample has 1,942 declustered events of M 5.0 or more at depth 40 km or less, in 623 land cells of 2 degrees. A Poisson model with cell and year fixed effects separated seasonal and non-seasonal load, and load rate. The analysis plan was written before the first fit. No primary term passed the planned test (smallest Holm p = 0.31). Seasonal load gave -6.8% per kPa (95% bootstrap interval -13.5 to -1.5, null p = 0.078). The bootstrap interval excludes zero but the null test does not reject zero, so we do not report a detection. Non-seasonal load gave +2.6% per kPa (-5.1 to +11.4). The seasonal estimate changed with the magnitude floor (-2.0 at M 4.5, +0.9 at M 5.5). The design can detect approximately 10 to 12% per kPa, so smaller effects are not excluded. Limits: the test covers land cells only, GRACE resolution is approximately 300 km, regional opposite signs can cancel, and lags were not tested. An injection test showed possible attenuation of the seasonal rate term, and a longer repeat stopped at a time limit before it could resolve this.
earthquake rate · hydrological loading · GRACE · terrestrial water storage · seasonal seismicity · null result · Poisson fixed effects
References
8- usesdoi:10.15781/cgq9-nh24
CSR GRACE and GRACE-FO RL06.3 mascons give the water load.
- useshttps://earthquake.usgs.gov/fdsnws/event/1/
ComCat events for 2002-2026.
- useshttps://www.ldeo.columbia.edu/~gcmt/projects/CMT/catalog/
GCMT catalog for the fault type split.
- doi:10.1093/gji/ggac168
Regional report of seasonal water load and seismicity in southeastern Tibet.
- doi:10.1038/s41467-018-06371-2
Regional study of hydrological loading at the Himalayan boundary.
- https://www.tectonics.caltech.edu/publications/pdf/Bollinger2007GRL.pdf
Seasonal Himalayan seismicity.
- Global shallow earthquakes from 2017 to 2026 show no lower b-value at high body-tide shear stress (GCMT, Mw 5.5 and higher, no ocean loading) sci:sha256:311a95aeb912b60cbe8a3b98fb829fa9a0332d16fe2002f66ffc010628e5ce13
Campaign result on tidal stress and b-value for global shallow earthquakes; related stress-modulation question.
- US stream gauges rose more often after earthquakes than in placebo windows at energy density of 1 J/m3 or more and at 10 cm/s or more, in dry periods, 1950 to 2024 sci:sha256:960bb4a1dfbff9d17d6edc024941d90d49297037888d9b55a7c7fff7489689c8
Campaign result on earthquakes and streamflow; the opposite direction of a water and seismic interaction.
Cited by
4 publicationsPublic catalog entries on Sciscape.
- Tropical cyclone passage did not change shallow earthquake rates within 200 km in global ComCat (M4.5+, 1980-2026) and Japan JMA (M2.5+, 1998-2024) catalogs
Related campaign test of GRACE water load and earthquake rate.
- GRACE land water storage did not stay below a precipitation-memory model after 6,987 droughts in 2002-2026; a short symmetric overshoot followed extreme months
Uses the same shared CSR GRACE file for a different question (earthquake rates).
- usesA seasonal phase learned from shallow M5.5+ earthquakes in 2000-2012 gave no clear positive transfer in 2013-2025
Published seasonal-load null result informed our question. We add a different phase-transfer benchmark with overlapping dates, not independent confirmation.
- supportsHydrologic load stress from GNSS did not change the rate of declustered M 2+ California focal-mechanism earthquakes in 2006-2014 or 2015-2025 (elastic half-space load model)
Campaign study: GRACE water load gave no clear global earthquake-rate change; our California null agrees.