paper

The 2016 Mw 7.0 Kumamoto Earthquake Sequence, Japan revisited: Insights from Spatio-Temporal Analysis of Seismicity Parameters

arXiv:2609.06911

Abstract

Understanding how crustal faults accumulate strain, nucleate ruptures, and redistribute post-seismic stress is fundamental to seismic hazard assessment. The 16 April 2016 7.0 Kumamoto earthquake, rupturing the Futagawa--Hinagu fault zone in the Beppu--Shimabara graben, central Kyushu, Japan, offers a premier dataset for tracking these processes across a full earthquake cycle. Using the JMA catalog (2014--2018) with dynamic completeness estimation (--), we examine the Gutenberg--Richter -value, the 3-D hypocentral fractal dimension (), and the seismicity-rate anomaly (-value). During the one- to two-year preparatory phase, declined from a baseline of to a precursory minimum of , while contracted from to --, recording stress concentration and microfracture coalescence onto a narrow nucleation zone. A coherent negative -value anomaly (-1.6 to -2.0) developed along the graben, strengthening with integration time---consistent with, though not proof of, progressive fault locking. Depth-sliced volumes show the low- locked core () was stratified at --~km depth and sharpened within the final four months before failure. Rupture reversed this within weeks: surged to -- and expanded to --, consistent with coseismic stress drop and aftershock activation, followed by recovery over 2.5 years. Cross-sections reveal sharp aftershock localization alongside two unrelaxed asperities (--) near the Hinagu termination and Mount Aso, consistent with positive Coulomb stress loading. This framework resolves asperity locking, release, and healing better than any single metric; however, since anomalies were identified retrospectively, they should be read as evidence of coherent behavior rather than a validated forecast tool alone.