Geology Reports⌕ Search

SEARCH · Geology Reports

Results for “Bulletin of Volcanology”

Search indexed USGS publications on groundwater, aquifers, geologic maps, mineral resources and earthquakes. Explore source records by subject and place.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 253 records · Page 14Linked to original sources
USGS2024

Pre-existing ground cracks as lava flow pathways at Kīlauea in 2014

In 2014, the Pāhoa lava flow at Kīlauea, on the Island of Hawaiʻi (USA), entered a string of pre-existing meter-width ground cracks in the volcano’s East Rift Zone. The ground cracks transported lava below the surface in a direction discordant to the slope of the landscape. The cracks, which were 100s of meters long and 10s to 100s of meters deep, also widened by up to several meters as they filled, probably in part at the expense of adjacent cracks, which likely closed. Widening of the cracks caused shallow crustal blocks on the volcano’s flank to shift—this deformation was captured by a nearby GPS station and a borehole tiltmeter. The GPS station moved away from the cracks in response, while the tiltmeter showed tilting toward the cracks, consistent with opening. Noting that the lava-filled cracks act as top-fed dikes, we adapt existing theory for the thermo-rheological evolution of dikes to analyze transport of lava captured by ground cracks and propose mechanisms for the exit of the lava back to the surface. This study shows that ground cracks as narrow as 50 cm wide can facilitate the transport of advancing lava flows and can carry lava in directions that differ from those expected based on surface topography, invalidating flow path projections based on the assumption of subaerial flow.

Hawaii↗
USGS2024

Tracking magma pathways and surface faulting in the Southwest Rift Zone and the Koaʻe fault system (Kīlauea volcano, Hawai ‘i) using photogrammetry and structural observations

Volcanic islands are often subject to flank instability, resulting from a combination of magmatic intrusions along rift zones and gravitational spreading causing extensional faulting at the surface. Here, we study the Koaʻe fault system (KFS), located south of the summit caldera of Kīlauea volcano in Hawaiʻi, one of the most active volcanoes on Earth, prone to active faulting, episodic dike intrusions, and flank instability. Two rift zones and the KFS are major structures controlling volcanic flank instability and magma propagation. Although several magmatic intrusions occurred over the KFS, the link between these faults, two nearby rift zones and the flank instability, is still poorly studied. To better characterize the KFS and its structural linkage with the surrounding fault and rift zones, we performed a detailed structural analysis of the extensional fault system, coupled with a helicopter photogrammetric survey, covering part of the south flank of Kīlauea. We generated a high-resolution DEM (~ 8 cm) and orthomosaic (~ 4 cm) to map the fracture field in detail. We also collected ~ 1000 ground structural measurements of extensional fractures during our three field missions (2019, 2022, and 2023). We observed many small, interconnected grabens, monoclines, rollover structures, and en-echelon fractures that were in part previously undocumented. We estimate the cumulative displacement rate across the KFS during the last 600 ~ 700 years and found a decrease toward the west of the horizontal component from 2 to 6 cm per year, consistent with GNSS data. Integrating morphology observations, fault mapping, and kinematic measurements, we propose a new kinematic model of the upper part of the Kīlauea’s south flank, suggesting a clockwise rotation and a translation of a triangular wedge. This wedge is bordered by the extensional structures (ERZ, SWRZ, and the KFS), largely influenced by gravitational spreading. These findings illustrate a structural linkage between the two rift zones and the KFS, the latter being episodically affected by dike intrusions.

Hawaii↗
USGS2024

Did steam boost the height and growth rate of the giant Hunga eruption plume?

