skip to main content


Title: To sink, or not to sink: The thermal and density structure of the modern northern Andean arc constrained by xenolith petrology
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.  more » « less
Award ID(s):
2131643 1926421
NSF-PAR ID:
10450535
Author(s) / Creator(s):
; ; ; ; ; ;
Date Published:
Journal Name:
Geology
Volume:
51
Issue:
6
ISSN:
0091-7613
Page Range / eLocation ID:
586 to 590
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. null (Ed.)
    Continental arcs in Cordilleran orogenic systems display episodic changes in magma production rate, alternating between flare ups (70–90 km3 km􀀀 1 Myr􀀀 1) and lulls (< 20 km3 km􀀀 1 Myr􀀀 1) on timescales of tens of millions of years. Arc segments or individual magmatic suites may have even higher rates, up several 100 s of km3 km􀀀 1 Myr􀀀 1, during flare ups. These rates are largely determined by estimating volumes of arc crust, but do not reflect melt production from the mantle. The bulk of mantle-derived magmas are recycled back into the mantle by delamination of arc roots after differentiation in the deep crust. Mantle-derived melt production rates for continental arcs are estimated to be 140–215 km3 km􀀀 1 Myr􀀀 1 during flare ups and ≤ 15 km3 km􀀀 1 Myr􀀀 1 during lulls. Melt production rates averaged over multiple magmatic cycles are consistent with independent estimates for partial melting of the mantle wedge in subduction zones, however, the rates during flare ups and lulls are both anomalously high and anomalously low, respectively. The difference in mantle-derived melt production between flare ups and lulls is larger than predicted by petrologic and numerical models that explore the range of globally observed subduction parameters (e.g., convergence rate, height of the mantle wedge). This suggests that other processes are required to increase magmatism during flare ups and suppress magmatism during lulls. There are many viable explanations, but one possibility is that crystallized melts from the asthenospheric mantle wedge are temporarily stored in the deep lithosphere during lulls and then remobilized during flare ups. Basaltic melts may stall in the mantle lithosphere in inactive parts of the arc system, like the back-arc, refertilizing the mantle lithosphere and suppressing melt delivery to the lower crust. Subsequent landward arc migration (i.e., toward the interior of the continent) may encounter such refertilized mantle lithosphere magma source regions, contributing to magmatic activity during a flare up. A review of continental arcs globally suggests that flare ups commonly coincide with landward arc migration and that this migration may start tens of millions of years before the flare up occurs. The region of magmatic activity, or arc width, can also expand significantly during a flare up. Arc migration or expansion into different mantle source regions and across lithospheric and crustal boundaries can cause temporal shifts in the radiogenic isotopic composition of magmatism. In the absence of arc migration, temporal shifts are more muted. Isotopic studies of mantle xenoliths and exposures of deep arc crust suggest that that primary, mantle-derived magmas generated during flare ups reflect substantial contributions from the subcontinental mantle lithosphere. Arc migration may be caused by a variety of mechanisms, including slab anchoring or slab folding in the mantle transition zone that could generate changes in slab dip. Episodic slab shallowing is associated with many tectonic processes in Cordilleran orogenic systems, like alternations between shortening and extension in the upper plate. Studies of arc migration may help to link irregular magmatic production in continental arcs with geodynamic models for orogenic cyclicity. 
    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. Abstract Carbon flux metasomatism in the subduction environment is an important process, but it remains poorly understood. The paucity of exposed lower crust and upper mantle rocks in continental arcs renders xenoliths a major target for studying the slab-derived carbon cycle. This study of the carbonate phases in volcanic rocks from three drill cores in Ulleung Island, South Korea, sheds light on the interaction of carbon flux in the upper mantle and lower crust in a back-arc setting. The volcanic rocks from Ulleung Island range in composition from trachybasalt to trachyte and contain abundant euhedral pseudomorphic carbonate grains, ulvöspinel-hosted and biotite-hosted carbonate-silicate melt inclusions, and irregular carbonate globules. Integrated petrographic and geochemical studies of a variety of phenocrysts, carbonate phases, and carbonate-silicate inclusions in biotite and ulvöspinel indicate that recharging of carbon flux affected magma evolution. Carbon and oxygen isotopes of the pseudomorphic carbonate grains overlap with mantle values, indicating a carbonatite-like origin of the carbonate phases. The (MgO, FeO, CaO)-rich silicates in ulvöspinel-hosted silicate inclusions and pseudomorphic carbonate grains likely represent a primary melt, which formed from the partial melting of carbonated eclogite of the subducted slab within the mantle wedge beneath Ulleung Island. A petrogenetic model is proposed to illustrate that the crystal mush in the magma chamber was intruded by carbonate-rich liquids and caused alteration of cumulate crystals to generate the euhedral pseudomorphic carbonate grains. The extrusive magma captured those pseudomorphic grains and erupted to form the trachybasalt-trachyte units. The observed carbonate phases and their geochemical characteristics indicate that carbon flux metasomatism played a fundamental role in this back-arc magmatism. 
    more » « less
  4. Abstract

    When continental rifting is accompanied by localized magmatism under extensional stress, the breakup duration can be short and the continent/ocean transition sharp, as mantle melts are thought to be efficient at heating and weakening the lithosphere. This mode of rifting has been invoked for the Eastern North American Margin (ENAM) based on the existing geophysical data. Here, we present results from multichannel seismic profiles from the ENAM Community Seismic Experiment offshore North Carolina, U.S. Our survey area encompasses both the East Coast Magnetic Anomaly (ECMA) and the Blake Spur Magnetic Anomaly (BSMA), which lies ~200‐km farther seaward. Our prestack depth‐migrated seismic images reveal major changes in the structure of the igneous crust across the BSMA. Between the ECMA and BSMA, we image a proto‐oceanic domain of rough, faulted, and thin igneous crust. The roughness of this oceanic crust is similar to modern ultraslow spreading environments which involve the continued presence of a pre‐existing lithospheric lid. Seaward of the BSMA the basement is smooth, and the crust is relatively thick, which is typical for Jurassic oceanic crust. Across the BSMA, we image a step up in basement and crustal root, which we interpret to represent complete lithospheric breakup and a transition to steady‐state seafloor spreading in agreement with coincident refraction results. Our results would also indicate low extension rates in the final stages of rifting that may have influenced the thermal structure of the lithosphere and could explain the delay for continental breakup. All of these observations show that although continental rifting between eastern North America and northwest Africa was assisted by magmatic activity, it did not lead to rapid localization of extensional strain as previously thought.

     
    more » « less
  5. Abstract The southern Coast Mountain batholith was episodically active from Jurassic to Eocene time and experienced four distinct high magmatic flux events during that period. Similar episodicity has been recognized in arcs worldwide, yet the mechanism(s) driving such punctuated magmatic behavior are debated. This study uses zircon Hf and O isotopes, with whole-rock and mineral geochemistry, to track spatiotemporal changes in southern Coast Mountains batholith melt sources and to evaluate models of flare-up behavior and crust formation in Cordilleran arc systems. Zircon Hf isotope analysis yielded consistently primitive values, with all zircon grains recording initial εHf between +6 and +16. The majority (97%) of zircons analyzed yielded δ18O values between 4.2‰ and 6.5‰, and only five grains recorded values of up to 8.3‰. These isotopic results are interpreted to reflect magmatism dominated by mantle melting during all time periods and across all areas of the southern batholith, which argues against the periodic input of more melt-fertile crustal materials as the driver of episodic arc magmatism. They also indicate that limited crustal recycling is needed to produce the large volumes of continental crust generated in the batholith. Although the isotopic character of intrusions is relatively invariant through time, magmas emplaced during flare-ups record higher Sr/Y and La/Yb(N) and lower zircon Ti and Yb concentrations, which is consistent with melting in thickened crust with garnet present as a fractionating phase. Flare-ups are also temporally associated with periods when the southern Coast Mountains batholith both widens and advances inboard. We suggest that the landward shift of the arc into more fertile lithospheric mantle domains triggers voluminous magmatism and is accompanied by magmatic and/or tectonic thickening. Overall, these results demonstrate that the magmatic growth of Cordilleran arcs can be spatially and temporally complex without requiring variability in the contributions of crust and/or mantle to the batholith. 
    more » « less