Geochemical Perspectives Letters
Geochemical
Perspectives Letters
Geochemical
Perspectives
  • Submit here
  • Track your paper
  • For authors
  • e-Alerts
  • Home
  • About
    • About the journal
    • Editorial Board
    • Publication Policy
    • Publication Ethics
  • Submission & Review
    • Copyright & Permissions
    • Information for Authors
    • Information for Reviewers
  • Current issue
  • All issues
  • Submit
Select Page Menu

by admin | Nov 11, 2022 | mainpost, vol24

C.E. Bucholz, J. Liebmann, C.J. Spencer

24

2240

5

August

2022

18

October

2022

11

November

2022

12

16

0

Next article >> << Previous article

Secular variability in zircon phosphorus concentrations prevents simple petrogenetic classification

C.E. Bucholz1,

1Division of Geological and Planetary Sciences, California Institute of Technology, Pasadena, CA 91125

J. Liebmann2,

2School of Earth and Planetary Sciences, Curtin University, 6845, Perth, Australia 6845

C.J. Spencer3

3Department of Geological Sciences and Geological Engineering, Queen’s University, Kingston, Canada K7L 3N6

Affiliations | Corresponding Author | Cite as | Funding information

C.E. Bucholz
Email: cbucholz@caltech.edu

1Division of Geological and Planetary Sciences, California Institute of Technology, Pasadena, CA 91125
2School of Earth and Planetary Sciences, Curtin University, 6845, Perth, Australia 6845
3Department of Geological Sciences and Geological Engineering, Queen’s University, Kingston, Canada K7L 3N6

Bucholz, C.E., Liebmann, J., Spencer, C.J. (2022) Secular variability in zircon phosphorus concentrations prevents simple petrogenetic classification. Geochem. Persp. Let. 24, 12–16. https://doi.org/10.7185/geochemlet.2240

National Science Foundation (EAR-1943629).

Geochemical Perspectives Letters v24 | https://doi.org/10.7185/geochemlet.2240
Received 5 August 2022 | Accepted 18 October 2022 | Published 11 November 2022

Copyright © 2022 The Authors

Published by the European Association of Geochemistry
under Creative Commons License CC BY-NC-ND 4.0

Keywords: Phosphorus, peraluminous granite, zircon

PDF PDF+SI
  • Share this article

  • Article views:
    801

    Cumulative count of HTML views and PDF downloads.

  • Download Citation
  • Rights & Permissions


top

Abstract

Abstract | Introduction | Studied Samples and Methods | Results | What Do Zircon P Concentrations Represent? | Implications and Conclusion | Acknowledgements | References | Supplementary Information

Phosphorus (P) concentrations in zircon have been used to discriminate their derivation from metaluminous versus strongly peraluminous granites (SPGs) based on the empirical observation of lower P concentrations in zircon from Phanerozoic metaluminous versus peraluminous granites. Higher P concentrations in zircon from Phanerozoic SPGs reflect enhanced apatite solubility in peraluminous melts and overall higher P concentrations in peraluminous granites. However, SPGs derived from partial melting of Precambrian sedimentary rocks have lower P concentrations compared to Phanerozoic metaluminous granites, reflecting lower P concentrations in Precambrian versus Phanerozoic sedimentary sources. We demonstrate that zircons from Precambrian SPGs also have lower P concentrations compared to Phanerozoic counterparts, likely reflecting lower P concentrations in their parental melts. Applying the P-in-zircon proxy to the detrital zircon record does not effectively discriminate between metaluminous and peraluminous sources and underestimates contributions from peraluminous granites. Although detrital zircons are an important early Earth archive, a uniformitarian perspective cannot always be applied when using trace element proxies developed on Phanerozoic samples.

Figures

Figure 1 (a) Cumulate frequency distributions of SPG bulk rock P2O5. (b) Histograms of SPG bulk rock P2O5 with source rock and crystallisation ages >720 Ma (top) and <720 Ma (bottom). Data in (a) and (b) from Bucholz (2022). (c) Cumulative frequency distributions of P in zircon for Proterozoic and Archean SPGs, as well as Phanerozoic SPGs, “I-type” granites, and continental arc volcanic rocks. (d) Histograms of zircon P concentrations for Proterozoic/Archean SPGs and Phanerozoic SPGs. Phanerozoic SPG and I-type granite data in (c) and (d) are from Burnham and Berry (2017) and Zhu et al. (2020) and continental arc volcanic zircon from Grimes et al. (2015).

Figure 2 REE + Y versus P (in μmol/g) in zircon. Phanerozoic granite data are the same as in Figure 1. Proterozoic and Archean SPGs are from this study (filled symbols) and Lu et al. (2021) (open symbols).

Figure 3 Model results for melt P2O5 (wt. %) versus temperature during crystallisation of a SPG melt for different initial melt P2O5 (wt. %). Melt trajectories end once apatite saturation is reached. Zircon saturation temperatures are indicated with vertical red lines.

Figure 1 Figure 2 Figure 3

View all figures and tables





top

Introduction

Abstract | Introduction | Studied Samples and Methods | Results | What Do Zircon P Concentrations Represent? | Implications and Conclusion | Acknowledgements | References | Supplementary Information


Trace element concentrations in zircon are used to “fingerprint” the magma from which they crystallised (e.g., Grimes et al., 2015

Grimes, C.B., Wooden, J.L., Cheadle, M.J., John, B.E. (2015) “Fingerprinting” tectono-magmatic provenance using trace elements in igneous zircon. Contributions to Mineralogy and Petrology 170, 46. https://doi.org/10.1007/s00410-015-1199-3

). Of particular interest is identification of zircon from “I-type” or metaluminous granites versus those from “S-type” or strongly peraluminous granites (SPGs1; Burnham and Berry, 2017

Burnham, A.D., Berry, A.J. (2017) Formation of Hadean granites by melting of igneous crust. Nature Geoscience 10, 457–461. https://doi.org/10.1038/ngeo2942

; Trail et al., 2017

Trail, D., Tailby, N., Wang, Y., Mark Harrison, T., Boehnke, P. (2017) Aluminum in zircon as evidence for peraluminous and metaluminous melts from the Hadean to present. Geochemistry, Geophysics, Geosystems 18, 1580–1593. https://doi.org/10.1002/2016GC006794

). Understanding the relative contribution of zircon from metaluminous versus peraluminous granites throughout Earth’s history yields information on lithologic diversity and tectonic environments through time. SPGs dominantly form in collisional orogenic environments through the partial melting of sedimentary rocks (Nabelek, 2020

Nabelek, P.I. (2020) Petrogenesis of leucogranites in collisional orogens. Geological Society, London, Special Publications 491, 179–207. https://doi.org/10.1144/SP491-2018-181

). Thus, their inferred presence (or absence) through time yields information on both tectonic regimes and the availability of sedimentary rocks to be recycled into magmas during orogenic cycles.