The eruption of Hunga volcano on 15 January 2022 produced a higher plume and faster-growing umbrella cloud than has ever been previously recorded. The plume height exceeded 58 km, and the umbrella grew to 450 km in diameter within 50 min. Assuming an umbrella thickness of 10 km, this growth rate implied an average volume injection rate into the umbrella of 330–500 km 3 s −1 . Conventional relationships between plume height, umbrella-growth rate, and mass eruption rate suggest that this period of activity should have injected a few to several cubic kilometers of rock particles (tephra) into the plume. Yet tephra fall deposits on neighboring islands are only a few centimeters thick and can be reproduced using ash transport simulations with only 0.1–0.2 km 3 erupted volume (dense-rock equivalent). How could such a powerful eruption contain so little tephra? Here, we propose that seawater mixing at the vent boosted the plume height and umbrella growth rate. Using the one-dimensional (1-D) steady plume model Plumeria, we find that a plume fed by ~90% water vapor at a temperature of 100 °C (referred to here as steam) could have exceeded 50 km height while keeping the injection rate of solids low enough to be consistent with Hunga’s modest tephra-fall deposit volume. Steam is envisaged to rise from intense phreatomagmatic jets or pyroclastic density currents entering the ocean. Overall, the height and expansion rate of Hunga’s giant plume is consistent with the total mass of fall deposits plus underwater density current deposits, even though most of the erupted mass decoupled from the high plume. This example represents a class of high (> 10 km), ash-poor, steam-driven plumes, that also includes Kīlauea (2020) and Fukutoku-oka-no-ba (2021). Their height is driven by heat flux following well-established relations; however, most of the heat is contained in steam rather than particles. As a result, the heights of these water-rich plumes do not follow well-known relations with the mass eruption rate of tephra.

Bulletin of Volcanology↗
USGS2024

Real-time lava flow forecasting during the 2022 Mauna Loa eruption response

On November 27, 2022, Mauna Loa (Hawai‘i) erupted for the first time in 38 years, initially producing lava flows that covered the floor of its summit caldera, Moku‘āweoweo. Over the first 12 h following the summit eruption, four main fissures opened on Mauna Loa’s Northeast Rift Zone, with “fissure 3” quickly becoming the dominant source of lava flows. For the next 12 days, fissure 3 produced a 19-km-long lava flow to the north, crossing the Mauna Loa Weather Observatory access road and coming within 2.8 km of inundating the Daniel K. Inouye Highway (Saddle Road). Within 40 min of fissure 3 opening, inundation modeling efforts had begun. For the duration of the eruption, the computational fluid dynamics model Lava2d was run in real time, using flow front locations and other field observations to make sequential forecast improvements, eventually producing a set of models which accurately predicted the routing and arrival times of lava from fissure 3. These models were used to inform timing estimates of possible future inundation of the Saddle Road. As the eruption progressed, almost 4000 Lava2d models were made using high-performance computing resources, providing critical information on uncertainty in multi-week forecasts. To our knowledge, this was the first ever real-time physics-based ensemble lava flow modeling and forecasting effort. Here, we present a chronology of these real-time efforts, focusing on the successes and limitations of this approach.

Hawaii↗
USGS2024

Deformation of Mauna Loa volcano before, during, and after its 2022 eruption

Mauna Loa on the Island of Hawaiʻi erupted on 27 November 2022, the first eruption since 1984, which marked the culmination of decades-long period of non-eruptive unrest and relative quiescence. We briefly describe the evolution of the USGS Hawaiian Volcano Observatory’s geodetic monitoring network at Mauna Loa and show patterns of deformation as measured by Global Navigation Satellite Systems (GNSS), interferometric synthetic aperture radar (InSAR), and borehole tilt. We highlight the long-term buildup and the imminent pre-eruptive geodetic signals, including subtle changes observed in early 2021 that suggested a potential eruption. We then describe the significant ramp up of activity in September 2022 that provided strong evidence of likely impending eruptive activity. Of particular importance are the first borehole tilt excursions related to magma movement measured at Mauna Loa’s summit, which began in 2021 and were accompanied by increased rates of seismicity. In addition to describing the evolution of surface displacements, we also model the co-eruption deformation, which can be fit by dike-like opening that matches the geometry of the surface eruptive fissures.

Hawaii↗
USGS2024

