An analytical model for gas overpressure in slug-driven explosions:insights into Strombolian volcanic eruptions

Strombolian eruptions, common at basaltic volcanoes, are mildly explosive events that are driven by a large bubble of magmatic gas (a slug) rising up the conduit and bursting at the surface. Gas overpressure within the bursting slug governs explosion dynamics and vigor and is the main factor controlling associated acoustic and seismic signals. We present a theoretical investigation of slug overpressure based on magma-static and geometric considerations and develop a set of equations that can be used to calculate the overpressure in a slug when it bursts, slug length at burst, and the depth at which the burst process begins. We find that burst overpressure is controlled by two dimensionless parameters: V', which represents the amount of gas in the slug, and A', which represents the thickness of the film of magma that falls around the rising slug. Burst overpressure increases nonlinearly as V' and A' increase. We consider two eruptive scenarios: (1) the "standard model," in which magma remains confined to the vent during slug expansion, and (2) the " overflow model," in which slug expansion is associated with lava effusion, as occasionally observed in the field. We find that slug overpressure is higher for the overflow model by a factor of 1.2-2.4. Applying our model to typical Strombolian eruptions at Stromboli, we find that the transition from passive degassing to explosive bursting occurs for slugs with volume >24-230 m(3), depending on magma viscosity and conduit diameter, and that at burst, a typical Strombolian slug (with a volume of 100-1000 m(3)) has an internal gas pressure of 1-5 bars and a length of 13-120 m. We compare model predictions with field data from Stromboli for low-energy " puffers," mildly explosive Strombolian eruptions, and the violently explosive 5 April 2003 paroxysm. We find that model predictions are consistent with field observations across this broad spectrum of eruptive styles, suggesting a common slug-driven mechanism; we propose that paroxysms are driven by unusually large slugs (large V')

Sequential Fragmentation / Transport Theory, Pyroclast Size-Density Relationships, and the Emplacement Dynamics of Pyroclastic Density Currents – A Case Study on the Mt. St. Helens (USA) 1980 Eruption

Pyroclastic density currents (PDCs) are the most dangerous hazard associated with explosive volcanic eruptions. Despite recent advancements in the general understanding of PDC dynamics, limited direct observation and/or outcrop scarcity often hinder the interpretation of specific transport and depositional processes at many volcanoes. This study explores the potential of sequential fragmentation / transport theory (SFT; cf. Wohletz et al. 1989), a modeling method capable of predicting particle mass distributions based on the physical principles of fragmentation and transport, to retrieve the transport and depositional dynamics of well-characterized PDCs from the size and density distributions of individual components within the deposits. The extensive vertical and lateral exposures through the May 18th, 1980 PDC deposits at Mt. St. Helens (MSH) provide constraints on PDC regimes and flow boundary conditions at specific locations across the depositional area. Application to MSH deposits suggests that SFT parameter distributions can be effectively used to characterize flow boundary conditions and emplacement processes for a variety of PDC lithofacies and deposit locations. Results demonstrate that (1) the SFT approach reflects particle fragmentation and transport mechanisms regardless of variations in initial component distributions, consistent with results from previous studies; (2) SFT analysis reveals changes in particle characteristics that are not directly observable in grain size and fabric data; (3) SFT parameters are more sensitive to regional transport conditions than local (outcrop-scale) depositional processes. The particle processing trends produced using SFT analysis are consistent with the degree of particle processing inferred from lithofacies architectures: for all lithofacies examined in this study, suspension sedimentation products exhibit much better processing than concentrated current deposits. Integrated field observations and SFT results provide evidence for increasing density segregation within the depositional region of the currents away from source, as well as for comparable density-segregation processes acting on lithic concentrations and pumice lenses within the current. These findings further define and reinforce the capability of SFT analysis to complement more conventional PDC study methods, significantly expanding the information gained regarding flow dynamics. Finally, this case study demonstrates that the SFT methodology has the potential to constrain regional flow conditions at volcanoes where outcrop exposures are limited

Maar-diatreme geometry and deposits: Subsurface blast experiments with variable explosion depth

Basaltic maar-diatreme volcanoes, which have craters cut into preeruption landscapes (maars) underlain by downward-tapering bodies of fragmental material commonly cut by hypabyssal intrusions (diatremes), are produced by multiple subsurface phreatomagmatic explosions. Although many maar-diatremes have been studied, the link between explosion dynamics and the resulting deposit architecture is still poorly understood. Scaled experiments employed multiple buried explosions of known energies and depths within layered aggregates in order to assess the effects of explosion depth, and the morphology and compaction of the host on the distribution of host materials in resulting ejecta, the development of subcrater structures and deposits, and the relationships between them. Experimental craters were 1–2 m wide. Analysis of high-speed video shows that explosion jets had heights and shapes that were strongly influenced by scaled depth (physical depth scaled against explosion energy) and by the presence or absence of a crater. Jet properties in turn controlled the distribution of ejecta deposits outside the craters, and we infer that this is also reflected in the diverse range of deposit types at natural maars. Ejecta were dominated by material that originated above the explosion site, and the shallowest material was dispersed the farthest. Subcrater deposits illustrate progressive vertical mixing of host materials through successive explosions. We conclude that the progressive appearance of deeper-seated material stratigraphically upward in deposits of natural maars probably records the length and time scale for upward mixing through multiple explosions with ejection by shallow blasts, rather than progressive deepening of explosion sites in response to draw down of aquifers

The electrical signature of mafic explosive eruptions at Stromboli volcano, Italy

Volcanic lightning is commonly observed in explosive volcanic eruptions of Volcanic Explosivity Index (VEI) > 2 and can be detected remotely providing real-time volcano monitoring information. However, little is known about the electrical activity accompanying the lower-magnitude spectrum of explosive eruptions, often involving mafic magmas. We narrow this gap in knowledge by presenting the electrical signature of the explosive activity (VEI ≤ 1) of Stromboli volcano (Italy) recorded by an electrostatic thunderstorm detector. The persistent eruptive activity of mild Strombolian explosions is occasionally interrupted by larger-scale major explosions and paroxysmal events. Here, we present electrical observations of three major explosions and unprecedented measurements of the 3 July 2019 paroxysm. The electrical signals of the major explosions show apparent similarities, with movements of charge and tens of electrical discharges, arising the question of whether these observations could be used to supplement the classification scheme of explosions on Stromboli. The electrical signals from the 3 July 2019 paroxysm exceed those from the major explosions in amplitude, discharge rate and complexity, showing characteristic variations during different phases of the eruption. These results show that also impulsive lower-magnitude explosions generate detectable electrical activity, which holds promise for monitoring low VEI activity at mafic volcanoes

Multiparametric approach to unravel the mechanism of Strombolian activity at a multivent system: Mt. Etna case study

On 5th July 2014 an eruptive fissure (hereafter referred to as EF) opened at the base of North-East Crater (NEC) of Mt. Etna. EF produced both Strombolian explosions and lava effusion. Thanks to the multiparametric experiment planned in the framework of MEDSUV project, we had the chance to acquire geophysical and volcanological data, in order to investigate the ongoing volcanic activity at EF. Temporary instruments (2 broadband seismometers, 2 microphones, 3-microphone arrays, a high-speed video camera and a thermal-camera) were deployed near the active vents during 15-16 July 2014 and were integrated with the data recorded by the permanent networks. Several kinds of studies are currently in progress, such as: frequency analysis by Fourier Transform and Short Time Fourier Transform to evaluate the spectral content of both seismic and acoustic signals; partitioning of seismic and acoustic energies, whose time variations could reflect changes in the volcanic dynamics; investigation on the intertimes between explosions to investigate their recurrence behaviour; classification of the waveforms, of infrasound events. Furthermore, joint analysis of video signals and seismic-acoustic wavefields outlined relationships between pyroclasts ejection velocity, total erupted mass, peak explosion pressure, and air-ground motion coupling. This multiparametric approach allowed distinguishing and characterizing individually the behavior of the two vents active along the eruptive fissure via their thermal, visible and infrasonic signatures and shed light in the eruptive dynamics.UnpublishedVienna (Austria)5V. Processi eruttivi e post-eruttiv

Effect of particle volume fraction on the settling velocity of volcanic ash particles: insights from joint experimental and numerical simulations

Most of the current ash transport and dispersion models neglect particle-fluid (two-way) and particle-fluid plus particle-particle (four-way) reciprocal interactions during particle fallout from volcanic plumes. These interactions, a function of particle concentration in the plume, could play an important role, explaining, for example, discrepancies between observed and modelled ash deposits. Aiming at a more accurate prediction of volcanic ash dispersal and sedimentation, the settling of ash particles at particle volume fractions (phi(p)) ranging 10(-7)-10(-3) was performed in laboratory experiments and reproduced by numerical simulations that take into account first the two-way and then the four-way coupling. Results show that the velocity of particles settling together can exceed the velocity of particles settling individually by up to 4 times for phi(p) similar to 10(-3). Comparisons between experimental and simulation results reveal that, during the sedimentation process, the settling velocity is largely enhanced by particle-fluid interactions but partly hindered by particle-particle interactions with increasing phi(p). Combining the experimental and numerical results, we provide an empirical model allowing correction of the settling velocity of particles of any size, density, and shape, as a function of phi(p). These corrections will impact volcanic plume modelling results as well as remote sensing retrieval techniques for plume parameters

Time-series analysis of fissure-fed multi-vent activity: a snapshot from the July 2014 eruption of Etna volcano (Italy)

The April to May 2010 eruption of Eyjafjallajökull (Iceland) volcano was characterized by a large compositional variability of erupted products. To contribute to the understanding of the plumbing system dynamics of this volcano, we present new EMPA and LA-ICP-MS data on groundmass glasses of ash particles and minerals erupted between April 15 and 22. The occurrence of disequilibrium textures in minerals, such as resorption and inverse zoning, indicate that open system processes were involved in determining the observed compositional variability. The variation of major and trace element data of glasses corroborates this hypothesis indicating that mixing between magma batches with different compositions interacted throughout the whole duration of the eruption. In particular, the arrival of new basaltic magma into the plumbing system of the volcano destabilized and remobilized magma batches of trachyandesite and rhyolite compositions that, according to geophysical data, might have intruded as sills over the past 20 years beneath the Eyjafjallajökull edifice. Two mixing processes are envisaged to explain the time variation of the compositions recorded by the erupted tephra. The first occurred between basaltic and trachyandesitic end-members. The second occurred between trachyandesite and rhyolites. Least-squares modeling of major elements supports this hypothesis. Furthermore, investi- gation of compositional histograms of trace elements allows us to estimate the initial proportions of melts that interacted to generate the compositional variability triggered by mixing of trachyandesites and rhyolites.Published515V. Dinamica dei processi eruttivi e post-eruttiviJCR Journa

Unoccupied Aircraft Systems (UASs) Reveal the Morphological Changes at Stromboli Volcano (Italy) before, between, and after the 3 July and 28 August 2019 Paroxysmal Eruptions

In July and August 2019, two paroxysmal eruptions dramatically changed the morphology of the crater terrace that hosts the active vents of Stromboli volcano (Italy). Here, we document these morphological changes, by using 2259 UAS-derived photographs from eight surveys and Structure-from-Motion (SfM) photogrammetric techniques, resulting in 3D point clouds, orthomosaics, and digital surface models (DSMs) with resolution ranging from 8.1 to 12.4 cm/pixel. We focus on the morphological evolution of volcanic features and volume changes in the crater terrace and the upper part of the underlying slope (Sciara del Fuoco). We identify both crater terrace and lava field variations, with vents shifting up to 47 m and the accumulation of tephra deposits. The maximum elevation changes related to the two paroxysmal eruptions (in between May and September 2019) range from +41.4 to −26.4 m at the lava field and N crater area, respectively. Throughout September 2018–June 2020, the total volume change in the surveyed area was +447,335 m3. Despite Stromboli being one of the best-studied volcanoes worldwide, the UAS-based photogrammetry products of this study provide unprecedented high spatiotemporal resolution observations of its entire summit area, in a period when volcanic activity made the classic field inspections and helicopter overflights too risky. Routinely applied UAS operations represent an effective and evolving tool for volcanic hazard assessment and to support decision-makers involved in volcanic surveillance and civil protection operations

Particle size-density relationships in pyroclastic deposits: inferences for emplacement processes

In moving pyroclastic systems, particles are sorted as a function of their sizes, densities and shapes. The analysis of the distribution of these characteristics in particle populations of pyroclastic deposits is a major tool in evaluating properties and regimes of parent transport systems. This paper reports a study of size-density relationships of component subpopulations from selected Quaternary pyroclastic deposits of the Latera Volcanic Complex (Vulsini, Italy), by means of integrated field and laboratory investigations, including grain-size and component analyses, density measurements and SEM crystal morphoscopy. Three kinds of particle size-density relationships characterize the studied deposits, as follows: (1) a significant correlation of decreasing sizes with increasing component densities exists for direct suspension sedimentation products, including Plinian fall and pyroclastic surge deposits; (2) a wide data dispersion with poor correlation normally occurs for the main zones of deposits from (high-concentration) pyroclastic flows; and (3) component subpopulations occur in a narrow size range irrespective of density in granular systems dominated by grain collisions, such as traction zones beneath (low-concentration) pyroclastic surges and sheared basal layers of pyroclastic flows. The regression coefficients of size-density relationships appear as a comparative measure of particle concentration in the parent systems, consistent with that inferred from the degree of sorting and abrasion of crystal subpopulations from sequential fragmentation/transport analysis of grain-size data, and quantitative morphoscopy. The crystal component becomes less sorted and more abraded as the fit to the size-density relationship decreases from Plinian fall to pyroclastic surge to pyroclastic flow, consistent with increasing particle concentration in the transporting current