One proposed proxy to discriminate zircon from metaluminous and peraluminous granites is phosphorus (P) concentrations in zircon (Burnham and Berry, 2017

Burnham, A.D., Berry, A.J. (2017) Formation of Hadean granites by melting of igneous crust. Nature Geoscience 10, 457–461. https://doi.org/10.1038/ngeo2942

). As the solubility of apatite increases with melt aluminium saturation index (Pichavant et al., 1992

Pichavant, M., Montel, J.-M., Richard, L.R. (1992) Apatite solubility in peraluminous liquids: Experimental data and an extension of the Harrison-Watson model. Geochimica et Cosmochimica Acta 56, 3855–3861. https://doi.org/10.1016/0016-7037(92)90178-L

; Wolf and London, 1994

Wolf, M.B., London, D. (1994) Apatite dissolution into peraluminous haplogranitic melts: An experimental study of solubilities and mechanisms. Geochimica et Cosmochimica Acta 58, 4127–4145. https://doi.org/10.1016/0016-7037(94)90269-0

), P concentrates in peraluminous melts during differentiation due to the lack of apatite precipitation (Bea et al., 1994

Bea, F., Pereira, M.D., Corretgé, L.G., Fershtater, G.B. (1994) Differentiation of strongly peraluminous, perphosphorus granites: The Pedrobernardo pluton, central Spain. Geochimica et Cosmochimica Acta 58, 2609–2627. https://doi.org/10.1016/0016-7037(94)90132-5

), a trend contrary to that observed in metaluminous melts (Lee and Bachmann, 2014

Lee, C.-T.A., Bachmann, O. (2014) How important is the role of crystal fractionation in making intermediate magmas? Insights from Zr and P systematics. Earth and Planetary Science Letters 393, 266–274. https://doi.org/10.1016/j.epsl.2014.02.044

). Indeed, Phanerozoic SPGs have significantly higher P2O5 (on average ∼0.25 wt. % and up to ∼1.5 wt. % P2O5; Bea et al., 1992

Bea, F., Fershtater, G., Corretgé, L.G. (1992) The geochemistry of phosphorus in granite rocks and the effect of aluminium. Lithos 29, 43–56. https://doi.org/10.1016/0024-4937(92)90033-U

; Bucholz, 2022

Bucholz, C.E. (2022) Coevolution of sedimentary and strongly peraluminous granite phosphorus records. Earth and Planetary Science Letters 596, 117795. https://doi.org/10.1016/j.epsl.2022.117795

; Fig. 1a,b) than metaluminous granites (which generally have <0.20 wt. % at bulk rock SiO2 > 65 wt. %). Further, upon apatite saturation, P concentrations are buffered at higher levels in peraluminous than metaluminous melts.


Figure 1 (a) Cumulate frequency distributions of SPG bulk rock P2O5. (b) Histograms of SPG bulk rock P2O5 with source rock and crystallisation ages >720 Ma (top) and <720 Ma (bottom). Data in (a) and (b) from Bucholz (2022)

Bucholz, C.E. (2022) Coevolution of sedimentary and strongly peraluminous granite phosphorus records. Earth and Planetary Science Letters 596, 117795. https://doi.org/10.1016/j.epsl.2022.117795

. (c) Cumulative frequency distributions of P in zircon for Proterozoic and Archean SPGs, as well as Phanerozoic SPGs, “I-type” granites, and continental arc volcanic rocks. (d) Histograms of zircon P concentrations for Proterozoic/Archean SPGs and Phanerozoic SPGs. Phanerozoic SPG and I-type granite data in (c) and (d) are from Burnham and Berry (2017)

Burnham, A.D., Berry, A.J. (2017) Formation of Hadean granites by melting of igneous crust. Nature Geoscience 10, 457–461. https://doi.org/10.1038/ngeo2942

and Zhu et al. (2020)

Zhu, Z., Campbell, I.H., Allen, C.M., Burnham, A.D. (2020) S-type granites: Their origin and distribution through time as determined from detrital zircons. Earth and Planetary Science Letters 536, 116140. https://doi.org/10.1016/j.epsl.2020.116140

and continental arc volcanic zircon from Grimes et al. (2015)

Grimes, C.B., Wooden, J.L., Cheadle, M.J., John, B.E. (2015) “Fingerprinting” tectono-magmatic provenance using trace elements in igneous zircon. Contributions to Mineralogy and Petrology 170, 46. https://doi.org/10.1007/s00410-015-1199-3

.
Full size image


Consequently, assuming equivalent zircon-melt P partition coefficients for metaluminous and peraluminous granitic systems, zircons from SPGs should have higher P contents than those from metaluminous granitic melts (Burnham and Berry, 2017

Burnham, A.D., Berry, A.J. (2017) Formation of Hadean granites by melting of igneous crust. Nature Geoscience 10, 457–461. https://doi.org/10.1038/ngeo2942

). Indeed, P in zircons from Phanerozoic peraluminous versus metaluminous granites have higher P concentrations (means of 1017 ± 53 versus 312 ± 48 ppm (2 s.e.); Burnham and Berry, 2017

Burnham, A.D., Berry, A.J. (2017) Formation of Hadean granites by melting of igneous crust. Nature Geoscience 10, 457–461. https://doi.org/10.1038/ngeo2942

; Zhu et al., 2020

Zhu, Z., Campbell, I.H., Allen, C.M., Burnham, A.D. (2020) S-type granites: Their origin and distribution through time as determined from detrital zircons. Earth and Planetary Science Letters 536, 116140. https://doi.org/10.1016/j.epsl.2020.116140

; Fig. 1c,d). Further, P and REE + Y concentrations are strongly correlated in zircon from Palaeozoic SPGs, suggesting the incorporation of P, REE, and Y into the zircon lattice via the coupled xenotime-type mechanism [(Y, REE)3+ + P5+ = Zr4+ + Si4+]. Based on these observations it has been proposed that P > 750 ppm and (REE + Y) < 1.15*P are robust criteria for identifying zircon from peraluminous granites (Burnham and Berry, 2017

Burnham, A.D., Berry, A.J. (2017) Formation of Hadean granites by melting of igneous crust. Nature Geoscience 10, 457–461. https://doi.org/10.1038/ngeo2942

).

However, P2O5 concentrations in SPGs are not elevated above those of metaluminous granites throughout Earth history (Bucholz, 2022

Bucholz, C.E. (2022) Coevolution of sedimentary and strongly peraluminous granite phosphorus records. Earth and Planetary Science Letters 596, 117795. https://doi.org/10.1016/j.epsl.2022.117795

). Maximum and average P2O5 concentrations in SPGs derived from sedimentary rocks deposited after 720 Ma are elevated above those derived from older source rocks (Fig. 1a,b). This temporal variation mirrors a similar increase in P concentrations in marine siliciclastic sediments deposited after 720 Ma, related to increased atmospheric and marine O2 levels and enhanced authigenic P burial (Reinhard et al., 2017

Reinhard, C.T., Planavsky, N.J., Gill, B.C., Ozaki, K., Robbins, L.J., Lyons, T.W., Fischer, W.W., Wang, C., Cole, D.B. and Konhauser, K.O. (2017) Evolution of the global phosphorus cycle. Nature 541, 386–389. https://doi.org/10.1038/nature20772

). After considering the factors affecting P concentrations in SPGs, Bucholz (2022)

Bucholz, C.E. (2022) Coevolution of sedimentary and strongly peraluminous granite phosphorus records. Earth and Planetary Science Letters 596, 117795. https://doi.org/10.1016/j.epsl.2022.117795

concluded that the temporal change in P concentrations is most likely explained as reflecting higher P concentrations in younger sedimentary source rocks. Elevated P concentrations in SPGs derived from sedimentary rocks deposited after 720 Ma could reflect both higher melt P concentrations and/or inherited P-rich restitic minerals (e.g., apatite, monazite, or P-bearing silicates).

If P concentrations in Archean and most Proterozoic SPGs are relatively low and if this, to some degree, represents lower melt P concentrations, zircon crystallising from such melts may record these low P concentrations. To test this hypothesis, we analysed trace element concentrations (including P) in zircon from Meso- to Paleoproterozoic and Archean SPGs.

