loading page

A statistical approach for spatial mapping and temporal forecasts of volcanic eruptions using monitoring data
  • +8
  • Abani Patra,
  • Andrea Bevilacqua,
  • Eric Bruce Pitman,
  • Marcus Bursik,
  • Barry Voight,
  • Augusto Neri,
  • Giovanni Macedonio,
  • Franco Flandoli,
  • Prospero De Martino,
  • Flora Giudicepietro,
  • Stefano Vitale
Abani Patra
Tufts University
Author Profile
Andrea Bevilacqua
Istituto Nazionale di Geofisica e Vulcanologia

Corresponding Author:andrea.bevilacqua@ingv.it

Author Profile
Eric Bruce Pitman
University at Buffalo
Author Profile
Marcus Bursik
University at Buffalo
Author Profile
Barry Voight
Pennsylvania State University
Author Profile
Augusto Neri
Istituto Nazionale di Geofisica e Vulcanologia
Author Profile
Giovanni Macedonio
Istituto Nazionale di Geofisica e Vulcanologia
Author Profile
Franco Flandoli
Scuola Normale Superiore di Pisa
Author Profile
Prospero De Martino
Istituto Nazionale di Geofisica e Vulcanologia
Author Profile
Flora Giudicepietro
Istituto Nazionale di Geofisica e Vulcanologia
Author Profile
Stefano Vitale
Università di Napoli Federico II
Author Profile

Abstract

We present two models using monitoring data in the production of volcanic eruption forecasts. The first model enhances the well-established failure forecast method introducing an SDE in its formulation. In particular, we developed new method for performing short-term eruption timing probability forecasts, when the eruption onset is well represented by a model of a significant rupture of materials. The method enhances the well-known failure forecast method equation. We allow random excursions from the classical solutions. This provides probabilistic forecasts instead of deterministic predictions, giving the user critical insight into a range of failure or eruption dates. Using the new method, we describe an assessment of failure time on present-day unrest signals at Campi Flegrei caldera (Italy) using either seismic count and ground deformation data. The new formulation enables the estimation on decade-long time windows of data, locally including the effects of variable dynamics. The second model establishes a simple method to update prior vent opening spatial maps. The prior reproduces the two-dimensional distribution of past vent distribution with a Gaussian Field. The likelihood relies on a one-dimensional variable characterizing the chance of material failure locally, based, for instance, on the horizontal ground deformation. In other terms, we introduce a new framework for performing short-term eruption spatial forecasts by assimilating monitoring signals into a prior (“background”) vent opening map. To describe the new approach, first we summarize the uncertainty affecting a vent opening map pdf of Campi Flegrei by defining an appropriate Gaussian random field that replicates it. Then we define a new interpolation method based on multiple points of central symmetry, and we apply it on discrete GPS data. Finally, we describe an application of the Bayes’ theorem that combines the prior vent opening map and the data-based likelihood product-wise. We provide examples based on either seismic count and interpolated ground deformation data collected in the Campi Flegrei volcanic area.