Operational Aftershock Forecasting for 2017-2018 Seismic Sequence in Western Iran

On Sunday November 12, 2017, at 18:18:16 UTC, (21:48:16 local time), a strong earthquake with Mw7.3 occurred in western Iran in the border region between Iran and Iraq in vicinity of the Sarpol-e Zahab town. Unfortunately, this catastrophic seismic event caused 572 causalities, thousands of injured and vast amounts of damage to the buildings, houses and infrastructures in the epicentral area. The mainshock of this seismic sequence was felt in the entire western and central provinces of Iran and surrounding areas. The main event was preceded by a foreshock with magnitude 4.5 about 43 minutes before the mainshock that warned the local residence to leave their home and possibly reduced the number of human casualties. More than 2500 aftershocks with magnitude greater than 2.5 have been reported up to January 2019 with the largest registered aftershock of Mw6.4. A novel and fully-probabilistic procedure is adopted for providing spatio-temporal predictions of aftershock occurrence in a prescribed forecasting time interval (in the order of hours or days). The procedure aims at exploiting the information provided by the ongoing seismic sequence in quasi-real time. The versatility of the Bayesian inference is exploited to adaptively update the forecasts based on the incoming information as it becomes available. The aftershock clustering in space and time is modelled based on an Epidemic Type Aftershock Sequence (ETAS). One of the main novelties of the proposed procedure is that it considers the uncertainties in the aftershock occurrence model and its model parameters. This is done by moving within a framework of robust reliability assessment which enables the treatment of uncertainties in an integrated manner. Pairing up the Bayesian robust reliability framework and the suitable simulation schemes (Markov Chain Monte Carlo Simulation) provides the possibility of performing the whole forecasting procedure with minimum (or no) need of human interference. The fully simulation-based procedure is examined for both Bayesian model updating of ETAS spatio-temporal model and robust operational forecasting of the number of events of interest expected to happen in various time intervals after main events within the sequence. The seismicity is predicted within a confidence interval from the mean estimate.