PYFLOW_2.0: a computer program for calculating flow properties and impact parameters of past dilute pyroclastic density currents based on field data

This paper presents PYFLOW_2.0, a hazard tool for the calculation of the impact parameters of dilute pyroclastic density currents (DPDCs). DPDCs represent the dilute turbulent type of gravity flows that occur during explosive volcanic eruptions; their hazard is the result of their mobility and the capability to laterally impact buildings and infrastructures and to transport variable amounts of volcanic ash along the path. Starting from data coming from the analysis of deposits formed by DPDCs, PYFLOW_2.0 calculates the flow properties (e.g., velocity, bulk density, thickness) and impact parameters (dynamic pressure, deposition time) at the location of the sampled outcrop. Given the inherent uncertainties related to sampling, laboratory analyses, and modeling assumptions, the program provides ranges of variations and probability density functions of the impact parameters rather than single specific values; from these functions, the user can interrogate the program to obtain the value of the computed impact parameter at any specified exceedance probability. In this paper, the sedimentological models implemented in PYFLOW_2.0 are presented, program functionalities are briefly introduced, and two application examples are discussed so as to show the capabilities of the software in quantifying the impact of the analyzed DPDCs in terms of dynamic pressure, volcanic ash concentration, and residence time in the atmosphere. The software and user’s manual are made available as a downloadable electronic supplement

Understanding the dynamics of unsteady buoyant Jets:an experimental analogue of Vulcanian and Strombolian style eruptions

Explosive volcanic eruptions, which are characterized by the discharge of ash and gas from the vent into the atmosphere, are an example of a naturally occurring buoyant jet. These buoyant jets can significantly impact the surrounding environment; for example, the presence of fine ash particles in the atmosphere can damage aircraft engines, potentially leading to engine failure. Therefore, during an explosive eruption, volcanic ash advisory centers (VAACs) consistently monitor the concentration of ash in the atmosphere using numerical models. These numerical models require the definition of a source term (i.e., source mass eruption rate, plume height and total grain size distribution), which is often obtained from simpler one-dimensional models. One dimensional models derived from well-established theories successfully replicate the dynamics of the initial buoyant jet; however, they assume time-averaged source conditions which are not observed in field-scale vulcanian and strombolian style eruptions. As such, there is a disconnectbetween these well-established theories assuming time averaged source conditions and reality. This disconnect may introduce uncertainties in ash concentration forecasts, potentially resulting in practical implications such as unnecessary airspace closures or flights operating in hazardous conditions. The present contribution utilizes scaled laboratory experiments to quantify theinfluence of source variability on the dynamics of buoyant jets and evaluates potential deviations from time-average assumptions

Understanding the dynamics of unsteady buoyant Jets:an experimental analogue of Vulcanian and Strombolian style eruptions

Explosive volcanic eruptions, which are characterized by the discharge of ash and gas from the vent into the atmosphere, are an example of a naturally occurring buoyant jet. These buoyant jets can significantly impact the surrounding environment; for example, the presence of fine ash particles in the atmosphere can damage aircraft engines, potentially leading to engine failure. Therefore, during an explosive eruption, volcanic ash advisory centers (VAACs) consistently monitor the concentration of ash in the atmosphere using numerical models. These numerical models require the definition of a source term (i.e., source mass eruption rate, plume height and total grain size distribution), which is often obtained from simpler one-dimensional models. One dimensional models derived from well-established theories successfully replicate the dynamics of the initial buoyant jet; however, they assume time-averaged source conditions which are not observed in field-scale vulcanian and strombolian style eruptions. As such, there is a disconnectbetween these well-established theories assuming time averaged source conditions and reality. This disconnect may introduce uncertainties in ash concentration forecasts, potentially resulting in practical implications such as unnecessary airspace closures or flights operating in hazardous conditions. The present contribution utilizes scaled laboratory experiments to quantify theinfluence of source variability on the dynamics of buoyant jets and evaluates potential deviations from time-average assumptions

Understanding the dynamics of unsteady buoyant Jets:an experimental analogue of Vulcanian and Strombolian style eruptions

