Probabilistic modeling of future volcanic eruptions at Mount Etna

Transcription

1 JOURNAL OF GEOPHYSICAL RESEARCH: SOLID EARTH, VOL. 118, , doi: /jgrb.50190, 2013 Probabilistic modeling of future volcanic eruptions at Mount Etna Annalisa Cappello, 1 Giuseppe Bilotta, 1,2 Marco Neri, 1 and Ciro Del Negro 1 Received 4 December 2012; revised 20 February 2013; accepted 11 April 2013; published 8 May [1] The statistical analysis of volcanic activity at Mt Etna was conducted with the twofold aim of (1) constructing a probability map for vent opening of future flank eruptions and (2) forecasting the expected number of eruptive events at the summit craters. The spatiotemporal map of new vent opening at Etna volcano is based on the analysis of spatial locations and frequency of flank eruptions starting from Thanks to the completeness and accuracy of historical data over the last four centuries, we examined in detail the spatial and temporal distribution of flank eruptions showing that effusive events follow a nonhomogenous Poisson process with space-time varying intensities. After demonstrating the spatial nonhomogeneity and the temporal nonstationarity of flank eruptions at Etna, we calculated the recurrence rates (events expected per unit area per unit time) and produced different spatiotemporal probability maps of new vent opening in the next 1, 10 and 50 years. These probabilistic maps have an immediate use in evaluating the future timing and areas of Etna prone to volcanic hazards. Finally, the results of the analysis of the persistent summit activity during the last 110 years indicate that the hazard rate for eruptive events is not constant with time, differs for each summit crater of Mt Etna, highlighting a general increase in the eruptive frequency starting from the middle of last century and particularly from 1971, when the SE crater was formed. Citation: Cappello, A., G. Bilotta, M. Neri, and C. Del Negro (2013), Probabilistic modeling of future volcanic eruptions at Mount Etna, J. Geophys. Res. Solid Earth, 118, , doi: /jgrb Introduction [2] Mt Etna (Italy) is one of the most active and hazardous volcanoes in the world [Behncke et al., 2005]. It is well known for the persistent activity from the summit craters, frequent lava flow-forming eruptions from vents situated on the flanks of the volcano, and the large population settled on and around the sides of the volcano that is at risk [Guest and Murray, 1979; Bisson et al., 2009]. The ever-expanding use of areas near the volcano increases the potential impact of future eruptions of Mt Etna on the regional economy and on the health and safety of its citizens (Figure 1). Eruptions can neither be prevented nor interrupted, but actions can be taken to minimize damage from them [Ganci et al., 2012b]. Reduction of risk to life and property can be undertaken by avoiding threatened areas and by taking protective measures to reduce the effects when and where vulnerable areas cannot be avoided [Scifoni et al., 2010]. Etna is continuously monitored using seismology, deformation, gravity, magnetism, gas emission studies, and petrology by Istituto Nazionale di Geofisica e Vulcanologia (INGV) and its observatory in Catania [Napoli 1 Istituto Nazionale di Geofisica e Vulcanologia, Sezione di Catania, Osservatorio Etneo, Catania, Italy. 2 Dipartimento di Matematica e Informatica, Università di Catania, Catania, Italy. Corresponding author: C. Del Negro, Istituto Nazionale di Geofisica e Vulcanologia, Sezione di Catania, Osservatorio Etneo, Italy. (ciro.delnegro@ct.ingv.it) American Geophysical Union. All Rights Reserved /13/ /jgrb et al., 2008; Greco et al., 2010; Bonaccorso et al., 2011]. Detection of volcanic precursors can generally identify the locality of impending volcanic activity, even though it often does not determine exactly the typology or timing of an eruption or even its certainty. Hazard-zonation maps can then be used for risk-based decision making in land-use planning and emergency management, as well as addressing the more general problem of quantitative volcanic hazards assessment [Wadge et al., 1994; Newhall, 2000; Sparks, 2003; Magill et al., 2005; Behncke et al., 2005; Favalli et al., 2009; Bisson et al., 2009; Crisci et al., 2010; Cappello et al., 2011a; Vicari et al., 2011a, 2011b]. Thus, effective monitoring of Etna volcano, combined with the preparation of emergency plans to face future eruptions, can help reduce the risk to lives and property in an area densely urbanized, where about 0.8 million people live [Behncke et al., 2005]. [3] Our purpose is to give a proper statistical treatment of the records of volcanic eruptions at Mt Etna to allow key atrisk areas to be rapidly and appropriately identified. In terms of hazard, persistent summit activity does not represent a significant threat to the towns located on Etna s slopes, although important tourist facilities and infrastructures close to the central craters have been repeatedly destroyed. However, periodic flank eruptions pose a serious danger to the lives and the property of the local populace [Behncke et al., 2005], especially when the eruptive vents open up at very low altitudes (up to m above sea level). In terms of volcanological data, the catalog of Etna flank eruptions is fully reliable only for the last ~400 years, while summit eruptions are sufficiently documented starting from the first decades of the 20th century [Mulargia et al., 1985; Behncke et al., 2005; Tanguy et al., 1925

2 CAPPELLO ET AL.: PROBABILISTIC HAZARD MAPPING Figure 1. Map of Etna, showing the spatial distribution of eruptive fissures (in dark gray) produced by the flank eruptions occurring from 1610 to The dashed black line delineates the area of 1170 km2 containing the 500 m spaced grid of potential vents. (a) Rose diagram representing the geometric orientation of the eruptive fissures; (b) frequency of the eruptive fissures during the last four centuries (the last column on right cumulates the fissures opened in the period ). 2007]. Therefore, the statistical analysis of volcanic activity of Mt Etna was conducted with the dual aim of (1) constructing a spatiotemporal probability map for vent opening of future flank eruptions and (2) forecasting the expected number of eruptive events at the summit craters. [4] Different approaches have been used to evaluate the probability of vent opening at Etna volcano by analyzing the distribution of spatial location of flank eruptions and their temporal frequency. [5] For the spatial distribution, the density of vents has been calculated to identify areas where the highest concentration is located, which correspond to the zones known as NE, S, and W Rifts [Guest and Murray, 1979; Acocella and Neri, 2003; Salvi et al., 2006; Crisci et al., 2010; Rongo et al., 2011; Neri et al., 2011; Cappello et al., 2012]. [6] For the temporal distribution of flank eruptions, Mulargia et al. [1985] concluded that Etna s lateral eruptions follow a stationary Poisson process with a constant intensity. Salvi et al. [2006] refuted this result, proving that flank eruptions at Etna do not follow a simple Poisson distribution but are more likely represented by a nonhomogeneous Poisson process with power law intensity. Recently, Bebbington [2007] applied the Hidden Markov Models (HMM) to volcanic occurrences to identify the activity state (hidden) of the volcano most consistent with the observations and to forecast the next eruptions. However, the stationary distribution assumed by HMM cannot properly model the increasing trend of flank eruptions observed over the last 40 years at Etna [Allard et al., 2006]. Finally, the temporal analysis performed by Smethurst et al. [2009] revealed that flank eruptions follow a nonhomogeneous Poisson process with a piecewise intensity increasing nearly linearly since the mid-1900s. [7] In this paper, we carefully assess the distribution of future vent locations using a probabilistic modeling that combines both the temporal and spatial analysis of flank eruptions in the past four centuries. Spatiotemporal probability maps for future vent opening on the flanks of Mt Etna in the next few years/decades are presented and discussed. We assume that the temporal patterns extracted from the records of past eruptions are useful to predict the future behavior of a volcano; likewise, we consider the regions with a high density of volcanic structures to be the most probable emission zones of future lava flows. Clearly the resulting map does not represent a deterministic forecast of future eruptions but a probabilistic estimate originating from a statistical analysis of the flank eruptive activity of Etna since Moreover, we show the results of the statistical modeling of the persistent summit activity during the last 110 years, considering Etna s summit craters both separately or as a whole. 2. Etna Volcanic Eruptions: Types, Frequencies, and Structural Features 2.1. Flank Eruptions in the Last Four Centuries [8] Before the beginning of the 17th century, the catalog of Etna flank eruptions is fragmentary and incomplete, but since 1926

