The Need to Quantify Hazard Related to Non-magmatic Unrest: From BET_EF to BET_UNREST

Most volcanic hazard studies focus on magmatic eruptions and their accompanying phenomena. However, hazardous volcanic events can also occur during non-magmatic unrest, defined as a state of volcanic unrest in which no migration of magma is recognised. Examples include tectonic unrest, and hydrothermal unrest that may lead to phreatic eruptions. Recent events (e.g. Ontake eruption, September 2014) have demonstrated that the successful forecasting of phreatic eruptions is still very difficult. It is therefore of paramount importance to identify indicators that define the state of non-magmatic unrest. Often, this type of unrest is driven by fluids-on-the-move, requiring alternative monitoring setups, beyond the classical seismic-geodetic-geochemical architectures. Here we present a new version of the probabilistic model BET (Bayesian Event Tree), called BET_UNREST, specifically developed to include the forecasting of non-magmatic unrest and related hazards. The structure of BET_UNREST differs from the previous BET_EF (BET for Eruption Forecasting) by adding a dedicated branch to detail non-magmatic unrest outcomes. Probabilities are calculated at each node by merging prior models and past data with new incoming monitoring data, and the results can be updated any time new data has been collected. Monitoring data are weighted through pre-defined thresholds of anomaly, as in BET_EF. The BET_UNREST model is introduced here, together with its software implementation PyBetUnrest, with the aim of creating a user-friendly, open-access, and straightforward tool to support short-term volcanic forecasting (already available on the VHub platform). The BET_UNREST model and PyBetUnrest tool are tested through three case studies in the frame of the EU VUELCO project.

data with new incoming monitoring data, and the results can be updated any time new data has been collected. Monitoring data are weighted through pre-defined thresholds of anomaly, as in BET_EF. The BET_UNREST model is introduced here, together with its software implementation PyBetUnrest, with the aim of creating a user-friendly, open-access, and straightforward tool to support short-term volcanic forecasting (already available on the VHub platform). The BET_UNREST model and PyBetUnrest tool are tested through three case studies in the frame of the EU VUELCO project.

