Search for: All records

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. 

Abstract The thermal and compositional structure of arcs influence magmatic differentiation and lower-crustal foundering, two key processes impacting the evolution of the continental crust. Although many studies have proposed time scales of lithospheric recycling based on convective downwelling calculations, these models depend on the composition, density (ρ), and thermal structure of the lower crust and mantle, which are difficult to quantify in active continental arcs. Here, we constrained these properties for the Andean Northern Volcanic Zone using direct petrologic observations from a unique suite of lower-crust and mantle xenoliths from Mercaderes, Colombia. Chemical abrasion–isotope dilution–thermal ionization mass spectrometry (CA-ID-TIMS) U-Pb dates for zircons within the host tuff indicate the xenoliths erupted no earlier than 238 (±19) ka and thus capture a recent snapshot of the arc and subarc mantle. Equilibrium pressure-temperature (P-T) estimates for 81 xenoliths define three distinct thermal domains, interpreted as (1) a steep conductive geothermal gradient in the lower arc crust; (2) a convecting mantle wedge; and (3) cooled mantle in proximity to the subducting slab. Our results indicate the presence of an ~10–14-km-thick, high-density lithospheric root that is ~0.1 g/cm3 denser than the underlying mantle. Unlike records from exhumed paleoarcs, Rayleigh-Taylor instability calculations using our P-T-ρ constraints are unrealistically short for the northern Andes. We suggest the presence of partial melts in this hot arc root as a potential source of buoyancy preventing or significantly slowing down foundering. 

Stable isotopes provide deep insights into processes across a wide range of scales, from micron- to cosmic-size systems. Here, we review how continued advances in mass-spectrometry have enabled the analysis of ever-smaller samples and brought the field of heavy stable isotope geochemistry to its next frontier: the single-crystal scale. Accessing this record can be as enlightening as it is challenging. Drawing on novel systematics at different stages of development (from well-established to nascent), we discuss how the isotopes of heavy elements, such as magnesium, iron, zirconium, or uranium, can be used at the single-crystal and subcrystal scales to reconstruct magma thermal histories, crystal growth timescales, or, possibly, magma redox conditions. 

The isotopic variability of the elements in our planet and Solar System is the end result of a complex mixture of processes, including variable production of isotopes in stars, ingrowth of daughter nuclides due to decay of radioactive parents, and selective incorporation of isotopes into solids, liquids, or gases as a function of their mass and/or nuclear volume. Interpreting the isotopic imprints that planetary formation and evolution have left in the rock and mineral record requires not only precise and accurate measurements but also an understanding of the drivers behind isotopic variability. Here, we introduce fundamental concepts needed to “read” the isotopic code, with particular emphasis on heavy stable isotope systems. 

Andean uplift played a fundamental role in shaping South American climate and species distribution, but the relationship between the rise of the Andes, plant composition, and local climatic evolution is poorly known. We investigated the fossil record (pollen, leaves, and wood) from the Neogene of the Central Andean Plateau and documented the earliest evidence of a puna-like ecosystem in the Pliocene and a montane ecosystem without modern analogs in the Miocene. In contrast to regional climate model simulations, our climate inferences based on fossil data suggest wetter than modern precipitation conditions during the Pliocene, when the area was near modern elevations, and even wetter conditions during the Miocene, when the cordillera was around ~1700 meters above sea level. Our empirical data highlight the importance of the plant fossil record in studying past, present, and future climates and underscore the dynamic nature of high elevation ecosystems.