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.Postgraduate scholarship from University of Bristol AXA Research Fund Royal Society Wolfson Merit Award Royal Society University Research Fellowship New Researchers Award from the Geologists’ Associatio

Contrasting mechanisms of magma fragmentation during coeval magmatic and hydromagmatic activity: the Hverfjall Fires fissure eruption, Iceland

Growing evidence for significant magmatic vesiculation prior to magma-water interaction (MWI) has brought into question the use of ‘diagnostic’ features, such as low vesicularities and blocky morphologies, to identify hydromagmatic pyroclasts. We address this question by quantifying co-variations in particle size, shape and texture in both magmatic and hydromagmatic deposits from the Hverfjall Fires fissure eruption, Iceland. Overlapping vesicularity and bubble number density distributions measured in rapidly quenched magmatic and hydromagmatic pyroclasts indicate a shared initial history of bubble nucleation and growth, with substantial vesiculation prior to MWI. Hydromagmatic fragmentation occurred principally by brittle mechanisms, where the length scale and geometry of fracturing was controlled by the bubble population. This suggests that the elevated fragmentation efficiency of hydromagmatic deposits is driven, at least in part, by brittle disintegration of vesicular pyroclasts due to high thermal stress generated during rapid cooling. In this way, the shape and size distributions of hydromagmatic pyroclasts, both critical input parameters for ash dispersion models, are strongly influenced by the dynamics of vesiculation prior to MWI. This result underlines the need to analyse multiple grain-size fractions to characterise the balance between magmatic and hydromagmatic processes. During the Hverfjall Fires eruption, the external water supply was sufficient to maintain MWI throughout the eruption, with no evidence for progressive exhaustion of a water reservoir. We suggest that both the longevity and the spatial distribution of MWI were determined by the pre-existing regional hydrology and represent continuous interaction between a propagating dike and a strong groundwater flow system hosted within permeable basalt lavas

Open-vent volcanoes fuelled by depth-integrated magma degassing

AbstractOpen-vent, persistently degassing volcanoes—such as Stromboli and Etna (Italy), Villarrica (Chile), Bagana and Manam (Papua New Guinea), Fuego and Pacaya (Guatemala) volcanoes—produce high gas fluxes and infrequent violent strombolian or ‘paroxysmal’ eruptions that erupt very little magma. Here we draw on examples of open-vent volcanic systems to highlight the principal characteristics of their degassing regimes and develop a generic model to explain open-vent degassing in both high and low viscosity magmas and across a range of tectonic settings. Importantly, gas fluxes from open-vent volcanoes are far higher than can be supplied by erupting magma and independent migration of exsolved volatiles is integral to the dynamics of such systems. The composition of volcanic gases emitted from open-vent volcanoes is consistent with its derivation from magma stored over a range of crustal depths that in general requires contributions from both magma decompression (magma ascent and/or convection) and iso- and polybaric second boiling processes. Prolonged crystallisation of water-rich basalts in crustal reservoirs produces a segregated exsolved hydrous volatile phase that may flux through overlying shallow magma reservoirs, modulating heat flux and generating overpressure in the shallow conduit. Small fraction water-rich melts generated in the lower and mid-crust may play an important role in advecting volatiles to subvolcanic reservoirs. Excessive gas fluxes at the surface are linked to extensive intrusive magmatic activity and endogenous crustal growth, aided in many cases by extensional tectonics in the crust, which may control the longevity and activity of open-vent volcanoes. There is emerging abundant geophysical evidence for the existence of a segregated exsolved magmatic volatile phase in magma storage regions in the crust. Here we provide a conceptual picture of gas-dominated volcanoes driven by magmatic intrusion and degassing throughout the crust.</jats:p

Contrasting mechanisms of magma fragmentation during coeval magmatic and hydromagmatic activity: the Hverfjall Fires fissure eruption, Iceland