Introduction
Monitoring activities represent the main source of information to understand the behaviour of volcanic systems on short time-scales and, possibly, during emergency crises. In this framework, one of the main challenges of volcano monitoring is the identification and characterisation of the phase defined as "unrest", which consists of a relevant physical or chemical change in the volcanic system with respect to its background behaviour, leading to cause for concern. Unrest can be due to several factors and depends on the local characteristics of each volcanic system, making it very difficult to find general features or patterns (Phillipson et al. 2013). Unrest may be followed by volcanic eruptions due to the movement of magma, but can also be associated with other dangerous phenomena: indeed, in addition to magma-related hazards (e.g., tephra fallout, lava flows, ballistics), hydrothermal and tectonic activities, without evidence for "magma-on-the-move", can also lead to dangerous outcomes (i.e., flank collapses, gas emissions, phreatic explosions, lahars). Such hazardous events related to non-magmatic unrest are not easy to track and, in volcanic hazard evaluations, are sometimes underestimated . For instance, many volcanoes pass through a phase of hydrothermal unrest for years, decades or even centuries. Due to this long-term behavioural similarity, it is often difficult to recognise how hydrothermal unrest can lead to related hazards in the short-term. Where the driving agent and the main eruptive product is not magma, but water (liquid or vapour) and occasionally liquid sulphur, or gas, this type of unrest can lead to non-magmatic eruptions. On the other hand, non-eruptive hydrothermal unrest can also promote volcanic hazards after prolonged gas emissions, acidic fluid infiltration into aquifers, soils and the hydrologic network, or deformation induced by a rising fluid front (see . In this light, although most volcanic hazard assessments focus only on magmatic eruptions as potential hazard sources, hazardous events can also occur during non-magmatic unrest, which in this chapter is defined as a state of volcanic unrest in which no migration of magma is recognised. Examples of non-magmatic unrest include the tectonic (which causes concern independently on how it evolves and eventually ends), and hydrothermal unrest types; the latter may eventually lead to phreatic eruptions. Recent occurrences of phreatic eruptions (e.g. Ontake eruption, September 2014, Japan) have demonstrated that they are still very hard to anticipate from classical observations based on seismic-geodetic-geochemical monitoring architectures. For these reasons, it is of paramount importance to identify indicators that define the state of non-magmatic unrest. Often, this type of unrest is driven by "fluids-on-the-move", requiring alternative and innovative monitoring setups, beyond the classical ones.
In the last decade it has become crucial to provide forecasts of the possible outcomes of volcanic unrest, to give quantitative support and scientific advice to decision makers (e.g., Woo 2008;Woo 2007, 2009). Because of this, event tree schemes have been proposed (e.g., Newhall and Hoblitt 2002;Marzocchi et al. 2004), and a few probabilistic tools based on event trees and Bayesian inference have been developed (e.g., BET_EF, Marzocchi et al. 2008;HASSET, Sobradelo et al. 2013) with the ability to quantify the probability of different possible outcomes related to magmatic unrest. However, the need for recognising and tracking the evolution of any type of volcanic unrest, and to quantify the probability linked to non-magmatic unrest as well, have led us, within the VUELCO project, to the development of a new probabilistic model, able to forecast both magmatic and non-magmatic hazardous events related to volcanic unrest: BET_UNREST. The BET_UNR-EST model is based on an event tree, whose structure is extended with respect to the previous schemes such as BET_EF (see the generalisation from BET_EF to BET_UNREST in Fig. 1, highlighted in red) by adding a specific branch to detail the track and outcome of non-magmatic unrest. Nonetheless, BET_UNREST adopts from BET_EF the Bayesian inferential paradigm and the ability to account both for long-term data (typically from the geological record) and short-term information from monitoring networks.
In this chapter, we briefly present the BET_UNREST model and its implementation in the PyBetUnrest software tool , made with the aim of providing a user-friendly, open-access, and straightforward tool to handle probabilistic forecasts and visualise results, and that has already been included on the Vhub platform (https://vhub.org/ resources/betunrest). The new event tree and tool are applied here as illustrative examples to the VUELCO target volcanoes Popocatépetl (Mexico), Cotopaxi (Ecuador) and Dominica (West Indies).

BET_UNREST Model and PyBetUnrest Tool
As with all the previous BET models (e.g., BET_EF, for short-and long-term eruption forecasting, Marzocchi et al. 2008; BET_VH, for long-term volcanic hazard associated to any potential hazardous phenomenon accompanying an eruption, Marzocchi et al. 2010;Tonini et al. 2015; BET_VHst, a model that merges the previous two, Selva et al. 2014), BET_UNREST performs probabilistic assessments in the frame of volcanic hazard analysis, based on an event tree scheme. The main novelty in the BET_UNREST event tree is the introduction, with respect to the BET_EF tree, of a new branch ( Fig. 1) for exploring and forecasting the outcomes of non-magmatic unrest ). Due to the resemblance of BET_UNR-EST to other BET models from a methodological and computational point of view, here we will only give a brief overview. The papers by Marzocchi et al. (2004Marzocchi et al. ( , 2008 provide a more detailed description. BET_UNREST probabilities are evaluated by a Bayesian inferential procedure, in order to quantify both the aleatory and epistemic uncertainty characterising the impact of volcanic eruptions in terms of eruption forecasting and/or hazard assessment. Such a procedure allows merging all the available information, such as models, a priori beliefs, past data from volcanic records and, when available, real-time monitoring data in order to include, in principle, all the knowledge about the considered volcanic system.
In general, the Bayesian inference procedure at the basis of BET_UNREST assigns a probability to each node, providing a framework where: probabilities are expressed through a probability density function (pdf), and not as a single number, to account for a best-evaluation value (for example the mean of the probability density function, representing a degree of aleatory uncertainty) and for a measure of the epistemic uncertainty (the dispersion of the pdf); the posterior pdf, at each node, is achieved by statistically combining, through Bayes' theorem, a prior probability distribution (usually coming from theoretical models and/or expert judgement) and information from the available data relevant for that node.
As in BET_EF, the probability ½h k at each node k is actually described by a statistical mixing of two pdfs, describing respectively the "so-called" long-term ½h f Mg k and short-term ½h fMg k regimes of the volcano as follows: where c k represents the weight in the interval [0,1] depending on the degree of unrest (Marzocchi et al. 2008). With such mixing, BET_UNREST switches between the two "regimes". In practice: • When anomalies with respect to the volcano's background activity are not observed at time t = t 0 , BET_UNREST relies on the so-called long-term information to assign the probabilities (hereinafter also referred to as background probabilities) at the various branches.
Such background probabilities (i.e., ½h f Mg k ) are based on theoretical models and information from the geological and eruptive record of the volcano studied, or of similar volcanoes, and describe the long-term frequencies of magmatic or non-magmatic unrest, and subsequent outcomes at these volcanic systems.
• When a clear state of unrest of whatever nature is detected at t = t 0 by BET_UNREST, the probabilistic assignment at all the successive nodes is based mainly on the monitoring information. In practice, monitoring data are transformed into subjective pdfs (i.e., ½h fMg k ) relative to the occurrence of magmatic or non-magmatic unrest and the following branches. Actually, at some nodes, monitoring data are not considered as relevant (for example, in forecasting the size of an eruption, magmatic or not), and here BET_UNREST continues to rely on theoretical models and long-term frequencies.
• When, at time t = t 0 , BET_UNREST observes a "degree of unrest" (of whatever nature) without it being completely clear, the statistical mixing provides a resulting pdf which accounts for both the regimes, giving the short-term regime a weight equal to the degree of unrest, and to the long-term regime its complement.
In this way, during a phase of unrest, the past data have less (null, in the case of complete unrest) importance. The short term hazard/ eruption forecasting depends exclusively on the translation of observed anomalies into pdfs describing all the branches of the event tree. This is done, separately at each node, by weighting monitoring data through pre-defined thresholds of anomaly (Marzocchi et al. 2008) and converting the resulting "degree of anomaly" into a best-evaluation probability, to which a degree of variance is associated (Fig. 2). This is a very simple and intuitive procedure, in which the basic assumptions are: 1. the first anomaly detected is the most informative 2. subsequent anomalies contribute less and less to the increase of the degree of anomaly 3. strong non-linear coupling among anomalies are neglected.
At each node, BET_UNREST evaluates the following probabilities (see also Fig. 1) by means of Bayesian inference (we give the acronyms used throughout the chapter to indicate the probability at each node in brackets): • Unrest: probability (P(U)) of unrest in the time period [t 0 ; t 0 + |], given the monitoring observations at time t = t 0 ; the time window | is defined by the user; • Magmatic unrest: probability (P(MU)) that the unrest is due to "magma-on-the-move", given the unrest; • Magmatic eruption: probability (P(MEr)) of a magmatic eruption, given magmatic unrest; the following sub-branches mirror the BET_EF structure, so we point the reader to Marzocchi et al. (Marzocchi et al. 2008) for them; • Non-magmatic unrest: this is the complementary of the Magmatic unrest branch, so by definition is the probability of non-magmatic unrest, given an unrest; • Hydrothermal unrest: probability (P(HU)) of hydrothermal unrest, given a non-magmatic unrest; • Tectonic unrest: this is the complementary of the Hydrothermal unrest branch, so it describes the probability (P(TU)) of a tectonic unrest, given a non-magmatic unrest; • Hydrothermal eruption: probability (P(HEr)) of a hydrothermal eruption, given a hydrothermal unrest; • Vent of hydrothermal eruption: here we explore the spatial probability of vent opening in a hydrothermal eruption, given a hydrothermal eruption occurring; this node is an extension with respect to the event tree proposed in ; • Size of hydrothermal eruption: probability of an explosive hydrothermal eruption, given a hydrothermal eruption occurring from a specific vent; its complementary branch is the effusive hydrothermal eruption.
In order to keep the structure of BET_UNR-EST as simple as possible, an effort has been made to maintain, where possible, a dichotomic branching into complementary (i.e., exhaustive and mutually exclusive) events. This is why the Unrest node does not branch directly into magmatic, hydrothermal and tectonic, but first it branches into magmatic-or-not. This allows a simplification in the evaluation of short-term probabilities. In particular, with this type of ramification, the user defines which monitoring measurements (plus thresholds and weight) affect the pdf of one of the two branches; the pdf of the complementary branch then comes automatically.
The new BET_UNREST model is applied here with its software implementation PyBe-tUnrest presented in Tonini et al. (2016), which aims to provide an open and usable tool to bridge between the scientific community and decision makers, with a graphical user interface which allows the exploration of the event tree and the visualisation of the results (see Fig. 1). This solution was also implemented in the VHub cyber-infrastructure (http://vhub.org/resources/ betunrest). In the present PyBetUnrest tool only one file needs to be adapted when new monitoring information is gathered. This structure makes PyBetUnrest extremely fast and user-friendly during crisis situations. More on the technical background of the BET_UNREST model and PyBetUnrest tool can be found in the VUELCO Deliverable 7.3 (at http://vhub.org) and in Tonini et al. (2016).
So far BET_UNREST and PyBetUnrest have not yet been blindly tested in real-time during an actual volcanic crisis, but only retrospectively  at Kawah Ijen (Indonesia), for the time period 2010-2012 (after a learning period based on the observations from 2000 to 2010). The term "blindly" signifies that the rules of BET_UNREST (the long-term pdfs, and the monitoring parameters, thresholds and weight at the different nodes) are set before the beginning of the application, on different data (the learning dataset), and then the model is applied untouched to new data (the voting dataset), typically covering a different time period (as in the case of Tonini et al. 2016). transformed into a best-evaluation probability at a given node of the event tree. First, a monitoring measure x i is translated in a degree of anomaly z i according to a selected anomaly function l(Á) (a). In the above example, a measure below x 1 is considered background, above x 2 is anomalous, and in between it has a certain degree of anomaly. After collecting the degree of anomaly for all parameters considered at the node, we combine them using a weighted average (xi is the weight of the i-th parameter) in order to obtain the total degree of anomaly (b). Then the total degree of anomaly is transformed into an average probability using a predefined function, in BET_UNR-EST, we use the function in (c). The parameters, weights, and thresholds are selected by the user, possibly through expert opinions' elicitation. Figure modified from (Marzocchi and Bebbington 2012) In the next section of this chapter, results and performances of the new model and tool will be discussed and validated by analysing the unrest crises for VUELCO target volcanoes Popocatépetl, Cotopaxi and Dominica through blind applications of BET_UNREST. The latter two applications show the results of the VUELCO crisis simulation exercises held in Quito (November 2014) and Dominica (May 2015). The time window, |, is defined as 1 month.

BET_UNREST Applications
In Table 1 we report the activity carried out 24/7 with regards to monitoring at Popocatépetl, available as short-term information for unrest, origin of unrest and eruption. However, for the time period 1994-2012, the available data (as listed in Mendoza-Rosas, VUELCO deliverable 5.1), are restricted mainly to seismicity (VT, tremor, number of events) and visual observations (i.e. number of eruptions, column height). No real-time SO 2 flux is available for our purpose, and deformation data would need further processing. Regarding past data (long-term information for unrest, origin of unrest, and eruption), there have been 13 unrest episodes, and constant unrest since December 1994 (so, a priori probability to be in a state of unrest for the next month is about 85%). Out of the 13 unrest episodes, 6 were due to magma-on-the-move (magmatic unrest), of which 3 lead into a magmatic major eruption. The monitoring parameters listed in Table 2, along with respective thresholds and weight, have been identified in the UNAM (Universidad Nacional Autónoma de México) database for the period 1997-2012, and used to set BET_UNREST for Popocatépetl. The volcano is a stratocone with a higher probability of an eruption to occur from the central vent. For the period of observation (1997)(1998)(1999)(2000)(2001)(2002)(2003)(2004)(2005)(2006)(2007)(2008)(2009)(2010)(2011)(2012) all eruptions were magmatic and occurred at the central crater. The a priori spatial distribution of vent opening is assigned as in Table 3. As a prior model to define the size/style of magmatic eruptions we take the power law from Simkin and Siebert (1994). As past data we take the Mendoza-Rosas and De la Cruz-Reyna (2008) catalog for the past 23,000 years, and assume it to be complete for VEI >= 2 (Table 3).
We retrospectively applied BET_UNREST for the voting period April-June 2013, in which respectively 10, 11 and 2 eruptions of 2, 3 and 4 km-high columns were observed. No major eruption occurred. Observed anomalies include ash eruptions up to 130/day (all with columns <4 km), seismic tremor, incandescence in the crater/dome, and VT events (but no shallow event with depth <5 km). There was no anomalous deformation, no dome growth, and no SO 2 data available. Results of P(MEr) for the retrospective application period (weekly updated) are presented in Fig. 3. For the whole period, P (MEr) of a major eruption (>8 km eruption column) was <1% per month.   Magmatic eruption -Dome growth Y/N -SO 2 >9000 t/d -Tectonic EQ > M5.5 (along the coast/arc Michoacán-Chiapas) Y/N -Incandescent debris Y/N -Change in # tremor Y/N -VT depth <5 km -Increase # VT >M2 Y/N -Duration tremor >30,000 s (inertia 2 months) Y/N -Increase # ash eruptions >2000-4000 (inertia 2 months) -Deformation Y/N rhyolitic composition, is one of the most active and hazardous volcanoes in Ecuador. Historic eruptions at Cotopaxi produced large lithic-rich pyroclastic flows, ash flows, lava flows as well as large lahars (Barberi et al. 1995;Hall and Mothes 2008;Biass and Bonadonna 2011). Some lahars reached the Pacific Ocean at >200 km distance (Aguilera et al. 2004;Pistolesi et al. 2013). Recent unrest periods at Cotopaxi occurred in 1975-1976and 2001-2002 and were characterised by increased fumarolic activity, elevated seismicity and edifice deformation (Molina et al. 2008). Fumarolic activity is a concern due to the heat transfer that may affect the ice cover resulting in non-eruptive debris flows or lahars.

Cotopaxi, Ecuador: Retrospective Application Inspired by the VUELCO Simulation Exercise in Quito
A still unstable version of PyBetUnrest was set up (along with parameters and thresholds at each node for Cotopaxi volcano derived from monitoring information) before the simulation exercise, based on the available data in the literature up to the beginning of the simulation (the learning period stopped with the beginning of the exercise), in order to preliminarily test its value in decision support by providing near-real time probabilities of (i) the occurrence of unrest, (ii) the origin and nature of unrest and (iii) eruptive activity. However, during the simulation, the reports from the "volcano team" did not reflect the real eruptive and unrest history of Cotopaxi, as the past activity for the simulation was  West flank 0.0025 VEI ! 5 0.0008 2 "invented". A different setting of BET_UNREST (and consequently of PyBetUnrest) on site was not possible due to the lack of time and the still premature customisability of the tool. This obliged us to set up and run the old BET_EF tool during the exercise (Constantinescu et al. 2015).
Obviously, this prevented us from providing probabilistic assessment of non-magmatic events during the exercise at Cotopaxi: this would have been possible with BET_UNREST, enabling the calculation of probabilities for hydrothermal unrest and hydrothermal eruptions (P(HU) and P (HEr)). Nevertheless, the unrest scenario proposed by the "volcano team" (Bulletins 1-5) did not emphasise a significant state of hydrothermal unrest, which, on the one hand, made our output less biased in not providing an evaluation for P (HU) and P(HEr); but on the other hand this simulation was probably not the best case to test BET_UNREST.
Here, we will re-run BET_UNREST and PyBetUnrest at Cotopaxi retrospectively for the unrest phases described in the five bulletins provided by the "volcano team" during the simulation exercise and using the BET_UNREST setup prepared prior to the simulation based on the real past activity of the volcano (Table 4). The time window | was set to 1 month. In Table 5 we show the probabilities resulting from the run of the code, after each bulletin: (1) Phase 0: The background activity of Cotopaxi (NO anomalies): results are based on the past activity of Cotopaxi, with all observation within background limits. (2) Phase 1 (Bulletin 1): the observed anomalies in this phase were limited to an increase in seismic activity compared to background level. Such an increase is indicative, according to pre-set parameters, of magma-on-themove (P(MU) = 0.68). The considerable uncertainty is summarised by the 10th to 90th percentiles confidence interval. (3) Phase 2 (Bulletin 2): the observed anomalies in this phase were: a drastic increase in seismicity, an increase in SO 2 emission (5 times background levels), and a crater thermal anomaly. As a consequence, the mean P (MU) increases, along with a decrease in the associated uncertainty. (4) Phase 3 (Bulletin 3): the observed anomalies in this phase were: an increase in VT and LP events, occurrence of tremor, appearance of new fumaroles, an increase in SO 2 emission, and an increase in the crater thermal anomaly. As a consequence, the P(MU) is similar to Bulletin 2, but the P(HU) increases slightly, due to the new fumaroles. (5) Phases 4 and 5 (Bulletins 4 and 5): the observed anomalies in these phases were similar, and included: intense fumarolic activity, occurrence of hybrid seismic events, an increase in SO 2 emission, and an increase in the crater thermal anomaly. As a consequence, P(MEr) increases from 0.21 (phase 3) to 0.57, combined with a lower uncertainty.

Dominica, West Indies, Lesser Antilles: VUELCO Simulation Exercise, Dominica, May 2015
Dominica is characterised by hydrothermal activity manifested as thermal springs (up to boiling temperature), boiling-temperature fumarolic emissions (e.g. Valley of Desolation) and a crater lake, known as 'Boiling Lake', with a particular hydrodynamic behaviour (Fournier et al. 2009;Joseph et al. 2011;Rouwet et al. 2017). No high-temperature manifestations occur on the island, so no clear evidence of active magmatic degassing exists at the present time.
The simulation exercise, and consequently the BET_UNREST application, for the VUELCO target island of Dominica mainly focused on an unrest scenario for the southern part of the island. The purpose of the exercise was to test the tracking/assessment of an unrest period, and the decision making process undertaken by the scientific advisory group and local authorities.
Due to the hydrothermal character of Dominica, the application of BET_UNREST is highly suited. Before the simulation exercise, the PyBetUnrest tool was set for Dominica, based on (1) existing literature of the past volcanic activity; (2) insights on the current hydrothermal activity; (3) discussion-based expert elicitation sessions (4 sessions at SRC and 1 at INGV-Bologna); and (4) exchanges with local experts in order to fine-tune the code with the monitoring parameters. We remark that all of this was done prior to the start of the simulation exercise (the learning period stopped at the beginning of the simulation exercise, as for Cotopaxi), and again no hindsight tuning was made. The long-term setup of PyBetUnrest is done by filling up a configuration file that includes the a priori and past data specifically for Dominica, whose main information is summarised in Table 6. The short-term information is listed in Table 7 (parameters and thresholds identified prior to the exercise onset, see above). Further details on the Dominica simulation exercise and on the BET_UNREST application are given in Constantinescu et al. (2016).
During the simulation exercise (May 14-15, 2015) three phases of changes in volcanic activity, each with a duration of six months, were distributed by the "volcano team" to the operators of the unrest crisis. The reports included four types of observations: (1) seismic bulletin, (2) GPS, (3) geothermal monitoring data, and (4) other observations.
The translation of the reported bulletins into the values for the selected parameters in the BET_UNREST for Dominica setup were reported back to the team of experts in real-time  during the simulation. In Table 8 we provide the probabilities resulting from the run of the code after each bulletin. In Fig. 4 we also provide the time evolution of some of the most relevant probability distributions, across all the time periods spanned by the simulation exercise in Dominica. For each bulletin, among the output information from PyBetUnrest, there were two Some of the data are too many to be listed (this is indicated by the label "file" in the table). They can be provided in the form of files on request maps of the spatial probability of vent opening: one for the case of magmatic eruption, and one for hydrothermal eruption (Fig. 4). We believe this could be particularly useful, for example in a volcanic system like Dominica, where there are numerous areas showing hydrothermal activity, thus increasing the uncertainty on the position of a possible phreatic event.
The parameter "detectable SO 2 , HCl, HF" created confusion and opened up a scientific discussion. For the sake of transparency, we provide the mean values of P(MU) and P(MEr) including, or not, the HCl anomaly (Table 8).
Beyond the scientific implications of this issue, this concern reflected the sensitivity of BET_UNREST to the interpretation of some parameters. When relatively few monitoring parameters are provided, the weight of a single anomaly can be high: this is somehow a measure of the epistemic uncertainty.

Discussion and Implications for Unrest Tracking
This chapter presents the need for an updated BET model and tool that is able to account for the non-magmatic nature of some volcanic unrest episodes, which can often go under-estimated, if not totally neglected. The new model (BET_UNREST) and tool (PyBetUnrest) allow the tracking of unrest phases at volcanic systems and enables short-term volcanic forecasts. It has been fully developed within the VUELCO project, during which time it has been applied to some of the project's target volcanoes. In general, when we are able to distinguish magma-on-the-move ) from the monitoring observations the new model basically "collapses" to BET_EF (or, better, the assessment of the probabilities related to  In bracket estimates of mean values without including HCl anomaly in Phase 3 magmatic outcomes provided by the two models coincide). On the other hand, if we are not able to identify a magmatic "active role" in the unrest (from the available monitoring observations), BET_UNREST is still able to provide the probabilities of hazardous events that accompany non-magmatic volcanic unrest, rather than neglecting them. As discussed in , a very difficult case is presented by phreatomagmatic eruptions that, sometimes, can occur without any precursors indicating magma movement. This is surely an important limit to overcome which requires further efforts to detect subtle changes in the very short-term (hours to minutes) by improving monitoring techniques.
The chapter illustrates the development and implementation of BET_UNREST model and PyBetUnrest tool through three different applications: (i) the pure retrospective analysis at Popocatépetl volcano, where there is no compelling need for a hydrothermal branch due to the current magmatic nature of the unrest episodes. Popocatépetl has remained in unrest from December 1994 to present and, for this application, BET_UNREST and PyBetUnrest were run using the UNAM Data Base for the learning period 1997-2012, with a On bottom left, a snapshot of PyBetUnrest tool shows the spatial probability of vent opening during Phase 1, localising the most probable position of the phreatic eruption retrospective application aiming to forecast major eruptions (column heights greater than 8 km) for the April-June 2013 volcanic activity. (ii) the application based on a simulation exercise at Cotopaxi. Here we tested the BET_UNREST retrospectively, but, this time, using the invented data provided during the VUELCO simulation exercise, in addition to data based on the real past history of the volcano. (iii) the almost real-time simulation exercise organised by the VUELCO project in Dominica (May 2015). The volcanic system of Dominica presents a "prototype" setting for BET_UNREST due to its hydrothermal character. Phreatic/ phreatomagmatic activity occurred during the simulation, coinciding with high associated probabilities from BET_UNREST (the average values P(HU) = 0.73 and P(HEr) = 0.32). We also positively tested the feasibility of providing different maps of the spatial probability of vent opening in case of magmatic or phreatic eruption.
As mentioned in previous sections, we implemented the BET_UNREST model into PyBetUnrest software tool using a graphical user interface aiming to provide a fast, open and user-friendly tool, which extends the usage of BET_UNREST to volcanologists with different expertise. The PyBetUnrest tool reached a mature and usable version during the Dominica simulation and its first stable release has been uploaded to Vhub cyber-infrastructure.
With these exercises we strongly believe we have brought BET a step closer to a full and proper implementation during a crisis situation. The PyBetUnrest tool eventually worked as expected, but it is important to take advantage of the lessons learned during these applications and pursue more tests that will improve its design and prove its usefulness in real-case scenarios.
As a final comment, we would like to remark that, as with any other event tree model (e.g. BET models by Marzocchi et al. 2004Marzocchi et al. , 2008Marzocchi et al. , 2010; HASSET model by Sobradelo et al. 2013), one can always apply and "populate" the BET_UNREST model in any "volcanic" circumstance. The uncertainty on the results provided by BET_UNREST, and consequently their practical use, will however be strongly dependent on the available information and data used to set up the models rules. If only a few pieces of evidence are available, the models results will be characterised by a large uncertainty, and thus might be not very helpful for decision-makers. As more and more knowledge is gathered, BET_UNREST output probabilities will become more attractive from a practical point of view, since their uncertainty will be increasingly small. This is an intrinsic feature of the Bayesian inferential procedure at the basis of the model.