Convective tipping point initiates localization of basaltic fissure eruptions

Basaltic fissure eruptions may evolve rapidly and unpredictably complicating hazard management. Localization of an elongate fissure to one or more focused vents may take days to months, and depends on fluid dynamic processes, such as thermally-driven viscous fingering, in the sub-volcanic plumbing system. However, fluid dynamics in a dyke geometry are poorly understood. We perform scaled analogue experiments to investigate convective magma exchange flow within a dyke-like conduit, and discover flow regimes ranging from chaotic mingling to stable, well-organized exchange, over the parameter space relevant for natural eruptions. Experiments are scaled via the Grashof number Gr, which is a Reynolds number for buoyancy-driven exchange flows. We propose that chaotic exchange at high Gr hinders thermally-driven localization by suppressing viscous fingering, whereas flow organization at low Gr enhances localization. Consequently, progressive decrease in Gr through increasing magma viscosity or decreasing dyke width pushes a fissure eruption towards a tipping point that results in rapid localization. Our findings indicate that current conceptual models for magma flow in a dyke require revision to account for this convective tipping point, and provide a quantitative framework for understanding the evolution of fissure eruptions

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')

The use of a shear thinning polymer as a bubbly magma analogue for scaled laboratory experiments

Analogue materials are commonly used in volcanology to perform scaled laboratory experiments. Analogue experiments inform on fundamental fluid dynamic, structural and mechanical processes that are typically very difficult to observe and quantify directly in the natural volcanic system. Here we investigate the suitability of an aqueous solution of hydroxyethyl cellulose polymer (HEC) for use as a lava/magma analogue, with a particular focus on its rheological behaviour. We characterize a range of physical properties as functions of the concentration and temperature of the solution: density; specific heat capacity; thermal diffusivity; thermal conductivity; surface tension; as well as rheology. HEC has a non-Newtonian, shear-thinning rheology that depends on the concentration and temperature of the solution. We show that the rheology is well described by the Cross model, which was originally developed for polymer solutions, but has also been applied to bubbly magmas. Using this similarity, an approach for scaling analogue experiments that use shear-thinning polymers, like HEC, to bubbly magma is presented. A detailed workflow and a spreadsheet are provided to allow experimentalists to investigate the effects of non-Newtonian behaviour in their existing laboratory set-ups. This contribution will allow for the more complex, but often more realistic case of bubble-bearing magmas to be rigorously studied in experimental volcanology

The α–β phase transition in volcanic cristobalite

Cristobalite is a common mineral in volcanic ash produced from dome-forming eruptions. Assessment of the respiratory hazard posed by volcanic ash requires understanding the nature of the cristobalite it contains. Volcanic cristobalite contains coupled substitutions of Al3+ and Na+ for Si4+; similar co-substitutions in synthetic cristobalite are known to modify the crystal structure, affecting the stability of the [alpha] and [beta] forms and the observed transition between them. Here, for the first time, the dynamics and energy changes associated with the [alpha]-[beta] phase transition in volcanic cristobalite are investigated using X-ray powder diffraction with simultaneous in situ heating and differential scanning calorimetry. At ambient temperature, volcanic cristobalite exists in the [alpha] form and has a larger cell volume than synthetic [alpha]-cristobalite; as a result, its diffraction pattern sits between ICDD [alpha]- and [beta]-cristobalite library patterns, which could cause ambiguity in phase identification. On heating from ambient temperature, volcanic cristobalite exhibits a lower degree of thermal expansion than synthetic cristobalite, and it also has a lower [alpha]-[beta] transition temperature (~473 K) compared with synthetic cristobalite (upwards of 543 K); these observations are discussed in relation to the presence of Al3+ and Na+ defects. The transition shows a stable and reproducible hysteresis loop with [alpha] and [beta] phases coexisting through the transition, suggesting that discrete crystals in the sample have different transition temperatures

A model for permeability evolution during volcanic welding

