Bubble nucleation is the critical first step during magma degassing. The resultant number density of bubbles provides a record of nucleation kinetics and underlying eruptive conditions. The rate of bubble nucleation is strongly dependent on the surface free energy associated with nucleus formation, making the use of bubble number density for the interpretation of eruptive conditions contingent upon a sound understanding of surface tension. Based on a suite of nucleation experiments with up to >1016bubbles per unit volume of melt, and using numerical simulations of bubble nucleation and growth during each experiment, we provide a new formulation for surface tension during homogeneous nucleation of H2O bubbles in rhyolitic melt. It is based on the Tolman correction with a Tolman length of
Bubbles in magmas drive explosive volcanic eruptions. The spatial distribution of bubble nucleation sites in an ascending, decompressing, and supersaturating magma is one of the primary controls on ash morphologies and volcanic hazards. The mechanism of bubble formation is important because it ultimately determines the spatial distribution of bubbles in the fragmenting magma. The initial nucleation of bubbles in a homogeneous magma is problematical because excessive surface tension pressure in very small, nascent bubbles should drive exsolved volatiles back into the melt. This thermodynamic barrier to bubble viability confounds understanding of homogeneous bubble nucleation, yet very small bubbles form, grow, and ultimately drive explosive volcanic eruptions. We refer to this as “the tiny bubble paradox.” Classical nucleation theory typically explains bubble formation and growth, but we propose that a spectrum of bubble‐forming mechanisms may include both homogeneous nucleation and spinodal decomposition (the spontaneous unmixing of phases by uphill diffusion) as end‐member processes. As spinodal decomposition progresses, regularly sized and regularly spaced quasi‐spherical zones form with increasingly high concentration of dissolved water at the centers. Bubble formation occurs as the concentration of water in the interior of the water‐rich zones approaches 100% and the concentration of melt approaches zero. The presence of a broad, diffuse, concentration gradient of water rather than a narrow water‐melt interface means that there is no surface, per se, for surface tension to arise. This is the crux of the solution of the tiny bubble paradox.
more » « less- NSF-PAR ID:
- 10456319
- Publisher / Repository:
- DOI PREFIX: 10.1029
- Date Published:
- Journal Name:
- Geochemistry, Geophysics, Geosystems
- Volume:
- 21
- Issue:
- 6
- ISSN:
- 1525-2027
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
more » « less
Abstract δ = 0.32 nm, which implies an increase in surface tension of bubbles with decreasing nucleus size. Our model results indicate that experiments encompass two distinct nucleation regimes, distinguishable by the ratio of the characteristic diffusion time of water,τ diff, to the decompression time,t d. Experiments with >1013 m−3bubbles are characterized byτ diff/t d≪ 1, wherein the nucleation rate predominantly depends on the interplay between decompression and diffusion rates. Nucleation occurs over a short time interval with nucleation rate peaks at high values. For experiments with comparatively low bubble number density the average distance between adjacent bubbles and the diffusion timescale is large. Consequently,τ diff/t d≫ 1 and nucleation is nearly unaffected by diffusion and independent of decompression rate, with bubbles nucleating at an approximately constant rate until the sample is quenched. -
more » « less
Abstract 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.
-
more » « less
Abstract Mafic volcanic activity is dominated by effusive to mildly explosive eruptions. Plinian and ignimbrite-forming mafic eruptions, while rare, are also possible; however, the conditions that promote such explosivity are still being explored. Eruption style is determined by the ability of gas to escape as magma ascends, which tends to be easier in low-viscosity, mafic magmas. If magma permeability is sufficiently high to reduce bubble overpressure during ascent, volatiles may escape from the magma, inhibiting violent explosive activity. In contrast, if the permeability is sufficiently low to retain the gas phase within the magma during ascent, bubble overpressure may drive magma fragmentation. Rapid ascent may induce disequilibrium crystallization, increasing viscosity and affecting the bubble network with consequences for permeability, and hence, explosivity. To explore the conditions that promote strongly explosive mafic volcanism, we combine microlite textural analyses with synchrotron x-ray computed microtomography of 10 pyroclasts from the 12.6 ka mafic Curacautín Ignimbrite (Llaima Volcano, Chile). We quantify microlite crystal size distributions (CSD), microlite number densities, porosity, bubble interconnectivity, bubble number density, and geometrical properties of the porous media to investigate the role of magma degassing processes at mafic explosive eruptions. We use an analytical technique to estimate permeability and tortuosity by combing the Kozeny-Carman relationship, tortuosity factor, and pyroclast vesicle textures. The groundmass of our samples is composed of up to 44% plagioclase microlites, > 85% of which are < 10 µm in length. In addition, we identify two populations of vesicles in our samples: (1) a convoluted interconnected vesicle network produced by extensive coalescence of smaller vesicles (> 99% of pore volume), and (2) a population of very small and completely isolated vesicles (< 1% of porosity). Computed permeability ranges from 3.0 × 10−13to 6.3 × 10−12m2, which are lower than the similarly explosive mafic eruptions of Tarawera (1886; New Zealand) and Etna (112 BC; Italy). The combination of our CSDs, microlite number densities, and 3D vesicle textures evidence rapid ascent that induced high disequilibrium conditions, promoting rapid syn-eruptive crystallization of microlites within the shallow conduit. We interpret that microlite crystallization increased viscosity while simultaneously forcing bubbles to deform as they grew together, resulting in the permeable by highly tortuous network of vesicles. Using the bubble number densities for the isolated vesicles (0.1-3−3 × 104 bubbles per mm3), we obtain a minimum average decompression rate of 1.4 MPa/s. Despite the textural evidence that the Curacautín magma reached the percolation threshold, we propose that rapid ascent suppressed outgassing and increased bubble overpressures, leading to explosive fragmentation. Further, using the porosity and permeability of our samples, we estimated that a bubble overpressure > 5 MPa could have been sufficient to fragment the Curacautín magma. Other mafic explosive eruptions report similar disequilibrium conditions induced by rapid ascent rate, implying that syn-eruptive disequilibrium conditions may control the explosivity of mafic eruptions more generally.
-
Nucleation of H2O vapor bubbles in magma requires surpassing a chemical supersaturation threshold via decompression. The threshold is minimized in the presence of a nucleation substrate (heterogeneous nucleation, <50 MPa), and maximized when no nucleation substrate is present (homogeneous nucleation, >100 MPa). The existence of explosively erupted aphyric rhyolite magma staged from shallow (<100 MPa) depths represents an apparent paradox that hints at the presence of a cryptic nucleation substrate. In a pair of studies focusing on Glass Mountain eruptive units from Medicine Lake, California, we characterize titanomagnetite nanolites and ultrananolites in pumice, obsidian, and vesicular obsidian (Brachfeld et al., 2024,more » « less
https://doi.org/10.1029/2023GC011336 ), calculate titanomagnetite crystal number densities, and compare titanomagnetite abundance with the physical properties of pumice to evaluate hypotheses on the timing of titanomagnetite crystallization. Titanomagnetite crystals with grain sizes of approximately 3–33 nm are identified in pumice samples from the thermal unblocking of low‐temperature thermoremanent magnetization. The titanomagnetite number densities for pumice are 10^18 to 10^20 m^−3, comparable to number densities in pumice and obsidian obtained from room temperature methods (Brachfeld et al., 2024,https://doi.org/10.1029/2023GC011336'>https://doi.org/10.1029/2023GC011336). This range exceeds reported bubble number densities (BND) within the pumice from the same eruptive units (average BND ∼4 × 10^14 m^−3). The similar abundances of nm‐scale titanomagnetite crystals in the effusive and explosive products of the same eruption, together with the lack of correlation between pumice permeability and titanomagnetite content, are consistent with titanomagnetite formation having preceded the bubble formation. Results suggest sub‐micron titanomagnetite crystals are responsible for heterogeneous bubble nucleation in this nominally aphyric rhyolite magma. -
more » « less
Abstract The most explosive basaltic scoria cone eruption yet documented (>20 km high plumes) occurred at Sunset Crater (Arizona) ca. 1085 AD by undetermined eruptive mechanisms. We present melt inclusion analysis, including bubble contents by Raman spectroscopy, yielding high total CO2(approaching 6000 ppm) and S (~2000 ppm) with moderate H2O (~1.25 wt%). Two groups of melt inclusions are evident, classified by bubble vol%. Modeling of post-entrapment modification indicates that the group with larger bubbles formed as a result of heterogeneous entrapment of melt and exsolved CO2and provides evidence for an exsolved CO2phase at magma storage depths of ~15 km. We argue that this exsolved CO2phase played a critical role in driving this explosive eruption, possibly analogous to H2O exsolution driving silicic caldera-forming eruptions. Because of their distinct gas compositions relative to silicic magmas (high S and CO2), even modest volume explosive basaltic eruptions could impact the atmosphere.
