On November 28th, 2022, following a record historical repose period of 38 years, Mauna Loa erupted about 145 × 106 m3 of lava and tephra over a 15-day period. The eruption was confined to the summit caldera region and the upper Northeast Rift Zone and is remarkably homogeneous in composition in both time and space. In these respects, it is typical of prior shallow summit reservoir magma bodies, recently estimated to be at a depth of around 1–2 km beneath the caldera. In contrast with these earlier magma bodies, which typically contain 6.7–7.1

More than 50 explosive eruptions occurred from Halemaʻumaʻu at Kīlauea volcano over 17 days from May 11 to 27, 1924. Ballistics weighing as much as 14,000 kg were ejected and most landed within 2 km of the vent. Fine ash made up a major component of the tephra and was dispersed tens of kilometers downwind. Draining of the Halemaʻumaʻu lava lake occurred in late February 1924, with the crater floor eventually subsiding by a further ∼70 m (to ∼180 m below the crater rim) by the time the first explosions took place during the night of May 10–11. The largest explosions occurred on May 17–18 and smaller explosions continued until May 27, at which point Halemaʻumaʻu had more than doubled in width and depth. The explosions generated plumes reaching up to ∼10 km high with ballistics ejected up to 2 km from the crater.Almost 100 years later, we investigate and characterize the preserved tephra deposits within ∼3 km of the 1924 crater rim. Grain size and shape analyses were performed on 202 samples collected from 34 tephra profiles using dynamic image analysis, with a subset of layers from nine tephra profiles used for componentry (200 grains per layer in the 0.5–1 mm size fraction). Additionally, we characterize the average diameters (using the five largest clasts) at 216 locations and measure the average diameters of 2291 ballistics (largest per ∼100 m2 area). Physical descriptions from fieldwork and grain size distributions were used to subdivide the tephra layers into five lithofacies: coarse homogeneous, fine homogenous, red ash, accretionary lapilli-bearing, and finely laminated. Grain size versus shape data show a range of values that demonstrate most grains are dense, smooth, and equant, in alignment with lithic clasts dominating the tephra componentry. The fine grained and accretionary lapilli-bearing nature of some of these lithofacies confirms that water influenced the style of the explosions. However, we also note juvenile clasts within many of the tephra layers, indicating that many of the layers were formed during phreatomagmatic explosions (sensu stricto), despite the eruptive mechanism being dominantly phreatic. Juvenile clasts are more abundant higher in the tephra profiles, suggesting that juvenile magma was more involved later in the explosive sequence. Thermal and hydrologic modeling indicate that groundwater inflow into a short-lived, small-diameter volcanic conduit (10-m to 120-m-diameter used for modeling) during the 78–85 days preceding the first explosion provides a physically plausible mechanism for this eruptive sequence.

More than 50 explosive eruptions occurred from Halema'uma'u at K & imacr;lauea volcano over 17 days from May 11 to 27, 1924. Ballistics weighing as much as 14,000 kg were ejected and most landed within 2 km of the vent. Fine ash made up a major component of the tephra and was dispersed tens of kilometers downwind. Draining of the Halema'uma'u lava lake occurred in late February 1924, with the crater floor eventually subsiding by a further similar to 70 m (to similar to 180 m below the crater rim) by the time the first explosions took place during the night of May 10-11. The largest explosions occurred on May 17-18 and smaller explosions continued until May 27, at which point Halema'uma'u had more than doubled in width and depth. The explosions generated plumes reaching up to similar to 10 km high with ballistics ejected up to 2 km from the crater. Almost 100 years later, we investigate and characterize the preserved tephra deposits within similar to 3 km of the 1924 crater rim. Grain size and shape analyses were performed on 202 samples collected from 34 tephra profiles using dynamic image analysis, with a subset of layers from nine tephra profiles used for componentry (200 grains per layer in the 0.5-1 mm size fraction). Additionally, we characterize the average diameters (using the five largest clasts) at 216 locations and measure the average diameters of 2291 ballistics (largest per similar to 100 m(2) area). Physical descriptions from fieldwork and grain size distributions were used to subdivide the tephra layers into five lithofacies: coarse homogeneous, fine homogenous, red ash, accretionary lapilli-bearing, and finely laminated. Grain size versus shape data show a range of values that demonstrate most grains are dense, smooth, and equant, in alignment with lithic clasts dominating the tephra componentry. The fine grained and accretionary lapilli-bearing nature of some of these lithofacies confirms that water influenced the style of the explosions. However, we also note juvenile clasts within many of the tephra layers, indicating that many of the layers were formed during phreatomagmatic explosions (sensu stricto), despite the eruptive mechanism being dominantly phreatic. Juvenile clasts are more abundant higher in the tephra profiles, suggesting that juvenile magma was more involved later in the explosive sequence. Thermal and hydrologic modeling indicate that groundwater inflow into a short-lived, small-diameter volcanic conduit (10-m to 120-m-diameter used for modeling) during the 78-85 days preceding the first explosion provides a physically plausible mechanism for this eruptive sequence.