3 Figure 2. Topographic map of Etna s summit area carried out in 2007 (from Neri et al. [2008], modified). Contour interval = 10 m. Dashed black line indicates the rim of the former central crater mapped in 1932 [Istituto Geografico Militare, 1934]. Light gray numbers indicate the age of the main lava flows in the period an almost complete account of eruptions is available, which has since been summarized by Sartorius von Waltershausen [1880], Imbo [1928, 1965], and more recently by Wadge [1977]. The eruptive behavior of Etna in the last 400 years is undoubtedly irregular, with important fluctuations in eruption frequencies, type of eruptions, and output rates, both in the long (centuries) and short terms (decades) [Hughes et al., 1990; Behncke and Neri, 2003; Allard et al., 2006]. [9] The period spanning from 1600 to 1669 was characterized by very high output rates (mean of ~1.4 m 3 s 1 ; Behncke and Neri, 2003) and represents the culminating stage of an earlier long cycle lasting several centuries [Tanguy et al., 2003; Clocchiatti et al., 2004]. [10] After 1669, a new long cycle began that continues today. This cycle is made up of a firstphase(uptomid s) characterized by very few, low-fed eruptions (mean effusion rate less than 0.03 m 3 s 1 ), occurring mostly at the volcano summit. Then, a second phase (up to 1865) marked an increase in the mean output rate (~0.2 m 3 s 1 )andthe resumption of a frequent flank eruptive activity (16 events). Flank eruptions and erupted volumes further increased in the third phase (mean output ~0.47 m 3 s 1 ), up to the present [Wadge et al., 1975; Guest and Murray, 1979; Allard et al., 2006]. [11] Moreover, short-term eruptive cycles have been observed since 1865 to date, a period during which flank eruptions have occurred in clusters separated by intervals of quiescence and summit activity [Behncke and Neri, 2003]. These short cycles consist of three phases too: (1) a period of eruptive quiescence lasting less than 3.5 years, (2) a period of continuous summit activity lasting 6 16 years, and finally (3) a period marked by flank eruptions alternating with summit activity, lasting 9 22 years. Each short cycle ends with ahighflux and/or voluminous flank eruption. [12] An evident increase in the lava output and frequency of eruptions began in 1950, although the behavior remained cyclic [Allard et al., 2006]. Effectively, more than 1.6 km 3 of volume of products have been erupted starting from 1950 at an average rate higher than 0.83 m 3 s 1. The last short eruptive cycle began after the flank eruption (~ m 3 volume of erupted products), evolving into a couple of years of quiescence, 6 years of gradually more frequent summit activity, followed by a series of flank (5) and summit (33) eruptions [Behncke et al., 2005; Allard et al., 2006; Neri et al., 2009, 2011; Vicari et al., 2011b]. [13] Based on numerous and detailed studies conducted on the structural features of flank eruptions at Etna [Neri et al., 2011, and references therein], we examined the eruptive fissure systems occurring between 1610 and 2012 (Figure 1). Specifically, we analyzed a total of 61 flank eruptions, which formed 130 segments of eruptive fissures (Figure 1; in some cases, the same eruption produced multiple eruptive systems with different orientations), propagated mainly southward and SSE-ward (38% of the fissures), NE-ward and ENEward (22%), eastward (9%), and westward (7%) directions (Figure 1a). During the 17th and 18th centuries, the eruptions generated 22 fissures; 27 eruptive fissures opened during the 19th century; and finally, 81 fissures formed during the 20th and 21st centuries (Figure 1b). The increased number of eruptive fissures in the 20th and 21st centuries is significant. [14] The length of eruptive systems ranged from a few hundred meters to about 11.5 km, with an average of ~3 km. About 60% of the eruptive fissures were located above 2400 m altitude; they were mainly radial with respect to the volcano summit, even if most follow a ~N-S trend, namely perpendicular to the direction of regional extension in the area of Mt Etna volcano [McGuire and Pullen, 1989; Lanzafame et al., 1996; Monaco and Tortorici, 2000; Solaro et al., 2010]. At lower altitude (below 2400 m), the eruptive fissures were clustered in SE, SSE, NE, ENE, east, and west alignments, following the rift zones. The minimum altitude (~650 m above sea level) of an eruptive vent was reached during the 1669 eruption, close to Nicolosi town (southern flank, Figure 1). In many cases, eruptive fissures do not erupt along their entire lengths (except at the very beginnings of the event); during the propagation of the fissure, the eruptive activity migrates downward, and the main vents are commonly located on the lower portions of the fissure itself. Table 1. Number and Type of Eruptions for Each Summit Crater a Crater Paroxysmal Episodes Strombolian Activities Lava Flows Central crater (since 1900) NE crater (since 1911) SE crater (since 1971) a The central crater includes the activity of Voragine and Bocca Nuova. The last column (Lava Flows) indicates the number of paroxysmal and strombolian events that were also characterized by (1) rheomorphic lavas formed during highly fed explosive activity; (2) lava overflows emerging from one of the summit crater; and (3) lavas erupted from short fissures opened along the external slope of the summit cone (above ~3000 m asl). 1927

