Search

In this study we present a novel general methodology for Probabilistic Volcanic Hazard Assessment (PVHA) for lahars. We apply the methodology to perform a probabilistic assessment in the Campanian Plain (Southern Italy), focussing on syn-eruptive lahars from a reference-size eruption from Somma-Vesuvius. We take advantage of new field data relative to volcaniclastic flow deposits in the target region and recent improvements in modelling lahars. We explore the effect of uncertainty on the flow initial conditions on the invasion of lahar in the target area, by sampling coherent sets of values for the input model parameters and performing a large number of simulations. We analyse the simulation outputs to produce hazard curves, hazard maps and probability maps for the maximum flow thickness, and hazard surface and probability maps for joint thresholds in flow thickness and dynamic pressure. We believe the latter hazard products represent, on one hand, a novel product in PVHA for lahars, and, on the other hand, a useful means for impact assessment.

We present the results of a meso-scale systematic structural analysis of fractures, faults and dykes exposed at the Somma-Vesuvius volcano (southern Italy). Observed fractures include: (i) radial and tangential (with respect the caldera axis), sub-metric to metric joints associated with the edifice load and volcano-tectonic activity (i.e. inflation, deflation and caldera collapse stages) and (ii) decameter-scale fractures related to volcano flank instabilities. For the Somma-Vesuvius volcano, preexisting radial joints were commonly reactivated as transfer faults during the caldera formation, allowing different blocks to move toward the center of the collapsing area. Dykes occur with different geometries, including en-echelon structures bounding structural depressions. The orientation analysis of all structures indicates that they are preferentially oriented. Furthermore, we provide a morphological lineament analysis using high-resolution Digital Terrain Models of Somma-Vesuvius. Azimuth and spatial distribution of dykes and morphological lineaments were analyzed for comparison with the old Somma Crater and Gran Cono axes, respectively. Results highlight the overprinting of radial and clustered strain patterns recorded in different volcano-tectonic evolution stages. We suggest a possible deformation evolution model in which structures develop along either radial or preferential trends, highlighting different volcanic conditions: (i) where radial patterns occur, the structures developed during volcanic inflation cycles with a closed magmatic conduit condition whereas (ii) clustered patterns are probably associated with a regional strain field that overcomes the local deformation field, a situation typical in the case of open-conduit activity.

This study presents a new method for producing long-term hazard maps for pyroclastic density currents (PDC) originating at Campi Flegrei caldera. Such method is based on a doubly stochastic approach and is able to combine the uncertainty assessments on the spatial location of the volcanic vent, the size of the flow and the expected time of such an event. The results are obtained by using a Monte Carlo approach and adopting a simplified invasion model based on the box model integral approximation. Temporal assessments are modelled through a Cox-type process including self-excitement effects, based on the eruptive record of the last 15 kyr. Mean and percentile maps of PDC invasion probability are produced, exploring their sensitivity to some sources of uncertainty and to the effects of the dependence between PDC scales and the caldera sector where they originated. Conditional maps representative of PDC originating inside limited zones of the caldera, or of PDC with a limited range of scales are also produced. Finally, the effect of assuming different time windows for the hazard estimates is explored, also including the potential occurrence of a sequence of multiple events. Assuming that the last eruption of Monte Nuovo (A.D. 1538) marked the beginning of a new epoch of activity similar to the previous ones, results of the statistical analysis indicate a mean probability of PDC invasion above 5% in the next 50 years on almost the entire caldera (with a probability peak of ~25% in the central part of the caldera). In contrast, probability values reduce by a factor of about 3 if the entire eruptive record is considered over the last 15 kyr, i.e. including both eruptive epochs and quiescent periods.

Campi Flegrei (CF) is an example of an active caldera containing densely populated settlements at very high risk of pyroclastic density currents (PDCs). We present here an innovative method for assessing background spatial PDC hazard in a caldera setting with probabilistic invasion maps conditional on the occurrence of an explosive event. The method encompasses the probabilistic assessment of potential vent opening positions, derived in the companion paper, combined with inferences about the spatial density distribution of PDC invasion areas from a simplified flow model, informed by reconstruction of deposits from eruptions in the last 15 ka. The flow model describes the PDC kinematics and accounts for main effects of topography on flow propagation. Structured expert elicitation is used to incorporate certain sources of epistemic uncertainty, and a Monte Carlo approach is adopted to produce a set of probabilistic hazard maps for the whole CF area. Our findings show that, in case of eruption, almost the entire caldera is exposed to invasion with a mean probability of at least 5%, with peaks greater than 50% in some central areas. Some areas outside the caldera are also exposed to this danger, with mean probabilities of invasion of the order of 5-10%. Our analysis suggests that these probability estimates have location-specific uncertainties which can be substantial. The results prove to be robust with respect to alternative elicitation models and allow the influence on hazard mapping of different sources of uncertainty, and of theoretical and numerical assumptions, to be quantified.

Campi Flegrei is an active volcanic area situated in the Campanian Plain (Italy) and dominated by a resurgent caldera. The great majority of past eruptions have been explosive, variable in magnitude, intensity, and in their vent locations. In this hazard assessment study we present a probabilistic analysis using a variety of volcanological data sets to map the background spatial probability of vent opening conditional on the occurrence of an event in the foreseeable future. The analysis focuses on the reconstruction of the location of past eruptive vents in the last 15 ka, including the distribution of faults and surface fractures as being representative of areas of crustal weakness. One of our key objectives was to incorporate some of the main sources of epistemic uncertainty about the volcanic system through a structured expert elicitation, thereby quantifying uncertainties for certain important model parameters and allowing outcomes from different expert weighting models to be evaluated. Results indicate that past vent locations are the most informative factors governing the probabilities of vent opening, followed by the locations of faults and then fractures. Our vent opening probability maps highlight the presence of a sizeable region in the central eastern part of the caldera where the likelihood of new vent opening per kilometer squared is about 6 times higher than the baseline value for the whole caldera. While these probability values have substantial uncertainties associated with them, our findings provide a rational basis for hazard mapping of the next eruption at Campi Flegrei caldera.