Vent Geometry Research Articles

Hot gas impingement and radiation on neighboring surfaces from venting and combustion in a package of 18650 cells

A quasi-steady, CFD-based modeling approach is employed to investigate the heat loading within a small package of twenty-five 18650 Li-ion cells. The quasi-steady approach allows for computationally efficient simulations to capture the compressible and turbulent flow field through the safety vent structure and out into the space surrounding a failing cell. Combustion of vent gases leads to high heat loading on neighboring cells and nearby surfaces. Heat transfer mechanisms within the enclosure include convection from hot gases, radiation from the participating medium, and radiation exchange between surfaces. Simulations provide insight into the magnitude of each heat transfer mechanism, and the spatial distribution of heat flux on nearby cells and surfaces within the pack. The complex geometry of the safety vent geometry resulted in an asymmetric jet flow pattern, which induces highly localized impingement heat transfer on specific cells within the enclosure. Radiation from hot surfaces was more significant than radiation from hot gases and soot to neighboring cells. The quasi-steady simulations may be used in the future to develop reduced-order heat transfer models that include the effects of venting and combustion on propagating failure.

Morphological transitions between lobate resurfacing and distal breakout lava flows in flood basalts: insights from analog experiments

Continental flood basalts (CFBs) are dominated by two characteristic lava morphologies. The first type, referred to as ‘compound’ or ‘hummocky pāhoehoe,’ exhibits pillow-like lava flow lobes with cross-sections of ~ 0.5–2 m and thin chilled margins. The second type, referred to as ‘simple’ or ‘sheet lobes’ preserves more massive, inflated flow interiors that are laterally continuous on scales of 100s of meters to kilometers. Previous hypotheses suggest that two factors may contribute to stratigraphic changes in morphology from ‘compound’ to ‘simple’: 1) increased eruption duration or 2) increased extrusion rate. We test the hypothesis that a large increase in extrusion rate would result in flow morphology transitioning from multiple small lobes to inflated sheet lobes due to a shift in flow propagation from intraflow resurfacing-dominated to marginal breakout-dominated. Using polyethylene glycol (PEG) wax extruded into a circular water-filled tank 130 cm in diameter, we produced larger, more complex experiments than previous studies. Our efforts simulated more complex lava fields which change flow morphology with distance from the eruptive vent, characteristic of CFBs. Whereas previous PEG studies linked extrusion rate to near-source surface morphologies, our experiments evaluated how flow propagation mechanisms change with variable extrusion rate and distance from the source. Two flow propagation styles were identified: 1) resurfacing, in which molten material breaks through the surface of a flow and covers the older crust and 2) marginal breakouts, in which molten material extends beyond the crust at the active distal margin of the flow. Flows that propagated via marginal breakouts were found to have lower proportions of resurfaced area and vice versa. We show that significant resurfacing is needed to preserve internal chilled boundaries within a flow and a low-extrusion-rate surface morphology, whereas marginal breakout-dominated flows tend to inflate the pillow-like surface morphology preserving a massive interior at great distances from the vent. Higher and more steady extrusion rates tend to decrease the extent of resurfacing and increase the distance between the source and preserved low-extrusion-rate surface morphologies. We find that an extrusion rate increase equivalent to a jump in the extrusion rate scaling factor, Ψ value, from < 1 to > 5 would be necessary to ensure a switch from resurfacing-dominated lobate morphologies to marginal breakout-dominated propagation style. This amounts to a factor of 125 increase in effusion rate for fissure eruptions and a factor of 625 for point source eruptions, assuming no change in vent geometry. This would be equivalent to an effusion rate of 0.2 m3/s, as documented in 1987–1990 Kīlauea eruptions, increasing to 125 m3/s, which was commonly measured during the 2014 Holuhraun eruption in Iceland and the 2018 eruption at Leilani Estates in Hawai‘i. Thus, we propose that continental flood basalts do not require unusually large effusion rates, but instead were active for a longer and more consistent time period than smaller-volume eruptions.

Assessment of eruption source parameters using infrasound and plume modelling: a case study from the 2021 eruption of Mt. Etna, Italy

