We undertook Zr isotope measurements on zircon, titanite, biotite, amphibole, and whole rocks from the La Posta pluton (Peninsular Ranges, southern California) together with trace element analyses and U-Pb age measurements to understand the controls on Zr isotope fractionation in igneous rocks, including temperature, crystallization sequence, and kinetic effects. We find large (>0.6‰) Zr isotope fractionations (expressed as δ94/90Zr) between titanite and zircon forming at approximately the same temperature. Using equilibrium fractionation factors calculated from ionic and ab initio models, we infer the controls on Zr isotope evolution to include the relative order in which phases appear on the liquidus, with titanite fractionation resulting in isotopically lighter melt and zircon fractionation resulting in isotopically heavier melt. While these models of Zr fractionation can explain δ94/90Zr variations in zircon of up to ∼1.5‰, crystallization order, temperature and presence of co-crystallizing phases do not explain all aspects of the intracrystalline Zr isotopic distribution in zircons in the La Posta pluton or the large range of Zr isotopic values among zircons (>2‰). Without additional constraints, such as knowledge of co-crystallizing phases and a better understand of the true causes of Zr isotope fractionation, Zr isotopes in zircon remains an ambiguous proxy of magmatic evolution.
more »
« less
Significant Zr isotope variations in single zircon grains recording magma evolution history
Zircons widely occur in magmatic rocks and often display internal zonation finely recording the magmatic history. Here, we presented in situ high-precision (2SD <0.15‰ for δ 94 Zr) and high–spatial-resolution (20 µm) stable Zr isotope compositions of magmatic zircons in a suite of calc-alkaline plutonic rocks from the juvenile part of the Gangdese arc, southern Tibet. These zircon grains are internally zoned with Zr isotopically light cores and increasingly heavier rims. Our data suggest the preferential incorporation of lighter Zr isotopes in zircon from the melt, which would drive the residual melt to heavier values. The Rayleigh distillation model can well explain the observed internal zoning in single zircon grains, and the best-fit models gave average zircon–melt fractionation factors for each sample ranging from 0.99955 to 0.99988. The average fractionation factors are positively correlated with the median Ti-in-zircon temperatures, indicating a strong temperature dependence of Zr isotopic fractionation. The results demonstrate that in situ Zr isotope analyses would be another powerful contribution to the geochemical toolbox related to zircon. The findings of this study solve the fundamental issue on how zircon fractionates Zr isotopes in calc-alkaline magmas, the major type of magmas that led to forming continental crust over time. The results also show the great potential of stable Zr isotopes in tracing magmatic thermal and chemical evolution and thus possibly continental crustal differentiation.
more »
« less
- Award ID(s):
- 1725002
- NSF-PAR ID:
- 10294215
- Date Published:
- Journal Name:
- Proceedings of the National Academy of Sciences
- Volume:
- 117
- Issue:
- 35
- ISSN:
- 0027-8424
- Page Range / eLocation ID:
- 21125 to 21131
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
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
-
Abstract Ignimbrite flare-ups are rare periods of intense silicic volcanism during which the pyroclastic volume and eruptive frequency is more than an order of magnitude higher than background activity. Investigating the compositional differences between flare-up and steady-state magmas provides critical constraints on the petrogenetic causes for the event and can offer unique opportunities to investigate the role of large-scale tectonic or geodynamic processes in arc magmatism. In this study, we focus on the bimodal Deschutes Formation ignimbrite flare-up of Central Oregon, which erupted unusually high volumes of pyroclastic material 6.25–5.45 Ma from a new axis of volcanism in the Cascades arc. This episode is marked by increased eruption rates and eruption of more silicic compositions relative to the Quaternary Cascade arc, which rarely erupts rhyolites. Ignimbrites are crystal-poor (<10%) dacite to rhyolites (mostly 65–77 wt.% SiO2) with anhydrous mineral assemblages and higher FeO/MgO, Y, Eu/Eu*, MREE and Zr/Sr, indicating drier magmatic evolution compared to the Quaternary arc, and are more similar to those from the rear-arc High Lava Plains (HLP) province that lies to the east. Magnetite-ilmenite oxybarometry indicates that Deschutes Formation felsic magmas tend to be hotter and more reduced (NNO-1 to NNO) than the Quaternary arc (NNO to NNO + 1.5). Rhyolite-MELTS geobarometry suggests complex storage of diverse Deschutes Formation magmas within the shallow crust (50–250 MPa), and the common co-eruption of multiple plagioclase populations, pumice compositions, and compositionally banded pumice suggest variable degrees of mixing and mingling of distinct magmas. Deschutes magmas also have low δ18Oplagioclase values that indicate partial melting and assimilation of hydrothermally altered shallow crust. Trace element systematics and rhyolite-MELTS modeling suggests that felsic pumice cannot be produced by simple fractionation of co-erupted mafic pumice or basaltic lavas, and requires a crustal melting origin, and trace elements and Pb isotopes suggest that young mafic crust may have been the primary protolith. We suggest that partial melting produced low-Si rhyolite melt (~72 wt.%) that acted as both a parent for the most evolved rhyolites, and as a mixing endmember to create the dacite to rhyodacite magmas with heterogenous plagioclase populations. Unlike the predominantly calc-alkaline basalts erupted in the Quaternary Cascade arc, Deschutes Formation primary basalts are mostly low-K tholeiites, indicative of decompression melting. These are similar to the compositions erupted during a contemporaneous pulse of low-K tholeiite volcanism across the whole HLP that reached into the Cascades rear-arc. We suggest that intra-arc extension focused decompression melts from the back-arc into the arc and that tensional stresses allowed this high flux of hot-dry-reduced basalt throughout the crustal column, causing partial melting of mafic protoliths and the production of hot-dry-reduced rhyolite melts. Depletion of incompatible elements in successive rhyolites implies progressive depletion in fertility of the protolith. Extension also allowed for the establishment of a robust hydrothermal system, and assimilation of hydrothermally-altered rocks by magmas residing in a shallow, complex storage network lead to low δ18O melts. Our findings suggest the integral role that extensional tectonics played in producing an unusual ignimbrite flare-up of hot-dry-reduced rhyolite magmas that are atypical of the Cascades arc and may be an important contributor to flare-ups at arcs worldwide.more » « less
-
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
-
Zirconium (Zr) plays a key role in the development of phases like zircon (ZrSiO 4 ) and baddeleyite (ZrO 2 ) in magmatic systems. These minerals are crucial for the study of geologic time and crustal evolution, and their high resistivity to weathering and erosion results in their preservation on timescales of billions of years. Although zircon and baddeleyite may also preserve a robust record of Zr isotope behavior in high-temperature terrestrial environments, little is known about the factors that control Zr isotope partitioning in magmatic systems, the petrogenetic significance of fractionated compositions, or how these variations are recorded in Zr-rich accessory phases. Here, we describe a new analytical protocol for accurately determining the Zr stable isotope composition of zircon by multicollector-inductively coupled plasma-mass spectrometry (MC-ICP-MS), using the double-spike method to correct for procedural and instrumental mass bias. We apply this technique to test whether zircon crystallization in carbonatite magmatic systems is a driver of Zr isotope fractionation by interrogating the internal zonation of a zircon megacryst from the Mud Tank carbonatite (MTUR1). We find the MTUR1 megacryst to lack internal zoning within analytical uncertainties with a mean μ 94/90 Zr NIST = −55 ± 28 ppm (2 SD, n = 151), which suggests that zircon crystallization is not a driver of Zr isotope fractionation in carbonatite magmas. This observation is in stark contrast with those made in silicate magmatic systems, raising the possibility that the bonding environment of Zr 4+ ions may be fundamentally different in carbonatite vs. silicate melts. Because of its remarkable homogeneity, the MTUR1 megacryst is an ideal natural reference material for Zr isotopic analysis of zircon using both solution and spatially resolved methods. The reproducibility of a pure Zr solution and our chemically purified zircon fractions indicate that the external reproducibility of our method is on the order of ±28 ppm for μ 94/90 Zr, or ±7 ppm per amu, at 95% confidence.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.1073/pnas.2002053117
