A Geodetically Constrained Petrogenetic Model for Evolved Lavas from the January 1997 Fissure Eruption of Kīlauea Volcano

Abstract

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.