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 | Dec 17, 2024 | mainpost, vol33

D. Trail, N.D. Tailby, M.R. Ackerson, M. Barboni, M. Newville, A. Lanzirotti, K.D. McKeegan

33

2448

12

August

2024

25

November

2024

17

December

2024

27

31

0

Next article >> << Previous article

Production of highly silicic 3.9 to 4.27 Ga crust on the Moon

D. Trail1,

1Department of Earth and Environmental Sciences, University of Rochester, Rochester NY USA

N.D. Tailby2,

2Division of Earth Sciences, School of Environmental and Rural Science, University of New England, Armidale, NSW, Australia

M.R. Ackerson3,

3Department of Mineral Sciences, National Museum of Natural History, Smithsonian Institution, Washington, DC

M. Barboni4,

4School of Earth and Space Exploration, Arizona State University, Tempe, AZ 85287, USA

M. Newville5,

5Center for Advanced Radiation Sources, University of Chicago, Chicago, Illinois, 60637, USA

A. Lanzirotti5,

5Center for Advanced Radiation Sources, University of Chicago, Chicago, Illinois, 60637, USA

K.D. McKeegan6

6Department of Earth, Planetary, and Space Sciences, University of California – Los Angeles, Los Angeles, CA 90095, USA

Affiliations | Corresponding Author | Cite as | Funding information

D. Trail
Email: dtrail@ur.rochester.edu

1Department of Earth and Environmental Sciences, University of Rochester, Rochester NY USA
2Division of Earth Sciences, School of Environmental and Rural Science, University of New England, Armidale, NSW, Australia
3Department of Mineral Sciences, National Museum of Natural History, Smithsonian Institution, Washington, DC
4School of Earth and Space Exploration, Arizona State University, Tempe, AZ 85287, USA
5Center for Advanced Radiation Sources, University of Chicago, Chicago, Illinois, 60637, USA
6Department of Earth, Planetary, and Space Sciences, University of California – Los Angeles, Los Angeles, CA 90095, USA

Trail, D., Tailby, N.D., Ackerson, M.R., Barboni, M., Newville, M., Lanzirotti, A., McKeegan, K.D. (2024) Production of highly silicic 3.9 to 4.27 Ga crust on the Moon. Geochem. Persp. Let. 33, 27–31. https://doi.org/10.7185/geochemlet.2448

NASA Solar System Workings grant 80NSSC20K1039. National Science Foundation – Earth Sciences (EAR – 1634415) NASA’s Planetary Science Enabling Facilities (PSEF) program (grant number 80NSSC23K0196). U.S. Department of Energy (DOE) Office of Science User Facility operated for the DOE Office of Science by Argonne National Laboratory under Contract No. DE-AC02-06CH11357.

Geochemical Perspectives Letters v33 | https://doi.org/10.7185/geochemlet.2448
Received 12 August 2024 | Accepted 25 November 2024 | Published 17 December 2024

Copyright © 2024 The Authors

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

Keywords: lunar, silicic, zircon, XANES

PDF PDF+SI
  • Share this article

  • Article views:
    2,683

    Cumulative count of HTML views and PDF downloads.

  • Download Citation
  • Rights & Permissions


top

Abstract

Abstract | Introduction | Results | Acknowledgements | References | Supplementary Information

The preserved remnants of a planetary body’s early crust represent a unique window into chemical and physical processes that shaped its beginning. While iron- and magnesium-rich magmas with low silica (SiO2) are common throughout a rocky body’s active lifetime, the timing of silica-rich crust production during the first 500 Myr of evolution remain uncertain. Silica-rich crust is relatively buoyant, influences the topography of rocky worlds, and concentrates heat-producing radioactive elements. However, identifying early silica-rich crust is challenging due to sample limitation and an inability to directly date silica-rich material like quartz and high-silica glass. Here, we show that the atomic environment around trace amounts of titanium in the mineral zircon is sensitive to the silica content of the magma it crystallised from, providing a new tool to recognise ancient silica-rich crust. By applying this method to 4.27 to 3.93 billion year old lunar zircons, we present evidence for early production of silica-rich crust on the Moon, suggesting this may be a common feature of planetary evolution.

Figures

Figure 1 XANES spectral analysis of Ti coordination in zircon. (a) Two standard spectra with Ti present in only 4-fold coordination (cristobalite); and a model spectrum for Ti in 8-fold coordination in zircon along the a axis (Tailby et al., 2011). (inset) Linear mixing of the pre-edge feature showing model spectra for different fractions of Ti in 4-fold coordination. (b) Example pre-edge features for experimental zircon for silica saturated and silica undersaturated (Tailby et al., 2011) experiments. The 1300 °C, 1 atm, SiO2 undersaturated experiment (srilankite (Ti,Zr)O4, rutile, zircon) has a pre-edge intensity similar to that of the 4-fold standard, whereas the intensity of the pre-edge feature for the two silica saturated experiments is muted meaning that a significant fraction of 8-fold Ti is present in zircon. The energy shift observed in the pre-edge spectra of the silica undersaturated samples is common, and has been shown to be related to crystal orientation (Tailby et al., 2011; Crisp et al., 2023). Differences in pre-edge intensity are due to orientation and/or sample heterogeneity (Tailby et al., 2011). Spectra were carefully evaluated for inclusions (see Supplementary Information and Fig. S-1). (c) Fraction of 4-fold Ti in zircon, each analysis from the three experiments calculated from a linear combination of the two standard spectra presented in (a) (Tailby et al., 2011; Crisp et al., 2023). When aqtzSiO2  ≈ 1, only ∼50 % of the Ti is occupying the 4-fold Si site, whereas with aqtzSiO2  ≈ 0.6 nearly 100 % of the Ti is in the 4-fold Si site. Error bars represent the 1 s.d. for each group. Experimental products can be found in Table S-1.

Figure 2 Ti XANES spectra and coordination in lunar zircon. (a, b) Two representative grain analyses showing XANES spectra before and after a 90° rotation. The top panel shows a significant change in pre-edge absorption intensity, while the lower panel, more typical, shows virtually no change. (c) Calculation of Ti in 4-fold coordination, with points averaging lunar zircon measurements in two crystal orientations, plotted against the previously reported 207Pb/206Pb age (Trail et al., 2020). Dashed lines show the average fraction of 4-fold Ti at 1300 °C, and shaded yellow areas represent the 1 s.d. of the fraction in synthetic zircon.

Figure 3 Th contents of highly silicic lunar systems. Circles represent calculated Th contents of igneous systems derived from zircon-melt Th experimental partitioning (see Supplementary Information; Fig. S-5 and Table S-5). Yellow squares represent Th contents for younger silicic systems investigated and quantified by remote sensing (Hagerty et al., 2006; Wilson et al., 2015). The ages for the silicic systems identified by remote sensing data are estimated through crater counting. The average Th contents from the PKT mare terrain (4.9 ppm) and other mare (2.2 ppm) are also shown for comparison (Jolliff et al., 2000). MMD1/2, first and second eruptions of Middle Mairan Dome; MTD, Mairan T Dome; GDδ/γ, Gruithuisen Dome (δ) and (γ); LM, Lassell Massif; HD, Hansteen Alpha; CBVC, Compton–Bel’kovich Volcanic Complex.

Figure 4 Simplified illustration (not to scale) of the lunar silicic systems whose formation was initiated by partial melting of fertile source rocks such as a KREEP basalt by basaltic underplating (Hagerty et al., 2006). High viscosity minimum melt (?) silicic systems may have developed not far from final emplacement location and/or made significant use of the higher porosity of the fractured crust/mega-regolith, including extensional features associated with impact basins (Valencia et al., 2024). Early silica-rich domes were likely destroyed by large basin-forming impacts, or covered by impact ejecta (Qiu et al., 2023).

Figure 1 Figure 2 Figure 3 Figure 4

View all figures and tables





top

Introduction

Abstract | Introduction | Results | Acknowledgements | References | Supplementary Information


The emergence of buoyant, silica-rich (>60 wt. %) crust on a rocky planetary body is a key marker of crustal reprocessing and/or extreme differentiation with crucial implications for a planet’s evolution and possibly the origin of life (Harrison, 2020

Harrison, T.M. (2020) Hadean Earth. Springer, Cham. https://doi.org/10.1007/978-3-030-46687-9

). Earth is the only known celestial body where silica-rich crust production occurs today, but whether silicic magma was present early or evolved gradually remains debated. Some models infer that there was continuous production and recycling of early high silica crust (Harrison, 2020

Harrison, T.M. (2020) Hadean Earth. Springer, Cham. https://doi.org/10.1007/978-3-030-46687-9

) while others favour the absence of a stable volume of buoyant continental material before 4 Ga (Cawood et al., 2022

Cawood, P.A., Chowdhury, P., Mulder, J.A., Hawkesworth, C.J., Capitanio, F.A., Gunawardana, P.M., Nebel, O. (2022) Secular Evolution of Continents and the Earth System. Reviews of Geophysics 60, e2022RG000789. https://doi.org/10.1029/2022RG000789

). Other solar system bodies – like the Moon, Mars, and Vesta – are composed predominantly of a primordial low silica (∼40–50 wt. %) crust but lack clear signs of early silicic magma production.