4 Figure 3. Kolmogorov-Smirnov test of stationarity for Etna s flank eruptions over the last four centuries. The green step function is the empirical distribution of the eruption onsets, the blue straight line is the theoretical uniform distribution, and the light blue stripes are the 95% confidence bands. Red dots indicate the repose times in years (see secondary Y axis on the right side) between the openings of successive eruptive fissures. The straight red line highlights the decreasing trend of repose times, i.e., eruptions occurring more frequently in recent years than in the past centuries Persistent Summit Activity During the Last 110 Years [15] The summit area of Mt Etna has frequently undergone major morphological changes due to its persistent eruptive activity, both effusive and explosive [Neri et al., 2008]. Four summit craters are active today: the Voragine, the NE crater, the Bocca Nuova, and the SE crater (Figure 2). Their eruptive activity consists of the following: (1) continuous degassing; (2) mild Strombolian explosions accompanied by low rate (~0.1 1m 3 s 1 ) lava effusion from vents inside or on the flanks of the summit craters; and (3) high-rate paroxysmal (up to m 3 s 1 ) eruptions, characterized by lava fountains, lava overflows, and tephra-rich eruption columns several kilometers high. [16] This morphology has not always been as it is today. Just 100 years ago, there was only one summit crater: the central crater (CC; Guest [1973]). Indeed, during the 17th 18th centuries, the summit of the volcano had a single deep and wide crater depression, which was formed because of the collapse of the summit area triggered by the huge 1669 eruption. Later on, during the 19th 20th centuries, the volcanic activity resumed within this crater, causing its gradual filling. [17] During the almost persistent activity of the central crater (CC), the NE crater (NEC) formed in At the beginning, it was a pit located ~3160 m above sea level (asl) on the northern slope of the CC cone [Guest, 1973] and characterized mainly by Strombolian activity. NEC became very active in the second half of the 20th century, producing numerous and occasionally violent and long eruptions (especially in , , and ), up to becoming a major pyroclastic cone as well as the highest point of the volcano (3330 m asl in 2007; Neri et al., 2008). [18] Starting from the middle of the 20th century, the CC changed its internal morphology frequently due to the intracrater volcanic activity alternating with subsidence and collapse phenomena. One of these events led to the creation of the Voragine (VOR; also known as The Chasm ) in VOR formed in the northeastern portion of the fairly flat floor of CC [Guest, 1973]. Major eruptive phases occurred at CC in , which culminated in the opening of eruptive fissures, several paroxysmal events accompanied by huge tephra emissions, continuous lava overflows, and the growth of some pyroclastic intracrater cones [Behncke et al., 2004]. Four years later, in 1968, the Bocca Nuova (BN) was formed, still inside the CC, a few tens of meters west of the VOR [Behncke et al., 2003]. BN widened and deepened gradually, always within the CC, becoming even larger than the VOR. Later, in 1998, the Strombolian activity resumed at the crater floor, generating several lava fountaining episodes (in 1999) that filled the crater and triggered numerous lava overflows [Behncke et al., 2003; Calvari et al., 2003]. After the 2001 flank eruption, BN and VOR were involved in subsidence and sporadic phreatomagmatic explosions, and progressively they collapsed, becoming almost a single crater depression whose edges are, in part, coincident with those of the CC [Giammanco et al., 2007]. Therefore, it is evident that VOR and BN are internal, minor vents of the CC itself (see dashed blacklineinfigure2). [19] The SE crater (SEC) formed during the 1971 eruption. Since its creation, SEC has been the most active 1928

5 Figure 4. Centroids of the eruptive fissures opened in the last 400 years, grouped in short periods of 50 years. A spatial shift is evident over time in the location of the eruptive activity, which becomes progressively closer to the summit of the volcano, on the high eastern side. Coordinates are in Universal Transverse Mercator (UTM) projection, zone 33 N. of the summit craters of the volcano. At first, it was a degassing pit located close to the southeast base of the CC cone, at ~3070 m asl [Behncke et al., 2006]. During the first 25 years of activity, SEC erupted quite frequently and built a cone about 100 m high. The cone has grown dramatically in the periods and , through over 143 paroxysmal episodes and effusive eruptions. Today it has reached an elevation of ~3300 m, forming also a huge new pyroclastic cone on its eastern slope [Alparone et al., 2005; Behncke et al., 2006; Neri et al., 2008; Vicari et al., 2011b; Ganci et al., 2012a]. [20] The progressive increase in eruptive activity characterizing the Mt Etna summit area during the last 110 years can be readily seen in Table 1. What is most striking is the fact that the appearance of the SEC marks a real change in the eruptive frequency. Over the past 41 years, the SEC has produced a number of eruptive events that correspond to more than double the events produced jointly by the CC and NEC in 110 years. Never in the recent past has Mt Etna been so active; to find similar levels of activity, we would probably have to go back several hundred years, when the historical chronicles, however, were certainly not as accurate as today. 3. Statistical Analysis of Flank Eruptions [21] Over the last four centuries, the completeness and accuracy of historical data at Mt Etna provides a robust basis for conducting a statistical analysis to detect spatial and temporal patterns in the distribution of past events and hence to formulate the appropriate spatiotemporal probabilistic model for future flank eruptions. [22] We analyzed the spatial pattern of eruptive fissures using the Clark and Evans [1954] aggregation index R that measures the ratio between the expected nearest-neighbor distances under complete spatial randomness (CSR) and the observed mean nearest-neighbor distances:! = R ¼ XN ðd min ðþ=n i Þ 0:5 ð A=N Þ1=2 (1) i¼1 where d min (i) is the distance between the ith eruptive fissure and its nearest neighbor, A is the area of the region, and N is the number of structures considered in the calculation. In order to consider the eruptive fissures in full, d min is calculated as the minimax distance, that is, the minimum value of the maximum distances between each end point of the ith eruptive fissure and all the end points of the jth eruptive fissure [Cappello et al., 2012]. Table 2. Estimates of the Parameters d and θ a Parameters Minimum Estimates Best Estimates Maximum Estimates d θ l(t 2 + 1) l(t ) l(t ) a Calculated by minimizing the square differences of residuals between the expected and actual number of eruptions over time (best estimates), with the 95% confidence intervals (minimum and maximum estimates).the evaluations of d and θ are then used to calculate the recurrence rate l(t 2 +p) with T 2 = 2012 and p = 1, 10,

