skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: Flow‐ and Fracture‐Driven Bubble Throat Growth Rates and Dynamic Permeability in Crystallizing Magma
Abstract Pyroclasts typically exhibit coalesced vesicle textures, which are the evidence of bubble coalescence and the incomplete bubble wall retraction in magma during volcanic eruptions. The sizes of bubble throats or inter‐bubble apertures in permeable networks control the extent of magma outgassing, and therefore, quantifying the growth rates of the bubble throats is important but has remained poorly constrained. Using dynamically similar experiments with spontaneous bursting of a single bubble in rheologically well‐characterized particulate suspensions, we investigate the growth rate of bubble throats for a range of particle volume fractions. For suspensions with 0.50 particle volume fraction, a circular hole (bubble throat) forms following bubble bursting, which after an initial fast growth starts plateauing at a throat‐bubble size ratio of 0.20. The throat growth time scale overall increases with increasing particle volume fraction due to the increase in suspension viscosity. On the other hand, bubbles in suspensions with particle volume fraction near the maximum packing fraction (0.64) exhibit a fracture‐like opening. Thus, our experimental results suggest that the plateauing of the bubble throat growth in crystal‐poor to crystal‐rich magma likely contributes to the wide occurrence of the incompletely retracted vesicle walls in pyroclasts. The implications of the flow‐ to fracture‐like growth of bubble throats on the development of dynamic permeability in magma are discussed.  more » « less
Award ID(s):
2202666
PAR ID:
10569216
Author(s) / Creator(s):
 ;  ;  ;  
Publisher / Repository:
DOI PREFIX: 10.1029
Date Published:
Journal Name:
Geochemistry, Geophysics, Geosystems
Volume:
26
Issue:
2
ISSN:
1525-2027
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Bubble and crystal textures evolve during magma ascent, altering properties that control ascent such as permeability and viscosity. Eruption style results from feedbacks between ascent, bubble nucleation and growth, microlite crystallization, and gas loss, all processes recorded in pyroclasts. We show that pyroclasts of the mafic Curacautín ignimbrite of Llaima volcano, Chile, record a history of repeated autobrecciation, fusing, and crystallization. We identified pyroclasts with domains of heterogeneous vesicle textures in sharp contact with one another that are overprinted by extensive microlite crystallization. Broken crystals with long axes (l) >10 μm record fragmentation events during the eruption. A second population of unbroken microlites with l ≤10 μm overprint sutures between fused domains, suggesting the highly crystalline groundmass formed at shallow depths after autobrecciation and fusing. Nearly all pyroclasts contain plutonic and ancestral Llaima lithics as inclusions, implying that fusing occurs from a few kilometers depth to as shallow as the surface. We propose that Curacautín ignimbrite magma autobrecciated during ascent and proto-pyroclasts remained melt rich enough to fuse together. Lithics from the conduit margins were entrained into the proto-pyroclasts before fusing. Autobrecciation broke existing phenocrysts and microlites; rapid post-fusing crystallization then generated the highly crystalline groundmass. This proposed conduit process has implications for interpreting the products of mafic explosive eruptions. 
    more » « less
  2. 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. 
    more » « less
  3. Laboratory experiments were conducted to study particle migration and flow properties of non-Brownian, noncolloidal suspensions ranging from 10% to 40% particle volume fraction in a pressure-driven flow over and through a porous structure at a low Reynolds number. Particle concentration maps, velocity maps, and corresponding profiles were acquired using a magnetic resonance imaging technique. The model porous medium consists of square arrays of circular rods oriented across the flow in a rectangular microchannel. It was observed that the square arrays of the circular rods modify the velocity profiles and result in heterogeneous concentration fields for various suspensions. As the bulk particle volume fraction of the suspension increases, particles tend to concentrate in the free channel relative to the porous medium while the centerline velocity profile along the lateral direction becomes increasingly blunted. Within the porous structure, concentrated suspensions exhibit smaller periodic axial velocity variations due to the geometry compared to semidilute suspensions (bulk volume fraction ranges from 10% to 20%) and show periodic concentration variations, where the average particle concentration is slightly greater between the rods than on top of the rods. For concentrated systems, high particle concentration pathways aligned with the flow direction are observed in regions that correspond to gaps between rods within the porous medium. 
    more » « less
  4. Jetting dynamics from bursting bubbles play a key role in mediating mass and momentum transport across the air–liquid interface, and have attracted widespread interest from researchers across disciplines. In marine environments, this phenomenon has drawn considerable attention due to its role in releasing biochemical contaminants, such as extracellular polymeric substances, into the atmosphere through aerosol production. These biocontaminants often exhibit non-Newtonian characteristics, yet the physics of bubble bursting with a rheologically complex layer at the bubble–liquid interface remains largely unexplored. In this study, we experimentally investigate the jetting dynamics of bubble bursting events in the presence of such a polymeric compound layer. Using bubbles coated by a polyethylene oxide solution, we document the cavity collapse and jetting dynamics produced by bubble bursting. At a fixed polymer concentration, the jet velocity increases while the jet radius decreases with an increasing compound layer volume fraction, as a result of stronger capillary wave damping due to capillary wave separation at the compound interface as well as the formation of smaller cavity cone angles during bubble cavity collapse. These dynamics produce smaller and more numerous jet drops. Meanwhile, as the polymer concentration increases, the jet velocity decreases while the jet radius increases for the same compound layer fraction due to the increasing viscoelastic stresses. In addition, fewer jet drops are ejected as the jets become slower and broader with increasing polymer concentration, as viscoelastic stresses persist throughout the jet formation and thinning process. We further obtain, for the first time, a regime map delineating the conditions for jet drop ejection versus no jet drop ejection in bursting bubbles coated with a polymeric compound layer. Our results may provide new insights into the mechanisms of mass transport of organic materials in bubble-mediated aerosolization processes, advancing our understanding of marine biology and environmental science. 
    more » « less
  5. We document the presence, composition, and number density (TND) of titanomagnetite nanolites and ultra‐nanolites in aphyric rhyolitic pumice, obsidian, and vesicular obsidian from the 1060 CE Glass Mountain volcanic eruption of Medicine Lake Volcano, California, using magnetic methods. Curie temperatures indicate compositions of Fe2.40Ti0.60O4 to Fe3O4. Rock‐magnetic parameters sensitive to domain state, which is dependent on grain volume, indicate a range of particle sizes spanning superparamagnetic (<50–80 nm) to multidomain (>10 μm) particles. Cylindrical cores drilled from the centers of individual pumice clasts display anisotropy of magnetic susceptibility with prolate fabrics, with the highest degree of anisotropy coinciding with the highest vesicularity. Fabrics within a pumice clast require particle alignment within a fluid, and are interpreted to result from the upward transport of magma driven by vesiculation, ensuing bubble growth, and shearing in the conduit. Titanomagnetite number density (TND) is calculated from titanomagnetite volume fraction, which is determined from ferromagnetic susceptibility. TND estimates for monospecific assemblages of 1,000 nm–10 nm cubes predict 10^12 to 10^20 m^−3 of solid material, respectively. TND estimates derived using a power law distribution of grain sizes predict 10^18 to 10^19  m^−3. These ranges agree well with TND determinations of 10^18 to 10^20  m^−3 made by McCartney et al. (2024), and are several orders of magnitude larger than the number density of bubbles in these materials. These observations are consistent with the hypothesis that titanomagnetite crystals already existed in extremely high number‐abundance at the time of magma ascent and bubble nucleation. 
    more » « less