Nonstationary ETAS models for nonstandard earthquakes

The conditional intensity function of a point process is a useful tool for generating probability forecasts of earthquakes. The epidemic-type aftershock sequence (ETAS) model is defined by a conditional intensity function, and the corresponding point process is equivalent to a branching process, assuming that an earthquake generates a cluster of offspring earthquakes (triggered earthquakes or so-called aftershocks). Further, the size of the first-generation cluster depends on the magnitude of the triggering (parent) earthquake. The ETAS model provides a good fit to standard earthquake occurrences. However, there are nonstandard earthquake series that appear under transient stress changes caused by aseismic forces such as volcanic magma or fluid intrusions. These events trigger transient nonstandard earthquake swarms, and they are poorly fitted by the stationary ETAS model. In this study, we examine nonstationary extensions of the ETAS model that cover nonstandard cases. These models allow the parameters to be time-dependent and can be estimated by the empirical Bayes method. The best model is selected among the competing models to provide the inversion solutions of nonstationary changes. To address issues of the uniqueness and robustness of the inversion procedure, this method is demonstrated on an inland swarm activity induced by the 2011 Tohoku-Oki, Japan earthquake of magnitude 9.0.

[1]  Y. Ogata.  Seismicity Analysis through Point-process Modeling: A Review , 1999 .

[2]  Nicholas Deichmann,et al.  High fluid pressure and triggered earthquakes in the enhanced geothermal system in Basel, Switzerland , 2011 .

[3]  Yosihiko Ogata,et al.  Seismicity anomaly scenario prior to the major recurrent earthquakes off the east coast of Miyagi Prefecture, northern Japan , 2006 .

[4]  Chihiro Hashimoto,et al.  Changes in seismic activity following the 2011 Tohoku-oki earthquake: Effects of pore fluid pressure , 2012 .

[5]  T. Utsu.  Aftershocks and Earthquake Statistics(2) : Further Investigation of Aftershocks and Other Earthquake Sequences Based on a New Classification of Earthquake Sequences , 1971 .

[6]  P. Reasenberg.  Second‐order moment of central California seismicity, 1969–1982 , 1985 .

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

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

[9]  Gerassimos A. Papadopoulos,et al.  OPERATIONAL EARTHQUAKE FORECASTING. State of Knowledge and Guidelines for Utilization , 2011 .

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

[11]  T. Utsu,et al.  A Relation between the Area of After-shock Region and the Energy of Main-shock , 1955 .

[12]  Ramesh Sharda,et al.  Modeling Brand Post Popularity in Online Social Networks , 2012 .

[13]  Mitsuhiro Matsu'ura,et al.  The 3‐D tectonic stress fields in and around Japan inverted from centroid moment tensor data of seismic events , 2010 .

[14]  S. Toda,et al.  Using the 2011 M = 9 . 0 Tohoku earthquake to test the Coulomb stress triggering hypothesis and to calculate faults brought closer to failure , 2011 .

[15]  Yosihiko Ogata,et al.  Modeling seismic swarms triggered by aseismic transients , 2009 .

[16]  V. Chavez-Demoulin,et al.  High-frequency financial data modeling using Hawkes processes , 2012 .

[17]  Lucile M. Jones,et al.  When and where the aftershock activity was depressed: Contrasting decay patterns of the proximate large earthquakes in southern California , 2003 .

[18]  Yosihiko Ogata,et al.  Statistical Models for Earthquake Occurrences and Residual Analysis for Point Processes , 1988 .

[19]  D. Kendall.  Stochastic Processes and Population Growth , 1949 .

[20]  Y. Ogata.  The asymptotic behaviour of maximum likelihood estimators for stationary point processes , 1978 .

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

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

[23]  Warner Marzocchi,et al.  Operational Earthquake Forecasting and Decision-Making , 2014 .

[24]  R. Herrera,et al.  Self-exciting Extreme Value Models for Stock Market Crashes , 2009 .

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

[26]  Jiancang Zhuang,et al.  Properties of the probability distribution associated with the largest event in an earthquake cluster and their implications to foreshocks. , 2006, Physical review. E, Statistical, nonlinear, and soft matter physics.

[27]  H. Akaike.  Factor analysis and AIC , 1987 .

[28]  Y. Ogata.  Significant improvements of the space-time ETAS model for forecasting of accurate baseline seismicity , 2011 .

[29]  I. Good,et al.  The Estimation of Probabilities: An Essay on Modern Bayesian Methods. , 1967 .

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

[31]  D. Vere-Jones,et al.  A statistical survey of earthquakes in the main seismic region of New Zealand: Part 2—Time series analyses , 1966 .

[32]  George E. Tita,et al.  Self-Exciting Point Process Modeling of Crime , 2011 .

