Search

This document summarizes the outcomes of the Modeling Collaboratory for Subduction Zone Science (MCS) Volcanic Systems Workshop and presents a vision for advancing collaborative modeling of volcano-magmatic systems. The U.S. Geological Survey (USGS) has identified 161 potentially active volcanoes in the United States and its territories, of which 57 are considered to be high or very high threats (Ewert et al., 2018). All western states, including Alaska and Hawaii, have potentially active volcanoes. Eruptions range from the quiet effusion of sluggish lava flows over hours to decades to immense explosive ejections of tephra which produce massive calderas.Understanding these volcanoes and assessing their threat to society requires the development of quantitative models, rooted in physics and chemistry, which can be used to interpret diverse observations including real-time monitoring data. Existing models have tremendously advanced our understanding of volcanic systems and have improved our ability to assess hazards and forecast future activity, contributing directly to reductions in the number of lives lost to volcanic eruptions and helping mitigate their costs to society. Magmatic system models also provide a quantitative framework for understanding processes that occur at depth beneath volcanoes, linking volcanic systems with a broad range of deeper processes associated with the production, transport, and storage of magma and associated fluids above subducting slabs.Despite this exciting progress much remains to be accomplished and workshop participants identified several important opportunities. First and foremost is the recognition that enhanced support for the development and dissemination of volcano-magmatic system models and associated methodologies will enable advances in ways not currently possible. A key outcome of the workshops is a recognition of the transformative potential of diverse groups of scientists working together on common problems. Support for collaborative working groups will enable communication across disciplines and between modelers and non-modelers, leveraging expertise from scientists studying different aspects of volcano-magmatic systems, and between geoscientists and outside experts from fields such as mathematics, statistics, and material sciences. Better support will also enable modelers to more fully verify, validate, benchmark, and document their codes, and also provide new training opportunities. Enhanced model sharing and interoperability will reduce the need for different groups to independently duplicate (re-invent) code and increase confidence in published results. This report lays out a proposal for a collaborative modeling environment that is centered in large part around community working groups manifested as workshops, summer schools, and sustained long-term research collaborations involving diverse groups of scientists working on common problems. Programmatic support is envisioned in the form of enhanced student and postdoc funding for model development, incentives and support for cross-disciplinary collaborative research projects, and related support for these activities. This support will fundamentally improve our ability to integrate and interpret observations using volcanic and magmatic system models.

NASA’s rapid response plan for gathering atmospheric data amid major volcanic eruptions, paired with efforts to improve eruption simulations, will offer better views of these events’ global effects.

Magma from Plinian volcanic eruptions contains an extraordinarily large numbers of bubbles. Nucleation of those bubbles occurs because pressure decreases as magma rises to the surface. As a consequence, dissolved magmatic volatiles, such as water, become supersaturated and cause bubbles to nucleate. At the same time, diffusion of volatiles into existing bubbles reduces supersaturation, resulting in a dynamical feedback between rates of nucleation due to magma decompression and volatile diffusion. Because nucleation rate increases with supersaturation, bubble number density (BND) provides a proxy record of decompression rate, and hence the intensity of eruption dynamics. Using numerical modeling of bubble nucleation, we reconcile a long-standing discrepancy in decompression rate estimated from BND and independent geospeedometers. We demonstrate that BND provides a record of the time-averaged decompression rate that is consistent with independent geospeedometers, if bubble nucleation is heterogeneous and facilitated by magnetite crystals., The authors simulate bubble nucleation in silica-rich magma with conditions appropriate for Plinian eruptions. They demonstrate that the gap between decompression rate estimates from bubble number density and independent geospeedometers can be largely closed if nucleation is heterogenous facilitated by magnetite crystals and decompression rate is calculated as time-averaged values.

Bubble nucleation is the first step in magma degassing. The kinetics of bubble nucleation during eruptive magma ascent determines the rate of magma degassing, which in turn, has a significant control on the explosivity of eruptions. Here, we examine if bubble nucleation in rhyolite can be continuous over the inferred timescale of magma travel time from reservoir to the surface in explosive silicic eruptions. We performed homogeneous bubble nucleation experiments using cylinders cored from natural rhyolite. Samples were hydrated at an initial pressure and temperature, then decompressed to a final pressure, and held at that pressure for various amounts of time before being quenched. Within a given suite samples had the same temperature as well as initial and final pressures, but were held for different annealing times. The resulting bubble number density of samples within each suite increases by up to 2 orders of magnitude as annealing time increases, indicating that nucleation continued during annealing. We estimate the nucleation rate and duration from measured bubble number densities and show that nucleation in the majority of samples was still ongoing at the time of quenching. Our findings confirm theoretical predictions that bubble nucleation can be a continuous process during eruptive magma ascent in explosive eruptions.

The fluxes of volatile elements from the mantle have long been used to understand mantle structure and evolution, and are critical controls on Earth's climate stability. Because of the ubiquity of magmatic degassing, inferring pre-degassing volatile concentrations from measured basalts requires the application of a degassing model. Such models, including the commonly-applied equilibrium Rayleigh distillation, typically assume equilibrium or solubility-based partitioning between melt and vapor. Here, we demonstrate that ratios of radiogenic isotopes of He, Ne, Ar and, especially, Xe measured in global mid-ocean ridge basalts (MORBs) are inconsistent with equilibrium degassing models, even when not considering He. We conclude that kinetic disequilibrium is a crucial process affecting volatile abundances during degassing. We present a simple disequilibrium Rayleigh distillation model to reconstruct pre-degassing MORB noble gas and carbon concentrations, which predicts that He and Ne achieve nearly equilibrium partitioning between melt and vesicles, but the slower-diffusing heavier noble gases are strongly affected by disequilibrium, resulting in non-equilibrium fractionation of noble gas elemental ratios, and other volatile element ratios like CO2/He. We apply our model to a large set of MORB data, and find average pre-degassing 3He, 22Ne, and 36Ar concentrations of 4.4 ± 0.9 × 10−10, 6.6 ± 1.4 × 10−11, and 6.8 ± 4.5 × 10−10 ccSTP/g (2σ), with variations of approximately 2 orders of magnitude, similar to other highly incompatible elements. Pre-degassing noble gas concentrations imply a mid-ocean ridge 3He flux of 800 ± 170 mol/yr and upper mantle 3He/22Ne and 3He/36Ar ratios of 6.6 ± 2.0 and 0.64 ± 0.44, but substantial variability in these ratios between samples. Applying our model to CO2, we calculate an average mantle CO2/3He molar ratio of 1.67 ± 0.21 × 109, which, when combined with our estimate of 3He flux, implies an upper mantle CO2 flux of 5.9 ± 1.0 × 1013 g/yr and a CO2 concentration of approximately 110 ppm. Our estimate of the mantle 3He flux is the first determined independently of oceanographic 3He measurements, and consequently represents a time-integrated flux substantially longer than ∼1000 years. And although at the high end of the range of previous estimates, our estimate does not resolve the long-standing heat-helium paradox. Additionally, our requirement for heterogeneous pre-degassing 3He/22Ne and 3He/36Ar ratios between samples is contrary to conclusions of previous applications of disequilibrium degassing models, which advocated uniform ratios. Furthermore, we find that CO2/Ba ratios are highly variable in MORB samples, but still consistent with an average mantle mass ratio of ∼100. However, estimating pre-degassing CO2 concentrations and the mantle CO2 flux depend strongly on the poorly-constrained carbon diffusivity. Consequently, our demonstration of the prevalence of disequilibrium during mid-ocean ridge degassing, and the potential for disequilibrium in other volcanic settings, highlights the need for better characterization of the physical parameters associated with volcanic degassing.

