Zircon of Triassic Age in the Stuttgart Formation (Schilfsandstein)—Witness of Tephra Fallout in the Central European Basin and New Constraints on the Mid-Carnian Episode

The Carnian Stuttgart-Formation (Schilfsandstein) of the Central European Basin contains relics of Triassic volcanic detritus in form of euhedral zircon grains and authigenic analcime. Multiple LA-ICP-MS spot analyses of single zircon crystals from an outcrop near Heilbronn (SW Germany) yielded weighted average 206Pb/238U ages between 250 and 230 Ma, providing evidence for tephra fallout in the southern part of the Central European Basin related to Olenekian, Anisian–Ladinian and Carnian volcanic activity. The tephra was probably transported by monsoonal circulations from volcanic centres of the NW Tethys to the Central European Basin. The four youngest zircon crystals gave a weighted average 206Pb/238U age of 231.1 ± 1.6 Ma (10 analyses), which is interpreted to date syn-depositional tephra fallout into the fluvial Lower Schilfsandstein Member of the Stuttgart Formation. This new maximum depositional age provides the first evidence that deposition of the Stuttgart Formation, which represents the type-example of the mid-Carnian episode, a global episode of enhanced flux of siliciclastic detritus and related environmental perturbations, occurred during the Tuvalian 2 substage at ca. 231 Ma, about 3 million years later than suggested by previous correlations. Zircon grains with weighted average 206Pb/238U ages of 236.0 ± 1.2 Ma (n = 17) and 238.6 ± 1.5 Ma (n = 6) and 206Pb/238U ages between 241 ± 6 and 250 ± 3 Ma point to the presence of tephra in early Carnian to Olenekian strata of the Keuper to Buntsandstein Groups. Traces of these reworked tephra were incorporated into the Stuttgart Formation due to fluvial erosion in the southern Central European Basin and at its margins.

Lahar risk assessment from source identification to potential impact analysis: the case of Vulcano Island, Italy

Lahars are rapid flows composed of water and volcaniclastic sediments, which have the potential to impact residential buildings and critical infrastructure as well as to disrupt critical services, especially in the absence of hazard-based land-use planning. Their destructive power is mostly associated with their velocity (related to internal flow properties and topographic interactions) and to their ability to bury buildings and structures (due to deposit thickness). The distance reached by lahars depends on their volume, on sediments/water ratio, as well as on the geometrical properties of the topography where they propagate. Here we present the assessment of risk associated with lahar using Vulcano island (Italy) as a case study. First, we estimated an initial lahar source volume considering the remobilisation by intense rain events of the tephra fallout on the slopes of the La Fossa cone (the active system on the island), where the tephra fallout is associated with the most likely scenario (e.g. long-lasting Vulcanian cycle). Second, we modelled and identified the potential syn-eruptive lahar impact areas on the northern sector of Vulcano, where residential and touristic facilities are located. We tested a range of parameters (e.g., entrainment capability, consolidation of tephra fallout deposit, friction angle) that can influence lahar propagation output both in terms of intensity of the event and extent of the inundation area. Finally, exposure and vulnerability surveys were carried out in order to compile exposure and risk maps for lahar-flow front velocity (semi-quantitative indicator-based risk assessment) and final lahar-deposit thickness (qualitative exposure-based risk assessment). Main outcomes show that the syn-eruptive lahar scenario with medium entrainment capability produces the highest impact associated with building burial by the final lahar deposit. Nonetheless, the syn-eruptive lahar scenario with low entrainment capacity is associated with higher runout and results in the highest impact associated with lahar-flow velocities. Based on our simulations, two critical infrastructures (telecommunication and power plant), as well as the main road crossing the island are exposed to potential lahar impacts (either due to lahar-flow velocity or lahar-deposit thickness or both). These results show that a risk-based spatial planning of the island could represent a valuable strategy to reduce the volcanic risk in the long term.

Volcanic hazard assessment for tephra fallout in Martinique

Mount Pelée (Martinique) is one of the most active volcanoes in the Lesser Antilles arc with more than 34 magmatic events in the last 24,000 years, including the deadliest eruption of the 20th century. The current volcanic hazard map used in the civil security plan puts the emphasis on the volcanic hazard close to the volcano. This map is however based on an incomplete eruptive history and does not take into account the variability of the expected source conditions (mass eruption rate, total erupted mass, and grain-size distribution) or the wind effect on ash dispersal. We propose here to refine the volcanic hazard map for tephra fallout by using the 2-D model of ash dispersal HAZMAP. We first simulate the maximum expected eruptive scenario at Mount Pelée (i.e., the P3 eruption) using a seasonal wind profile. Building upon the good agreement with field data, we compute probability maps based on this maximum expected scenario, which show that tephra fallout hazard could threaten not only areas close to the volcano but also the southern part of Martinique. We then use a comprehensive approach based on 16 eruptive scenarios that include new field constraints obtained in the recent years on the past Plinian eruptions of Mount Pelée volcano. Each eruptive scenario considers different values of total erupted mass and mass eruption rate, and is characterized by a given probability of occurrence estimated from the refined eruptive history of the volcano. The 1979-2019 meteorological ERA-5 database is used to further take into account the daily variability of winds. These new probability maps show that the area of probable total destruction is wider when considering the 16 scenarios compared to the maximum expected scenario. The southern part of Martinique, although less threatened than when considering the maximum expected scenario, would still be impacted both by tephra fallout and by its high dependence on the water and electrical network carried from the northern part of the island. Finally, we show that key infrastructures in Martinique (such as the international airport) have a non-negligible probability of being impacted by a future Plinian eruption of the Mount Pelée. These results provide strong arguments for and will support significant and timely reconceiving of the emergency procedures as the local authorities have now placed Mount Pelée volcano on alert level yellow (vigilance) based on increased seismicity and tremor-type signals.

