A way to synchronize models with seismic faults for earthquake forecasting: Insights from a simple stochastic model

Abstract

Numerical models are starting to be used for determining the future behaviour of seismic faults and fault networks. Their final goal would be to forecast future large earthquakes. In order to use them for this task, it is necessary to synchronize each model with the current status of the actual fault or fault network it simulates (just as, for example, meteorologists synchronize their models with the atmosphere by incorporating current atmospheric data in them). However, lithospheric dynamics is largely unobservable: important parameters cannot (or can rarely) be measured in Nature. Earthquakes, though, provide indirect but measurable clues of the stress and strain status in the lithosphere, which should be helpful for the synchronization of the models.The rupture area is one of the measurable parameters of earthquakes. Here we explore how it can be used to at least synchronize fault models between themselves and forecast synthetic earthquakes. Our purpose here is to forecast synthetic earthquakes in a simple but stochastic (random) fault model. By imposing the rupture area of the synthetic earthquakes of this model on other models, the latter become partially synchronized with the first one. We use these partially synchronized models to successfully forecast most of the largest earthquakes generated by the first model. This forecasting strategy outperforms others that only take into account the earthquake series. Our results suggest that probably a good way to synchronize more detailed models with real faults is to force them to reproduce the sequence of previous earthquake ruptures on the faults. This hypothesis could be tested in the future with more detailed models and actual seismic data.

Predicting Solar Flares by Data Assimilation in Avalanche Models

“…However, as pointed out by Isliker, Anastasiadis, and Vlahos (2000), it is still possible to compute derivatives and thus operators. To the best of our knowledge, at the present time the only other area in which data assimilation techniques are being developed in conjunction with a cellular automaton is in seismic data assimilation of a stochastic random fault model (Rundle et al, 2003;González et al, 2006). We are also using a time series of a global, model-produced variable to define the error, as opposed to the spatial state of the system measured at some time interval beyond the initial condition in which data is being assimilated.…”

Brittle Failure Mechanism of Multiple Locked Patches in a Seismogenic Fault and Exploration of a New Approach to Earthquake Prediction

“…(1) Based on seismological observations, statistical physics, experiments and mechanical theory, the mechanism of instability evolution is analyzed which allows event prediction to be made. However, due to strong nonlinearity of the seismogenic process, the kinetic equations [11] describing such a process cannot be accurately established. Even if these equations were written out, it would be still rather difficult to determine many geometrical and mechanical parameters.…”