Seismic velocity changes from repetitive seismicity at Mauna Loa prior to and during its 2022 eruption

Mauna Loa’s short-lived eruption from late November to early December 2022 marked the culmination of nearly a decade of elevated seismic activity and geodetic inflation. The volcano has been monitored by a network of permanent, short period and broadband seismometers. I used the continuous waveform data from that network starting in 2012 to generate a catalog of seismicity that enhances the US Geological Survey Hawaiian Volcano Observatory’s public seismic catalog with four times the number of earthquakes, which were then grouped by waveform similarity. Analysis of subtle delays in the timing of arrivals of scattered waves between pairs of earthquakes in this catalog yields a history of small changes in the shallow seismic velocity structure of the volcano. Seismic velocities have been shown at other volcanoes to change during unrest and eruption. My results show a decrease in seismic velocity centered on the summit beginning in September 2022, corresponding to the onset of a vigorous precursory swarm of seismic activity and shallow inflation. During the eruption itself, I observe large changes due likely to dike opening along the northeast rift zone and deflation of the summit reservoir. However, seismic velocity changes associated with non-volcanic sources such as ground shaking from large earthquakes and meteorological influences at seasonal and diurnal time scales are also observed, and these dominate the velocity changes prior to the eruption. Proper accounting of these effects will be a requirement for use in real-time monitoring, and this work serves as a starting point in that endeavor for Mauna Loa.

Hawaii↗
USGS2025

Local, regional, and distal recordings of seismic unrest at Tau Island volcano, American Samoa

A seismic swarm near Taʻū Island, a volcanic island in eastern American Samoa, occurred from July to October 2022. The earliest unrest was noted as felt shaking reports in late July, and instrumentation varied in the beginning of the sequence as the U.S. Geological Survey (USGS) Hawaiian Volcano Observatory responded by installing temporary and then permanent seismometers to monitor the activity. This network variability made it difficult to characterize the earliest seismicity and contextualize the entire sequence to discriminate between an underlying tectonic or volcanic source. Here, we present results analyzing hydroacoustic detections from an International Monitoring System hydrophone array near Wake Island, 4500 km northwest of Taʻū Island volcano. Using least-squares beamforming analysis, we create a catalog of T-wave detections from the direction of Taʻū Island to track the earthquakes, some of which were located by the USGS National EarthquakeInformation Center. Both the rate and hydroacoustic pressures, which we interpret as a proxy for earthquake size, gradually increased from late July to August, peaking on August 19 (rate) and August 24 (size), before decreasing to background in late September. Minutes-long bursts of tremor were also contemporaneously recorded as local network data became avail-able on August 20. Tremor activity continued throughout the rest of August, peaking on August 25, before ending in earlySeptember. These tremor bursts were band-limited to ~ 1–5 Hz and recorded as S waves at a regional station on the island of Upolu in Samoa, 250 km to the west of Taʻū Island. Our results do not constrain the tremor locations, but comparisons of earthquake and tremor reduced displacements recorded locally and regionally suggest a deeper tremor source. We interpret the increase in earthquake size and rate, together with the occurrence, characteristics, and relative depth of the tremor to be the result of magmatic activity beneath Taʻū Island volcano.

American Samoa, Taʻū Island↗
USGS2025

Complex staged emplacement of a basaltic lava: The example of the July 1974 flow of Kīlauea

