Fluid inclusion constraints on the geometry of the magmatic plumbing system beneath Mauna Loa – Part 2: Xenoliths

Penny E. Wieser; Matthew Gleeson; Berenise Rangel; Charlotte DeVitre; Alexander T. Bearden; Kendra J. Lynn; Paula Antoshechkina; Amy Gaffney; Brian Monteleone

doi:10.1007/s00445-025-01869-2

USGS · 70271918

Fluid inclusion constraints on the geometry of the magmatic plumbing system beneath Mauna Loa – Part 2: Xenoliths

Penny E. Wieser·Matthew Gleeson·Berenise Rangel·Charlotte DeVitre·Alexander T. Bearden·Kendra J. Lynn·Paula Antoshechkina·Amy Gaffney·Brian Monteleone

View original source ↗Export metadata ↓

Abstract

Mauna Loa volcano erupts crystal-poor material at its summit and more crystal-rich material on its rift zones. Some of the more olivine-rich lava flows contain xenoliths with diverse mineralogy, including cumulate harzburgites with high-Mg# orthopyroxenes and high-Fo olivines (both > 84). Previous experimental work and thermodynamic modelling has proposed that high-Mg# orthopyroxenes only crystallize from Mauna Loa melts at high pressures (> 6 kbar, > 20 km), leading to suggestions that there is a region of sub-Moho magma storage at Mauna Loa in addition to the geophysically imaged magma reservoir at 2–5 km depth below the summit. We use melt and fluid inclusion barometry combined with thermodynamic models to further investigate this suggestion. Fluid inclusion data from harzburgites and dunitic xenoliths yield storage depths remarkably similar to those found in non-xenolithic crystals from lavas and tephras, with a clear peak at ~ 2–3 km (below the summit). Depths from melt inclusions in these xenoliths overlap with fluid inclusion pressures, ruling out the possibility of fluid inclusion re-equilibration during a period of stalling in a shallower reservoir. We examine five different thermodynamic models and find that the minimum pressure of olivine-orthopyroxene co-saturation varies by ~ 4 kbar (~ 12 km). These models also fail to predict that orthopyroxene is stable in ~ 15–80% of compositionally relevant experimental charges which grew orthopyroxene. Overall, this shows that phase stability modelling is an unreliable method of determining magma storage depth at Mauna Loa. We suggest that model discrepancies reflect a lack of experimental constraints on orthopyroxene stability at > 1200 ℃ and 0.01–5 kbar. Based on the presence of large oikocrystic orthopyroxenes completely enclosing rounded olivine chadacrysts, we suggest that these harzburgitic xenoliths formed through the reaction of intruding melts with olivine mush piles within the Mauna Loa edifice at ~ 3 km depth below the summit, with no need for a deeper storage reservoir. The predominance of pre-eruptive shallow storage means that there is more chance of detecting reservoir destabilization with geophysical monitoring techniques compared to a scenario where melts are supplied from sub-Moho reservoirs.

Explore related subjects

Bulletin of Volcanology Article

Source-reported bounding extent: 18.997816607636395° to 19.719087637337836° latitude; -155.9431663120125° to -155.33750669977726° longitude. This indicates report coverage, not an exact sampling location. View area on OpenStreetMap.

Keep this discovery

Explore connections, maps & timelines

BibTeX RIS

Penny E. Wieser, Matthew Gleeson, Berenise Rangel, Charlotte DeVitre, Alexander T. Bearden, Kendra J. Lynn, Paula Antoshechkina, Amy Gaffney, Brian Monteleone. 2025-09-23. Fluid inclusion constraints on the geometry of the magmatic plumbing system beneath Mauna Loa – Part 2: Xenoliths. https://doi.org/10.1007/s00445-025-01869-2

Cite the original work for its findings. Save a collection to share your selection of sources.

Discover connections

Hawaii place United STates place Article subject Bulletin of Volcanology subject Mauna Loa place

Connections use source metadata and explicit phrase matches, not verified experimental comparisons.

KEEP EXPLORING

Related USGS reports

USGS2026

Geophysical insights into Kīlauea’s eruptive reawakening after its 2018 eruption and caldera collapse

