Experiments on conduit flow and eruption behavior of basaltic volcanic eruptions

Multiphase flow in basaltic volcanic conduits is investigated using analog experiments and theoretical approaches. Depending on gas supply, large gas bubbles (gas slugs) may rise through basaltic magma in regimes of distinct fluid‐dynamical behavior: ascent of single slugs, supplied slugs fed from the gas source during ascent, and periodic slug flow. An annular flow regime commences at the highest gas supply rates. A first set of experiments demonstrates that the growth of gas slugs due to hydrostatic decompression does not affect their ascent velocity and that excess pressure in the slugs remain negligible. The applicability of theoretical formulae describing slug ascent velocity as a function of liquid and conduit properties is evaluated in a second set of experiments. A third set of experiments with continuous gas supply into a cylindrical conduit are scaled to basaltic conditions over Morton, Eotvös, Reynolds, and Froude numbers. Gas flow rate and liquid viscosity are varied over the whole range of flow regimes to observe flow dynamics and to measure gas and liquid eruption rates. Foam generation by slug bursting at the surface and partial slug disruption by wake turbulence can modify the bubble content and size distribution of the magma. At the transition from slug to annular flow, when the liquid bridges between the gas slugs disappear, pressure at the conduit entrance drops by ∼60% from the hydrostatic value to the dynamic‐flow resistance of the annular flow, which may trigger further degassing in a stored magma to maintain the annular flow regime until the gas supply is exhausted and the eruption ends abruptly. Magma discharge may also terminate when magma ascent is hindered by wall friction in long volcanic conduits and the annular gas flow erodes all magma from the conduit. Supplied slugs are found to reach much higher rise velocities than unsupplied slugs and to collapse to turbulent annular flow upon bursting at the surface. A fourth set of experiments uses a conduit partially blocked by built‐in obstacles providing traps for gas pockets. Once gas pockets are filled, rising gas slugs deform but remain intact as they move around obstacles without coalescence or significant velocity changes. Bursting of bubbles coalescing with trapped gas pockets causes pressure signals at least 3 orders of magnitude more powerful than gas pocket oscillation induced by passing liquid. Our experiments suggest a refined classification of Strombolian and Hawaiian eruptions according to time‐dependant behavior into sporadically pulsating lava fountains (driven by stochastic rise of single slugs), periodically pulsating lava fountains (resulting from slug flow), and quasi‐steady lava fountains (oscillating at the frequency of annular‐flow turbulence)

A criterion for the fragmentation of bubbly magma based on brittle failure theory

The fragmentation of bubbly magma is a defining point in a volcanic eruption-before fragmentation the magma flows relatively slowly, during fragmentation the bubbles break up to release compressed gas and, afterwards, the eruption becomes a violent gas flow carrying suspended magma particles. Seemingly benign lava flows or domes can suddenly fragment into deadly pyroclastic flows(1-3). Several criteria have been proposed to define the point of magma fragmentation or foam stability(4-7). The criterion of Papale(7) is based on melt relaxation theory and equates magma strain rate with the rate of increase of flow velocity with distance. It ignores, however, the role of bubble pressure in causing fragmentation. Two empirical approaches(4,5) consider the role of high bubble pressure in causing fragmentation but do not address the underlying physics of magma fragmentation. Here I develop a fragmentation criterion for bubbly magma based on brittle failure theory and apply it to the fragmentation of lava domes and flows. On the basis of this theory, a bubbly magma will fragment when the tensile stress at the inner walls of bubbles exceeds the tensile strength of the magma. The fragmentation conditions depend strongly on initial water content, with calculated vesicularity and final water levels coinciding reasonably well with those in observed pumices. This suggests that the proposed criterion captures the essence of the fragmentation process in bubbly magma.Peer Reviewedhttp://deepblue.lib.umich.edu/bitstream/2027.42/62724/1/402648a0.pd

The size range of bubbles that produce ash during explosive volcanic eruptions

