skip to main content


Title: Crystal and melt inclusion timescales reveal the evolution of magma migration before eruption
Abstract

Volatile element concentrations measured in melt inclusions are a key tool used to understand magma migration and degassing, although their original values may be affected by different re-equilibration processes. Additionally, the inclusion-bearing crystals can have a wide range of origins and ages, further complicating the interpretation of magmatic processes. To clarify some of these issues, here we combined olivine diffusion chronometry and melt inclusion data from the 2008 eruption of Llaima volcano (Chile). We found that magma intrusion occurred about 4 years before the eruption at a minimum depth of approximately 8 km. Magma migration and reaction became shallower with time, and about 6 months before the eruption magma reached 3–4 km depth. This can be linked to reported seismicity and ash emissions. Although some ambiguities of interpretation still remain, crystal zoning and melt inclusion studies allow a more complete understanding of magma ascent, degassing, and volcano monitoring data.

 
more » « less
Award ID(s):
1664308
NSF-PAR ID:
10154353
Author(s) / Creator(s):
; ; ; ; ;
Publisher / Repository:
Nature Publishing Group
Date Published:
Journal Name:
Nature Communications
Volume:
9
Issue:
1
ISSN:
2041-1723
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. 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
  2. Abstract Magmatic systems are composed of melt accumulations and crystal mush that evolve with melt transport, contributing to igneous processes, volcano dynamics, and eruption triggering. Geophysical studies of active volcanoes have revealed details of shallow-level melt reservoirs, but little is known about fine-scale melt distribution at deeper levels dominated by crystal mush. Here, we present new seismic reflection images from Axial Seamount, northeastern Pacific Ocean, revealing a 3–5-km-wide conduit of vertically stacked melt lenses, with near-regular spacing of 300–450 m extending into the inferred mush zone of the mid-to-lower crust. This column of lenses underlies the shallowest melt-rich portion of the upper-crustal magma reservoir, where three dike intrusion and eruption events initiated. The pipe-like zone is similar in geometry and depth extent to the volcano inflation source modeled from geodetic records, and we infer that melt ascent by porous flow focused within the melt lens conduit led to the inflation-triggered eruptions. The multiple near-horizontal lenses are interpreted as melt-rich layers formed via mush compaction, an interpretation supported by one-dimensional numerical models of porous flow in a viscoelastic matrix. 
    more » « less
  3. null (Ed.)
    Abstract Despite multidisciplinary evidence for crustal magma accumulation below Santorini volcano, Greece, the structure and melt content of the shallow magmatic system remain poorly constrained. We use three-dimensional (3-D) velocity models from tomographic inversions of active-source seismic P-wave travel times to identify a pronounced low-velocity anomaly (–21%) from 2.8 km to 5 km depth localized below the northern caldera basin. This anomaly is consistent with depth estimates of pre-eruptive storage and a recent inflation episode, supporting the interpretation of a shallow magma body that causes seismic attenuation and ray bending. A suite of synthetic tests shows that the geometry is well recovered while a range of melt contents (4%–13% to fully molten) are allowable. A thin mush region (2%–7% to 3%–10% melt) extends from the main magma body toward the northeast, observed as low velocities confined by tectono-magmatic lineaments. This anomaly terminates northwest of Kolumbo; little to no melt underlies the seamount from 3 to 5 km depth. These structural constraints suggest that crustal extension and edifice loads control the geometry of magma accumulation and emphasize that the shallow crust remains conducive to melt storage shortly after a caldera-forming eruption. 
    more » « less
  4. Abstract

    Classical mechanisms of volcanic eruptions mostly involve pressure buildup and magma ascent towards the surface1. Such processes produce geophysical and geochemical signals that may be detected and interpreted as eruption precursors1–3. On 22 May 2021, Mount Nyiragongo (Democratic Republic of the Congo), an open-vent volcano with a persistent lava lake perched within its summit crater, shook up this interpretation by producing an approximately six-hour-long flank eruption without apparent precursors, followed—rather than preceded—by lateral magma motion into the crust. Here we show that this reversed sequence was most likely initiated by a rupture of the edifice, producing deadly lava flows and triggering a voluminous 25-km-long dyke intrusion. The dyke propagated southwards at very shallow depth (less than 500 m) underneath the cities of Goma (Democratic Republic of the Congo) and Gisenyi (Rwanda), as well as Lake Kivu. This volcanic crisis raises new questions about the mechanisms controlling such eruptions and the possibility of facing substantially more hazardous events, such as effusions within densely urbanized areas, phreato-magmatism or a limnic eruption from the gas-rich Lake Kivu. It also more generally highlights the challenges faced with open-vent volcanoes for monitoring, early detection and risk management when a significant volume of magma is stored close to the surface.

     
    more » « less
  5. Abstract

    Mount Cleveland is one of Alaska's most active volcanoes, yet little is known about the magmatic system driving persistent and dynamic volcanic activity. Volcanic gas and melt inclusion (MI) data from 2016 were combined to investigate shallow magmatic processes. SO2emission rates were between 166 and 324 t/day and the H2O/SO2was 600 ± 53, whereas CO2and H2S were below detection. Olivine‐, clinopyroxene‐, and plagioclase‐hosted MIs have up to 3.8 wt.% H2O, 514 ppm CO2, and 2,320 ppm S. Equilibration depths, based on MI H2O contents, suggest that a magmatic column extended from 0.5 to 3.0 km (~10–60 MPa). We used MI data to empirically model open‐system H‐C‐S degassing from 0 to 12 km and found that a column of magma between 0.5 and 3 km could produce the measured gas H2O/SO2ratio. However, additional magma deeper than 3 km is required to sustain emissions over periods greater than days to weeks, if the observed vent dimension is a valid proxy for the conduit. Assuming an initial S content of 2,320 ppm, the total magma supply needed to sustain the annual SO2flux was 5 to 9.8 Mm3/yr, suggesting a maximum intrusive‐to‐extrusive ratio of 13:1. The model predicts degassing of <50 t/day CO2for July 2016, which corresponds to a maximum predicted CO2/SO2of 0.2. Ultimately, frequent recharge from deeper, less degassed magma is required to drive the continuous activity observed over multiple years. During periods of recharge we would expect lower H2O/SO2and measurable volcanic CO2.

     
    more » « less