Atmospheric injection of volcanic ash during eruptions is a threat to aviation. Reliable forecast of airborne ash dispersal relies on empirical and numerical models. Key inputs into these models are so-called eruption source parameters such as the rate at which pyroclastic material is ejected from the vent and the maximum height of eruptive columns. Here, we use infrasound data recorded during eruptive activity in June 2021 at Mt. Etna, Italy, to demonstrate its potential for assessment of eruption rates in near-real time. We calculate a time series of flow velocity at the vent using data corrected for topographic scattering, and the effect of vent geometry on the acoustic source radiation. We obtain values of flow velocity of 50–125 m/s during a period of sustained, paroxysmal, activity. We use independent estimates from other ground-based remote sensing data to validate our results. Further, we use the infrasound-derived flow velocities as input into a 1D plume model to estimate the maximum height of the eruption column. Our results suggest that infrasound technology holds promise to assess eruption rates and inform modelling of volcanic plumes. We anticipate that implementation of real-time operational workflows based on infrasound data and plume modelling will impact decision-making and risk mitigation at active volcanoes.

Crystals and inclined conduits: analogue experiments for slug-driven volcanism

Basaltic volcanism is the dominant mode of volcanism on Earth and exhibits a range of activity, from passive degassing to the most common explosive style: strombolian volcanism. Strombolian volcanism is driven by gas slugs, making it vital to consider the effects of variable magmatic rheology and internal vent geometry on slug flow dynamics. Emerging experimental technologies play a major role in developing our understanding of the natural complexity of such basaltic systems. This study examines slug ascent within particle-free and particle-containing media experimentally across a range of inclinations. Dimensionless parameters are derived to describe specific flow characteristics at laboratory and volcanic scales, and to demonstrate the viability of current theoretical framework. Slug ascent is shown to be dependent on its morphology, which is a function of inclination, liquid viscosity, and related controlling characteristics i.e. particle fraction. Maxima for ascent velocities and associated dimensionless parameters occur within the range 40—60°.

Morphometric trends and implications for the formation of araneiform clusters

Araneiforms are an exotic type of feature exclusive to the Martian south polar terrain which have no known direct Earth analog. The dilation and recession of the seasonal CO2 ice deposit has long been hypothesized to be responsible for their gradual formation. However, high south polar latitude araneiforms have not been observed to form or extend in the present day and it is unclear whether they are growing below the limits of high resolution imagery, or whether they formed in a paleoclimate. Furthermore, it is unclear what physical factors govern the distinct araneiform morphologies, which often appear in discrete ‘clusters’, showing non-random spatial distribution. While qualitative observations of morphology have been reported, there has been a paucity of quantitative measurements of the many physical controls that govern araneiform activity and morphology. Here we present the results of a survey designed to initiate a quantitative morphometric approach to exploring some of the factors which may influence araneiform morphologies. We survey distinct clusters within 3 sites in the Martian south polar regions representing unique araneiform morphologies ranging from ‘fat’ to ‘thin’ to starburst' in order to test the relationship between vent spacing and diameter to morphology found in experimental results. Using the diameters of their central pits as proxies for their starting vent diameters, we find that the level of branching of these dendritic features appears to be correlated with (i) inferred vent spacing and (ii) inferred vent diameter. We suggest that the individual araneiforms within such distinct clusters may have formed contemporaneously, and that vent geometry reflects at least one control on their distinct morphologies. We surveyed and measured araneiforms in three HiRISE images and we have observed and characterized new morphologies at eight sites in HiRISE and CTX images between lat=-87.521°N and lon=162.39°E and lat=-74.7°N, lon=168.41°E in order to investigate if there are any trends in araneiform morphology with latitude. For the sites we studied, we find no discernable relationship between morphology and latitude, though we caution that a more extended survey is needed to confirm this observation on the broader scale. Additionally, from this survey, we report six new araneiform morphologies; ‘fan’, ‘linear’, ‘spiky’, ‘segmented’, ‘flowing’ and ‘sunflower’.

The initial phase of the 2021 Cumbre Vieja ridge eruption (Canary Islands): Products and dynamics controlling edifice growth and collapse

