Bayesian forecasting of recurrent earthquakes and predictive performance for a small sample size

[1] This paper presents a Bayesian method of probability forecasting for a renewal of earthquakes. When only limited records of characteristic earthquakes on a fault are available, relevant prior distributions for renewal model parameters are essential to computing unbiased, stable time-dependent earthquake probabilities. We also use event slip and geological slip rate data combined with historical earthquake records to improve our forecast model. We apply the Brownian Passage Time (BPT) model and make use of the best fit prior distribution for its coefficient of variation (the shape parameter, alpha) relative to the mean recurrence time because the Earthquake Research Committee (ERC) of Japan uses the BPT model for long-term forecasting. Currently, more than 110 active faults have been evaluated by the ERC, but most include very few paleoseismic events. We objectively select the prior distribution with the Akaike Bayesian Information Criterion using all available recurrence data including the ERC datasets. These data also include mean recurrence times estimated from slip per event divided by long-term slip rate. By comparing the goodness of fit to the historical record and simulated data, we show that the proposed predictor provides more stable performance than plug-in predictors, such as maximum likelihood estimates and the predictor currently adopted by the ERC.

[1]  Yosihiko Ogata.  Biases and Uncertainties when Estimating the Hazard of the Next Nankai Earthquake , 2001 .

[2]  Yosihiko Ogata,et al.  Estimating the hazard of rupture using uncertain occurrence times of paleoearthquakes , 1999 .

[3]  Yosihiko Ogata.  Slip‐size‐dependent renewal processes and Bayesian inferences for uncertainties , 2002 .

[4]  André Jalobeanu,et al.  Long‐term slip history discriminates among occurrence models for seismic hazard assessment , 2010 .

[5]  Tom Parsons.  Monte Carlo method for determining earthquake recurrence parameters from short paleoseismic catalogs: Example calculations for California , 2008 .

[6]  Tom Parsons,et al.  Earthquake recurrence on the south Hayward fault is most consistent with a time dependent, renewal process , 2008 .

[7]  H. Akaike.  Prediction and Entropy , 1985 .

[8]  David A. Rhoades,et al.  On the handling of uncertainties in estimating the hazard of rupture on a fault segment , 1994 .

[9]  J. Aitchison.  Goodness of prediction fit , 1975 .

[10]  Takashi Nakata,et al.  Time‐predictable recurrence model for large earthquakes , 1980 .

[11]  Gert Zöller,et al.  Recurrence Time Distributions of Large Earthquakes in a Stochastic Model for Coupled Fault Systems: The Role of Fault Interaction , 2007 .

[12]  D. P. Schwartz,et al.  Fault behavior and characteristic earthquakes: Examples from the Wasatch and San Andreas Fault Zones , 1984 .

[13]  Yan Y. Kagan,et al.  Random stress and earthquake statistics: spatial dependence , 1987 .

[14]  David D. Jackson,et al.  Seismic hazards in southern California: probable earthquakes, 1994 to 2024 , 1996 .

[15]  W. H. Bakun,et al.  The Parkfield, California, Earthquake Prediction Experiment , 1985, Science.

[16]  Glenn P. Biasi,et al.  Quasi-periodic recurrence of large earthquakes on the southern San Andreas fault , 2010 .

[17]  Hirotugu Akaike,et al.  On the Likelihood of a Time Series Model , 1978 .

[18]  Stuart P. Nishenko,et al.  Circum-Pacific seismic potential: 1989–1999 , 1991 .

[19]  Domenico Giardini,et al.  Multiarchive paleoseismic record of late Pleistocene and Holocene strong earthquakes in Switzerland , 2005 .

[20]  Hirotugu Akaike,et al.  Likelihood and the Bayes procedure , 1980 .

[21]  R. A. Leibler,et al.  On Information and Sufficiency , 1951 .

[22]  S. P. Nishenko,et al.  A generic recurrence interval distribution for earthquake forecasting , 1987 .

[23]  John J. Clague,et al.  Evidence for large earthquakes at the Cascadia Subduction Zone , 1997 .

[24]  H. Akaike.  A new look at the statistical model identification , 1974 .

[25]  Masami Okada,et al.  A Prior Distribution of the Parameters in the Renewal Model with Lognormal Distribution Used for Estimating the Probability of Recurrent Earthquakes , 2007 .

[26]  Akira Hasegawa,et al.  Characteristic small‐earthquake sequence off Sanriku, northeastern Honshu, Japan , 2002 .

[27]  Tetsuo Hashimoto,et al.  Detection of the Characteristic Earthquake Activity of the Middle Class Earthquakes in the Tohoku District, Japan , 2005 .

[28]  Tom Parsons,et al.  Significance of stress transfer in time‐dependent earthquake probability calculations , 2004 .

[29]  M. Wyss,et al.  Regular Intervals Between Hawaiian Earthquakes: Implications for Predicting the Next Event , 1986, Science.

[30]  David Vere-Jones,et al.  Earthquake prediction - a statistician's view. , 1978 .

[31]  Rodolfo Console,et al.  Stress interaction effect on the occurrence probability of characteristic earthquakes in Central Apennines , 2008 .

[32]  Yehuda Ben-Zion,et al.  Estimating recurrence times and seismic hazard of large earthquakes on an individual fault , 2007 .

[33]  William L. Ellsworth,et al.  A Brownian Model for Recurrent Earthquakes , 2002 .