Volcanic eruptions can produce ash particles with a range of sizes and morphologies. Here we morphologically distinguish two textural types: Simple (generally smaller) ash particles, where the observable surface displays a single measureable bubble because there is at most one vesicle imprint preserved on each facet of the particle; and complex ash particles, which display multiple vesicle imprints on their surfaces for measurement and may contain complete, unfragmented vesicles in their interiors. Digital elevation models from stereo-scanning electron microscopic images of complex ash particles from the 14 October 1974 sub-Plinian eruption of Volcán Fuego, Guatemala and the 18 May 1980 Plinian eruption of Mount St. Helens, Washington, U.S.A. reveal size distributions of bubbles that burst during magma fragmentation. Results were compared between these two well-characterized eruptions of different explosivities and magma compositions and indicate that bubble size distributions (BSDs) are bimodal, suggesting a minimum of two nucleation events during both eruptions. The larger size mode has a much lower bubble number density (BND) than the smaller size mode, yet these few larger bubbles represent the bulk of the total bubble volume. We infer that the larger bubbles reflect an earlier nucleation event (at depth within the conduit) with subsequent diffusive and decompressive bubble growth and possible coalescence during magma ascent, while the smaller bubbles reflect a relatively later nucleation event occurring closer in time to the point of fragmentation. Bubbles in the Mount St. Helens complex ash particles are generally smaller, but have a total number density roughly one order of magnitude higher, compared to the Fuego samples. Results demonstrate that because ash from explosive eruptions preserves the size of bubbles that nucleated in the magma, grew, and then burst during fragmentation, the analysis of the ash-sized component of tephra can provide insights into the spatial distribution of bubbles in the magma prior to fragmentation, enabling better parameterization of numerical eruption models and improved understanding of ash transport phenomena that result in pyroclastic volcanic hazards. Additionally, the fact that the ash-sized component of tephra preserves BSDs and BNDs consistent with those preserved in larger pyroclasts indicates that these values can be obtained in cases where only distal ash samples from particular eruptions are obtainable

Submarine deposits from pumiceous pyroclastic density currents traveling over water: an outstanding example from offshore Montserrat (IODP 340)

Pyroclastic density currents have been observed to both enter the sea, and to travel over water for tens of kilometers. Here, we identified a 1.2-m-thick, stratified pumice lapilli-ash cored at Site U1396 offshore Montserrat (Integrated Ocean Drilling Program [IODP] Expedition 340) as being the first deposit to provide evidence that it was formed by submarine deposition from pumice-rich pyroclastic density currents that traveled above the water surface. The age of the submarine deposit is ca. 4 Ma, and its magma source is similar to those for much younger Soufrière Hills deposits, indicating that the island experienced large-magnitude, subaerial caldera-forming explosive eruptions much earlier than recorded in land deposits. The deposit’s combined sedimentological characteristics are incompatible with deposition from a submarine eruption, pyroclastic fall over water, or a submarine seafloor-hugging turbidity current derived from a subaerial pyroclastic density current that entered water at the shoreline. The stratified pumice lapilli-ash unit can be subdivided into at least three depositional units, with the lowermost one being clast supported. The unit contains grains in five separate size modes and has a >12 phi range. Particles are chiefly subrounded pumice clasts, lithic clasts, crystal fragments, and glass shards. Pumice clasts are very poorly segregated from other particle types, and lithic clasts occur throughout the deposit; fine particles are weakly density graded. We interpret the unit to record multiple closely spaced (<2 d) hot pyroclastic density currents that flowed over the ocean, releasing pyroclasts onto the water surface, and settling of the various pyroclasts into the water column. Our settling and hot and cold flotation experiments show that waterlogging of pumice clasts at the water surface would have been immediate. The overall poor hydraulic sorting of the deposit resulted from mixing of particles from multiple pulses of vertical settling in the water column, attesting to complex sedimentation. Slow-settling particles were deposited on the seafloor together with faster-descending particles that were delivered at the water surface by subsequent pyroclastic flows. The final sediment pulses were eventually deflected upon their arrival on the seafloor and were deposited in laterally continuous facies. This study emphasizes the interaction between products of explosive volcanism and the ocean and discusses sedimentological complexities and hydrodynamics associated with particle delivery to water

Dynamic observations of vesiculation reveal the role of silicate crystals in bubble nucleation and growth in andesitic magmas

Bubble nucleation and growth control the explosivity of volcanic eruptions, and the kinetics of these processes are generally determined from examinations of natural samples and quenched experimental run products. These samples, however, only provide a view of the final state, from which the initial conditions of a time-evolving magmatic system are then inferred. The interpretations that follow are inexact due to the inability of determining the exact conditions of nucleation and the potential detachment of bubbles from their nucleation sites, an uncertainty that can obscure their nucleation location \u2013 either homogeneously within the melt or heterogeneously at the interface between crystals and melts. We present results of a series of dynamic, real-time 4D X-ray tomographic microscopy experiments where we observed the development of bubbles in crystal bearing silicate magmas. Experimentally synthesized andesitic glasses with 0.25\u20130.5 wt% H2O and seed silicate crystals were heated at 1 atm to induce bubble nucleation and track bubble growth and movement. In contrast to previous studies on natural and experimentally produced samples, we found that bubbles readily nucleated on plagioclase and clinopyroxene crystals, that their contact angle changes during growth and that they can grow to sizes many times that of the silicate on whose surface they originated. The rapid heterogeneous nucleation of bubbles at low degrees of supersaturation in the presence of silicate crystals demonstrates that silicates can affect when vesiculation ensues, influencing subsequent permeability development and effusive vs. explosive transition in volcanic eruptions

On the application of scattering matrix measurements to detection and identification of major types of airborne aerosol particles: Volcanic ash, desert dust and pollen