Tajogaite cone in the Cumbre Vieja ridge (La Palma, Canary Islands) erupted between 19 September and 13 December 2021. The tephra and lava sourced from the newly formed fissure rapidly built a pyroclastic cone. During the early days of eruption and after several small-scale landslides, the west flank of the edifice partially collapsed on 25 September, breaching the cone and emplacing a prominent raft-bearing lava flow. Our research combines direct observations, digital elevation models, thermal and visible imaging, and textural and compositional investigation of the explosive products to describe and characterize the edifice growth and collapse. The cone built over a steep slope (26°) and its failure occurred after an intense phase of lava fountaining (up to 30 m3 s−1) that produced rapid pyroclastic accumulation. We suggest that an increased magma supply, to an ascent rate of 0.30 m s−1, led to the rapid growth of the cone (at 2.4 × 106 m3 day−1). Simultaneously, the SW lava flow reactivated and formed a lava ‘seep’ that undercut the flank of the cone, triggering a lateral collapse via rotational rockslide that moved at minimum speeds of 34–70 m h−1. The lateral collapse formed a ~ 200 m wide scar, involving 5.5 × 106 m3 of material, and covered 1.17 km2 with decametric edifice portions and raft-bearing lava. The collapse produced a modest change in the vent geometry, but did not affect eruptive activity long term. A short pause in the eruption after the collapse may have been favored by rapid emptying of the shallower magma system, reducing ascent rates and increasing crystallization times. These results reveal the complex chain of events related to the growth and destruction of newly formed volcanic cones and highlight hazards when situated close to inhabited areas.

Complex geometry of volcanic vents and asymmetric particle ejection: experimental insights

Explosive volcanic eruptions eject a gas-particle mixture into the atmosphere. The characteristics of this mixture in the near-vent region are a direct consequence of the underlying initial conditions at fragmentation and the geometry of the shallow plumbing system. Yet, it is not possible to observe directly the sub-surface parameters that drive such eruptions. Here, we use scaled shock-tube experiments mimicking volcanic explosions in order to elucidate the effects of a number of initial conditions. As volcanic vents can be expected to possess an irregular geometry, we utilise three vent designs, two “complex” vents and a vent with a “real” volcanic geometry. The defining geometry elements of the “complex” vents are a bilateral symmetry with a slanted top plane. The “real” geometry is based on a photogrammetric 3D model of an active volcanic vent with a steep and a diverging vent side. Particle size and density as well as experimental pressure are varied. Our results reveal a strong influence of the vent geometry, on both the direction and the magnitude of particle spreading and the velocity of particles. The overpressure at the vent herby controls the direction of the asymmetry of the gas-particle jet. These findings have implications for the distribution of volcanic ejecta and resulting areas at risk.

The impact of vent geometry on the growth of lava domes

SUMMARY Thick lava flows that are a feature of many volcanic fields on the Earth and Venus vary from sheet-like to nearly perfect axisymmetric domes. Here, we investigate how these geometrical characteristics depend on the shape of the feeder vent. We study the gravitational spreading of viscous lava erupting from elliptical vents onto a flat surface using 3-D numerical models. The aspect ratio of the vent, defined to be the major to minor axes ratio, varies between 1 and 25. In the limit of an aspect ratio of one, the vent is circular and spreading is axisymmetric. In the limit where the ratio is large, the vent behaves as a fissure. The numerical models rely on an isoviscous lava rheology and a constant volumetric eruption rate. In all cases, the initial phase of the dome’s evolution is in a lava-discharge dominated regime such that spreading is insignificant and the height of the dome increases at a constant rate over the vent area. For vent aspect ratios greater than five, three successive regimes of spreading are identified: 2-D spreading in the direction perpendicular to the major axis of the vent, a transient phase such that the dome shape evolves towards that of a circular dome and a late axisymmetric spreading phase that does not depend on the vent shape. These regimes are delimited by the times required for the flow thickness above the vent to reach a given height and for the flow to spread axisymmetrically up to a length equal to the semi-major axis of the vent. Numerical results for the flow height and runout length tend towards the similarity solutions in the 2-D and axisymmetric regimes. Two main implications for highly viscous (rhyolitic) fissure eruptions can be drawn. First, the fissure length determines the flow regimes. The longer the vent fissure length, the longer the early lava discharge regime and 2-D spreading perpendicular to the length of the fissure. Second, the aspect ratio of fissure-fed lava flows can be used as an indicator of the fissure length and the duration of lava discharge. The ellipticity of some terrestrial fissure-fed flows provides evidence for viscous gravity-driven spreading terminated before the onset of the axisymmetric regime. On the other hand, the circular domes on Venus appear to be the result of fissure-fed eruptions sustained enough for the spreading to reach the axisymmetric regime. We propose relationships providing estimates of the fissure length and the duration of lava discharge based on fossil dome dimensions.

