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.

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.

No abstract

We provide a database of the coseismic geological surface effects following the Mw 6.5 Norcia earthquake that hit central Italy on 30 October 2016. This was one of the strongest seismic events to occur in Europe in the past thirty years, causing complex surface ruptures over an area of >400 km2. The database originated from the collaboration of several European teams (Open EMERGEO Working Group; about 130 researchers) coordinated by the Istituto Nazionale di Geofisica e Vulcanologia. The observations were collected by performing detailed field surveys in the epicentral region in order to describe the geometry and kinematics of surface faulting, and subsequently of landslides and other secondary coseismic effects. The resulting database consists of homogeneous georeferenced records identifying 7323 observation points, each of which contains 18 numeric and string fields of relevant information. This database will impact future earthquake studies focused on modelling of the seismic processes in active extensional settings, updating probabilistic estimates of slip distribution, and assessing the hazard of surface faulting.

Morphostructural data derived from Lidar (Light detection and ranging) surveys carried out on Mount Etna in 2005 and 2007 are compared with earlier aerophotogrammetric surveys in 1986 and 1998. These data render an unprecedentedly clear and quantitative image of morphostructural and volumetric changes that have affected the summit area of the volcano in the past two decades and permit the production of a new topographic map. The computed volume gain during the 1986–2007 period amounts to 112 ± 12 × 10 6 m 3, at a mean annual rate of ∼5.3 × 10 6 m 3. The comparison of the various surveys furthermore emphasizes the levels of accuracy and resolution of the different techniques applied. The Lidar technology used in 2007 allows production of high‐precision maps in near‐real‐time, facilitating work concerning environmental hazards such as numerical simulations of, e.g., lava flows.

Volcanic islands pose several major types of natural hazards, often interconnected and concentrated in relatively small areas. The quantification of these hazards must be framed from a multi-hazard perspective whilst building on existing single-hazard analyses. Ischia is a densely inhabited volcanic island with a long eruptive history lasting more than 150 ka (last in 1302 AD) characterized by the significant asymmetric resurgence of a caldera block. Here, we review the state-of-art of the natural hazards of Ischia, aiming at building a solid base for future holistic multihazard quantifications. We frame our analysis in three steps: i) review of geological, historical and current activity; ii) review of available hazard models and analyses; iii) development of an interpretative framework for the interdependent hazards. The results highlight that volcanic activity has been quite intense and many volcanorelated hazardous phenomena have affected the island including in very recent times, both for eruptive (phreatic or magmatic eruptions) and non-eruptive (earthquakes, landslides, and tsunamis) phenomena. The effects of some of them (e.g. tsunamis, tephra) are also relevant beyond the island territory. Quantitative hazard assessments are almost absent and should be developed in the future considering the evident interconnections between hazards. To this end, we propose a conceptual interpretative multi-hazard framework that highlights the fundamental role played by the resurgent block in controlling and connecting the different hazards, in terms of both spatial distribution of the sources and temporal clustering.

This study presents new and revised data sets about the spatial distribution of past volcanic vents, eruptive fissures, and regional/local structures of the Somma‐Vesuvio volcanic system (Italy). The innovative features of the study are the identification and quantification of important sources of uncertainty affecting interpretations of the data sets. In this regard, the spatial uncertainty of each feature is modeled by an uncertainty area, i.e., a geometric element typically represented by a polygon drawn around points or lines. The new data sets have been assembled as an updatable geodatabase that integrates and complements existing databases for Somma‐Vesuvio. The data are organized into 4 data sets and stored as 11 feature classes (points and lines for feature locations and polygons for the associated uncertainty areas), totaling more than 1700 elements. More specifically, volcanic vent and eruptive fissure elements are subdivided into feature classes according to their associated eruptive styles: (i) Plinian and sub‐Plinian eruptions (i.e., large‐ or medium‐scale explosive activity); (ii) violent Strombolian and continuous ash emission eruptions (i.e., small‐scale explosive activity); and (iii) effusive eruptions (including eruptions from both parasitic vents and eruptive fissures). Regional and local structures (i.e., deep faults) are represented as linear feature classes. To support interpretation of the eruption data, additional data sets are provided for Somma‐Vesuvio geological units and caldera morphological features. In the companion paper, the data presented here, and the associated uncertainties, are used to develop a first vent opening probability map for the Somma‐Vesuvio caldera, with specific attention focused on large or medium explosive events.

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 modeled 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.