FALL3D-8.0: a computational model for atmospheric transport and deposition of particles, aerosols and radionuclides – Part 2: Model validation

This paper presents model validation results for the latest version release of the FALL3D atmospheric transport model. The code has been redesigned from scratch to incorporate different categories of species and to overcome legacy issues that precluded its preparation towards extreme-scale computing. The model validation is based on the new FALL3D-8.0 test suite, which comprises a set of four real case studies that encapsulate the major features of the model; namely, the simulation of long-range fine volcanic ash dispersal, volcanic SO2 dispersal, tephra fallout deposits and the dispersal and deposition of radionuclides. The first two test suite cases (i.e. the June 2011 Puyehue-Cordón Caulle ash cloud and the June 2019 Raikoke SO2 cloud) are validated against geostationary satellite retrievals and demonstrate the new FALL3D data insertion scheme. The metrics used to validate the volcanic ash and SO2 simulations are the structure, amplitude and location (SAL) metric and the figure of merit in space (FMS). The other two test suite cases (i.e. the February 2013 Mt. Etna ash cloud and associated tephra fallout deposit, and the dispersal of radionuclides resulting from the 1986 Chernobyl nuclear accident) are validated with scattered ground-based observations of deposit load and local particle grain size distributions and with measurements from the Radioactivity Environmental Monitoring database. For validation of tephra deposit loads and radionuclides, we use two variants of the normalised root-mean-square error metric. We find that FALL3D-8.0 simulations initialised with data insertion consistently improve agreement with satellite retrievals at all lead times up to 48 h for both volcanic ash and SO2 simulations. In general, SAL scores lower than 1.5 and FMS scores greater than 0.40 indicate acceptable agreement with satellite retrievals of volcanic ash and SO2. In addition, we show very good agreement, across several orders of magnitude, between the model and observations for the 2013 Mt. Etna and 1986 Chernobyl case studies. Our results, along with the validation datasets provided in the publicly available test suite, form the basis for future improvements to FALL3D (version 8 or later) and also allow for model intercomparison studies.

Development of the muographic tephra deposit monitoring system

Measurements of volcanic tephra fallout deposits provide useful information about the magnitude and intensity of explosive volcanic eruptions and potential for remobilization of deposits as dangerous volcanic flows. However, gathering information in the vicinity of erupting craters is extremely dangerous, and moreover, it is often quite difficult to determine deposit thickness proximal to volcanic craters because the thickness of the deposit is too great to easily measure; thus, airborne remote sensing technologies have generally been utilized during the intermission between eruptions. As an alternative tool, a muographic tephra deposit monitoring system was developed in this work. Here we report the performance of this system by applying the muographic data acquired at Sakurajima volcano, Japan as an example. By assuming the average density of the deposit was 2.0 g cm−3, the deposit thicknesses measured with muography were in agreement with the airborne results, indicating that volcanic fallout built up within the upper river basin, showed its potential for monitoring the episodic tephra fallouts even during eruptions.

Near-Real-Time Tephra Fallout Assessment at Mt. Etna, Italy

During explosive eruptions, emergency responders and government agencies need to make fast decisions that should be based on an accurate forecast of tephra dispersal and assessment of the expected impact. Here, we propose a new operational tephra fallout monitoring and forecasting system based on quantitative volcanological observations and modelling. The new system runs at the Istituto Nazionale di Geofisica e Vulcanologia, Osservatorio Etneo (INGV-OE) and is able to provide a reliable hazard assessment to the National Department of Civil Protection (DPC) during explosive eruptions. The new operational system combines data from low-cost calibrated visible cameras and satellite images to estimate the variation of column height with time and model volcanic plume and fallout in near-real-time (NRT). The new system has three main objectives: (i) to determine column height in NRT using multiple sensors (calibrated cameras and satellite images); (ii) to compute isomass and isopleth maps of tephra deposits in NRT; (iii) to help the DPC to best select the eruption scenarios run daily by INGV-OE every three hours. A particular novel feature of the new system is the computation of an isopleth map, which helps to identify the region of sedimentation of large clasts (≥5 cm) that could cause injuries to tourists, hikers, guides, and scientists, as well as damage buildings in the proximity of the summit craters. The proposed system could be easily adapted to other volcano observatories worldwide.