1We prefer the use of strongly peraluminous granite over “S-type” granites as the former is purely a geochemical definition that avoids confusion associated with the “alphabet” granite classification scheme (c.f., Frost et al., 2001

Frost, C.D., Da Prat, F.A. (2021) Petrogenetic and tectonic interpretation of strongly peraluminous granitic rocks and their significance in the Archean rock record. American Mineralogist 106, 1195–1208. https://doi.org/10.2138/am-2022-8001

). However, we emphasise that the SPGs discussed here have all been previously demonstrated to have formed via the partial melting of metasedimentary rocks.


top

Studied Samples and Methods

Abstract | Introduction | Studied Samples and Methods | Results | What Do Zircon P Concentrations Represent? | Implications and Conclusion | Acknowledgements | References | Supplementary Information


We analysed trace elements in zircon (n = 201 analyses) from 16 previously characterised SPGs (Liebmann et al., 2021a

Liebmann, J., Spencer, C.J., Kirkland, C.L., Bucholz, C., He, X., Santosh, M., Xia, X.P., Martin, L., Evans, N.J. (2021a) Emergence of continents above sea‐level influences sediment melt composition. Terra Nova 33, 465–474. https://doi.org/10.1111/ter.12531

,b

Liebmann, J., Spencer, C.J., Kirkland, C.L., Bucholz, C.E., Xia, X.-P., Martin, L., Kitchen, N., Shumlyanskyy, L. (2021b) Coupling sulfur and oxygen isotope ratios in sediment melts across the Archean-Proterozoic transition. Geochimica et Cosmochimica Acta 307, 242–257. https://doi.org/10.1016/j.gca.2021.05.045

) from the southwestern USA (ca. 1.4 Ga), the North China Craton (ca. 1.9 and 2.5 Ga), Finland (ca. 1.8–1.9 Ga), and Ghana (ca. 2.2 Ga) (see Tables S-1 and S-2 for locality, lithologic, and bulk rock major elements). The average bulk rock P2O5 contents of these SPGs are 0.088 ± 0.036 wt. % (2 s.e.). U-Pb isotopes and trace elements (Ca, P, Y, Hf, and REEs) were collected simultaneously using split stream laser ablation inductively coupled plasma mass spectrometry at UC Santa Barbara. Full analytical details are given in the Supplementary Information (SI) and sample and standard analyses are given in Tables S-4 and S-5, respectively. Cathodoluminescence imaging and U-Pb ages were used to ensure analysis of magmatic zircon (i.e. related to the crystallisation of the SPG) and not inherited cores/grains (Fig. S-1). Any analyses of inherited cores (with older U-Pb ages) or with >5 % discordance of U-Pb ages were excluded. No apatite inclusions were observed. We are aware of only one other locality of >720 Ma SPGs where trace elements in magmatic zircon were analysed (Lu et al., 2021

Lu, G.-M., Spencer, C.J., Tian, Y., Wang, W. (2021) Significant Increase of Continental Freeboard During the Early Paleoproterozoic: Insights From Metasediment‐Derived Granites. Geophysical Research Letters 48, e2021GL096049. https://doi.org/10.1029/2021GL096049

) and include these data in our discussion.

top

Results

Abstract | Introduction | Studied Samples and Methods | Results | What Do Zircon P Concentrations Represent? | Implications and Conclusion | Acknowledgements | References | Supplementary Information


The P concentrations in Archean and Proterozoic SPGs are positively skewed (average = 132 ± 10 (2 s.e.) ppm; median =110 ppm; and including Lu et al. (2021)

Lu, G.-M., Spencer, C.J., Tian, Y., Wang, W. (2021) Significant Increase of Continental Freeboard During the Early Paleoproterozoic: Insights From Metasediment‐Derived Granites. Geophysical Research Letters 48, e2021GL096049. https://doi.org/10.1029/2021GL096049

data, average =196 ± 17 ppm; median = 138 ppm). In comparison, P concentrations in zircon from Phanerozoic SPGs are normally distributed with an average of 1017 ± 53 (2 s.e.) (Fig. 1c,d). Cumulative frequency distributions of P concentrations highlight this difference (Fig. 1c) and the mean P concentration in zircon from Precambrian and Phanerozoic SPGs are statistically distinguishable (P < 0.00001, Wilcoxon rank sum test). Zircon REE + Y concentrations from Precambrian SPGs, although weakly correlated with P (R2 = 0.63), are defined by a slope of >1.15, distinct from the Phanerozoic SPG trend (Fig. 2). Neither P nor total REE concentrations correlate with Ca concentrations, suggesting that P concentrations are not reflecting the incorporation of apatite inclusions in the ablation volume.


Figure 2 REE + Y versus P (in μmol/g) in zircon. Phanerozoic granite data are the same as in Figure 1. Proterozoic and Archean SPGs are from this study (filled symbols) and Lu et al. (2021)

Lu, G.-M., Spencer, C.J., Tian, Y., Wang, W. (2021) Significant Increase of Continental Freeboard During the Early Paleoproterozoic: Insights From Metasediment‐Derived Granites. Geophysical Research Letters 48, e2021GL096049. https://doi.org/10.1029/2021GL096049

(open symbols).
Full size image


top

What Do Zircon P Concentrations Represent?

Abstract | Introduction | Studied Samples and Methods | Results | What Do Zircon P Concentrations Represent? | Implications and Conclusion | Acknowledgements | References | Supplementary Information


The P concentration of zircon reflects equilibrium and/or disequilibrium processes during crystallisation from melt, including (1) equilibrium partitioning dependent on temperature, pressure, and melt composition, (2) partitioning reflecting melt-crystal boundary layer concentrations controlled by diffusion of P in the melt, or (3) surface controlled non-equilibrium growth (e.g., Watson and Liang, 1995

Watson, E.B., Liang, Y. (1995) A simple model for sector zoning in slowly grown crystals; implications for growth rate and lattice diffusion, with emphasis on accessory minerals in crustal rocks. American Mineralogist 80, 1179–1187. https://doi.org/10.2138/am-1995-11-1209

; Hofmann et al., 2009

Hofmann, A.E., Valley, J.W., Watson, E.B., Cavosie, A.J., Eiler, J.M. (2009) Sub-micron scale distributions of trace elements in zircon. Contributions to Mineralogy and Petrology 158, 317–335. https://doi.org/10.1007/s00410-009-0385-6

). In the first scenario, assuming that P partitioning is both Henrian and that equilibrium between zircon and melt occurred, then lower P concentrations in zircon from Archean/Proterozoic versus Phanerozoic SPGs could reflect either lower P concentrations in their parental melts or different conditions of crystallisation (e.g., melt composition, pressure, and temperature) which resulted in variable zircon-melt P partition coefficients. Phosphorus partitioning between melt and zircon as a function of melt composition, temperature, and pressure is poorly understood with limited experimental studies available (Rubatto and Hermann, 2007

Rubatto, D., Hermann, J. (2007) Experimental zircon/melt and zircon/garnet trace element partitioning and implications for the geochronology of crustal rocks. Chemical Geology 241, 38–61. https://doi.org/10.1016/j.chemgeo.2007.01.027

; Taylor et al., 2015

Taylor, R.J.M., Harley, S.L., Hinton, R.W., Elphick, S., Clark, C., Kelly, N.M. (2015) Experimental determination of REE partition coefficients between zircon, garnet and melt: a key to understanding high- T crustal processes. Journal of Metamorphic Geology 33, 231–248. https://doi.org/10.1111/jmg.12118

). However, Archean/Proterozoic SPGs are thought to have formed under similar conditions (e.g., pressure, temperature) to those in the Phanerozoic (Bucholz and Spencer, 2019

Bucholz, C.E., Spencer, C.J. (2019) Strongly Peraluminous Granites across the Archean–Proterozoic Transition. Journal of Petrology 60, 1299–1348. https://doi.org/10.1093/petrology/egz033

). Thus, in the first scenario of equilibrium partitioning, the simplest explanation is that lower P concentrations in zircon from Archean/Proterozoic SPGs, as compared to Phanerozoic ones, reflect lower melt P concentrations.