6 Figure 5. (a) Cumulative number of eruptions (black curve); best estimate calculated by minimizing the square differences of residuals between the expected and actual number of eruptions over time (red curve); and the 95% confidence interval (blue curves) of the best estimate. (b) The corresponding power intensity functions l(t). The future rates (listed in detail in Table 2) are obtained by extending the curves for the different forecasting periods highlighted by the dashed lines over the gray bar, i.e., 2013 and the next 10 and 50 years. [23] R values less than 1 highlight a cluster tendency, since the observed mean distance between neighboring fissures is less than that expected in random patterns. R equal to 1 points out that casualness is present among data, while R greater than 1 indicates uniformity. The theoretical maximum of R is occurring when data are maximally dispersed. With historical data at Etna, we obtained an R value equal to 0.027, which proves the nonhomogeneous spatial distribution of eruptive fissures. [24] In addition to the spatial analysis, we examined the historical records during the last four centuries to detect the presence of trends or patterns characterizing the temporal distribution of eruption occurrences. We calculated the cumulative number of eruptive fissure systems linked to flank eruptions over time and the inter-event periods, i.e., the repose times between subsequent eruptions (Figure 3). Both distributions underline the temporal nonstationarity of eruption occurrences, since a growing trend is evident for the number of flank eruptions, while a mainly decreasing tendency stands out for repose times during the last 30 years. This highlights that volcanic eruptions occurred with a higher frequency in this last period than in the previous ~370 years of activity. To formally prove that the sequence of eruptive events is not a stationary process, we applied the Kolmogorov-Smirnov test statistic D n [Mulargia et al., 1985], which is based on the maximum difference between the theoretical (uniform) and the empirical cumulative distribution functions. Since the historical records of Etna s flank eruptions provided D n = 0.267, which is higher than the critical value at a 5% significance level (0.174 = 1.36/ n for a large number n of events), we can reject the null hypothesis of uniform distribution (Figure 3). [25] In order to evaluate the possible shift during time of the location of the eruptive activity, we also calculated the centroids closer to the eruptive fractures in the last 400 years, grouping them every 50 years (Figure 4). The result indicates that the increasing trend in volcanic activity is accompanied by a migration of the centroids that became progressively closer to the summit of the volcano, on the high eastern side. This fact is testified by the rapid morphostructural evolution of the summit area during the last decades: the former central crater (the single summit crater existing up to early 1900) was modified by the growth of the NE crater (formed in 1911 on the northern side), the Bocca Nuova (formed in 1968 inside the central crater and close to the its western rim, see Figure 2), and the SE crater (formed in 1971 on the southeastern side), which represent the connections of the three main volcanic rifts in the summit area [Neri et al., 2011]. 1930

7 Figure 6. Spatiotemporal probability maps of flank vent opening for the next (a) 1, (b) 10, and (c) 50 years. A gradual increase in the probabilities is evident as the forecasting period increases. The highest value estimated with our model is very low for 2013 ( ), reaches a medium peak for the next 10 years (0.018), and attains the top value in the spatiotemporal probability map for the next 50 years (0.091). Summit craters are masked because their activity was separately investigated. [26] The statistical analysis conducted on flank eruptions in the last four centuries demonstrates the spatial nonhomogeneity and the temporal nonstationarity of flank activity, hence suggesting that the most appropriate space-time probabilistic model for future flank eruptions at Etna volcano is a nonhomogenous Poisson process (NHPP) with a timevarying intensity function. 4. The Space-Time Probabilistic Modeling of Flank Eruptions [27] Probabilistic modeling provides a powerful tool to estimate the timing and location of future eruptions. Typically, modeling exploits statistical analysis on historical information to calculate probability estimates over a particular area of interest and construct a spatiotemporal map furnishing detailed recurrence rates, i.e., events expected per unit area per unit time. [28] Here the area of interest on Etna volcano is 1170 km 2 (rectangular dotted area of km in Figure 1), over which we superimposed a regular grid of potential vents yielding a resolution of 0.5 km. For each potential vent (x, y), the probability of at least a new vent opening in the unit area ΔxΔy (0.5 km 0.5 km) within time Δt is calculated with the NHPP model as follows: p a ðx; y; Δt Þ ¼ 1 exp l xyt ΔxΔyΔt with l xyt the space-time varying intensity function calculated as l xyt ¼ XN i¼1 exp di 2 = 2h2 exp ð kti Þ where d i is the distance between the point (x, y) andtheith eruptive fissure, h is the smoothing factor, t i is the elapsed time since the eruption at the ith fissure, N is the total number of fissures, and k is a constant that assigns larger weights to more recent events with respect to older ones. (2) (3) [29] The smoothing factor h (also called bandwidth or window width) strictly influences the result of equation (3), controlling how l xyt varies with distance from historical eruptive fissures. We assigned it a value of 1 km by applying the explicit version of the Least Squares Cross-Validation (LSCV) proposed by Worton [1995], which is based on minimizing the integrated square error between the true and the estimated distribution: LSCVðhÞ ¼ 1 ph 2 N þ 1 4ph 2 N 2 XN X N exp dij 2 = 4h2 4 exp dij 2 = 2h2 i¼1 j¼1 where N is the total number of structures considered in the calculation and d ij is the minimax distance. [30] In order to calculate the constant k, we conducted a retrospective search to quantify how the distribution of the oldest eruptive fissures has affected the formation of the most recent ones. The approach used, called back analysis [Cappello et al., 2011b], allowed us to empirically retrieve the best value of the parameter k using only the known past data. We divided all eruptive fissures according to the age of their formation in two parts and used the first one as a learning Table 3. Estimates of the Best Values for d, θ (Best Estimates), and the 95% Confidence Intervals (Minimum and Maximum Estimates) and the Relative Recurrence rates l(t 2 +p)witht 2 = 2012 and Forecasting Periods p = 1, 10 Years a Parameters Minimum Estimates Best Estimates Maximum Estimates d θ l(t 2 + 1) l(t ) a The estimates are calculated considering the activity of summit craters as a whole. (4) 1931

