Suchergebnisse

Abstract. The study of volcanic flow hazards in a probabilistic framework centers around systematic experimental numerical modeling of the hazardous phenomenon and the subsequent generation and interpretation of a probabilistic hazard map (PHM). For a given volcanic flow (e.g., lava flow, lahar, pyroclastic flow, ash cloud), the PHM is typically interpreted as the point-wise probability of inundation by flow material. In the current work, we present new methods for calculating spatial representations of the mean, standard deviation, median, and modal locations of the hazard's boundary as ensembles of many deterministic runs of a physical model. By formalizing its generation and properties, we show that a PHM may be used to construct these statistical measures of the hazard boundary which have been unrecognized in previous probabilistic hazard analyses. Our formalism shows that a typical PHM for a volcanic flow not only gives the point-wise inundation probability, but also represents a set of cumulative distribution functions for the location of the inundation boundary with a corresponding set of probability density functions. These distributions run over curves of steepest probability gradient ascent on the PHM. Consequently, 2-D space curves can be constructed on the map which represents the mean, median, and modal locations of the likely inundation boundary. These curves give well-defined answers to the question of the likely boundary location of the area impacted by the hazard. Additionally, methods of calculation for higher moments including the standard deviation are presented, which take the form of map regions surrounding the mean boundary location. These measures of central tendency and variance add significant value to spatial probabilistic hazard analyses, giving a new statistical description of the probability distributions underlying PHMs. The theory presented here may be used to aid construction of improved hazard maps, which could prove useful for planning and emergency management purposes. This formalism also allows for application to simplified processes describable by analytic solutions. In that context, the connection between the PHM, its moments, and the underlying parameter variation is explicit, allowing for better source parameter estimation from natural data, yielding insights about natural controls on those parameters.

Abstract. We describe a method for calculating the probability that a distal geographic location is impacted by a pyroclastic density current (PDC) of a given size, considering the key related uncertainties. Specifically, we evaluate the minimum volume and mass of a PDC generated at the Aso caldera (Japan) that might affect each of five distal infrastructure (marker) sites, with model input parameter uncertainties derived from expert judgment. The 5 marker sites are all located 115–145 km from the caldera; as these lie in well-separated directions, we can test the effects of the different topographic shielding effects in each case. To inform our probabilistic analysis, we apply alternative kinetic energy assessment approaches, i.e., rock avalanche and density current dynamics. In the latter formulation, the minimum mass needed to reach the markers ranges between median values of ∼153×1012 and ∼465×1012 kg (M 7.2–7.7), depending on the site. Rock avalanche dynamics modeling indicates that a ∼3-times greater mass would be required to reach the marker sites with 50 % probability, while the hypothetical scenario of a relatively dilute distal ash cloud would require ∼3-times less mass. We compare our results with the largest recorded Aso eruption, showing that a catastrophic eruption, similar to Aso-4, ≈ M8, would present a significant conditional probability of PDCs reaching the marker sites, in the density current formulation and contingent on uncertainty in the erupted mass and on marker site direction.

Abstract. The San Salvador volcanic complex (El Salvador) and Nejapa-Chiltepe volcanic complex (Nicaragua) have been characterized by a significant variability in eruption style and vent location. Densely inhabited cities are built on them and their surroundings, including the metropolitan areas of San Salvador (∼2.4 million people) and Managua (∼1.4 million people), respectively. In this study we present novel vent opening probability maps for these volcanic complexes, which are based on a multi-model approach that relies on kernel density estimators. In particular, we present thematic vent opening maps, i.e., we consider different hazardous phenomena separately, including lava emission, small-scale pyroclastic density currents, ejection of ballistic projectiles, and low-intensity pyroclastic fallout. Our volcanological dataset includes: (1) the location of past vents, (2) the mapping of the main fault structures, and (3) the eruption styles of past events, obtained from critical analysis of the literature and/or inferred from volcanic deposits and morphological features observed remotely and in the field. To illustrate the effects of considering the expected eruption style in the construction of vent opening maps, we focus on the analysis of small-scale pyroclastic density currents derived from phreatomagmatic activity or from low-intensity magmatic volcanism. For the numerical simulation of these phenomena we adopted the recently developed branching energy cone model by using the program ECMapProb. Our results show that the implementation of thematic vent opening maps can produce significantly different hazard levels from those estimated with traditional, non-thematic maps.

Abstract. We detail a new prediction-oriented procedure aimed at volcanic hazard
assessment based on geophysical mass flow models constrained with
heterogeneous and poorly defined data. Our method relies on an itemized
application of the empirical falsification principle over an arbitrarily wide
envelope of possible input conditions. We thus provide a first step towards a
objective and partially automated experimental design construction. In
particular, instead of fully calibrating model inputs on past observations,
we create and explore more general requirements of consistency, and then we
separately use each piece of empirical data to remove those input values that
are not compatible with it. Hence, partial solutions are defined to the inverse
problem. This has several advantages compared to a traditionally posed
inverse problem: (i) the potentially nonempty inverse images of partial
solutions of multiple possible forward models characterize the solutions to
the inverse problem; (ii) the partial solutions can provide hazard estimates
under weaker constraints, potentially including extreme cases that are
important for hazard analysis; (iii) if multiple models are applicable,
specific performance scores against each piece of empirical information can
be calculated. We apply our procedure to the case study of the Atenquique
volcaniclastic debris flow, which occurred on the flanks of Nevado de Colima
volcano (Mexico), 1955. We adopt and compare three depth-averaged models
currently implemented in the TITAN2D solver, available from https://vhub.org
(Version 4.0.0 – last access: 23 June 2016). The associated inverse problem
is not well-posed if approached in a traditional way. We show that our procedure
can extract valuable information for hazard assessment, allowing the exploration
of the impact of synthetic flows that are similar to those that occurred in the
past but different in plausible ways. The implementation of multiple models is
thus a crucial aspect of our approach, as they can allow the covering of other
plausible flows. We also observe that model selection is inherently linked to
the inversion problem.

Formalised elicitation of expert judgements has been used to help tackle several problematic societal issues, including volcanic crises and pandemic threats. We present an expert elicitation exercise for Piton de la Fournaise volcano, La Réunion island, held remotely in April 2021. This involved 28 experts from nine countries who considered a hypothetical effusive eruption crisis involving a new vent opening in a high-risk area. The tele-elicitation presented several challenges, but is a promising and workable option for application to future volcanic crises. Our exercise considered an "uncommon" eruptive scenario with a vent outside the present caldera and within inhabited areas, and provided uncertainty ranges for several hazard-related questions for such a scenario (e.g. probability of eruption within a defined timeframe; elapsed time until lava flow reaches a critical location, and other hazard management issues). Our exercise indicated that such a scenario would probably present very different characteristics compared to recent eruptions, and that it is fundamental to include well-prepared expert elicitations in updated civil protection evacuation plans to improve disaster response procedures. ; This work was funded by the Agence National de la Recherche (ANR) through project Lava Advance into Vulnerable Areas (LAVA; ANR program: DS0902 2016; project: ANR-16CE39–0009). This is ANR-LAVA con- tribution n°21. This contribution is part of the European Commis- sion grant EVE (DG ECHO Ref: 826292). J.M. received funding from the IMAGINE ERC Grant No 804162 to support the development of this paper. A.T. was funded by the ClerVolc project - Programme 1 "Detection and characterization of volcanic plumes and ash clouds" funded by the French government 'Laboratory of Excellence' initiative. This is ClerVolc contribution n°532. A.T. was also partially funded by the by the French government IDEX-ISITE initiative 16- IDEX-0001 (CAP 20-25). ; Peer reviewed