In the second scenario, the growth rate of zircon approaches (or exceeds) the diffusivity of P (and/or REE) in the adjacent melt. For hydrous rhyolites (∼6 wt. % H2O), experimentally calibrated P and Zr diffusivities (D) are similar with log(DP) and log(DZr) between −13.5 to −14.5 m2/s at 750–850 °C (see review of Zhang and Gan, 2022

Zhang, Y., Gan, T. (2022) Diffusion in Melts and Magmas. Reviews in Mineralogy and Geochemistry 87, 283–337. https://doi.org/10.2138/rmg.2022.87.07

). However, when zircon growth is sufficiently fast, P due to its general incompatibility is excluded from the zircon and builds up in a melt boundary layer. Thus, the growth of new zircon in equilibrium with this boundary layer acquires higher P concentrations elevated above that predicted by equilibrium partitioning with the bulk (i.e. far field) melt. However, all else being equal (e.g., P diffusivity and partition coefficients), melts with higher bulk P concentrations would produce boundary layers with higher P concentrations which would be reflected in higher zircon P concentrations.

Finally, in the third scenario, unequivocable evidence for surface controlled, non-equilibrium growth manifests as sector zoning (Watson and Liang, 1995

Watson, E.B., Liang, Y. (1995) A simple model for sector zoning in slowly grown crystals; implications for growth rate and lattice diffusion, with emphasis on accessory minerals in crustal rocks. American Mineralogist 80, 1179–1187. https://doi.org/10.2138/am-1995-11-1209

), however (sub-)μm scale oscillatory variations P concentrations have also been attributed to non-equilibrium growth (Hofmann et al., 2009

Hofmann, A.E., Valley, J.W., Watson, E.B., Cavosie, A.J., Eiler, J.M. (2009) Sub-micron scale distributions of trace elements in zircon. Contributions to Mineralogy and Petrology 158, 317–335. https://doi.org/10.1007/s00410-009-0385-6

). We generally avoided analysing portions of zircon that exhibited sector zoning (Fig. S-1). However, when analysed, sector zoned zircon had similar P concentrations to other zircon from the sample (c.f., 18IM15b; Fig. S-1). Further, the 25 μm laser spot used homogenised any visible oscillatory zoning (Fig. S-1), providing an integrated P concentration of the ablation volume. In zircon where multiple (2–3) spots were analysed, P concentrations agree within 25 ppm for 70 % and within 50 ppm for 80 % of zircon (Fig. S-3), suggesting that potential local variability induced by non-equilibrium growth was mostly homogenised.

No matter the P incorporation mechanism, our study demonstrates unequivocally that zircons from Archean and Proterozoic SPGs have lower P concentrations on average than their Phanerozoic SPGs (Figs. 1, 2). A parsimonious explanation for this observation is that lower P concentrations in zircons from Precambrian SPGs reflect lower P concentrations of their granitic parental melts. To explore this idea, we modelled melt P2O5 and zircon and apatite saturation during crystallisation of a SPG melt using alphaMELTS (Pichavant et al., 1992

Pichavant, M., Montel, J.-M., Richard, L.R. (1992) Apatite solubility in peraluminous liquids: Experimental data and an extension of the Harrison-Watson model. Geochimica et Cosmochimica Acta 56, 3855–3861. https://doi.org/10.1016/0016-7037(92)90178-L

; Smith and Asimow, 2005

Smith, P.M., Asimow, P.D. (2005) Adiabat_1ph: A new public front-end to the MELTS, pMELTS, and pHMELTS models: ADIABAT_1PH FRONT-END. Geochemistry, Geophysics, Geosystems 6. https://doi.org/10.1029/2004GC000816

; Boehnke et al., 2013

Boehnke, P., Watson, E.B., Trail, D., Harrison, T.M., Schmitt, A.K. (2013) Zircon saturation re-revisited. Chemical Geology 351, 324–334. https://doi.org/10.1016/j.chemgeo.2013.05.028

). Full modelling details are given in the SI. Calculations were performed for initial melt P2O5 contents typical of SPGs across Earth history (0.05, 0.10, 0.20, and 0.30 wt. %; Fig. 1a,b).

Due to compositional proximity to the granitic eutectic, zircon saturation temperatures were primarily controlled by increasing Zr concentrations in the melt with progressive differentiation. For initial Zr concentrations in the melt of 50, 100, 150, and 200 ppm, calculated zircon saturation temperatures are ∼727, 772, 800, and 822 °C, respectively. [Average Zr concentrations of pre-720 Ma SPGs from Bucholz (2022)

Bucholz, C.E. (2022) Coevolution of sedimentary and strongly peraluminous granite phosphorus records. Earth and Planetary Science Letters 596, 117795. https://doi.org/10.1016/j.epsl.2022.117795

are 120 ± 7 ppm (2 s.e.).] The calculated saturation temperatures agree with calculated Ti-in-zircon temperatures which on average are 810 ± 119 °C and 741 ± 105 °C (2 s.d.) assuming a melt TiO2 activity of 0.5 and 1, respectively (see SI for calculation details). For 100 ppm Zr in the initial melt, modelled melt P2O5 concentrations at zircon saturation are 0.08, 0.15, 0.30, and 0.45 wt. % for initial melt P2O5 contents of 0.05, 0.10, 0.20, and 0.30 wt. %, respectively (Fig. 3). Importantly, at zircon saturation, the melt P2O5 concentration has not significantly increased from its initial value as the melt fraction is still high (∼75 %). For a melt with initially low P2O5 (0.05–0.10 wt. %) typical for Precambrian SPGs, melt P2O5 increases significantly only within 20–30 °C of the granitic eutectic when the melt fraction decreases dramatically. Calculated melts ultimately reach apatite saturation 70–110 °C below that of the zircon saturation for a melt with 100 ppm initial Zr (Fig. 3). Therefore, although melt P2O5 increases due to delayed apatite saturation in these peraluminous melts, it remains low at zircon saturation in SPG melts with initially low P contents.


Figure 3 Model results for melt P2O5 (wt. %) versus temperature during crystallisation of a SPG melt for different initial melt P2O5 (wt. %). Melt trajectories end once apatite saturation is reached. Zircon saturation temperatures are indicated with vertical red lines.
Full size image


top

Implications and Conclusion

Abstract | Introduction | Studied Samples and Methods | Results | What Do Zircon P Concentrations Represent? | Implications and Conclusion | Acknowledgements | References | Supplementary Information