Basaltic lava flows can be highly destructive. Forecasting the future path and/or behavior of an active lava flow is challenging because topography is often poorly constrained and lava has a complex rheology and emplacement history. Preserved lavas are an important source of information which, combined with observations of active flows, underpins conceptual models of lava flow emplacement. However, the value of preserved lavas is limited because pre-eruptive topography and, thus, syn-eruptive lava flow geometry are usually not known. Here, we use tree-mold data to constrain pre-eruptive topography and syn-eruptive lava flow geometry of the July 1974 flow of Kīlauea (USA). Tree molds, which are formed after advancing lava encloses standing trees, preserve the lava inundation height and the final preserved thickness of lava. We used data from 282 tree molds to reconstruct the temporal and spatial evolution of the ~ 2.1 km-long July 1974 flow. The tree mold dataset yields a detailed dynamic picture of staged emplacement, separated by intervals of ponding. In some ponded areas, flow depth during emplacement (~ 5 m) was twice the preserved thickness of the final lava (2–3 m). Drainage of the ponds led to episodic surges in flow advancement, decoupled from fluctuations in vent discharge rate. We infer that the final breakout occurred after the cessation of fountaining. Such complex emplacement histories may be common for pāhoehoe lavas at Kīlauea and elsewhere in situations where the terrain is of variable slope, and/or where lava is temporarily perched and stored.

Hawaii↗
USGS2025

Characteristics of volcanic ash reveal changes in fragmentation and eruption dynamics at Poás volcano, Costa Rica, 2016–2019

Determining whether fresh magma has reached the surface during a volcanic eruption can provide important information for forecasts of future activity, especially in the early stages of an eruption. However, identifying fresh, juvenile pyroclasts in tephra fall deposits can be challenging and inconclusive. We studied the products of explosions at Poás volcano, Costa Rica, in 2016–2019, a period during which the volcano transitioned from a pressurized, hydrothermally sealed state to an open conduit with increased degassing to the atmosphere. The activity consisted of semi-continuous explosions producing <500-m-high plumes, with the exception of explosions on April 14 and 22, 2017, that produced 4-km-high plumes. We analyzed the grain size distribution, componentry, and particle density of the products of twenty explosions, and collected groundmass glass composition on juvenile particles for three of them. Our work demonstrates varying degrees of magma involvement with the hydrothermal system through time, with juvenile material representing a wide range of abundance (~10–70 vol.%) in deposits of individual explosions. Before early April 2017, we infer that small phreatomagmatic explosions were triggered by contact between magmatic fluids and/or magma and the hydrothermal system, based on the presence of abundant hydrothermal fragments and minor juvenile magma. Concurrent with decreasing hydrothermal component in the deposits, explosions in April–May 2017 eroded the walls of the shallow plumbing system, evidenced by an increase in wall-rock lithics in the deposits. These changes coincided with drying of the crater lake, leading to several magmatic explosions in April-September 2017, whose juvenile-rich deposits are consistent with primary fragmentation of fresh magma in the conduit. The eruptive activity changed after this magmatic phase, and in 2019, small explosions mostly recycled the heterogeneous deposits of previous events, producing fine material with high proportions of recycled particles. All explosions from 2016–2019 remobilized already-emplaced magma from the 1953–1955 eruptive period, although distinct glass compositions between explosions suggest difference in crystallization or they tapped different portions of this magma body. Our work sheds light on the eruption dynamics and shallow plumbing system of this persistently active volcano and provides a case study for understanding the variable efficiency of phreatomagmatic fragmentation.

Poás volcano↗
USGS2025

Rhenium-osmium and oxygen isotope homogeneity during the 2022 Mauna Loa eruption and implications for basaltic magma storage

