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.


This content will become publicly available on October 1, 2026

Title: Origins and fluxes of gas emissions from the Central Volcanic Zone of the Andes
We present geochemical data from gas samples from ~1200 km of arc in the Central Volcanic Zone of the Andes (CVZA), the volcanic arc with the thickest (~70 km) continental crust globally. The primary goals of this study are to characterize and understand how magmatic gases interact with hydrothermal systems, assess the origins of the major gas species, and constrain gas emission rates. To this end, we use gas chemistry, isotope compositions of H, O, He, C, and S, and SO2 fluxes from the CVZA. Gas and isotope ratios (CO2/ST, CO2/CH4, H2O/ST, δ13C, δ34S, 3He/4He) vary dramatically as magmatic gases are progressively affected by hydrothermal processes, reflecting removal and crustal sequestration of reactive species (e.g., S) and addition of less reactive meteoric and crustal components (e.g., He). The observed variations are similar in magnitude to those expected during the magmatic reactivation of volcanoes with hydrothermal systems. Carbon and sulfur isotope compositions of the highest temperature emissions (97–408 ◦C) are typical of arc magmatic gases. Helium isotope compositions reach values similar to upper mantle in some volcanic gases indicating that transcustal magma systems are effective conduits for volatiles, even through very thick continental crust. However, He isotopes are highly sensitive to even low degrees of hydrothermal interaction and radiogenic overprinting. Previous work has significantly underestimated volatile fluxes from the CVZA; however, emission rates from this study also appear to be lower than typical arcs, which may be related to crustal thickness.  more » « less
Award ID(s):
2151015
PAR ID:
10629128
Author(s) / Creator(s):
; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; more » ; ; « less
Publisher / Repository:
Journal of Volcanology and Geothermal Research
Date Published:
Journal Name:
Journal of Volcanology and Geothermal Research
Volume:
466
Issue:
C
ISSN:
0377-0273
Page Range / eLocation ID:
108382
Subject(s) / Keyword(s):
Volcanic gas geochemistry Stable isotope geochemistry Volatile fluxes Hydrothermal-magmatic transition
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Subduction zones represent the interface between Earth’s interior (crust and mantle) and exterior (atmosphere and oceans), where carbon and other volatile elements are actively cycled between Earth reservoirs by plate tectonics. Helium is a sensitive tracer of volatile sources and can be used to deconvolute mantle and crustal sources in arcs; however it is not thought to be recycled into the mantle by subduction processes. In contrast, carbon is readily recycled, mostly in the form of carbon-rich sediments, and can thus be used to understand volatile delivery via subduction. Further, carbon is chemically-reactive and isotope fractionation can be used to determine the main processes controlling volatile movements within arc systems. Here, we report helium isotope and abundance data for 42 deeply-sourced fluid and gas samples from the Central Volcanic Zone (CVZ) and Southern Volcanic Zone (SVZ) of the Andean Convergent Margin (ACM). Data are used to assess the influence of subduction parameters (e.g., crustal thickness, subduction inputs, and convergence rate) on the composition of volatiles in surface volcanic fluid and gas emissions. He isotopes from the CVZ backarc range from 0.1 to 2.6 R A ( n = 23), with the highest values in the Puna and the lowest in the Sub-Andean foreland fold-and-thrust belt. Atmosphere-corrected He isotopes from the SVZ range from 0.7 to 5.0 R A ( n = 19). Taken together, these data reveal a clear southeastward increase in 3 He/ 4 He, with the highest values (in the SVZ) falling below the nominal range associated with pure upper mantle helium (8 ± 1 R A ), approaching the mean He isotope value for arc gases of (5.4 ± 1.9 R A ). Notably, the lowest values are found in the CVZ, suggesting more significant crustal inputs (i.e., assimilation of 4 He) to the helium budget. The crustal thickness in the CVZ (up to 70 km) is significantly larger than in the SVZ, where it is just ∼40 km. We suggest that crustal thickness exerts a primary control on the extent of fluid-crust interaction, as helium and other volatiles rise through the upper plate in the ACM. We also report carbon isotopes from ( n = 11) sites in the CVZ, where δ 13 C varies between −15.3‰ and −1.2‰ [vs. Vienna Pee Dee Belemnite (VPDB)] and CO 2 / 3 He values that vary by over two orders of magnitude (6.9 × 10 8 –1.7 × 10 11 ). In the SVZ, carbon isotope ratios are also reported from ( n = 13) sites and vary between −17.2‰ and −4.1‰. CO 2 / 3 He values vary by over four orders of magnitude (4.7 × 10 7 –1.7 × 10 12 ). Low δ 13 C and CO 2 / 3 He values are consistent with CO 2 removal (e.g., calcite precipitation and gas dissolution) in shallow hydrothermal systems. Carbon isotope fractionation modeling suggests that calcite precipitation occurs at temperatures coincident with the upper temperature limit for life (122°C), suggesting that biology may play a role in C-He systematics of arc-related volcanic fluid and gas emissions. 
    more » « less
  2. Zirconium (Zr) stable isotope variations occur among co-existing Zr-rich accessory phases as well as at the bulk-rock scale, but the petrologic mechanism(s) responsible for Zr isotope fractionation during magmatic differentiation remain unclear. Juvenile magma generation and intra-crustal differentiation in convergent continental margins may play a crucial role in developing Zr isotope variations, and the Northern Volcanic Zone of the Andes is an ideal setting to test this hypothesis. To investigate the influence of these processes on Zr stable isotope compositions, we report δ94/90ZrNIST of whole rock samples from: 1) juvenile arc basalts from the Quaternary Granatifera Tuff, Colombia; 2) lower crust-derived garnet pyroxenites (i.e., arclogites), hornblendites, and gabbroic cumulates found in the same unit; and 3) felsic volcanic products from the Doña Juana Volcanic Complex, a dacitic composite volcano in close proximity to and partially covering the Granatifera Tuff. The basalts have δ94/90ZrNIST values ranging from −0.025 ± 0.018 ‰ to +0.003 ± 0.015 ‰ (n = 8), within the range of mid-ocean ridge basalts. The dacites have δ94/90ZrNIST values ranging from +0.008 ± 0.013 ‰ to +0.043 ± 0.015 ‰ (n = 14), slightly positive relative to the Granatifera and mid-ocean ridge basalts. In contrast, the (ultra)mafic cumulates have highly variable, predominantly positive δ94/90ZrNIST values, ranging from −0.134 ± 0.012 ‰ to +0.428 ± 0.012 ‰ (n = 15). Individual grains and mineral fractions of major rock-forming phases, including garnet (n = 21), amphibole (n = 9), and clinopyroxene (n = 18), were analyzed from 8 (ultra)mafic cumulates. The mineral fractions record highly variable Zr isotopic compositions, with inter-mineral fractionation (Δ94/90Zrgarnet-amphibole) up to 2.067 ‰. Recent ab initio calculations of Zr–O bond force constants in rock-forming phases predict limited inter-mineral Zr isotope fractionation in high-temperature environments, suggesting that the large fractionations we observe are not the product of vibrational equilibrium processes. Instead, we propose a scenario in which large Zr isotopic fractionations develop kinetically, induced by sub-solidus Zr diffusion between coexisting phases via changes in Zr distribution coefficients that arise from changes in temperature. Altogether, Zr isotope variability in this calc-alkaline continental arc setting exhibits no correlation with indices of magmatic differentiation (e.g., Mg#, SiO2), and is not a simple function of fractional crystallization. Furthermore, the garnet clinopyroxenite cumulates studied here represent density-unstable lower arc crust material; consequently, material with isotopically variable δ94/90Zr can be recycled into the mantle as a consequence of lower crustal foundering. 
    more » « less
  3. The Mineral King pendant in the Sierra Nevada batholith (California, USA) contains at least four rhyolite units that record high-silica volcanism during magmatic lulls in the Sierran magmatic arc. U-Th-Pb, trace element (single crystal spot analyses via sensitive high-resolution ion microprobe–reverse geometry, SHRIMP-RG), and bulk oxygen isotope analyses of zircon from these units provide a record of the age and compositional properties of the magmas that is not available from whole-rock analysis because of intense hydrothermal alteration of the pendant. U-Pb spot ages reveal that the Mineral King rhyolites are from two periods, the Early Jurassic (197 Ma) and the Early Cretaceous (134–136 Ma). These two rhyolite packages have zircons with distinct compositional trends for trace elements and δ18O; the Early Jurassic rhyolite shows less evidence of crustal influences on the rhyolites and the Early Cretaceous rhyolite shows evidence of increasing crustal influences and crystal recycling. These rhyolites capture evidence of magmatism during two periods of low magmatic flux in the Sierran Arc; however, they still show that magmas were derived from interactions of maturing continental crust, increasing from the Early to Late Jurassic. This finding likely reflects the transition of the North America margin from one of docking island arcs in the Early Jurassic to one of a more mature continental arc in the Early Cretaceous. This also shows the utility in examining zircon spot ages combined with trace element and bulk isotopic composition to unlock the petrogenetic history of altered volcanic rocks. 
    more » « less
  4. These data correspond to the article “Deep Nitrogen Fluxes and Sources Constrained by Arc Lava Phenocrysts” by Hudak et al. submitted to Geophysical Research Letters. Table S1 includes N-He-Ar data for FIs in phenocrysts from mafic are lavas and tephras. Table S2 contains the corrected N2/3He data used for volcanic arc N flux calculations and the arc-averaged mean N arc flux. Table S3 summarizes previous literature estimates of N fluxes and the data used for those calculations. Table S4 provides the N concentrations, He concentrations, N isotope compositions of the mantle, sediments, and altered oceanic crust, as well as sediment thicknesses. Finally, Table S5 gives information about the sources of the mineral separates used for these analyses. 
    more » « less
  5. null (Ed.)
    Volcanic arcs are the surface expression of magmatic systems that result from subduction of mostly oceanic lithosphere at convergent plate boundaries. Arcs with a submarine component include intraoceanic arcs and island arcs that span almost 22,000 km on Earth’s surface, and the vast majority of them are located in the Pacific region. Hydrothermal systems hosted by submarine arc volcanoes commonly contain a large component of magmatic fluid. This magmatic-hydrothermal signature, coupled with the shallow water depths of arc volcanoes and their high volatile contents, strongly influences the chemistry of the fluids and resulting mineralization and likely has important consequences for the biota associated with these systems. The high metal content and very acidic fluids in these hydrothermal systems are thought to be important analogs to numerous porphyry copper and epithermal gold deposits mined today on land. During International Ocean Discovery Program (IODP) Expedition 376 (5 May–5 July 2018), a series of five sites was drilled on Brothers volcano in the Kermadec arc. The expedition was designed to provide the missing link (i.e., the third dimension) in our understanding of hydrothermal activity and mineral deposit formation at submarine arc volcanoes and the relationship between the discharge of magmatic fluids and the deep biosphere. Brothers volcano hosts two active and distinct hydrothermal systems: one is seawater influenced and the other is affected by magmatic fluids (largely gases). In total, 222.4 m of volcaniclastics and lavas were recovered from the five sites drilled, which include Sites U1527 and U1530 in the Northwest (NW) Caldera seawater-influenced hydrothermal field; Sites U1528 and U1531 in the magmatic fluid-influenced hydrothermal fields of the Upper and Lower Cones, respectively; and Site U1529, located within an area of low crustal magnetization that marks the West (W) Caldera upflow zone on the caldera floor. Downhole logging and borehole fluid sampling were completed at two sites, and two tests of a prototype turbine-driven coring system (designed by the Center for Deep Earth Exploration [CDEX] at Japan Agency for Marine-Earth Science and Technology [JAMSTEC]) for drilling and coring hard rocks were conducted. Core recovered from all five sites consists of dacitic volcaniclastics and lava flows with only limited chemical variability relative to the overall range in composition of dacites in the Kermadec arc. Pervasive alteration with complex and variable mineral assemblages attest to a highly dynamic hydrothermal system. The upper parts of several drill holes at the NW Caldera hydrothermal field are characterized by secondary mineral assemblages of goethite + opal + zeolites that result from low-temperature (<150°C) reaction of rock with seawater. At depth, NW Caldera Site U1527 exhibits a higher temperature (~250°C) secondary mineral assemblage dominated by chlorite + quartz + illite + pyrite. An older mineral assemblage dominated by diaspore + quartz + pyrophyllite + rutile at the bottom of Hole U1530A is indicative of acidic fluids with temperatures of ~230°–320°C. In contrast, the alteration assemblage at Site U1528 on the Upper Cone is dominated by illite + natroalunite + pyrophyllite + quartz + opal + pyrite, which attests to high-temperature reaction of rocks with acid-sulfate fluids derived from degassed magmatic volatiles and the disproportionation of magmatic SO2. These intensely altered rocks exhibit extreme depletion of major cation oxides, such as MgO, K2O, CaO, MnO, and Na2O. Furthermore, very acidic (as low as pH 1.8), relatively hot (≤236°C) fluids collected at 160, 279, and 313 meters below seafloor in Hole U1528D have chemical compositions indicative of magmatic gas input. In addition, preliminary fluid inclusion data provide evidence for involvement of two distinct fluids: phase-separated (modified) seawater and a ~360°C hypersaline brine, which alters the volcanic rock and potentially transports metals in the system. The material and data recovered during Expedition 376 provide new stratigraphic, lithologic, and geochemical constraints on the development and evolution of Brothers volcano and its hydrothermal systems. Insights into the consequences of the different types of fluid–rock reactions for the microbiological ecosystem elucidated by drilling at Brothers volcano await shore-based studies. 
    more » « less