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
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:
- 10544681
- 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
-
-
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
-
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
-
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
-
Abstract Olivine‐hosted melt inclusions are an important archive of pre‐eruptive processes such as magma storage, mixing and subsequent ascent through the crust. However, this record can be modified by post‐entrapment diffusion of H+through the olivine lattice. Existing studies often use spherical or 1D models to track melt inclusion dehydration that fail to account for complexities in geometry, diffusive anisotropy and sectioning effects. Here we develop a finite element 3D multiphase diffusion model for the dehydration of olivine‐hosted melt inclusions that includes natural crystal geometries and multiple melt inclusions. We use our 3D model to test the reliability of simplified analytical and numerical models (1D and 2D) using magma ascent conditions from the 1977 eruption of Seguam volcano, Alaska. We find that 1D models underestimate melt inclusion water loss, typically by ∼30%, and thus underestimate magma decompression rates, by up to a factor of 5, when compared to the 3D models. An anisotropic analytical solution that we present performs well and recovers decompression rates within a factor of 2, in the situations in which it is valid. 3D models that include multiple melt inclusions show that inclusions can shield each other and reduce the amount of water loss upon ascent. This shielding effect depends on decompression rate, melt inclusion size, and crystallographic direction. Our modeling approach shows that factors such as 3D crystal geometry and melt inclusion configuration can play an important role in constraining accurate decompression rates and recovering water contents in natural magmatic systems.more » « less
An official website of the United States government