Mauna Loa is one of the largest and most active volcanoes on Earth. The most recent eruption of Mauna Loa started on 27 November 2022, lasted for 13 days, and was preceded by the longest repose time of 38 years in its modern history. In this contribution, new trace- and highly siderophile-element (HSE: Os, Ir, Ru, Pt, Pd, Re) abundances, 187 Re- 187 Os, and 18 O/ 16 O data are reported for the 2022 lavas. These lavas have a limited range of MgO (6.2 ± 0.1 wt.%) and Ni (83 ± 2 µg/g), with a broader range of Re (0.3 to 1.3 ng/g) and consistent Os (0.031 to 0.080 ng/g) contents. They have 187 Os/ 188 Os ratios (0.1345 to 0.1385) which are, on average, more radiogenic than Mauna Loa picrites (0.1331 to 0.1349) and are similar in composition to more differentiated Mauna Loa tholeiite lavas (0.1340 to 0.1381). The oxygen isotope compositions of glassy samples are 5.35 ± 0.15‰ ( n = 13) and span a range in δ 18 O of + 5.0 to + 5.5‰, with an average composition 0.2 to 0.3‰ lower than MORB. The δD value is − 81 ± 11‰ ( n = 5) at very low (0.03 ± 0.015 wt.%) H 2 O concentrations. The 2022 Mauna Loa eruption is similar in terms of δ 18 O but contrasts in terms of 187 Os/ 188 Os variability, with the recent longer-lived eruptions on La Palma (Canary Islands; 85 days) in 2021 and on the Reykjanes Peninsula (Iceland) that began in 2021 and are still ongoing. Initial lavas were more fractionated for both the Canary Islands and Iceland eruptions, producing more radiogenic Os isotope compositions than later erupted products. The 2022 Mauna Loa eruption showed no such trends. The limited range in isotope compositions of the 2022 Mauna Loa lavas and their strongly fractioned HSE patterns reflect long-term storage, crystal fractionation, and assimilation of related basaltic volcanic edifice materials by the parent magma beneath the volcano prior to eruption triggering. Eruption of differentiated and homogeneous tholeiite lavas at the summit caldera and high on the volcano’s flank, with emplacement of accumulative picrites lower on the volcano, are consistent with neutral buoyancy arguments.

Hawaii↗
USGS2025

Onset and evolution of summit lava fountaining during the Mauna Loa 2022 eruption

The start of the Mauna Loa 2022 eruption in the Mokuʻāweoweo summit caldera was entirely captured through webcam videos. We analyzed footage from the ~ 7-h summit episode, processing > 87,000 frames using a newly automated method to measure fountain heights, fissure lengths, and inflight ejecta volumes. The summit episode comprised four phases. In Phase 1 (~ 1 h ), a ~ 1 km long fissure propagated from southwest to northeast, with steady fountain heights reaching 125 m. In Phase 2 (~ 50 min ), two further fissure segments extended the total length to ~ 2.4 km, with fountains focused into point sources reaching up to 70 m. In Phase 3 (~ 70 min ), the eruption was steady with no major changes in eruptive behavior. In Phase 4 (~ 4 h ), fountains became unsteady and weak (< 30 m). Fissure growth rates varied from ~ 20 to -9 m/min, and bulk magma flux peaked at 2.7 × 10 4 m 3 /s. An inverse relationship between maximum fountain heights and fissure lengths suggests that total magma flux was nearly constant until Phase 4. We propose an interconnected feeding system for the summit episode with a preferential pathway for gas-rich magma in the southwest portion of the caldera. The sustained nature of the episode was likely driven by the fast rise of melt coupled with small bubbles, with trains of decoupled larger bubbles producing pulsations. Webcam videography was an effective tool to quantify early-stage eruption parameters and could be further explored for rapid response and fine-tuning of hazard mitigation strategies on erupting basaltic volcanoes.

Hawaii↗
USGS2025

New insights into gas-driven phase segregation in andesitic enclaves from Mt. Mazama (Crater Lake), USA

A key process in active magmatic systems is the “recharge” of deep-sourced mafic magma into cooler, more evolved, and crystal-rich shallow reservoirs; recharge may be the cause of, or response to, eruptive activity. Although compositional evidence for recharge has been extensively documented, physical models of recharge are limited, particularly processes that separate exsolving volatiles and melts from rapidly growing crystals. To improve constraints on phase separation behaviors, we re-examine andesitic enclaves in silicic andesite lava flows of Mt. Mazama (Crater Lake), USA, that provided early evidence of gas-driven filter pressing (Bacon, 1986). 2D and 3D imaging shows that enclaves have a sample-spanning crystal framework that is disrupted by melt patches, indicating that initially deformable crystal networks were subject to early phase reorganization. Small enclaves are poorly vesicular and require early gas loss. Large enclaves have porous cores with angular (diktytaxitic) voids that are well-connected in 3D and denser rinds with isolated pores. Large enclave rinds have similar bulk compositions to small enclaves but their less evolved cores require ~ 20% melt removal. In the large enclave, diktytaxitic core textures and gas fingering structures at the core–rind boundary suggest relatively slow late-stage outward gas migration. Both scaling arguments and evidence of outward gas/melt migration require a resistant rind. Rind formation is best explained by differential cooling and demonstrates the importance of thermal gradients for gas-driven filter pressing. A corollary is a limited time scale of recharge, enclave formation, and vesiculation to produce diktytaxitic textures, suggesting that recharge was (near) synchronous with eruption.