Linking gas and particle ejection dynamics to boundary conditions in scaled shock-tube experiments

Predicting the onset, style and duration of explosive volcanic eruptions remains a great challenge. While the fundamental underlying processes are thought to be known, a clear correlation between eruptive features observable above Earth’s surface and conditions and properties in the immediate subsurface is far from complete. Furthermore, the highly dynamic nature and inaccessibility of explosive events means that progress in the field investigation of such events remains slow. Scaled experimental investigations represent an opportunity to study individual volcanic processes separately and, despite their highly dynamic nature, to quantify them systematically. Here, impulsively generated vertical gas-particle jets were generated using rapid decompression shock-tube experiments. The angular deviation from the vertical, defined as the “spreading angle”, has been quantified for gas and particles on both sides of the jets at different time steps using high-speed video analysis. The experimental variables investigated are 1) vent geometry, 2) tube length, 3) particle load, 4) particle size, and 5) temperature. Immediately prior to the first above-vent observations, gas expansion accommodates the initial gas overpressure. All experimental jets inevitably start with a particle-free gas phase (gas-only), which is typically clearly visible due to expansion-induced cooling and condensation. We record that the gas spreading angle is directly influenced by 1) vent geometry and 2) the duration of the initial gas-only phase. After some delay, whose length depends on the experimental conditions, the jet incorporates particles becoming a gas-particle jet. Below we quantify how our experimental conditions affect the temporal evolution of these two phases (gas-only and gas-particle) of each jet. As expected, the gas spreading angle is always at least as large as the particle spreading angle. The latter is positively correlated with particle load and negatively correlated with particle size. Such empirical experimentally derived relationships between the observable features of the gas-particle jets and known initial conditions can serve as input for the parameterisation of equivalent observations at active volcanoes, alleviating the circumstances where an a priori knowledge of magma textures and ascent rate, temperature and gas overpressure and/or the geometry of the shallow plumbing system is typically chronically lacking. The generation of experimental parameterisations raises the possibility that detailed field investigations on gas-particle jets at frequently erupting volcanoes might be used for elucidating subsurface parameters and their temporal variability, with all the implications that may have for better defining hazard assessment.

Characterising vent and crater shape changes at Stromboli: implications for risk areas

Active volcanoes are typically subject to frequent substantial topographic changes as well as variable eruption intensity, style and/or directionality. Gravitational instabilities and local accumulation of pyroclasts affect conditions at the active vents, through which gas-particle jets are released. In turn, the vent geometry strongly impacts the eruption characteristics. Here, we compare five high-resolution topographic data sets (&lt;4 cm/pixel) of volcanic craters and vents from Stromboli volcano, Italy, that were acquired by unoccupied aerial vehicle (UAV) during five field campaigns between May 2019 and January 2020. This period includes two paroxysmal explosions (3 July and 28 August 2019) and exhibited significant changes on day-to-month timescales. Our results highlight changes to vent geometry and their strong control on the directionality of explosions. Recurrent UAV surveys enable the monitoring of temporal morphologic changes and aid the interpretation of observed changes in eruption style. Ultimately, this may contribute to repeatedly revised risk areas on permanently active volcanoes, especially those that are important tourist destinations.

Onset and evolution of Kīlauea's 2018 flank eruption and summit collapse from continuous gravity

