An official website of the United States government
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
Mafic explosive volcanism at Llaima Volcano: 3D x-ray microtomography reconstruction of pyroclasts to constrain shallow conduit processes
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
- Award ID(s):
- 1831213
- PAR ID:
- 10361584
- Publisher / Repository:
- Springer Science + Business Media
- Date Published:
- Journal Name:
- Bulletin of Volcanology
- Volume:
- 84
- Issue:
- 1
- ISSN:
- 0258-8900
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Abstract After more than three months of lava dome extrusion, La Soufrière (St Vincent) transitioned to a series of explosive eruptions in April 2021. Here we present a time-series petrologic analysis of the phenocryst and microlite populations during the firstc.48 h of explosivity to constrain ascent conditions and processes that drove changes in behaviour. Primary eruptive products were crystal-rich (45–50 vol%) basaltic andesites with similar phenocryst phase assemblages and compositions. The change in eruptive style is consistent with overpressurization as a consequence of second boiling from anhydrous microlite crystallization. The microlites display variation between the explosive phases, with two populations: (1) ‘inherited’ − normally zoned high-An plagioclase (>An70) + olivine (Fo62–79) + clinopyroxene + titanomagnetite, inferred to have crystallized at depths >15 km and high water pressures; (2) ‘juvenile’ − unzoned plagioclase (An45–65) + clinopyroxene + orthopyroxene + intermediate pyroxene (Wo12–38) + titanomagnetite, inferred to have crystallized upon ascent due to decompression and degassing. Scoria from the first explosions featured extensive groundmass crystallization and a significant ‘inherited’ microlite population. Later explosions had a more abundant ‘juvenile’ microlite population and lower crystallinity, consistent with more rapid ascent from depth, initiated by decompression following initial blasts and destruction of the lava dome.more » « less
-
We investigate whether decompression rates derived from three often-disparate petrologic techniques (microlites, bubbles, and melt embayments) can be reconciled or integrated for a more complete understanding of magma ascent in the conduit. We focus on the well-studied and -documented earliest Plinian eruptions (June 12, 1991) of Mount Pinatubo. Using a newly developed two-stage decompression-diffusion model, volatile profiles in quartz-hosted embayments reveal an initial stage of decompression nearly two orders of magnitude slower than final rates. In applying time-integrated models of microlite and bubble nucleation and growth, initial decompression rates from embayments are supported by microlite modeling results, whereas final rates are in close agreement with bubble number densities. This consistency and continuity between speedometers supports the sensitivity of different petrologic recorders to specific regions of the conduit system and highlights the fidelity of embayments as recorders of decompression throughout the entire conduit. Ascent timescales derived from Pinatubo embayments range from hours to days, coinciding with the visual onset of lava effusion leading to explosive activity.more » « less
-
null (Ed.)Silicic volcanic activity has long been framed as either violently explosive or gently effusive. However, recent ob- servations demonstrate that explosive and effusive behavior can occur simultaneously. Here, we propose that rhyolitic magma feeding subaerial eruptions generally fragments during ascent through the upper crust and that effusive eruptions result from conduit blockage and sintering of the pyroclastic products of deeper cryptic frag- mentation. Our proposal is supported by (i) rhyolitic lavas are volatile depleted; (ii) textural evidence supports a pyroclastic origin for effusive products; (iii) numerical models show that small ash particles !10−5 m can diffusive- ly degas, stick, and sinter to low porosity, in the time available between fragmentation and the surface; and (iv) inferred ascent rates from both explosive and apparently effusive eruptions can overlap. Our model reconciles previously paradoxical observations and offers a new framework in which to evaluate physical, numerical, and geochemical models of Earth’s most violent volcanic eruptions.more » « less
-
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