Oregon↗
USGS2025

Magma fragmentation and tephra deposition from a small-volume phreatomagmatic eruption: Blue Lake crater, Oregon, USA

Maars pose considerable hazards due to their more explosive nature (compared with more common scoria cones) and likelihood that eruptions produce pyroclastic surges. Blue Lake crater is a maar in the Oregon High Cascades that erupted within the last 3000 years, making it one of the youngest eruptions in the Oregon Cascades. Its young, unaltered deposits make it an excellent site to examine the relationship between fragmentation processes and ash characteristics. This paper presents an extensive data set of grain size, componentry, texture, particle morphology, and surface features for 23 samples from 17 layers from Blue Lake crater to better understand fragmentation style and eruptive dynamics over the course of the eruption. We present detailed stratigraphy from 22 tephra pits and analyze tephra samples following a standardized method. An improved isopach map and a new isopleth map show the extensive, ENE-trending fallout and surge deposits. Based on the tephra sheet, the eruption can be divided into three phases, starting with a phreatomagmatic phase that produced laterally extensive, lithic-rich fallout deposits and excavated the initial crater. The middle phase of the eruption produced phreatomagmatically fragmented intercalated fallout and surge deposits. The eruption closed with coarse grained fallout deposits with a declining lithic content, indicating a shift towards a hybrid or phreato-Strombolian style. This detailed examination of the deposits leads to a more nuanced explanation of the eruption and fragmentation dynamics, which contribute to a better understanding of maar eruption processes and hazards.

Oregon↗
USGS2025

A decade of precursory seismicity at Mauna Loa: 2013–2022

The 2022 eruption of Mauna Loa was preceded by nearly 10 years of seismic unrest. Using the earthquake catalog of the Hawaiian Volcano Observatory, we describe and explore many aspects of unrest between January 1, 2013, and the start of the eruption on November 28, 2022 (UTC). Earthquakes associated with unrest occurred above an inflating dike at shallow depths under and to the south of Moku‘āweoweo and in the cluster of seismicity to the west of Moku‘āweoweo associated with a buried rift zone. Earthquakes also occurred in a broad region at the base of the crust, likely associated with a zone of magma accumulation. Unrest is broken down into four phases of increasing intensity with time and into individual geographic regions. While activity occurring at shallow depths under the caldera often alternated with activity on the west flank in the first two phases of unrest, the third phase of unrest featured increased rates of earthquake activity in both areas acting concurrently. We also document a slow shift in seismicity across multiple years from the south of Moku‘āweoweo to an area under the caldera that likely represents a change in the location of the intrusion through time. Broadly speaking, unrest observed in advance of the 2022 eruption occurred in the same areas as the unrest preceding eruptions in 1975 and 1984. Furthermore, precursory unrest prior to eruptions in 1975, 1984 and 2022 lasted years and culminated in an intense phase of seismic activity lasting on the order of hours that represented the final dike propagation to the surface. Concurrent earthquake activity under Moku‘āweoweo and the west flank of Mauna Loa, paired with a shift in earthquakes under Moku‘āweoweo were the two most notable long-term precursors to an eruption; however, neither provided an actionable timeline of an impending eruption. The decade of precursory seismicity prior to the 2022 eruption is similar to long precursory phases prior to eruptions in 1975 and 1984 and underscores the challenges of informing on eruption hazards in a timely manner. By combining seismic activity with other geophysical fields such as geodesy and gravity, we present a conceptual model to aid in the interpretation of future unrest at Mauna Loa.