Prior to the 2018 lower East Rift Zone (ERZ) eruption and summit collapse of Kīlauea Volcano, Hawai‘i, continuous gravimeters operated on the vent rims of ongoing eruptions at both the summit and Pu‘u ‘Ō‘ō. These instruments captured the onset of the 2018 lower ERZ eruption and the effects of lava withdrawal from both locales, providing constraints on the timing and style of activity and the physical properties of the lava lakes at both locations. At the summit, combining gravity, lava level, and a three-dimensional model of the vent indicates that the upper ∼200 m of the lava lake had a density of about 1700 kgm−3, slightly greater than estimates from 2011–2015 and possibly indicating a gradual densification over time. At Pu‘u ‘Ō‘ō, gravity and vent geometry were used to model both the density and the rate of crater collapse, which was unknown owing to a lack of visual observations. Results suggest the withdrawal of at least 11×106 m3 of lava over the course of two hours, and a material density of 1800–1900 kgm−3. In addition, gravity data at Pu‘u ‘Ō‘ō captured a transient decrease and increase about an hour prior to crater collapse and that was probably related to a small, short-lived fissure eruption on the west flank of the cone and possibly to dike intrusion beneath Pu‘u ‘Ō‘ō. The fissure was the first event in the subsequent cascade that ultimately led to the extrusion of over 1 km3 of lava from lower ERZ vents, collapse of the summit caldera floor by more than 500 m, and the destruction of over 700 homes and other structures. These results emphasize the importance of continuous gravity in operational monitoring of active volcanoes.

Release characteristics of overpressurised gas from complex vents: implications for volcanic hazards

Many explosive volcanic eruptions produce underexpanded starting gas-particle jets. The dynamics of the accompanying pyroclast ejection can be affected by several parameters, including magma texture, gas overpressure, erupted volume and geometry. With respect to the latter, volcanic craters and vents are often highly asymmetrical. Here, we experimentally evaluate the effect of vent asymmetry on gas expansion behaviour and gas jet dynamics directly above the vent. The vent geometries chosen for this study are based on field observations. The novel element of the vent geometry investigated herein is an inclined exit plane (5, 15, 30° slant angle) in combination with cylindrical and diverging inner geometries. In a vertical setup, these modifications yield both laterally variable spreading angles as well as a diversion of the jets, where inner geometry (cylindrical/diverging) controls the direction of the inclination. Both the spreading angle and the inclination of the jet are highly sensitive to reservoir (conduit) pressure and slant angle. Increasing starting reservoir pressure and slant angle yield (1) a maximum spreading angle (up to 62°) and (2) a maximum jet inclination for cylindrical vents (up to 13°). Our experiments thus constrain geometric contributions to the mechanisms controlling eruption jet dynamics with implications for the generation of asymmetrical distributions of proximal hazards around volcanic vents.

Control of Vent Geometry on the Fluid Dynamics of Volcanic Plumes: Insights From Numerical Simulations

AbstractWe present three‐dimensional numerical simulations of eruption clouds from circular to linear fissure vents to investigate the control of vent shape on the height and stability of volcanic plumes during large explosive eruptions. Our results show that clouds ejected from circular or low‐aspect‐ratio (nearly square‐like) fissure vents can be associated with radially suspended flow (RSF) at the top of the jet region, whereas those emitted from narrow‐fissure vents are not. Non‐RSF plumes are more stable than those associated with RSF because the highly concentrated parts of the ejected mixture are easily dissipated and mixed with air near the vent. Plume height in the RSF regime decreases while that in the non‐RSF regime increases with increasing aspect ratio, even for a fixed magma flow rate. These observations suggest that the efficiency of air entrainment is influenced by the vent shape, which in turn controls the dynamics of eruption plumes.

Seismic, morphologic and scale variabilities of subsurface pipes and vent complexes in a magma-rich margin