Quantitative metrics for eruption rates at mid-ocean ridges (MORs) would improve our understanding of the structure and formation of the uppermost oceanic crust and would provide a means to link volcanic processes with the conditions of the underlying magmatic system. However, these metrics remain elusive because no MOR eruptions have been directly observed. The possibility of disequilibrium degassing in mid-ocean ridge basalts (MORB), due to high eruptive depressurization rates, makes the analysis of volatile concentrations in MORB glass a promising method for evaluating eruption rates. In this study, we estimate magma ascent and lava flow emplacement rates during the 2011 eruption of Axial Seamount based on numerical modeling of diffusion-controlled bubble growth and new measurements of dissolved volatiles, vesicularity, and vesicle size distributions in erupted basalts. This dataset provides a unique view of the variability in magma ascent (∼0.02–1.2 m/s) and lava flow rates (∼0.1–0.7 m/s) during a submarine MOR eruption based on 50 samples collected from a >10 km long fissure system and three individual lava flow lobes. Samples from the 2011 eruption display an unprecedented range in dissolved CO2 concentrations, nearly spanning the full range observed on the global MOR system. The variable vesicularity and dissolved CO2 concentrations in these samples can be explained by differences in the extent of degassing, dictated by flow lengths and velocities during both vertical ascent and horizontal flow along the seafloor. Our results document, for the first time, the variability in magma ascent rates during a submarine eruption (∼0.02–1.2 m/s), which spans the global range previously proposed based on CO2 degassing. The slowest ascent rates are associated with hummocky flows while faster ascent rates produce channelized sheet flows. This study corroborates degassing-based models for eruption rates using comparisons with independent methods and documents the relationship between eruption dynamics, magma ascent rates, and the morphology of eruptive products. Globally, this approach allows interrogation of the processes that govern mid-ocean ridge eruptions and influence the formation of the oceanic crust.