Hawaii↗
USGS2025

Growth of a passive lava lake during the 2020–2021 eruption of Kīlauea volcano, Hawaii

We investigate the growth of a passive lava lake in Halemaʻumaʻu crater during the December 2020 to May 2021 eruption of Kīlauea volcano, Hawaii. Fed by vents above their surfaces, the formation of passive lava lakes in topographic lows is an important process in the growth of basaltic volcanoes. We captured visible and thermal images during 17 helicopter overflights and applied structure-from-motion photogrammetry to create digital elevation models and orthomosaics of Halemaʻumaʻu. These data products allowed us to track eruptive activity and processes. The bulk time-averaged discharge rate (TADR) in December 2020 initially exceeded 100 m 3 s –1 but decreased to < 10 m 3 s –1 within seven days. By February 2021, TADR was < 2 m 3 s –1 and continued to decrease until the eruption ended in May 2021. A total volume of 40.6 ± 0.5 × 10 6 m 3 of lava filled Halemaʻumaʻu to a depth of 225 m. As TADR decreased, the lake progressively developed an immobile, solidified crust, beginning with surfaces farthest from the vent. This immobile surface rose endogenously, whereas exogenous surface rise occurred near the vent. Eruptive activity at a vent ended when the level of the lake surface exceeded that of the vent, which we attribute to the effects of lava sitting above the vent on ascending magma. Regular helicopter overflights, combined with field observations and the extensive monitoring network at Kīlauea, generated an unprecedented density of observations that provide insights into the emplacement of passive lava lakes and how these eruptions wane and end.

Hawaii↗
USGS2025

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

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.

Hawaii↗
USGS2025

Fluid inclusion constraints on the geometry of the magmatic plumbing system beneath Mauna Loa – Part I: Lavas and tephras

There are few petrological constraints on magma storage depths at Mauna Loa, Hawai‘i. Yet understanding the geometry of the magmatic plumbing system is critical for interpreting geophysical signals of unrest at this very high-threat volcano. We address this gap by examining CO₂-rich fluid inclusions (FI) in lava and tephra from seven eruptions (8870 ± 56 14C yr BP, 1852, 1855, 1868, 1949, 1950, and 1984), supplemented with published data from 2022. Raman spectroscopy was used to determine FI densities, from which entrapment pressures were calculated using a CO₂-H₂O equation of state. Most FI record pressures of ~ 0.25–1.25 kbar (~ 2–5 km depth below the summit), consistent with geophysical estimates from the past 40 years. In summit eruptions, FI hosted in more evolved olivine and orthopyroxene clots (Fo and Mg# < 84) record slightly shallower pressures than those in more primitive olivines (Fo > 84) from rift zone eruptions, suggesting a crystal-poor evolved cap near the top of the reservoir (~ 2 km). The similarity in storage depths across all eight eruptions indicates that Mauna Loa’s magmas have tapped a quasi-stable reservoir over the past two centuries, and possibly over 10 kyr. Electron backscatter diffraction reveals deformations to the crystal lattice in Fo82-83 olivines, likely due to deformation during storage in mush piles. The intensity of deformation is comparable to that seen at Kīlauea, implying that mush pile stress may be decoupled from edifice size or longevity. Finally, SO₂ contents in FI increase from ~ 2 mol% at 2 kbar to ~ 15 mol% at 0.5 kbar, suggesting sulphur degassing begins far deeper than the 0.2–0.3 kbar commonly assumed for Hawaiian systems. This validates the newest generation of S degassing models (e.g., Sulfur_X), and explains precursory SO2 emissions in the ~3 hours prior to the onset of the 2022 eruption (Esse et al. 2025).

Hawaii↗