Subsurface pipes related to fluid flow are observed in seismic reflection data as vertical to sub-vertical features of low-amplitude reflections, linking chaotic reflections at the base (root zones) of volcanic constructs to their summits which can comprise craters, mounds or eye-shaped vents. To date, uncertainties remain regarding the linear correlation and scaling relationship between pipes and their overlying vent structures. Using 3D seismic reflection data from a magma-rich basin, the Vigrid Syncline in the Voring Basin (offshore Norway), this study provides a seismic, morphometric and statistical description of forty (40) pipes and their associated structures. These pipes include conical, bifurcated, fault-controlled and columnar types, which are a consequence of the intrusion of two mappable magmatic sills of Early Eocene age. The heights (av. 1098 m), widths (av. 1740 m) and slenderness ratios (i.e. the heights/widths (Ω)) of the pipes (av. 0.6) show low to moderate correlation coefficients with vent parameters such as height (av. 258 m), area (av. 2.3 km2) and ellipticity (av. 1.5). The morpho-statistical analyses provided here show that ‘pipe’ formation mechanisms are mutually exclusive and that vent geometries are not particularly systematic in nature.

Ice-confined construction of a large basaltic volcano—Austurfjöll massif, Askja, Iceland

Austurfjoll is the largest basaltic glaciovolcanic massif at Askja volcano (Central Iceland), and through detailed studies of its volcanological and geochemical characteristics, we provide a detailed account of the sequence and structure of the ice-confined construction of a large Icelandic basaltic volcano. In particular, Austurfjoll represents a geometry of vents, and resulting glaciovolcanic morphology, not previously documented in ice-confined basaltic volcanoes. Austurfjoll was constructed during two major phases of basaltic volcanism, via seven eruptive episodes through disperse fissure-dominated eruptions. The earliest episode involved a rare and poorly exposed example of subaerial activity, and this was succeeded by six episodes involving the eruption of ice-confined pillow lavas and numerous overlapping fissure eruptions of phreatomagmatic tephra. Evidence of local subaerial lavas and tephras indicates the local growth of eruptive centers above englacial lake levels, and subsequent flooding, but no prolonged subaerial activity. Localized ice-contact facies, paleowater levels, and diamictons indicate the position and thickness of the ice was variable during the construction of Austurfjoll, and eruptive activity likely occurred in multiple and variable level meltwater lakes during the last glacial period. Lithofacies evidence including gradational transitions from effusive to explosive deposits, superposition of fragmental facies above coherent facies, and drainage channels suggest that changes in eruptive style were driven largely by external factors such as drainage and the increasing elevation of the massif. This study emphasizes the unique character of Austurfjoll, being composed of large pillow lava sheets, numerous (> 40) overlapping glaciovolcanic tindars, and only localized emergent deposits, as a product of its prolonged ice-confined eruptive history, contrasts with previous descriptions of tuyas and tindars.

Parameterizing multi-vent activity at Stromboli Volcano (Aeolian Islands, Italy)

The crater terrace of Stromboli Volcano (Italy) hosts several active vents which have evolved and migrated through time within three main vent areas: south-west (SW), central (C), and north-east (NE). Frequent, jet-like explosions typically take place, episodically interrupted by larger-scale paroxysms, which can substantially modify the morphology of the crater terrace and vent geometries. However, the link between the time-space evolution of vent activity and the shallow conduit system are still a matter of debate. In this work, we analyze the vent position and explosion parameters (jet duration and geometry) of 4296 events at Stromboli in five 72-h-long time-windows between 2005 and 2009, as recorded by an infrared surveillance camera. Vent locations illustrate the resilience of the shallow conduit system, which controls explosive activity at different time scales and depths. At the shallowest depth, where slugs burst, conduit branching and merging determines the evolution of simultaneous or alternating twin vents, while vent shape and slug size control local explosion parameters. These processes show variability on an hourly to daily time scale. Below the depth of the slug burst, the conduit system feeding each vent area controls which specific vent will host the explosions and also, possibly, the size of the slugs. Several observations suggest that the C and SW vent areas may be connected at this depth. The deeper conduit system, common to all vent areas, sets the overall explosion rate of the volcano and maintains a balance of this rate between the NE and the combined SW and C vent areas.

Insights into the dynamics of mafic magmatic-hydromagmatic eruptions from volatile degassing behaviour: The Hverfjall Fires, Iceland