The Moon exhibits evidence of younger (∼3.4–4 Ga) silica-rich systems, identified through remote sensing (Glotch et al., 2010

Glotch, T.D., Lucey, P.G., Bandfield, J.L., Greenhagen, B.T., Thomas, I.R., Elphic, R.C., Bowles, N., Wyatt, M.B., Allen, C.C., Hanna, K.D., Paige, D.A. (2010) Highly Silicic Compositions on the Moon. Science 329, 1510–1513. https://doi.org/10.1126/science.1192148

; Siegler et al., 2023

Siegler, M.A., Feng, J., Lehman-Franco, K., Andrews-Hanna, J.C., Economos, R.C., St. Clair, M., Million, C., Head, J.W., Glotch, T.D., White, M.N. (2023) Remote detection of a lunar granitic batholith at Compton–Belkovich. Nature 620, 116–121. https://doi.org/10.1038/s41586-023-06183-5

), although detecting such features earlier in lunar history remains challenging. Apollo breccias contain ancient “granite” clasts with 4–4.3 Ga zircon (Meyer et al., 1996

Meyer, C., Williams, I.S., Compston, W. (1996) Uranium-lead ages for lunar zircons: Evidence for a prolonged period of granophyre formation from 4.32 to 3.88 Ga. Meteoritics & Planetary Science 31, 370–387. https://doi.org/10.1111/j.1945-5100.1996.tb02075.x

), but their complex history of crust formation (Su et al., 2023

Su, B., Chen, Y., Yue, Z., Chen, L., Mitchell, R.N., Tang, M., Yang, W., Huang, G., Guo, J., Li, X.-H., Wu, F.-Y. (2023) Crustal remelting origin of highly silicic magmatism on the Moon. Communications Earth & Environment 4, 228. https://doi.org/10.1038/s43247-023-00900-8

), meteorite impacts (Thiessen et al., 2018

Thiessen, F., Nemchin, A.A., Snape, J.F., Bellucci, J.J., Whitehouse, M.J. (2018) Apollo 12 breccia 12013: Impact-induced partial Pb loss in zircon and its implications for lunar geochronology. Geochimica et Cosmochimica Acta 230, 94–111. https://doi.org/10.1016/j.gca.2018.03.023

; Valencia et al., 2024

Valencia, S.N., Korotev, R.L., Jolliff, B.L. (2024) Compositional analysis of Apollo 12 Granitic Breccia 12013: Insights into protoliths and formation. Geochimica et Cosmochimica Acta 367, 189–212. https://doi.org/10.1016/j.gca.2023.12.034

) and disequilibrium assemblages (Seddio et al., 2013

Seddio, S.M., Jolliff, B.L., Korotev, R.L., Zeigler, R.A. (2013) Petrology and geochemistry of lunar granite 12032,366-19 and implications for lunar granite petrogenesis. American Mineralogist 98, 1697–1713. https://doi.org/10.2138/am.2013.4330

; Bellucci et al., 2019

Bellucci, J.J., Nemchin, A.A., Grange, M., Robinson, K.L., Collins, G., Whitehouse, M.J., Snape, J.F., Norman, M.D., Kring, D.A. (2019) Terrestrial-like zircon in a clast from an Apollo 14 breccia. Earth and Planetary Science Letters 510, 173–185. https://doi.org/10.1016/j.epsl.2019.01.010

) obscure primary chemistry. Large variations in accessory mineral ages within individual clasts (Grange et al., 2013

Grange, M.L., Pidgeon, R.T., Nemchin, A.A., Timms, N.E., Meyer, C. (2013) Interpreting U–Pb data from primary and secondary features in lunar zircon. Geochimica et Cosmochimica Acta 101, 112–132. https://doi.org/10.1016/j.gca.2012.10.013

; Thiessen et al., 2018

Thiessen, F., Nemchin, A.A., Snape, J.F., Bellucci, J.J., Whitehouse, M.J. (2018) Apollo 12 breccia 12013: Impact-induced partial Pb loss in zircon and its implications for lunar geochronology. Geochimica et Cosmochimica Acta 230, 94–111. https://doi.org/10.1016/j.gca.2018.03.023

) further complicate interpretations. To circumvent these challenges and better understand the early lunar crust’s composition, we need a tool that can probe the silica component of the system, independent of the samples’ complex histories. To accomplish this, we combine lunar zircon (ZrSiO4) ages with a new proxy for melt silica activity – a strong indicator of silicic magmatism – to characterise ancient crust missing from the lunar record. From first principles in thermodynamics, zircon can crystallise with a silica activity from ∼0.5, like a mafic or anorthositic system, to 1 in quartz-saturated silicic (granitic) environments (see Supplementary Information). We also calculate the Th content of the melts from which the zircons crystallised and find strong evidence that they were enriched in radioactive elements relative to the Procellarum KREEP Terrane. These data reveal persistent, highly silicic, Th-rich crust in the Apollo 14 lunar breccia record.

Our approach utilises X-ray absorption near edge structure spectroscopy (XANES) to investigate Ti coordination in zircon. The Ti content in zircon is often used as a tool to determine crystallisation temperature (Ferry and Watson, 2007

Ferry, J.M., Watson, E.B. (2007) New thermodynamic models and revised calibrations for the Ti-in-zircon and Zr-in-rutile thermometers. Contributions to Mineralogy and Petrology 154, 429–437. https://doi.org/10.1007/s00410-007-0201-0

; Trail et al., 2020

Trail, D., Barboni, M., McKeegan, K.D. (2020) Evidence for diverse lunar melt compositions and mixing of the pre-3.9 Ga crust from zircon chemistry. Geochimica et Cosmochimica Acta 284, 173–195. https://doi.org/10.1016/j.gca.2020.06.018

). A critical component of the thermodynamic theory that dictates the application of the thermometer is that Ti substitutes overwhelmingly for Si on the tetrahedral site in the zircon structure, as opposed to the Zr dodecahedral site, at crustal pressures (Tailby et al., 2011

Tailby, N.D., Walker, A.M., Berry, A.J., Hermann, J., Evans, K.A., Mavrogenes, J.A., O’Neill, H.St.C., Rodina, I.S., Soldatov, A.V., Rubatto, D., Sutton, S.R. (2011) Ti site occupancy in zircon. Geochimica et Cosmochimica Acta 75, 905–921. https://doi.org/10.1016/j.gca.2010.11.004

) (see Supplementary Information). We demonstrate that this long held view requires significant modification. We present experiments in which zircon co-crystallised with rutile but where the silica activity of each experimental system was buffered and calculable (see Methods and Supplementary Information). In one case, the silica activity was buffered by tridymite, akin to the crystallisation environment of a granite (aSiO2 ≈ 1, relative to quartz). The other experiment was undersaturated in SiO2 with a silica activity buffered at ∼0.6, representative of a crystallising mafic/anorthositic rock. Both experiments were conducted at 1300 °C and 1 atm (Table S-1).

top

Results

Abstract | Introduction | Results | Acknowledgements | References | Supplementary Information


We report 43 XANES spectra on randomly oriented zircons from these two experimental products to determine the fraction of Ti residing in the Si tetrahedral site (in 4-fold coordination) vs. Zr dodecahedral site (in 8-fold coordination) of zircon (see Methods). Contrary to expectation (Ferry and Watson, 2007

Ferry, J.M., Watson, E.B. (2007) New thermodynamic models and revised calibrations for the Ti-in-zircon and Zr-in-rutile thermometers. Contributions to Mineralogy and Petrology 154, 429–437. https://doi.org/10.1007/s00410-007-0201-0

; Tailby et al., 2011

Tailby, N.D., Walker, A.M., Berry, A.J., Hermann, J., Evans, K.A., Mavrogenes, J.A., O’Neill, H.St.C., Rodina, I.S., Soldatov, A.V., Rubatto, D., Sutton, S.R. (2011) Ti site occupancy in zircon. Geochimica et Cosmochimica Acta 75, 905–921. https://doi.org/10.1016/j.gca.2010.11.004

; Crisp et al., 2023

Crisp, L.J., Berry, A.J., Burnham, A.D., Miller, L.A., Newville, M. (2023) The Ti-in-zircon thermometer revised: The effect of pressure on the Ti site in zircon. Geochimica et Cosmochimica Acta 360, 241–258. https://doi.org/10.1016/j.gca.2023.04.031

), a significant fraction of the Ti was found to be in the zircon 8-fold site (Tables S-2, S-3). To further confirm our observations, we performed a separate experiment at different P/T conditions. While the previous experiments were conducted at 1300 °C/1 atm, lunar zircons likely crystallised at a lower T and moderately higher P (Crow et al., 2017

Crow, C.A., McKeegan, K.D., Moser, D.E. (2017) Coordinated U–Pb geochronology, trace element, Ti-in-zircon thermometry and microstructural analysis of Apollo zircons. Geochimica et Cosmochimica Acta 202, 264–284. https://doi.org/10.1016/j.gca.2016.12.019

; Bellucci et al., 2019

Bellucci, J.J., Nemchin, A.A., Grange, M., Robinson, K.L., Collins, G., Whitehouse, M.J., Snape, J.F., Norman, M.D., Kring, D.A. (2019) Terrestrial-like zircon in a clast from an Apollo 14 breccia. Earth and Planetary Science Letters 510, 173–185. https://doi.org/10.1016/j.epsl.2019.01.010

; Trail et al., 2020

Trail, D., Barboni, M., McKeegan, K.D. (2020) Evidence for diverse lunar melt compositions and mixing of the pre-3.9 Ga crust from zircon chemistry. Geochimica et Cosmochimica Acta 284, 173–195. https://doi.org/10.1016/j.gca.2020.06.018

). Our follow up experiment was performed at 800 °C and 1 GPa – i.e. at lower T and higher P than expected for lunar zircons (see Methods), with SiO2 and TiO2 activities buffered by quartz and rutile, respectively. These zircons from this experiment also yield a low fraction of Ti in the Si tetrahedral site (∼0.3). Considering these results alongside the other silica saturated experiment at 1300 °C and 1 atm, the differences in Ti site occupancy of zircon represent a drastic change from past studies, which interpreted ∼95 % of the Ti to be present in the Si site (Tailby et al., 2011

Tailby, N.D., Walker, A.M., Berry, A.J., Hermann, J., Evans, K.A., Mavrogenes, J.A., O’Neill, H.St.C., Rodina, I.S., Soldatov, A.V., Rubatto, D., Sutton, S.R. (2011) Ti site occupancy in zircon. Geochimica et Cosmochimica Acta 75, 905–921. https://doi.org/10.1016/j.gca.2010.11.004

; Crisp et al., 2023

Crisp, L.J., Berry, A.J., Burnham, A.D., Miller, L.A., Newville, M. (2023) The Ti-in-zircon thermometer revised: The effect of pressure on the Ti site in zircon. Geochimica et Cosmochimica Acta 360, 241–258. https://doi.org/10.1016/j.gca.2023.04.031

). These findings require a revision of the thermodynamic model for Ti in zircon and its application to thermometry (see Supplementary Information). However, they also present a significant opportunity, as the Ti coordination in zircon separated from the parent rock can serve as proxy for identifying high silica systems (Fig. 1).


Figure 1 XANES spectral analysis of Ti coordination in zircon. (a) Two standard spectra with Ti present in only 4-fold coordination (cristobalite); and a model spectrum for Ti in 8-fold coordination in zircon along the a axis (Tailby et al., 2011

Tailby, N.D., Walker, A.M., Berry, A.J., Hermann, J., Evans, K.A., Mavrogenes, J.A., O’Neill, H.St.C., Rodina, I.S., Soldatov, A.V., Rubatto, D., Sutton, S.R. (2011) Ti site occupancy in zircon. Geochimica et Cosmochimica Acta 75, 905–921. https://doi.org/10.1016/j.gca.2010.11.004

). (inset) Linear mixing of the pre-edge feature showing model spectra for different fractions of Ti in 4-fold coordination. (b) Example pre-edge features for experimental zircon for silica saturated and silica undersaturated (Tailby et al., 2011

Tailby, N.D., Walker, A.M., Berry, A.J., Hermann, J., Evans, K.A., Mavrogenes, J.A., O’Neill, H.St.C., Rodina, I.S., Soldatov, A.V., Rubatto, D., Sutton, S.R. (2011) Ti site occupancy in zircon. Geochimica et Cosmochimica Acta 75, 905–921. https://doi.org/10.1016/j.gca.2010.11.004

) experiments. The 1300 °C, 1 atm, SiO2 undersaturated experiment (srilankite (Ti,Zr)O4, rutile, zircon) has a pre-edge intensity similar to that of the 4-fold standard, whereas the intensity of the pre-edge feature for the two silica saturated experiments is muted meaning that a significant fraction of 8-fold Ti is present in zircon. The energy shift observed in the pre-edge spectra of the silica undersaturated samples is common, and has been shown to be related to crystal orientation (Tailby et al., 2011

Tailby, N.D., Walker, A.M., Berry, A.J., Hermann, J., Evans, K.A., Mavrogenes, J.A., O’Neill, H.St.C., Rodina, I.S., Soldatov, A.V., Rubatto, D., Sutton, S.R. (2011) Ti site occupancy in zircon. Geochimica et Cosmochimica Acta 75, 905–921. https://doi.org/10.1016/j.gca.2010.11.004

; Crisp et al., 2023

Crisp, L.J., Berry, A.J., Burnham, A.D., Miller, L.A., Newville, M. (2023) The Ti-in-zircon thermometer revised: The effect of pressure on the Ti site in zircon. Geochimica et Cosmochimica Acta 360, 241–258. https://doi.org/10.1016/j.gca.2023.04.031

). Differences in pre-edge intensity are due to orientation and/or sample heterogeneity (Tailby et al., 2011

Tailby, N.D., Walker, A.M., Berry, A.J., Hermann, J., Evans, K.A., Mavrogenes, J.A., O’Neill, H.St.C., Rodina, I.S., Soldatov, A.V., Rubatto, D., Sutton, S.R. (2011) Ti site occupancy in zircon. Geochimica et Cosmochimica Acta 75, 905–921. https://doi.org/10.1016/j.gca.2010.11.004

). Spectra were carefully evaluated for inclusions (see Supplementary Information and Fig. S-1). (c) Fraction of 4-fold Ti in zircon, each analysis from the three experiments calculated from a linear combination of the two standard spectra presented in (a) (Tailby et al., 2011

Tailby, N.D., Walker, A.M., Berry, A.J., Hermann, J., Evans, K.A., Mavrogenes, J.A., O’Neill, H.St.C., Rodina, I.S., Soldatov, A.V., Rubatto, D., Sutton, S.R. (2011) Ti site occupancy in zircon. Geochimica et Cosmochimica Acta 75, 905–921. https://doi.org/10.1016/j.gca.2010.11.004

; Crisp et al., 2023

Crisp, L.J., Berry, A.J., Burnham, A.D., Miller, L.A., Newville, M. (2023) The Ti-in-zircon thermometer revised: The effect of pressure on the Ti site in zircon. Geochimica et Cosmochimica Acta 360, 241–258. https://doi.org/10.1016/j.gca.2023.04.031

). When aqtzSiO2 ≈ 1, only ∼50 % of the Ti is occupying the 4-fold Si site, whereas with aqtzSiO2 ≈ 0.6 nearly 100 % of the Ti is in the 4-fold Si site. Error bars represent the 1 s.d. for each group. Experimental products can be found in Table S-1.
Full size image


