67 research outputs found

    Influence of fluctuating supply on the emplacement dynamics of channelized lava flows

    Get PDF
    The evolution of lava flows emplaced on Mount Etna (Italy) in September 2004 is examined in detail through the analysis ofmorphometricmeasurements of flow units. The growth of the main channelized flow is consistent with a layering of lava blankets, which maintains the initial geometry of the channel (although levees are widened and raised), and is here explicitly related to the repeated overflow of lava pulses. A simple analytical model is introduced describing the evolution of the flow level in a channelized flow unit fed by a fluctuating supply. The model, named FLOWPULSE, shows that a fluctuation in the velocity of lava extrusion at the vent triggers the formation of pulses, which become increasingly high the farther they are from the vent, and are invariably destined to overflow within a given distance. The FLOWPULSE simulations are in accordance with the observed morphology, characterized by a very flat initial profile followed by a massive increase in flow unit cross-section area between 600 and 700 m downflow. The modeled emplacement dynamics provides also an explanation for the observed substantial “loss” of the original flowing mass with increasing distance from the vent

    Uncertainty quantification and sensitivity analysis of volcanic columns models: results from the integral model PLUME-MoM

    Get PDF
    The behavior of plumes associated with explosive volcanic eruptions is complex and dependent on eruptive source parameters (e.g. exit velocity, gas fraction, temperature and grain-size distribution). It is also well known that the atmospheric environment interacts with volcanic plumes produced by explosive eruptions in a number of ways. The wind field can bend the plume but also affect atmospheric air entrainment into the column, enhancing its buoyancy and in some cases, preventing column collapse. In recent years, several numerical simulation tools and observational systems have investigated the action of eruption parameters and wind field on volcanic column height and column trajectory, revealing an important influence of these variables on plume behavior. In this study, we assess these dependencies using the integral model PLUME-MoM, whereby the continuous polydispersity of pyroclastic particles is described using a quadrature-based moment method, an innovative approach in volcanology well-suited for the description of the multiphase nature of magmatic mixtures. Application of formalized uncertainty quantification and sensitivity analysis techniques enables statistical exploration of the model, providing information on the extent to which uncertainty in the input or model parameters propagates to model output uncertainty. In particular, in the framework of the IAVCEI Commission on tephra hazard modeling inter-comparison study, PLUME-MoM is used to investigate the parameters exerting a major control on plume height, applying it to a weak plume scenario based on 26 January 2011 Shinmoe-dake eruptive conditions and a strong plume scenario based on the climatic phase of the 15 June 1991 Pinatubo eruption

    Insights into the formation and dynamics of coignimbrite plumes from one-dimensional models

    Get PDF
    Coignimbrite plumes provide a common and effective mechanism by which large volumes of fine-grained ash are injected into the atmosphere. Nevertheless, controls on formation of these plumes as a function of eruptive conditions are still poorly constrained. Herein, two 1-D axysymmetric steady state models were coupled, the first describing the parent pyroclastic density current and the second describing plume rise. Global sensitivity analysis is applied to investigate controls on coignimbrite plume formation and describe coignimbrite source and the maximum plume height attained. For a range of initial mass flow rates between 108 and 1010 kg/s, modeled liftoff distance (the distance at which neutral buoyancy is attained), assuming radial supercritical flow, is controlled by the initial flow radius, gas mass fraction, flow thickness, and temperature. The predicted decrease in median grain size between flow initiation and plume liftoff is negligible. Calculated initial plume vertical velocities, assuming uniform liftoff velocity over the pyroclastic density current invasion area, are much greater (several tens of m/s) than those previously used in modeling coignimbrite plumes (1 m/s). Such velocities are inconsistent with the fine grain size of particles lofted into coignimbrite plumes, highlighting an unavailability of large clasts, possibly due to particle segregation within the flow, prior to plume formation. Source radius and initial vertical velocity have the largest effect on maximum plume height, closely followed by initial temperature. Modeled plume heights are between 25 and 47 km, comparable with Plinian eruption columns, highlighting the potential of such events for distributing fine-grained ash over significant areas

    Evolution of Conduit Geometry and Eruptive Parameters During Effusive Events

    Get PDF
    The dynamics of effusive events is controlled by the interplay between conduit geometry and source conditions. Dyke-like geometries have been traditionally assumed for describing conduits during effusive eruptions, but their depth-dependent and temporal modifications are largely unknown. We present a novel model which describes the evolution of conduit geometry during effusive eruptions by using a quasi steady state approach based on a 1-D conduit model and appropriate criteria for describing fluid shear stress and elastic deformation. This approach provides time-dependent trends for effusion rate, conduit geometry, exit velocity, and gas flow. Fluid shear stress leads to upward widening conduits, whereas elastic deformation becomes relevant only during final phases of effusive eruptions. Simulations can reproduce different trends of effusion rate, showing the effect of magma source conditions and country rock properties on the eruptive dynamics. This model can be potentially applied for data inversion in order to study specific case studies

    Effusion Rate Evolution During Small-Volume Basaltic Eruptions: Insights From Numerical Modeling

    Get PDF
    The temporal evolution of effusion rate is the main controlling factor of lava spreading and emplacement conditions. Therefore, it represents the most relevant parameter for characterizing the dynamics of effusive eruptions and thus for assessing the volcanic hazard associated with this type of volcanism. Since the effusion rate curves can provide important insights into the properties of the magma feeding system, several efforts have been performed for their classification and interpretation. Here, a recently published numerical model is employed for studying the effects of magma source and feeding dike properties on the main characteristics (e.g., duration, erupted mass, and effusion rate trend) of small-volume effusive eruptions, in the absence of syn-eruptive magma injection from deeper storages. We show that the total erupted mass is mainly controlled by magma reservoir conditions (i.e., dimensions and overpressure) prior to the eruption, whereas conduit processes along with reservoir properties can significantly affect mean effusion rate, and thus, they dramatically influence eruption duration. Simulations reproduce a wide variety of effusion rate trends, whose occurrence is controlled by the complex competition between conduit enlargement and overpressure decrease due to magma withdrawal. These effusion rate curves were classified in four groups, which were associated with the different types described in the literature. Results agree with the traditional explanation of effusion rate curves and provide new insights for interpreting them, highlighting the importance of magma reservoir size, initial overpressure, and initial width of the feeding dike in controlling the nature of the resulting effusion rate curve

    Lahar events in the last 2000 years from Vesuvius eruptions – Part 2: Formulation and validation of a computational model based on a shallow layer approach

    Get PDF
    In this paper we present a new model for the simulation of lahars based on the depth-averaged code IMEX-SfloW2D with new governing and constitutive equations introduced to better describe the dynamics of lahars. A thorough sensitivity analysis is carried out to identify the critical processes (such as erosion and deposition) and parameters (both numerical and physical) controlling lahar runout using both synthetic and real case topographies. In particular, an application of the model to a syn-eruptive lahar from a reference size eruption from Somma–Vesuvius, affecting the Campanian Plain (southern Italy), described in Di Vito et al. (2024), is used in this work for the sensitivity analysis. Effects of erosion and deposition are investigated by comparing simulations with and without these processes. By comparing flow thickness and area covered by the flow and their evolution with time, we show that the modelling of both the processes is important to properly simulate the effects of the bulking and debulking as well as the associated changes in rheology. From a computational point of view, the comparisons of simulations obtained for different numerical grids (from 25 to 100 m), scheme order, and grain size discretization were useful to find a good compromise between resolution and computational speed. The companion paper by Sandri et al. (2024) shows an application of the presented model for probabilistic volcanic hazard assessment for lahars from Vesuvius deposits in the Neapolitan area.</p

    Eruption type probability and eruption source parameters at Cotopaxi and Guagua Pichincha volcanoes (Ecuador) with uncertainty quantification

    Get PDF
    Future occurrence of explosive eruptive activity at Cotopaxi and Guagua Pichincha volcanoes, Ecuador, is assessed probabilistically, utilizing expert elicitation. Eight eruption types were considered for each volcano. Type event probabilities were evaluated for the next eruption at each volcano and for at least one of each type within the next 100&nbsp;years. For each type, we elicited relevant eruption source parameters (duration, average plume height, and total tephra mass). We investigated the robustness of these elicited evaluations by deriving probability uncertainties using three expert scoring methods. For Cotopaxi, we considered both rhyolitic and andesitic magmas. Elicitation findings indicate that the most probable next eruption type is an andesitic hydrovolcanic/ash-emission (~ 26–44% median probability), which has also the highest median probability of recurring over the next 100&nbsp;years. However, for the next eruption at Cotopaxi, the average joint probabilities for sub-Plinian or Plinian type eruption is of order 30–40%—a significant chance of a violent explosive event. It is inferred that any Cotopaxi rhyolitic eruption could involve a longer duration and greater erupted mass than an andesitic event, likely producing a prolonged emergency. For Guagua Pichincha, future eruption types are expected to be andesitic/dacitic, and a vulcanian event is judged most probable for the next eruption (median probability ~40–55%); this type is expected to be most frequent over the next 100&nbsp;years, too. However, there is a substantial probability (possibly &gt;40% in average) that the next eruption could be sub-Plinian or Plinian, with all that implies for hazard levels
    • …
    corecore