The style and intensity of hydromagmatic activity is governed by a complex interplay between the relative volumes of magma and water that interact, their relative viscosities, the depth of subsurface explosions, the substrate properties, and the vent geometry. Fundamental questions remain, however, regarding the role of magmatic vesiculation in determining the dynamics of magma-water interaction (MWI). Petrological reconstructions of magmatic degassing histories are commonly employed to interpret the pre- and syn-eruptive conditions during ‘dry’ magmatic eruptions, but the application of similar techniques to hydromagmatic activity has not yet been fully explored. In this study, we integrate glass volatile measurements (S, Cl, H2O and CO2) with field observations and microtextural measurements to examine the relationship between degassing and eruptive style during the Hverfjall Fires fissure eruption, Iceland. Here, coeval fissure vents produced both ‘dry’ magmatic (Jarðbaðshólar scoria cone complex) and variably wet hydromagmatic (Hverfjall tuff ring) activity, generating physically distinct pyroclastic deposits with contrasting volatile signatures. Matrix glass volatile concentrations in hydromagmatic ash (883 ± 172 [1σ] ppm S; 0.45 ± 0.03 [1σ] wt% H2O; ≤20 ppm CO2) are consistently elevated relative to magmatic ash and scoria lapilli (418 ± 93 [1σ] ppm S; 0.12 ± 0.48 [1σ] wt% H2O; CO2 below detection) and overlap with the range for co-erupted phenocryst-hosted melt inclusions (1522 ± 127 [1σ] ppm S; 165 ± 27 [1σ] ppm Cl). Measurements of hydromagmatic glasses indicate that the magma has degassed between 17 and 70% of its initial sulfur prior to premature quenching at variably elevated confining pressures.By comparing volatile saturation pressures for both magmatic and hydromagmatic glasses, and how these vary through the eruptive stratigraphy, we place constraints on the conditions of MWI. Crucially, our data demonstrate that the magma was already vesiculating when it encountered groundwater at depths of 100–200 m, and that the external water supply was sufficient to maintain MWI throughout the eruption with no significant vertical or lateral migration of the fragmentation surface. We propose that development of an in-vent water-sediment slurry provides a mechanism through which the elevated confining pressures of ~1.6–2.6 MPa (or up to 6 MPa accounting for uncertainty in CO2 below analytical detection) could be maintained and buffered throughout the eruption, whilst enabling vertical mixing and ejection of fragmented juvenile and lithic material from a range of depths. Importantly, these results demonstrate that the volatile contents of hydromagmatic deposits provide valuable records of (1) the environment of MWI (e.g., groundwater versus surface water, vertical migration of the fragmentation level) and (2) the state of the magma at the time of fragmentation and quenching. We further suggest that the volatile content of tephra glasses provides a reliable alternative (or additional) indicator of a hydromagmatic origin, particularly for reduced Ocean Island Basalts where late-stage volatile saturation and degassing (S, H2O) occurs over a pressure range relevant to typical MWI environments.

Magma wagging and whirling in volcanic conduits

Seismic tremor characterized by 0.5–7Hz ground oscillations commonly occur before and during eruptions at silicic volcanoes with widely ranging vent geometries and edifice structures. The ubiquitous characteristics of this tremor imply that its causes are potentially common to silicic volcanoes. Here we revisit and extend to three dimensions the magma-wagging model for tremor (Jellinek and Bercovici, 2011; Bercovici et al., 2013), wherein a stiff magma column rising in a vertical conduit oscillates against a surrounding foamy annulus of bubbly magma, giving rise to tremor. While prior studies were restricted to two-dimensional lateral oscillations, here we explore three-dimensional motion and additional modes of oscillations. In the absence of viscous damping, the magma column undergoes ‘whirling’ motion: the center of each horizontal section of the column traces an elliptical trajectory. In the presence of viscous effect we identify new ‘coiling’ and ‘uncoiling’ column bending shapes with relatively higher and comparable rates of dissipation to the original two-dimensional magma wagging model. We also calculate the seismic P-wave response of the crustal material around the volcanic conduit to the new whirling motions and propose seismic diagnostics for different wagging patterns using the time-lag between seismic stations. We test our model by analyzing pre-eruptive seismic data from the 2009 eruption of Redoubt Volcano. In addition to suggesting that the occurrence of elliptical whirling motion more than 1week before the eruption, our analysis of seismic time-lags also implies that the 2009 eruption was accompanied by qualitative changes in the magma wagging behavior including fluctuations in eccentricity and a reversal in the direction of elliptical whirling motion when the eruption was immediately impending.

