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.

Development of an Inverse Plume Model for Mass Eruption Rate in Unsteady Conditions

Following volcanic eruptions, forecasters need accurate estimates of mass eruption rate (MER) to appropriately predict the downstream effects. Most analyses use simple correlations or models based on large eruptions at steady conditions, even though many volcanoes feature significant unsteadiness. To address this, a superposition model is developed based on a technique used for spray injection applications, which predicts plume height as a function of the time-varying exit velocity. This model can be inverted, providing estimates of MER using field observations of a plume. The model parameters are optimized using laboratory data for plumes with physically-relevant exit profiles and Reynolds numbers, resulting in predictions that agree to within 10% of measured exit velocities. The model performance is examined using a historic eruption from Stromboli with well-documented unsteadiness, again providing MER estimates of the correct order of magnitude. This method can provide a rapid alternative for real-time forecasting of small, unsteady eruptions.

Modelling pyroclastic density currents from a subplinian eruption at La Soufrière de Guadeloupe (West Indies, France)

We have used a three-dimensional, non-equilibrium multiphase flow numerical model to simulate subplinian eruption scenarios at La Soufrière de Guadeloupe (Lesser Antilles, France). Initial and boundary conditions for computer simulations were set on the basis of independent estimates of eruption source parameters (i.e. mass eruption rate, volatile content, temperature, grain size distribution) from a field reconstruction of the 1530 CE subplinian eruption. This event is here taken as a reference scenario for hazard assessment at La Soufrière de Guadeloupe. A parametric study on eruption source parameters allowed us to quantify their influence on the simulated dynamics and, in particular, the increase of the percentage of column collapse and pyroclastic density current (PDC) intensity, at constant mass eruption rate, with variable vent diameter. Numerical results enabled us to quantify the effects of the proximal morphology on distributing the collapsing mass around the volcano and into deep and long valleys and to estimate the areas invaded by PDCs, their associated temperature and dynamic pressure. Significant impact (temperature > 300 °C and dynamic pressure > 1 kPa) in the inhabited region around the volcano is expected for fully collapsing conditions and mass eruption rates > 2 × 107 kg/s. We thus combine this spatial distribution of temperature and dynamic pressure with an objective consideration of model-related uncertainty to produce preliminary PDC hazard maps for the reference scenario. In such a representation, we identify three areas of varying degree of susceptibility to invasion by PDCs—very likely to be invaded (and highly impacted), susceptible to invasion (and moderately impacted), and unlikely to be invaded (or marginally impacted). The study also raises some key questions about the use of deterministic scenario simulations for hazard assessment, where probability distributions and uncertainties are difficult to estimate. Use of high-performance computing techniques will in part allow us to overcome such difficulties, but the problem remains open in a scientific context where validation of numerical models is still, necessarily, an incomplete and ongoing process. Nevertheless, our findings provide an important contribution to the quantitative assessment of volcanic hazard and risk at La Soufrière de Guadeloupe particularly in the context of the current unrest of the volcano and the need to prepare for a possible future reawakening of the volcano that could culminate in a magmatic explosive eruption.

A new database of independently estimated eruption source parameters devoted to eruptive column model evaluation

<p>Eruptive column models are crucial for managing volcanic crises, forecasting future events, and reconstructing past eruptions. Given their central role in volcanology and the large uncertainties weakening their predictions, the evaluation and improvement of these models is critical. Such evaluation is challenging as it requires independent estimates of the main model inputs (e.g. mass eruption rate) and outputs (e.g. column height). Despite recent efforts to extend datasets of independently estimated eruption source parameters (ESP) (e.g. Mastin 2014, Aubry et al. 2017), there is no standardized, maintained, and community-based ESP database devoted to the evaluation of eruptive column models.</p><p>Here we present a new ESP database designed to respond to the needs of the plume modelling community, and which will also be valuable to observatories, field volcanologists, and volcanic ash advisory centers. We compiled data for over 130 eruptive events with independent estimates of: i) the mass eruption rate; ii) the height reached by the column; and iii) atmospheric conditions during the eruption. In contrast with previous ESP datasets, we distinguish estimates of column height that relate to different phases (ash and SO2) and parts of the column (plume top or umbrella). We additionally provide the total grain size distribution, uncertainties in eruption parameters, and multiple sources for atmospheric profiles for events where these parameters are available. The database also includes a wealth of additional information which will enable modelers to distinguish between different eruptions when evaluating or calibrating models. This includes the type of eruption, the morphology of the plume (weak/transitional/strong), and the occurrence and mass entrained within pyroclastic density currents.</p><p>We will apply the new database to revisit empirical scaling relationships between the mass eruption rate and “plume height”. In particular, we will show how such relationships depend on the type of height (e.g. SO2 height vs. ash top height) and eruption (e.g. magmatic vs. phreatomagmatic) considered. We will also discuss the difficulties and limitations of compiling ESP estimates from the literature as well as characterizing fundamentally unsteady volcanic events by a single value for each ESP.</p>