The experimental observations are crucial because they can be applied to evaluate the silica activity in the igneous systems that crystallised zircon. Lunar zircons are exactly that: they are commonly found in breccias and soils as fragments of larger crystals yet they are high value analytical targets because U-Pb ages, typically with 10 Ma precision, can be measured in situ (Trail et al., 2020

Trail, D., Barboni, M., McKeegan, K.D. (2020) Evidence for diverse lunar melt compositions and mixing of the pre-3.9 Ga crust from zircon chemistry. Geochimica et Cosmochimica Acta 284, 173–195. https://doi.org/10.1016/j.gca.2020.06.018

). With age information already available for a suite of samples (Trail et al., 2020

Trail, D., Barboni, M., McKeegan, K.D. (2020) Evidence for diverse lunar melt compositions and mixing of the pre-3.9 Ga crust from zircon chemistry. Geochimica et Cosmochimica Acta 284, 173–195. https://doi.org/10.1016/j.gca.2020.06.018

), we collected Ti XANES spectra on zircons with ages spanning ∼400 Myr (Table S-4). The majority of zircons (n = 22) are from Apollo sample 14311,58 whose parent rock is a polymict impact melt breccia. In addition, zircons from Apollo 14 soil samples 14259 and 14163 were also analysed (n = 4).

To evaluate the sensitivity of grain orientation to spectral features, we collected XANES spectra for two different grain orientations, separated by 90° (Fig. 2). Some grains exhibit differences in absorption characteristics depending on the orientation whereas others do not (Fig. 2a,b). We attribute this to a combination of beam polarisation, grain orientation, and mineral anisotropy. The difference in the pre-edge intensity is, on average, 0.07 normalised units, about a factor of seven smaller than the effects caused by a change in the silica activity of the system. Thus, while grain orientation does change the pre-edge intensities, it is a second order effect (see Supplementary Information and Figs. S-2, S-3). We compare the fraction of Ti in 4-fold coordination in lunar zircon with that from 1300 °C experiments, which provide strong evidence for a high silica activity (ameltSiO2 ≈ 1) for zircon crystallisation in the lunar environment (Fig. 2c). This differs from previous thermodynamic calculations for lunar melt evolution, which used the MELTS algorithm and zircon saturation thermometry to estimate ameltSiO2 in igneous systems at zircon saturation ranging from 0.5 to 0.7 (Crow et al., 2017

Crow, C.A., McKeegan, K.D., Moser, D.E. (2017) Coordinated U–Pb geochronology, trace element, Ti-in-zircon thermometry and microstructural analysis of Apollo zircons. Geochimica et Cosmochimica Acta 202, 264–284. https://doi.org/10.1016/j.gca.2016.12.019

).


Figure 2 Ti XANES spectra and coordination in lunar zircon. (a, b) Two representative grain analyses showing XANES spectra before and after a 90° rotation. The top panel shows a significant change in pre-edge absorption intensity, while the lower panel, more typical, shows virtually no change. (c) Calculation of Ti in 4-fold coordination, with points averaging lunar zircon measurements in two crystal orientations, plotted against the previously reported 207Pb/206Pb age (Trail et al., 2020

Trail, D., Barboni, M., McKeegan, K.D. (2020) Evidence for diverse lunar melt compositions and mixing of the pre-3.9 Ga crust from zircon chemistry. Geochimica et Cosmochimica Acta 284, 173–195. https://doi.org/10.1016/j.gca.2020.06.018

). Dashed lines show the average fraction of 4-fold Ti at 1300 °C, and shaded yellow areas represent the 1 s.d. of the fraction in synthetic zircon.
Full size image


Pressure influences the Ti site occupancy in zircon at >2 GPa (Ferriss et al., 2008

Ferriss, E.D.A., Essene, E.J., Becker, U. (2008) Computational study of the effect of pressure on the Ti-in-zircon geothermometer. European Journal of Mineralogy 20, 745–755. https://doi.org/10.1127/0935-1221/2008/0020-1860

; Crisp et al., 2023

Crisp, L.J., Berry, A.J., Burnham, A.D., Miller, L.A., Newville, M. (2023) The Ti-in-zircon thermometer revised: The effect of pressure on the Ti site in zircon. Geochimica et Cosmochimica Acta 360, 241–258. https://doi.org/10.1016/j.gca.2023.04.031

). Although these pressures exceed typical lunar crustal formation conditions, many surface samples likely experienced transient pressures >2 GPa due to impact shock metamorphism (Thiessen et al., 2018

Thiessen, F., Nemchin, A.A., Snape, J.F., Bellucci, J.J., Whitehouse, M.J. (2018) Apollo 12 breccia 12013: Impact-induced partial Pb loss in zircon and its implications for lunar geochronology. Geochimica et Cosmochimica Acta 230, 94–111. https://doi.org/10.1016/j.gca.2018.03.023

; Valencia et al., 2024

Valencia, S.N., Korotev, R.L., Jolliff, B.L. (2024) Compositional analysis of Apollo 12 Granitic Breccia 12013: Insights into protoliths and formation. Geochimica et Cosmochimica Acta 367, 189–212. https://doi.org/10.1016/j.gca.2023.12.034

). To investigate the effect of shock on Ti coordination in zircon, we collected Ti XANES spectra on ZrSiO4 synthetically shocked to ∼30 GPa (Szumila et al., 2023

Szumila, I., Trail, D., Erickson, T., Simon, J.I., Wielicki, M.M., Lapen, T., Nakajima, M., Fries, M., Bell, E.A. (2023) Microstructural changes and Pb mobility during the zircon to reidite transformation: Implications for planetary impact chronology. American Mineralogist 108, 1516–1529. https://doi.org/10.2138/am-2022-8604

). The products include shocked zircon, amorphous ZrSiO4, and reidite (Szumila et al., 2023

Szumila, I., Trail, D., Erickson, T., Simon, J.I., Wielicki, M.M., Lapen, T., Nakajima, M., Fries, M., Bell, E.A. (2023) Microstructural changes and Pb mobility during the zircon to reidite transformation: Implications for planetary impact chronology. American Mineralogist 108, 1516–1529. https://doi.org/10.2138/am-2022-8604

) (a high pressure ZrSiO4 polymorph). The fraction of 4-fold coordinated Ti is equal to or slightly higher in the shocked ZrSiO4 products when compared to unshocked zircon starting material (Fig. S-4 and Table S-5). Results indicate either no impact shock modification of the Ti XANES spectral features, or shock modification that would shift the calculated silica activity to lower apparent values. Thus, there is no or very limited Ti site exchange between Si and Zr sites in zircon during high pressure shock events, confirming that lunar zircons in this study formed in highly silicic systems. This interpretation is further supported by the relatively limited evidence of impact shock twining among the zircons investigated here, as determined by electron backscatter diffraction (see Fig. S-3). With these experimental constraints in mind, we now explore the implications of high silica systems within the broader context of early lunar crustal evolution and, more generally, the evolution of planetary crusts.

Remote sensing has identified younger ∼3.4–4.0 Ga plutonic and volcanic silicic rock on the near side of the Moon (Glotch et al., 2010

Glotch, T.D., Lucey, P.G., Bandfield, J.L., Greenhagen, B.T., Thomas, I.R., Elphic, R.C., Bowles, N., Wyatt, M.B., Allen, C.C., Hanna, K.D., Paige, D.A. (2010) Highly Silicic Compositions on the Moon. Science 329, 1510–1513. https://doi.org/10.1126/science.1192148

), also characterised by high Th contents (∼10–55 ppm; Hagerty et al., 2006

Hagerty, J.J., Lawrence, D.J., Hawke, B.R., Vaniman, D.T., Elphic, R.C., Feldman, W.C. (2006) Refined thorium abundances for lunar red spots: Implications for evolved, nonmare volcanism on the Moon. Journal of Geophysical Research: Planets 111, E06002. https://doi.org/10.1029/2005JE002592

; Wilson et al., 2015

Wilson, J.T., Eke, V.R., Massey, R.J., Elphic, R.C., Jolliff, B.L., Lawrence, D.J., Llewellin, E.W., McElwaine, J.N., Teodoro, L.F.A. (2015) Evidence for explosive silicic volcanism on the Moon from the extended distribution of thorium near the Compton-Belkovich Volcanic Complex. Journal of Geophysical Research: Planets 120, 92–108. https://doi.org/10.1002/2014JE004719

) relative to the Procellarum KREEP Terrane (∼5 ppm; Jolliff et al., 2000

Jolliff, B.L., Gillis, J.J., Haskin, L.A., Korotev, R.L., Wieczorek, M.A. (2000) Major lunar crustal terranes: Surface expressions and crust-mantle origins. Journal of Geophysical Research: Planets 105, 4197–4216. https://doi.org/10.1029/1999JE001103

). The Compton–Bel’kovich Volcanic Complex represents a similar highly silicic Th-rich feature on the far side of the Moon, that erupted at ∼3.5 Ga and possibly earlier at ∼3.8 Ga, based on crater size frequency distributions (Jolliff et al., 2011

Jolliff, B.L., Wiseman, S.A., Lawrence, S.J., Tran, T.N., Robinson, M.S., Sato, H., Hawke, B.R., Scholten, F., Oberst, J., Hiesinger, H., van der Bogert, C.H., Greenhagen, B.T., Glotch, T.D., Paige, D.A. (2011) Non-mare silicic volcanism on the lunar farside at Compton–Belkovich. Nature Geoscience 4, 566–571. https://doi.org/10.1038/ngeo1212

; Wilson et al., 2015

Wilson, J.T., Eke, V.R., Massey, R.J., Elphic, R.C., Jolliff, B.L., Lawrence, D.J., Llewellin, E.W., McElwaine, J.N., Teodoro, L.F.A. (2015) Evidence for explosive silicic volcanism on the Moon from the extended distribution of thorium near the Compton-Belkovich Volcanic Complex. Journal of Geophysical Research: Planets 120, 92–108. https://doi.org/10.1002/2014JE004719

; Siegler et al., 2023

Siegler, M.A., Feng, J., Lehman-Franco, K., Andrews-Hanna, J.C., Economos, R.C., St. Clair, M., Million, C., Head, J.W., Glotch, T.D., White, M.N. (2023) Remote detection of a lunar granitic batholith at Compton–Belkovich. Nature 620, 116–121. https://doi.org/10.1038/s41586-023-06183-5

). The lunar zircons from highly silicic systems identified here are older, with U-Pb ages ranging from 3.93 to 4.27 Ga. To compare these older zircons with the younger silicic systems mentioned above, we apply zircon-melt Th partition coefficients to calculate Th melt contents of ∼5–40 ppm for the 3.93 to 4.27 Ga zircon saturated systems. These values are in broad agreement with the younger silicic systems but extend the record of highly silicic, Th-enriched magmatism back to the early Moon (Fig. 3).


