Abstract
The composition and origin of Earth’s early crust remains hotly debated. Here we use partition coefficients to invert the trace element composition of 4.3–3.3 Gyr Jack Hills zircons to calculate the composition of the melts from which they crystallised. Using this approach, the average SiO2 content of these melts was 59 ± 6 wt. % with Th/Nb, Dy/Yb and Sr/Y ratios of 2.7 ± 1.9, 0.9 ± 0.2 and 1.6 ± 0.7, respectively. Such features strongly indicate that the protolith for the Jack Hills zircons was not an intra-plate mafic rock, nor a TTG (tondjhemite-tonalite-granodiorite) or a Sudbury-like impact melt. Instead, the inferred equilibrium melts are much more similar to andesites formed in modern subduction settings. We find no evidence for any secular variation between 4.3 and 3.3 Gyr implying little change in the composition or tectonic affinity of the Earth’s early crust from the Hadean to Mesoarchaean.
Similar content being viewed by others
Introduction
The composition and tectonic affiliation of the Earths earliest crust are two of the most fundamental questions in Earth science1,2,3,4 but challenging to constrain due to the lack of suitably aged rocks. Many studies have concentrated on 4.4–3.3 Gyr old aged zircons extracted from metasediments in the Jack Hills in Western Australia1,5,6,7,8,9 but consensus on the nature of their protolith(s) remains elusive. Some inferences about the bulk composition of their protoliths have come from the Hf isotope composition of the zircons5,8,9. One thing that is not in dispute is that these zircons derive from igneous source rocks7,10,11. Most crystallised at temperatures around 700 ± 50 °C12,13 but the parental magmas would likely have been hotter. It is generally accepted that elevated U/Th ratios (>0.1) indicate that the zircons derive from igneous source rocks7,10,11 but thereafter opinions differ. Some advocate low temperature, hydrous granites as the protoliths, implying conditions similar to those on Earth today14,15. Others infer more mafic, possibly intra-plate basaltic source rocks9,16 or even that the zircons crystallised from impact melts17,18,19. The latter interpretations do not require plate tectonics. Other studies have suggested that plate tectonics was in a stop—start mode on the early Earth3,20 which would predict that key geochemical signatures (e.g., SiO2 content, Th/Nb ratios) alternate between within-plate and subduction like over time during the Archaean. As many of these conclusions are based on the Hf isotope signature of the zircons1,5,8,9 there is a need for a complementary approach. Whilst mostly concentrating on rare earth element (REE) abundances, zircon trace element data have been collected in several different studies11,21 and fingerprinting of these also indicate the Jack Hills zircons most likely crystallised from continental igneous rocks rather than fractionates of mafic intra-plate melts5,21. However, additional information can be extracted by exploiting a wider range of trace element data from the Jack Hills zircons. These suggest that they crystallised from melts with trace element signatures akin to modern arc andesites making an early onset of modern style plate tectonics entirely permissible.
Results
We have prepared a new aliquot of zircons from Jack Hills metasediment and obtained their U-Pb ages using standard SHRIMP methodology (see Methods section and Supplementary Data 1 and Supplementary Data 2). The ages range from 4.3 to 3.3 Ga consistent with earlier work5,6,7. We then used a standard laser-extraction technique (LA-ICP-MS) to obtain trace element data from the same location from which the ages were obtained (see Methods section and Supplementary Table S2). Next, we used zircon/melt partition coefficients to calculate the trace element composition of the melts from which the zircons crystallised (i.e., their protoliths). For this purpose the choice of partition coefficients is critical. The majority of studies do not report large ion lithophile and high field strength element data that are especially critical here. Accordingly, we have used LA-ICP-MS partition coefficients from the only experimental study that analysed both a wide range of trace elements and that also showed these to be consistent with the lattice strain model22. These experiments were conducted on an andesitic composition that, as we demonstrate below seem highly applicable to the Jack Hills zircons. Moreover, although a temperature dependence on absolute partition coefficient values has been well demonstrated, the ratios of partition coefficient (as we use throughout) are far less prone to this problem22,23. Moreover, zircons from plutonic rocks often contain micro-inclusions and their hosts are notoriously heterogeneous, rendering most empirical partition coefficients significantly inaccurate.
The new results are illustrated on a series of time-sequence diagrams in Fig. 1. Th/Y ratios were used to calculate the SiO2 content of the equilibrium melts (see Methods section) and these range from 51 to 68 wt.% with an average of 59 wt.% (Fig. 1a). These overlap estimates for the present day continental crust that has a broadly andesitic composition24. Th/Nb ratios provide one measure of the presence or absence of a negative Nb anomaly that is absent in mid-ocean ridge and intra-plate rocks. The melt Th/Nb ratios we have estimated from the Jack Hills zircons are 2.7 ± 1.9, significantly higher than any modern oceanic rocks and extend from values akin to typical continental crust to even higher values (Fig. 1b). Finally, we also calculated the temperatures at which these zircons crystallised yielding 700 ± 50 °C, consistent with previous observations12,13. However, zircon is typically a late crystallised phase in magmatic systems and so these must be regarded as minimum magmatic temperatures if the melts from which the zircons formed evolved from more mafic parental melts. In respect to this it is interesting that the two oldest zircons appear to record higher crystallisation temperatures than the remaining zircons. A key observation is that neither the calculated melt SiO2 contents and Th/Nb ratios or the zircon crystallisation temperatures show any statistically meaningful (R2 > 0.5) secular change from the Hadean to Mesoarchaean (i.e., between 4.3 and 3.3 Ga).
a Melt SiO2 calculated from inferred melt Th/Y is similar to average bulk crust. b Melt Th/Nb ratios are distinct from mid-ocean ridge and ocean island basalts but similar to or higher than average crust. c Zircon crystallisation temperatures are similar to the granite minimum though the parental magmas would have been hotter. Errors are smaller than the symbol size. There is no secular change in any of these parameters and all are features of typical subduction zone andesites. Average compositions of continental crustal24 and MORB and OIB43 are also shown. Error bar in a is 1σ (see Methods section).
Discussion
The combination of elevated SiO2 contents and elevated Th/Nb ratios are characteristics of andesitic melts formed in a modern subduction setting. Although some remain sceptical, plots like the Th/Yb versus Nb/Yb diagram shown in Fig. 2a have been repeatedly appraised to help distinguish the tectonic setting in which ancient rocks formed25,26. As can be seen the protoliths of the Jack Hills zircons clearly plot well above the fields for oceanic rocks and in or above the field for modern arc rocks3. Trondhjemite-tonalite-granodiorite suites (TTGs) characterise much of the felsic portion of the Archaean geological record and might be considered an obvious potential candidate source for the Jack Hills zircons. Many TTGs have elevated Th/Nb (and Sr/Y) ratios and that has sparked much debate as to whether these rock associations formed via subduction26. However, the Jack Hills protolihts do not overlap the field of TTGs on Fig. 2a. They are also distinct from the composition of rocks from Sudbury that are inferred to reflect crystallisation from an impact-derived melt (Fig. 2a).
a Trace element discrimination diagram25 showing comparison of melts in equilibrium with Jack Hills zircons with the fields for modern arc rocks, mid-ocean ridge and ocean island basalts. b Plot of Th/Nb versus Dy/Yb for melts in equilibrium with the Jack Hills zircons along with vectors indicating the contrasting effects of subduction and melting (or crystallisation) in the presence of garnet. Black circles assume experimental zircon/melt partition coefficients22. For comparison, the grey circles assume empirical zircon/melt partition coefficients28. Also shown are the composition of Sudbury rocks (white squares) inferred to have formed from impact-induced melting44 and the TTG field for trondhjemite-tonalite-granodiorite suites (data compilation from GEOROC). Average crustal composition24 (CC) and MORB43 are also shown. Error bars in b are 1σ (see Methods section).
A caveat to the preceding discussion is the choice of zircon/melt partition coefficients as noted earlier. In order to estimate the effect of differing partition coefficients we also applied the partition coefficients from a recent empirical study of a high silica rhyolite glass27. These should avoid many of the problems with whole rock estimates and were also shown to be consistent with the lattice strain model. Unlike Th and the rare earth elements (REE), Nb is especially vulnerable to the choice of zircon/melt partition coefficient because it is quite incompatible in zircon and there can be analytical difficulties in measuring Nb22. For example, Nb is significantly more incompatible in the empirical partitioning data27. We note that the latter were obtained from a high silica rhyolite whereas the broadly andesitic SiO2 content inferred for the Jack Hills zircon protoliths is more consistent with the andesitic starting composition used in the experimental study22. Nevertheless, if we use the empirical values27 the calculated melts in equilibrium with the Jack Hills zircons shift to higher Nb/Yb such that they can overlap the TTG field on Fig. 2a. Irrespective, they remain displaced well above the fields for oceanic rocks (i.e., at elevated Th/Nb).
Another well-established feature of the TTGs is that they have steep heavy REE patterns and elevated Sr/Y ratios (>10) indicative of an important role for garnet either during partial melting or crystal fractionation27. REE partitioning into garnet is well constrained and ratios of the partition coefficients for the heavy REE (which are highly compatible in zircon) appear particularly robust22,28. Thus, a second striking result of our investigation is that the inferred Jack Hills protoliths had flat REE patterns and so their Dy/Yb ratios do not overlap the TTG field on Fig. 2b and they have an average Sr/Y ratio of only 1.6 ± 0.7 (not shown). Their Dy/Yb ratios are also distinct from the composition of rocks from the Sudbury intrusion (Fig. 2b).
What then, can be said about the nature of the Jack Hills zircon protoliths? The most robust conclusions are that they had intermediate SiO2 contents and elevated Th/Nb ratios coupled with low Dy/Yb ratios. These results do not change even if a very different set of partition coefficients are used (see Fig. 2). Since impact melts have the composition of their target rock it is unsurprising that the Sudbury rocks lie close to average continental crust (Fig. 2b). As shown on Fig. 2b, the melts in equilibrium with the Jack Hills zircons are most similar to modern day arc lavas. Interestingly the average composition of the inferred Jack Hills protolith is very similar to the average low-Ti enriched basaltic andesites from the Nuvvuagittuq Greenstone Belt (SiO2 ~ 55 wt.%; Th/Nb ~ 0.8) that may be 3.8 or 4.3 Gyr in age29 and have been argued to derive from a stratigraphic sequence similar to those observed in present day subduction initiation sequences2. Experiments have shown that partial melting of those rocks can reproduce the associated TTG rocks30. Partial melting of the Jack Hills protoliths in the presence of residual garnet provides a viable explanation for the origin of TTGs. Finally, we find no evidence for within-plate geochemical signatures or alternation between with-plate and subduction like signatures over time. The presence of a subduction signature does not prove plate tectonics operated in the early Archaean. However, our findings are entirely permissible of an onset of modern-style plate tectonics very early in Earth history consistent with recent geodynamic modelling31 and Si and Mo isotope studies32,33.
Methods
Zircons and SHRIMP age dating
Grains for this study were extracted from one of the vials of non-magnetic zircon prepared from the original Jack Hills conglomerate sample obtained from the W74 site34. Two mounts were prepared with a total of approximately 50 zircons on each, along with several pieces of the CZ3 standard with a 206Pb/238U age of 564 Ma35. The zircons were cast in epoxy resin, ground, polished and cleaned prior to gold coating. Images were taken in both reflected and transmitted light in order to identify any imperfections or inclusions, and the mounts were then imaged in cathodoluminescence (CL). Based on these images, not all zircons were found to be suitable for analysis due to cracks and inclusions, and a total of 37 grains were analysed on mount 1 and 24 on mount 2. Initial analyses involved single-cycle runs on a SHRIMP II at Curtin University to determine the approximate age of the selected grains and to ascertain how many were Hadean. It was established that the grains were representative of the W74 zircon population and all were then re-analysed using 6-cycle runs through the mass stations following published methods36. A mass resolution of ~5000 was recorded during measurement of the Pb/Pb and U/Pb isotopic ratios and the latter were normalised to those measured on the standard zircon [CZ3–(206Pb/238U = 0.0914)]. The uncertainty associated with measurement of U/Pb for the standard at 1σ was 1.52% for the mount 1 session and 1.32% for the mount 2 session. The common lead correction was modelled on the composition of Broken Hill lead and the data were processed using Squid 137 and Isoplot 3.7538. The data are contained in Supplementary Data 1. Discordant zircons (>5%) and any indicating chemical evidence of alteration39 were removed from dataset prior to further processing and interpretation.
Trace element determinations
Trace-element concentrations of zircons were measured using laser ablation inductively coupled plasma mass spectrometry (LA-ICP-MS) using an Agilent 7700× quadrupole ICP-MS coupled with a Photon Machines Excite Analyte 193 nm excimer laser ablation system with HelEx sample cell at Macquarie GeoAnalytical (MQGA), Macquarie University. The data were integrated near the surface to be as close as possible to where the U-Pb ages were obtained. A fluence of 7.59 J/cm2 and a repetition rate of 5 Hz were adopted for the analytical conditions for the laser. A laser beam of 50 µm was used for the zircon analysis. Each measurement consists of 60 s of background and 120 s of ablation. NIST-610 glass standards were analysed at the beginning and the end of each analysis series to correct the machine drift. Concentrations for the trace elements were obtained through calibration of relative element sensitivities using the NIST-610 standard as the external calibration standard and Zr content was used for internal calibration. Basaltic reference material, BCR-2G, and zircon reference materials, GJ and 91500, were run before the unknown samples to monitor the accuracy and reproducibility of the measurements. Data reduction was carried out online using the GLITTER software40. The spectra were trimmed to avoid both surface contamination and clear depth related zoning although some unavoidable variation in U and Th concentrations remains. P and REE data were the used to identify and discard any zircons containing evidence for alteration39. Trace element concentrations of BCR-2G and 91500 obtained in this study are in good agreement with recommended values41,42. The data are contained in Supplementary Data 2. Errors for calculated melt composition trace element ratios have been propagated from the external reproducibility of the laser ablation analyses and incorporate the uncertainties on the partition coefficients22.
Calculation of protolith SiO2 contents
In order to derive the SiO2 contents of the equilibrium melts we regressed 41186 analyses between 50 and 70 wt.% SiO2 from an initial data set of >100,000 whole rock obtained from the GEOROC database and unfiltered for age (see below). The data was binned into 1 wt.% SiO2 intervals and a linear regression was used to determine the form Th/Y = 0.0269 × SiO2–1.3169 with an R2 value of 0.9812. Since both Th and Y are compatible in zircon with Dzircon/melt = 1.1 and 2.1, respectively22 and neither have redox sensitive partitioning their inferred melt compositions can be regarded as robust. The standard error about the mean was then subsequently calculated for each 1% bin interval. Each bin contained a minimum of n = 527 (66 wt.% SiO2) to a maximum of 4936 (50 wt.% SiO2) analyses. Further regressions were then conducted on the maximum (Th/Y = 0.0263 × SiO2–1.3066, R2 = 0.9812) and minimum (Th/Y = 0.0257 × SiO2–1.2963; R2 = 0.9808) values for each bin interval in order to calculate the maximum variation of SiO2 within individual bin intervals. This corresponded to ±1.1 wt.% SiO2 at 69 wt.% SiO2 and was a minimum of 0.75 at 51 wt.% SiO2 (Supplementary Fig. 1). Hence we conservatively report the purely statistical errors on the calculated SiO2 to be ±1.1 wt.% SiO2 (1 SE) and when this is combined with the analytical errors from the trace element analyses and errors for the partition coefficients the combined errors on calculated melt SiO2 is ±3.5 wt.%. The errors on calculated SiO2 are significantly less than the overall calculated variation in SiO2 and thus do not impinge upon the major conclusions drawn from the calculations.
As noted above, the dataset used for the regression was unfiltered for age. If we use only Archaean magmatic rocks we obtain a regression that is Th/Y = 0.1 × SiO2–6. This leads to even higher inferred SiO2 contents for the Jack Hills protoliths for 2 reasons we prefered to retain the original calculations. First, our original regression gives a “conservative” estimate of protolith SiO2. Second, The Archaean-only regression is naturally dominated by TTGs yet our Fig. 2b indicates that TTGs at not an appropriate the source of the Jack Hills zircons.
References
Harrison, T. M. et al. Heterogeneous Hadean Hafnium: evidence for continental crust at 4.4 to 4.5 Ga. Science 310, 1947 (2005).
Turner, S., Rushmer, T., Reagan, M. & Moyen, J.-F. Heading down early on? Start of subduction on Earth. Geology 42, 139–142 (2014).
O’Neill, C., Turner, S., Rushmer, T. The inception of plate tectonics: a record of failure. Philos. Trans. Roy. Soc. A376, https://doi.org/10.1098/rsta.2017.0414 (2018).
Johnson, T. E., Brown, M., Gardiner, N. J., Kirkland, C. J. & Smithies, R. H. Earth’s first stable continents did not form by subduction. Nature 543, 239–242 (2017).
Wilde, S. A., Valley, J. W., Peck, W. H. & Graham, C. M. Evidence from detrital zircons for the existence of continental crust and oceans on the Earth 4.4 Ga ago. Nature 409, 175–178 (2001).
Cavoisie, et al. (eds) Precambrian Ophiolites and Related Rocks, Vol. 15, Ch. 2.5, Developments in Precambrian Geology (Elsevier B. V., 2007).
Wang, Q. & Wilde, S. A. New constraints on the Hadean to Proterozoic history of the Jack Hills belt, Western Australia. Gondw. Res. 55, 74–91 (2018).
Blichert-Toft, J. & Albarède, F. Hafnium isotopes in Jack Hills zircons and the formation of the Hadean crust. Earth Planet. Sci. Lett. 265, 686–702 (2008).
Kemp, A. I. S. et al. Hadean crustal evolution revisited: new constraints from Pb-Hf isotope systematics of the Jack Hills zircons. Earth Planet. Sci. Lett. 296, 45–56 (2010).
Burnham, A. D. & Berry, A. J. Formation of Hadean granites by melting of igneous crust. Nat. Geosci. 10, 457–462 (2017).
Crowley, J. L., Myers, J. S., Sylvester, P. J. & Cox, R. A. Detrital zircon from the Jack Hills and Mount Narryer, Western Australia: evidence for diverse > 4.0 Ga source rocks. J. Geol. 113, 239–263 (2005).
Harrison, T. M., Watson, E. B. & Aikman, A. B. Temperature spectra of zircon crystallization in plutonic rocks. Geology 35, 635–638 (2007).
Fu, B. et al. Ti-in-zircon thermometry: applications and limitations. Contrib. Mineral. Petrol. 156, 197–215 (2008).
Mojzsis, S. J., Harrison, T. M. & Pigeon, R. T. Oxygen-isotope evidence from ancient zircons for liquid water at the Earth’s surface 4300 Myr ago. Nature 409, 178–181 (2001).
Trail, D., Watson, E. B. & Tailby, N. D. The oxidation state of Hdean magmas and implications for early Earth’s atmosphere. Nature 480, 79–82 (2011).
Dhuime, B., Wuestefeld, A. & Hawkesworth, C. J. Emergence of modern continental crust about 3 billion years ago. Nat. Geosci. 8, 552–555 (2015).
Darling, J., Storey, C. & Hawkesworth, C. Impact melt sheet zircons and their implications for the Hadean crust. Geology 37, 927–930 (2009).
Kenny, G. G., Whitehouse, M. J. & Kamber, B. S. Differentiated impact melt sheets may be a potential source of Hadean detrital zircon. Geology 44, 435–438 (2016).
Marchi, S. et al. Widespread mixing and burial of Earth’s Hadean crust by asteroid impacts. Nature 511, 578–582 (2014).
Moyen, J.-F. & Jvan Hunen., J. Short term episodicity of Archean subduction. Geology 40, 451–454 (2012).
Grimes, C. B. et al. Trace element chemistry of zircons from oceanic crust: a method for distinguishing detrital zircon provenance. Geology 35, 643–646 (2007).
Burnham, A. D. & Berry, A. J. 2012. An experimental study of trace element partitioning between zircon and melt as a function of oxygen fugacity. Geochim. Cosmochim. Acta 95, 196–212 (2012).
Blundy, J. D. & Wood, B. J. Prediction of crystal-melt partition coefficients from elastic moduli. Nature 372, 452–454 (1994).
Rudnick, R. L. & Gao, S. Composition of the continental crust. Treatise Geochem. 3, 1–64 (2003).
Pearce, J. A. Geochemical fingerprinting of oceanic basalts with application to ophiolite classification and the search for Archean oceanic crust. Lithos 100, 14–48 (2008).
Pearce, J. A. & Reagan, M. K. identification, classification, and interpretation of boninites from Anthropocene to Eoarchean using Si-Mg-Ti systematics. Geosphere 15, 1008–1037 (2019).
Moyen, J.-F. & Martin, H. Forty years of TTG research. Lithos 148, 312–336 (2012).
Padilla, A. J. & Gualda, A. R. Crystal-melt elemental partitioning in silicic magmatic systems: an example from the Peach Spring Tuff high-silica rhyolite, Southwest USA. Chem. Geol. 440, 326–344 (2016).
O’Neil, J., Carlson, R. W., Paquette, J.-L. & Francis, D. Formation age and metamorphic history of the Nuvvuagittuq Greenstone Belt. Precamb. Res. 220, 23–44 (2012).
Adam, J., Rushmer, T., O’Neil, J. & Francis, D. Hadean greenstones from the Nuvvuagittuq fold belt and the origin of the Earth’s early continental crust. Geology 40, 363–366 (2012).
Weller, O. M., Copley, A., Miller, W. G. R., Palin, R. M. & Dyck, B. The relationship between mantle potential temperature and oceanic lithospheric buoyancy. Earth Planet. Sci. Lett. 518, 86–99 (2019).
Deng, Z., et al. An oceanic subduction origin of Archaean granitoids revealed by silicon isotopes. Nat. Geosci. 12, 774–779 (2019).
McCoy-West, A. J., et al. Extensive crustal extraction in Earth’s early history inferred from molybdenum isotopes. Nat. Geosci. 12, 946–951 (2019).
Compston, W. & Pidgeon, R. T. Jack Hills, evidence of more very old detrital zircons in Western Australia. Nature 321, 766–769 (1986).
Nelson, D. R. Compilation of SHRIMP U–Pb geochronological data 1999. Geological Survey of Western Australia Record 2000/2 (2002).
Williams, I. S. U–Th–Pb geochronology by ion microprobe. (eds. McKibben, M. A., Shanks III, W. C. & Ridley, W. I. Applications of microanalytical techniques to understanding mineralizing processes. Rev. Econom. Geol. 7, 1–35 (1998).
Ludwig, K. R. Users manual for Squid 1.02. Berkeley Geochronol. Cent. Spec. Publ. 1a, 19 (2001).
Ludwig, K. R. Isoplot 3.75: a geochronological toolkit for Microsoft Excel. Berkeley Geochronology Center Special Publication 5 (2012).
Bell, E. A., Boehnke, P. & Harrison, T. M. Recovering the primary geochemistry of Jack Hills zircons through quantitative estimates of chemical alteration. Geochim. Cosmochim. Acta 191, 187–202 (2016).
Griffin, W. L., Powell, W. J., Pearson, N. J., O’Reilly, S. Y. GLITTER: data reduction software for laser ablation ICP-MS. (ed. Sylvester, P.), Laser Ablation-ICP-MS in the Earth Sciences. Mineralogical Association of Canada Short Course Series 40, 204–207 (Appendix 2) (2008).
Jochum, K. P. et al. GeoReM: A new geochemical database for reference materials and isotopic standards. Geostand. Geoanalyt. Res. 29, 333–338 (2007).
Wiedenbeck, M. et al. Further Characterisation of the 91500 Zircon Crystal. Geostand. Geoanalyt. Res. 28, 9–39 (2004).
Sun, S.-S. & McDonough, W. F. Chemical and isotopic systematics of oceanic basalts: implications for mantle composition and processes. Geol. Soc. Lond. Spec. Publ. 42, 313–345 (1989).
Lightfoot, P. C., Keays, R. R., Morrison, G. G., Bite, A. & Farrell, K. P. Geochemical relationships in the Sudbury Igneous Complex: origin of the main mass and offset dykes. Econom. Geol. 92, 289–307 (1997).
Acknowledgements
S.T. thanks the Akademie der Wissenschaften zu Göttingen for a Gauss Professorship during which the paper was written. We thank Craig O’Neill and Chris Hawkesworth for stimulating discussions about the onset of plate tectonics.
Author information
Authors and Affiliations
Contributions
S.W. selected the zircons and obtained the U-Pb age data, S.T. obtained the trace element data with Y.-J. Lai. S.T., G.W. and B.S. calculated and interpreted the equilibrium melt compositions and all authors contributed to the writing of the paper.
Corresponding author
Ethics declarations
Competing interests
The authors declare no competing interests.
Additional information
Peer review information Nature Communications thanks the anonymous reviewer(s) for their contribution to the peer review of this work. Peer reviewer reports are available.
Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Turner, S., Wilde, S., Wörner, G. et al. An andesitic source for Jack Hills zircon supports onset of plate tectonics in the Hadean. Nat Commun 11, 1241 (2020). https://doi.org/10.1038/s41467-020-14857-1
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/s41467-020-14857-1