Understanding the evolution of a small-volume silicic fissure eruption: Puketerata Volcanic Complex, Taupo Volcanic Zone, New Zealand

The Puketerata Volcanic Complex formed as a result of a series of rhyolitic eruptions along a 2.5km long NE-trending eruptive fissure at the southern margin of Maroa Volcanic Centre (Whakamaru caldera) at 16.5ka. Initial deep-seated explosions were excavated and widen the pre-existing volcanic structures associated to earlier basaltic activity. The following emplacement of two lava domes and associated shallow-seated phreatomagmatic activity produced a widespread tephra blanket characterized by alternating surge and fall dominated pyroclastic units. The silicic magma (73% SiO2) erupted during the Puketerata activity was highly degassed, where explosive fragmentation was driven only by the interaction with external water. The lava domes and ejecta ring makes up more than two-thirds (98×106m3 DRE) of the total eruptive volume (~143×106m3 DRE). The compound ring structures surrounding the larger dome totally overlap the underlying volcanic landforms that formed during the initial phase. The ejecta ring around the larger dome represents a complex landform, which predominantly consists of outwardly thinning pyroclastic breccias forming fan-shape bodies relating to the series of destruction of the larger dome. Despite the relatively large combined dispersal area of the pyroclastic deposits generated by the entire activity (235km2), the volume of single fall beds do not suggest individual explosive events exceeded VEI 2–3. The volcanic hazards associated with such small-volume volcanoes are largely ignored despite these eruptions being more frequent than eruption characterized by caldera collapses. The example of Puketerata suggests that small-volume volcanism often exhibits activity producing fissures, whereas the varied advancement of dykes influenced by external factors makes the course of the eruptions unpredictable in terms of geometry of vents and eruption styles.

Proximal lava drainage controls on basaltic fissure eruption dynamics

Hawaiian basaltic eruptions commonly initiate as a fissure, producing fountains, spattering, and clastogenic lava flows. Most fissures rapidly localize to form a small number of eruptive vents, the location of which may influence the subsequent distribution of lava flows and associated hazards. We present results from a detailed field investigation of the proximal deposits of episode 1 of the 1969 fissure eruption of Mauna Ulu, Kīlauea, Hawai‘i. Exceptional preservation of the deposits allows us to reconstruct vent-proximal lava drainage patterns and to assess the role that drainage played in constraining vent localization. Through detailed field mapping, including measurements of the height and internal depth of lava tree moulds, we reconstruct high-resolution topographic maps of the pre-eruption ground surface, the lava high-stand surface and the post-eruption ground surface. We calculate the difference in elevation between pairs of maps to estimate the lava inundation depth and lava drainage depth over the field area and along different segments of fissure. Aerial photographs collected during episode 1 of the eruption allow us to locate those parts of the fissure that are no longer exposed at the surface. By comparing with the inundation and drainage maps, we find that fissure segments that were inundated with lava to greater depths (typically 1–6 m) during the eruption later became foci of lava drainage back into the fissure (internal drain-back). We infer that, in these areas, lava ponding over the fissure suppressed discharge of magma, thereby favouring drain-back and stagnation. By contrast, segments with relatively shallow inundation (typically less than ~ 1 m), such as where the fissure intersects pre-eruptive topographic highs, or where flow away from the vent (outflow) was efficient, are often associated with sub-circular vent geometries in the post-eruption ground surface. We infer that these parts of the fissure became localization points for ongoing magma ascent and discharge. We conclude that lava inundation and drainage processes in basaltic fissure eruptions can play an important role in controlling their localization and longevity.