Figure 3 Th contents of highly silicic lunar systems. Circles represent calculated Th contents of igneous systems derived from zircon-melt Th experimental partitioning (see Supplementary Information; Fig. S-5 and Table S-5). Yellow squares represent Th contents for younger silicic systems investigated and quantified by remote sensing (Hagerty et al., 2006

Hagerty, J.J., Lawrence, D.J., Hawke, B.R., Vaniman, D.T., Elphic, R.C., Feldman, W.C. (2006) Refined thorium abundances for lunar red spots: Implications for evolved, nonmare volcanism on the Moon. Journal of Geophysical Research: Planets 111, E06002. https://doi.org/10.1029/2005JE002592

; Wilson et al., 2015

Wilson, J.T., Eke, V.R., Massey, R.J., Elphic, R.C., Jolliff, B.L., Lawrence, D.J., Llewellin, E.W., McElwaine, J.N., Teodoro, L.F.A. (2015) Evidence for explosive silicic volcanism on the Moon from the extended distribution of thorium near the Compton-Belkovich Volcanic Complex. Journal of Geophysical Research: Planets 120, 92–108. https://doi.org/10.1002/2014JE004719

). The ages for the silicic systems identified by remote sensing data are estimated through crater counting. The average Th contents from the PKT mare terrain (4.9 ppm) and other mare (2.2 ppm) are also shown for comparison (Jolliff et al., 2000

Jolliff, B.L., Gillis, J.J., Haskin, L.A., Korotev, R.L., Wieczorek, M.A. (2000) Major lunar crustal terranes: Surface expressions and crust-mantle origins. Journal of Geophysical Research: Planets 105, 4197–4216. https://doi.org/10.1029/1999JE001103

). MMD1/2, first and second eruptions of Middle Mairan Dome; MTD, Mairan T Dome; GDδ/γ, Gruithuisen Dome (δ) and (γ); LM, Lassell Massif; HD, Hansteen Alpha; CBVC, Compton–Bel’kovich Volcanic Complex.
Full size image


The XANES data also allow us to evaluate whether these zircons record significant changes to the silica activities of the igneous systems during crystallisation. For instance, if zircons crystallised from the extreme fractionation of mafic melts, evolving silica activity should produce variable pre-edge intensities of the XANES spectra. However, for ∼90 % of the analysed lunar grains, the pre-edge intensity is within 0.05 normalised units across different analytical locations on each grain (see Supplementary Information). We interpret this as evidence that most zircons crystallised from high silica igneous bodies with approximately constant silica activity (e.g., quartz saturated), rather than from fractionating mafic systems. A few exceptional grains show large pre-edge intensity variations of ∼0.3, which we interpret to represent evolving silica conditions along the liquid line of descent (Figs. S-6, S-7). But overall, the zircon data represent highly silicic systems that are under-represented in sample collections from Apollo, Luna, meteorites, and Chang’e 5 missions (Neal, 2009

Neal, C.R. (2009) The Moon 35 years after Apollo: What’s left to learn? Geochemistry 69, 3–43. https://doi.org/10.1016/j.chemer.2008.07.002

; Wu et al., 2024

Wu, F.-Y., Li, Q.-L., Chen, Y., Hu, S., Yue, Z.-Y., Zhou, Q., Wang, H., Yang, W., Tian, H.-C., Zhang, C., Li, J.-H., Li, L.-X., Hui, H.-J., Li, C.-L., Lin, Y.-T., Li, X.-H., Delano, J.W. (2024) Lunar Evolution in Light of the Chang’e-5 Returned Samples. Annual Review of Earth and Planetary Sciences 52, 159–194. https://doi.org/10.1146/annurev-earth-040722-100453

).

Next, we examine models that could lead to high silica activity in lunar systems, including felsic mesostasis, liquid immiscibility, and basaltic underplating of the crust accompanied by partial melting. Based on textural evidence, phase assemblage, and partition coefficients that conflict with our zircon data, we consider felsic mesostasis and liquid immiscibility between mafic and felsic systems to be unlikely (Seddio et al., 2013

Seddio, S.M., Jolliff, B.L., Korotev, R.L., Zeigler, R.A. (2013) Petrology and geochemistry of lunar granite 12032,366-19 and implications for lunar granite petrogenesis. American Mineralogist 98, 1697–1713. https://doi.org/10.2138/am.2013.4330

; Gullikson et al., 2016

Gullikson, A.L., Hagerty, J.J., Reid, M.R., Rapp, J.F., Draper, D.S. (2016) Silicic lunar volcanism: Testing the crustal melting model. American Mineralogist 101, 2312–2321. https://doi.org/10.2138/am-2016-5619

) (see Supplementary Information). Instead, we favour a model similar to that ascribed to the origin of the red spot features (Hagerty et al., 2006

Hagerty, J.J., Lawrence, D.J., Hawke, B.R., Vaniman, D.T., Elphic, R.C., Feldman, W.C. (2006) Refined thorium abundances for lunar red spots: Implications for evolved, nonmare volcanism on the Moon. Journal of Geophysical Research: Planets 111, E06002. https://doi.org/10.1029/2005JE002592

), in which basaltic underplating and partial melting of pre-existing lunar crust produces highly silicic, high Th lunar melts (Hagerty et al., 2006

Hagerty, J.J., Lawrence, D.J., Hawke, B.R., Vaniman, D.T., Elphic, R.C., Feldman, W.C. (2006) Refined thorium abundances for lunar red spots: Implications for evolved, nonmare volcanism on the Moon. Journal of Geophysical Research: Planets 111, E06002. https://doi.org/10.1029/2005JE002592

; Glotch et al., 2010

Glotch, T.D., Lucey, P.G., Bandfield, J.L., Greenhagen, B.T., Thomas, I.R., Elphic, R.C., Bowles, N., Wyatt, M.B., Allen, C.C., Hanna, K.D., Paige, D.A. (2010) Highly Silicic Compositions on the Moon. Science 329, 1510–1513. https://doi.org/10.1126/science.1192148

; Gullikson et al., 2016

Gullikson, A.L., Hagerty, J.J., Reid, M.R., Rapp, J.F., Draper, D.S. (2016) Silicic lunar volcanism: Testing the crustal melting model. American Mineralogist 101, 2312–2321. https://doi.org/10.2138/am-2016-5619

; Valencia et al., 2024

Valencia, S.N., Korotev, R.L., Jolliff, B.L. (2024) Compositional analysis of Apollo 12 Granitic Breccia 12013: Insights into protoliths and formation. Geochimica et Cosmochimica Acta 367, 189–212. https://doi.org/10.1016/j.gca.2023.12.034

) which is consistent with our data (Fig. 4).


Figure 4 Simplified illustration (not to scale) of the lunar silicic systems whose formation was initiated by partial melting of fertile source rocks such as a KREEP basalt by basaltic underplating (Hagerty et al., 2006

Hagerty, J.J., Lawrence, D.J., Hawke, B.R., Vaniman, D.T., Elphic, R.C., Feldman, W.C. (2006) Refined thorium abundances for lunar red spots: Implications for evolved, nonmare volcanism on the Moon. Journal of Geophysical Research: Planets 111, E06002. https://doi.org/10.1029/2005JE002592

). High viscosity minimum melt (?) silicic systems may have developed not far from final emplacement location and/or made significant use of the higher porosity of the fractured crust/mega-regolith, including extensional features associated with impact basins (Valencia et al., 2024

Valencia, S.N., Korotev, R.L., Jolliff, B.L. (2024) Compositional analysis of Apollo 12 Granitic Breccia 12013: Insights into protoliths and formation. Geochimica et Cosmochimica Acta 367, 189–212. https://doi.org/10.1016/j.gca.2023.12.034

). Early silica-rich domes were likely destroyed by large basin-forming impacts, or covered by impact ejecta (Qiu et al., 2023

Qiu, D., Sasaki, S., Yan, J., Ye, M., Deng, Q., Liang, F., Liu, L., Li, F. (2023) Buried Silicic Volcanoes Discovered in the Gruithuisen Region on the Moon. Geophysical Research Letters 50, e2023GL103336. https://doi.org/10.1029/2023GL103336

).
Full size image