Explosive volcanic eruptions, which are characterized by the discharge of ash and gas from the vent into the atmosphere, are an example of a naturally occurring buoyant jet. These buoyant jets can significantly impact the surrounding environment; for example, the presence of fine ash particles in the atmosphere can damage aircraft engines, potentially leading to engine failure. Therefore, during an explosive eruption, volcanic ash advisory centers (VAACs) consistently monitor the concentration of ash in the atmosphere using numerical models. These numerical models require the definition of a source term (i.e., source mass eruption rate, plume height and total grain size distribution), which is often obtained from simpler one-dimensional models. One dimensional models derived from well-established theories successfully replicate the dynamics of the initial buoyant jet; however, they assume time-averaged source conditions which are not observed in field-scale vulcanian and strombolian style eruptions. As such, there is a disconnectbetween these well-established theories assuming time averaged source conditions and reality. This disconnect may introduce uncertainties in ash concentration forecasts, potentially resulting in practical implications such as unnecessary airspace closures or flights operating in hazardous conditions. The present contribution utilizes scaled laboratory experiments to quantify theinfluence of source variability on the dynamics of buoyant jets and evaluates potential deviations from time-average assumptions

The Impact of Eruption Source Parameter Uncertainties on Ash Dispersion Forecasts During Explosive Volcanic Eruptions

Volcanic ash in the atmosphere is a hazard to aviation. To predict which areas of airspace are most likely to be affected by the presence of ash, Volcanic Ash Advisory Centers (VAACs) use observations and atmospheric dispersion models. These models are initialized with, among other parameters, a mass eruption rate (MER), which quantifies the emission rate into the atmosphere at the source. This influences the predicted spatial–temporal evolution and concentration of the ash cloud. Different models are available to estimate MER from the volcanic plume height and some models also include the weather conditions (e.g., wind speed). The REFIR software tool uses time‐series of plume height estimated from observations and weather data to provide estimates of MER through time. Here we present an updated version of REFIR that can now be used also to calculate MER for past eruptions and produce output parameters in a format suitable for use with the NAME dispersion model (UK Met Office—London VAAC). We also investigate how uncertainty in input parameters is propagated through to dispersion model output. Our results show that a +/−1 km uncertainty on a 6 km high plume can result in the affected area ranging by a factor of three between the minimum and maximum estimates. Additionally, we show that using wind‐affected plume models results in affected areas that are five times larger than using no‐wind‐affected models. This demonstrates the sensitivity of MER to the type of plume model chosen (no‐wind‐ vs. wind‐affected)

Inverting sediment bedforms for evaluating the hazard of dilute pyroclastic density currents in the field

Pyroclastic density currents are ground hugging gas-particle flows associated to explosive volcanic eruptions and moving down a volcano's slope, causing devastation and deaths. Because of the hostile nature they cannot be analyzed directly and most of their fluid dynamic behavior is reconstructed by the deposits left in the geological record, which frequently show peculiar structures such as ripples and dune bedforms. Here, a set of equations is simplified to link flow behavior to particle motion and deposition. This allows to construct a phase diagram by which impact parameters of dilute pyroclastic density currents, representing important factors of hazard, can be calculated by inverting bedforms wavelength and grain size, without the need of more complex models that require extensive work in the laboratory

Drag forces at the ice-sheet bed and resistance of hard-rock obstacles:The physics of glacial ripping

Glacial ripping involves glaciotectonic disintegration of rock hills and extensive removal of rock at the ice-sheet bed, triggered by hydraulic jacking caused by fluctuating water pressures. Evidence from eastern Sweden shows that glacial ripping caused significant subglacial erosion during the final deglaciation of the Fennoscandian ice sheet, distinct from abrasion and plucking (quarrying). Here we analyse the ice drag forces exerted onto rock obstacles at the base of an ice sheet, and the resisting forces of such rock obstacles: glaciotectonic disintegration requires that ice drag forces exceed the resisting forces of the rock obstacle. We consider rock obstacles of different sizes, shapes and fracture patterns, informed by natural examples from eastern Sweden. Our analysis shows that limited overpressure events, unfavourable fracture patterns, low-Transmissivity fractures, slow ice and streamlined rock hamper rock hill disintegration. Conversely, under fast ice flow and fluctuating water pressures, disintegration is possible if the rock hill contains subhorizontal, transmissive fractures. Rock steps on previously smooth, abraded surfaces, caused by hydraulic jacking, also enhance drag forces and can cause disintegration of a rock hill. Glacial ripping is a physically plausible erosion mechanism, under realistic glaciological conditions prevalent near ice margins.</p