As SPGs predominantly form during collisional orogenesis, an accurate understanding of their temporal distribution provides information on the tectonic evolution and crustal recycling through time (Bucholz and Spencer, 2019

Bucholz, C.E., Spencer, C.J. (2019) Strongly Peraluminous Granites across the Archean–Proterozoic Transition. Journal of Petrology 60, 1299–1348. https://doi.org/10.1093/petrology/egz033

; Frost and Da Prat, 2021

Frost, C.D., Da Prat, F.A. (2021) Petrogenetic and tectonic interpretation of strongly peraluminous granitic rocks and their significance in the Archean rock record. American Mineralogist 106, 1195–1208. https://doi.org/10.2138/am-2022-8001

). Although detrital zircon can help characterise igneous rocks throughout Earth history, a uniformitarian approach cannot always be applied when using trace element proxies developed on Phanerozoic samples. Specifically, this study provides an example of how temporal variations in sedimentary rocks affect our interpretation of the igneous rock (and mineral) record. For example, the use of zircon P concentrations to identify metaluminous versus peraluminous source rocks has been applied to Archean and Hadean detrital zircon (Burnham and Berry, 2017

Burnham, A.D., Berry, A.J. (2017) Formation of Hadean granites by melting of igneous crust. Nature Geoscience 10, 457–461. https://doi.org/10.1038/ngeo2942

; Zhu et al., 2020

Zhu, Z., Campbell, I.H., Allen, C.M., Burnham, A.D. (2020) S-type granites: Their origin and distribution through time as determined from detrital zircons. Earth and Planetary Science Letters 536, 116140. https://doi.org/10.1016/j.epsl.2020.116140

). In particular, P and the correlation between P and REE + Y in Jack Hills zircon are more similar to Phanerozoic metaluminous granites and have been used to infer that the Jack Hills zircon crystallised from “TTG-like” magmas (Burnham and Berry, 2017

Burnham, A.D., Berry, A.J. (2017) Formation of Hadean granites by melting of igneous crust. Nature Geoscience 10, 457–461. https://doi.org/10.1038/ngeo2942

). However, the nature of the source rocks for the Jack Hills zircon is debated with elevated 18O/16O ratios (Cavosie et al., 2005

Cavosie, A.J., Valley, J.W., Wilde, S.A., E.I.M.F. (2005) Magmatic δ18O in 4400–3900 Ma detrital zircons: A record of the alteration and recycling of crust in the Early Archean. Earth and Planetary Science Letters 235, 663–681. https://doi.org/10.1016/j.epsl.2005.04.028

; Trail et al., 2007

Trail, D., Mojzsis, S.J., Harrison, T.M., Schmitt, A.K., Watson, E.B., Young, E.D. (2007) Constraints on Hadean zircon protoliths from oxygen isotopes, Ti-thermometry, and rare earth elements. Geochemistry, Geophysics, Geosystems 8, Q06014. https://doi.org/10.1029/2006GC001449

) and Al concentrations in Jack Hills zircon (Ackerson et al., 2021

Ackerson, M.R., Trail, D., Buettner, J. (2021) Emergence of peraluminous crustal magmas and implications for the early Earth. Geochemical Perspectives Letters 17, 50–54. https://doi.org/10.7185/geochemlet.2114

), potentially indicating that some Jack Hills zircons are from granites with a sedimentary source. Our results demonstrate that low P contents and P/(Y + REE) values <1 in Precambrian zircons cannot rule out crystallisation from a peraluminous granitic melt. Therefore, low P contents in Jack Hills zircon is consistent with crystallisation from both metaluminous and peraluminous melts, supporting inferences that the Jack Hills zircon could be sourced from a variety of granitic magmas (Bell, 2017

Bell, E. (2017) Ancient magma sources revealed. Nature Geoscience 10, 397–398. https://doi.org/10.1038/ngeo2955

; Ackerson et al., 2021

Ackerson, M.R., Trail, D., Buettner, J. (2021) Emergence of peraluminous crustal magmas and implications for the early Earth. Geochemical Perspectives Letters 17, 50–54. https://doi.org/10.7185/geochemlet.2114

).

Similarly, both cumulative distributions of P contents and P correlations with REE + Y in detrital zircons have been used to infer that SPG formation in the Archean was negligible and limited until the Neoproterozoic to Palaeozoic (Zhu et al., 2020

Zhu, Z., Campbell, I.H., Allen, C.M., Burnham, A.D. (2020) S-type granites: Their origin and distribution through time as determined from detrital zircons. Earth and Planetary Science Letters 536, 116140. https://doi.org/10.1016/j.epsl.2020.116140

). However, we again urge caution about extrapolating these metrics to Precambrian zircon when sedimentary source regions of SPGs had low P concentrations. Further, Phanerozoic SPGs with sedimentary source rocks deposited before 720 Ma also have low P concentrations (Fig. 1a) and should also contain magmatic zircon with low P. However, P concentrations of detrital zircon may be useful in local studies of young terranes to identify contributions from Phanerozoic P-rich SPGs. Identification of detrital zircon from SPGs with Precambrian source rocks is perhaps best done through a combination of O isotopes with new developing zircon proxies such as Al concentrations (Ackerson et al., 2021

Ackerson, M.R., Trail, D., Buettner, J. (2021) Emergence of peraluminous crustal magmas and implications for the early Earth. Geochemical Perspectives Letters 17, 50–54. https://doi.org/10.7185/geochemlet.2114

) or Si isotopes (Trail et al., 2018

Trail, D., Boehnke, P., Savage, P.S., Liu, M.-C., Miller, M.L., Bindeman, I. (2018) Origin and significance of Si and O isotope heterogeneities in Phanerozoic, Archean, and Hadean zircon. Proceedings of the National Academy of Sciences 115, 10287–10292. https://doi.org/10.1073/pnas.1808335115

).

top

Acknowledgements

Abstract | Introduction | Studied Samples and Methods | Results | What Do Zircon P Concentrations Represent? | Implications and Conclusion | Acknowledgements | References | Supplementary Information


We thank A. Kylander-Clark for his assistance with LA-ICPMS analyses. The constructive reviews of Jesse Reimink and an anonymous reviewer clarified and strengthened our manuscript. This work was supported by NSF grant EAR-1943629 to CB.

Editor: Maud Boyet

top

References

Abstract | Introduction | Studied Samples and Methods | Results | What Do Zircon P Concentrations Represent? | Implications and Conclusion | Acknowledgements | References | Supplementary Information

Ackerson, M.R., Trail, D., Buettner, J. (2021) Emergence of peraluminous crustal magmas and implications for the early Earth. Geochemical Perspectives Letters 17, 50–54. https://doi.org/10.7185/geochemlet.2114
Show in context

However, the nature of the source rocks for the Jack Hills zircon is debated with elevated 18O/16O ratios (Cavosie et al., 2005; Trail et al., 2007) and Al concentrations in Jack Hills zircon (Ackerson et al., 2021), potentially indicating that some Jack Hills zircons are from granites with a sedimentary source.
View in article
Therefore, low P contents in Jack Hills zircon is consistent with crystallisation from both metaluminous and peraluminous melts, supporting inferences that the Jack Hills zircon could be sourced from a variety of granitic magmas (Bell, 2017; Ackerson et al., 2021).
View in article
Identification of detrital zircon from SPGs with Precambrian source rocks is perhaps best done through a combination of O isotopes with new developing zircon proxies such as Al concentrations (Ackerson et al., 2021) or Si isotopes (Trail et al., 2018).
View in article


Bea, F., Fershtater, G., Corretgé, L.G. (1992) The geochemistry of phosphorus in granite rocks and the effect of aluminium. Lithos 29, 43–56. https://doi.org/10.1016/0024-4937(92)90033-U
Show in context

Indeed, Phanerozoic SPGs have significantly higher P2O5 (on average ∼0.25 wt. % and up to ∼1.5 wt. % P2O5; Bea et al., 1992; Bucholz, 2022; Fig. 1a,b) than metaluminous granites (which generally have <0.20 wt. % at bulk rock SiO2 > 65 wt. %).
View in article


Bea, F., Pereira, M.D., Corretgé, L.G., Fershtater, G.B. (1994) Differentiation of strongly peraluminous, perphosphorus granites: The Pedrobernardo pluton, central Spain. Geochimica et Cosmochimica Acta 58, 2609–2627. https://doi.org/10.1016/0016-7037(94)90132-5
Show in context