In summary, the high silica activities and Th contents recorded by these zircons, provide compelling evidence for production of a silica-rich component of the lunar crust from at least 4.27 to 3.4 Ga. We propose that these zircons represent vestiges of early highly silicic rocks that have been largely obliterated by basin-forming impacts, buried in impact ejecta, or covered by later mare basalts. Similar processes are observed in younger silicic localities such as the Gruithuisen region (Qiu et al., 2023

Qiu, D., Sasaki, S., Yan, J., Ye, M., Deng, Q., Liang, F., Liu, L., Li, F. (2023) Buried Silicic Volcanoes Discovered in the Gruithuisen Region on the Moon. Geophysical Research Letters 50, e2023GL103336. https://doi.org/10.1029/2023GL103336

), Aristarchus, and Lassell regions, where silicic lithologies are exposed through later impact cratering (Glotch et al., 2010

Glotch, T.D., Lucey, P.G., Bandfield, J.L., Greenhagen, B.T., Thomas, I.R., Elphic, R.C., Bowles, N., Wyatt, M.B., Allen, C.C., Hanna, K.D., Paige, D.A. (2010) Highly Silicic Compositions on the Moon. Science 329, 1510–1513. https://doi.org/10.1126/science.1192148

; Clegg-Watkins et al., 2017

Clegg-Watkins, R.N., Jolliff, B.L., Watkins, M.J., Coman, E., Giguere, T.A., Stopar, J.D., Lawrence, S.J. (2017) Nonmare volcanism on the Moon: Photometric evidence for the presence of evolved silicic materials. Icarus 285, 169–184. https://doi.org/10.1016/j.icarus.2016.12.004

; Yang et al., 2023

Yang, J., Ju, D., Pang, R., Li, R., Liu, J., Du, W. (2023) Significance of silicate liquid immiscibility for the origin of young highly evolved lithic clasts in Chang’E-5 regolith. Geochimica et Cosmochimica Acta 340, 189–205. https://doi.org/10.1016/j.gca.2022.11.008

). This aligns with the occurrence of 1–3 mm felsic clasts found in lunar breccia (Meyer et al., 1996

Meyer, C., Williams, I.S., Compston, W. (1996) Uranium-lead ages for lunar zircons: Evidence for a prolonged period of granophyre formation from 4.32 to 3.88 Ga. Meteoritics & Planetary Science 31, 370–387. https://doi.org/10.1111/j.1945-5100.1996.tb02075.x

; Seddio et al., 2013

Seddio, S.M., Jolliff, B.L., Korotev, R.L., Zeigler, R.A. (2013) Petrology and geochemistry of lunar granite 12032,366-19 and implications for lunar granite petrogenesis. American Mineralogist 98, 1697–1713. https://doi.org/10.2138/am.2013.4330

) containing impact-generated glass and brecciated material (Warren et al., 1983

Warren, P.H., Taylor, G.J., Keil, K., Shirley, D.N., Wasson, J.T. (1983) Petrology and chemistry of two “large” granite clasts from the moon. Earth and Planetary Science Letters 64, 175–185. https://doi.org/10.1016/0012-821X(83)90202-9

; Jolliff, 1991

Jolliff, B.L. (1991) Fragments of quartz monzodiorite and felsite in Apollo 14 soil particles. Proceedings of Lunar and Planetary Science 21, 101–118.

).

Our zircon record demonstrates that, despite the broadly mafic/anorthositic character of the lunar crust, silicic lithologies were an integral and persistent component of the early Moon. Crucially, this variability is not confined to a specific period in lunar history, nor do we uncover clear evidence for distinct intervals of silicic melt initiation and cessation (Meyer et al., 1996

Meyer, C., Williams, I.S., Compston, W. (1996) Uranium-lead ages for lunar zircons: Evidence for a prolonged period of granophyre formation from 4.32 to 3.88 Ga. Meteoritics & Planetary Science 31, 370–387. https://doi.org/10.1111/j.1945-5100.1996.tb02075.x

). Further confirmation will require additional measurements similar to those presented here for time intervals not represented including the pre-4.3 Ga lunar zircon record. We conclude that even small single plate planetary bodies acquire a heterogeneous crust of mafic and highly silicic rocks, and that the ongoing production of silicic systems is likely a common trait of planetary bodies, even in the first billion years of development.

top

Acknowledgements

Abstract | Introduction | Results | Acknowledgements | References | Supplementary Information


We are grateful for the constructive and thorough reviews of Timmons Erickson and Steve Parman. We thank Hugh O’Neill for advice and input. This work was funded by the NASA Solar System Workings grant 80NSSC20K1039. We thank Axel Wittman for assistance with EBSD analyses. Portions of this work were performed at GeoSoilEnviroCARS (The University of Chicago, Sector 13), Advanced Photon Source (APS), Argonne National Laboratory. GeoSoilEnviroCARS is supported by the National Science Foundation – Earth Sciences (EAR – 1634415) and NASA’s Planetary Science Enabling Facilities (PSEF) program (grant number 80NSSC23K0196). This research used resources of the Advanced Photon Source, a U.S. Department of Energy (DOE) Office of Science User Facility operated for the DOE Office of Science by Argonne National Laboratory under Contract No. DE-AC02-06CH11357.

Editor: Francis McCubbin

top

References

Abstract | Introduction | Results | Acknowledgements | References | Supplementary Information

Bellucci, J.J., Nemchin, A.A., Grange, M., Robinson, K.L., Collins, G., Whitehouse, M.J., Snape, J.F., Norman, M.D., Kring, D.A. (2019) Terrestrial-like zircon in a clast from an Apollo 14 breccia. Earth and Planetary Science Letters 510, 173–185. https://doi.org/10.1016/j.epsl.2019.01.010
Show in context

The Moon exhibits evidence of younger (∼3.4–4 Ga) silica-rich systems, identified through remote sensing (Glotch et al., 2010; Siegler et al., 2023), although detecting such features earlier in lunar history remains challenging. Apollo breccias contain ancient “granite” clasts with 4–4.3 Ga zircon (Meyer et al., 1996), but their complex history of crust formation (Su et al., 2023), meteorite impacts (Thiessen et al., 2018; Valencia et al., 2024) and disequilibrium assemblages (Seddio et al., 2013; Bellucci et al., 2019) obscure primary chemistry. Large variations in accessory mineral ages within individual clasts (Grange et al., 2013; Thiessen et al., 2018) further complicate interpretations.
View in article
To further confirm our observations, we performed a separate experiment at different P/T conditions. While the previous experiments were conducted at 1300 °C/1 atm, lunar zircons likely crystallised at a lower T and moderately higher P (Crow et al., 2017; Bellucci et al., 2019; Trail et al., 2020).
View in article


Cawood, P.A., Chowdhury, P., Mulder, J.A., Hawkesworth, C.J., Capitanio, F.A., Gunawardana, P.M., Nebel, O. (2022) Secular Evolution of Continents and the Earth System. Reviews of Geophysics 60, e2022RG000789. https://doi.org/10.1029/2022RG000789
Show in context

Some models infer that there was continuous production and recycling of early high silica crust (Harrison, 2020) while others favour the absence of a stable volume of buoyant continental material before 4 Ga (Cawood et al., 2022).
View in article


Clegg-Watkins, R.N., Jolliff, B.L., Watkins, M.J., Coman, E., Giguere, T.A., Stopar, J.D., Lawrence, S.J. (2017) Nonmare volcanism on the Moon: Photometric evidence for the presence of evolved silicic materials. Icarus 285, 169–184. https://doi.org/10.1016/j.icarus.2016.12.004
Show in context

Similar processes are observed in younger silicic localities such as the Gruithuisen region (Qiu et al., 2023), Aristarchus, and Lassell regions, where silicic lithologies are exposed through later impact cratering (Glotch et al., 2010; Clegg-Watkins et al., 2017; Yang et al., 2023).
View in article


Crisp, L.J., Berry, A.J., Burnham, A.D., Miller, L.A., Newville, M. (2023) The Ti-in-zircon thermometer revised: The effect of pressure on the Ti site in zircon. Geochimica et Cosmochimica Acta 360, 241–258. https://doi.org/10.1016/j.gca.2023.04.031
Show in context

We report 43 XANES spectra on randomly oriented zircons from these two experimental products to determine the fraction of Ti residing in the Si tetrahedral site (in 4-fold coordination) vs. Zr dodecahedral site (in 8-fold coordination) of zircon (see Methods). Contrary to expectation (Ferry and Watson, 2007; Tailby et al., 2011; Crisp et al., 2023), a significant fraction of the Ti was found to be in the zircon 8-fold site (Tables S-2, S-3).
View in article
Considering these results alongside the other silica saturated experiment at 1300 °C and 1 atm, the differences in Ti site occupancy of zircon represent a drastic change from past studies, which interpreted ∼95 % of the Ti to be present in the Si site (Tailby et al., 2011; Crisp et al., 2023).
View in article
The energy shift observed in the pre-edge spectra of the silica undersaturated samples is common, and has been shown to be related to crystal orientation (Tailby et al., 2011; Crisp et al., 2023).
View in article
(c) Fraction of 4-fold Ti in zircon, each analysis from the three experiments calculated from a linear combination of the two standard spectra presented in (a) (Tailby et al., 2011; Crisp et al., 2023).
View in article
Pressure influences the Ti site occupancy in zircon at >2 GPa (Ferriss et al., 2008; Crisp et al., 2023).
View in article


Crow, C.A., McKeegan, K.D., Moser, D.E. (2017) Coordinated U–Pb geochronology, trace element, Ti-in-zircon thermometry and microstructural analysis of Apollo zircons. Geochimica et Cosmochimica Acta 202, 264–284. https://doi.org/10.1016/j.gca.2016.12.019
Show in context

To further confirm our observations, we performed a separate experiment at different P/T conditions. While the previous experiments were conducted at 1300 °C/1 atm, lunar zircons likely crystallised at a lower T and moderately higher P (Crow et al., 2017; Bellucci et al., 2019; Trail et al., 2020).
View in article
This differs from previous thermodynamic calculations for lunar melt evolution, which used the MELTS algorithm and zircon saturation thermometry to estimate ameltSiO2 in igneous systems at zircon saturation ranging from 0.5 to 0.7 (Crow et al., 2017).
View in article


Ferriss, E.D.A., Essene, E.J., Becker, U. (2008) Computational study of the effect of pressure on the Ti-in-zircon geothermometer. European Journal of Mineralogy 20, 745–755. https://doi.org/10.1127/0935-1221/2008/0020-1860
Show in context

Pressure influences the Ti site occupancy in zircon at >2 GPa (Ferriss et al., 2008; Crisp et al., 2023).
View in article


Ferry, J.M., Watson, E.B. (2007) New thermodynamic models and revised calibrations for the Ti-in-zircon and Zr-in-rutile thermometers. Contributions to Mineralogy and Petrology 154, 429–437. https://doi.org/10.1007/s00410-007-0201-0
Show in context

The Ti content in zircon is often used as a tool to determine crystallisation temperature (Ferry and Watson, 2007; Trail et al., 2020).
View in article
We report 43 XANES spectra on randomly oriented zircons from these two experimental products to determine the fraction of Ti residing in the Si tetrahedral site (in 4-fold coordination) vs. Zr dodecahedral site (in 8-fold coordination) of zircon (see Methods). Contrary to expectation (Ferry and Watson, 2007; Tailby et al., 2011; Crisp et al., 2023), a significant fraction of the Ti was found to be in the zircon 8-fold site (Tables S-2, S-3).
View in article


Glotch, T.D., Lucey, P.G., Bandfield, J.L., Greenhagen, B.T., Thomas, I.R., Elphic, R.C., Bowles, N., Wyatt, M.B., Allen, C.C., Hanna, K.D., Paige, D.A. (2010) Highly Silicic Compositions on the Moon. Science 329, 1510–1513. https://doi.org/10.1126/science.1192148
Show in context

The Moon exhibits evidence of younger (∼3.4–4 Ga) silica-rich systems, identified through remote sensing (Glotch et al., 2010; Siegler et al., 2023), although detecting such features earlier in lunar history remains challenging. Apollo breccias contain ancient “granite” clasts with 4–4.3 Ga zircon (Meyer et al., 1996), but their complex history of crust formation (Su et al., 2023), meteorite impacts (Thiessen et al., 2018; Valencia et al., 2024) and disequilibrium assemblages (Seddio et al., 2013; Bellucci et al., 2019) obscure primary chemistry. Large variations in accessory mineral ages within individual clasts (Grange et al., 2013; Thiessen et al., 2018) further complicate interpretations.
View in article
Remote sensing has identified younger ∼3.4–4.0 Ga plutonic and volcanic silicic rock on the near side of the Moon (Glotch et al., 2010), also characterised by high Th contents (∼10–55 ppm; Hagerty et al., 2006; Wilson et al., 2015) relative to the Procellarum KREEP Terrane (∼5 ppm; Jolliff et al., 2000).
View in article
Instead, we favour a model similar to that ascribed to the origin of the red spot features (Hagerty et al., 2006), in which basaltic underplating and partial melting of pre-existing lunar crust produces highly silicic, high Th lunar melts (Hagerty et al., 2006; Glotch et al., 2010; Gullikson et al., 2016; Valencia et al., 2024) which is consistent with our data (Fig. 4).
View in article
Similar processes are observed in younger silicic localities such as the Gruithuisen region (Qiu et al., 2023), Aristarchus, and Lassell regions, where silicic lithologies are exposed through later impact cratering (Glotch et al., 2010; Clegg-Watkins et al., 2017; Yang et al., 2023).
View in article


