Next-day earthquake forecasts for the Japan region generated by the ETAS model

This paper gives the technical solutions of implementing the space-time epidemic-type aftershock sequence (ETAS) model for short-term (1-day) earthquake forecasts for the all-Japan region in the Collaboratory for the Study of Earthquake Predictability (CSEP) project in Japan. For illustration, a retrospective forecasting experiment is carried out to forecast the seismicity in the Japan region before and after the Tokachi-Oki earthquake (M 8.0) at 19:50:07 (UTC) on 25 September 2003, in the format of contour images. The optimal model parameters used for the forecasts are estimated by fitting the model to the observation records up to the starting time of the forecasting period, and the probabilities of earthquake occurrences are obtained through simulations. To tackle the difficulty of heavy computations in fitting a complicated point-process to a huge dataset, an “off-line optimization” and “online forecasting” scheme is proposed to keep both the estimates of model parameters and forecasts updated according to the most recent observations. The results show that the forecasts have captured the spatial distribution and temporal evolution of the features of future seismicity. These forecasts are tested against the reference Poisson model that is stationary in time but spatially inhomogeneous.

[1]  D. Vere-Jones,et al.  Analyzing earthquake clustering features by using stochastic reconstruction , 2004 .

[2]  Annemarie Christophersen,et al.  Differences between spontaneous and triggered earthquakes: Their influences on foreshock probabilities , 2008 .

[3]  Warner Marzocchi,et al.  On the Increase of Background Seismicity Rate during the 1997-1998 Umbria-Marche, Central Italy, Sequence: Apparent Variation or Fluid-Driven Triggering? , 2010 .

[4]  Missing links in earthquake clustering models , 2010 .

[5]  Yosihiko Ogata,et al.  On Lewis' simulation method for point processes , 1981, IEEE Trans. Inf. Theory.

[6]  Foreshocks explained by cascades of triggered seismicity , 2002, physics/0210130.

[7]  Y. Kagan,et al.  High-resolution Time-independent Grid-based Forecast for M ≥ 5 Earthquakes in California , 2007 .

[8]  D. Vere-Jones,et al.  A Variable-Grid Algorithm for Smoothing Clustered Data , 1986 .

[9]  Jiancang Zhuang,et al.  Gambling scores for earthquake predictions and forecasts , 2010 .

[10]  Yan Y. Kagan,et al.  Comparison of short-term and long-term earthquake forecast models for southern California , 2006 .

[11]  Yosihiko Ogata,et al.  Detecting fluid signals in seismicity data through statistical earthquake modeling , 2005 .

[12]  Yosihiko Ogata,et al.  Space‐time model for regional seismicity and detection of crustal stress changes , 2004 .

[13]  Jiancang Zhuang,et al.  Second‐order residual analysis of spatiotemporal point processes and applications in model evaluation , 2006 .

[14]  Jiancang Zhuang,et al.  Space–time ETAS models and an improved extension , 2006 .

[15]  D. Vere-Jones,et al.  Stochastic Declustering of Space-Time Earthquake Occurrences , 2002 .

[16]  Y. Kagan,et al.  Comparison of Short-Term and Time-Independent Earthquake Forecast Models for Southern California , 2006 .

[17]  Y. Ogata.  Space-Time Point-Process Models for Earthquake Occurrences , 1998 .

[18]  G. M. Molchan,et al.  Strategies in strong earthquake prediction , 1990 .

[19]  D. Sornette,et al.  Are aftershocks of large Californian earthquakes diffusing , 2003, physics/0303075.

[20]  Temporal behavior of the background seismicity rate in central Japan, 1998 to mid-2003 , 2006 .

[21]  Rodolfo Console,et al.  Real Time Forecasts through an Earthquake Clustering Model Constrained by the Rate-and- State Constitutive Law: Comparison with a Purely Stochastic ETAS Model , 2007 .

[22]  Warner Marzocchi,et al.  Real‐time forecasting following a damaging earthquake , 2009 .

[23]  R. Dennis Cook,et al.  Cross-Validation of Regression Models , 1984 .

[24]  Rodolfo Console,et al.  Refining earthquake clustering models , 2003 .

[25]  A. Hawkes.  Point Spectra of Some Mutually Exciting Point Processes , 1971 .

[26]  Naoshi Hirata,et al.  Overview of the first earthquake forecast testing experiment in Japan , 2011 .

[27]  Yan Y. Kagan,et al.  Statistical search for non-random features of the seismicity of strong earthquakes , 1976 .

[28]  D. Vere-Jones.  Stochastic Models for Earthquake Occurrence , 1970 .

[29]  Rodolfo Console,et al.  Probability gains of an epidemic-type aftershock sequence model in retrospective forecasting of M ≥ 5 earthquakes in Italy , 2009 .

[30]  Chung-Pai Chang,et al.  A study on the background and clustering seismicity in the Taiwan region by using point process models : Stress transfer, earthquake triggering, and time-dependent seismic hazard , 2005 .

[31]  Roger Fletcher,et al.  A Rapidly Convergent Descent Method for Minimization , 1963, Comput. J..

[32]  Jie Liu,et al.  An evaluation of Chinese annual earthquake predictions, 1990–1998 , 2001, Journal of Applied Probability.

[33]  Richard S. Bird,et al.  An introduction to the theory of lists , 1987 .

[34]  Y. Ogata.  Seismicity quiescence and activation in western Japan associated with the 1944 and 1946 great earthquakes near the Nankai trough , 2004 .

[35]  David Vere-Jones,et al.  The Entropy Score and its Uses in Earthquake Forecasting , 2005 .

[36]  K. Murase.  A Characteristic Change in Fractal Dimension Prior to the 2003 Tokachi-oki Earthquake (MJ = 8.0), Hokkaido, Northern Japan , 2004 .

[37]  David Vere-Jones.  Probabilities and Information Gain for Earthquake Forecasting , 2013 .

[38]  A. Hawkes.  Spectra of some self-exciting and mutually exciting point processes , 1971 .

[39]  Rodolfo Console,et al.  Real time earthquake forecasting in Italy , 2009 .