In this study, using the tephra dispersal model HAZMAP, we investigate the effect of using different meteorological datasets and eruption source parameters on tephra fallout hazard assessment for a sub-Plinian eruption of Vesuvius, which is considered as a reference case for hazard assessment analysis. We analyze the effect of using different meteorological data, from: i) radio-sounding carried out at the meteorological station of Brindisi (Italy) between 1962 and 1976 and between 1996 and 2012, and at Pratica di Mare (Rome, Italy) between 1995 and 2013; ii) meteorological models of the National Oceanic and Atmospheric Administration (NOAA), and of the European Centre for Medium-Range Weather Forecasts (ECMWF). Furthermore, we consider the effects of perturbing reference eruptive source parameters. In particular, we vary the total mass, the total grain-size distribution, the column height, and the effective atmospheric diffusion coefficient to evaluate how these parameters affect the hazard probability maps. Moreover, the effect of the seasonal variation of the wind field and the effect of the rain on the deposit loading are considered. Results show that the parameter that mostly affects hazard maps is, as expected, the total erupted mass; furthermore, keeping constant the erupted mass, the most important control on hazard is due to the particle terminal settling velocity distribution which is a function of the total grain-size distribution, particle density and shape. Within the considered range variations, the hazard depends less on the use of different meteorological datasets, column height and effective diffusion coefficient.
In this paper we present the main outcomes of a detailed investigation of the effects of meteorological datasets and key eruption source parameters on the fallout hazard maps for the reference scenario adopted for the Emergency Plan of Vesuvius (DPC 2001; 2015; Regione Campania 2015). The goal is to highlight the role of the uncertainty and the natural variability on the probability maps commonly used for hazard assessment.
In order to analyze the effect that the use of different wind datasets and the variation of the main eruptive source parameters have on the fallout hazard maps, we used the HAZMAP model (Macedonio et al. 2005). The HAZMAP model (freely available at ) is based on the solution of a simplified 2D equations of diffusion, transport and sedimentation of small particles, describing the dispersion of tephra generated by a volcanic convective column. The HAZMAP code was previously used to produce hazard maps at Vesuvius (Macedonio et al. 2008) and Campi Flegrei (Costa et al. 2009).
The investigation of the effects of tephra transport model on the hazard maps is outside the scope of this work, which is mainly focused on the analysis of the effects of the meteorological and the volcanological source parameters. In a previous work, Scollo et al. (2008) found that, at Etna volcano, larger differences between the output of tephra transport models are related to topographic effects, which strongly affect the dispersal of fine particles for column heights lower than 12 km; however these difference decrease in simulating higher eruption columns (Scollo et al. 2008). With respect to Scollo et al. (2008), who deal with a sensitivity analysis of tephra transport at Etna for small to moderate eruptions, the present work is focused on a sub-Plinian eruption at Vesuvius.
We performed a statistical analysis of wind speeds and directions in the Vesuvius area using data obtained from different datasets derived from radio-sounding stations and from reanalysis. Radio-sounding launches are routinely performed by the Aeronautica Militare (the Italian Air-Force); these data can be downloaded from the web-site of the University of Wyoming (USA) or the NOAA data center either in ASCII or NetCDF formats (Pincus and Rew 2010; Rew et al. 2010). Our investigation considered the following meteorological data sets:
radio-soundings at Brindisi (LIBR station) in the period between 1 July 1996 and 31 December 2012. This dataset consists of 18,482 records. Only records containing data between ground level and 25,000 m above sea level, and with more than 8 data in this altitude range, were selected for this work, resulting in a dataset containing 12,788 records.
radio-soundings at Pratica di Mare, Rome, about 190 km North-West of Vesuvius (LIRE station) in the period between 1 November 1995 and 21 December 2013. This dataset consists of 21211 records. Only records containing data between ground level and 25,000 m above sea level, and with more than 8 data in this altitude range, were selected for this work, resulting in a dataset containing 14,976 records.
Polar diagrams of wind direction and intensity recorded at the Brindisi station (Aeronautica Militare), about 300 km from Vesuvius, in the period between 1962 and 1976. In the analysis we consider data in the period between 1962 and 1976 (3125 records). Please note that, differently from the convention usually adopted in meteorology, these plots do not represent the direction of provenance of the wind, but the direction towards which the wind blows
Figure 2 shows the polar diagrams of the winds, averaged on different horizontal layers 5 km thick, for the LIBR, LIRE, NOAA1, NOAA2 and ERAI datasets used in the present work. Again, here, the plots show the direction towards which the wind blows. A comparison between the plots shows no major differences among the datasets.
Wind direction and speed show a moderate seasonal variation, similar to that described by Costa et al. (2009). A statistics on this effect is reported in Fig. 3 for the ECMWF ERA-Interim (ERAI) dataset. The other datasets do not show significant difference with respect the ECMWF dataset.
The total erupted mass is the most important input parameter having a first order effect. The estimation of the total mass is commonly derived by field data analysis (Bonadonna and Costa 2012; Bonadonna and Houghton 2005; Pyle 1989) or by best-fit procedure (Bonasia et al. 2010; 2012; Connor and Connor 2006; Costa et al. 2009; Pfeiffer et al. 2005; Scollo et al. 2008). Keeping fixed the other parameters, the effect of the total mass on the deposit loading, neglecting or simplifying other non-linear phenomena (e.g. the aggregation process, the variability of the wind, etc.), follows a linear relationship.
A total mass of 2 1011 kg was used by Barberi et al. (1990) for assessing volcanic hazard from sub-Plinian eruptions at Vesuvius. This value derived from an estimate of the mass erupted during the main fallout phase of the Vesuvius A.D. 1631 sub-Plinian eruption obtained by Rosi et al. (1993), who found values in the range 4.5 1010 and 1.35 1011 kg, using the method of deposit thinning. However this kind of estimation can have a large uncertainty, depending on the quality and spatial distribution of the data, up to a factor of four (Bonadonna and Costa 2012).
In this paper we use, as reference, the same total mass adopted by Cioni et al. (2003), who proposed a value of 5 1011 kg as representative of sub-Plinian eruptions at Vesuvius. This value was also adopted by the Italian Department of Civil Protection for the reference scenario at Vesuvius (DPC 2012). Moreover, we performed a sensitivity analysis on the extension of the hazard zone by varying the total mass between 1011 and 1012 kg, in agreement with the corresponding volume range proposed by Cioni et al. (2008) for the sub-Plinian (type I) eruptions.
An analysis of particle types characterizing the TGSD of the Vesuvius 79 A.D. eruption can be found in Macedonio et al. (1988). Lithic particles usually show constant density whereas juvenile particles show a density variable with the size (e.g. lower density for particles of greater size) due to the presence of larger bubbles inside larger clasts. The variation of the density of juvenile particles with their size was analyzed by different authors (e.g. Bonadonna and Phillips 2003; Walker 1971). For sake of consistency, we used the same density variation as Macedonio et al. (1990), who adopted a simplified function for describing the particle density variations. However, as shown by Pfeiffer et al. (2005), particle density has a second order effect on the settling velocity with respect to the particle dimensions.
The evaluation of particle terminal settling velocity from the Total Grain Size Distribution (TGSD) is not trivial because of the irregular and variable shape of the real pyroclast particles (Armienti et al. 1988; Pfeiffer et al. 2005; Wilson and Huang 1979). Usually, it is obtained combining both theoretical models and experimental measurements (e.g. Arastoopour et al. 1982; Bagheri et al. 2015; Dellino et al. 2005; Ganser 1993; Wilson and Huang 1979). The terminal settling velocity depends on particle size, shape, and density and plays an important role on tephra dispersal (e.g. Pfeiffer et al. 2005). Normally, the densest and greatest particles fall near the volcanic vent, less dense and finest particles are instead dispersed farther from the volcanic vent. This natural trend is not valid in presence of significant volcanic ash aggregation, as discussed above.
It is worth noting that, due to the variation of the air density and viscosity with the altitude, the particle settling velocity has different values at different heights and changes, therefore, during the particle fall. For this reason, in this work, the terminal velocity is evaluated at each height by adopting the same approach used by Pfeiffer et al. (2005), who considered an empirical settling velocity-altitude relationship for each particle size.
The effective horizontal atmospheric diffusion coefficient is an empirical parameter adopted in plume models, that describes the effective spreading due to atmospheric turbulence and gravitational spreading at the Neutral Buoyancy Level (Costa et al. 2013). It is usually obtained by best-fit procedure (e.g. Bonadonna et al. 2002; Bonasia et al. 2010; Macedonio et al. 1988) and typically ranges between 102 and 104 m2/s. Studies on Vesuvius eruptions have shown that the horizontal diffusion coefficient, for sub-Plinian and Plinian eruptions, for the target area, ranges between 1000 and 10,000 m2/s (Bonasia et al. 2010; Macedonio et al. 1988).
7fc3f7cf58