A Machine Learning Approach for Parameter Screening in Earthquake Simulation

Earthquakes are the result of rupture in the Earth's crust. The rupture process is difficult to model deterministically due to the number of unmeasurable parameters involved and poorly constrained physical conditions, as well as the very diverse scales involved in their nucleation (meters) and complete failure (up to hundreds of kilometers). In this research work we focus on synthetic seismic catalogs generated with a stochastic modeling technique called Fiber Bundle Model (FBM). These catalogs can be readily compared with statistical measures computed from real earthquake series, but the link between the FBM parameters and the characteristics of the obtained earthquake series is difficult to assess. Furthermore, the stochastic nature of the model requires a large amount of realizations in order to attain statistical robustness. The aim of this work is to estimate the FBM parameters that generate aftershock sequences that are similar to those generated by real seismic events. In order to estimate the optimal combination of parameters that generate such sequences, we executed a large number of simulations with different combinations of parameters using High-Performance Computing (HPC) resources to reduce compute time. Lastly, the synthetic datasets were used to train a supervised Machine Learning (ML) model that analyzes and extracts statistical patterns that reproduce the observations regarding aftershock occurrence and its spatio-temporal distribution in real seismic events.

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

[2]  Y. Morenoa,et al.  Self-organized criticality in a bre-bundle-type model , 1999 .

[3]  E. H. Lloyd,et al.  Long-Term Storage: An Experimental Study. , 1966 .

[4]  Trevor Hastie,et al.  The Elements of Statistical Learning , 2001 .

[5]  X. Lana,et al.  Aftershock sequences of three seismic crises at southern California, USA, simulated by a cellular automata model based on self-organized criticality , 2015, Geosciences Journal.

[6]  Yoshida Hiroaki,et al.  Rapid Feature Selection Based on Random Forests for High-Dimensional Data , 2012 .

[7]  Kurt Hornik,et al.  Misc Functions of the Department of Statistics (e1071), TU Wien , 2014 .

[8]  R. Tibshirani,et al.  Flexible Discriminant Analysis by Optimal Scoring , 1994 .

[9]  Luciano Telesca,et al.  Tsallis-Based Nonextensive Analysis of the Southern California Seismicity , 2011, Entropy.

[10]  Eric R. Ziegel,et al.  The Elements of Statistical Learning , 2003, Technometrics.

[11]  R. Tibshirani,et al.  Penalized Discriminant Analysis , 1995 .

[12]  Ron Kohavi,et al.  Guest Editors' Introduction: On Applied Research in Machine Learning , 1998, Machine Learning.

[13]  Leo Breiman,et al.  Random Forests , 2001, Machine Learning.

[14]  B. Gutenberg,et al.  Seismicity of the Earth and associated phenomena , 1950, MAUSAM.

[15]  Matteo Hessel,et al.  Machine Learning for Parameter Screening in Computer Simulations , 2014, MESAS.

[16]  G. Atkinson,et al.  Aftershock Statistics for Earthquakes in the St. Lawrence Valley , 2014 .

[17]  Pier Luca Lanzi,et al.  Automatic Tuning of Computational Models , 2014, SIMULTECH.

[18]  Á. Figueroa-Soto,et al.  Simulation of spatial and temporal properties of aftershocks by means of the fiber bundle model , 2017, Journal of Seismology.

[19]  Charles F. Richter,et al.  Earthquake magnitude, intensity, energy, and acceleration , 1942 .

[20]  D. Turcotte,et al.  Aftershock Statistics , 2005 .

[21]  G. King,et al.  STATIC STRESS CHANGES AND THE TRIGGERING OF EARTHQUAKES , 1994 .