[33]  Y. Ogata,et al.  Modelling heterogeneous space–time occurrences of earthquakes and its residual analysis , 2003 .

[34]  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 .

[35]  Yosihiko Ogata,et al.  Exploratory analysis of earthquake clusters by likelihood-based trigger models , 2001, Journal of Applied Probability.

[36]  Yosihiko Ogata,et al.  Statistical model for standard seismicity and detection of anomalies by residual analysis , 1989 .

[37]  T. Utsu.  A statistical study on the occurrence of aftershocks. , 1961 .

[38]  Frederic Paik Schoenberg,et al.  Application of Branching Models in the Study of Invasive Species , 2012 .

[39]  T. Utsu.  Aftershocks and Earthquake Statistics (4) : Analyses of the Distribution of Earthquakes in Magnitude, Time and Space with Special Consideration to Clustering Characteristics of Earthquake Occurrence(2) , 1972 .

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

[41]  Y. Ogata,et al.  Precursory seismic anomalies and transient crustal deformation prior to the 2008 Mw = 6.9 Iwate‐Miyagi Nairiku, Japan, earthquake , 2010 .

[42]  Roger D. Peng,et al.  On the distribution of wildfire sizes , 2003 .

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

[44]  Yosihiko Ogata,et al.  A non‐stationary epidemic type aftershock sequence model for seismicity prior to the December 26, 2004 M 9.1 Sumatra‐Andaman Islands mega‐earthquake , 2013 .

[45]  Yosihiko Ogata,et al.  Detection of anomalous seismicity as a stress change sensor , 2003 .

[46]  Y. Ogata.  Anomalies of Seismic Activity and Transient Crustal Deformations Preceding the 2005 M 7.0 Earthquake West of Fukuoka , 2010 .

[47]  Y. Ogata,et al.  The Centenary of the Omori Formula for a Decay Law of Aftershock Activity , 1995 .

[48]  C. Lomnitz.  Global tectonics and earthquake risk , 1974 .

[49]  T. Utsu.  Aftershocks and Earthquake Statistics(1) : Some Parameters Which Characterize an Aftershock Sequence and Their Interrelations , 1970 .

[50]  Roger D. Peng,et al.  A Space–Time Conditional Intensity Model for Evaluating a Wildfire Hazard Index , 2005 .

[51]  Daryl J. Daley,et al.  An Introduction to the Theory of Point Processes , 2013 .

[52]  Shinji Toda,et al.  Widespread seismicity excitation throughout central Japan following the 2011 M=9.0 Tohoku earthquake and its interpretation by Coulomb stress transfer , 2011 .

[53]  T. Utsu.  On the nature of three Alaskan aftershock sequences of 1957 and 1958 , 1962 .

[54]  Satoshi Ide,et al.  Statistic analysis of swarm activities around the Boso Peninsula, Japan: Slow slip events beneath Tokyo Bay? , 2011 .

[55]  L. Adamopoulos.  Cluster models for earthquakes: Regional comparisons , 1976 .

[56]  Yosihiko Ogata,et al.  Immediate and updated forecasting of aftershock hazard , 2006 .

[57]  Robert Tibshirani,et al.  An Introduction to the Bootstrap , 1994 .

[58]  R. A. Gaskins,et al.  Nonparametric roughness penalties for probability densities , 2022 .

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

[60]  Kazuyuki Aihara,et al.  Forecasting large aftershocks within one day after the main shock , 2013, Scientific reports.

[61]  Yosihiko Ogata,et al.  Detection of precursory relative quiescence before great earthquakes through a statistical model , 1992 .

[62]  Y. Ogata,et al.  Analysis of temporal and spatial heterogeneity of magnitude frequency distribution inferred from earthquake catalogues , 1993 .

[63]  H. Akaike,et al.  Information Theory and an Extension of the Maximum Likelihood Principle , 1973 .

[64]  Yosihiko Ogata,et al.  Seismicity and geodetic anomalies in a wide area preceding the Niigata-Ken-Chuetsu earthquake of 23 October 2004, central Japan , 2007 .

[65]  Hirotugu Akaike,et al.  On entropy maximization principle , 1977 .

[66]  L. Knopoff,et al.  Statistical Short-Term Earthquake Prediction , 1987, Science.

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

[68]  P. Laplace.  Memoir on the Probability of the Causes of Events , 1986 .

[69]  Giada Adelfio,et al.  Hybrid kernel estimates of space–time earthquake occurrence rates using the epidemic-type aftershock sequence model , 2009 .

[70]  Yoshiyuki Yamamoto,et al.  Change detection method for landslide area using RGB color composite image of SRTM and ALOS/PALSAR InSAR DEM: A case study of the Iwate-Miyagi Nairiku Earthquake in 2008 , 2011, 2011 IEEE International Geoscience and Remote Sensing Symposium.