Grange, M.L., Pidgeon, R.T., Nemchin, A.A., Timms, N.E., Meyer, C. (2013) Interpreting U–Pb data from primary and secondary features in lunar zircon. Geochimica et Cosmochimica Acta 101, 112–132. https://doi.org/10.1016/j.gca.2012.10.013
Show in context

The Moon exhibits evidence of younger (∼3.4–4 Ga) silica-rich systems, identified through remote sensing (Glotch et al., 2010; Siegler et al., 2023), although detecting such features earlier in lunar history remains challenging. Apollo breccias contain ancient “granite” clasts with 4–4.3 Ga zircon (Meyer et al., 1996), but their complex history of crust formation (Su et al., 2023), meteorite impacts (Thiessen et al., 2018; Valencia et al., 2024) and disequilibrium assemblages (Seddio et al., 2013; Bellucci et al., 2019) obscure primary chemistry. Large variations in accessory mineral ages within individual clasts (Grange et al., 2013; Thiessen et al., 2018) further complicate interpretations.
View in article


Gullikson, A.L., Hagerty, J.J., Reid, M.R., Rapp, J.F., Draper, D.S. (2016) Silicic lunar volcanism: Testing the crustal melting model. American Mineralogist 101, 2312–2321. https://doi.org/10.2138/am-2016-5619
Show in context

Next, we examine models that could lead to high silica activity in lunar systems, including felsic mesostasis, liquid immiscibility, and basaltic underplating of the crust accompanied by partial melting. Based on textural evidence, phase assemblage, and partition coefficients that conflict with our zircon data, we consider felsic mesostasis and liquid immiscibility between mafic and felsic systems to be unlikely (Seddio et al., 2013; Gullikson et al., 2016) (see Supplementary Information).
View in article
Instead, we favour a model similar to that ascribed to the origin of the red spot features (Hagerty et al., 2006), in which basaltic underplating and partial melting of pre-existing lunar crust produces highly silicic, high Th lunar melts (Hagerty et al., 2006; Glotch et al., 2010; Gullikson et al., 2016; Valencia et al., 2024) which is consistent with our data (Fig. 4).
View in article


Hagerty, J.J., Lawrence, D.J., Hawke, B.R., Vaniman, D.T., Elphic, R.C., Feldman, W.C. (2006) Refined thorium abundances for lunar red spots: Implications for evolved, nonmare volcanism on the Moon. Journal of Geophysical Research: Planets 111, E06002. https://doi.org/10.1029/2005JE002592
Show in context

Remote sensing has identified younger ∼3.4–4.0 Ga plutonic and volcanic silicic rock on the near side of the Moon (Glotch et al., 2010), also characterised by high Th contents (∼10–55 ppm; Hagerty et al., 2006; Wilson et al., 2015) relative to the Procellarum KREEP Terrane (∼5 ppm; Jolliff et al., 2000).
View in article
Yellow squares represent Th contents for younger silicic systems investigated and quantified by remote sensing (Hagerty et al., 2006; Wilson et al., 2015).
View in article
Instead, we favour a model similar to that ascribed to the origin of the red spot features (Hagerty et al., 2006), in which basaltic underplating and partial melting of pre-existing lunar crust produces highly silicic, high Th lunar melts (Hagerty et al., 2006; Glotch et al., 2010; Gullikson et al., 2016; Valencia et al., 2024) which is consistent with our data (Fig. 4).
View in article
Simplified illustration (not to scale) of the lunar silicic systems whose formation was initiated by partial melting of fertile source rocks such as a KREEP basalt by basaltic underplating (Hagerty et al., 2006).
View in article


Harrison, T.M. (2020) Hadean Earth. Springer, Cham. https://doi.org/10.1007/978-3-030-46687-9
Show in context

The emergence of buoyant, silica-rich (>60 wt. %) crust on a rocky planetary body is a key marker of crustal reprocessing and/or extreme differentiation with crucial implications for a planet’s evolution and possibly the origin of life (Harrison, 2020).
View in article
Some models infer that there was continuous production and recycling of early high silica crust (Harrison, 2020) while others favour the absence of a stable volume of buoyant continental material before 4 Ga (Cawood et al., 2022).
View in article


Jolliff, B.L. (1991) Fragments of quartz monzodiorite and felsite in Apollo 14 soil particles. Proceedings of Lunar and Planetary Science 21, 101–118.
Show in context

This aligns with the occurrence of 1–3 mm felsic clasts found in lunar breccia (Meyer et al., 1996; Seddio et al., 2013) containing impact-generated glass and brecciated material (Warren et al., 1983; Jolliff, 1991).
View in article


Jolliff, B.L., Gillis, J.J., Haskin, L.A., Korotev, R.L., Wieczorek, M.A. (2000) Major lunar crustal terranes: Surface expressions and crust-mantle origins. Journal of Geophysical Research: Planets 105, 4197–4216. https://doi.org/10.1029/1999JE001103
Show in context

Remote sensing has identified younger ∼3.4–4.0 Ga plutonic and volcanic silicic rock on the near side of the Moon (Glotch et al., 2010), also characterised by high Th contents (∼10–55 ppm; Hagerty et al., 2006; Wilson et al., 2015) relative to the Procellarum KREEP Terrane (∼5 ppm; Jolliff et al., 2000).
View in article
The ages for the silicic systems identified by remote sensing data are estimated through crater counting. The average Th contents from the PKT mare terrain (4.9 ppm) and other mare (2.2 ppm) are also shown for comparison (Jolliff et al., 2000).
View in article


Jolliff, B.L., Wiseman, S.A., Lawrence, S.J., Tran, T.N., Robinson, M.S., Sato, H., Hawke, B.R., Scholten, F., Oberst, J., Hiesinger, H., van der Bogert, C.H., Greenhagen, B.T., Glotch, T.D., Paige, D.A. (2011) Non-mare silicic volcanism on the lunar farside at Compton–Belkovich. Nature Geoscience 4, 566–571. https://doi.org/10.1038/ngeo1212
Show in context

The Compton–Bel’kovich Volcanic Complex represents a similar highly silicic Th-rich feature on the far side of the Moon, that erupted at ∼3.5 Ga and possibly earlier at ∼3.8 Ga, based on crater size frequency distributions (Jolliff et al., 2011; Wilson et al., 2015; Siegler et al., 2023).
View in article


Meyer, C., Williams, I.S., Compston, W. (1996) Uranium-lead ages for lunar zircons: Evidence for a prolonged period of granophyre formation from 4.32 to 3.88 Ga. Meteoritics & Planetary Science 31, 370–387. https://doi.org/10.1111/j.1945-5100.1996.tb02075.x
Show in context

The Moon exhibits evidence of younger (∼3.4–4 Ga) silica-rich systems, identified through remote sensing (Glotch et al., 2010; Siegler et al., 2023), although detecting such features earlier in lunar history remains challenging. Apollo breccias contain ancient “granite” clasts with 4–4.3 Ga zircon (Meyer et al., 1996), but their complex history of crust formation (Su et al., 2023), meteorite impacts (Thiessen et al., 2018; Valencia et al., 2024) and disequilibrium assemblages (Seddio et al., 2013; Bellucci et al., 2019) obscure primary chemistry. Large variations in accessory mineral ages within individual clasts (Grange et al., 2013; Thiessen et al., 2018) further complicate interpretations.
View in article
This aligns with the occurrence of 1–3 mm felsic clasts found in lunar breccia (Meyer et al., 1996; Seddio et al., 2013) containing impact-generated glass and brecciated material (Warren et al., 1983; Jolliff, 1991).
View in article
Crucially, this variability is not confined to a specific period in lunar history, nor do we uncover clear evidence for distinct intervals of silicic melt initiation and cessation (Meyer et al., 1996).
View in article


Neal, C.R. (2009) The Moon 35 years after Apollo: What’s left to learn? Geochemistry 69, 3–43. https://doi.org/10.1016/j.chemer.2008.07.002
Show in context

But overall, the zircon data represent highly silicic systems that are under-represented in sample collections from Apollo, Luna, meteorites, and Chang’e 5 missions (Neal, 2009; Wu et al., 2024).
View in article


Qiu, D., Sasaki, S., Yan, J., Ye, M., Deng, Q., Liang, F., Liu, L., Li, F. (2023) Buried Silicic Volcanoes Discovered in the Gruithuisen Region on the Moon. Geophysical Research Letters 50, e2023GL103336. https://doi.org/10.1029/2023GL103336
Show in context

Early silica-rich domes were likely destroyed by large basin-forming impacts, or covered by impact ejecta (Qiu et al., 2023).
View in article
Similar processes are observed in younger silicic localities such as the Gruithuisen region (Qiu et al., 2023), Aristarchus, and Lassell regions, where silicic lithologies are exposed through later impact cratering (Glotch et al., 2010; Clegg-Watkins et al., 2017; Yang et al., 2023).
View in article


Seddio, S.M., Jolliff, B.L., Korotev, R.L., Zeigler, R.A. (2013) Petrology and geochemistry of lunar granite 12032,366-19 and implications for lunar granite petrogenesis. American Mineralogist 98, 1697–1713. https://doi.org/10.2138/am.2013.4330
Show in context

The Moon exhibits evidence of younger (∼3.4–4 Ga) silica-rich systems, identified through remote sensing (Glotch et al., 2010; Siegler et al., 2023), although detecting such features earlier in lunar history remains challenging. Apollo breccias contain ancient “granite” clasts with 4–4.3 Ga zircon (Meyer et al., 1996), but their complex history of crust formation (Su et al., 2023), meteorite impacts (Thiessen et al., 2018; Valencia et al., 2024) and disequilibrium assemblages (Seddio et al., 2013; Bellucci et al., 2019) obscure primary chemistry. Large variations in accessory mineral ages within individual clasts (Grange et al., 2013; Thiessen et al., 2018) further complicate interpretations.
View in article
Next, we examine models that could lead to high silica activity in lunar systems, including felsic mesostasis, liquid immiscibility, and basaltic underplating of the crust accompanied by partial melting. Based on textural evidence, phase assemblage, and partition coefficients that conflict with our zircon data, we consider felsic mesostasis and liquid immiscibility between mafic and felsic systems to be unlikely (Seddio et al., 2013; Gullikson et al., 2016) (see Supplementary Information).
View in article
This aligns with the occurrence of 1–3 mm felsic clasts found in lunar breccia (Meyer et al., 1996; Seddio et al., 2013) containing impact-generated glass and brecciated material (Warren et al., 1983; Jolliff, 1991).
View in article