Kīlauea volcano underwent dramatic morphological changes in 2018. That year recorded the end of the 35-year-long eruption of Puʻuʻōʻō (1983–2018) and 10-year-long (2008–2018) Halemaʻumaʻu lava lake and emplacement of the 4-month-long lower East Rift Zone lava flows that coincided with 500 m of summit caldera collapse. Starting on December 20, 2020, eruptions resumed at Kīlauea’s summit. There were five summit eruptions between December 2020 and September 2023, which ranged in duration from more than a year to as short as a week. Following these summit eruptions, seismicity and deformation increased in the upper Southwest Rift Zone in 2024, culminating in a 8.5-h-long eruption in this region on June 3, 2024. Increased seismicity and deformation then shifted to the upper and middle East Rift Zone and after several months culminated in an eruption just west of, and within, Nāpau Crater in the middle East Rift Zone from September 15 to 20, 2024. Despite vast morphological changes at Kīlauea’s summit, the geochemical compositions (i.e., whole rock and glass) that erupted from December 2020 to September 2023 are all remarkably similar to each other. Whole-rock compositions appear distinct from the preceding 2008–2018 Halemaʻumaʻu lava lake and phase 3 (i.e., summit or uprift-derived mafic lavas) of the 2018 lower East Rift Zone lava flows, although glass compositions appear to have more overlap with 2018 lower East Rift Zone glasses. The June 3, 2024, upper Southwest Rift Zone spatter and lava flows exhibit a dramatic enrichment in whole-rock MgO that is not recorded in glass, which reflects accumulation of olivine (e.g., antecrysts or xenocrysts) during dike emplacement, and is consistent with the abundance of olivine in the lava flows (5–10

Differentiated magmas stored in the rift zones of Kilauea have received more attention in recent years following eruption of andesite during the early phase of 2018 lower East Rift Zone activity. Despite this growing interest, some of the most voluminous eruptions of differentiated rift zone magmas remain poorly studied. One such eruption, and the most voluminous exposed differentiated flow field at Kilauea, is the Kamakai'a Hills. This eruption took place in the Southwest Rift Zone of Kilauea, a region that is hypothesized to contain a long-lived rift zone reservoir. The Kamakai'a Hills flow field encompasses >250 & times; 10(6) m(3) of basaltic andesite and basalt compositions with a mineral assemblage of orthopyroxene + clinopyroxene + plagioclase during its early 'a'a phase and clinopyroxene + plagioclase + olivine during its late pahoehoe phase. To better understand storage conditions and magma accumulation, this study focuses on major, minor, and trace elements from the mineral assemblage present within the early 'a'a and late pahoehoe phases. The diversity of clinopyroxene and plagioclase compositions within the early 'a'a and late pahoehoe phases, as well as diverse compositions of plagioclase and orthopyroxene within the early 'a'a phase, suggest multiple magma bodies and limited pre-eruption magma mixing within the broader Kamakai'a Hills reservoir. Oscillatory zoning patterns (particularly in clinopyroxene) imply processes such as recharge events, magma mixing or mingling, or convection within a differentially cooling, chemically stratified reservoir over protracted time intervals, whereas only limited resorbed mineral textures indicate incomplete mixing of heat and chemically distinct magmas during the dike intrusion that triggered the eruption. Mineral-mineral and mineral-melt thermobarometry indicate predominantly shallow (<= 2.5 km depth) crustal storage conditions of the cooled, differentiated magma (similar to 1100 degrees C and cooler for the basaltic andesites) to hotter temperatures for the basalts (all >1100 degrees C). Despite the known large standard errors estimated for mineral-melt and mineral-mineral barometry (10s to >100 MPa), the calculated pressures and depths broadly correspond with earthquake swarm depths beneath the Kamakai'a Hills, and drill core and fluid inclusion barometry storage depths of differentiated magmas within the lower East Rift Zone. The Kamakai'a Hills differentiated magmas have H2O contents (similar to 0.5 wt%, using plagioclase-melt hygrometry) equivalent to typical Kilauea basalts. Our data and interpretations demonstrate a complex, long-lived rift zone storage system that consisted of multiple magma bodies and was mobilized into eruption through intrusion of a hotter and more primitive summit-derived (uprift) magma.

Volcano-geologic mapping and stratigraphic reconstructions provide important information toward understanding patterns of edifice construction and destruction of volcanic systems, their eruptive histories and recurrence rates, and magmatic evolution and plumbing systems, all of which are required to make informed hazard assessments. Geologic mapping of volcanic terrains also provides context in the search for and identification of natural resources, including geothermal reservoirs and magmatic-related ore deposits, and can provide useful background for communication and outreach. Most volcano-geologic maps and chrono-stratigraphic frameworks are based on a systematic lithostratigraphic approach and are constructed using an iterative process whereby field observations and mapping alternate with acquisition of rock compositional and geochronological data. Modern technological and analytical tools have greatly advanced geochemical and geochronological data acquisition and accuracy, while some technologies have provided new useful tools for fieldwork, including in poorly accessible environments.

Details of Mauna Loa’s plumbing system and magma storage dynamics are relatively poorly constrained, relative to nearby Kīlauea volcano. Recent studies have begun to address this gap in knowledge using geophysical imaging or olivine-focused petrology. Here, we use the crystal cargo from the 2022 summit and Northeast Rift Zone eruption to advance our understanding of Mauna Loa’s storage regions and pre˗eruptive processes. We identify textural and geochemical populations of olivine, orthopyroxene, clinopyroxene, and plagioclase macrocrysts as well as evolved glomerocrysts consisting of orthopyroxene, clinopyroxene, and plagioclase. Summit versus rift zone trends can be distinguished between specific crystal populations and are used to infer magma transfer pathways. Melt inclusion barometry from three orthopyroxene populations indicate crystallization depths of 2–3 km which coincide with previous olivine˗hosted fluid inclusion barometry across Mauna Loa’s recent 200˗year eruptive history. We combine orthopyroxene˗hosted melt inclusion barometry with diffusion chronometry on the same crystals using Fe˗Mg interdiffusion to estimate timescales of mixing and recharge prior to the eruption. Summit orthopyroxene yield timescales of weeks to months, aligning with geophysically recognized pre˗eruptive intrusion events and previously determined olivine timescales. Northeast Rift Zone orthopyroxene crystals have resorbed cores and yield short pre˗eruptive mixing timescales on the order of hours to days. Chemistry and timescales from the orthopyroxene populations suggest that some of the 2022 magmas may have bypassed the shallowest portions of the summit reservoir system and were transferred directly into the rift zone syn˗eruptively.

From July to October 2022, a noneruptive volcanic earthquake swarm occurred within ∼15 km of Ta‘ū Island, located in eastern American Samoa. Felt reports from residents were the only available information about the swarm when it started, as American Samoa lacked a seismic monitoring network. We developed a consistent single-station catalog for the entire swarm, using seismic data from the nearest station IU.AFI, ∼250 km away. We applied the EQTransformer deep-learning model (Mousavi et al., 2020), automatically picking Pn and Sn arrivals on IU.AFI continuous data. We retained only events with Sn–Pn times of 22.5–25 s, consistent with the expected locations based on felt reports, then detected smaller swarm events with subsequent template-matching. This single-station catalog characterized the swarm’s onset and escalation to peak activity before a multiagency field response team installed a local seismic network in mid-August 2022. This permanent seismic network captured the swarm’s decline. EQTransformer identified short S–P times on the first two locally deployed seismometers, both Raspberry Shake sensors, to constrain the swarm’s distance from Ta‘ū Island. Modern seismological processing methods, combined with basic observations such as felt reports, can quickly contribute useful information during an earthquake response in a poorly monitored region.

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, 187Re-187Os, and 18O/16O data are reported for the 2022 lavas. These lavas have a limited range of MgO (6.2 ± 0.1 wt.

The Keaiwa Lava Flow of 1823 in the Southwest Rift Zone of Kilauea volcano is unusual for its expansive pahoehoe sheet flow morphology and lack of constructive vent topography, despite having a similar tholeiitic basalt composition to other lavas erupted from Kilauea. This lava flow issued from a similar to 10-km-long continuous fissure now known as the Great Crack, and has an unusually thin sheet flow morphology with margin thicknesses of similar to 15-110 cm (average of 42 cm). Based on field observations of the lava flow at its fissure vent (e.g., drainback features), we propose that the Great Crack formed, or at least significantly widened, just prior to and syneruptively with this 1823 eruption. The absence of pyroclastic cones or spatter ramparts indicates that the eruption consisted of a rapid outpouring of relatively degassed lava as the fissure unzipped. The rapidly moving lava flow overtopped pre-existing tumuli and scoria cones (e.g., Lava Plastered Cones) up to similar to 10 m tall. Glass and whole-rock chemistry yield homogeneous compositions for the lavas erupted from the Great Crack, with glass compositions of 6.40 +/- 0.10 wt% MgO and whole-rock compositions of 7.39 +/- 0.07 wt% MgO. Lava pads erupted from a short western fissure system are richer in mafic minerals (e.g., olivine and clinopyroxene), and show slightly more MgO-rich whole-rock compositions (7.79 +/- 0.05 wt%). MgO-in-glass thermometry on juvenile spatter yield eruption temperatures of 1153 +/- 13 degrees C that are typical of Kilauea lavas. Thus, the extensive sheet-like lava flow morphology is not a direct consequence of unusual magmatic or rheological conditions (i.e., low viscosity). Instead, the flow morphology is associated with high effusion rates caused by sudden drainage of uprift magma as it erupted from the Great Crack. Lava flow modeling on a 2-m-resolution digital elevation model indicates that a minimum bulk effusion rate of similar to 5800 m(3)/s (similar to 3500 m3/s dense rock equivalent) and a minimum flow velocity of similar to 11 m/s are required for the lava flow to overcome the topography of the Lava Plastered Cones. This effusion rate is among the highest inferred for eruptions in Hawai'i and around the world. This study highlights a less frequent eruption style at Hawaiian volcanoes characterized by a sudden outpouring of lava from an unusual fissure system. Local eyewitness accounts indicate that the 1823 eruption was preceded by seismicity. Given the complex magmatic-volcanic-tectonic relations across Kilauea, we speculate that the south flank could have slipped over one or more events that ultimately triggered unzipping of the Great Crack and passive release of briefly stored uprift magma. An eruption similar to 1823 at Kilauea or Mauna Loa, with an eruptive timeframe that could be as short as an hour, with high effusion rates and rapid flow front velocities, would not easily allow for a timely response.

We integrate new mineral geochemistry with field relations, 36Cl surface-exposure ages, whole-rock geochemistry, and paleomagnetic directions to understand the petrogenesis and eruption duration of three young basaltic lava flows in northern Harrat Rahat, a distributed volcanic field in central-western Saudi Arabia. These lavas, informally referred to as the Five Fingers, erupted at 24 ka from three elongate vent complexes that define a common lineament. Paleomagnetic data from the northern and central fingers are consistent with penecontemporaneous eruption, while the southern fingers erupted at a different time. Based on paleomagnetic data, the eruptive interval for this set of lava flows is estimated to be 160 ± 70 years. Mineral and whole-rock chemistry suggest that the southern fingers magma was deeply sourced and near-primitive, whereas the northern fingers magma had fractionated a gabbroic assemblage in the mid-crust prior to eruption. The central finger has a whole-rock composition nearly identical to the northern fingers but contains crystals compositionally similar to those present in both the southern and northern fingers, suggesting mixing between a deeply sourced magma and a magma stalled in the mid-crust. These data suggest that the Five Fingers eruptive episode reflects renewed injection of mantle-derived magma into the crust ca. 24 ka, in which basalts from different crustal levels interacted and erupted in two discrete events separated by 160 ± 70 years. Compared with similar datasets from other distributed volcanic fields, our results show that the duration of eruption episodes involving multiple eruption events can span from months to > 100 years.

Syn-eruption geochemical monitoring during volcanic activity is an important component of integrated volcanic monitoring. Volcanoes on the Island of Hawai‘i are primarily monitored by the U.S. Geological Survey’s Hawaiian Volcano Observatory using instrumental networks, field surveys, satellite observations, and petrologic monitoring. In collaboration with the University of Hawaiʻi at Hilo, an important component of this monitoring is near-real-time (generally within 24 h of sample collection) bulk-rock geochemistry using benchtop energy-dispersive X-ray fluorescence (EDXRF). Bulk-rock geochemistry, coupled with additional petrologic and geophysical analyses, are used to track eruption dynamics to help facilitate identification and understanding of changes that signal the onset of potentially hazardous activity. In this paper, we present EDXRF methodology for bulk-rock analysis used to evaluate the 2018 Kīlauea and 2022 Mauna Loa eruptions and compare this technique to other geochemical methods. Syn-eruption geochemical monitoring using EDXRF provides reliable quantitative data that is important for understanding eruption dynamics. Critical to this process are calibrating the instrument for analysis over the same range of compositions expected during volcanic activity, compiling a relevant database from past eruptions on the same instrument for comparison, and having the physical infrastructure and trained personnel to complete the analyses in a timely fashion. In addition to the scientific goals of petrologic monitoring, near-real-time geochemical analysis in Hawaii provides important information about the nature of on-going eruptions and informs decision-makers and the public about associated evolving hazards.

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 m3 s–1 but decreased to < 10 m3 s–1 within seven days. By February 2021, TADR was < 2 m3 s–1 and continued to decrease until the eruption ended in May 2021. A total volume of 40.6 ± 0.5 × 106 m3 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.

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.

Abstract Vesicularity of individual pyroclasts from airfall tephra deposits is an important parameter that is commonly measured at basaltic volcanoes. Conventional methods used to determine pyroclast vesicularity on a large number of clasts has the potential to be time consuming, particularly when rapid analysis is required. Here we propose dynamic image analysis on two-dimensional (2D) projection shapes of crushed pyroclasts from tephra deposits as a new method to estimate vesicularity. This method relies on the influence of vesicles and uses grain morphology as a proxy for vesicle size and abundance. Pyroclasts from a variety of basaltic tephra deposits from the volcanoes of Mauna Loa and Kīlauea were analyzed. Vesicularities between 52–98% were measured via nitrogen-gas pycnometry. The same pyroclasts were then crushed and sieved, and their grain shapes measured using dynamic image analysis on a CAMSIZER®. This yields values for the mean sphericity, elongation, compactness, and Krumbein roundness of the grains. Our data show that grains become increasingly irregular with increasing vesicularity, with the degree of correlation between shape parameters and vesicularity depending on the size of measured grains. Shape irregularities in small grains (60–250 µm) are mostly area-based, with elongation being the best vesicularity indicator, whereas shape irregularities in large grains (250–700 µm) are mostly perimeter-based, with Krumbein roundness as the best vesicularity indicator. Using mean shape parameter values with all grain sizes included, grain elongation is the most well-correlated shape parameter with vesicularity, with the best fitted model explaining 76% of variation in the observations. Microscope images of thin sections of intact pyroclasts, as well as from crushed pyroclasts, were analyzed using CSDCorrections 1.6 software in ImageJ to find local vesicularity, vesicle size, grain size, grain elongation, and vesicle spatial distribution by stereological conversion. Observed correlation between grain shape and vesicularity can be explained by the local effect of vesicles on the shape of the solid structure in between those vesicles. Grain shape depends not only on vesicularity, but also on vesicle to grain size ratio and the spatial distribution of vesicles. The influence of vesicles on grain shape is best captured by grains with the size of the solid structure in between vesicles, which generally increases with decreasing vesicularity. Dynamic image analysis is a useful tool to quickly gauge vesicularity, which could be used in near-real-time during an eruption response. However, this method is best suited for highly vesicular (> 80%) basaltic pyroclasts from tephra deposits with few microlites and phenocrysts. Further research on crushing techniques, optimum grain size for shape measurements, and Krumbein roundness measurements for the grain size range of 250–700 µm might enable application of this method to lower vesicularity pyroclasts.

Abstract Distinguishing periods of intermittent unrest from the run-up to eruption is a major challenge at volcanoes around the globe. Comparing multidisciplinary monitoring data with mineral chemistry that records the physical and spatio-temporal evolution of magmas fundamentally advances our ability to forecast eruptions. The recent eruption of Mauna Loa, Earth’s largest active volcano, provides a unique opportunity to differentiate unrest from run-up and improve forecasting of future eruptions. After decades of intermittent seismic and geodetic activity over 38 years of repose, Mauna Loa began erupting on 27 November 2022. Here we present a multidisciplinary synthesis that tracks the spatio-temporal evolution of precursory activity by integrating mineral and melt chemistry, fluid inclusion barometry, numerical modeling of mineral zoning, syn-eruptive gas plume measurements, the distribution and frequency of earthquake hypocenters, seismic velocity changes, and ground deformation. These diverse data indicate that the eruption occurred following a 2-month period of sustained magma intrusion from depths of 3–5 km up to 1–2 km beneath the summit caldera, providing a new model of the plumbing system at this very high threat volcano. Careful correlation of both the geochemistry and instrumental monitoring data improves our ability to distinguish unrest from the run-up to eruption by providing deeper understanding of the both the monitoring data and the magmatic system—an approach that could be applied at other volcanic systems worldwide.