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: Crystals and Melt Inclusions Record Deep Storage of Superhydrous Magma Prior to the Largest Known Eruption of Cerro Machín Volcano, Colombia
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
Award ID(s):
2237994
PAR ID:
10571554
Author(s) / Creator(s):
; ; ;
Publisher / Repository:
Oxford University Press
Date Published:
Journal Name:
Journal of Petrology
Volume:
65
Issue:
9
ISSN:
0022-3530
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. The origin of gaps or zoning in the composition of erupted products is critical to understanding how sub-volcanic reservoirs operate. We characterize the compositionally zoned magma that produced the 2053 ± 50 cal. yr BP Paso Puyehue Tephra from the Antillanca Volcanic Complex in the Andean Southern Volcanic Zone (SVZ). The 3.7 km3 Paso Puyehue Tephra is zoned from dacite (69 wt% SiO2) lapilli and ash comprising the lowermost 80% of the deposit that abruptly transitions upward into basaltic andesite scoria (54 wt% SiO2) making up the remaining ~20%. Variations in whole-rock, matrix glass, and mineral compositions through the deposit allow us to estimate pre-eruptive magma storage conditions and to develop a model of how this magma body was generated. Our findings suggest that amphibole-bearing basaltic andesitic magma stored at ~8.0 ± 1.3 km depth fractionally crystallized and cooled from 1048 ± 1.1 to 811 ± 28.6 ◦C under highly oxidizing conditions to produce silicic a melt that upon extraction and rise, pooled at ~6.4 ± 1.2 km depth at temperatures as low as 810 ◦C before eruption. MELTS models suggest that crystallization of a basaltic andesite parent magma with 4 wt% dissolved H2O can produce the dacite under conditions predicted by mineral thermobarometers with phase compositions comparable to those measured in minerals. Pervasive normal zoning at the rims of plagioclase crystals—most pronounced at the transition between dacite and basaltic andesite, and compatible vs. incompatible trace element concentrations, suggest that magma mixing was limited and likely occurred at the interface between the dacitic and basaltic andesitic magmas during ascent within the conduit upon eruption. Compositionally bimodal tephras are increasingly recognized throughout the SVZ with several interpreted to reflect basaltic recharge and mixing into extant rhyolitic reservoirs. In further contrast to other SVZ rhyolitic products, e.g., from the nearby Cord´on Callue and Mocho Choshuenco volcanoes, the Paso Puyehue magma was highly oxidized. This may reflect enhanced delivery of H2O from the subducting plate into the mantle wedge, which in turn may facilitate efficient extraction and separation of buoyant, low-viscosity rhyolitic melt from crystal-rich basaltic andesitic parent magmas and the co-eruption of both end members. 
    more » « less
  2. The 2021 La Palma eruption (Tajogaite) was unprecedented in magnitude, duration, and degree of monitoring compared to historical volcanism on La Palma. Here, we provide data on melt inclusions in samples from the beginning and end of the eruption to compare the utility of both melt and fluid inclusions as recorders of magma storage. We also investigated compositional heterogeneities within the magmatic plumbing system. We found two populations of olivine crystals: a low Mg# (78–82) population present at the beginning and end of eruption, recording the maximum volatile contents (2.5 wt % H2O, 1,800 ppm F, 700 ppm Cl, 3,800 ppm S) and a higher Mg# (83–86) population sampled toward the end of the eruption, with lower volatile contents. Despite their host composition, melt inclusions share the same maximum range of CO2 concentrations (1.2–1.4 wt %), indicating olivine growth and inclusion capture at similar depths. Overall, both melt and fluid inclusions record similar pressures (450–850 MPa, ∼15–30 km), and when hosted in the same olivine crystal pressures are indistinguishable within error. At these mantle pressures, CO2 is expected to be an exsolved phase explaining the similar range of CO2 between the two samples, but other volatile species (F, Cl, S) behave incompatibly, and thus, the increase between the two olivine populations can be explained by fractional crystallization prior to eruption. Finally, based on our new data, we provide estimates on the total volatile emission of the eruption. 
    more » « less
  3. Inclusions of basaltic melt trapped inside of olivine phenocrysts during igneous crystallization provide a rich, crystal-scale record of magmatic processes ranging from mantle melting to ascent, eruption, and quenching of magma during volcanic eruptions. Melt inclusions are particularly valuable for retaining information on volatiles such as H 2 O and CO 2 that are normally lost by vesiculation and degassing as magma ascends and erupts. However, the record preserved in melt inclusions can be variably obscured by postentrapment processes, and thus melt inclusion research requires careful evaluation of the effects of such processes. Here we review processes by which melt inclusions are trapped and modified after trapping, describe new opportunities for studying the rates of magmatic and volcanic processes over a range of timescales using the kinetics of post-trapping processes, and describe recent developments in the use of volatile contents of melt inclusions to improve our understanding of how volcanoes work. ▪  Inclusions of silicate melt (magma) trapped inside of crystals formed by magma crystallization provide a rich, detailed record of what happens beneath volcanoes. ▪  These inclusions record information ranging from how magma forms deep inside Earth to its final hours as it ascends to the surface and erupts. ▪  The melt inclusion record, however, is complex and hazy because of many processes that modify the inclusions after they become trapped in crystals. ▪  Melt inclusions provide a primary archive of dissolved gases in magma, which are the key ingredients that make volcanoes erupt explosively. 
    more » « less
  4. 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
  5. Abstract Whole rock compositions at Buldir Volcano, western Aleutian arc, record a strong, continuous trend of iron depletion with decreasing MgO, classically interpreted as a calc-alkaline liquid line of descent. In contrast, olivine-hosted melt inclusions have higher total iron (FeO*) than whole rocks and show little change in FeO* with decreasing MgO. To investigate this discrepancy and determine the conditions required for strong iron depletion, we conducted oxygen fugacity (ƒO2) buffered, water-saturated crystallization experiments at 800 MPa and ƒO2 = QFM + 1.6 ± 0.4 (1$$\sigma$$) (where QFM refers to the quartz-fayalite-magnetite buffer) on a high-Al, basaltic starting material modeled after a Buldir lava. Experimental conditions were informed by olivine-hosted melt inclusions that record minimum entrapment pressures as high as 570 MPa, >6 wt % H2O, and ƒO2 of QFM + 1.4 (±0.2), making Buldir one of the most oxidized and wettest arc volcanoes documented globally. The experiments produce melts with Si-enrichment and Fe-depletion signatures characteristic of evolved, calc-alkaline magmas at the lowest MgO, although FeO* remains roughly constant over most of the experimental temperature range. Experiments saturate CrAl-spinel and olivine at 1160°C, followed by clinopyroxene and Al-spinel at 1085°C, hornblende at 1060°C, and, finally, plagioclase and magnetite between 1040°C and 960°C. Hornblende crystallization, not magnetite, generates the largest increase in SiO2 and largest decrease in FeO* in coexisting melts. Compositions of melt inclusions are consistent with experimental melts and reflect crystallization of a basaltic parent magma at high PH2O. In contrast, the whole rock compositional trends are influenced by magma mixing and phenocryst redistribution and accumulation. The crystallization experiments and natural liquids (melt inclusions and groundmass glass) from Buldir suggest that for an oxidized, hydrous primary basalt starting composition, significant Fe depletion from the melt will not occur until intermediate to late stages of magma crystallization (< ~4.5 wt % MgO). We conclude that the Buldir whole rock trend cannot be reproduced by crystallization at arc-relevant oxygen fugacities and is not a true liquid line of descent, warranting caution when interpreting volcanic trends globally. 
    more » « less