Growing evidence for significant magmatic vesiculation prior to magma-water interaction (MWI) has brought into question the use of ‘diagnostic’ features, such as low vesicularities and blocky morphologies, to identify hydromagmatic pyroclasts. We address this question by quantifying co-variations in particle size, shape and texture in both magmatic and hydromagmatic deposits from the Hverfjall Fires fissure eruption, Iceland. Overlapping vesicularity and bubble number density distributions measured in rapidly quenched magmatic and hydromagmatic pyroclasts indicate a shared initial history of bubble nucleation and growth, with substantial vesiculation prior to MWI. Hydromagmatic fragmentation occurred principally by brittle mechanisms, where the length-scale and geometry of fracturing was controlled by the bubble population. This suggests that the elevated fragmentation efficiency of hydromagmatic deposits is driven, at least in part, by brittle disintegration of vesicular pyroclasts due to high thermal stress generated during rapid cooling. In this way, the shape and size distributions of hydromagmatic pyroclasts, both critical input parameters for ash dispersion models, are strongly influenced by the dynamics of vesiculation prior to MWI. This result underlines the need to analyse multiple grain-size fractions to characterise the balance between magmatic and hydromagmatic processes. During the Hverfjall Fires eruption, the external water supply was sufficient to maintain MWI throughout the eruption, with no evidence for progressive exhaustion of a water reservoir. We suggest that both the longevity and the spatial distribution of MWI was determined by the pre-existing regional hydrology, and represents continuous interaction between a propagating dike and a strong groundwater flow system hosted within permeable basalt lavas.University of Bristol postgraduate scholarship New Researchers Award, Geologists’ Association AXA Research Grant Royal Society Wolfson Merit Award Royal Society UR

Petrologic monitoring at Volcán de Fuego, Guatemala

Paroxysmal activity represents an end-member in the common range of activity at mafic arc volcanoes, characterised by rapid transitions across the effusive-explosive interface and thus posing significant challenges to hazard assessment. Conceptual models to explain changes in the frequency and magnitude of these paroxysmal events are based either on magma recharge or an increase in gas flux, largely framed in the context of two-phase flow. Gas- and magma-driven models are both viable mechanisms to explain the varying styles of paroxysmal behaviour observed in mafic systems; however, each has different implications for future activity. We present time series petrologic data for ash and lava samples collected at Volcán de Fuego, Guatemala, during paroxysmal eruptions between 2011 and 2018. We show that a step-change in glass composition occurred between 2015 and 2016, reflecting an increase in magma temperature and a reduction in pre-eruptive crystallisation, concurrent with an escalation in the frequency of paroxysmal activity. There was no change in the bulk or phase compositions during this period. To explain these observations, we propose that the increase in frequency of paroxysmal eruptions is modulated by the supply of exsolved volatiles from lower crustal degassing magmas, without invoking repeated transfer of new, primitive magma to a shallow reservoir. Protracted lava effusion, accompanied by more vigorous and more frequent Strombolian explosions and gas ‘chugging’, prior to the transition to sustained fountaining suggests that gas retention in crystal-rich magma may modulate the height of the magma column as gas supply increases. Slow decompression associated with effusion may determine the timing of effusive to explosive transitions in mafic arc systems more generally. A large paroxysmal eruption of Fuego on 3 June 2018, notable for the rapid escalation in eruptive intensity several hours into the eruption, produced ash with a range of textures and glass compositions consistent with magma evacuation over a range of depths and decompression rates. Given the protracted repose time between paroxysms before this event, we suggest that a shallow crystallised plug degraded, and ultimately failed, several hours into the eruption of 3 June 2018, triggering top-down decompression of magma in the conduit synchronous with the observed rapid acceleration in eruption rate. Ultimately, we propose that the frequency of paroxysms at Fuego is broadly proportional to the gas supply rate, while the range in glass compositions is related to the repose time prior to eruptive activity. Our data illustrate the potential of petrologic monitoring to distinguish between gas- and magma-driven paroxysm triggers and to anticipate future events, especially when interpreted in the context of geophysical observations and implemented within community-based ash collection initiatives

Crystal and melt inclusion timescales reveal the evolution of magma migration before eruption

Volatile element concentrations measured in melt inclusions are a key tool used to understand magma migration and degassing, although their original values may be affected by different re-equilibration processes. Additionally, the inclusion-bearing crystals can have a wide range of origins and ages, further complicating the interpretation of magmatic processes. To clarify some of these issues, here we combined olivine diffusion chronometry and melt inclusion data from the 2008 eruption of Llaima volcano (Chile). We found that magma intrusion occurred about 4 years before the eruption at a minimum depth of approximately 8 km. Magma migration and reaction became shallower with time, and about 6 months before the eruption magma reached 3–4 km depth. This can be linked to reported seismicity and ash emissions. Although some ambiguities of interpretation still remain, crystal zoning and melt inclusion studies allow a more complete understanding of magma ascent, degassing, and volcano monitoring data.NRF (Natl Research Foundation, S’pore)MOE (Min. of Education, S’pore)Published versio

Textural and geochemical constraints on andesitic plug emplacement prior to the 2004–2010 vulcanian explosions at Galeras volcano, Colombia

