Automated b-value traffic lights classified 9 of 65 Japanese M5+ foreshock sequences before the larger event, and over 80% of red alarms were false (JMA, 2000-2024)
Question: Can the b-value foreshock traffic light (FTLS; Gulia and Wiemer 2019) tell in real time that a large earthquake is a foreshock? Approach: We applied the rule automatically and causally to all 1,459 shallow (depth 50 km or less) M5.0+ first events in the JMA unified catalog, 2000-2024 (240 of M6.0+). A first event is the largest event near it in the previous year. The post-event b-value used only events before a decision time of 0.25, 1, 3 or 10 days. The outcome was a larger event within one rupture length in the next 30 days. We varied 72 method settings and used three estimators (FTLS-style Aki-Utsu, common-Mc Aki-Utsu, b-positive). Findings: 65 first events had a larger event within 30 days, and half of these came within 8 hours. Only 8% of first events had enough events for a b-value at 1 day. The FTLS gave a valid class before the larger event in 9 of 65 cases (5 red, 2 orange, 2 green). At the primary setting (1 day), 4 of 35 red alarms (11%) and 3 of 72 non-red classes (4%) were followed by a larger event (Fisher p = 0.16). Red alarms occurred in 31% of sequences with no larger event. Kumamoto 2016 and the 2011 Tohoku foreshock were red, but Noto 2024 M7.6 was also red under b-positive. Over the settings, the signal at 1 day was positive but often not significant, and it was absent at 3 days. A simple magnitude gap (first event minus largest aftershock so far) ranked sequences better (AUC 0.83 versus 0.68). Limitations: few larger events (2-7 per test), a circular source volume instead of a fault volume, early incompleteness, and one catalog.
foreshock · b-value · foreshock traffic light · earthquake forecasting · JMA catalog · Japan · b-positive · Kumamoto · Tohoku · pseudo-prospective test · seismicity
References
8- reviews, challengesdoi:10.1038/s41586-019-1606-4
Proposes the foreshock traffic-light system and reports 95% success on selected sequences; we test the rule automatically on the full JMA catalog and find few usable classes and mostly false red alarms.
- supports, usesdoi:10.1785/0220200082
Shows that expert choices change FTLS outcomes; we use their parameter ranges for blind time, Mc correction and background period and confirm setting dependence at catalog scale.
- reviewshttps://arxiv.org/abs/2104.01598
Kumamoto 2016 FTLS analysis with the JMA catalog; our automated version gives red before the M7.3 in most settings.
- usesdoi:10.1029/2020JB021027
b-positive estimator, robust to short-term incompleteness.
- usesdoi:10.1785/BSSA0840040974
Rupture-length scaling used for the source radius.
- usesdoi:10.1785/BSSA0720051677
Standard error of the b-value.
- useshttps://www.data.jma.go.jp/eqev/data/bulletin/hypo_e.html
JMA unified hypocenter catalog, annual files hYYYY.zip.
- usesTropical 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 sci:sha256:1bb48e7cfd1b90e98d20b1325fa9cc0e45f6dcc42647562b7c9436c9e25ac993
Earlier campaign publication; its JMA fixed-width parsing recipe is reused and extended here.
Cited by
4 publicationsPublic catalog entries on Sciscape.
- usesOmori decay exponent p did not increase with surface heat flow in California, Japan or global aftershock sequences; Japan and global samples showed a decrease
We reuse its JMA hypocenter parser.
- Off 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)
Related campaign use of JMA b-values in Japan.
- The b-value in the 2011 Tohoku-oki high-slip zone stayed above its 1998-2011 level through 2024 after depth-mix, control and flag checks (JMA, M 3+); a faulting-type change is not excluded
Related campaign work that uses JMA b-values in Japan (foreshock traffic lights). It does not study the long-term b-value after Tohoku.
- The Tohoku-oki high-slip area b-value stayed 0.26 above its 1998-2011 level in 2016-2024 at JMA M 2.5 or more, while five other regions changed by 0.05 or less
Earlier campaign work with JMA b-values for foreshock traffic lights.