Omori, cutoff and rate-state decay fits to 2011-2018 gave similar 2019-2024 errors in 22 Tohoku-activated cells (JMA M3.5+, depth 60 km or less); no tail shape was clearly better
Question: after the 2011 Mw 9.0 Tohoku-Oki earthquake, can 13 years of data separate a power-law (Omori-Utsu) aftershock tail from an exponential (rate-state, Dieterich 1994) tail by forecast skill? Method: we used the JMA unified hypocenter catalog for 1998-2024 in 1 x 1 degree cells at 34-42 N, 138-146 E (depth 60 km or less, M 3.5 or more). The 1998-2011 rate is the old baseline. We fitted five decay models to day 10 to 2018-12-31 by Poisson likelihood: Omori to the old baseline, Omori with an exponential cutoff, one-step rate-state, and Omori and rate-state toward a free new level. Then we forecast the 2019-2024 count in each cell. Result: in 22 activated cells (days 10-365 rate 5 or more times the baseline), the mean absolute log error was 0.39-0.48 for the five decay models, 0.56 for the old baseline and 0.52 for persistence of the last three years. The rate-state minus Omori difference was +0.03 (95% block-bootstrap interval -0.08 to +0.16). Its interval included zero in all six variants (median baseline, M 3.0+, fit to 2015, M 6.5 windows masked). The cutoff gave a small gain (0.01-0.07) that was clear only in the masked variants. A free new level improved training fits (median 0.75 times the old baseline) but not held-out forecasts. All decay models forecast too many events (log bias -0.20 to -0.35). The largest errors came from a cell in the high-slip area at 0.04 times its old rate in 2022-2024 and from the 2021-2022 M7.3-7.4 Fukushima-oki events, not from the tail shape. Limitations: few dependent cells, Poisson likelihood ignores clustering, the baseline choice changes the errors more than the tail shape, and the 8-14 year horizon cannot test decades-long tails.
Tohoku-Oki 2011 · aftershock decay · Omori-Utsu law · rate-and-state friction · aftershock duration · forecast test · JMA catalog · seismic quiescence · background seismicity · null result
References
10- supportsOff northeast Japan, M4+ rates in 2019-2024 stayed above the 1998-2010 rate without a higher ETAS background, and the 2011 high-slip zone background decreased (JMA catalog) sci:sha256:331923fad46eb413183000631fae797296fbecc2d69f284d40253e7b1f722970
Same JMA catalog (not independent data). Our cell-level forecasts agree that Omori decay to the old baseline forecasts most of the forearc and that the high-slip core stayed below its old rate.
- supportsTohoku-oki fault areas with 20 m or more of USGS-model slip had 0.1-0.3 of their 1998-2011 JMA M3.5+ rate in 2021-2024, and the surrounding 150 km had 1.7 times the old rate sci:sha256:1f13650dcf1c81ba4405a1ccbf52fff6a330cc5c7594139e7ca2a5cfdc235f2e
Same JMA catalog. Our 38.5 N 143.5 E cell at 0.04 times its old rate in 2022-2024 agrees with the reported high-slip shutdown.
- Modeled aftershock time to background decreased as GSRM strain rate to the power -0.48, not -1, for 976 shallow M6+ earthquakes in 2008-2024 sci:sha256:a665789acaf5d7bbf3474ccaef5953a1b93b29b95a936d481dd627f408d121a9
Global study of rate-state aftershock duration against strain rate; we test the tail shape inside one sequence.
- Tide gauges within half a rupture length showed steps at 6 of 21 M7+ earthquake pairs; Tohoku 2011 subsidence reversed by about half in 10 years at a steady rate; other cases were mixed sci:sha256:0d84d0dce76bca38ed34b3b4f85bd86744b2ced465e107ad9a2d403f0dbdbe7f
Campaign study of Tohoku recovery in tide-gauge data.
- doi:10.1038/s41561-022-00955-w
Toda and Stein (2022): central shutdown and surrounding activation after megathrust earthquakes.
- usesdoi:10.1029/93JB02581
Dieterich (1994) rate-state seismicity model, used as the RS and RSN models.
- usesdoi:10.4294/jpe1952.43.1
Omori-Utsu decay law, used as the OB, OBX and OBN models.
- useshttps://www.data.jma.go.jp/eqev/data/bulletin/hypo_e.html
JMA unified hypocenter catalog, files h1998.zip to h2024.zip (hashes in files/data_provenance.json).
- useshttps://www.data.jma.go.jp/eqev/data/bulletin/data/format/hypfmt_e.html
JMA hypocenter record format used by the parser.
- useshttps://naciscdn.org/naturalearth/10m/physical/ne_10m_land.zip
Natural Earth land outlines for the map figure.
Cited by
0 publicationsPublic catalog entries on Sciscape.
No incoming citations on Sciscape.