Siegler, M.A., Feng, J., Lehman-Franco, K., Andrews-Hanna, J.C., Economos, R.C., St. Clair, M., Million, C., Head, J.W., Glotch, T.D., White, M.N. (2023) Remote detection of a lunar granitic batholith at Compton–Belkovich. Nature 620, 116–121. https://doi.org/10.1038/s41586-023-06183-5
Show in context

The Moon exhibits evidence of younger (∼3.4–4 Ga) silica-rich systems, identified through remote sensing (Glotch et al., 2010; Siegler et al., 2023), although detecting such features earlier in lunar history remains challenging. Apollo breccias contain ancient “granite” clasts with 4–4.3 Ga zircon (Meyer et al., 1996), but their complex history of crust formation (Su et al., 2023), meteorite impacts (Thiessen et al., 2018; Valencia et al., 2024) and disequilibrium assemblages (Seddio et al., 2013; Bellucci et al., 2019) obscure primary chemistry. Large variations in accessory mineral ages within individual clasts (Grange et al., 2013; Thiessen et al., 2018) further complicate interpretations.
View in article
The Compton–Bel’kovich Volcanic Complex represents a similar highly silicic Th-rich feature on the far side of the Moon, that erupted at ∼3.5 Ga and possibly earlier at ∼3.8 Ga, based on crater size frequency distributions (Jolliff et al., 2011; Wilson et al., 2015; Siegler et al., 2023).
View in article


Su, B., Chen, Y., Yue, Z., Chen, L., Mitchell, R.N., Tang, M., Yang, W., Huang, G., Guo, J., Li, X.-H., Wu, F.-Y. (2023) Crustal remelting origin of highly silicic magmatism on the Moon. Communications Earth & Environment 4, 228. https://doi.org/10.1038/s43247-023-00900-8
Show in context

The Moon exhibits evidence of younger (∼3.4–4 Ga) silica-rich systems, identified through remote sensing (Glotch et al., 2010; Siegler et al., 2023), although detecting such features earlier in lunar history remains challenging. Apollo breccias contain ancient “granite” clasts with 4–4.3 Ga zircon (Meyer et al., 1996), but their complex history of crust formation (Su et al., 2023), meteorite impacts (Thiessen et al., 2018; Valencia et al., 2024) and disequilibrium assemblages (Seddio et al., 2013; Bellucci et al., 2019) obscure primary chemistry. Large variations in accessory mineral ages within individual clasts (Grange et al., 2013; Thiessen et al., 2018) further complicate interpretations.
View in article


Szumila, I., Trail, D., Erickson, T., Simon, J.I., Wielicki, M.M., Lapen, T., Nakajima, M., Fries, M., Bell, E.A. (2023) Microstructural changes and Pb mobility during the zircon to reidite transformation: Implications for planetary impact chronology. American Mineralogist 108, 1516–1529. https://doi.org/10.2138/am-2022-8604
Show in context

To investigate the effect of shock on Ti coordination in zircon, we collected Ti XANES spectra on ZrSiO4 synthetically shocked to ∼30 GPa (Szumila et al., 2023).
View in article
The products include shocked zircon, amorphous ZrSiO4, and reidite (Szumila et al., 2023) (a high pressure ZrSiO4 polymorph).
View in article


Tailby, N.D., Walker, A.M., Berry, A.J., Hermann, J., Evans, K.A., Mavrogenes, J.A., O’Neill, H.St.C., Rodina, I.S., Soldatov, A.V., Rubatto, D., Sutton, S.R. (2011) Ti site occupancy in zircon. Geochimica et Cosmochimica Acta 75, 905–921. https://doi.org/10.1016/j.gca.2010.11.004
Show in context

A critical component of the thermodynamic theory that dictates the application of the thermometer is that Ti substitutes overwhelmingly for Si on the tetrahedral site in the zircon structure, as opposed to the Zr dodecahedral site, at crustal pressures (Tailby et al., 2011) (see Supplementary Information).
View in article
We report 43 XANES spectra on randomly oriented zircons from these two experimental products to determine the fraction of Ti residing in the Si tetrahedral site (in 4-fold coordination) vs. Zr dodecahedral site (in 8-fold coordination) of zircon (see Methods). Contrary to expectation (Ferry and Watson, 2007; Tailby et al., 2011; Crisp et al., 2023), a significant fraction of the Ti was found to be in the zircon 8-fold site (Tables S-2, S-3).
View in article
Considering these results alongside the other silica saturated experiment at 1300 °C and 1 atm, the differences in Ti site occupancy of zircon represent a drastic change from past studies, which interpreted ∼95 % of the Ti to be present in the Si site (Tailby et al., 2011; Crisp et al., 2023).
View in article
(a) Two standard spectra with Ti present in only 4-fold coordination (cristobalite); and a model spectrum for Ti in 8-fold coordination in zircon along the a axis (Tailby et al., 2011).
View in article
(b) Example pre-edge features for experimental zircon for silica saturated and silica undersaturated (Tailby et al., 2011) experiments.
View in article
The energy shift observed in the pre-edge spectra of the silica undersaturated samples is common, and has been shown to be related to crystal orientation (Tailby et al., 2011; Crisp et al., 2023).
View in article
Differences in pre-edge intensity are due to orientation and/or sample heterogeneity (Tailby et al., 2011).
View in article
(c) Fraction of 4-fold Ti in zircon, each analysis from the three experiments calculated from a linear combination of the two standard spectra presented in (a) (Tailby et al., 2011; Crisp et al., 2023).
View in article


Thiessen, F., Nemchin, A.A., Snape, J.F., Bellucci, J.J., Whitehouse, M.J. (2018) Apollo 12 breccia 12013: Impact-induced partial Pb loss in zircon and its implications for lunar geochronology. Geochimica et Cosmochimica Acta 230, 94–111. https://doi.org/10.1016/j.gca.2018.03.023
Show in context

The Moon exhibits evidence of younger (∼3.4–4 Ga) silica-rich systems, identified through remote sensing (Glotch et al., 2010; Siegler et al., 2023), although detecting such features earlier in lunar history remains challenging. Apollo breccias contain ancient “granite” clasts with 4–4.3 Ga zircon (Meyer et al., 1996), but their complex history of crust formation (Su et al., 2023), meteorite impacts (Thiessen et al., 2018; Valencia et al., 2024) and disequilibrium assemblages (Seddio et al., 2013; Bellucci et al., 2019) obscure primary chemistry. Large variations in accessory mineral ages within individual clasts (Grange et al., 2013; Thiessen et al., 2018) further complicate interpretations.
View in article
Although these pressures exceed typical lunar crustal formation conditions, many surface samples likely experienced transient pressures >2 GPa due to impact shock metamorphism (Thiessen et al., 2018; Valencia et al., 2024).
View in article


Trail, D., Barboni, M., McKeegan, K.D. (2020) Evidence for diverse lunar melt compositions and mixing of the pre-3.9 Ga crust from zircon chemistry. Geochimica et Cosmochimica Acta 284, 173–195. https://doi.org/10.1016/j.gca.2020.06.018
Show in context

The Ti content in zircon is often used as a tool to determine crystallisation temperature (Ferry and Watson, 2007; Trail et al., 2020).
View in article
To further confirm our observations, we performed a separate experiment at different P/T conditions. While the previous experiments were conducted at 1300 °C/1 atm, lunar zircons likely crystallised at a lower T and moderately higher P (Crow et al., 2017; Bellucci et al., 2019; Trail et al., 2020).
View in article
Lunar zircons are exactly that: they are commonly found in breccias and soils as fragments of larger crystals yet they are high value analytical targets because U-Pb ages, typically with 10 Ma precision, can be measured in situ (Trail et al., 2020).
View in article
With age information already available for a suite of samples (Trail et al., 2020), we collected Ti XANES spectra on zircons with ages spanning ∼400 Myr (Table S-4).
View in article
(c) Calculation of Ti in 4-fold coordination, with points averaging lunar zircon measurements in two crystal orientations, plotted against the previously reported 207Pb/206Pb age (Trail et al., 2020).
View in article


Valencia, S.N., Korotev, R.L., Jolliff, B.L. (2024) Compositional analysis of Apollo 12 Granitic Breccia 12013: Insights into protoliths and formation. Geochimica et Cosmochimica Acta 367, 189–212. https://doi.org/10.1016/j.gca.2023.12.034
Show in context

The Moon exhibits evidence of younger (∼3.4–4 Ga) silica-rich systems, identified through remote sensing (Glotch et al., 2010; Siegler et al., 2023), although detecting such features earlier in lunar history remains challenging. Apollo breccias contain ancient “granite” clasts with 4–4.3 Ga zircon (Meyer et al., 1996), but their complex history of crust formation (Su et al., 2023), meteorite impacts (Thiessen et al., 2018; Valencia et al., 2024) and disequilibrium assemblages (Seddio et al., 2013; Bellucci et al., 2019) obscure primary chemistry. Large variations in accessory mineral ages within individual clasts (Grange et al., 2013; Thiessen et al., 2018) further complicate interpretations.
View in article
Although these pressures exceed typical lunar crustal formation conditions, many surface samples likely experienced transient pressures >2 GPa due to impact shock metamorphism (Thiessen et al., 2018; Valencia et al., 2024).
View in article
Instead, we favour a model similar to that ascribed to the origin of the red spot features (Hagerty et al., 2006), in which basaltic underplating and partial melting of pre-existing lunar crust produces highly silicic, high Th lunar melts (Hagerty et al., 2006; Glotch et al., 2010; Gullikson et al., 2016; Valencia et al., 2024) which is consistent with our data (Fig. 4).
View in article
High viscosity minimum melt (?) silicic systems may have developed not far from final emplacement location and/or made significant use of the higher porosity of the fractured crust/mega-regolith, including extensional features associated with impact basins (Valencia et al., 2024).
View in article


Warren, P.H., Taylor, G.J., Keil, K., Shirley, D.N., Wasson, J.T. (1983) Petrology and chemistry of two “large” granite clasts from the moon. Earth and Planetary Science Letters 64, 175–185. https://doi.org/10.1016/0012-821X(83)90202-9
Show in context

This aligns with the occurrence of 1–3 mm felsic clasts found in lunar breccia (Meyer et al., 1996; Seddio et al., 2013) containing impact-generated glass and brecciated material (Warren et al., 1983; Jolliff, 1991).
View in article


Wilson, J.T., Eke, V.R., Massey, R.J., Elphic, R.C., Jolliff, B.L., Lawrence, D.J., Llewellin, E.W., McElwaine, J.N., Teodoro, L.F.A. (2015) Evidence for explosive silicic volcanism on the Moon from the extended distribution of thorium near the Compton-Belkovich Volcanic Complex. Journal of Geophysical Research: Planets 120, 92–108. https://doi.org/10.1002/2014JE004719
Show in context

Remote sensing has identified younger ∼3.4–4.0 Ga plutonic and volcanic silicic rock on the near side of the Moon (Glotch et al., 2010), also characterised by high Th contents (∼10–55 ppm; Hagerty et al., 2006; Wilson et al., 2015) relative to the Procellarum KREEP Terrane (∼5 ppm; Jolliff et al., 2000).
View in article
The Compton–Bel’kovich Volcanic Complex represents a similar highly silicic Th-rich feature on the far side of the Moon, that erupted at ∼3.5 Ga and possibly earlier at ∼3.8 Ga, based on crater size frequency distributions (Jolliff et al., 2011; Wilson et al., 2015; Siegler et al., 2023).
View in article
Yellow squares represent Th contents for younger silicic systems investigated and quantified by remote sensing (Hagerty et al., 2006; Wilson et al., 2015).
View in article


