{"status":"ok","message-type":"work","message-version":"1.0.0","message":{"indexed":{"date-parts":[[2026,4,22]],"date-time":"2026-04-22T07:14:10Z","timestamp":1776842050875,"version":"3.51.2"},"reference-count":31,"publisher":"Frontiers Media SA","license":[{"start":{"date-parts":[[2023,3,23]],"date-time":"2023-03-23T00:00:00Z","timestamp":1679529600000},"content-version":"vor","delay-in-days":0,"URL":"https:\/\/creativecommons.org\/licenses\/by\/4.0\/"}],"funder":[{"DOI":"10.13039\/100010661","name":"Horizon 2020 Framework Programme","doi-asserted-by":"publisher","id":[{"id":"10.13039\/100010661","id-type":"DOI","asserted-by":"publisher"}]},{"DOI":"10.13039\/501100000270","name":"Natural Environment Research Council","doi-asserted-by":"publisher","id":[{"id":"10.13039\/501100000270","id-type":"DOI","asserted-by":"publisher"}]}],"content-domain":{"domain":["frontiersin.org"],"crossmark-restriction":true},"short-container-title":["Front. Appl. Math. Stat."],"abstract":"<jats:p>The epidemic type aftershock sequence (ETAS) model is widely used to model seismic sequences and underpins operational earthquake forecasting (OEF). However, it remains challenging to assess the reliability of inverted ETAS parameters for numerous reasons. For example, the most common algorithms just return point estimates with little quantification of uncertainty. At the same time, Bayesian Markov chain Monte Carlo implementations remain slow to run and do not scale well, and few have been extended to include spatial structure. This makes it difficult to explore the effects of stochastic uncertainty. Here, we present a new approach to ETAS modeling using an alternative Bayesian method, the integrated nested Laplace approximation (INLA). We have implemented this model in a new R-Package called <jats:monospace>ETAS.inlabru<\/jats:monospace>, which is built on the R packages R-INLA and <jats:monospace>inlabru<\/jats:monospace>. Our study has included extending these packages, which provided tools for modeling log-Gaussian Cox processes, to include the self-exciting Hawkes process that ETAS is a special case of. While we just present the temporal component here, the model scales to a spatio-temporal model and may include a variety of spatial covariates. This is a fast method that returns joint posteriors on the ETAS background and triggering parameters. Using a series of synthetic case studies, we explore the robustness of ETAS inversions using this method of inversion. We also included runnable notebooks to reproduce the figures in this article as part of the package's GitHub repository. We demonstrate that reliable estimates of the model parameters require that the catalog data contain periods of relative quiescence, as well as triggered sequences. We explore the robustness of the method under stochastic uncertainty in the training data and show that the method is robust to a wide range of starting conditions. We show how the inclusion of historic earthquakes prior to the modeled time window affects the quality of the inversion. Finally, we show that rate-dependent incompleteness of earthquake catalogs after large earthquakes have a significant and detrimental effect on the ETAS posteriors. We believe that the speed of the <jats:monospace>inlabru<\/jats:monospace> inversion, which includes a rigorous estimation of uncertainty, will enable a deeper exploration of how to use ETAS robustly for seismicity modeling and operational earthquake forecasting.<\/jats:p>","DOI":"10.3389\/fams.2023.1126759","type":"journal-article","created":{"date-parts":[[2023,3,23]],"date-time":"2023-03-23T06:18:03Z","timestamp":1679552283000},"update-policy":"https:\/\/doi.org\/10.3389\/crossmark-policy","source":"Crossref","is-referenced-by-count":9,"title":["Bayesian modeling of the temporal evolution of seismicity using the ETAS.inlabru package"],"prefix":"10.3389","volume":"9","author":[{"given":"Mark","family":"Naylor","sequence":"first","affiliation":[],"role":[{"role":"author","vocabulary":"crossref"}]},{"given":"Francesco","family":"Serafini","sequence":"additional","affiliation":[],"role":[{"role":"author","vocabulary":"crossref"}]},{"given":"Finn","family":"Lindgren","sequence":"additional","affiliation":[],"role":[{"role":"author","vocabulary":"crossref"}]},{"given":"Ian G.","family":"Main","sequence":"additional","affiliation":[],"role":[{"role":"author","vocabulary":"crossref"}]}],"member":"1965","published-online":{"date-parts":[[2023,3,23]]},"reference":[{"key":"B1","doi-asserted-by":"publisher","first-page":"9","DOI":"10.1080\/01621459.1988.10478560","article-title":"Statistical models for earthquake occurrences and residual analysis for point processes","volume":"83","author":"Ogata","year":"1988","journal-title":"J Am Stat Assoc"},{"key":"B2","doi-asserted-by":"publisher","first-page":"13","DOI":"10.1016\/j.tecto.2005.10.016","article-title":"Space-time ETAS models and an improved extension","volume":"413","author":"Ogata","year":"2006","journal-title":"Tectonophysics"},{"key":"B3","doi-asserted-by":"publisher","first-page":"217","DOI":"10.5047\/eps.2010.09.001","article-title":"Significant improvements of the space-time ETAS model for forecasting of accurate baseline seismicity","volume":"63","author":"Ogata","year":"2011","journal-title":"Earth Planets Space"},{"key":"B4","doi-asserted-by":"publisher","first-page":"83","DOI":"10.1093\/biomet\/58.1.83","article-title":"Spectra of some self-exciting and mutually exciting point processes","volume":"58","author":"Hawkes","year":"1971","journal-title":"Biometrika"},{"key":"B5","doi-asserted-by":"publisher","first-page":"1","DOI":"10.18637\/jss.v088.c01","article-title":"ETAS: an R package for fitting the space-time ETAS model to earthquake data","volume":"88","author":"Jalilian","year":"2019","journal-title":"J Stat Software Code Snippets"},{"key":"B6","doi-asserted-by":"publisher","first-page":"1473","DOI":"10.1785\/0120200198","article-title":"Bayesian estimation of the ETAS model for earthquake occurrences","volume":"111","author":"Ross","year":"2021","journal-title":"Bull Seismol Soc Am"},{"key":"B7","doi-asserted-by":"publisher","first-page":"395","DOI":"10.1146\/annurev-statistics-060116-054045","article-title":"Bayesian computing with INLA: a review","volume":"4","author":"Rue","year":"2017","journal-title":"Ann Rev Stat Appl"},{"key":"B8","doi-asserted-by":"publisher","first-page":"760","DOI":"10.1111\/2041-210X.13168","article-title":"inlabru: an R package for Bayesian spatial modelling from ecological survey data","volume":"10","author":"Bachl","year":"2019","journal-title":"Methods Ecol Evol"},{"key":"B9","doi-asserted-by":"publisher","first-page":"e2020JB020226","DOI":"10.1029\/2020JB020226","article-title":"Data-Driven optimization of seismicity models using diverse data sets: generation, evaluation, and ranking using Inlabru","volume":"125","author":"Bayliss","year":"2020","journal-title":"J Geophys Res Solid Earth"},{"key":"B10","doi-asserted-by":"publisher","first-page":"3231","DOI":"10.5194\/nhess-22-3231-2022","article-title":"Pseudo-prospective testing of 5-year earthquake forecasts for California using inlabru","volume":"22","author":"Bayliss","year":"2022","journal-title":"Natural Hazards Earth Syst Sci"},{"key":"B11","doi-asserted-by":"publisher","first-page":"e2644","DOI":"10.1002\/env.2644","article-title":"A joint Bayesian space-time model to integrate spatially misaligned air pollution data in R-INLA","volume":"31","author":"Forlani","year":"2020","journal-title":"Environmetrics"},{"key":"B12","doi-asserted-by":"publisher","first-page":"1145","DOI":"10.1177\/0962280216660421","article-title":"An intuitive Bayesian spatial model for disease mapping that accounts for scaling","volume":"25","author":"Riebler","year":"2016","journal-title":"Stat Methods Med Res"},{"key":"B13","doi-asserted-by":"publisher","first-page":"e0147172","DOI":"10.1371\/journal.pone.0147172","article-title":"Spatiotemporal evolution of Ebola virus disease at sub-national level during the 2014 West Africa epidemic: model scrutiny and data meagreness","volume":"11","author":"Santermans","year":"2016","journal-title":"PLoS ONE"},{"key":"B14","doi-asserted-by":"publisher","first-page":"241","DOI":"10.1007\/s00180-010-0208-2","article-title":"A primer on disease mapping and ecological regression using INLA","volume":"26","author":"Schr\u00f6dle","year":"2011","journal-title":"Comput Stat"},{"key":"B15","doi-asserted-by":"publisher","first-page":"725","DOI":"10.1002\/env.1065","article-title":"Spatio-temporal disease mapping using INLA","volume":"22","author":"Schr\u00f6dle","year":"2011","journal-title":"Environmetrics"},{"key":"B16","doi-asserted-by":"publisher","first-page":"1095","DOI":"10.1093\/jxb\/erv453","article-title":"Extensive tissue-specific transcriptomic plasticity in maize primary roots upon water deficit","volume":"67","author":"Opitz","year":"2016","journal-title":"J Exp Bot"},{"key":"B17","doi-asserted-by":"publisher","first-page":"2653","DOI":"10.1093\/eurheartj\/ehv216","article-title":"Road traffic noise is associated with increased cardiovascular morbidity and mortality and all-cause mortality in London","volume":"36","author":"Halonen","year":"2015","journal-title":"Eur Heart J"},{"key":"B18","doi-asserted-by":"publisher","first-page":"92","DOI":"10.1016\/j.marenvres.2015.08.005","article-title":"Modeling sensitive parrotfish (Labridae: Scarini) habitats along the Brazilian coast","volume":"110","author":"Roos","year":"2015","journal-title":"Mar Environ Res"},{"key":"B19","doi-asserted-by":"publisher","first-page":"12999","DOI":"10.1111\/zph.12999","article-title":"Bayesian spatiotemporal modelling analysis of hemorrhagic fever with renal syndrome outbreaks in China using R-INLA","volume":"2022","author":"Teng","year":"2022","journal-title":"Zoonoses Public Health"},{"key":"B20","doi-asserted-by":"publisher","first-page":"e1443","DOI":"10.1002\/wics.1443","article-title":"Spatial modeling with R-INLA: a review","volume":"10","author":"Bakka","year":"2018","journal-title":"Wiley Interdisc Rev Comput Stat"},{"key":"B21","doi-asserted-by":"publisher","first-page":"33","DOI":"10.1016\/j.sste.2012.12.001","article-title":"Spatial and spatio-temporal models with R-INLA","volume":"4","author":"Blangiardo","year":"2013","journal-title":"Spat Spatio Temporal Epidemiol"},{"key":"B22","doi-asserted-by":"crossref","DOI":"10.1201\/9781315175584","author":"G\u00f3mez-Rubio","year":"2020","journal-title":"Bayesian Inference With INLA"},{"key":"B23","doi-asserted-by":"publisher","first-page":"2266","DOI":"10.1080\/00949655.2013.788653","article-title":"INLA or MCMC? A tutorial and comparative evaluation for spatial prediction in log-Gaussian Cox processes","volume":"84","author":"Taylor","year":"2014","journal-title":"J Stat Comput Simulat"},{"key":"B24","doi-asserted-by":"publisher","DOI":"10.48550\/arXiv.2206.13360","article-title":"Approximation of bayesian Hawkes process models with Inlabru","author":"Serafini","year":"2022","journal-title":"arXiv preprint"},{"key":"B25","doi-asserted-by":"publisher","first-page":"438","DOI":"10.1111\/j.2517-6161.1971.tb01530.x","article-title":"Point spectra of some mutually exciting point processes","volume":"33","author":"Hawkes","year":"1971","journal-title":"J R Stat Soc B"},{"key":"B26","volume-title":"EnglishAdvanced Spatial Modeling with Stochastic Partial Differential Equations Using R and INLA","author":"Krainski","year":"2019"},{"key":"B27","doi-asserted-by":"publisher","first-page":"168501","DOI":"10.1103\/PhysRevLett.102.168501","article-title":"Origin and nonuniversality of the earthquake interevent time distribution","volume":"102","author":"Touati","year":"2009","journal-title":"Phys Rev Lett"},{"key":"B28","doi-asserted-by":"publisher","first-page":"90","DOI":"10.1785\/0120050067","article-title":"Comparison of short-term and time-independent earthquake forecast models for Southern California","volume":"96","author":"Helmstetter","year":"2006","journal-title":"Bull Seismol Soc Am"},{"key":"B29","doi-asserted-by":"publisher","first-page":"e2020JB021027","DOI":"10.1029\/2020JB021027","article-title":"B-Positive: a robust estimator of aftershock magnitude distribution in transiently incomplete catalogs","volume":"126","author":"van der Elst","year":"2021","journal-title":"J Geophy ResSolid Earth"},{"key":"B30","doi-asserted-by":"publisher","DOI":"10.1201\/9781351165761","article-title":"Bayesian Regression Modeling with INLA","author":"Wang","year":"2018","journal-title":"UK United Kingdom"},{"key":"B31","doi-asserted-by":"publisher","first-page":"885","DOI":"10.1785\/0120130270","article-title":"Statistical modeling of the 1997\u20131998 Colfiorito earthquake sequence: locating a stationary solution within parameter uncertainty","volume":"104","author":"Touati","year":"2014","journal-title":"Bull Seismol Soc Am"}],"container-title":["Frontiers in Applied Mathematics and Statistics"],"original-title":[],"link":[{"URL":"https:\/\/www.frontiersin.org\/articles\/10.3389\/fams.2023.1126759\/full","content-type":"unspecified","content-version":"vor","intended-application":"similarity-checking"}],"deposited":{"date-parts":[[2023,3,23]],"date-time":"2023-03-23T06:18:11Z","timestamp":1679552291000},"score":1,"resource":{"primary":{"URL":"https:\/\/www.frontiersin.org\/articles\/10.3389\/fams.2023.1126759\/full"}},"subtitle":[],"short-title":[],"issued":{"date-parts":[[2023,3,23]]},"references-count":31,"alternative-id":["10.3389\/fams.2023.1126759"],"URL":"https:\/\/doi.org\/10.3389\/fams.2023.1126759","relation":{},"ISSN":["2297-4687"],"issn-type":[{"value":"2297-4687","type":"electronic"}],"subject":[],"published":{"date-parts":[[2023,3,23]]},"article-number":"1126759"}}