skip to main content

Title: Climbing the crustal ladder: Magma storage-depth evolution during a volcanic flare-up
Very large eruptions (>50 km 3 ) and supereruptions (>450 km 3 ) reveal Earth’s capacity to produce and store enormous quantities (>1000 km 3 ) of crystal-poor, eruptible magma in the shallow crust. We explore the interplay between crustal evolution and volcanism during a volcanic flare-up in the Taupo Volcanic Zone (TVZ, New Zealand) using a combination of quartz-feldspar-melt equilibration pressures and time scales of quartz crystallization. Over the course of the flare-up, crystallization depths became progressively shallower, showing the gradual conditioning of the crust. Yet, quartz crystallization times were invariably very short (<100 years), demonstrating that very large reservoirs of eruptible magma were transient crustal features. We conclude that the dynamic nature of the TVZ crust favored magma eruption over storage. Episodic tapping of eruptible magmas likely prevented a supereruption. Instead, multiple very large bodies of eruptible magma were assembled and erupted in decadal time scales.  more » « less
Award ID(s):
Author(s) / Creator(s):
; ; ; ; ; ; ;
Date Published:
Journal Name:
Science Advances
Page Range / eLocation ID:
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract The long-term thermochemical conditions at which large bodies of silicic magma are stored in the crust is integral to our understanding of the timing, frequency, and intensity of volcanic eruptions and provides important context for interpreting volcano monitoring data. Despite this, however, individual magmatic systems may exhibit a range of time–temperature paths, or thermal histories, that are the result of many complex and, in some cases, competing processes. This complexity contributes to an incomplete understanding of the long-term thermal evolution of magma stored within the Earth’s crust. Of recent interest to the volcanology community is the length of time large volumes of rheologically eruptible and geophysically detectable magma exist within the crust prior to their eruption. Here we use a combination of diffusion chronometry, trace element, and thermodynamic modeling to quantify the long-term thermal evolution of the 2.08 Ma, 630 km3 Cerro Galán Ignimbrite (CGI) in NW Argentina; one of the largest explosive volcanic eruptions in the recent geologic record. We find that diffusion of both Mg and Sr in plagioclase indicate that erupted magmatic material only spent decades to centuries at or above temperatures (~750°C) required to maintain significant volumes of stored eruptible magma. Calculated plagioclase equilibrium compositions reveal an array of liquids that is controlled overall by fractionation of plagioclase + biotite + sanidine, although high-resolution trace element transects record a diversity of fractionation pathways. Overall, we suggest that there is compelling evidence that the magma erupted from the CGI magmatic system spent most of its upper crustal residence in a largely uneruptible state and was rapidly remobilized shortly before eruption. 
    more » « less
  2. Abstract

    Interpreting unrest at silicic volcanoes requires knowledge of the magma storage conditions and dynamics that precede eruptions. The Laguna del Maule volcanic field, Chile, has erupted ~40 km3of rhyolite over the last 20 ka. Astonishing rates of sustained surface inflation at >25 cm/year for >12 years reveal a large, restless system. Integration of geochronologic, petrologic, geomorphic, and geophysical observations provides an unusually rich context to interpret ongoing and prehistoric processes. We present new volatile (H2O, CO2, S, F, and Cl), trace element, and major element concentrations from 109 melt inclusions hosted in quartz, plagioclase, and olivine from seven eruptions. Silicic melts contain up to 8.0 wt. % H2O and 570 ppm CO2. In rhyolites melt inclusions track decompression‐driven fractional crystallization as magma ascended from ~14 to 4 km. This mirrors teleseismic tomography and magnetotelluric findings that reveal a domain containing partial melt spanning from 14 to 4 km. Ce and Cl contents of rhyolites support the generation of compositionally distinct domains of eruptible rhyolite within the larger reservoir. Heat, volatiles, and melt derived from episodic mafic recharge likely incubate and grow the shallow reservoir. Olivine‐hosted melt inclusions in mafic tephra contain up to 2.5 wt. % H2O and 1,140 ppm CO2and proxy for the volatile load delivered via recharge into the base of the silicic mush at ~14 to 8 km. We propose that mafic recharge flushes deeper reaches of the magma reservoir with CO2that propels H2O exsolution, upward accumulation of fluid, pressurization, and triggering of rhyolitic eruptions.

    more » « less
  3. Interpretation of erupted products we observe on the seafloor requires that we understand the petrogenesis of melts in the oceanic crust and where crystallization initially takes place. Our work focuses on estimating depth of crystallization of the plagioclase megacrysts using CO2 and H2O concentrations from plagioclase ultraphyric basalts (PUBs). Samples were analyzed from the Lucky Strike segment on the Mid-Atlantic Ridge and from three locations on the Juan de Fuca Ridge (West Valley, Endeavor Segment, and Axial Segment). Melt inclusions were re-homogenized to remove the effects of post-entrapment crystallization. The CO2 in the vapor bubbles present in the melt inclusions were analyzed at Virginia Tech using Raman spectroscopy, and associated glassy melt inclusions were analyzed at WHOI using the ion microprobe for CO2 and H2O. Vapor-saturation pressures calculated from these volatiles stored in melt inclusions and vapor bubbles range from 359-3994 bars, corresponding to depths of 1.0-11.4 km below the sea floor. The proportion of CO2 partitioned in the bubbles range from 11-98%. In summary, about 14% of the melt inclusions from Lucky Strike record crystallization depths of 3-4 km, consistent with the depth of the seismically imaged melt lens, whereas ~55% of melt inclusions crystallized at depths >4 km with a maximum at 9.8 km. These data are similar to depths of formation determined through olivine-hosted melt inclusions from the same segment (Wanless et al., 2015), although a greater portion of plagioclase-hosted melt inclusions record crystallization below the melt lens. At the Juan de Fuca ridge, ~24% of the melt inclusions record crystallization depths of 2-3 km, consistent with a seismically imaged mid-crustal magma chamber at the Endeavor Segment, while an additional ~62% crystallize at depths >3 km with a maximum at 11.4 km. This suggests that while crystallization can be focused within the melt lenses and magma chambers at these ridge localities, a significant and greater proportion of the megacrysts were sampled from the lower crust or upper mantle. 
    more » « less
  4. Abstract Mount Elbrus, Europe's tallest and largely glaciated volcano, is made of silicic lavas and is known for Holocene eruptions, but the size and state of its magma chamber remain poorly constrained. We report high spatial resolution U–Th–Pb zircon ages, co-registered with oxygen and hafnium isotopic values, span ~ 0.6 Ma in each lava, documenting magmatic initiation that forms the current edifice. The best-fit thermochemical modeling constrains magmatic fluxes at 1.2 km 3 /1000 year by hot (900 °C), initially zircon-undersaturated dacite into a vertically extensive magma body since ~ 0.6 Ma, whereas a volcanic episode with eruptible magma only extends over the past 0.2 Ma, matching the age of oldest lavas. Simulations explain the total magma volume of ~ 180 km 3 , temporally oscillating δ 18 O and εHf values, and a wide range of zircon age distributions in each sample. These data provide insights into the current state (~ 200 km 3 of melt in a vertically extensive system) and the potential for future activity of Elbrus calling for much-needed seismic imaging. Similar zircon records worldwide require continuous intrusive activity by magmatic accretion of silicic magmas generated at depths, and that zircon ages do not reflect eruption ages but predate them by ~ 10 3 to 10 5  years reflecting protracted dissolution–crystallization histories. 
    more » « less
  5. Abstract

    Cordilleran arcs are built by long periods of steady-state magmatism punctuated by transient high-flux magmatic episodes or flare-ups. Such flare-ups, manifested as periods of prodigious silicic volcanism and magmatism, result from geodynamic perturbations that cause elevated rates of magma addition to the crust. Questions remain, however, about how magmatic addition rates quantitatively compare between steady-state and flare-up modes of arc magmatism, and how long after the major geodynamic perturbation the flare-up begins. Here, we compute new estimates of erupted volumes over the last 35 Myr for the 22.5–29°S segment of the Central Andes based on a new volcanic geospatial database. These yield magmatic addition rates at least an order of magnitude higher during flare-up compared to steady-state conditions. A lag time of ~8–12 Myr between ocean ridge subduction (the major geodynamic perturbation in the Central Andean arc) and the onset of flare-up conditions is estimated.

    more » « less