skip to main content


Title: The formation of high-Sr/Y plutons in cordilleran-arc crust by crystal accumulation and melt loss
Bulk-rock data are commonly used in geochemical studies as a proxy for melt compositions in order to understand the evolution of crustal melts. However, processes of crystal accumulation and melt migration out of deep-crustal, crystal-rich mush zones to shallower storage regions raise questions about how faithfully bulk-rock compositions in plutons approximate melt compositions. This problem is particularly acute in the lower crust of arcs, where melt reservoirs are subject to periodic melt extraction that leaves behind a cumulate residue. Here, we examine bulk-rock data from the perspective of high-Sr/Y plutonic rocks in the lower crust of a well-exposed Early Cretaceous cordilleran-arc system in Fiordland, New Zealand. We test the validity of using high-Sr/Y bulk-rock compositions as proxies for melts by comparing bulk-rock compositions to melts modeled from >100 major- and trace-element analyses of 23 magmatic clinopyroxene grains from the same samples. The sampling locations of the igneous clinopyroxenes and encompassing bulk rocks are distributed across ~550 km2 of exhumed lower crust and are representative of Mesozoic lower-crustal arc rocks in the Median batholith. We confirm that bulk-rock data have characteristics of high-Sr/Y plutons (Sr/Y >50, Na2O >3.5 wt%, Sr >1000 ppm, and Y <20 ppm), features that have been previously interpreted to indicate the presence of garnet as a residual or fractionating phase. In contrast to bulk rocks, igneous clinopyroxenes have low Sr (<100 ppm), high Y (25–100 ppm), and low molar Mg# [100 × Mg/(Mg + Fe)] values (60–70), which are consistent with derivation from fractionated, low-Sr/Y melts. Chondrite-normalized rare-earth-element patterns and Sm/Yb values in clinopyroxenes also show little to no evidence for involvement of garnet in the source or in differentiation processes. Fe-Mg partitioning relationships indicate that clinopyroxenes are not in equilibrium with their encompassing bulk rocks but could have been in equilibrium with melt compositions determined from chemometry of coexisting igneous hornblendes. Moho-depth calculations based on bulk-rock Sr/Y values also yield Moho depths (average = 69 km) that are inconsistent with Moho depths based on bulk-rock Ce/Y, contact aureole studies, Al-in-horn- blende crystallization pressures, and our modeled clinopyroxene crystallization pressures. These data indicate that most Mesozoic high-Sr/Y bulk rocks in the lower crust of Fiordland are cumulates formed by plagioclase + amphibole + clinopyroxene accumulation and interstitial melt loss from crystal-rich mush zones. Our data do not support widespread fractionation of igneous garnet nor partial melting of a garnet-bearing source in the petrogenesis of these melts. We speculate that melt extraction and the production of voluminous cumulates in the lower crust were aided by unusually high heat flow and high magma addition rates associated with an Early Cretaceous arc flareup. We conclude that bulk-rock compositions are poor proxies for melt compositions in the lower crust of the Median batholith, and geochemical modeling of these high-Sr/Y bulk rocks would overemphasize the role of garnet in their petrogenesis.  more » « less
Award ID(s):
1901827
NSF-PAR ID:
10351519
Author(s) / Creator(s):
;
Date Published:
Journal Name:
Geosphere
ISSN:
1553-040X
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract We present microbeam major- and trace-element data from 14 monzodiorites collected from the Malaspina Pluton (Fiordland, New Zealand) with the goal of evaluating processes involved in the production of andesites in lower arc crust. We focus on relict igneous assemblages consisting of plagioclase and amphibole with lesser amounts of clinopyroxene, orthopyroxene, biotite and quartz. These relict igneous assemblages are heterogeneously preserved in the lower crust within sheeted intrusions that display hypersolidus fabrics defined by alignment of unstrained plagioclase and amphibole. Trace-element data from relict igneous amphiboles in these rocks reveal two distinct groups: one relatively enriched in high field strength element concentrations and one relatively depleted. The enriched amphibole group has Zr values in the range of ∼25–110 ppm, Nb values of ∼5–32 ppm, and Th values up to 2·4 ppm. The depleted group, in contrast, shows Zr values <35 ppm and Nb values <0·25 ppm, and Th is generally below the level of detection. Amphibole crystallization temperatures calculated from major elements range from ∼960 to 830 °C for all samples in the pluton; however, we do not observe significant differences in the range of crystallization temperatures between enriched (∼960–840 °C) and depleted groups (∼940–830 °C). Bulk-rock Sr and Nd isotopes are also remarkably homogeneous and show no apparent difference between enriched (εNdi = 0·1 to –0·1; 87Sr/86Sri = 0·70420–0·70413) and depleted groups (εNdi = 0·3 to –0·4; 87Sr/86Sri = 0·70424–0·70411). Calculated amphibole-equilibrium melt compositions using chemometric equations indicate that melts were highly fractionated (molar Mg# <50), andesitic to dacitic in composition, and were much more evolved than bulk lower continental crust or primitive basalts and andesites predicted to have formed from hydrous melting of mantle-wedge peridotite beneath an arc. We suggest that melts originated from a common, isotopically homogeneous source beneath the Malaspina Pluton, and differences between enriched and depleted trace-element groups reflect varying contributions from subducted sediment-derived melt and sediment-derived fluid, respectively. Our data demonstrate that andesites and dacites were the dominant melts that intruded the lower crust, and their compositions mirror middle and upper bulk-continental crust estimates. Continental crust-like geochemical signatures were acquired in the source region from interaction between hydrous mantle-wedge melts and recycled subducted sediment rather than assimilation and/or remelting of pre-existing lower continental crust. 
    more » « less
  2. Abstract

    Water is key to plate tectonics on Earth, which, in turn, is vital to the production of continental crust. Although arc lavas erupt in a volatile‐rich state and calc‐alkaline arc plutons are distinguished by the presence of hydrous minerals such as hornblende and biotite, the water content of arc magmas earlier in their evolution—in the deep crust—remains poorly constrained. Here, we report H2O contents in nominally anhydrous minerals measured in situ on petrographic thin sections by secondary ion mass spectrometry of Proterozoic deep crustal xenoliths from Colorado, USA. Clinopyroxene, orthopyroxene, and garnet contain average H2O contents ranging from 75–760, 233–410, and 42–139 ppm, respectively. Reconstructed bulk rock H2O contents range from ~60 to ~650 ppm. Intermineral H2O ratios overlap experimental mineral/meltDvalues and are used to calculate H2O of melts last in equilibrium with the xenoliths. We propose that these xenoliths represent cumulates fractionated from a primitive, hydrous (≥1 wt.% H2O) melt at high (~1 GPa) pressures, similar to conditions in modern subduction zones and potentially associated with widespread arc accretion that formed the core of North America in the Precambrian.

     
    more » « less
  3. 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
  4. Interpretation of geochronological and petrological data from partially-melted granulite is challenging. However, integration of multiple chronometers and mineral assemblage diagrams (MAD) can be used to estimate the nature and duration of processes. Excellent lower-crustal exposures of garnet granulite from the Malaspina Pluton, Fiordland New Zealand provide an ideal place to employ this kitchen sink approach. We use zircon U-Pb ages from LA-ICPMS, SHRIMP-RG, and CA-TIMS, garnet Lu-Hf and Sm-Nd ages, and MAD in order to evaluate local partial melting vs. melt injection, equilibrium volumes, P-T conditions, and the duration of lower crustal thermal events. Host diorite (H), garnet-clinopyroxene reaction zones (GRZ), coarse garnet selvages, and tonalite veins provide a record of intrusion and granulite facies partial melting. Zircon U-Pb ages range from 123 to 107 Ma (all); LA-ICPMS ages contain the entire range; CA-TIMS ages range from 118.30±0.13 to 115.7±0.18 Ma; and SHRIMP-RG ages range from 121.4±2 to 109.8±1.8 Ma. The latter two techniques are interpreted to indicate primary igneous crystallization from ~119 to ~116 Ma and the youngest ~110 Ma ages are interpreted as metamorphic zircon growth. Garnet ages for ~1 cm grains are ~113 Ma (Lu-Hf & Sm-Nd) and record metamorphic growth, and <0.3 mm grains with Sm-Nd ages from 113 to 104 Ma reflect high temperature intracrystalline diffusion and isotopic closure during cooling to amphibolite facies. Zircon trace-element compositions indicate 2 distinct crystallization trends reflecting evolution of primary magma batches. MAD indicate that garnet was not in equilibrium with sampled rock compositions. Instead, garnet shows apparent equilibrium with a modeled mixture of the GRZ and the H and grew in equilibrium with an effective bulk composition that shifted toward the leucosome. This would produce the observed increase in garnet grossular content. We conclude that: Malaspina rocks from Crooked Arm preserve evidence for 2 igneous layers which evolved as discrete magmas, igneous crystallization lasted 2 to 3 m.y., granulite metamorphism peaked ~ 3 m.y. after intrusion, metamorphism lasted ≥3 m.y., cooling occurred at ~20°C/m.y., and granulite minerals equilibrated with a mixture of solid phases and melt at ~14 kbar and 920°C (based on garnet compositions and MAD). 
    more » « less
  5. null (Ed.)
    Abstract The 119 Ma Dinkey Dome pluton in the central Sierra Nevada Batholith is a peraluminous granite and contains magmatic garnet and zircon that are complexly zoned with respect to oxygen isotope ratios. Intracrystalline SIMS analysis tests the relative importance of magmatic differentiation processes vs. partial melting of metasedimentary rocks. Whereas δ18O values of bulk zircon concentrates are uniform across the entire pluton (7.7‰ VSMOW), zircon crystals are zoned in δ18O by up to 1.8‰, and when compared to late garnet, show evidence of changing magma chemistry during multiple interactions of the magma with wall rock during crustal transit. The evolution from an early high-δ18O magma [δ18O(WR) = 9.8‰] toward lower values is shown by high-δ18O zircon cores (7.8‰) and lower δ 18O rims (6.8‰). Garnets from the northwest side of the pluton show a final increase in δ18O with rims reaching 8.1‰. In situ REE measurements show zircon is magmatic and grew before garnets. Additionally, δ18O in garnets from the western side of the pluton are consistently higher (avg = 7.3‰) relative to the west (avg = 5.9‰). These δ18O variations in zircon and garnet record different stages of assimilation and fractional crystallization whereby an initially high-δ18O magma partially melted low-δ18O wallrock and was subsequently contaminated near the current level of emplacement by higher δ18O melts. Collectively, the comparison of δ18O zoning in garnet and zircon shows how a peraluminous pluton can be constructed from multiple batches of variably contaminated melts, especially in early stages of arc magmatism where magmas encounter significant heterogeneity of wall-rock assemblages. Collectively, peraluminous magmas in the Sierran arc are limited to small <100 km2 plutons that are intimately associated with metasedimentary wall rocks and often surrounded by later and larger metaluminous tonalite and granodiorite plutons. The general associations suggest that early-stage arc magmas sample crustal heterogeneities in small melt batches, but that with progressive invigoration of the arc, such compositions are more effectively blended with mantle melts in source regions. Thus, peraluminous magmas provide important details of the nascent Sierran arc and pre-batholithic crustal structure. 
    more » « less