8 Figure 7. Actual cumulative number of summit eruptions (black curve) compared with the best estimate (red curve) and its 95% confidence interval (blue curves). The formations of the NE crater in 1911 and of the SE crater in 1971 are highlighted. (b) Power intensity functions l(t) calculated with the temporal rates for the next 1 and 10 years reported in Table 3. data set and the second one as a testing data set. The cutoff time we chose to divide the eruptive fissures is 1981, since it determines two meaningful sets from a statistical point of view: 50 fissures were used for learning and 11 for testing (i.e., the 1983, March 1985, December 1985, , 1989, , 2001, , , 2006, and eruptions). Moreover, from a volcanological point of view, activity at Etna since 1981 shows a considerably increasing trend that further justified our choice [Behncke and Neri, 2003; Behncke et al., 2005; Allard et al., 2006; Salvi et al., 2006; Neri et al., 2008, 2009; Cappello et al., 2011b]. [31] The best fit for the constant k was retrieved by evaluating the mean square error of each probability forecast with the Brier score [Brier, 1950]: BS ¼ 1 N X N i¼1 ðp i o i Þ 2 (5) where p i is the forecast probability of each grid cell and o i indicates whether an eruptive fissure has been formed in that cell after 1981 (o i = 1 if at least a fissure has occurred, and o i = 0 if not). The Brier score runs from 0 to 1, with smaller values indicating better forecasts. We found that the value for k that minimizes BS is and used it in equation (3) for all eruptive fissures. [32] The function l xyt is then rescaled so that its integral across the whole area equals the temporal intensity calculated with the power intensity function [Smethurst et al., 2009; and reference therein]: lðþ¼ t d t T d 1 1 for t T 1 (6) θ θ where T 1 is the year of the oldest eruptive fissure, t is the year of interest, and d and θ are unknown positive parameters. The parameter d determines the shape of the power function: d > 1 (resp. d < 1) corresponds to an increasing (resp. decreasing) intensity with time, while d = 1 reduces equation (6) to a homogenous process. The best estimates of d and θ (Table 2) are found by minimizing the square differences of residuals between the expected amount (computed as the integral of l from T 1 to t) and the actual number of eruptions over time (estimated as the cumulative number of events observed up to t) on the entire catalog (from T 1 =1610 to T 2 = 2012; see Figure 5a). In order to indicate the reliability of the estimates of d and θ, we also listed the 95% confidence intervals (minimum and maximum estimates in Table 2). Using the obtained values of d and θ, we extrapolated the curves of annual recurrence rates for different forecasting periods p (1, 10 or 50 years) and calculated the minimum, best, and maximum temporal intensities evaluating equation (6) 1932

9 Figure 8. Actual (blue curve) and modeled (red curve) cumulative number of eruptions for the (a) central crater (CC), (b) the NE crater (NEC), and (c) the SE crater (SEC). The power intensity functions calculated with the best rates for the next 1 and 10 years are also shown respectively for (d) CC, (e) NEC, and (f ) SEC. for t = T 2 + p (Figure 5b). Obviously, the estimated values of d > 1 (Table 2) are only representative of the increasing phase of the current eruptive cycle. Indeed, it is not physically plausible that the intensity will increase continuously in the next years due to the cyclic variations in output rate and frequency of Etna s eruptive activity [Wadge et al., 1975; Hughes et al., 1990; Behncke and Neri, 2003; Allard et al., 2006]. 5. Spatiotemporal Probability Maps for Future Vent Opening [33] The spatiotemporal probability hazard maps for future vent opening for the three different forecasting periods (1, 10 or 50 years) are reported in Figure 6. The temporal recurrence rates are those obtained with the best estimates of d and θ (Table 2). The highest value reached by the probability map of vent opening for the coming year is fairly low, at This is bounded by the temporal recurrence rate (0.216), which predicts about one eruption every 5 years. As is evident in Figures 6b and 6c, the probability estimates gradually increase when the forecasting periods are 10 years (maximum value of 0.018) and 50 years (maximum value of 0.091). Nevertheless, the global trend of the probability distributions is kept, showing that the areas closer to the summit of the volcano, above 2500 m above sea level, are where it is very probable that a new eruption will take place. Dissimilarity in the probability estimates can be noticed also in the areas at lower altitude, between 2500 and 1800 m above sea level. Below 1000 m, the expectation of at least a new eruption is very low, except in the South Rift, where the probabilities slowly decrease until the altitude of ~600 m above sea level. This effect is caused by the presence of the 1669 eruptive fissure. 6. Statistical Modeling of Summit Eruptions [34] The persistent summit activity documented at Etna from the first decades of the 20th century until today represents a good source for a statistical analysis. Thus, here we provide a more comprehensive and accurate study in forecasting volcanic eruptions, estimating also the space-time annual rate of eruptive events at the summit craters (Figure 2): the CC (comprising VOR and BN), the NEC, and the SEC (including the new cone recently formed on its eastern slope). [35] The data set for the statistical analysis comprises the various types of eruptions occurring since 1900 at the summit craters of Etna, including paroxysmal episodes and Strombolian activities (see Table 1). We conducted two kinds of analysis on this data set. First we calculated the recurrence rate of the summit activity considering the whole catalog of eruptions. Next we divided the data set into three and estimated space-time forecasts for each summit crater (CC, NEC, and SEC). [36] As for the recurrence rate of the summit activity as a whole (without distinguishing the activity of each summit crater), we estimated the best values for the parameters d and θ (Table 3) minimizing the square differences of residuals between the expected and the actual number of eruptions over time (Figure 7a). These best estimates were then used in the power intensity function described by equation (6) to obtain the curves of recurrence rates with the expected number of eruptions for the next 1 and 10 years (Figure 7b). Unlike the analysis of flank eruptions, we decided not to use forecasting periods longer than 10 years since they might Table 4. Best Values Obtained for d and θ and the Relative Recurrence Rates l(t 2 +p) with T 2 = 2012 and Forecasting Periods p = 1, 10 Years a Parameters Central Crater NE Crater SE Crater d θ l(t 2 + 1) l(t ) a The estimates are calculated considering the activity of the central crater (since 1900), of the NE crater (since 1911), and of the SE crater (since 1971), separately. 1933