A study on the influence of internal structures on the shape of pyroclastic particles by X-ray microtomography investigations

1.X-Ray computed microtomography is a non-destructive 3D imaging technique that can be used for the investigation of both the morphology and internal structures of a solid object. Thanks to its versatility, it is currently of common use in many research fields and applications, from medical science to geosciences. The latter include volcanology, where this analytical technique is becoming increasingly popular, in particular for quantifying the shape as well as the internal structure of particles constituting tephra deposits. Particle morphology plays a major role in controlling the mobility of pyroclastic material in the atmosphere and particle-laden flows, while the internal structure (e.g. vesicles and crystal content) is of importance for constraining the processes occurring in magmatic chamber or volcanic conduits. In this paper, we present results of X-Ray microtomography morphological and textural analyses on volcanic particles carried out to study how particle shape is influenced by their internal structures. Particles were selected from tephra generated during explosive eruptions of different magnitude and composition. Results show how particle morphology is strongly influenced by their internal structure, which is characterized by textural features like vesicularity, vesicle and solid structure distribution, vesicles inter-connection and distance between adjacent vesicles. These have been found to vary with magma composition, vesiculation and crystallization history. Furthermore, our results confirm that X-Ray microtomography is a powerful tool for investigating shape and internal structure of particles. It both allows us to characterize the particle shape by means of tridimensional shape parameters and to relate them to their internal structures

Probabilistic hazard analysis of the gas emission of Mefite d'Ansanto, southern Italy

The emission of gas species dangerous to human health and life is a widespread source of hazard in various natural contexts. These mainly include volcanic areas but also non-volcanic geological contexts. A notable example of the latter occurrence is the Mefite d'Ansanto area in the southern Apennines in Italy. Here, large emissions of carbon dioxide (CO2) occur at rates that make this the largest non-volcanic CO2 gas emissions area in Italy and probably on Earth. Given the topography of the area, in certain meteorological conditions a cold-gas stream forms in the valleys surrounding the emission zone, which has proved to be potentially lethal to humans and animals in the past. In this study, we present a gas hazard modelling study that considers the main species, CO2, and the potential effect of another notable species, hydrogen sulfide (H2S). For these purposes, we used VolcanIc Gas dIspersion modeLling v1.3.7 (VIGIL), a tool that manages the workflow of gas dispersion simulations in both the dense- and dilute-gas regimes and is specifically optimised for probabilistic hazard applications. In its latest version, VIGIL can automatically detect the most appropriate regime to simulate based on the gas emission properties and meteorological conditions at the source. Results are discussed and presented in the form of maps of CO2 and H2S concentration and persistence at various exceedance probabilities, which consider the gas emission rates and their possible ranges of variation defined in previous studies. The effect of seasonal variations is also presented and discussed

Investigating Source Conditions and Controlling Parameters of Explosive Eruptions: Some Experimental-Observational- Modelling Case Studies

Explosive volcanic eruptions are complex systems that can generate a variety of hazardous phenomena, for example, the injection of volcanic ash into the atmosphere or the generation of pyroclastic density currents. Explosive eruptions occur when a turbulent multiphase mixture, initially predominantly composedf of fragmented magma and gases, is injected from the volcanic vent into the atmosphere. For plume modelling purposes, a specific volcanic eruption scenario based on eruption type, style or magnitude is strictly linked to magmatic and vent conditions, despite the subsequent evolution of the plume being influenced by the interaction of the erupted material with the atmosphere. In this chapter, different methodologies for investigating eruptive source conditions and the subsequent evolution of the eruptive plumes are presented. The methodologies range from observational techniques to large-scale experiments and numerical models. Results confirm the relevance of measuring and observing source conditions, as such studies can improve predictions of the hazards of eruptive columns. The results also demonstrate the need for fundamental future research specifically tailored to answer some of the still open questions: the effect of unsteady flow conditions at the source on the eruptive column dynamics and the interaction between a convective plume and wind