Attention:The NSF Public Access Repository (PAR) system and access will be unavailable from 5:00 PM ET until 8:00 PM ET on Friday, September 11 due to maintenance. We apologize for the inconvenience.


Title: Advances in seismic imaging of magma and crystal mush
Seismic imaging methods have provided detailed three-dimensional constraints on the physical properties of magmatic systems leading to invaluable insight into the storage, differentiation and dynamics of magma. These constraints have been crucial to the development of our modern understanding of magmatic systems. However, there are still outstanding knowledge gaps resulting from the challenges inherent in seismic imaging of volcanoes. These challenges stem from the complex physics of wave propagation across highly heterogeneous low-velocity anomalies associated with magma reservoirs. Ray-based seismic imaging methods such as travel-time and surface-wave tomography lead to under-recovery of such velocity anomalies and to under-estimation of melt fractions. This review aims to help the volcanologist to fully utilize the insights gained from seismic imaging and account for the resolution limits. We summarize the advantages and limitations of the most common imaging methods and propose best practices for their implementation and the quantitative interpretation of low-velocity anomalies. We constructed and analysed a database of 277 seismic imaging studies at 78 arc, hotspot and continental rift volcanoes. Each study is accompanied by information about the seismic source, part of the wavefield used, imaging method, any detected low-velocity zones, and estimated melt fraction. Thirty nine studies attempted to estimate melt fractions at 22 different volcanoes. Only five studies have found evidence of melt storage at melt fractions above the critical porosity that separates crystal mush from mobile magma. The median reported melt fraction is 13% suggesting that magma storage is dominated by low-melt fraction crystal mush. However, due to the limits of seismic resolution, the seismological evidence does not rule out the presence of small (<10 km 3 ) and medium-sized (<100 km 3 ) high-melt fraction magma chambers at many of the studied volcanoes. The combination of multiple tomographic imaging methods and the wider adoption of methods that use more of the seismic wavefield than the first arriving travel-times, promise to overcome some of the limitations of seismic tomography and provide more reliable constraints on melt fractions. Wider adoption of these new methods and advances in data collection are needed to enable a revolution in imaging magma reservoirs.  more » « less
Award ID(s):
2023338
PAR ID:
10440842
Author(s) / Creator(s):
; ; ; ; ;
Publisher / Repository:
Frontiers in Earth Science
Date Published:
Journal Name:
Frontiers in Earth Science
Volume:
10
ISSN:
2296-6463
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Quantitative estimates of magma storage are fundamental to evaluating volcanic dynamics and hazards. Yet our understanding of subvolcanic magmatic plumbing systems and their variability remains limited. There is ongoing debate regarding the ephemerality of shallow magma storage and its volume relative to eruptive output, and so whether an upper-crustal magma body could be a sign of imminent eruption. Here we present seismic imaging of subvolcanic magmatic systems along the Cascade Range arc from systematically modelling the three-dimensional scattered wavefield of teleseismic body waves. This reveals compelling evidence of low-seismic-velocity bodies indicative of partial melt between 5 and 15 km depth beneath most Cascade Range volcanoes. The magma reservoirs beneath these volcanoes vary in depth, size and complexity, but upper-crustal magma bodies are widespread, irrespective of the eruptive flux or time since the last eruption of the associated volcano. This indicates that large volumes of melts can persist at shallow depth throughout eruption cycles beneath large volcanoes. 
    more » « less
  2. 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
  3. 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
  4. 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
  5. Abstract From 2019 to 2024, gravity surveys were conducted at the Three Sisters volcanic cluster (TSVC), measuring 246 gravity sites using a spring relative gravimeter. We calculated the residual Bouguer anomaly and identified three main zones with negative anomalies, ranging from −4 to −8 mGal, located southwest and west of South Sister, within an area that has been uplifting for the past two decades. After inversion, we obtain a 3D density model of the subsurface and identify low‐density bodies extending from the surface down to 3 km. We estimate a total of 15 k of crustal bodies with density close to 2 g/ that could store up to 5 k of water, forming an extensive hydrothermal system beneath the TSVC. We explore the possible combinations of melt compositions and temperatures that could create a bulk density close to our reference crustal density (2.5 g/) using MELTS thermodynamic simulations. Our results indicate that a magmatic mush with as little as 15% partial melt of bulk rhyolitic composition or as much as 52%–57% partial melt of a bulk dacitic composition could be stored in a magmatic system under TSVC without generating a detectable gravity anomaly. Episodic magma injections at the base of the magmatic system, such as the 1998–2000 intrusion at 6 km depth, would bring heat and gas to the hydrothermal system while maintaining a low melt fraction in the magmatic mush, as imaged at other Cascade volcanoes. 
    more » « less