10 provide inconsistent results starting from the short time series of documented summit eruptions (~110 years). The best values of recurrence rates for the coming year ( eruptions/year) and the next 10 years ( eruptions/year) estimate that slightly more than one eruption per month is expected in the entire summit area. In the worst cases, more than two (25.530) and about three eruptions (35.268) per month could occur, respectively, in the next 1 and 10 years. These estimates are strongly influenced by the extraordinary activity that took place at SEC in 2000 (66 short-term lava fountains) and recently between 2011 and 2012 (25 midterm lava fountains). [37] For the space-time forecasts, we divided the catalog in three data sets selecting the eruptions occurring at CC, NEC, and SEC. With these three data sets, we used the same mathematical modeling, calculating the best expected number of eruptions on the basis of the historical catalogs (Figures 8a 8c), from which we derived the best estimates for parameters d and θ (Table 4). The recurrence rates were again estimated for the coming year and the next 10 years (Figures 8d 8f). The spatiotemporal analysis suggests that the expected number of eruptions at the CC (from to 0.484) is comparable with the rates obtained for the flank eruptions (from to 0.219). Unlike those of the NEC and of the SEC, the curve of the annual rate for the CC shows a decrease in the speed of growth. At the NEC, the expected number of eruptions reaches annual rates of and 1.641, respectively, for the next 1 and 10 years. But the summit crater showing the most exceptional activity is undoubtedly at the SEC. On the basis of the historical eruptions since its formation in 1971, we expected a number of eruptions equal to for the coming year, that is, more than one eruption per month, and eruptions in the next 10 years, namely slightly less than two events per month. 7. Conclusive Remarks [38] We have demonstrated a new way to construct the spatiotemporal probability map of vent opening over a specific time period. The spatiotemporal probability map furnishes detailed recurrence rates, estimating the number of expected events per unit area per unit time. In order to find the areas where it is more likely that a new eruption will occur, the best approach consists in analyzing the past eruptive history of a volcano to deduce the possible location of future eruptions. [39] We based the construction of the spatiotemporal probability map of vent opening on an exhaustive statistical analysis and a NHPP model applied to the Etna flank eruptions of the last 400 years. The spatial distribution of past eruptive fissures is studied, as well as the spatial shift in the effusive eruptions, which appears evident in the last 30 years of Etna activity. An innovation of our approach consists in forecasting the location of new eruptive vents by weighting the most recent events more heavily than the older ones. If a new eruption occurs, the map can simply be updated by merging past information with the timing and location of the new vent. In order to make major changes to the map, a statistically significant number of eruptions would have to occur at different unexpected locations in the study area. The spatiotemporal probability maps of vent opening are quite different in content and meaning from the susceptibility map produced by Cappello et al. [2012], which estimates only the most likely emission zones. Here, we also considered the frequency of flank eruptions to evaluate the future timing and areas of Etna prone to vent opening for the coming year and the next 10 and 50 years. We chose three different forecasting periods to highlight how the probability estimates change over time. Since the catalog of flank eruptions at Etna is well studied and documented only for the last four centuries, forecasting periods longer than 50 years would provide results lacking in significance from a statistical point of view. [40] Even though the focus of this study is Mt Etna, the approach used to assess the most exposed areas to be affected by new eventual eruptions is meant to be general and applicable to any volcanic region, if available data are significant from a statistical point of view. The spatiotemporal probability map is extremely important in volcanic hazards assessment because the degree of danger presented by a hazardous event strongly depends on the eruption vent location. Hence, the more rigorous the evaluation of the spatiotemporal probability map of vent opening, then the greater the accuracy of the hazard map. [41] Finally, the results of the analysis of the persistent summit activity during the last 110 years highlight that the hazard rate for eruptive events is not constant with time and differs for each summit crater of Mt Etna. The variations in the recurrence rates at the summit craters indicate that in the last century the main eruptive activity was initially located at the central crater, then migrated to the NE crater and later to the SE crater, as the new summit craters were formed. Most likely, the SE crater will remain the most active crater in the summit area for the near future. [42] Acknowledgments. This work was developed in the frame of the TecnoLab, the Laboratory for the Technological Advance in Volcano Geophysics organized by INGV-CT, DIEES-UNICT, and DMI-UNICT. We thank Geoff Wadge and an anonymous reviewer for constructive and supportive comments that helped to improve the manuscript. References Acocella, V., and M. Neri (2003), What makes flank eruptions? The 2001 Etna eruption and its possible triggering mechanisms, Bull. Volcanol., 65, , doi: /s Allard, P., B. Behncke, S. D Amico, M. Neri, and S. Gambino (2006), Mount Etna : Anatomy of an evolving eruptive cycle, Earth- Sci. Rev., 78, , doi: /j.earscirev Alparone, S., B. Behncke, S. Giammanco, M. Neri, and E. Privitera (2005), Paroxysmal summit activity at Mt. Etna monitored through continuous soil radon measurements. Geophys. Res. Lett., 32, L16307, doi: / 2005GL Bebbington, M. S. (2007), Identifying volcanic regimes using hidden Markov models, Geophys. J. Int., 171, , doi: / j x x. Behncke, B., and M. Neri (2003), Cycles and trends in the recent eruptive behaviour of Mount Etna (Italy), Can. J. Earth Sci., 40, , doi: /e Behncke, B., M. Neri, and R. Carniel (2003), An exceptional case of lava dome growth spawning pyroclastic avalanches at Mt. Etna (Italy): The 1999 Bocca Nuova eruption, J. Volcanol. Geotherm. Res., 124/1 2, , doi: /s (03) Behncke, B., M. Neri, and G. Sturiale (2004), Rapid morphological changes at the summit of an active volcano: Reappraisal of the poorly documented 1964 eruption of Mount Etna (Italy), Geomorphology, 63/3 4, , doi: /j.geomorph Behncke, B., M. Neri, and A. Nagay (2005), Lava flow hazard at Mount Etna (Italy): New data from a GIS-based study, in Kinematics and dynamics of lava flows, edited by M. Manga and G. Ventura, Spec. Pap. Geol. Soc. Am., 396, , doi: / (13). 1934