As the solubility of apatite increases with melt aluminium saturation index (Pichavant et al., 1992; Wolf and London, 1994), P concentrates in peraluminous melts during differentiation due to the lack of apatite precipitation (Bea et al., 1994), a trend contrary to that observed in metaluminous melts (Lee and Bachmann, 2014).
View in article


Bell, E. (2017) Ancient magma sources revealed. Nature Geoscience 10, 397–398. https://doi.org/10.1038/ngeo2955
Show in context

Therefore, low P contents in Jack Hills zircon is consistent with crystallisation from both metaluminous and peraluminous melts, supporting inferences that the Jack Hills zircon could be sourced from a variety of granitic magmas (Bell, 2017; Ackerson et al., 2021).
View in article


Boehnke, P., Watson, E.B., Trail, D., Harrison, T.M., Schmitt, A.K. (2013) Zircon saturation re-revisited. Chemical Geology 351, 324–334. https://doi.org/10.1016/j.chemgeo.2013.05.028
Show in context

To explore this idea, we modelled melt P2O5 and zircon and apatite saturation during crystallisation of a SPG melt using alphaMELTS (Pichavant et al., 1992; Smith and Asimow, 2005; Boehnke et al., 2013).
View in article


Bucholz, C.E. (2022) Coevolution of sedimentary and strongly peraluminous granite phosphorus records. Earth and Planetary Science Letters 596, 117795. https://doi.org/10.1016/j.epsl.2022.117795
Show in context

Indeed, Phanerozoic SPGs have significantly higher P2O5 (on average ∼0.25 wt. % and up to ∼1.5 wt. % P2O5; Bea et al., 1992; Bucholz, 2022; Fig. 1a,b) than metaluminous granites (which generally have <0.20 wt. % at bulk rock SiO2 > 65 wt. %).
View in article
Data in (a) and (b) from Bucholz (2022).
View in article
However, P2O5 concentrations in SPGs are not elevated above those of metaluminous granites throughout Earth history (Bucholz, 2022).
View in article
After considering the factors affecting P concentrations in SPGs, Bucholz (2022) concluded that the temporal change in P concentrations is most likely explained as reflecting higher P concentrations in younger sedimentary source rocks.
View in article
[Average Zr concentrations of pre-720 Ma SPGs from Bucholz (2022) are 120 ± 7 ppm (2 s.e.).]
View in article


Bucholz, C.E., Spencer, C.J. (2019) Strongly Peraluminous Granites across the Archean–Proterozoic Transition. Journal of Petrology 60, 1299–1348. https://doi.org/10.1093/petrology/egz033
Show in context

However, Archean/Proterozoic SPGs are thought to have formed under similar conditions (e.g., pressure, temperature) to those in the Phanerozoic (Bucholz and Spencer, 2019).
View in article
As SPGs predominantly form during collisional orogenesis, an accurate understanding of their temporal distribution provides information on the tectonic evolution and crustal recycling through time (Bucholz and Spencer, 2019; Frost and Da Prat, 2021).
View in article


Burnham, A.D., Berry, A.J. (2017) Formation of Hadean granites by melting of igneous crust. Nature Geoscience 10, 457–461. https://doi.org/10.1038/ngeo2942
Show in context

Of particular interest is identification of zircon from “I-type” or metaluminous granites versus those from “S-type” or strongly peraluminous granites (SPGs1; Burnham and Berry, 2017; Trail et al., 2017).
View in article
One proposed proxy to discriminate zircon from metaluminous and peraluminous granites is phosphorus (P) concentrations in zircon (Burnham and Berry, 2017).
View in article
Phanerozoic SPG and I-type granite data in (c) and (d) are from Burnham and Berry (2017) and Zhu et al. (2020) and continental arc volcanic zircon from Grimes et al. (2015).
View in article
Consequently, assuming equivalent zircon-melt P partition coefficients for metaluminous and peraluminous granitic systems, zircons from SPGs should have higher P contents than those from metaluminous granitic melts (Burnham and Berry, 2017).
View in article
Indeed, P in zircons from Phanerozoic peraluminous versus metaluminous granites have higher P concentrations (means of 1017 ± 53 versus 312 ± 48 ppm (2 s.e.); Burnham and Berry, 2017; Zhu et al., 2020; Fig. 1c,d).
View in article
Based on these observations it has been proposed that P > 750 ppm and (REE + Y) < 1.15*P are robust criteria for identifying zircon from peraluminous granites (Burnham and Berry, 2017).
View in article
For example, the use of zircon P concentrations to identify metaluminous versus peraluminous source rocks has been applied to Archean and Hadean detrital zircon (Burnham and Berry, 2017; Zhu et al., 2020).
View in article
In particular, P and the correlation between P and REE + Y in Jack Hills zircon are more similar to Phanerozoic metaluminous granites and have been used to infer that the Jack Hills zircon crystallised from “TTG-like” magmas (Burnham and Berry, 2017).
View in article


Cavosie, A.J., Valley, J.W., Wilde, S.A., E.I.M.F. (2005) Magmatic δ18O in 4400–3900 Ma detrital zircons: A record of the alteration and recycling of crust in the Early Archean. Earth and Planetary Science Letters 235, 663–681. https://doi.org/10.1016/j.epsl.2005.04.028
Show in context

However, the nature of the source rocks for the Jack Hills zircon is debated with elevated 18O/16O ratios (Cavosie et al., 2005; Trail et al., 2007) and Al concentrations in Jack Hills zircon (Ackerson et al., 2021), potentially indicating that some Jack Hills zircons are from granites with a sedimentary source.
View in article


Frost, C.D., Da Prat, F.A. (2021) Petrogenetic and tectonic interpretation of strongly peraluminous granitic rocks and their significance in the Archean rock record. American Mineralogist 106, 1195–1208. https://doi.org/10.2138/am-2022-8001
Show in context

As SPGs predominantly form during collisional orogenesis, an accurate understanding of their temporal distribution provides information on the tectonic evolution and crustal recycling through time (Bucholz and Spencer, 2019; Frost and Da Prat, 2021).
View in article


Grimes, C.B., Wooden, J.L., Cheadle, M.J., John, B.E. (2015) “Fingerprinting” tectono-magmatic provenance using trace elements in igneous zircon. Contributions to Mineralogy and Petrology 170, 46. https://doi.org/10.1007/s00410-015-1199-3
Show in context

Trace element concentrations in zircon are used to “fingerprint” the magma from which they crystallised (e.g., Grimes et al., 2015).
View in article
Phanerozoic SPG and I-type granite data in (c) and (d) are from Burnham and Berry (2017) and Zhu et al. (2020) and continental arc volcanic zircon from Grimes et al. (2015).
View in article


Hofmann, A.E., Valley, J.W., Watson, E.B., Cavosie, A.J., Eiler, J.M. (2009) Sub-micron scale distributions of trace elements in zircon. Contributions to Mineralogy and Petrology 158, 317–335. https://doi.org/10.1007/s00410-009-0385-6
Show in context

The P concentration of zircon reflects equilibrium and/or disequilibrium processes during crystallisation from melt, including (1) equilibrium partitioning dependent on temperature, pressure, and melt composition, (2) partitioning reflecting melt-crystal boundary layer concentrations controlled by diffusion of P in the melt, or (3) surface controlled non-equilibrium growth (e.g., Watson and Liang, 1995; Hofmann et al., 2009).
View in article
Finally, in the third scenario, unequivocable evidence for surface controlled, non-equilibrium growth manifests as sector zoning (Watson and Liang, 1995), however (sub-)μm scale oscillatory variations P concentrations have also been attributed to non-equilibrium growth (Hofmann et al., 2009).
View in article