This work has been funded by the excellence research program of the Andalusian Regional Government, grant number P18RT-1854, the National Plan of Scientific and Technical Research and Innovation of the Spanish Ministry of Science and Innovation, grant number RTI2018-095330-B-100 (LEONIDAS), and the "Center of Excellence Severo Ochoa" award to the Instituto de Astrofisica de Andalucia (SEV-2017-0709) by the Spanish State Agency for Research. J.C.G.M acknowledges financial support from the Ramon y Cajal Program of the Spanish Ministry of Science and Innovation (RYC-2016-19570). JoseLuis de la Rosa, JoseAntonio Ruiz and Shi Zongbo are acknowledged for collecting the Sahara-OSN and GobiBeijing desert dust samples.Atmospheric aerosols play key roles in climate and have important impacts on human activities and health. Hence, much effort is directed towards developing methods of improved detection and discrimina- tion of different types of aerosols. Among these, light scattering-based detection of aerosol offers several advantages including applications in both in situ and remote sensing devices. In this work, new scat- tering matrix measurements for two samples of airborne desert dust collected in Spain and China are reported. The average extrapolated scattering matrices of airborne desert dust and of volcanic ash at two wavelengths have been calculated and compared with the aim of finding criteria to distinguish these two types of aerosol. Additionally, the scattering matrix of cypress pollen has been measured and extrapo- lated to explore differences with mineral dust that can be exploited in atmospheric detection. Field mea- surements of the backscattering linear depolarization ratio δL (180 °) are used to obtain information about non-sphericity and discrimination between fine and coarse aerosol. However, the average δL (180 °) for the three types of aerosols considered in this work in the visible spectral range is δL (180 °) = 0.40 ±0.05. This shows that δL (180 °) is not informative about the composition or morphology of irregular particles. By contrast, measurements of scattering matrix elements or depolarization ratios at different scattering angles may provide information about the structural differences of particles, and in particular may en- able to differentiate airborne volcanic ash from desert dust, which are otherwise similar in terms of size and optical constants. Cypress pollen shows a characteristic degree of linear polarization curve that is very different from that of polydisperse irregular mineral dust. Light scattering field instruments and re- mote sensing methods could extract more information about the characteristics of aerosol particles if modifications were introduced to measure the phase curves of several scattering matrix elements or de- polarization ratios.excellence research program of the Andalusian Regional Government P18RT-1854National Plan of Scientific and Technical Research and Innovation of the Spanish Ministry of Science and Innovation RTI2018-095330-B-100Spanish State Agency for Research SEV-2017-0709Spanish Government RYC-2016-1957

Radar Doppler velocity of volcanic eruptions: Theoretical considerations and quantitative documentation of changes in eruptive behaviour at Stromboli volcano, Italy

The use of radar Doppler velocimetry for the observation of volcanic activity is new. We used this method to continuously observe the activity of one vent of Stromboli volcano, Italy, from the end of 2000 April until early May. During this period we recorded 702 eruptions, 132 of which occurred before a strong rain storm passed over the island on April 29. In order to interpret the recorded Doppler data we developed a program that simulates different strombolian eruption scenarios, for which we then calculate the theoretical Doppler spectra. Comparing our theoretical data with the observed data we are able to show that most of the eruptions are nearly vertical, although we did observe only one component of the eruption vector with our Doppler radar. One of the most interesting features of the data set is a significant change in eruptive behaviour that correlates with the occurrence of the rain storm: we find that on average the eruption duration increased by a factor of 2, eruptive velocities were much higher and indirect evidence indicates that the average particle diameter of the erupted material decreased. This change may have several causes, but the coincidence with the rain storm may be evidence of magma–water interaction and feedback on the volcanic activity. If the fluid source (rain) changing the eruptive style is at the surface and in near-surface layers then the main control on final eruption dynamics at Stromboli volcano must also be in rather shallow regions

Silicate melt properties and volcanic eruptions

Peer Reviewedhttp://deepblue.lib.umich.edu/bitstream/2027.42/95246/1/rog1657.pd

Confocal microscopy 3D imaging of diesel particulate matter

To date, diesel particulate matter (DPM) has been described as aggregates of spherule particles with a smooth appearing surface. We have used a new colour confocal microscope imaging method to study the 3D shape of diesel particulate matter (DPM); we observed that the particles can have sharp jagged appearing edges and consistent with these findings, 2D light microscopy demonstrated that DPM adheres to human lung epithelial cells. Importantly, the slide preparation and confocal microscopy method applied avoids possible alteration to the particles' surfaces and enables colour 3D visualisation of the particles. From twenty-one PM particles, the mean (standard deviation) major axis length was 5.6 (2.25) μm with corresponding values for the minor axis length of 3.8 (1.25) μm. These new findings may help explain why air pollution particulate matter (PM) has the ability to infiltrate human airway cells, potentially leading to respiratory tract, cardiovascular and neurological disease