11 Behncke, B., M. Neri, E. Pecora, and V. Zanon (2006), The exceptional activity and growth of the Southeast Crater, Mount Etna (Italy), between 1996 and 2001, Bull. Volcanol., 69, , doi: /s x. Bisson, M., B. Behncke, A. Fornaciai, and M. Neri (2009), Lidar-based digital terrain analysis of an area exposed to the risk of lava flow invasion: The Zafferana Etnea territory, Mt. Etna (Italy), Natural Hazards, 50(2), , doi: /s Bonaccorso, A., A. Bonforte, S. Calvari, C. Del Negro, G. Di Grazia, G. Ganci, M. Neri, A. Vicari, and E. Boschi (2011), The initial phases of the Mt. Etna eruption: A multi-disciplinary approach for hazard assessment, J. Geophys. Res., 116, B03203, doi: /2010jb Brier, G. W. (1950), Verification of forecasts expressed in terms of probabilities, Mon. Wea. Rev., 78, 1 3, doi: / (1950)078 < 0001: VOFEIT > 2.0.CO;2. Calvari, S., M. Neri, and H. Pinkerton (2003), Effusion rate estimation during the 1999 summit eruption on Mt. Etna, and growth of two distinct lava flow fields, J. Volcanol. Geotherm. Res., 119/1 4, , doi: /s (02) Cappello, A., A. Vicari, and C. Del Negro (2011a), Assessment and modeling of lava flow hazard on Etna volcano, Boll. Geof. Teor. Appl., 52(2), , doi: /bgta0003. Cappello, A., A. Vicari, and C. Del Negro (2011b), Retrospective validation of a lava flow hazard map for Mount Etna volcano, Ann. Geophys., 54(5), , doi: /ag Cappello, A., M. Neri, V. Acocella, G. Gallo, A. Vicari, and C. Del Negro (2012), Spatial vent opening probability map of Mt Etna volcano (Sicily, Italy), Bull. Volcanol., 74, , doi: /s Clark, J., and F. C. Evans (1954), Distance to nearest neighbour as a measure of spatial relationships in populations, Ecology, 35, , doi: / Clocchiatti, R., M. Condomines, N. Guenot, and J. C. Tanguy (2004), Magma changes at Mount Etna: The 2001 and eruptions, Earth Planet. Sci. Lett., 226, , doi: /j.epsl Crisci, G. M., M. V. Avolio, B. Behncke, D. D Ambrosio, S. Di Gregorio, V. Lupiano, M. Neri, R. Rongo, and W. Spataro (2010), Predicting the impact of lava flows at Mount Etna, Italy, J. Geophys. Res., 115, B04203, doi: /2009jb Favalli, M., S. Tarquini, A. Fornaciai, and E. Boschi (2009), A new approach to risk assessment of lava flow at Mount Etna, Geology, 37(12), , doi: /g30187a. Ganci, G., A. J. L. Harris, C. Del Negro, Y. Guehenneux, A. Cappello, P. Labazuy, S. Calvari, and M. Gouhier (2012a), A year of lava fountaining at Etna: Volumes from SEVIRI, Geophys. Res. Lett. 39, L06305, doi: /2012gl Ganci, G., A. Vicari, A. Cappello, and C. Del Negro (2012b), An emergent strategy for volcano hazard assessment: From thermal satellite monitoring to lava flow modeling, Remote Sens. Environ., 119, , doi: /j.rse Giammanco, S., K. W. W. Sims, and M. Neri (2007), Measurements of 220 Rn and 222 Rn and CO 2 emissions in soil and fumarole gases on Mt. Etna volcano (Italy): Implications for gas transport and shallow ground fracture, Geochem. Geophys. Geosyst., 8, Q10001, doi: / 2007GC Greco, F., G. Currenti, C. Del Negro, R. Napoli, G. Budetta, M. Fedi, and E. Boschi (2010), Spatiotemporal gravity variations to look deep into the southern flank of Etna volcano, J. Geophys. Res., 115, B11411, doi: /2009jb Guest, J. E. (1973), The summit area of Mount Etna prior to the 1971 eruption, Phil. Trans. Roy. Soc. London, A, 274, Guest, J. E., and J. B. Murray (1979), An analysis of hazard from Mount Etna volcano, J. Geol. Soc. London, 136, , doi: / gsjgs Hughes, J. W., J. E. Guest, and A. M. Duncan (1990), Changing styles of effusive eruption on Mount Etna since AD 1600, in Magma Storage and Ascent, edited by M. P. Ryan, pp , John Wiley, Chichester. Imbo, G. (1928), Sistemi eruttivi etnei, Bulletin Volcanologique, Series I, 5, Imbo, G. (1965), Catalog of Active Volcanoes of the World, 28, Italy, IAVCEI, 18, Istituto Geografico Militare (1934), Cono Terminale del Monte Etna, Settembre 1932, Anno X. scala 1:5000. Lanzafame, G., M. Neri, and D. Rust (1996), Active Tectonics Affecting the Eastern Flank of Mount Etna: Structural Interactions at Regional and Local Scale, in Etna: Fifteen Years On, edited by P. H. Gravestock and W. J. McGuire, pp , Cheltenham and Gloucester College of Higher Education, Centre for Volcanic Research, Cheltenham, U.K. Magill, C. R., K. J. McAneney, and I. E. M. Smith (2005), Probabilistic assessment of vent locations for the next Auckland volcanic field event, Mathematical Geol., 37(3), , doi: /s McGuire, W. J., and A. D. Pullen (1989), Location and orientation of eruptive fissures and feeder dykes at Mount Etna; influence of gravitational and regional tectonic stress regimes, J. Volcanol. Geotherm. Res., 38, , doi: / (89) Monaco, C., and L. Tortorici (2000), Active faulting in the Calabrian arc and eastern Sicily, J. Geodynamics, 29, , doi: /s (99) Mulargia, F., S. Tinti, and E. Boschi (1985), A statistical analysis of flank eruptions on Etna Volcano, J. Volcanol. Geotherm. Res., 23, Napoli, R., G. Currenti, C. Del Negro, F. Greco, and D. Scandura (2008), Volcanomagnetic evidence of the magmatic intrusion on 13th May 2008 Etna eruption, Geophys. Res. Lett., 35, L22301, doi: / 2008GL Neri, M., F. Mazzarini, S. Tarquini, M. Bisson, I. Isola, B. Behncke, and M. T. Pareschi (2008), The changing face of Mount Etna s summit area documented with LiDAR technology, Geophys. Res. Lett., 35, L09305, doi: /2008gl Neri, M., F. Casu, V. Acocella, G. Solaro, S. Pepe, P. Berardino, E. Sansosti, T. Caltabiano, P. Lundgren, and R. Lanari (2009), Deformation and eruptions at Mt. Etna (Italy): A lesson from 15 years of observations, Geophys. Res. Lett., 36, L02309, doi: /2008gl Neri, M., V. Acocella, B. Behncke, S. Giammanco, F. Mazzarini, and D. Rust (2011), Structural analysis of the eruptive fissures at Mount Etna (Italy), Ann. Geophys., 54(5), , doi: /ag Newhall, C. G. (2000), Volcano warnings, in Encyclopedia of Volcanoes, edited by H. Sigurdsson, pp , Academic Press, San Diego. Rongo, R., M. V. Avolio, B. Behncke, D. D Ambrosio, S. Di Gregorio, V. Lupiano, M. Neri, W. Spataro, and G. M. Crisci (2011), Defining high detailed hazard maps by a cellular automata approach: Application to Mt. Etna (Italy), Ann. Geophys., 54(5), , doi: /ag Salvi, F., R. Scandone, and C. Palma (2006), Statistical analysis of the historical activity of Mount Etna, aimed at the evaluation of volcanic hazard, J. Volcanol. Geotherm. Res., 154, Sartorius von Waltershausen, W. (1880), Der Aetna, edited by A. von Lasaulx. W. Engelmann, Leipzig, Germany. Scifoni, S., M. Coltelli, M. Marsella, C. Proietti, Q. Napoleoni, A. Vicari, and C. Del Negro (2010), Mitigation of lava flow invasion hazard through optimized barrier configuration aided by numerical simulation: The case of the 2001 Etna eruption, J. Volcanol. Geotherm. Res., 192, 16 26, doi: /j.jvolgeores Smethurst, L., M. R. James, H. Pinkerton, and J. A. Tawn (2009), A statistical analysis of eruptive activity on Mount Etna, Sicily, Geophys. J. Int., 179, , doi: /j x x. Solaro, G., V. Acocella, S. Pepe, J. Ruch, M. Neri, and E. Sansosti (2010), Anatomy of an unstable volcano through InSAR data: Multiple processes affecting flank instability at Mt. Etna in , J. Geophys. Res., 115, B10405, doi: /2009jb Sparks, R. S. J. (2003), Forecasting volcanic eruptions, Earth Planet, Sci. Lett., 210, 1 15, doi: /s x(03) Tanguy, J. C., M. Le Goff, C. Principe, S. Arrighi, V. Chillemi, A. Paiotti, S. La Delfa, and G. Patanè (2003), Archeomagnetic dating of Mediterranean volcanics of the last 2100 years: Validity and limits, Earth Planet, Sci. Lett., 211, , doi: /s x(03) Tanguy, J.-C., M. Condomines, M. Le Goff, V. Chillemi, S. Le Delfa, and G. Patanè (2007), Mount Etna eruptions of the last 2,750 years: Revised chronology and location through archeomagnetic and 226Ra-230Th dating, Bull. Volcanol., 70, 55 83, doi: /s x. Vicari, A., G. Bilotta, S. Bonfiglio, A. Cappello, G. Ganci, A. Hérault, E. Rustico, G. Gallo, and C. Del Negro (2011a), LAV@HAZARD: A Web-Gis interface for volcanic hazard assessment, Ann. Geophys., 54, 5, doi: /ag Vicari, A., G. Ganci, B. Behncke, A. Cappello, M. Neri, and C. Del Negro (2011b), Near-real-time forecasting of lava flow hazards during the January 2011 Etna eruption, Geophys. Res. Lett., 38, L13317, doi: /2011gl Wadge, G. (1977), The storage and release of magma on Mount Etna, J. Volcanol. Geotherm. Res., 2, Wadge, G., G. P. L. Walker, and J. E. Guest (1975), The output of the Etna volcano, Nature, 255, , doi: /255385a0. Wadge, G., P. A. V. Young, and I. J. Mckendrick (1994), Mapping lava flow hazards using computer simulation, J. Geophys. Res., 99, , doi: /93jb Worton, B. J. (1995), Using Monte Carlo simulation to evaluate kernel based home range estimators, J. Wild. Manag., 59,