Abstract:The Epidemic-Type Aftershock Sequences (ETAS) model and its variants effectively capture the space-time clustering of seismicity, setting the standard for earthquake forecasting. Accurate unbiased ETAS calibration is thus crucial. But we identify three sources of bias, (i) boundary effects, (ii) finite-size effects, and (iii) censorship, which are often overlooked or misinterpreted, causing errors in seismic analysis and predictions. By employing an ETAS model variant with variable spatial background rates, we propose a method to correct for these biases, focusing on the branching ratio n, a key indicator of earthquake triggering potential. Our approach quantifies the variation in the apparent branching ratio (napp) with increased cut-off magnitude (Mco) above the optimal cut-off (Mcobest). The napp(Mco) function yields insights superior to traditional point estimates. We validate our method using synthetic earthquake catalogs, accurately recovering the true branching ratio (ntrue) after correcting biases with napp(Mco). Additionally, our method introduces a refined estimation of the minimum triggering magnitude (m0), a crucial parameter in the ETAS model. Applying our framework to the earthquake catalogs of California, New Zealand, and the China Seismic Experimental Site (CSES) in Sichuan and Yunnan provinces, we find that seismicity hovers away from the critical point, nc = 1, remaining distinctly subcritical, however with values tending to be larger than recent reports that do not consider the above biases. It is interesting that, m0 is found around 4 for California, 3 for New Zealand and 2 for CSES, suggesting that many small triggered earthquakes may not be fertile. Understanding seismicity's critical state significantly enhances our comprehension of seismic patterns, aftershock predictability, and informs earthquake risk mitigation and management strategies.

Estimating Real Seismicity Rates from Nonhomogeneously Recorded Earthquake Catalogs

Aftershock forecasts based on incomplete earthquake catalogues: ETASI model application to the 2023 SE Türkiye earthquake sequence

Embracing Data Incompleteness for Better Earthquake Forecasting

Nonstationary ETAS models for nonstandard earthquakes

Nowcasting ETAS Earthquakes: Information Entropy of Earthquake Catalogs

Heterogeneity of Direct Aftershock Productivity of the Mainshock Rupture

Inverse analysis of seismic activity rate changes for severely incomplete sequences: comparison of aftershock activity patterns immediately following the 2023 M6.5 and 2024 M7.6 Noto Peninsula earthquakes

Algorithms for Fitting the Space-Time ETAS Model to Earthquake Catalog Data: A Comparative Study

Increased earthquake rate prior to mainshocks

Bayesian earthquake forecasting approach based on the epidemic type aftershock sequence model

High Resolution Long- and Short-Term Earthquake Forecasts for California

Nowcasting Earthquakes With Stochastic Simulations: Information Entropy of Earthquake Catalogs

Reduced seismic activity after mega earthquakes

Possible Origin of Memory in Earthquakes: Real catalogs and ETAS model

Adaptive estimation of the Gutenberg–Richter b value using a state space model and particle filtering

Number of aftershocks in epidemic-type seismicity models

Revisiting Seismicity Criticality: A New Framework for Bias Correction of Statistical Seismology Model Calibrations

Onset of Aftershocks: Constraints on the Rate‐and‐State Model

A Hypocentral Version of the Space–time ETAS Model

A Test on Methods for MC Estimation Based on Earthquake Catalog

A 3D Spherical ETAS Model: A Case Study in the Alaska-Aleutian Region