skip to main content


Title: Autobrecciation and fusing of mafic magma preceding explosive eruptions
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
Award ID(s):
1831213
NSF-PAR ID:
10345564
Author(s) / Creator(s):
; ; ;
Date Published:
Journal Name:
Geology
ISSN:
0091-7613
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. 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
  2. Abstract

    Plagioclase microlites in a magma nucleate and grow in response to melt supersaturation (Δϕplag). The resultant frozen plagioclase crystal size distribution (CSD) preserves the history of decompression pathways (dP/dt). SNGPlag is a numerical model that calculates the equilibrium composition of a decompressing magma and nucleates and grows plagioclase in response to an imposed Δϕplag. Here, we test a new version of SNGPlag calibrated for use with basaltic andesite magmas and modeldP/dtfor the ca. 12.6 ka Curacautín eruption of Llaima volcano, Chile. Instantaneous nucleation (Nplag) and growth (Gplag) rates of plagioclase were computed using the experimental results of Shea and Hammer (J Volcanol Geotherm Res 260:127–145, 10.1016/j.jvolgeores.2013.04.018, 2013) and used for SNGPlag modeling of basaltic andesite composition. MaximumNplagof 6.1 × 105 cm h−1is achieved at a Δϕplagof 44% and the maximumGplagof 27.4 μm h−1is achieved at a Δϕplagof 29%. Our modeled logdP/dtavgrange from 2.69 ± 0.09 to 6.89 ± 0.96 MPa h−1(1σ) with an average duration of decompression from 0.87 ± 0.25 to 16.13 ± 0.29 h assuming a starting pressurePiof 110–150 MPa. These rates are similar to those derived from mafic decompression experiments for other explosive eruptions. Using assumptions for lithostatic pressure gradients (dP/dz), we calculate ascent rates of < 1–6 m s−1. We conducted a second set of Monte Carlo simulations usingPiof 15–30 MPa to investigate the influence of shallower decompression, resulting in logdP/dtavgfrom 2.86 ± 0.49 to 6.00 ± 0.86 MPa h−1. ThedP/dtmodeled here is two orders of magnitude lower than those calculated by Valdivia et al. (Bull Volcanol, 10.1007/s00445-021-01514-8, 2022) for the same eruption using a bubble number density meter, and suggests homogeneous nucleation raisesdP/dtby orders of magnitude in the shallow conduit. Our modeling further supports the rapid-ascent hypothesis for driving highly explosive mafic eruptions.

     
    more » « less
  3. null (Ed.)
    Dense, vitric, dacitic pyroclasts (dacite lithics) from the 1991 preclimactic explosions of Mt. Pinatubo were analyzed for their vesicular and crystal textures and dissolved H2O and CO2 contents. Micron-scale heterogeneities in groundmass glass volatile contents (0.9 wt% differences in H2O within 500 μm) are observed and argue that parts of the dacite lithics equilibrated at different depths before finally being constructed. Greater vesicularities and larger and greater number densities of vesicles are observed in groundmass glass around phenocrysts compared to groundmass glass away from phenocrysts, similar to textures produced in experiments that sintered bimodal distributions of particles. Furthermore, increasingly greater proportions of stretched and distorted vesicles are observed in lithics from the later explosions, which parallels the increasingly shorter reposes between explosions. Finally, micron-sized crystal fragments are ubiquitous in groundmass glass of all dacite lithics. The textures, together with the variable volatile contents, lead us to propose a model that the dacite lithics formed by rapid and repetitive sintering of ash particles derived from a variety of depths on the conduit walls above the fragmentation level. We speculate that sintering of conduit material produced impermeable layers that retarded gas flow through the conduit, causing pressure to build until the cap failed and the next explosion occurred. 
    more » « less
  4. Puyehue-Cordon Caulle (PCC) is an active volcanic complex located in the SVZ of the Andes that has had three major historic rhyodacitic eruptions with the most recent event in 2011-12. We provide petrologic and geochemical evidence that PCC is underlain by a crystal mush using recently identified basaltic mafic enclaves that highlights the involvement of distinct mafic magma components during the 2011-12 eruption. We suggest the mafic enclaves represent remnants of the crystal-rich mush that get entrained during eruption of the crystal-poor rhyodacite melt lens cap. This architecture requires the basaltic mush to produce rhyodacite through efficient fractionation. The dominant population of enclaves are equigranular, crystal-rich (45-55%), vesiculated (10-20%), and display interlocking grains between phases. Vesicles have complex shapes filling the irregular interlocking textures, while phenocrysts show stepwise normal zoning (uniform plagioclase cores, ~An90, overgrown with weakly zoned rims, ~An60). A second porphyritic population may represent mafic recharge into the system that bypasses the mush unperturbed. The porphyritic enclaves have spherical vesicles and tightly bound primitive mineral compositions (Fo80-86 vs Fo70-86 in the equigranular enclaves). Published geothermobarometry from the 2011-12 rhyodacite suggests shallow magma storage (5-7 km, 100-140 MPa, 895°C), which we compare against newly determined mineral-mineral trace-element partitioning based thermometry. Our thermometry indicates the equigranular enclaves were stored at ~900-1000°C at the time of eruption suggesting both a compositionally and thermally zoned magma system. We combine this temperature information with trace element data and mass balance calculations from various minerals phases and melt to substantiate our previous hypothesis that the basaltic enclaves can produce rhyodacite given their crystallinity. These estimates may support a spatially connected basaltic crystal-mush underlying a rhyodacite melt lens cap further proving highly efficient rhyolite formation at PCC. PCC’s enclaves present one of the largest compositional gaps on record globally. We compare them to other enclave-bearing systems and how PCC is an important end-member to understand enclaves as well as rhyolite formation. 
    more » « less
  5. Abstract The first known occurrence of rhyolite along the submarine segments of the mid-ocean ridge (MOR) system was discovered on Alarcon Rise, the northernmost segment of the East Pacific Rise (EPR), by the Monterey Bay Aquarium Research Institute in 2012. Zircon trace element and Hf and O isotope patterns indicate that the rhyolite formed by extreme crystal fractionation of primary mid-ocean ridge basalt (MORB) sourced from normal to enriched MOR mantle with little to no addition of continental lithosphere or hydrated oceanic crust. A large range in zircon ɛHf spanning 11 ɛ units is comparable to the range of whole rock ɛHf from the entire EPR. This variability is comparable to continental granitoids that develop over long periods of time from multiple sources. Zircon geochronology from Alarcon Rise suggests that at least 20 kyr was needed for rhyolite petrogenesis. Grain-scale textural discontinuities and trace element trends from zircon cores and rims are consistent with crystal fractionation from a MORB magma with possible perturbations associated with mixing or replenishment events. Comparison of whole rock and zircon oxygen isotopes with modeled fractionation and zircon-melt patterns suggests that, after they formed, rhyolite magmas entrained hydrated mafic crust from conduit walls during ascent and/or were hydrated by seawater in the vent during eruption. These data do not support a model where rhyolites formed directly from partial melts of hydrated oceanic crust or do they require assimilation of such crust during fractional crystallization, both models being commonly invoked for the formation of oceanic plagiogranites and dacites. A spatial association of highly evolved lavas (rhyolites) with an increased number of fault scarps on the northern Alarcon Rise might suggest that low magma flux for ~20 kyr facilitated extended magma residence necessary to generate rhyolite from MORB. 
    more » « less