Lee, C.-T.A., Bachmann, O. (2014) How important is the role of crystal fractionation in making intermediate magmas? Insights from Zr and P systematics. Earth and Planetary Science Letters 393, 266–274. https://doi.org/10.1016/j.epsl.2014.02.044
Show in context

As the solubility of apatite increases with melt aluminium saturation index (Pichavant et al., 1992; Wolf and London, 1994), P concentrates in peraluminous melts during differentiation due to the lack of apatite precipitation (Bea et al., 1994), a trend contrary to that observed in metaluminous melts (Lee and Bachmann, 2014).
View in article


Liebmann, J., Spencer, C.J., Kirkland, C.L., Bucholz, C., He, X., Santosh, M., Xia, X.P., Martin, L., Evans, N.J. (2021a) Emergence of continents above sea‐level influences sediment melt composition. Terra Nova 33, 465–474. https://doi.org/10.1111/ter.12531
Show in context

We analysed trace elements in zircon (n = 201 analyses) from 16 previously characterised SPGs (Liebmann et al., 2021a,b) from the southwestern USA (ca. 1.4 Ga), the North China Craton (ca. 1.9 and 2.5 Ga), Finland (ca. 1.8–1.9 Ga), and Ghana (ca. 2.2 Ga) (see Tables S-1 and S-2 for locality, lithologic, and bulk rock major elements).
View in article


Liebmann, J., Spencer, C.J., Kirkland, C.L., Bucholz, C.E., Xia, X.-P., Martin, L., Kitchen, N., Shumlyanskyy, L. (2021b) Coupling sulfur and oxygen isotope ratios in sediment melts across the Archean-Proterozoic transition. Geochimica et Cosmochimica Acta 307, 242–257. https://doi.org/10.1016/j.gca.2021.05.045
Show in context

We analysed trace elements in zircon (n = 201 analyses) from 16 previously characterised SPGs (Liebmann et al., 2021a,b) from the southwestern USA (ca. 1.4 Ga), the North China Craton (ca. 1.9 and 2.5 Ga), Finland (ca. 1.8–1.9 Ga), and Ghana (ca. 2.2 Ga) (see Tables S-1 and S-2 for locality, lithologic, and bulk rock major elements).
View in article


Lu, G.-M., Spencer, C.J., Tian, Y., Wang, W. (2021) Significant Increase of Continental Freeboard During the Early Paleoproterozoic: Insights From Metasediment‐Derived Granites. Geophysical Research Letters 48, e2021GL096049. https://doi.org/10.1029/2021GL096049
Show in context

We are aware of only one other locality of >720 Ma SPGs where trace elements in magmatic zircon were analysed (Lu et al., 2021) and include these data in our discussion.
View in article
The P concentrations in Archean and Proterozoic SPGs are positively skewed (average = 132 ± 10 (2 s.e.) ppm; median =110 ppm; and including Lu et al. (2021) data, average =196 ± 17 ppm; median = 138 ppm).
View in article
Proterozoic and Archean SPGs are from this study (filled symbols) and Lu et al. (2021) (open symbols).
View in article


Nabelek, P.I. (2020) Petrogenesis of leucogranites in collisional orogens. Geological Society, London, Special Publications 491, 179–207. https://doi.org/10.1144/SP491-2018-181
Show in context

SPGs dominantly form in collisional orogenic environments through the partial melting of sedimentary rocks (Nabelek, 2020).
View in article


Pichavant, M., Montel, J.-M., Richard, L.R. (1992) Apatite solubility in peraluminous liquids: Experimental data and an extension of the Harrison-Watson model. Geochimica et Cosmochimica Acta 56, 3855–3861. https://doi.org/10.1016/0016-7037(92)90178-L
Show in context

As the solubility of apatite increases with melt aluminium saturation index (Pichavant et al., 1992; Wolf and London, 1994), P concentrates in peraluminous melts during differentiation due to the lack of apatite precipitation (Bea et al., 1994), a trend contrary to that observed in metaluminous melts (Lee and Bachmann, 2014).
View in article
To explore this idea, we modelled melt P2O5 and zircon and apatite saturation during crystallisation of a SPG melt using alphaMELTS (Pichavant et al., 1992; Smith and Asimow, 2005; Boehnke et al., 2013).
View in article


Reinhard, C.T., Planavsky, N.J., Gill, B.C., Ozaki, K., Robbins, L.J., Lyons, T.W., Fischer, W.W., Wang, C., Cole, D.B. and Konhauser, K.O. (2017) Evolution of the global phosphorus cycle. Nature 541, 386–389. https://doi.org/10.1038/nature20772
Show in context

This temporal variation mirrors a similar increase in P concentrations in marine siliciclastic sediments deposited after 720 Ma, related to increased atmospheric and marine O2 levels and enhanced authigenic P burial (Reinhard et al., 2017).
View in article


Rubatto, D., Hermann, J. (2007) Experimental zircon/melt and zircon/garnet trace element partitioning and implications for the geochronology of crustal rocks. Chemical Geology 241, 38–61. https://doi.org/10.1016/j.chemgeo.2007.01.027
Show in context

Phosphorus partitioning between melt and zircon as a function of melt composition, temperature, and pressure is poorly understood with limited experimental studies available (Rubatto and Hermann, 2007; Taylor et al., 2015).
View in article


Smith, P.M., Asimow, P.D. (2005) Adiabat_1ph: A new public front-end to the MELTS, pMELTS, and pHMELTS models: ADIABAT_1PH FRONT-END. Geochemistry, Geophysics, Geosystems 6. https://doi.org/10.1029/2004GC000816
Show in context

To explore this idea, we modelled melt P2O5 and zircon and apatite saturation during crystallisation of a SPG melt using alphaMELTS (Pichavant et al., 1992; Smith and Asimow, 2005; Boehnke et al., 2013).
View in article


Taylor, R.J.M., Harley, S.L., Hinton, R.W., Elphick, S., Clark, C., Kelly, N.M. (2015) Experimental determination of REE partition coefficients between zircon, garnet and melt: a key to understanding high- T crustal processes. Journal of Metamorphic Geology 33, 231–248. https://doi.org/10.1111/jmg.12118
Show in context

Phosphorus partitioning between melt and zircon as a function of melt composition, temperature, and pressure is poorly understood with limited experimental studies available (Rubatto and Hermann, 2007; Taylor et al., 2015).
View in article


Trail, D., Mojzsis, S.J., Harrison, T.M., Schmitt, A.K., Watson, E.B., Young, E.D. (2007) Constraints on Hadean zircon protoliths from oxygen isotopes, Ti-thermometry, and rare earth elements. Geochemistry, Geophysics, Geosystems 8, Q06014. https://doi.org/10.1029/2006GC001449
Show in context

However, the nature of the source rocks for the Jack Hills zircon is debated with elevated 18O/16O ratios (Cavosie et al., 2005; Trail et al., 2007) and Al concentrations in Jack Hills zircon (Ackerson et al., 2021), potentially indicating that some Jack Hills zircons are from granites with a sedimentary source.
View in article


Trail, D., Tailby, N., Wang, Y., Mark Harrison, T., Boehnke, P. (2017) Aluminum in zircon as evidence for peraluminous and metaluminous melts from the Hadean to present. Geochemistry, Geophysics, Geosystems 18, 1580–1593. https://doi.org/10.1002/2016GC006794
Show in context