Hazardous sequences of vulcanian explosions are thought to result from the repeated emplacement and destruction of degassed, highly crystalline magma plugs in the shallow conduit of arc volcanoes. The processes governing the timing and magnitude of these explosions are thought to be related to magma ascent rate and efficiency of degassing and crystallisation. We study a rare suite of time-constrained ballistic bombs from the 2004–2010 period of activity of Galeras volcano to reconstruct magma plug architecture prior to six individual vulcanian explosions. We find that each plug was vertically stratified with respect to crystallinity, vesicularity and melt volatile content, melt composition and viscosity. We interpret this structure as resulting from multiple bubble nucleation events and degassing-driven crystallisation during multi-step ascent of the magma forming the plug, followed by spatially variable crystallisation within the plug under contrasting conditions of effective undercooling created by degassing. We propose that the shallow conduit evolved from more open degassing conditions during 2004–2008 to more closed conditions during 2009–2010. This resulted in explosions becoming smaller and less frequent over time during 2004–2008, then larger and more frequent over time during 2009–2010. This evolution was controlled by changing average ascent rates and is recorded by systematic changes in plagioclase microlite textures. Our results suggest that small volume vulcanian explosions (~ 105 m3) should generally be associated with longer repose times (hundreds of days) and produce ballistics characterised by small numbers of large, prismatic plagioclase microlites. Larger volume vulcanian explosions (~ 106 m3) should be associated with shorter repose times (tens of days) and produce ballistics characterised by high numbers of small, more tabular plagioclase microlites

Transitions between explosive and effusive phases during the cataclysmic 2010 eruption of Merapi volcano, Java, Indonesia

Transitions between explosive and effusive activity are commonly observed during dome-forming eruptions and may be linked to factors such as magma influx, ascent rate and degassing. However, the interplay between these factors is complex and the resulting eruptive behaviour often unpredictable. This paper focuses on the driving forces behind the explosive and effusive activity during the well-documented 2010 eruption of Merapi, the volcano’s largest eruption since 1872. Time-controlled samples were collected from the 2010 deposits, linked to eruption stage and style of activity. These include scoria and pumice from the initial explosions, dense and scoriaceous dome samples formed via effusive activity, as well as scoria and pumice samples deposited during subplinian column collapse. Quantitative textural analysis of groundmass feldspar microlites, including measurements of areal number density, mean microlite size, crystal aspect ratio, groundmass crystallinity and crystal size distribution analysis, reveal that shallow pre- and syn-eruptive magmatic processes acted to govern the changing behaviour during the eruption. High-An (up to ∼80 mol% An) microlites from early erupted samples reveal that the eruption was likely preceded by an influx of hotter or more mafic magma. Transitions between explosive and effusive activity in 2010 were driven primarily by the dynamics of magma ascent in the conduit, with degassing and crystallisation acting via feedback mechanisms, resulting in cycles of effusive and explosive activity. Explosivity during the 2010 eruption was enhanced by the presence of a ‘plug’ of cooled magma within the shallow magma plumbing system, which acted to hinder degassing, leading to overpressure prior to initial explosive activity

A transmission electron microscope study of white mica crystallite size distribution in a mudstone to slate transitional sequence, North Wales, UK

High-resolution transmission electron microscopy (HRTEM) measurements of the thickness of white mica crystallites were made on three pelite samples that represented a prograde transition from diagenetic mudstone though anchizonal slate to epizonal slate. Crystallite thickness, measured normal to (001), increases as grade increases, whereas the XRD measured 10 Å peak-profile, the Kubler index, decreases. The mode of the TEM-measured size population can be correlated with the effective crystallite size N (001) determined by XRD. The results indicate that the Kubler index of white mica crystallinity measures changes in the crystallite size population that result from prograde increases in the size of coherent X-ray scattering domains. These changes conform to the Scherrer relationship between XRD peak broadening and small crystallite size. Lattice ‘strain’ broadening is relatively unimportant, and is confined to white mica populations in the diagenetic mudstone. Rapid increases in crystallite size occur in the anchizone, coincident with cleavage development. Changes in the distribution of crystallite thickness with advancing grade and cleavage development are characteristic of grain-growth by Ostwald ripening. The Kubler index rapidly loses sensitivity as an indicator of metapelitic grade within the epizone.Peer Reviewedhttp://deepblue.lib.umich.edu/bitstream/2027.42/47293/1/410_2004_Article_BF00306406.pd