An official website of the United States government
Abstract Magmatic systems below volcanoes are often dominated by partially crystalline magma over the long term. Rejuvenation of these systems during eruptive events can impact lava composition and eruption style—sometimes resulting in more violent or explosive activity than would be expected, as was the case at Fissure 17 during Kīlauea’s 2018 eruption. Here, we explore how the crystallinity of unerupted intrusion magmas affect hybrid magma compositions and petrological signatures by constructing phase-equilibria models to evaluate mineral and melt compositions of low-MgO lavas erupted along the East Rift Zone of Kīlauea volcano on 30 to 31 January 1997 (Episode 54, Fissures A-F). We then compare calculated mixing proportions and petrologically derived magma volumes to GPS-based geodetic inversions of ground deformation and intrusion growth in an attempt to reconcile geodetic and petrologically estimated magma volumes. Open-system phase-equilibria thermodynamic models were used to constrain the composition, degree of differentiation, and thermodynamic state of a rift-stored, two pyroxene + plagioclase saturated low-MgO magma body immediately preceding its mixing with high-MgO recharge and degassed drainback (lava lake) magma from Pu‘u‘ō‘ō‘, shortly before fissure activity within Nāpau Crater began on 29 January 1997. Mixing models constructed using the Magma Chamber Simulator reproduce the mineralogy and compositions of Episode 54 lavas within uncertainties and suggest that the identity of the low-MgO magma body may be either variably differentiated remnants of un-erupted magmas intruded into Nāpau Crater in October 1968, or another spatially and compositionally similar magma body. We find that magmas derived from a single, compositionally stratified magma emplaced beneath Nāpau Crater in 1968 can mix with mafic Kīlauea magmas to reproduce average Episode 54 bulk lava, mineralogy and mineral compositions without necessitating the interaction of multiple, low-MgO rift-stored magma bodies to produce Episode 54 lava compositions. Further, by constructing phase equilibria-based mixing models of Episode 54, we can better define the pre-eruptive state of the magmatic system. The resultant mineral assemblages and compositions are consistent with the possibility that the now-fractionated, rift-stored magma body was compositionally stratified and ~ 40% to 50% crystalline at the time of mixing. Finally, we estimate the volume of the low-MgO magma body to be ~7.51 Mm3. Phase-equilibria model results corroborate field and geochemical relationships demonstrating how shallow intrusions at intraplate shield volcanoes can crystallize, evolve, and then be remobilized by new, later batches of mafic magma. Most notably, our MCS models demonstrate that the pre-eruptive conditions of an intrusive body may be recovered by examining mineral compositions within mixed lavas. Discrepancies between the geodetic constraints on volumes of stored rift versus newly intruded (recharge) magma and our best-fit results produced by MCS mixing models (which respectively are mmafic:mlow-MgO ≈ 2 vs. mmafic:mlow-MgO ≈ 0.75) are interpreted to highlight the complex nature of incomplete mixing on more localized scales as reflected in erupted lavas, compared to geodetically constrained volumes that likely reflect large spatial scale contributions to a magmatic system. These dissimilar volume relationships may also help to constrain eruptive versus unerupted volumes in magmatic systems undergoing mixing. By demonstrating the usefulness of MCS in modeling past eruptions, we highlight the potential to use it as a tool to aid in petrologic monitoring of ongoing activity.
more »
« less
This content will become publicly available on May 24, 2025
Magma mingling during the 1959 eruption of Kīlauea Iki, Hawaiʻi
Magma mingling and mixing are common processes at basaltic volcanoes and play a fundamental role in magma petrogenesis and eruption dynamics. Mingling occurs most commonly when hot primitive magma is introduced into cooler magma. Here, we investigate a scenario whereby cool, partially degassed lava is drained back into a conduit, where it mingles with hotter, less degassed magma. The 1959 eruption of Kīlauea Iki, Hawaiʻi involved 16 high fountaining episodes. During each episode, fountains fed a lava lake in a pit crater, which then partially drained back into the conduit during and after each episode. We infer highly crystalline tachylite inclusions and streaks in the erupted crystal-poor scoria to be the result of the recycling of this drain-back lava. The crystal phases present are dendrites of plagioclase, augite and magnetite/ilmenite, at sizes of up to 10 μm. Host sideromelane glass contains 7–8 wt% MgO and the tachylite glass (up to 0.5% by area) contains 2.5–6 wt% MgO. The vesicle population in the tachylite is depleted in the smallest size classes (< 0.5 mm) and has overall lower vesicle number densities and a higher degree of vesicle coalescence than the sideromelane component. The tachylite exhibits increasingly complex ‘stretching and folding’ mingling textures through the episodes, with discrete blocky tachylite inclusions in episodes 1 and 3 giving way to complex, folded, thin filaments of tachylite in pyroclasts erupted in episodes 15 and 16. We calculate that a lava lake crust 8–35 cm thick may have formed in the repose times between episodes, and then foundered and been entrained into the conduit during drain-back. The recycled fragments of crust would have been reheated in the conduit, inducing glass devitrification and crystallisation of pyroxene, magnetite and plagioclase dendrites and eventually undergoing ductile flow as the temperature of the fragments approached the host magma temperature. We use simple models of magma mingling to establish that stretching and folding of recycled, ductile lava could involve thinning of the clasts by up to a factor of 10 during the timescale of the eruption, consistent with observations of streaks and filaments of tachylite erupted during episodes 15 and 16, which may have undergone multiple cycles of eruption, drain-back and reheating.
more »
« less
- Award ID(s):
- 2119838
- PAR ID:
- 10540039
- Publisher / Repository:
- Springer-Verlag
- Date Published:
- Journal Name:
- Bulletin of volcanology
- ISSN:
- 1432-0819
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Abstract Ignimbrite flare-ups are rare periods of intense silicic volcanism during which the pyroclastic volume and eruptive frequency is more than an order of magnitude higher than background activity. Investigating the compositional differences between flare-up and steady-state magmas provides critical constraints on the petrogenetic causes for the event and can offer unique opportunities to investigate the role of large-scale tectonic or geodynamic processes in arc magmatism. In this study, we focus on the bimodal Deschutes Formation ignimbrite flare-up of Central Oregon, which erupted unusually high volumes of pyroclastic material 6.25–5.45 Ma from a new axis of volcanism in the Cascades arc. This episode is marked by increased eruption rates and eruption of more silicic compositions relative to the Quaternary Cascade arc, which rarely erupts rhyolites. Ignimbrites are crystal-poor (<10%) dacite to rhyolites (mostly 65–77 wt.% SiO2) with anhydrous mineral assemblages and higher FeO/MgO, Y, Eu/Eu*, MREE and Zr/Sr, indicating drier magmatic evolution compared to the Quaternary arc, and are more similar to those from the rear-arc High Lava Plains (HLP) province that lies to the east. Magnetite-ilmenite oxybarometry indicates that Deschutes Formation felsic magmas tend to be hotter and more reduced (NNO-1 to NNO) than the Quaternary arc (NNO to NNO + 1.5). Rhyolite-MELTS geobarometry suggests complex storage of diverse Deschutes Formation magmas within the shallow crust (50–250 MPa), and the common co-eruption of multiple plagioclase populations, pumice compositions, and compositionally banded pumice suggest variable degrees of mixing and mingling of distinct magmas. Deschutes magmas also have low δ18Oplagioclase values that indicate partial melting and assimilation of hydrothermally altered shallow crust. Trace element systematics and rhyolite-MELTS modeling suggests that felsic pumice cannot be produced by simple fractionation of co-erupted mafic pumice or basaltic lavas, and requires a crustal melting origin, and trace elements and Pb isotopes suggest that young mafic crust may have been the primary protolith. We suggest that partial melting produced low-Si rhyolite melt (~72 wt.%) that acted as both a parent for the most evolved rhyolites, and as a mixing endmember to create the dacite to rhyodacite magmas with heterogenous plagioclase populations. Unlike the predominantly calc-alkaline basalts erupted in the Quaternary Cascade arc, Deschutes Formation primary basalts are mostly low-K tholeiites, indicative of decompression melting. These are similar to the compositions erupted during a contemporaneous pulse of low-K tholeiite volcanism across the whole HLP that reached into the Cascades rear-arc. We suggest that intra-arc extension focused decompression melts from the back-arc into the arc and that tensional stresses allowed this high flux of hot-dry-reduced basalt throughout the crustal column, causing partial melting of mafic protoliths and the production of hot-dry-reduced rhyolite melts. Depletion of incompatible elements in successive rhyolites implies progressive depletion in fertility of the protolith. Extension also allowed for the establishment of a robust hydrothermal system, and assimilation of hydrothermally-altered rocks by magmas residing in a shallow, complex storage network lead to low δ18O melts. Our findings suggest the integral role that extensional tectonics played in producing an unusual ignimbrite flare-up of hot-dry-reduced rhyolite magmas that are atypical of the Cascades arc and may be an important contributor to flare-ups at arcs worldwide.more » « less
-
Abstract Two distinct types of rare crystal-rich mafic enclaves have been identified in the rhyolite lava flow from the 2011–12 Cordón Caulle eruption (Southern Andean Volcanic Zone, SVZ). The majority of mafic enclaves are coarsely crystalline with interlocking olivine-clinopyroxene-plagioclase textures and irregular shaped vesicles filling the crystal framework. These enclaves are interpreted as pieces of crystal-rich magma mush underlying a crystal-poor rhyolitic magma body that has fed recent silicic eruptions at Cordón Caulle. A second type of porphyritic enclaves, with restricted mineral chemistry and spherical vesicles, represents small-volume injections into the rhyolite magma. Both types of enclaves are basaltic end-members (up to 9.3 wt% MgO and 50–53 wt% SiO 2 ) in comparison to enclaves erupted globally. The Cordón Caulle enclaves also have one of the largest compositional gaps on record between the basaltic enclaves and the rhyolite host at 17 wt% SiO 2 . Interstitial melt in the coarsely-crystalline enclaves is compositionally identical to their rhyolitic host, suggesting that the crystal-poor rhyolite magma was derived directly from the underlying basaltic magma mush through efficient melt extraction. We suggest the 2011–12 rhyolitic eruption was generated from a primitive basaltic crystal-rich mush that short-circuited the typical full range of magmatic differentiation in a single step.more » « less
-
Abstract High-threat explosive silicic eruptions commonly contain banded pumice, reflecting magma mingling in the conduit prior to or during eruption. Heterogeneities in tuffs have been attributed to the draw-up of compositionally distinct magmas, in which low-viscosity magmas ascend more quickly than high-viscosity magmas. The Rattlesnake Tuff of the High Lava Plains in Oregon (northwestern United States) represents a zoned magma reservoir where at least five different rhyolite compositions are preserved in banded pumice samples in variable mingled combinations. Geochemical gradients recorded across band boundaries in pumice were modeled using a Monte Carlo least-square minimization procedure to find the complementary error function that best fit observed Si and Ba diffusion profiles by iteratively varying the concentration of each plateau (i.e., the concentration on either side of the band boundary), the center and spacing of the diffusion profile, diffusion length scale, and temperature. Modeling indicates maximum time scales between mingling and conduit ascent from minutes to hours. Viscosity calculations for each rhyolite composition confirm that highly viscous rhyolites have longer ascent times than low-viscosity magmas, strongly supporting a model of sequential tapping of a zoned chamber controlled by viscosity.more » « less
-
Abstract Magmas with matrix glass compositions ranging from basalt to dacite erupted from a series of 24 fissures in the first 2 weeks of the 2018 Lower East Rift Zone (LERZ) eruption of Kīlauea Volcano. Eruption styles ranged from low spattering and fountaining to strombolian activity. Major element trajectories in matrix glasses and melt inclusions hosted by olivine, pyroxene and plagioclase are consistent with variable amounts of fractional crystallization, with incompatible elements (e.g., Cl, F, and H2O) becoming enriched by 4–5 times as melt MgO contents evolve from 6 to 0.5 wt%. The high viscosity and high H2O contents (∼2 wt%) of the dacitic melts erupting at Fissure 17 account for the explosive Strombolian behavior exhibited by this fissure, in contrast to the low fountaining and spattering observed at fissures erupting basaltic to basaltic‐andesite melts. Saturation pressures calculated from melt inclusion CO2‐H2O contents indicate that the magma reservoir(s) supplying these fissures was located at ∼2–3 km depth, which is in agreement with the depth of a dacitic magma body intercepted during drilling in 2005 (∼2.5 km) and a seismically imaged lowVp/Vsanomaly (∼2 km depth). Nb/Y ratios in erupted products are similar to lavas erupted between 1955 and 1960, indicating that melts were stored and underwent variable amounts of crystallization in the LERZ for >60 years before being remobilized by a dike intrusion in 2018. We demonstrate that extensive fractional crystallization generates viscous and volatile‐rich magma with potential for hazardous explosive eruptions, which may be lurking undetected at many ocean island volcanoes.more » « less
