Abstract
Lunar chronology models are built by associating the radiometric ages of samples returned by the Apollo and Luna missions measured in the laboratory with compiled crater distributions of those sites. Such models have not only been widely used to determine the absolute ages of various regions on the Moon1,2,3,4,5,6, but have also been generalized to date the surfaces of the rocky bodies of the inner Solar System7,8,9,10,11,12. However, there is a gap in lunar samples ages between 3.0 Gyr ago and 1.0 Gyr ago13, which occupies almost half of the history of the Moon. Chang’e-5, the first lunar sample return mission since the Luna 24 lander in 1976, brought back basalt material from a young mare area that has been dated to the centre of this gap at 2.030 ± 0.004 Gyr old14. Using this radiometric age, we updated the most widely used chronology models, focusing in particular on the Neukum model13. We found that the updated model is consistent with a combination of an exponential decrease and a linear rate. The updated chronology gives older ages with respect to the Neukum model for most of the lunar history, with a maximum difference of 0.24 Gyr at 2.55 Gyr ago. Differences from other models are of comparable magnitude or greater. These results have important implications for the chronology and impact history of the inner Solar System.
Similar content being viewed by others
Main
On 1 December 2020, Chang’e-5 successfully landed in Northern Oceanus Procellarum at (51.92° W, 43.06° N)15. The landing site is located in the geologic unit Em4 mapped by Qian et al.16 (see their fig. 1b), which belongs to the Eratosthenian-aged mare. The landing area is characterized as a young phase of volcanism16 with the basalt thickness varying from ∼39 to ∼63 m overlying the old Imbrian-aged mare units17. The thickness of the lunar regolith in Em4 is ~3–15 m (ref. 18). Before landing and sampling, our team identified and mapped craters >200 m in diameter in Em4 (Fig. 1b) and other units in the landing area19.
Automated sampling was the core task of the Chang’e-5 mission. Two sampling methods were employed: surficial sample collection by a robotic arm and underground sample collection by a drill. The Chang’e-5 samples that have been radiometrically measured were scooped from the lunar regolith surface20. The samples were sealed in a storage device to ensure that they were unaffected by the external environment when returning to Earth. The samples were carried by the ascender and transferred to the returner in lunar orbit. The returner was subsequently separated from the orbiter and finally carried the samples back to Earth on 17 December 2020. The total weight of the returned samples was 1,731 g. The samples collected by the robotic arm have been radiometrically measured at the Institute of Geology and Geophysics, Chinese Academy of Sciences14,20.
Mineral and radiometric analyses were conducted at the Institute of Geology and Geophysics, Chinese Academy of Sciences, using 3 g of the surficial sample. Basalt clasts used for radiometric measurements have similar major mineral constituents and formed during the same final crystallization stage of the magma14. Radiometric measurements using the Pb–Pb dating technique were taken on both Zr-rich minerals and phosphates, and reveal the same age of 2.030 ± 0.004 Gyr ago (Ga), which represents the absolute age of the basalt within the landing area of Chang’e-520. This age has been used here to update the lunar cratering chronology function that is widely used at present.
We updated the chronology function established by Neukum13 using the the radiometric age and corresponding terrain N(1) value; that is, the total number of craters with diameters ≥1 km per square kilometre. The crater size–frequency distribution in Em4 has been studied previously and different N(1) values have been derived17,19,21,22. Here we use the result from Jia et al.19 (that is, N(1) = (1.74 ± 0.32) × 10−3 km−2), which was derived from crater measurements in a high-resolution orthophoto map of the Chang’e-5 landing area generated from more than 700 Narrow Angle Camera images from the Lunar Reconnaissance Orbiter Camera. In studying the size–frequency distribution with high-resolution images, Hiesinger et al.23 updated several N(1) values of the young units sampled. Based on the data used in Neukum13, the new data from Chang’e-5 and the updated data for Copernicus crater in Hiesinger et al.23 were used to fit the chronology function model. As shown in Fig. 2, the data point from Chang’e-5 does not change the chronology function model greatly from that of Neukum13; thus a function with an exponential decline and a linear rate is again used to fit the data presently available. Our new lunar cratering chronology function is as follows, where t is the model age:
Figure 2 also shows the model age difference (that is, new model − Neukum model) with respect to N(1). The new chronology model gives older ages for all N(1) values except 0.055 km−2 to 0.0065 km−2 (with corresponding model ages from 3.98 Ga to 3.59 Ga), for which the new chronology model gives slightly younger ages. The maximum difference appears at 2.55 Ga, with a corresponding N(1) of 0.0021 km−2, and the maximum model age difference is 0.24 Gyr. It is worth noting that the norm of the residuals decreases quickly with the iteration times in the least-squares fitting (see details in the Methods), demonstrating that the residuals have been stabilized with the iteration. A detailed comparison with other N(1) values and other chronologies is presented in the Methods.
The impact rate resulting in craters of 1 km in diameter and larger is given by the derivative of the impact chronology function with respect to t; that is, \(\phi = \frac{{\partial N}}{{\partial t}}\). From the new lunar cratering chronology model (equation (1)), the derived impact rate for diameter D ≥ 1 km is as follows:
Similar to the rate from Neukum13, the impact rate resulting in craters of 1 km in diameter and larger from the new model also exponentially decreases before ~3.0 Ga, and then an almost constant value is obtained. The difference in the impact rate resulting in 1 km and larger craters between the two models is also shown in Fig. 3. Before 3.866 Ga, the cumulative cratering rate from the new model is always smaller than that from Neukum13, with the difference becoming larger for older ages. From 3.866 Ga to ~3.156 Ga, the new model gives a slightly higher cumulative impact rate. From that time until today, the two models are almost same.
The differences in ages and cratering rates between the new and previous models are shown in Fig. 4. The models from Neukum13, Wagner et al.24, Le Feuvre and Wieczorek25 and our new model are generally close, with maximum differences of ~0.2 Gyr. These four models have significant differences from those of Marchi et al.26 and Robbins27, with maximum differences of up to 1.2 Gyr. So far, the model from Neukum13 has been the most widely used in lunar and planetary research. Our new model, which was updated from Neukum13 with the new radiometric age of the Chang’e-5 sample, should be more reliable and applicable in future research.
In summary, a new chronology model based on the radiometric age of Chang’e-5 samples is presented and compared with previous models. The radiometric age from Chang’e-5 samples falls within the gap between 3.0 and 1.0 Ga and serves as a ‘golden spike’ for lunar chronology. The result confirms that a summation of an exponential equation and a linear equation is applicable to describe the chronology model. The new model refines the previous model and will enhance the reliability of ages estimated with the crater size–frequency distribution method. We propose that the new chronology model is used in lunar surface dating, and that the chronology models for inner Solar System bodies (for example, Mars, Mercury and asteroids) are revised on the basis of this new lunar model while accounting for differences in their impact rate, impact velocity, gravity and target properties compared with the Moon28,29,30.
Methods
Crater counting has long been used to determine the model age of units of the lunar surface31. The rationale of the approach is to fit the observed crater size–frequency distribution of a given surface unit to a known crater production function1,32, which is further used to derive the absolute model ages along with a chronology function chronology function calibrated to radiometric dating from lunar samples13,33.
The production function refers to the standard crater size–frequency distribution of the lunar surface, which was established by combining a number of crater size–frequency distributions of measured craters on geologic units of different ages and in overlapping crater diameter ranges10,13,33,34. The production function widely used at present for the Moon was derived by Neukum13, Hartmann34 and Neukum et al.10 based on the hypothesis that the craters are randomly distributed on the lunar surface. Possible asymmetries in the distribution of impact craters on the Moon were proposed and their effects on the production function were modelled in several studies24,25,35,36,37. However, these models are mutually inconsistent, given the asymmetries in crater distribution, and they are not necessarily widely accepted at present23.
The chronology function is usually expressed as the relationship between the radiometric ages of the returned samples and the cumulative frequency of craters larger than some particular diameter, such as 1 km, that is, N(1) in the corresponding sampled area. However, the values of N(1) were often calculated through extrapolation using the production function instead of direct measurements, thus the production function also influences the chronology function. The current widely used chronology function is based on the values of N(1) from Neukum13, and corresponding radiometric ages. Hiesinger et al.23 updated N(1) for young ages with the highest-resolution lunar orbital images so far and indicated that the updated values (especially for Copernicus crater) are more consistent with the chronology function from Neukum13. Therefore, N(1) values from Neukum13 and the updated N(1) for Copernicus crater are used in this study. However, the updated N(1) from Hiesinger et al.23 was calculated using the production function from Neukum et al.10. It is better to derive the values of N(1) using the same production function. The production functions of Neukum13 and Neukum et al.10 have a similar relationship with different coefficients. If the diameter of the craters used to derive the N(1) is D1, the ratio of N(1) values from the two production functions can be expressed as follows:
where N(1)* is from Neukum13 and N(1) is from Neukum et al.10; the coefficients of ai and \(a_i^ \ast\) are for Neukum et al.10 and Neukum13, respectively. Supplementary Fig. 1 shows the above equation for crater diameters ranging from 10 m to 300 km, which is also the recommended diameter interval for the two chronology functions. It is clear that for crater diameter less than 1.0 km the values of N(1) from the two models are very similar. For Copernicus crater in Hiesinger et al.23 (See Fig. 10 therein), the above ratio is between ~1.01 and 1.02. Thus N(1) is also very similar for the two production functions and was directly used here, as the accurate ratio is unknown and the difference plays a small role in generating the new chronology function.
Through the distribution of all of the radiometric ages, including that from the Chang’e-5 samples, the chronology model can be expressed as the summation of an exponential equation and a linear equation. Thus a function similar to Neukum13 was adopted. The trust-region-reflective algorithm38,39 of the nonlinear least-squares fitting was used to solve the coefficients of the chronology function by minimizing the sum of squared residuals through iteration. The uncertainties of the ages and N(1) values were not incorporated in the model fitting for simplicity and easy comparison with previous models. The coefficients of the Neukum model13 were used as the initial values, and the iteration was executed 1,250 times, when the norm of the residuals stabilized at the level of 1.42242 × 10−5.
In addition to the N(1) value provided by Jia et al.27 of (1.74 ± 0.32)×10−3 km−2, Wu et al.22 and Qian et al.17 gave N(1) values for Em4 of 1.24 × 10−3 km−2 and 1.28 × 10−3 km−2, respectively. The production function of Neukum et al.10 was used for both studies. If the above values are directly used to fit the new chronology function, as we have done here, the results are as follows:
Equations (4) and (5) were derived with the N(1) of Em4 unit from Wu et al.22 and Qian et al.17, respectively. Comparing the two equations with equation (1), the differences are not large. Supplementary Fig. 2 shows the differences between equations (1), (4) and (5) and the chronology function from Neukum13. These curves have similar shapes; that is, the maximum difference appears at a model age of about 2.52 Ga, and the modifications are small at young and old ages for all three models. The differences between the three models result from the different N(1) values in the Chang’e-5 landing area. This indicates the importance of utilizing high-resolution images and the subjectivity of identifying craters, including the careful identification and elimination of secondary craters when applying the crater size–frequency distribution method to estimate the age of a lunar surface unit.
Besides the model from Neukum13, there are other chronology function models24,25,28,29 and their differences from the new model in this research are shown in Fig. 4. In summary, the largest differences in the ages and cumulative cratering rates are between the new model and that from Robbins27; the new model is most similar to that from Neukum13.
Data availability
The data points used to fit the new chronology model include those from Neukum13, Hiesinger et al.23, Jia et al.19, Wu et al.22, Qian et al.17 and the radiometric age (2.030 ± 0.004 Ga) from Li et al.14. The mapped craters in the Chang'e-5 landing area can be found at https://zenodo.org/record/5615501#.YXvMT5pByF4.
Code availability
The code for fitting the new lunar cratering chronology function can be found at https://zenodo.org/record/5615459#.YXvHJppByF4.
Change history
11 March 2022
A Correction to this paper has been published: https://doi.org/10.1038/s41550-022-01649-4
References
Hiesinger, H. et al. Ages of mare basalts on the lunar nearside. J. Geophys. Res. 105, 29239–29275 (2000).
Hiesinger, H. et al. Ages and stratigraphy of mare basalts in Oceanus Procellarum, Mare Nubium, Mare Cognitum, and Mare Insularum. J. Geophys. Res. 108, 5065 (2003).
Hiesinger, H. et al. Ages and stratigraphy of lunar mare basalts in Mare Frigoris and other nearside maria based on crater size-frequency distribution measurements. J. Geophys. Res. 115, E03003 (2010).
Whitten, J. et al. Lunar mare deposits associated with the Orientale impact basin: new insights into mineralogy, history, mode of emplacement, and relation to Orientale Basin evolution from Moon Mineralogy Mapper (M3) data from Chandrayaan-1. J. Geophys. Res. 116, E00G9 (2011).
Kirchoff, M. R. et al. Ages of large lunar impact craters and implications for bombardment during the Moon’s middle age. Icarus 225, 325–341 (2013).
Yue, Z. et al. Refined model age for Orientale Basin derived from zonal crater dating of its ejecta. Icarus 346, 113804 (2020).
Ivanov, B. A. Mars/Moon cratering rate ratio estimates. Space Sci. Rev. 96, 87–104 (2001).
Michael, G. G. Planetary surface dating from crater size–frequency distribution measurements: multiple resurfacing episodes and differential isochron fitting. Icarus 226, 885–890 (2013).
Hartmann, W. K. & Neukum, G. Cratering chronology and the evolution of Mars. Space Sci. Rev. 96, 165–194 (2001).
Neukum, G., Ivanov, B. A. & Hartmann, W. K. Cratering records in the Inner Solar System in relation to thelunar reference. Space Sci. Rev. 96, 55–86 (2001).
Hartmann, W. K. Martian cratering 8: isochron refinement and the chronology of Mars. Icarus 174, 294–320 (2005).
Hartmann, W. K. & Daubar, I. J. Martian cratering 11. Utilizing decameter scale crater populations to study Martian history. Meteorit. Planet. Sci. 52, 493–510 (2017).
Neukum, G. Meteorite Bombardment and Dating of Planetary Surfaces. Habilitation thesis, Univ. Munich (1983).
Li, Q. L. et al. Two billion-year-old volcanism on the Moon from Chang’e-5 basalts. Nature 600, 54–58 (2021).
Wang, J. et al. Localization of the Chang’e-5 lander using radio-tracking and image-based methods. Remote Sens. 13, 590–601 (2021).
Qian, Y. et al. The regolith properties of the Chang’e-5 landing region and the ground drilling experiments using lunar regolith simulants. Icarus 337, 113508 (2020).
Qian, Y. et al. Young lunar mare basalts in the Chang’e-5 sample return region, northern Oceanus Procellarum. Earth Planet. Sci. Lett. 555, 116702 (2021).
Yue, Z. et al. Lunar regolith thickness deduced from concentric craters in the CE-5 landing area. Icarus 329, 46–54 (2019).
Jia, M. et al. A catalogue of impact craters larger than 200 m and surface age analysis in the Chang’e-5 landing area. Earth Planet. Sci. Lett. 541, 116272 (2020).
Hu, S. et al. A dry lunar mantle reservoir for young mare basalts of Chang’e-5. Nature 600, 49–53 (2021).
Qian, Y. Q. et al. Geology and scientific significance of the Rümker region in northern Oceanus Procellarum: China’s Chang’E-5 landing region. J. Geophys. Res. Planets 123, 1407–1430 (2018).
Wu, B. et al. Rock abundance and crater density in the candidate Chang’E-5 landing region on the Moon. J. Geophys. Res. Planets 123, 3256–3272 (2018).
Hiesinger, H. et al. How old are young lunar craters? J. Geophys. Res. 117, E00H10 (2012).
Wagner, R., Head, J. W. III, Wolf, U. & Neukum, G. Stratigraphic sequence and ages of volcanic units in the Gruithuisen region of the Moon. J. Geophys. Res. 107, E00H10 (2002).
Le Feuvre, M. & Wieczorek, M. A. Nonuniform cratering of the Moon and a revised crater chronology of the inner Solar System. Icarus 214, 1–20 (2011).
Marchi, S. et al. A New chronology for the Moon and Mercury. Astrophys. J. 137, 4936–4948 (2009).
Robbins, S. J. New crater calibrations for the lunar crater-age chronology. Earth Planet. Sci. Lett. 403, 188–198 (2014).
Marchi, S. et al. Small crater populations on Vesta. Planet. Space Sci. 103, 96–103 (2014).
Schmedemann, N. et al. The cratering record, chronology and surface ages of (4) Vesta in comparison to smaller asteroids and the ages of HED meteorites. Planet. Space Sci. 103, 104–130 (2014).
Chapman, C. R. et al. Cratering on Ida. Icarus 120, 77–86 (1996).
Shoemaker, E. M., Hackman, R. J. & Eggleton, R. E. Interplanetary correlation of geologic time. Adv. Astronaut. Sci. 8, 70–79 (1962).
Michael, G. G. & Neukum, G. Planetary surface dating from crater size–frequency distribution measurements: partial resurfacing events and statistical age uncertainty. Earth Planet. Sci. Lett. 294, 223–229 (2010).
Stöffler, D. & Ryder, G. Stratigraphy and isotope ages of lunar geologic units: chronological standard for the inner solar system. Space Sci. Rev. 96, 9–54 (2001).
Hartmann, W. K. Martian cratering VI: crater count isochrons and evidence for recent volcanism from Mars Global Surveyor. Meteorit. Planet. Sci. 34, 167–177 (1999).
Morota, T., Ukai, T. & Furumoto, M. Influence of the asymmetrical cratering rate on the lunar cratering chronology. Icarus 173, 322–332 (2005).
Le Feuvre, M. & Wieczorek, M. A. Nonuniform cratering of the terrestrial planets. Icarus 197, 291–306 (2008).
Gallant, J., Gladman, B. & Cuk, M. Current bombardment of the Earth-Moon system: emphasis on cratering asymmetries. Icarus 202, 371–382 (2009).
Coleman, T. F. & Li, Y. An interior trust region approach for nonlinear minimization subject to bounds. SIAM J. Optim. 6, 418–445 (1996).
Coleman, T. F. & Li, Y. On the convergence of reflective Newton methods for large-scale nonlinear minimization subject to bounds. Math. Program. 67, 189–224 (1994).
Acknowledgements
The Change'e-5 mission was carried out by the Chinese Lunar Exploration Program. This work was supported by the Strategic Priority Program of the Chinese Academy of Sciences (grant no. XDB41000000), the National Natural Science Foundation of China (grant nos 41941003, 41773065 and 41972321) and the Key Research Program of Frontier Sciences, CAS (grant no. QYZDY-SSW-DQC028).
Author information
Authors and Affiliations
Contributions
Z.Y. and K.D. planned the research and prepared the manuscript with the help of Z.O. K.D., W.W., Z.L., S.G., B.L., M.P., Y.W., M.J. and J. L. processed the data.
Corresponding author
Ethics declarations
Competing interests
The authors declare no competing interests.
Peer review
Peer review information
Nature Astronomy thanks Nicolle Zellner, Michelle Kirkhoff and the other, anonymous, reviewer(s) for their contribution to the peer review of this work.
Additional information
Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary information
Supplementary Information
Supplementary Figs. 1 and 2.
Rights and permissions
About this article
Cite this article
Yue, Z., Di, K., Wan, W. et al. Updated lunar cratering chronology model with the radiometric age of Chang’e-5 samples. Nat Astron 6, 541–545 (2022). https://doi.org/10.1038/s41550-022-01604-3
Received:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1038/s41550-022-01604-3
This article is cited by
-
Mineralogy and chronology of the young mare volcanism in the Procellarum-KREEP-Terrane
Nature Astronomy (2023)
-
Progresses and prospects of impact crater studies
Science China Earth Sciences (2023)
-
Innovative developments in lunar and planetary science promoted by China’s lunar exploration
Science China Earth Sciences (2023)
-
Spectral interpretation of late-stage mare basalt mineralogy unveiled by Chang’E-5 samples
Nature Communications (2022)