A little more than two years after the destructive 2018 lower East Rift Zone eruption and summit collapse ended, renewed volcanic activity began at Kīlauea volcano, Hawaiʻi. During 2020–2025, activity at the volcano included six eruptions in the summit caldera, one in the Southwest Rift Zone, and one in the middle East Rift Zone, as well as multiple non-eruptive intrusions in all three locations. We categorize the volcano’s activity post-2018 into four phases based on eruption locations and style, as well as patterns of magma storage and migration inferred from geophysical data. In some cases, these observations illuminated parts of the volcano’s magmatic system for the first time since the advent of modern geophysical monitoring—particularly deformation datasets that include tilt, continuous GNSS, and InSAR. Analyses of these geodetic data and complementary seismic, geological, and geochemical observations provide an opportunity to explore summit and rift zone connectivity, rift-zone magma storage conditions, the variable geophysical manifestations of eruption onsets, magma-fault interactions, and other long-standing questions associated with eruptions and intrusions at Kīlauea.

Hawaii↗

USGS2026

An extensive and short-lived episode of mid-Holocene fayalite- bearing high-silica rhyolite eruptions at Mono Craters, California

Having last erupted ~ 700 years ago, the Mono Craters volcanic chain is an active system with a history of explosive and effusive high-silica rhyolite eruptions. Approximately ~ 30 lava domes and flows range from ~ 42 to ~ 0.7 ka, form a ~ 15 km curvilinear trend, and change mineralogy (biotite-, orthopyroxene-, or fayalite-bearing) and crystallinity through time. However, the eruptive chronology of this system is incompletely understood. We report new 238 U- 230 Th isochron ages for nine fayalite-bearing domes from paired zircon and allanite crystal surfaces that reveal a major shift in eruptive frequency and volume over the history of the chain, with extrusion of > 0.45 km 3 of high-silica rhyolite lava between 7.1 and 5.2 ka. This contrasts with older biotite- (~ 42–20 ka) and orthopyroxene-bearing (~ 12–9 ka) domes, each of which erupted smaller volumes over a more protracted period. Whole-rock compositions suggest these older domes tapped separate domains within a shared magmatic system because the former are higher in Sr and Ba whereas the latter are lower in SiO 2 , higher in Zr, and more enriched in light rare-earth elements. Titanomagnetite chemistry and zircon saturation temperatures also reflect changing magma compositions and conditions associated with these domes. Fayalite-bearing domes, in contrast, have uniform whole-rock and titanomagnetite compositions which support rapid extrusion of magma over an interval short enough to prevent Fe-Ti oxides from recording any changes in magmatic conditions. The youngest lavas were also produced in a short-lived episode ~ 700 years ago, showing that multivent, geologically brief eruptive episodes occur repeatedly at Mono Craters.

California↗

USGS2026

Satellite observations, plume dynamics and volcanic lightning from the 2023 eruption of Shishaldin Volcano, Alaska

The 2023 eruption of Shishaldin Volcano, Alaska, provided an unusually well-monitored example of a “dry” (non-hydromagmatic) explosive basaltic eruption with sustained, buoyant plumes, and pyroclastic density currents. Here, we examine the 13 basaltic paroxysms that occurred between July and November 2023 using geostationary satellite observations, plume modeling, and volcanic lightning detection to characterize mass eruption rates and plume dynamics. We show that seven volcanic plumes reached 10–14 km above sea level, sustained for minutes to tens of minutes. Translating plume heights into mass eruption rates, using both empirical scaling and one-dimensional plume modeling, gives peak rates of 10 4 –10 7 kg·s −1 and a total erupted volume of ~ 0.05 km 3 (dense rock equivalent) across the eruption sequence. Despite the moderate intensity of this eruption, automated satellite detection algorithms generally failed to detect volcanic ash. We attribute this persistent challenge to the fines-poor grain size distribution of the eruption. Volcanic lightning provides an additional window into the eruption dynamics. A total of 491 volcanic lightning flashes were detected, concentrated during seven of the most energetic phases. The timing of this electrical activity suggests that lightning only became detectable by the global lightning monitoring network when plumes ascended to atmospheric levels favorable for ice formation. Based on the available data, we infer that the low fine ash content of the plumes, combined with minimal water and ice, led to the overall modest rates of lightning (1–13 flashes per minute) compared to “wet” eruptions of similar intensity. These results show that combining satellite observations, plume modeling, and lightning detection provides a valuable timeline of eruptive intensity, even for “dry” basaltic eruptions in remote locations.

Alaska↗