The source conditions of volcanic plumes and collapsing fountains are investigated by means of large-scale experiments. In the experiments, gas-particle jets issuing from a cylindrical conduit are forced into the atmosphere at different mass flow rates. Dense jets (high particle volumetric concentration, e.g., C 0 > 0.01) generate collapsing fountains, whose height scales with the squared exit velocity. This is consistent with Bernoulli's equation, which is a good approximation if air entrainment is negligible. In this case, kinetic energy is transformed into potential energy without any significant loss by friction with the atmosphere. The dense collapsing fountain, on hitting the ground, generates an intense shear flow similar to a pyroclastic density current. Dilute hot jets (low particle volumetric concentration, e.g., C 0 < 0.01) dissipate their initial kinetic energy at much smaller heights than those predicted by Bernoulli's equation. This is an indication that part of the total mechanical energy is lost by friction with the atmosphere. Significant air entrainment results in this case, leading to the formation of a buoyant column (plume) from which particles settle similarly to pyroclastic fallout. The direct measurement of entrainment coefficient in the experiments suggests that dense collapsing fountains form only when air entrainment is not significant. This is a consequence of the large density difference between the jet and the atmosphere. Cold dilute experiments result in an entrainment coefficient of about 0.06, which is typical of pure jets of fluid dynamics. Hot dilute experiments result in an entrainment coefficient of about 0.11, which is typical of thermally buoyant plumes. The entrainment coefficients obtained by experiments were used as input data in numerical simulations of fountains and plumes. A numerical model was used to solve the classic top-hat system of governing equations, which averages the field variables (e.g., column velocity and density) across the column. The maximum heights calculated with the model agree well with those observed experimentally, showing that our entrainment coefficients are compatible with a top-hat model. Dimensional analysis of the experimental data shows that a value of 3 for the source densimetric Froude number characterizes the transition between dense collapsing fountains and dilute plumes. This value delimits the source conditions (exit velocity, conduit radius, and particle volumetric concentration) for pyroclastic flow (<3) and fallout (>3). © 2014 Springer-Verlag Berlin Heidelberg

[1] It is currently impractical to measure what happens in a volcano during an explosive eruption, and up to now much of our knowledge depends on theoretical models. Here we show, by means of large-scale experiments, that the regime of explosive events can be constrained on the basis of the characteristics of magma at the point of fragmentation and conduit geometry. Our model, whose results are consistent with the literature, is a simple tool for defining the conditions at conduit exit that control the most hazardous volcanic regimes. Besides the well-known convective plume regime, which generates pyroclastic fallout, and the vertically collapsing column regime, which leads to pyroclastic flows, we introduce an additional regime of radially expanding columns, which form when the eruptive gas-particle mixture exits from the vent at overpressure with respect to atmosphere. As a consequence of the radial expansion, a dilute collapse occurs, which favors the formation of density currents resembling natural base surges. We conclude that a quantitative knowledge of magma fragmentation, i.e., particle size, fragmentation energy, and fragmentation speed, is critical for determining the eruption regime.

No abstract

The aspect ratio of ignimbrites is a commonly used parameter that has been related to the energy of the parent pyroclastic density currents (PDCs). However this parameter, calculated as the ratio between the average thickness and the average lateral extent of ignimbrites, does not capture fundamental differences in pyroclastic flow mobility nor relates to lithofacies variations of the final deposits. We herein introduce the “topological aspect ratio” (ARt) as the ratio of the local deposit thickness (Ht) to the distance between the local site and the maximum runout distance (Lt), where Ht is a proxy for the PDC tendency to deposit, and Lt a proxy for the PDC mobility or its tendency to further transport the pyroclastic material. The positive versus negative spatial gradient d(ARt)/dx along flow paths discriminate zones where PDCs are forced (i.e. where they transport the total energy under the action of mass discharge rate) from zones where they are inertial (i.e. where they transport the total energy under the action of viscous or turbulent fluidization). Though simple to apply, the topological aspect ratio and its spatial gradient are powerful descriptors of the interplay between sedimentation and mobility of PDCs, and of the resulting lithofacies variations.

No abstract

Large‐scale experiments generating ground‐hugging multiphase flows were carried out with the aim of modelling the rate of sedimentation, of pyroclastic density currents. The current was initiated by the impact on the ground of a dense gas‐particle fountain issuing from a vertical conduit. On impact, a thick massive deposit was formed. The grain size of the massive deposit was almost identical to that of the mixture feeding the fountain, suggesting that similar layers formed at the impact of a natural volcanic fountain should be representative of the parent grain‐size distribution of the eruption. The flow evolved laterally into a turbulent suspension current that sedimented a thin, tractive layer. A good correlation was found between the ratio of transported/sedimented load and the normalized Rouse number of the turbulent current. A model of the sedimentation rate was developed, which shows a relationship between grain size and flow runout. A current fed with coarser particles has a higher sedimentation rate, a larger grain‐size selectivity and runs shorter than a current fed with finer particles. Application of the model to pyroclastic deposits of Vesuvius and Campi Flegrei of Southern Italy resulted in sedimentation rates falling inside the range of experiments and allowed definition of the duration of pyroclastic density currents which add important information on the hazard of such dangerous flows. The model could possibly be extended, in the future, to other geological density currents as, for example, turbidity currents.

Explosive activity and lava dome collapse at stratovolcanoes can lead to pyroclastic density currents (PDCs; mixtures of volcanic gas, air, and volcanic particles) that produce complex deposits and pose a hazard to surrounding populations. Two-dimensional computer simulations of dilute PDCs (characterized by a turbulent suspended load and deposition through a bed load) show that PDC transport, deposition, and hazard potential are sensitive to the shape of the volcano slope (profile) down which they flow. We focus on three generic volcano profiles: straight, concave-upward, and convex-upward. Dilute PDCs that flow down a constant slope gradually decelerate over the simulated run-out distance (5 km in the horizontal direction) due to a combination of sedimentation, which reduces the density of the PDC, and mixing with the atmosphere. However, dilute PDCs down a concave-upward slope accelerate high on the volcano flanks and have less sedimentation until they begin to decelerate over the shallow lower slopes. A convex-upward slope causes dilute PDCs to lose relatively more of their pyroclast load on the upper slopes of a volcano, and although they accelerate as they reach the lower, steeper slopes, the acceleration is reduced because of the upstream loss of pyroclasts (lower density contrast with the atmosphere). Dynamic pressure, a measure of the damage that can be caused by PDCs, reflects these complex relations

No abstract