Surface tension of hydrous rhyolitic melt is high enough that large degrees of supersaturation are needed to homogeneously nucleate H2O bubbles during eruptive magma ascent. This study examines whether dissolved fluorine lowers surface tension of hydrous rhyolite, and thus lowers the supersaturation required for bubble nucleation. Fluorine was targeted because it, like H2O, changes melt properties and is highly soluble, unlike all other common magmatic volatiles. Rhyolite melts were saturated at Ps = 245 MPa with H2O fluid that contained F, generating rhyolite with 6.7 ± 0.4 wt.% H2O and 1.1–1.3 wt.% F. When these melts were decompressed rapidly to Pf = 149–202 MPa and quenched after 60 s, bubbles nucleated at supersaturations of ΔP = Ps − Pf ≥52 MPa, and reached bubble number densities of NB = 1012–13 m−3 at ΔP = 78–101 MPa. In comparison, rhyolite saturated with 6.34 ± 0.09 wt.% H2O, but only 0.25 wt.% F, did not nucleate bubbles until ΔP ≥ 100–116 MPa, and even then, at significantly lower NB (

Plinian lapilli from the 1060 Common Era Glass Mountain rhyolitic eruption of Medicine Lake Volcano, California, were collected and analyzed for vesicularity and permeability. A subset of the samples were deformed at a temperature of 975°, under shear and normal stress, and postdeformation porosities and permeabilities were measured. Almost all undeformed samples fall within a narrow range of vesicularity (0.7–0.9), encompassing permeabilities between approximately 10−15 m2 and 10−10 m2. A percolation threshold of approximately 0.7 is required to fit the data by a power law, whereas a percolation threshold of approximately 0.5 is estimated by fitting connected and total vesicularity using percolation modeling. The Glass Mountain samples completely overlap with a range of explosively erupted silicic samples, and it remains unclear whether the erupting magmas became permeable at porosities of approximately 0.7 or at lower values. Sample deformation resulted in compaction and vesicle connectivity either increased or decreased. At small strains permeability of some samples increased, but at higher strains permeability decreased. Samples remain permeable down to vesicularities of less than 0.2, consistent with a potential hysteresis in permeability-porosity between expansion (vesiculation) and compaction (outgassing). We attribute this to retention of vesicle interconnectivity, albeit at reduced vesicle size, as well as bubble coalescence during shear deformation. We provide an equation that approximates the change in permeability during compaction. Based on a comparison with data from effusively erupted silicic samples, we propose that this equation can be used to model the change in permeability during compaction of effusively erupting magmas.

The East Antarctic ice sheet was larger than present during past cold periods. Seafloor geophysical data show that in the Ross Sea, the extended ice sheet was underlain by an active hydrologic system during the glacial termination.

International audience; Coalescence during bubble nucleation and growth in crystal‐free rhyolitic melt was experimentally investigated, and the percolation threshold, defined as the porosity at which the vesicular melt first becomes permeable, was estimated. Experiments with bubble number densities between 1014 and 1015 m−3 were compared to four suites of rhyolitic Plinian pumices, which have approximately equal bubble number densities. At the same total porosity, Plinian samples have a higher percentage of coalesced bubbles compared to their experimental counterparts. Percolation modeling of the experimental samples indicates that all of them are impermeable and have percolation thresholds of approximately 80–90%, irrespective of their porosity. Percolation modeling of the Plinian pumices, all of which have been shown to be permeable, gives a percolation threshold of approximately 60%. The experimental samples fall on a distinct trend in terms of connected versus total porosity relative to the Plinian samples, which also have a greater melt‐bubble structural complexity. The same holds true for experimental samples of lower bubble number densities. We interpret the comparatively higher coalescence within the Plinian samples to be a consequence of shear deformation of the erupting magma, together with an inherently greater structural complexity resulting from a more complex nucleation process.

Eruption rates at the mid-ocean ridges (MORs) are believed to strongly control the morphology and length of lava flows emplaced along the ridge axis, and thus the structure and porosity of the upper oceanic crust. Eruption rate also represents one of the few tools to gain insight into the driving pressures within sub-ridge magmatic systems. As eruption rate is inferred to vary systematically along the global mid-ocean ridge system, understanding of how to assess eruption rate in submarine systems and how it maps to observable features of the ridge axis would provide a powerful tool to understand Earth's largest magmatic system. Eruption rates at MORs are poorly constrained, however, because of a lack of direct observations, preventing the duration of an eruption to be quantified. This study uses decompression experiments of MORB samples and numerical modeling of CO 2 degassing to reconstruct the timescales for magma ascent and lava emplacement during the 2005–06 eruption of the East Pacific Rise at ca. 9°51′N. Samples collected from the lava flow are all supersaturated in dissolved CO 2 contents, but CO 2 decreases with distance from the vent, presumably as a consequence of progressive CO 2 diffusion into growing bubbles. Samples collected at the vent contain ∼10 5 vesicles per cm 3 . Pieces of these samples were experimentally heated to 1225 °C at high pressure and then decompressed at controlled rates. Results, plus those from numerical modeling of diffusive bubble growth, indicate that magma rose from the axial magma chamber to the seafloor in ≤1 h and at a rate of ≥2–3 km h −1 . Our modeling, as validated by experimental decompression of MORB samples with ∼10 6 vesicles cm −3 , also suggests that CO 2 degassed from the melt within ∼10–100 min as the vesicular lava traveled ∼1.7 km along the seafloor, implying a volumetric flow rate on order of 10 3–4 m 3 s −1 . Given an ascent rate of ≥0.2 m s −1 , the width of a rectangular dike feeding the lava would have been ∼1–2 m wide. MORB samples from the Pacific ridge are generally more supersaturated in dissolved CO 2 than those from slower spreading Atlantic and Indian ridges. Our results suggest that Pacific MORBs ascend to the seafloor faster than Atlantic or Indian MORBs.

Sinters are siliceous, sedimentary deposits that form in geothermal areas. Formation occurs in two steps. Hot water circulates in the subsurface and dissolves silica from the host rock, usually rhyolites. Silica then precipitates after hot water is discharged and cools. Extensive sinter formations are linked to up-flow areas of fluids originating from high temperature (> 175 °C) deep reservoirs. Fluid geochemistry, microbial communities, and environmental conditions of deposition determine the texture of sinter and pore framework. Porosity strongly influences physical and hydraulic properties of rocks. To better understand the properties controlling the transport of fluids, and interpret geophysical observations in geothermal systems, we studied 17 samples of modern geyserite sinter deposits (

Summary The amendment of soil with biochar can sequester carbon and alter hydrologic properties by changing physical and chemical characteristics of soil. To understand the effect of biochar amendment on soil hydrology, we measured the hydraulic conductivity ( K ) of biochar–sand mixtures as well as dissolved organic carbon (DOC) in leachate. Specifically, we assessed the effects of biochar concentration and particle size on K and amount of DOC in the soil leachate. To better understand how physical properties influenced K , we also measured the skeletal density of biochars and sand, and the bulk density, the water saturation, and the porosity of biochar–sand mixtures. Our model soil was sand (0.251–0.853 mm) with biochar rates from 2 to 10 wt% (g biochar/g total soil × 100%). As biochar ( K decreased by 72 ± 3%. When biochar particle size was equal to, greater than, and less than particle size of sand, we found that biochar in different particle sizes have different effects on K . For a 2 wt% biochar rate, K decreased by 72 ± 2% when biochar particles were finer than sand particles, and decreased by 15 ± 2% when biochar particles were coarser than sand particles. When biochar and sand particle size were comparable, we observed no significant effect on K . We propose that the decrease of K through the addition of fine biochar was because finer biochar particles filled spaces between sand particles, which increased tortuosity and reduced pore throat size of the mixture. The decrease of K associated with coarser biochar was caused by the bimodal particle size distribution, resulting in more compact packing and increased tortuosity. The loss of biochar C as DOC was related to both biochar rate and particle size. The cumulative DOC loss was 1350% higher from 10 wt% biochar compared to pure sand. This large increase reflected the very small DOC yield from pure sand. In addition, DOC in the leachate decreased as biochar particle size increased. For all treatments, the fraction of carbon lost as DOC ranged from 0.06 to 0.18 wt% of biochar. These experiments suggest that mixing sandy soils with biochar is likely to reduce infiltration rates, holding water near the surface longer with little loss of biochar-derived carbon to groundwater and streams.

The weakest explosive volcanic eruptions globally, Strombolian explosions and Hawaiian fountaining, are also the most common. Yet, despite over a hundred years of observations, no classifications have offered a convincing, quantitative way of demarcating these two styles. New observations show that the two styles are distinct in their eruptive time scale, with the duration of Hawaiian fountaining exceeding Strombolian explosions by ∼300–10,000 s. This reflects the underlying process of whether shallow-exsolved gas remains trapped in the erupting magma or is decoupled from it. We propose here a classification scheme based on the duration of events (brief explosions versus prolonged fountains) with a cutoff at 300 s that separates transient Strombolian explosions from sustained Hawaiian fountains.

Despite the strong theoretical basis for the importance of magma decompression rate in modulating the dynamics of explosive basaltic eruptions, few measurements exist. Here, we quantify magma decompression rate using the incomplete diffusive re-equilibration of water through olivine for two contrasting eruptions at Cerro Negro volcano in Nicaragua: the 1992 VEI 3 and 1995 VEI 2 eruptions. Both eruptions have the same basaltic composition (SiO2 ∼ 50 wt%) and maximum volatile concentrations (H2O ∼ 4.7 wt%). SIMS and NanoSIMS measurements of water in olivine-hosted melt inclusions (MIs) from both eruptions show a negative correlation between MI size and water loss, as predicted for diffusive re-equilibration. MIs from the 1995 eruption appear to have experienced more water loss compared to those from the 1992 eruption (up to 80–90% water loss in MIs Quantifying the timescale for water diffusion during magma decompression requires an accurate constraint on H+ diffusivity through olivine. We perform dehydration experiments on oriented cuboids of olivine from Cerro Negro in a 1 atm. furnace at 800 and 1000 °C and at an oxygen fugacity 0.25 log units above the nickel-nickel oxide buffer (NNO + 0.25) in a sequence of heating steps. Between each heating step we measure the H+ concentration along each crystallographic direction by FTIR and model these profiles to determine the diffusivity. We find the fastest diffusion direction along a and develop an Arrhenius relationship for the diffusivity along this direction: D a ( m 2 / s ) = 9.6 × 1 0 − 6 exp − 125 R T , where R is the gas constant 8.314 J/mol K; T is the temperature in K and the activation energy is 125 kJ/mol. These are the fastest rates ever measured for olivine dehydration, which may relate to the low forsterite content of the olivines (Fo ∼ 79). We have developed a 1D model for H+ diffusion along the fast direction since this dominates water loss (Da > 10 × Db or Dc). Using our new H+ diffusivity, we model water diffusion from MIs during magma decompression to determine ascent timescales. We find that magma from the VEI 3 eruption of 1992 ascended at a median rate of 0.007 MPa/s, which is ∼5 × faster than the median for the 1995 VEI 2 eruption, although both eruptions show a ∼ two orders of magnitude spread in decompression rate. Assuming linear decompression, this places a minimum bound of ∼6.5 h on magma ascent during the 1992 eruption and ∼35 h for the 1995 eruption. Monte Carlo simulations give a 2σ uncertainty on log 10(dP/dt) of ∼ 10% — much smaller than the order of magnitude spread in decompression rate within each eruption, which suggests that there is real variation in decompression rate for a given deposit. Our results support the hypothesis that higher decompression rates are correlated with more explosive eruptions and show that diffusion modelling of MI water loss is a powerful tool that is sensitive to changes in decompression between VEI 2 and 3 eruption intensities.

Magma fragmentation is the breakup of a continuous volume of molten rock into discrete pieces, called pyroclasts. Because magma contains bubbles of compressible magmatic volatiles, decompression of low-viscosity magma leads to rapid expansion. The magma is torn into fragments, as it is stretched into hydrodynamically unstable sheets and filaments. If the magma is highly viscous, resistance to bubble growth will instead lead to excess gas pressure and the magma will deform viscoelastically by fracturing like a glassy solid, resulting in the formation of a violently expanding gas-pyroclast mixture. In either case, fragmentation represents the conversion of potential energy into the surface energy of the newly created fragments and the kinetic energy of the expanding gas-pyroclast mixture. If magma comes into contact with external water, the conversion of thermal energy will vaporize water and quench magma at the melt-water interface, thus creating dynamic stresses that cause fragmentation and the release of kinetic energy. Lastly, shear deformation of highly viscous magma may cause brittle fractures and release seismic energy.

Erupting magma often contains crystals over a wide range of sizes and shapes, potentially affecting magma viscosity over many orders of magnitude. A robust relation between viscosity and the modality of crystal sizes and shapes remains lacking, principally because of the dimensional complexity and size of the governing parameter space. We have performed a suite of shear viscosity measurements on liquid-particle suspensions of dynamical similarity to crystal-bearing magma. Our experiments encompass five suspension types, each consisting of unique mixtures of two different particle sizes and shapes. The experiments span two orthogonal subspaces of particle concentration, as well as particle size and shape for each suspension type, thereby providing insight into the topology of parameter space. For each suspension type, we determined the dry maximum packing fraction and measured shear rates across a range of applied shear stresses. The results were fitted using a Herschel-Bulkley model and augment existing predictive capabilities. We demonstrate that our results are consistent with previous work, including friction-based constitutive laws for granular materials. We conclude that predictions for ascent rates of crystal-rich magmas must take the shear-rate dependence of viscosity into account. Shear-rate dependence depends first and foremost on the volume fraction of crystals, relative to the maximum packing fraction, which in turn depends on crystal size and shape distribution.

Pyroclasts from explosive eruptions, such as the 1060 CE explosive Glass Mountain eruption of Medicine Lake volcano, California, contain large amounts of water. This may be the consequence of diffusive rehydration of the volcanic glass by meteoric (secondary) water after the eruption. Discriminating between magmatic and secondary water in the matrix glass of pyroclasts is important, because the degassing of magmatic water affects the intensity of volcanic eruptions. Such discrimination has remained a challenging problem, especially because some aspects of water diffusion in silicate glasses at low temperatures and atmospheric pressure remain poorly constrained. We used thermogravimetry to analyze the loss of water from natural volcanic glasses and glasses that were hydrated in the laboratory at magmatic temperatures and pressures. Numerical modeling of diffusive water loss during thermogravimetric analyses accounted for the interconversion of molecular water (H 2 O m ) and hydroxyls groups (OH), and indicates that Glass Mountain pumices contain 0.2–0.5 wt% primary water, but gained 1–2 wt% of meteoric water by diffusive rehydration during the past 950 years. These results confirm that the majority of magmatic water is lost from the magma during explosive eruptions. Furthermore, the integration of thermogravimetric analysis and numerical modeling facilitates discrimination between the magmatic and secondary water content of volcanic glasses.

Many studies report that, under some circumstances, amending soil with biochar can improve field capacity and plant-available water. However, little is known about the mechanisms that control these improvements, making it challenging to predict when biochar will improve soil water properties. To develop a conceptual model explaining biochar's effects on soil hydrologic processes, we conducted a series of well constrained laboratory experiments using a sand matrix to test the effects of biochar particle size and porosity on soil water retention curves. We showed that biochar particle size affects soil water storage through changing pore space between particles (interpores) and by adding pores that are part of the biochar (intrapores). We used these experimental results to better understand how biochar intrapores and biochar particle shape control the observed changes in water retention when capillary pressure is the main component of soil water potential. We propose that biochar's intrapores increase water content of biochar-sand mixtures when soils are drier. When biochar-sand mixtures are wetter, biochar particles' elongated shape disrupts the packing of grains in the sandy matrix, increasing the volume between grains (interpores) available for water storage. These results imply that biochars with a high intraporosity and irregular shapes will most effectively increase water storage in coarse soils.

Silicic volcanic eruptions commonly show abrupt shifts between powerful and dangerous (Plinian) explosive episodes and gentle effusion of lava. Whether the onset of magma permeability and ensuing gas loss controls these transitions has been a subject of debate. We measured porosities and permeabilities in samples from the A.D. 1912 eruption of Novarupta volcano, Alaska, and analyzed them within the context of a well-constrained eruptive sequence that encompasses sustained explosive and effusive activity. For the explosive samples, we find that the degree of vesicle interconnectivity, measured as the ratio of connected to total porosity, decreases with phenocryst content and with increasing eruption intensity. Permeabilities of explosive samples show a weak dependence on porosity. Dome samples are not significantly different in permeability, but are of lower porosity, which together with abundant flattened vesicles is consistent with bubble collapse by permeable outgassing. Quantitative analysis indicates that outgassing alone was insufficient to affect the transition to effusive activity. Rather, the change from explosive to effusive activity was probably a consequence of high versus low magma ascent rates.

It is clear that the density and porosity of biochar will impact its mobility in the environment, its interaction with the soil hydrologic cycle, and its suitability as an ecological niche for soil microorganisms. However, the wide range of biochar pore sizes complicates biochar porosity characterization, making it challenging to find methods appropriate to connect the fundamental physical properties of density and porosity to environmental outcomes. Here, we report the use of two fast, simple density measurement techniques to characterize biochar density and porosity. We measured biochar skeletal density by helium pycnometry and envelope density by displacement of a dry granular suspension. We found that biochar skeletal density ranged from 1.34 g cm−3 to 1.96 g cm−3, and increased with pyrolysis temperature. Biochar envelope density ranged from 0.25 g cm−3 to 0.60 g cm−3, and was higher for wood biochars than grass biochars—a difference we attribute to plant cell structures preserved during pyrolysis. We compared the pore volumes measured by pycnometry with those measured by nitrogen gas sorption and mercury porosimetry. We show that biochar pore volumes measured by pycnometry are comparable to the values obtained by mercury porosimetry, the current benchmark method. We also show that the majority of biochar pore volume is in macropores, and thus, is not measured by gas sorption analysis. These fast, simple techniques can now be used to study the relationship between biochar's physical properties and its environmental behaviors.

We propose a new bubble dynamics model to study the evolution of a suspension of bubbles over a wide range of vesicularity, and that accounts for hydrodynamical interactions between bubbles while they grow, deform under shear flow conditions, and exchange mass by diffusion coarsening. The model is based on a lattice Boltzmann method for free surface flows. As such, it assumes an infinite viscosity contrast between the exsolved volatiles and the melt. Our model allows for coalescence when two bubbles approach each other because of growth or deformation. The parameter (disjoining pressure) that controls the coalescence efficiency, i.e., drainage time for the fluid film between the bubbles, can be set arbitrarily in our calculations. We calibrated this parameter by matching the measured time for the drainage of the melt film across a range of Bond numbers (ratio of buoyancy to surface tension stresses) with laboratory experiments of a bubble rising to a free surface. The model is then used successfully to model Ostwald ripening and bubble deformation under simple shear flow conditions. The results we obtain for the deformation of a single bubble are in excellent agreement with previous experimental and theoretical studies. For a suspension, we observe that the collective effect of bubbles is different depending on the relative magnitude of viscous and interfacial stresses (capillary number). At low capillary number, we find that bubbles deform more readily in a suspension than for the case of a single bubble, whereas the opposite is observed at high capillary number.

[1] The transfer of volatiles from the Earth's interior to the atmosphere occurs through degassing of magma, the dynamics of which assert a significant control on volcanic eruptions. The first and most critical step in degassing is the nucleation of gas bubbles, which requires that a sufficient number of volatile molecules cluster together to overcome the free energy associated with the formation of a new interface between nucleus and surrounding melt. This free energy is a function of surface tension, typically assumed to equate to the macroscopically measurable value. Surface tension estimates inferred from bubble nucleation experiments in silicate melts are, however, lower than direct macroscopic measurements, making it difficult to accurately predict magma ascent and decompression rates from measured bubble number densities in pyroclasts. We provide a potential resolution to this problem through an integrated study of bubble nucleation experiments and modeling thereof, based on nonclassical nucleation theory. We find that surface tension between critical bubble nuclei and the surrounding melt depends on the degree of supersaturation and is lower than the macroscopically measured value. This is consistent with the view that far from equilibrium the interface between a nucleus and surrounding metastable bulk phase is diffuse instead of sharp. As a consequence, the increase in nucleation rate with supersaturation is significantly larger at high supersaturations than predicted by classical nucleation theory.

[1] We present the results of new laboratory experiments that provide constraints on inter bubble film thinning and bubble coalescence as a consequence of liquid expulsion by gravitational and capillary forces. To ensure dynamic similarity to magmatic systems, the experiments are at small Reynolds numbers (Re≪1) and cover a wide range of Bond numbers (10−3 ≤ Bo ≤ 102). Results indicate that at Bo 0.25 gravitational forces result in film thinning. The film drainage time scale is given by t ∼ C ln (α) τ and is orders of magnitude faster than often assumed for magmatic systems. Here, C ∼ 10 is an empirical constant and α is the ratio of initial film thickness to film thickness at the time of rupture and τ is the characteristic capillary or buoyancy time scale at values of Bo 0.25, respectively.

[1] Volcanic eruptions are often accompanied by spatiotemporal migration of ground deformation, a consequence of pressure changes within magma reservoirs and pathways. We modeled the propagation of pressure variations through the east rift zone (ERZ) of Klauea Volcano, Hawai‘i, caused by magma withdrawal during the early eruptive episodes (1983–1985) of the ongoing Pu‘u ‘Ō‘ō-Kupaianaha eruption. Eruptive activity at the Pu‘u ‘Ō‘ō vent was typically accompanied by abrupt deflation that lasted for several hours and was followed by a sudden onset of gradual inflation once the eruptive episode had ended. Similar patterns of deflation and inflation were recorded at Klauea's summit, approximately 15 km to the northwest, albeit with time delays of hours. These delay times can be reproduced by modeling the spatiotemporal changes in magma pressure and flow rate within an elastic-walled dike that traverses Klauea's ERZ. Key parameters that affect the behavior of the magma-dike system are the dike dimensions, the elasticity of the wall rock, the magma viscosity, and to a lesser degree the magnitude and duration of the pressure variations themselves. Combinations of these parameters define a transport efficiency and a pressure diffusivity, which vary somewhat from episode to episode, resulting in variations in delay times. The observed variations in transport efficiency are most easily explained by small, localized changes to the geometry of the magma pathway.

The matrix-glass water concentrations in samples from volcanic eruptions of intermediate to highly silicic magmas were measured and compiled. They range from 0.1 wt% to more than 3.5 wt% and show a positive correlation with vesicles surface area over glass volume ratio. Modeling of water diffusion suggests that most of this correlation can be explained by the post-eruptive diffusion of external water at atmospheric temperature and pressure into the matrix-glass, a process referred to as rehydration. Although the precise proportion of primary (magmatic) to secondary (meteoric) water is not determined by our analysis, we find that most samples can be modeled by progressive rehydration of an initially ‘dry’ sample during the time interval between deposition and sample collection at an average rehydration diffusivity of approximately 10−23 m2 s−1. This diffusivity estimate is consistent with values provided in the literature on obsidian hydration dating and with the extrapolation of diffusivity formulations for silicic melts to lower temperatures and pressures.

Microtextures of juvenile pyroclasts from Kīlauea’s (Hawai‘i) early A.D. 2008 explosive activity record the velocity and depth of convection within the basaltic magma-filled conduit. We use X-ray microtomography (μXRT) to document the spatial distribution of bubbles. We find small bubbles (radii from 5 μm to 70 μm) in a halo surrounding larger millimeter-size bubbles. This suggests that dissolved water was enriched around the larger bubbles—the opposite of what is expected if bubbles grow as water diffuses into the bubble. Such volatile enrichment implies that the volatiles within the large bubbles were redissolving into the melt as they descended into the conduit by the downward motion of convecting magma within the lava lake. The thickness of the small bubble halo is ∼100–150 μm, consistent with water diffusing into the melt on time scales on the order of 103 s. Eruptions, triggered by rockfall, rapidly exposed this magma to lower pressures, and the haloes of melt with re-dissolved water became sufficiently supersaturated to cause nucleation of the population of smaller bubbles. The required supersaturation pressures are consistent with a depth of a few hundred meters and convection velocities of the order of 0.1 m s−1, similar to the circulation velocity observed on the surface of the Halema‘uma‘u lava lake.

The decompression rate of magma as it ascends during volcanic eruptions is an important but poorly constrained parameter that controls many of the processes that influence eruptive behavior. In this study, we quantify decompression rates for basaltic magmas using volatile diffusion in olivine-hosted melt tubes (embayments) for three contrasting eruptions of Kīlauea volcano, Hawaii. Incomplete exsolution of H2O, CO2, and S from the embayment melts during eruptive ascent creates diffusion profiles that can be measured using microanalytical techniques, and then modeled to infer the average decompression rate. We obtain average rates of ~0.05–0.45 MPa s−1 for eruptions ranging from Hawaiian style fountains to basaltic subplinian, with the more intense eruptions having higher rates. The ascent timescales for these magmas vary from around ~5 to ~36 min from depths of ~2 to ~4 km, respectively. Decompression-exsolution models based on the embayment data also allow for an estimate of the mass fraction of pre-existing exsolved volatiles within the magma body. In the eruptions studied, this varies from 0.1 to 3.2 wt% but does not appear to be the key control on eruptive intensity. Our results do not support a direct link between the concentration of pre-eruptive volatiles and eruptive intensity; rather, they suggest that for these eruptions, decompression rates are proportional to independent estimates of mass discharge rate. Although the intensity of eruptions is defined by the discharge rate, based on the currently available dataset of embayment analyses, it does not appear to scale linearly with average decompression rate. This study demonstrates the utility of the embayment method for providing quantitative constraints on magma ascent during explosive basaltic eruptions.

Over the past decade, Kīlauea and Mauna Loa—adjacent volcanoes in Hawai‘i—have inflated and deflated in tandem, yet their shallow volcanic plumbing systems are separate. Numerical modelling of the volcanoes shows that dynamic stress transfer by asthenospheric pore pressure is a viable mechanism for volcano coupling in Hawai‘i.

Hawaiian eruptions are characterized by fountains of gas and ejecta, sustained for hours to days that reach tens to hundreds of meters in height. Quantitative analysis of the pyroclastic products from the 1959 eruption of Kīlauea Iki, Kīlauea volcano, Hawai‘i, provides insights into the processes occurring during typical Hawaiian fountaining activity. This short-lived but powerful eruption contained 17 fountaining episodes and produced a cone and tephra blanket as well as a lava lake that interacted with the vent and fountain during all but the first episode of the eruption, the focus of this paper. Microtextural analysis of Hawaiian fountaining products from this opening episode is used to infer vesiculation processes within the fountain and shallow conduit. Vesicle number densities for all clasts are high (106–107 cm−3). Post-fragmentation expansion of bubbles within the thermally-insulated fountain overprints the pre-fragmentation bubble populations, leading to a reduction in vesicle number density and increase in mean vesicle size. However, early quenched rims of some clasts, with vesicle number densities approaching 107 cm−3, are probably a valid approximation to magma conditions near fragmentation. The extent of clast evolution from low vesicle-to-melt ratio and corresponding high vesicle number density to higher vesicle-to-melt ratio and lower vesicle-number density corresponds to the length of residence time within the fountain.

Biochar may improve soil hydrology by altering soil porosity, density, hydraulic conductivity, and water-holding capacity. These properties are associated with the grain size distributions of both soil and biochar, and therefore may change as biochar weathers. Here we report how freeze-thaw (F-T) cycling impacts the grain size of pine, mesquite, miscanthus, and sewage waste biochars under two drainage conditions: undrained (all biochars) and a gravity-drained experiment (mesquite biochar only). In the undrained experiment plant biochars showed a decrease in median grain size and a change in grain-size distribution consistent with the flaking off of thin layers from the biochar surface. Biochar grain size distribution changed from unimodal to bimodal, with lower peaks and wider distributions. For plant biochars the median grain size decreased by up to 45.8% and the grain aspect ratio increased by up to 22.4% after 20 F-T cycles. F-T cycling did not change the grain size or aspect ratio of sewage waste biochar. We also observed changes in the skeletal density of biochars (maximum increase of 1.3%), envelope density (maximum decrease of 12.2%), and intraporosity (porosity inside particles, maximum increase of 3.2%). In the drained experiment, mesquite biochar exhibited a decrease of median grain size (up to 4.2%) and no change of aspect ratio after 10 F-T cycles. We also document a positive relationship between grain size decrease and initial water content, suggesting that, biochar properties that increase water content, like high intraporosity and pore connectivity large intrapores, and hydrophilicity, combined with undrained conditions and frequent F-T cycles may increase biochar breakdown. The observed changes in biochar particle size and shape can be expected to alter hydrologic properties, and thus may impact both plant growth and the hydrologic cycle.

High 3He/4He ratios sampled at many ocean islands are traditionally attributed to an essentially undegassed lower mantle reservoir with high primordial 3He concentrations. A large and mostly undegassed mantle reservoir is also thought to be required to balance the Earth's 40Ar budget. On the other hand, geophysical and geochemical observations suggest slab subduction into the lower mantle, implying that most of the Earth's mantle should have been processed by partial melting over the history of the Earth, leaving noble gases in both upper and lower mantle extensively outgassed. Helge Gonnermann and Sujoy Mukhopadhyay propose a simple solution to this contradiction, which lies with the recycling and mixing of noble-gas-depleted slabs. Such mixing dilutes the concentrations of noble gases in the mantle, thereby decreasing the rate of mantle degassing and leaving significant amounts of noble gases in the processed mantle. As a result, over the 4.5 billion years of Earth's history, the lower mantle can preserve high 3He contents, high 3He/4He ratios and 40Ar concentrations high enough to satisfy the 40Ar mass balance of the Earth, even with a substantial mass flux from the upper to lower mantle. An essentially undegassed mantle reservoir is thought to be required to balance the Earth's 40Ar budget and an undegassed lower mantle reservoir is traditionally believed to be responsible for the high 3He/4He ratios sampled at many ocean islands. On the other hand, geophysical and geochemical observations suggest slab subduction into the lower mantle; the subsequent process of partial melting over the history of the Earth would leave noble gases in both upper and lower mantle extensively outgassed. A simple solution to this contradiction is now proposed, which lies in the recycling and mixing of noble-gas-depleted slabs. High 3He/4He ratios sampled at many ocean islands are usually attributed to an essentially undegassed lower-mantle reservoir with high 3He concentrations1,2,3,4. A large and mostly undegassed mantle reservoir is also required to balance the Earth’s 40Ar budget, because only half of the 40Ar produced from the radioactive decay of 40K is accounted for by the atmosphere and upper mantle5. However, geophysical6,7 and geochemical observations8 suggest slab subduction into the lower mantle, implying that most or all of Earth’s mantle should have been processed by partial melting beneath mid-ocean ridges and hotspot volcanoes. This should have left noble gases in both the upper and the lower mantle extensively outgassed, contrary to expectations from 3He/4He ratios and the Earth’s 40Ar budget. Here we suggest a simple solution: recycling and mixing of noble-gas-depleted slabs dilutes the concentrations of noble gases in the mantle, thereby decreasing the rate of mantle degassing and leaving significant amounts of noble gases in the processed mantle. As a result, even when the mass flux across the 660-km seismic discontinuity is equivalent to approximately one lower-mantle mass over the Earth’s history, high 3He contents, high 3He/4He ratios and 40Ar concentrations high enough to satisfy the 40Ar mass balance of the Earth can be preserved in the lower mantle. The differences in 3He/4He ratios between mid-ocean-ridge basalts and ocean island basalts, as well as high concentrations of 3He and 40Ar in the mantle source of ocean island basalts4, can be explained within the framework of different processing rates for the upper and the lower mantle. Hence, to preserve primitive noble gas signatures, we find no need for hidden reservoirs or convective isolation of the lower mantle for any length of time.

Basaltic pyroclastic volcanism takes place over a range of scales and styles, from weak discrete Strombolian explosions (∼10 2 –10 3 kg s −1 ) to Plinian eruptions of moderate intensity (10 7 –10 8 kg s −1 ). Recent well-documented historical eruptions from Etna, Kīlauea and Stromboli typify this diversity. Etna is Europe's largest and most voluminously productive volcano with an extraordinary level and diversity of Strombolian to subplinian activity since 1990. Kīlauea, the reference volcano for Hawaiian fountaining, has four recent eruptions with high fountaining (>400 m) activity in 1959, 1960, 1969 (–1974) and 1983–1986 (–2008); other summit (1971, 1974, 1982) and flank eruptions have been characterized by low fountaining activity. Stromboli is the type location for mildly explosive Strombolian eruptions, and from 1999 to 2008 these persisted at a rate of ca. 9 per hour, briefly interrupted in 2003 and 2007 by vigorous paroxysmal eruptions. Several properties of basaltic pyroclastic deposits described here, such as bed geometry, grain size, clast morphology and vesicularity, and crystal content are keys to understand the dynamics of the parent eruptions. The lack of clear correlations between eruption rate and style, as well as observed rapid fluctuations in eruptive behavior, point to the likelihood of eruption style being moderated by differences in the fluid dynamics of magma and gas ascent and the mechanism by which the erupting magma fragments. In all cases, the erupting magma consists of a mixture of melt and gaseous bubbles. The depth and rate of degassing, melt rheology, bubble rise and coalescence rates, and extent of syn-eruptive microlite growth define complex feedbacks that permit reversible shifts between fragmentation mechanisms and in eruption style and intensity. However, many basaltic explosive eruptions end after an irreversible shift to open-system outgassing and microlite crystallization in melt within the conduit. Clearer understanding of the factors promoting this diversity of basaltic pyroclastic eruptions is of fundamental importance in order to improve understanding of the range of behaviors of these volcanoes and assess hazards of future explosive events at basaltic volcanoes. The three volcanoes used for this review are the sites of large and growing volcano-tourism operations and there is a public need both for better knowledge of the volcanoes’ behavior and improved forecasting of the likely course of future eruptions.

Radioactive decay of uranium and thorium produces 4He, whereas 3He in the Earth's mantle is not produced by radioactive decay and was only incorporated during accretion-that is, it is primordial. 3He/4He ratios in many ocean-island basalts (OIBs) that erupt at hotspot volcanoes, such as Hawaii and Iceland, can be up to sixfold higher than in mid-ocean ridge basalts (MORBs). This is inferred to be the result of outgassing by melt production at mid-ocean ridges in conjunction with radiogenic ingrowth of 4He, which has led to a volatile-depleted upper mantle (MORB source) with low 3He concentrations and low 3He/4He ratios. Consequently, high 3He/4He ratios in OIBs are conventionally viewed as evidence for an undegassed, primitive mantle source, which is sampled by hot, buoyantly upwelling deep-mantle plumes. However, this conventional model provides no viable explanation of why helium concentrations and elemental ratios of He/Ne and He/Ar in OIBs are an order of magnitude lower than in MORBs. This has been described as the 'helium concentration paradox' and has contributed to a long-standing controversy about the structure and dynamics of the Earth's mantle. Here we show that the helium concentration paradox, as well as the full range of noble-gas concentrations observed in MORB and OIB glasses, can self-consistently be explained by disequilibrium open-system degassing of the erupting magma. We show that a higher CO2 content in OIBs than in MORBs leads to more extensive degassing of helium in OIB magmas and that noble gases in OIB lavas can be derived from a largely undegassed primitive mantle source.

The style and evolution of volcanic eruptions are dictated by the fluid mechanics governing magma ascent. Decompression during ascent causes dissolved volatile species, such as water and carbon dioxide, to exsolve from the melt to form bubbles, thus providing a driving force for the eruption. Ascent is influenced not only by the nucleation and growth of gas bubbles, but also magma rheology and brittle deformation (fragmentation). In fact, all processes and magma properties within the conduit interact and are coupled. Ultimately, it is the ability of gas trapped within growing bubbles to expand or to be lost by permeable gas flow, which determines whether ascending magmas can erupt nonexplosively. We review and integrate models of the primary conduit processes to show when each process or property dominates and how these interact within a conduit. In particular, we illustrate how and why ascent rate may control eruptive behavior: slowly ascending magmas erupt effusively and rapidly ascending magmas erupt explosively.

Dikes are vertical sheet-like bodies of igneous rock that were emplaced as magma within preexisting rock and represent the principal mode of magma transport within and through planetary lithospheres. They represent the critical link between regions of melt production and regions of melt accumulation and/or volcanic eruptions. Theoretical treatments of dike propagation include the fracture mechanics at the dike tip, magma flow within the dike, deformation of the surrounding wall rock, as well as cooling of the flowing magma within the dike. Two end-member cases are a finite volume of buoyant magma forming a self-propagating dike, and dikes that are connected to and fed by a magma reservoir.

We present results from a numerical conduit model of nonequilibrium magma degassing. We show that CO2/H2O concentration ratios in pyroclastic obsidian from the ca. 1340 A.D. Mono Craters eruption may record nonequilibrium degassing during magma ascent. Our results also indicate that permeability-controlled, open-system gas loss is consistent with obsidian formation and promotes nonequilibrium degassing at shallow depths. Because of the low diffusivity of CO2 relative to H2O, we find that CO2 concentrations can remain above equilibrium during magma ascent. To reproduce Mono Craters volatile concentrations, our model does not require volatile supersaturation, or an exsolved gas phase, prior to magma eruption. We estimate magmatic volatile contents of approximately 400 ppm CO2 and 4.6 wt.% H2O. However, nonequilibrium is contingent upon low rates of bubble nucleation to shallow depths, so that average bubble number densities are of order 10 11 m -3 during most of the magma ascent. If bubble number densities during ascent of the Mono Craters magma were greater than 1011 m � 3

We perform a quantitative textural analysis of banding in obsidian from Big Glass Mountain (BGM), Medicine Lake Volcano, California and from Mayor Island (MI), New Zealand. The samples are compositionally homogeneous, with banding defined by variable microlite content (BGM) or vesicularity (MI). The spatial distribution of alternating light and dark bands has a 1/wavenumber power-spectral density, as well as multifractal characteristics. Through the use of simple mathematical models, we demonstrate that formation of the observed banding structure is consistent with a continuous deformational reworking of magma. Our results indicate that changes in crystallinity (vesicularity) were concurrent with magma deformation and are random multiples of the total amount of crystallinity (vesicularity) that was already present at any given time. Banded obsidian from BGM and MI can be found in close spatial proximity to brittle deformational textures. We therefore propose that a key element of banding formation in our samples may be magma brecciation or fragmentation. Fragmentation can result in variable rates of degassing of fragments, which in turn can lead to nonuniform crystal (vesicle) nucleation and growth rates. Repeated fragmentation, welding and subsequent viscous deformation can thus lead to banding formation.

We report results from analog laboratory experiments, in which a large-scale flow is imposed upon natural convection from a hot boundary layer at the base of a large tank of corn syrup. The experiments show that the subdivision of the convective flow into four regions provides a reasonable conceptual framework for interpreting the effects of large-scale flow on plumes. Region I includes the area of the hot thermal boundary layer (TBL) that is thinned by the large-scale flow, thereby suppressing plumes. Region II encompasses the critically unstable boundary layer where plumes form. Region III is the area above the boundary layer that is devoid of plumes. Region IV comprises the area of hot upwelling and plume conduits. Quantitative analysis of our experiments results in a scaling law for heat flux from the hot boundary and for the spatial extent of plume suppression. When applied to the Earth’s core–mantle boundary (CMB), our results suggest that large-scale mantle flow, due to sinking lithospheric plates, can locally thin the TBL and suppress plume formation over large fractions of the CMB. Approximately 30% of heat flow from the core may be due to increased heat flux from plate-scale flow. Furthermore, CMB heat flux is non-uniformly distributed along the CMB, with large areas where heat flux is increased on average by a factor of 2. As a consequence, the convective flow pattern in the outer core may be affected by CMB heat-flux heterogeneity and sensitive to changes in plate-scale mantle flow. BecauseofplumesuppressionanddfocusingTofhotmantlefromtheCMBintozonesofupwellingflow,plumeconduits(hotspots)

The fragmentation of magma, containing abundant gas bubbles, is thought to be the defining characteristic of explosive eruptions. When viscous stresses associated with the growth of bubbles and the flow of the ascending magma exceed the strength of the melt, the magma breaks into disconnected fragments suspended within an expanding gas phase. Although repeated effusive and explosive eruptions for individual volcanoes are common, the dynamics governing the transition between explosive and effusive eruptions remain unclear. Magmas for both types of eruptions originate from sources with similar volatile content, yet effusive lavas erupt considerably more degassed than their explosive counterparts. One mechanism for degassing during magma ascent, consistent with observations, is the generation of intermittent permeable fracture networks generated by non-explosive fragmentation near the conduit walls. Here we show that such fragmentation can occur by viscous shear in both effusive and explosive eruptions. Moreover, we suggest that such fragmentation may be important for magma degassing and the inhibition of explosive behaviour. This implies that, contrary to conventional views, explosive volcanism is not an inevitable consequence of magma fragmentation.

The coexistence of stationary mantle plumes with plate-scale flow is problematic in geodynamics. We present results from laboratory experiments aimed at understanding the effects of an imposed large-scale circulation on thermal convection at high Rayleigh number (106≤Ra≤109) in a fluid with a temperature-dependent viscosity. In a large tank, a layer of corn syrup is heated from below while being stirred by large-scale flow due to the opposing motions of a pair of conveyor belts immersed in the syrup at the top of the tank. Three regimes are observed, depending on the ratio V of the imposed horizontal flow velocity to the rise velocity of plumes ascending from the hot boundary, and on the ratio λ of the viscosity of the interior fluid to the viscosity of the hottest fluid in contact with the bottom boundary. When V≪1 and λ≥1, large-scale circulation has a negligible effect on convection and the heat flux is due to the formation and rise of randomly spaced plumes. When V>10 and λ>100, plume formation is suppressed entirely, and the heat flux is carried by a sheet-like upwelling located in the center of the tank. At intermediate V, and depending on λ, established plume conduits are advected along the bottom boundary and ascending plumes are focused towards the central upwelling. Heat transfer across the layer occurs through a combination of ascending plumes and large-scale flow. Scaling analyses show that the bottom boundary layer thickness and, in turn, the basal heat flux q depend on the Peclet number, Pe, and λ. When λ>10, q∝Pe1/2 and when λ→1, q∝(Peλ)1/3, consistent with classical scalings. When applied to the Earth, our results suggest that plate-driven mantle flow focuses ascending plumes towards upwellings in the central Pacific and Africa as well as into mid-ocean ridges. Furthermore, plumes may be captured by strong upwelling flow beneath fast-spreading ridges. This behavior may explain why hotspots are more abundant near slow-spreading ridges than fast-spreading ridges and may also explain some observed variations of mid-ocean ridge basalt (MORB) geochemistry with spreading rate. Moreover, our results suggest that a potentially significant fraction of the core heat flux is due to plumes that are drawn into upwelling flows beneath ridges and not observed as hotspots.