Wu, F.-Y., Li, Q.-L., Chen, Y., Hu, S., Yue, Z.-Y., Zhou, Q., Wang, H., Yang, W., Tian, H.-C., Zhang, C., Li, J.-H., Li, L.-X., Hui, H.-J., Li, C.-L., Lin, Y.-T., Li, X.-H., Delano, J.W. (2024) Lunar Evolution in Light of the Chang’e-5 Returned Samples. Annual Review of Earth and Planetary Sciences 52, 159–194. https://doi.org/10.1146/annurev-earth-040722-100453
Show in context

But overall, the zircon data represent highly silicic systems that are under-represented in sample collections from Apollo, Luna, meteorites, and Chang’e 5 missions (Neal, 2009; Wu et al., 2024).
View in article


Yang, J., Ju, D., Pang, R., Li, R., Liu, J., Du, W. (2023) Significance of silicate liquid immiscibility for the origin of young highly evolved lithic clasts in Chang’E-5 regolith. Geochimica et Cosmochimica Acta 340, 189–205. https://doi.org/10.1016/j.gca.2022.11.008
Show in context

Similar processes are observed in younger silicic localities such as the Gruithuisen region (Qiu et al., 2023), Aristarchus, and Lassell regions, where silicic lithologies are exposed through later impact cratering (Glotch et al., 2010; Clegg-Watkins et al., 2017; Yang et al., 2023).
View in article



top

Supplementary Information

Abstract | Introduction | Results | Acknowledgements | References | Supplementary Information


The Supplementary Information includes:
  • Methods
  • Grain Orientation and Intensity of the Pre-edge Feature
  • Impact Modification of ZrSiO4
  • Determining the Th Melt Contents of Zircon Bearing Igneous Systems
  • Implications for Ti-in-Zircon Thermometry
  • Explanation for Ti Zircon XANES Differences from Previous Studies
  • Silica Activity Calculation
  • Other High Silica Activity Lunar Environments
  • Correlation with Other Trace Elements
  • Figures S-1 to S-9
  • Tables S-1 to S-5
  • Supplementary Information References


Download the Supplementary Information (PDF)

Download Tables S-1 to S-5 (.xlsx)
top

Figures



Figure 1 XANES spectral analysis of Ti coordination in zircon. (a) Two standard spectra with Ti present in only 4-fold coordination (cristobalite); and a model spectrum for Ti in 8-fold coordination in zircon along the a axis (Tailby et al., 2011

Tailby, N.D., Walker, A.M., Berry, A.J., Hermann, J., Evans, K.A., Mavrogenes, J.A., O’Neill, H.St.C., Rodina, I.S., Soldatov, A.V., Rubatto, D., Sutton, S.R. (2011) Ti site occupancy in zircon. Geochimica et Cosmochimica Acta 75, 905–921. https://doi.org/10.1016/j.gca.2010.11.004

). (inset) Linear mixing of the pre-edge feature showing model spectra for different fractions of Ti in 4-fold coordination. (b) Example pre-edge features for experimental zircon for silica saturated and silica undersaturated (Tailby et al., 2011

Tailby, N.D., Walker, A.M., Berry, A.J., Hermann, J., Evans, K.A., Mavrogenes, J.A., O’Neill, H.St.C., Rodina, I.S., Soldatov, A.V., Rubatto, D., Sutton, S.R. (2011) Ti site occupancy in zircon. Geochimica et Cosmochimica Acta 75, 905–921. https://doi.org/10.1016/j.gca.2010.11.004

) experiments. The 1300 °C, 1 atm, SiO2 undersaturated experiment (srilankite (Ti,Zr)O4, rutile, zircon) has a pre-edge intensity similar to that of the 4-fold standard, whereas the intensity of the pre-edge feature for the two silica saturated experiments is muted meaning that a significant fraction of 8-fold Ti is present in zircon. The energy shift observed in the pre-edge spectra of the silica undersaturated samples is common, and has been shown to be related to crystal orientation (Tailby et al., 2011

Tailby, N.D., Walker, A.M., Berry, A.J., Hermann, J., Evans, K.A., Mavrogenes, J.A., O’Neill, H.St.C., Rodina, I.S., Soldatov, A.V., Rubatto, D., Sutton, S.R. (2011) Ti site occupancy in zircon. Geochimica et Cosmochimica Acta 75, 905–921. https://doi.org/10.1016/j.gca.2010.11.004

; Crisp et al., 2023

Crisp, L.J., Berry, A.J., Burnham, A.D., Miller, L.A., Newville, M. (2023) The Ti-in-zircon thermometer revised: The effect of pressure on the Ti site in zircon. Geochimica et Cosmochimica Acta 360, 241–258. https://doi.org/10.1016/j.gca.2023.04.031

). Differences in pre-edge intensity are due to orientation and/or sample heterogeneity (Tailby et al., 2011

Tailby, N.D., Walker, A.M., Berry, A.J., Hermann, J., Evans, K.A., Mavrogenes, J.A., O’Neill, H.St.C., Rodina, I.S., Soldatov, A.V., Rubatto, D., Sutton, S.R. (2011) Ti site occupancy in zircon. Geochimica et Cosmochimica Acta 75, 905–921. https://doi.org/10.1016/j.gca.2010.11.004

). Spectra were carefully evaluated for inclusions (see Supplementary Information and Fig. S-1). (c) Fraction of 4-fold Ti in zircon, each analysis from the three experiments calculated from a linear combination of the two standard spectra presented in (a) (Tailby et al., 2011

Tailby, N.D., Walker, A.M., Berry, A.J., Hermann, J., Evans, K.A., Mavrogenes, J.A., O’Neill, H.St.C., Rodina, I.S., Soldatov, A.V., Rubatto, D., Sutton, S.R. (2011) Ti site occupancy in zircon. Geochimica et Cosmochimica Acta 75, 905–921. https://doi.org/10.1016/j.gca.2010.11.004

; Crisp et al., 2023

Crisp, L.J., Berry, A.J., Burnham, A.D., Miller, L.A., Newville, M. (2023) The Ti-in-zircon thermometer revised: The effect of pressure on the Ti site in zircon. Geochimica et Cosmochimica Acta 360, 241–258. https://doi.org/10.1016/j.gca.2023.04.031

). When aqtzSiO2, only ∼50 % of the Ti is occupying the 4-fold Si site, whereas with aqtzSiO2 nearly 100 % of the Ti is in the 4-fold Si site. Error bars represent the 1 s.d. for each group. Experimental products can be found in Table S-1.
Back to article


Figure 2 Ti XANES spectra and coordination in lunar zircon. (a, b) Two representative grain analyses showing XANES spectra before and after a 90° rotation. The top panel shows a significant change in pre-edge absorption intensity, while the lower panel, more typical, shows virtually no change. (c) Calculation of Ti in 4-fold coordination, with points averaging lunar zircon measurements in two crystal orientations, plotted against the previously reported 207Pb/206Pb age (Trail et al., 2020

Trail, D., Barboni, M., McKeegan, K.D. (2020) Evidence for diverse lunar melt compositions and mixing of the pre-3.9 Ga crust from zircon chemistry. Geochimica et Cosmochimica Acta 284, 173–195. https://doi.org/10.1016/j.gca.2020.06.018

). Dashed lines show the average fraction of 4-fold Ti at 1300 °C, and shaded yellow areas represent the 1 s.d. of the fraction in synthetic zircon.
Back to article


Figure 3 Th contents of highly silicic lunar systems. Circles represent calculated Th contents of igneous systems derived from zircon-melt Th experimental partitioning (see Supplementary Information; Fig. S-5 and Table S-5). Yellow squares represent Th contents for younger silicic systems investigated and quantified by remote sensing (Hagerty et al., 2006

Hagerty, J.J., Lawrence, D.J., Hawke, B.R., Vaniman, D.T., Elphic, R.C., Feldman, W.C. (2006) Refined thorium abundances for lunar red spots: Implications for evolved, nonmare volcanism on the Moon. Journal of Geophysical Research: Planets 111, E06002. https://doi.org/10.1029/2005JE002592

; Wilson et al., 2015

Wilson, J.T., Eke, V.R., Massey, R.J., Elphic, R.C., Jolliff, B.L., Lawrence, D.J., Llewellin, E.W., McElwaine, J.N., Teodoro, L.F.A. (2015) Evidence for explosive silicic volcanism on the Moon from the extended distribution of thorium near the Compton-Belkovich Volcanic Complex. Journal of Geophysical Research: Planets 120, 92–108. https://doi.org/10.1002/2014JE004719

). The ages for the silicic systems identified by remote sensing data are estimated through crater counting. The average Th contents from the PKT mare terrain (4.9 ppm) and other mare (2.2 ppm) are also shown for comparison (Jolliff et al., 2000

Jolliff, B.L., Gillis, J.J., Haskin, L.A., Korotev, R.L., Wieczorek, M.A. (2000) Major lunar crustal terranes: Surface expressions and crust-mantle origins. Journal of Geophysical Research: Planets 105, 4197–4216. https://doi.org/10.1029/1999JE001103

). MMD1/2, first and second eruptions of Middle Mairan Dome; MTD, Mairan T Dome; GDδ/γ, Gruithuisen Dome (δ) and (γ); LM, Lassell Massif; HD, Hansteen Alpha; CBVC, Compton–Bel’kovich Volcanic Complex.
Back to article


Figure 4 Simplified illustration (not to scale) of the lunar silicic systems whose formation was initiated by partial melting of fertile source rocks such as a KREEP basalt by basaltic underplating (Hagerty et al., 2006

Hagerty, J.J., Lawrence, D.J., Hawke, B.R., Vaniman, D.T., Elphic, R.C., Feldman, W.C. (2006) Refined thorium abundances for lunar red spots: Implications for evolved, nonmare volcanism on the Moon. Journal of Geophysical Research: Planets 111, E06002. https://doi.org/10.1029/2005JE002592

). High viscosity minimum melt (?) silicic systems may have developed not far from final emplacement location and/or made significant use of the higher porosity of the fractured crust/mega-regolith, including extensional features associated with impact basins (Valencia et al., 2024

Valencia, S.N., Korotev, R.L., Jolliff, B.L. (2024) Compositional analysis of Apollo 12 Granitic Breccia 12013: Insights into protoliths and formation. Geochimica et Cosmochimica Acta 367, 189–212. https://doi.org/10.1016/j.gca.2023.12.034

). Early silica-rich domes were likely destroyed by large basin-forming impacts, or covered by impact ejecta (Qiu et al., 2023

Qiu, D., Sasaki, S., Yan, J., Ye, M., Deng, Q., Liang, F., Liu, L., Li, F. (2023) Buried Silicic Volcanoes Discovered in the Gruithuisen Region on the Moon. Geophysical Research Letters 50, e2023GL103336. https://doi.org/10.1029/2023GL103336

).
Back to article

  • Contact us
  • |
  • Subscribe
  • |
  • Sign up to the EAG newsletter
  • Connect with us
  • Bluesky
  • facebook
  • Linkedin
  • 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