skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: Thermal Budgets of Magma Storage Constrained by Diffusion Chronometry: the Cerro Galán Ignimbrite
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
Award ID(s):
1948862 1949173
PAR ID:
10428481
Author(s) / Creator(s):
; ;
Date Published:
Journal Name:
Journal of Petrology
Volume:
63
Issue:
7
ISSN:
0022-3530
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. 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
  2. Abstract We present a model for a coupled magma chamber–dike system to investigate the conditions required to initiate volcanic eruptions and to determine what controls the size of eruptions. The model combines the mechanics of dike propagation with internal chamber dynamics including crystallization, volatile exsolution, and the elastic response of the magma and surrounding crust to pressure changes within the chamber. We find three regimes for dike growth and eruptions: (1) below a critical magma chamber size, eruptions are suppressed because chamber pressure drops to lithostatic before a dike reaches the surface; (2) at an intermediate chamber size, the erupted volume is less than the dike volume (“dike-limited” eruption regime); and (3) above a certain chamber size, dikes can easily reach the surface and the erupted volume follows a classic scaling law, which depends on the attributes of the magma chamber (“chamber-limited” eruption regime). The critical chamber volume for an eruption ranges from ∼0.01 km3 to 10 km3 depending on the water content in the magma, depth of the chamber, and initial overpressure. This implies that the first eruptions at a volcano likely are preceded by a protracted history of magma chamber growth at depth, and that the crust above the magma chamber may have trapped several intrusions or “failed eruptions.” Model results can be combined with field observations of erupted volume, pressure, and crystal and volatile content to provide tighter constraints on parameters such as the eruptible chamber size. 
    more » « less
  3. Mocho-Choshuenco volcano (39.9°S, 72.0°W) produced ∼75 explosive eruptions following retreat of the >1.5-km-thick Patagonian Ice Sheet associated with the local Last Glacial Maximum (LGM, from 35 to 18 ka). Here, we extend this record of volcanic evolution to include pre- and syn-LGM lavas that erupted during the Pleistocene. We establish a long-term chronology of magmatic and volcanic evolution and evaluate the relationship between volcanism and loading/unloading of the Patagonian Ice Sheet via twenty-four 40Ar/39Ar and two 3He age determinations integrated with stratigraphy and whole-rock compositions of lava flows and glass compositions of tephra. Our findings reveal that the edifice is much younger than previously thought and preserves 106 km3 of eruptive products, of which 50% were emplaced immediately following the end of the penultimate glaciation and 20% after the end of the LGM. A period of volcanic inactivity between 37 and 26 ka, when glaciers expanded, was followed by the eruption of incompatible element-rich basaltic andesites. Several of these syn-LGM lavas dated between 26 and 16 ka, which crop out at 1500−1700 m above sea level, show ice contact features that are consistent with emplacement against a 1400- to 1600-m-thick Patagonian Ice Sheet. Small volume dacitic eruptions and two explosive rhyolitic eruptions dominate the volcanic output from 18 to 8 ka, when the Patagonian Ice Sheet began to retreat rapidly. We hypothesize that increased lithostatic loading as the Patagonian Ice Sheet grew prohibited dike propagation, thus stalling the ascent of magma, promoting growth of at least three discrete magma reservoirs, and enhancing minor crustal assimilation to generate incompatible element-rich basaltic andesitic to dacitic magmas that erupted between 26 and 17 ka. From an adjacent reservoir, incompatible element-poor dacites erupted from 17 to 12 ka. These lava flows were followed by the caldera-forming eruption at 11.5 ka of 5.3 km3 of rhyolite from a deeper reservoir atop which a silicic melt lens had formed and expanded. Subsequent eruptions of oxidized dacitic magmas from the Choshuenco cone from 11.5 to 8 ka were followed by andesitic to dacitic eruptions at the more southerly Mocho cone, as well as small flank vent eruptions of basaltic andesite at 2.5 and 0.5 ka. This complex history reflects a multi-reservoir plumbing system beneath Mocho-Choshuenco, which is characterized by depths of magma storage, oxidation states, and trace element compositions that vary over short periods of time (<2 k.y.). 
    more » « less
  4. Abstract Cerro Machín, a volcano located in the northern segment of the Andes, is considered one of the most dangerous volcanoes in Colombia with an explosive record that involves at least five plinian events. Prior studies focused on the last dome-building eruption have suggested the presence of a water-rich mid-crustal magma reservoir. However, no direct volatile measurements have been published and little work has been completed on the explosive products of the volcano. Here, we study the largest known eruption of Cerro Machín volcano which occurred 3600 years BP producing dacitic pyroclastic fall deposits that can be traced up to 40 km from the vent. Lapilli pumice clasts have a mineral assemblage of plagioclase, amphibole, quartz, and biotite phenocrysts, with accessory olivine, Fe–Ti oxides, and apatite. The occurrence of Fo89–92 olivine rimmed by high Mg# amphibole and the established high-water contents in the magma imply the presence of magma near or at water saturation at pressures > ~ 500 MPa. Measurements of up to 10.7 wt % H2O in melt inclusions hosted in plagioclase and quartz in the 3600 years BP eruption products support the idea that Cerro Machín is a remarkably water-rich volcanic system. Moreover, this is supported by measurements of ~103 to 161 ppm H2O in plagioclase phenocrysts. The application of two parameterizations of water partitioning between plagioclase and silicate melt allows us to use our water in plagioclase measurements to estimate equilibrium melt water contents of 5 ± 1 wt % to 11 ± 2 wt % H2O, which are in good agreement with the water contents we measured in melt inclusions. Results of amphibole geobarometry are consistent with a magma reservoir stored in the mid-to-lower crust at a modal pressure of 700 ± 250 MPa, corresponding to a depth of ~25 km. Minor crystallization in the shallow crust is also recorded by amphibole barometry and calculated entrapment pressures in melt inclusions. Amphibole is present as unzoned and zoned crystals. Two populations of unzoned amphibole crystals are present, the most abundant indicate crystallization conditions of 853 ± 26°C (1 se; standard error), and the less abundant crystallized at an average temperature of 944 ± 24°C (1 se). Approximately 18% of the amphibole crystals are normally or reversely zoned, providing evidence for a minor recharge event that could have been the trigger mechanism for the explosive eruption. Plagioclase crystals also show normal and reverse zoning. The moderate Ni concentrations (<1600 μg/g) in the high-Fo olivine xenocrysts suggest that Cerro Machín primary magmas are generated by inefficient interaction of mantle peridotite with a high-silica melt produced by slab melting of basaltic material. Some sediment input is also suggested by the high Pb/Th (>2.2) and Th/La (0.3–0.4) ratios. Whole rock chemistry reveals heavy rare earth element (HREE) depletion and Sr enrichment that likely formed during the crystallization of garnet and amphibole in the upper part of the mantle or lower portion of the crust, promoting the formation of water-rich dacitic magma that was then injected into the middle-to-lower crust. Textural and compositional differences in the crystal cargo that erupted during dome-building and plinian events support the idea that large volumes of magma recharge lead to effusive eruptions, while only small recharge events are needed to trigger plinian eruptions at Cerro Machín. 
    more » « less
  5. Abstract Magmatic systems below volcanoes are often dominated by partially crystalline magma over the long term. Rejuvenation of these systems during eruptive events can impact lava composition and eruption style—sometimes resulting in more violent or explosive activity than would be expected, as was the case at Fissure 17 during Kīlauea’s 2018 eruption. Here, we explore how the crystallinity of unerupted intrusion magmas affect hybrid magma compositions and petrological signatures by constructing phase-equilibria models to evaluate mineral and melt compositions of low-MgO lavas erupted along the East Rift Zone of Kīlauea volcano on 30 to 31 January 1997 (Episode 54, Fissures A-F). We then compare calculated mixing proportions and petrologically derived magma volumes to GPS-based geodetic inversions of ground deformation and intrusion growth in an attempt to reconcile geodetic and petrologically estimated magma volumes. Open-system phase-equilibria thermodynamic models were used to constrain the composition, degree of differentiation, and thermodynamic state of a rift-stored, two pyroxene + plagioclase saturated low-MgO magma body immediately preceding its mixing with high-MgO recharge and degassed drainback (lava lake) magma from Pu‘u‘ō‘ō‘, shortly before fissure activity within Nāpau Crater began on 29 January 1997. Mixing models constructed using the Magma Chamber Simulator reproduce the mineralogy and compositions of Episode 54 lavas within uncertainties and suggest that the identity of the low-MgO magma body may be either variably differentiated remnants of un-erupted magmas intruded into Nāpau Crater in October 1968, or another spatially and compositionally similar magma body. We find that magmas derived from a single, compositionally stratified magma emplaced beneath Nāpau Crater in 1968 can mix with mafic Kīlauea magmas to reproduce average Episode 54 bulk lava, mineralogy and mineral compositions without necessitating the interaction of multiple, low-MgO rift-stored magma bodies to produce Episode 54 lava compositions. Further, by constructing phase equilibria-based mixing models of Episode 54, we can better define the pre-eruptive state of the magmatic system. The resultant mineral assemblages and compositions are consistent with the possibility that the now-fractionated, rift-stored magma body was compositionally stratified and ~ 40% to 50% crystalline at the time of mixing. Finally, we estimate the volume of the low-MgO magma body to be ~7.51 Mm3. Phase-equilibria model results corroborate field and geochemical relationships demonstrating how shallow intrusions at intraplate shield volcanoes can crystallize, evolve, and then be remobilized by new, later batches of mafic magma. Most notably, our MCS models demonstrate that the pre-eruptive conditions of an intrusive body may be recovered by examining mineral compositions within mixed lavas. Discrepancies between the geodetic constraints on volumes of stored rift versus newly intruded (recharge) magma and our best-fit results produced by MCS mixing models (which respectively are mmafic:mlow-MgO ≈ 2 vs. mmafic:mlow-MgO ≈ 0.75) are interpreted to highlight the complex nature of incomplete mixing on more localized scales as reflected in erupted lavas, compared to geodetically constrained volumes that likely reflect large spatial scale contributions to a magmatic system. These dissimilar volume relationships may also help to constrain eruptive versus unerupted volumes in magmatic systems undergoing mixing. By demonstrating the usefulness of MCS in modeling past eruptions, we highlight the potential to use it as a tool to aid in petrologic monitoring of ongoing activity. 
    more » « less