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: CO2 and H2O in Plagioclase-Hosted Melt Inclusions from Ocean-Ridge Lavas: An Indicator of Crystallization in the Lower Oceanic Crust
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
Award ID(s):
1634206
PAR ID:
10191142
Author(s) / Creator(s):
; ; ;
Date Published:
Journal Name:
American Geophysical Union, Fall Meeting 2019, abstract #V13C-0172
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract Recent multi‐channel seismic studies of fast spreading and hot‐spot influenced mid‐ocean ridges reveal magma bodies located beneath the mid‐crustal Axial Magma Lens (AML), embedded within the underlying crustal mush zone. We here present new seismic images from the Juan de Fuca Ridge that show reflections interpreted to be from vertically stacked magma lenses in a number of locations beneath this intermediate‐spreading ridge. The brightest reflections are beneath Northern Symmetric segment, from ∼46°42′‐52′N and Split Seamount, where a small magma body at local Moho depths is also detected, inferred to be a source reservoir for the stacked magma lenses in the crust above. The imaged magma bodies are sub‐horizontal, extend continuously for along‐axis lengths of ∼1–8 km, with the shallowest located at depths of ∼100–1,200 m below the AML, and are similar to sub‐AML bodies found at the East Pacific Rise. At both ridges, stacked sill‐like lenses are detected beneath only a small fraction of the ridge length examined and are inferred to mark local sites of higher melt flux and active replenishment from depth. The imaged magma lenses are focused in the upper part of the lower crust, which coincides with the most melt rich part of the crystal mush zone detected in other geophysical studies and where sub‐vertical fabrics are observed in geologic exposures of oceanic crust. We infer that the multi‐level magma accumulations are ephemeral and may result from porous flow and mush compaction, and that they can be tapped and drained during dike intrusion and eruption events. 
    more » « less
  2. 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
  3. 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
  4. 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
  5. We present DiadFit—an open-source Python3 tool for efficient processing of Raman spectroscopy data collected from fluid inclusions, melt inclusions and silicate melts. DiadFit is optimized to fit the characteristic peaks from CO2 fluids (Fermi diads, hot bands, 13C), gas species such as SO2, N2, solid precipitates (e.g. carbonates), and Ne emission lines with easily tweakable background positions and peak shapes. DiadFit's peak fitting functions are used as part of a number of workflows optimized for quantification of CO2 in melt inclusion vapour bubbles and fluid inclusions. DiadFit can also convert between temperature, pressure and density using various CO2 and CO2-H2O equations of state (EOS), allowing calculation of fluid inclusion pressures (and depths in the crust), conversion of homogenization temperatures from microthermometry to CO2 density, and propagation of uncertainties associated with EOS calculations using Monte Carlo methods. There are also functions to quantify the area ratio of the silicate vs. H2O region of spectra collected on silicate glasses to determine H2O contents in glasses and melt inclusions. Documentation and worked examples are available (https://bit.ly/DiadFitRTD, https://bit.ly/DiadFitYouTube). 
    more » « less