Volcanic ash and pyroclasts can weld when deposited hot by pyroclastic density currents, in near-vent fall deposits, or in fractures in volcano interiors. Welding progressively decreases the permeability of the particle packs, influencing a range of magmatic and volcanic processes, including magma outgassing, which is an important control on eruption dynamics. Consequently, there is a need for a quantitative model for permeability evolution during welding of ash and pyroclasts under the range of conditions encountered in nature. Here we present in situ experiments in which hydrous, crystal-free, glassy pyroclasts are imaged via x-ray tomography during welding at high temperature. For each 3D dataset acquired, we determine the porosity, Darcian gas permeability, specific surface area, and pore connectivity. We find that all of these quantities decrease as a critical percolation threshold is approached. We develop a constitutive mathematical model for the evolution of permeability in welding volcanic systems based on percolation theory, and validate the model against our experimental data. Importantly, our model accounts for polydispersivity of the grainsize in the particle pack, the pressures acting on the pack, and changes in particle viscosity arising from degassing of dissolved H2O during welding. Our model is theoretically grounded and has no fitting parameters, hence it should be valid across all magma compositions. The model can be used to predict whether a cooling pyroclast pack will have sufficient time to weld and to degas, the scenarios under which a final deposit will retain a permeable network, the timescales over which sealing occurs, and whether a welded deposit will have disequilibrium or equilibrium H2O content. A user-friendly implementation of the model is provided

In Vulcan's Forge

Investigations on the correlations between glass and magma involving members of staff from the Glass Department at the University of Sunderland and volcanologists, Dr Fabian Wadsworth and Dr Ed Llewellin from the Natural Sciences Department of Durham University. From a physico-chemical perspective, both magmas and the hot glass manipulated by glass artists are one and the same. The proposed the research question; ‘what lessons can geoscientists learn from knowledge-exchange and experimentation with glass artists?

A reappraisal of explosive–effusive silicic eruption dynamics: syn-eruptive assembly of lava from the products of cryptic fragmentation

Silicic volcanic eruptions range in style from gently effusive to highly explosive, and may switch style unpredictably during a single eruption. Direct observations of subaerial rhyolitic eruptions (Chaiten 2008, Cordón Caulle 2011–2012, Chile) challenged long-standing paradigms of explosive and effusive eruptive styles and led to the formulation of new models of hybrid activity. However, the processes that govern such hybrid explosive–effusive activity remain poorly understood. Here, we bring together observations of the well-studied 2011–2012 Cordón Caulle eruption with new textural and petrologic data on erupted products, and video and still imagery of the eruption. We infer that all of the activity – explosive, effusive, and hybrid – was fed by explosive fragmentation at depth, and that effusive behaviour arose from sticking and sintering, in the shallow vent region, of the clastic products of deeper, cryptic fragmentation. We use a scaling approach to determine that there is sufficient time available, during emplacement, for diffusive pyroclast degassing and sintering to produce a degassed plug that occludes the shallow conduit, feeding clastogenic, apparently effusive, lava-like deposits. Based on evidence from Cordón Caulle, and from other similar eruptions, we further argue that hybrid explosive–effusive activity is driven by episodic gas-fracking of the occluding lava plug, fed by the underlying pressurized ash- and pyroclast-laden region. The presence of a pressurized pocket of ash-laden gas within the conduit provides a mechanism for generation of harmonic tremor, and for syn-eruptive laccolith intrusion, both of which were features of the Cordón Caulle eruption. We conclude that the cryptic fragmentation models is more consistent with available evidence than the prevailing model for effusion of silicic lava that assume coherent non-fragmental rise of magma from depth to the surface without wholesale explosive fragmentation

Experimental observations of bubbling regimes at in-line multi-orifice bubblers

Bubble formation and bubbling regimes are well-characterized for the cases of single-orifice bubblers and industrial perforated plates. However, bubbling regimes from bubblers with multiple in-line orifices remain poorly described. Here, we investigate the dynamics of bubble formation at both single-orifice and multi-orifice bubblers, with one, three, five and nine in-line orifices in an 80-cm-long bubbler. We use high-speed videography and image processing to identify the effects of bubbler volume, and the number, spacing, and diameter of orifices, on bubbling regimes, bubble period, and bubble formation time. We identify five main bubbling regimes based on synchronization among orifices, and discuss the parameters affecting the bubbling dynamics. Decreasing bubbler volume leads to a decrease in bubble volume and bubble period, and enhances synchronization. Increasing orifice diameter leads to an increase in bubble volume and enhances synchronization. Spacing between orifices doesn't play an important role in determining the bubbling regime. Based on the experimental observations, we develop new bubbling regime maps constructed using the dimensionless Capacitance number and Weber number