The project deals with the assimilation of data with different degrees of uncertainty into the space-time Epidemic Type Aftershock Sequences Model (ETAS). More specifically, we will use instrumental earthquake catalogs as well as historic and paleoseismic data. Because earthquake data are sparse and inhomogeneous, the frequently used maximum likelihood estimates need to be augmented by a quantification of uncertainties. Here we will follow a Bayesian approach that allows one to quantify the relevant parameters along with their uncertainties in space and time. The project starts with the design of appropriate prior distributions with a spatial correlation structure based on known geological features. In this context, the choice of the spatial discretization in relation to the data availability will be a key issue. The ETAS parameters will then be estimated using instrumental catalogs.
Here Markov chain Monte Carlo (MCMC) and sequential data assimilation techniques will be employed. Finally, pre-instrumental data (historic and paleo earthquakes) will also be considered. This requires an error model and its implementation into the ETAS framework. As a result, we will obtain a flexible tool that allows one to adjust the ETAS model to a wide range of seismic data sets and to continuously update the ETAS model through new incoming data. An important application is the forecast of future earthquake rates based on the earthquake history as performed by the Collaboratory for the Study of Earthquake Predictability (CSEP) and used in seismic hazard assessment. In particular, we will be able to add confidence bounds to the forecasted earthquake numbers and to evaluate to what degree existing forecasts and model comparisons are stable.
Reich, S. and Weissmann, S. (2019). Fokker-Planck particle systems for Bayesian inference: Computational approaches.arXiv:1911.10832
Nuesken, N. and Reich, S. (2019). Note on Interacting Langevin diffusions: Gradient structure and ensemble Kalman sampler by Garbuno-Inigo, Hoffmann, Li and Stuart. arXiv:1908.10890
Pathiraja, S. and Reich, S. (2019). Discrete gradients for computational Bayesian inference. Journal of Computational Dynamics (in press) arXiv:1901.06300
Somogyvári, M. and Reich, S. (2019). Convergence tests for transdimensional Markov chains in geoscience imaging, Math Geosci, 2019. doi: 10.1007/s11004-019-09811-x
Lontsi, A. M., García-Jerez, A., Molina-Villegas, J. C., Sánchez-Sesma, F. J., Molkenthin, C., Ohrnberger, M., Krüger, F., Wang, R. and Fäh, D. (2019). A generalized theory for full microtremor horizontal-to-vertical [H/V(z, f)] spectral ratio interpretation in offshore and onshore environments, Geophysical Journal International, 218(2), 1276–1297. doi: 10.1093/gji/ggz223 arXiv: 1907.04606
Zöller, G. (2018). A Statistical Model for Earthquake Recurrence Based on the Assimilation of Paleoseismicity, Historic Seismicity, and Instrumental Seismicity. Journal of Geophysical Research: Solid Earth, 123, 4906-4921. doi: 10.1029/2017JB015099
Fiedler, B., Zöller, G., Holschneider, M. and Hainzl, S. (2018). Multiple Change‐Point Detection in Spatiotemporal Seismicity Data, Bulletin of the Seismological Society of America. 108 (3A): 1147-1159. doi: 10.1785/0120170236