Numerical shake prediction for Earthquake Early Warning: data assimilation, real-time shake-mapping, and simulation of wave propagation

Many of the present earthquake early warning (EEW) systems quickly determine an event’s hypocenter and magnitude and then predict strengths of ground motions. The M w 9.0 Tohoku earthquake, however, revealed some technical issues with such methods: (1) underprediction at large distances due to the large extent of the fault rupture and (2) overprediction because the system was confused by multiple aftershocks that occurred simultaneously. To address these issues, we propose a new concept for EEW, in which the distribution of the present wavefield is estimated precisely in real time (real‐time shake mapping) by applying a data assimilation technique, and then the future wavefield is predicted time evolutionally by simulation of seismic‐wave propagation. Information on the hypocenter location and magnitude is not necessarily required. We call this method, in which physical processes are simulated from the precisely estimated present condition, numerical shake prediction because of its analogy to numerical weather prediction in meteorology. By applying the proposed method to the 2011 Tohoku earthquake and the 2004 Mid‐Niigata Prefecture earthquake ( M w 6.7), we show that numerical shake prediction can precisely and rapidly predict ground motion in real time. Online Material: Animations as examples of numerical shake prediction.

[1]  Mitsuyuki Hoshiba,et al.  Separation of scattering attenuation and intrinsic absorption in Japan using the multiple lapse time window analysis of full seismogram envelope , 1993 .

[2]  H. Nakamura,et al.  How precisely can we anticipate seismic intensities? A study of uncertainty of anticipated seismic intensities for the Earthquake Early Warning method in Japan , 2010 .

[3]  A. Uribe,et al.  The Mexican Seismic Alert System (SASMEX): Its Alert Signals, Broadcast Results and Performance During the M 7.4 Punta Maldonado Earthquake of March 20th, 2012 , 2014 .

[4]  Mitsuyuki Hoshiba,et al.  Real‐time prediction of ground motion by Kirchhoff‐Fresnel boundary integral equation method: Extended front detection method for Earthquake Early Warning , 2013 .

[5]  Mustafa Erdik,et al.  A study on warning algorithms for Istanbul earthquake early warning system , 2009 .

[6]  Hiroyuki Fujiwara,et al.  A Real-Time Processing of Seismic Intensity , 2008 .

[7]  G. Manfredi,et al.  Development of Earthquake Early Warning Systems in the European Union , 2014 .

[8]  T. Ozaki,et al.  Earthquake Early Warning and Tsunami Warning of the Japan Meteorological Agency, and Their Performance in the 2011 off the Pacific Coast of Tohoku Earthquake ({M}_{{\mathrm{w}}} 9.0) , 2014 .

[9]  Scattering attenuation and intrinsic absorption using uniform and depth dependent model – Application to full seismogram envelope recorded in Northern Chile , 2001 .

[10]  M. Hoshiba Simulation of multiple-scattered coda wave excitation based on the energy conservation law , 1991 .

[11]  Masumi Yamada,et al.  Estimation of Fault Rupture Extent Using Near-Source Records for Earthquake Early Warning , 2014 .

[12]  Richard M. Allen,et al.  Designing a Network‐Based Earthquake Early Warning Algorithm for California: ElarmS‐2 , 2014 .

[13]  Michael Fehler,et al.  Seismic Wave Propagation and Scattering in the Heterogeneous Earth , 2012 .

[14]  Mitsuyuki Hoshiba Seismic coda wave envelope in depth-dependent S wave velocity structure , 1997 .

[15]  Kazuo Yoshimoto,et al.  Monte Carlo simulation of seismogram envelopes in scattering media , 2000 .

[16]  Masumi Yamada,et al.  Bayesian Approach for Identification of Multiple Events in an Early Warning System , 2014 .

[17]  Mitsuyuki Hoshiba,et al.  Outline of the 2011 off the Pacific coast of Tohoku Earthquake (Mw 9.0) —Earthquake Early Warning and observed seismic intensity— , 2011 .

[18]  Shigeki Aoki,et al.  Real-time site correction based on evaluating relative responses to common reference station for wide area network , 2014 .

[19]  Mitsuyuki Hoshiba Real‐Time Correction of Frequency‐Dependent Site Amplification Factors for Application to Earthquake Early Warning , 2013 .

[20]  Eugenia Kalnay,et al.  Atmospheric Modeling, Data Assimilation and Predictability , 2002 .

[21]  M. Fehler,et al.  Separation of scattering and intrinsic attenuation for the Kanto-Tokai region, Japan, using measurements of S-wave energy versus hypocentral distance , 1992 .

[22]  Walter H. F. Smith,et al.  New version of the generic mapping tools , 1995 .

[23]  Thomas H. Heaton,et al.  Real-time Finite Fault Rupture Detector (FinDer) for large earthquakes , 2012 .

[24]  Nai-Chi Hsiao,et al.  Development of earthquake early warning system in Taiwan , 2009 .

[25]  Makoto Saito,et al.  Earthquake Early Warning Starts Nationwide in Japan , 2008 .

[26]  Mustafa Erdik,et al.  PreSEIS: A Neural Network-Based Approach to Earthquake Early Warning for Finite Faults , 2008 .

[27]  Akira Ishimaru,et al.  Wave propagation and scattering in random media , 1997 .

[28]  Source model for strong ground motion generation in the frequency range 0.1–10 Hz during the 2011 Tohoku earthquake , 2012, Earth, Planets and Space.

[29]  Mitsuyuki Hoshiba,et al.  Estimation of nonisotropic scattering in western Japan using coda wave envelopes : application of a multiple nonisotropic scattering model , 1995 .

[30]  A. A. Gusev,et al.  Monte-Carlo simulation of record envelope of a near earthquake , 1987 .

[31]  K. Irikura,et al.  Short‐Period Source Model of the 2011 Mw 9.0 Off the Pacific Coast of Tohoku Earthquake , 2013 .