Of particular interest is identification of zircon from “I-type” or metaluminous granites versus those from “S-type” or strongly peraluminous granites (SPGs1; Burnham and Berry, 2017; Trail et al., 2017).
View in article


Trail, D., Boehnke, P., Savage, P.S., Liu, M.-C., Miller, M.L., Bindeman, I. (2018) Origin and significance of Si and O isotope heterogeneities in Phanerozoic, Archean, and Hadean zircon. Proceedings of the National Academy of Sciences 115, 10287–10292. https://doi.org/10.1073/pnas.1808335115
Show in context

Identification of detrital zircon from SPGs with Precambrian source rocks is perhaps best done through a combination of O isotopes with new developing zircon proxies such as Al concentrations (Ackerson et al., 2021) or Si isotopes (Trail et al., 2018).
View in article


Watson, E.B., Liang, Y. (1995) A simple model for sector zoning in slowly grown crystals; implications for growth rate and lattice diffusion, with emphasis on accessory minerals in crustal rocks. American Mineralogist 80, 1179–1187. https://doi.org/10.2138/am-1995-11-1209
Show in context

The P concentration of zircon reflects equilibrium and/or disequilibrium processes during crystallisation from melt, including (1) equilibrium partitioning dependent on temperature, pressure, and melt composition, (2) partitioning reflecting melt-crystal boundary layer concentrations controlled by diffusion of P in the melt, or (3) surface controlled non-equilibrium growth (e.g., Watson and Liang, 1995; Hofmann et al., 2009).
View in article
Finally, in the third scenario, unequivocable evidence for surface controlled, non-equilibrium growth manifests as sector zoning (Watson and Liang, 1995), however (sub-)μm scale oscillatory variations P concentrations have also been attributed to non-equilibrium growth (Hofmann et al., 2009).
View in article


Wolf, M.B., London, D. (1994) Apatite dissolution into peraluminous haplogranitic melts: An experimental study of solubilities and mechanisms. Geochimica et Cosmochimica Acta 58, 4127–4145. https://doi.org/10.1016/0016-7037(94)90269-0
Show in context

As the solubility of apatite increases with melt aluminium saturation index (Pichavant et al., 1992; Wolf and London, 1994), P concentrates in peraluminous melts during differentiation due to the lack of apatite precipitation (Bea et al., 1994), a trend contrary to that observed in metaluminous melts (Lee and Bachmann, 2014).
View in article


Zhang, Y., Gan, T. (2022) Diffusion in Melts and Magmas. Reviews in Mineralogy and Geochemistry 87, 283–337. https://doi.org/10.2138/rmg.2022.87.07
Show in context

In the second scenario, the growth rate of zircon approaches (or exceeds) the diffusivity of P (and/or REE) in the adjacent melt. For hydrous rhyolites (∼6 wt. % H2O), experimentally calibrated P and Zr diffusivities (D) are similar with log(DP) and log(DZr) between −13.5 to −14.5 m2/s at 750–850 °C (see review of Zhang and Gan, 2022).
View in article


Zhu, Z., Campbell, I.H., Allen, C.M., Burnham, A.D. (2020) S-type granites: Their origin and distribution through time as determined from detrital zircons. Earth and Planetary Science Letters 536, 116140. https://doi.org/10.1016/j.epsl.2020.116140
Show in context

Phanerozoic SPG and I-type granite data in (c) and (d) are from Burnham and Berry (2017) and Zhu et al. (2020) and continental arc volcanic zircon from Grimes et al. (2015).
View in article
Indeed, P in zircons from Phanerozoic peraluminous versus metaluminous granites have higher P concentrations (means of 1017 ± 53 versus 312 ± 48 ppm (2 s.e.); Burnham and Berry, 2017; Zhu et al., 2020; Fig. 1c,d).
View in article
For example, the use of zircon P concentrations to identify metaluminous versus peraluminous source rocks has been applied to Archean and Hadean detrital zircon (Burnham and Berry, 2017; Zhu et al., 2020).
View in article
Similarly, both cumulative distributions of P contents and P correlations with REE + Y in detrital zircons have been used to infer that SPG formation in the Archean was negligible and limited until the Neoproterozoic to Palaeozoic (Zhu et al., 2020).
View in article



top

Supplementary Information

Abstract | Introduction | Studied Samples and Methods | Results | What Do Zircon P Concentrations Represent? | Implications and Conclusion | Acknowledgements | References | Supplementary Information


The Supplementary Information includes:
  • Analytical methods
  • Ti-in-zircon Thermometry
  • Zircon and Apatite Saturation Modelling
  • Tables S-1 to S-5
  • Figures S-1 to S-4
  • Supplementary Information References


  • Download Table S-4 (xlsx)

    Download Table S-5 (xlsx)

    Download the Supplementary Information (PDF)
    top

    Figures



    Figure 1 (a) Cumulate frequency distributions of SPG bulk rock P2O5. (b) Histograms of SPG bulk rock P2O5 with source rock and crystallisation ages >720 Ma (top) and <720 Ma (bottom). Data in (a) and (b) from Bucholz (2022)

    Bucholz, C.E. (2022) Coevolution of sedimentary and strongly peraluminous granite phosphorus records. Earth and Planetary Science Letters 596, 117795. https://doi.org/10.1016/j.epsl.2022.117795

    . (c) Cumulative frequency distributions of P in zircon for Proterozoic and Archean SPGs, as well as Phanerozoic SPGs, “I-type” granites, and continental arc volcanic rocks. (d) Histograms of zircon P concentrations for Proterozoic/Archean SPGs and Phanerozoic SPGs. Phanerozoic SPG and I-type granite data in (c) and (d) are from Burnham and Berry (2017)

    Burnham, A.D., Berry, A.J. (2017) Formation of Hadean granites by melting of igneous crust. Nature Geoscience 10, 457–461. https://doi.org/10.1038/ngeo2942

    and Zhu et al. (2020)

    Zhu, Z., Campbell, I.H., Allen, C.M., Burnham, A.D. (2020) S-type granites: Their origin and distribution through time as determined from detrital zircons. Earth and Planetary Science Letters 536, 116140. https://doi.org/10.1016/j.epsl.2020.116140

    and continental arc volcanic zircon from Grimes et al. (2015)

    Grimes, C.B., Wooden, J.L., Cheadle, M.J., John, B.E. (2015) “Fingerprinting” tectono-magmatic provenance using trace elements in igneous zircon. Contributions to Mineralogy and Petrology 170, 46. https://doi.org/10.1007/s00410-015-1199-3

    .
    Back to article


    Figure 2 REE + Y versus P (in μmol/g) in zircon. Phanerozoic granite data are the same as in Figure 1. Proterozoic and Archean SPGs are from this study (filled symbols) and Lu et al. (2021)

    Lu, G.-M., Spencer, C.J., Tian, Y., Wang, W. (2021) Significant Increase of Continental Freeboard During the Early Paleoproterozoic: Insights From Metasediment‐Derived Granites. Geophysical Research Letters 48, e2021GL096049. https://doi.org/10.1029/2021GL096049

    (open symbols).
    Back to article


    Figure 3 Model results for melt P2O5 (wt. %) versus temperature during crystallisation of a SPG melt for different initial melt P2O5 (wt. %). Melt trajectories end once apatite saturation is reached. Zircon saturation temperatures are indicated with vertical red lines.
    Back to article

    • Contact us
    • |
    • Subscribe
    Sign up for our newsletter
    • Connect with us
    • facebook
    • twitter
    • youtube

    Geochemical Perspectives Letters is a registered trademark of the European Association of Geochemistry
    ISSN 2410-339X (print) | ISSN 2410-3403 (online)
    EAG Privacy Policy