Main

It has long been argued that there could be water and other volatile species at the surface of the Moon1,2,3. Renewed lunar exploration and advances in remote-sensing measurements in the late 1990s allowed the neutron spectrometer on board the Lunar Prospector mission to confirm the existence of water ice at the lunar poles4. Following this, the Moon mineralogy mapper instrument on board the Chandrayaan-1 spacecraft detected the diagnostic absorption bands of hydroxyl and/or water at 2.8–3.0 μm on the surface of the Moon5. Detailed analyses indicated that hydroxyl/water could be present all over the Moon’s surface, which also displayed temporal and spatial variations, with equivalent water abundances ranging between ~10 and 1,000 μg g−1 (refs. 6,7,8,9,10). Furthermore, the Lunar Crater Observation and Sensing Satellite impact experiment carried out in 2009 provided direct evidence for high water-ice abundances (5.6 ± 2.9 wt%) in permanently shadowed regions within Cabeus crater11. Elevated water-ice abundance in lunar polar regions was further supported by the neutron flux measurements performed by the Lunar Exploration Neutron Detector on board the Lunar Reconnaissance Orbiter spacecraft12. Recently, the neutral mass spectrometer on the Lunar Atmosphere and Dust Environment Explorer detected exospheric water liberated by meteoroid impacts13, and ground-based telescope observations detected molecular water on the lunar surface14. Today, there is little doubt that most of the Moon’s surface harbours water in one form or another. However, the origin(s) of this lunar surface water and its spatial distribution and evolution during regolith gardening remain largely unknown, despite key implications for future lunar surface exploration and for better understanding the (sub)surface water reservoir and processing on Solar System airless bodies.

There are several potential sources and processes that could have contributed to the water inventory at the surface of the Moon15, such as (1) solar wind implantation16,17,18,19,20, (2) outgassing of volatiles during lunar volcanism21,22,23, (3) deposition of volatile-bearing pyroclastic deposits24 and minerals25 and (4) delivery by impacts of comets and asteroids26,27,28,29. It is generally thought that solar wind hydrogen-ion implantation could react with surface minerals to produce hydroxyl or water in lunar soils30,31. The surface water produced and/or delivered on the Moon’s equatorial regions may migrate to polar regions, driven by temperature oscillations16,19,32,33, and could also be partially released to space13. Hence, a scenario of lunar water cycle was further proposed to describe the retention, release and replenishment of water on the surface of the Moon13. To sustain a water cycle at the surface of the Moon, there should be a hydrated layer (reservoir) at depth in lunar soils13. However, finding this water reservoir has remained elusive, despite several studies having investigated the water inventory of fine mineral grains in lunar soils (for example, refs. 34,35), impact-produced agglutinates27, volcanic rocks (for example, ref. 26) and pyroclastic glass beads21,36. Consequently, there must be a yet-unidentified water reservoir in lunar soils that has the capacity to buffer a lunar surface water cycle.

The lunar soils returned by the Apollo, Luna and Chang’e-5 (CE5) missions comprise mainly rock and mineral fragments, pyroclastic glass beads, impact-produced agglutinates and impact glass beads37,38. Among these components, the water inventory of impact glass beads has not yet been investigated in detail, despite these glassy beads being potential candidates for playing a significant role in a lunar surface water cycle. To investigate this possibility, we carried out a systematic characterization of the petrography, major element composition, Raman characteristics, water abundance and hydrogen isotope composition on the impact glass beads returned by the CE5 mission, aiming to identify and characterize the missing water reservoir on the Moon’s surface. A recent geochronological study of CE5 impact glass beads has shown that they formed more or less continuously for the past 2 Gyr, with prominent peaks in formation ages at ~575 million years ago (Ma), 380 Ma, 68 Ma and 35 Ma, originating predominantly from a few 1- to 5-km-diameter impact craters in the Em4 (Eratosthenian-aged mare) basaltic unit where CE5 landed39.

Water abundance and hydrogen isotope composition of impact glasses

The studied impact glass beads were hand picked from a lunar soil sample (sample CE5C0100YJFM00103; ~1 g) scooped by the CE5 lander robotic arm. A total of 117 individual spherical glass beads were characterized using field-emission scanning electron microscopy (SEM), electron probe microanalyser (EPMA) and Raman spectroscopy (Fig. 1, Supplementary Figs. 1 and 2, Supplementary Table 1 and Extended Data Fig. 1). Among them, 32 impact glass beads characterized by smooth exposed surfaces and chemical compositions consistent with those of bulk CE5 basalts were chosen for in situ water abundance and H isotope analysis using a nanoscale secondary ion mass spectrometer (NanoSIMS) (Fig. 1, Extended Data Figs. 2 and 3 and Supplementary Table 1). The CE5 impact glass beads contain 0–1,909 μg g−1 equivalent H2O, with H isotope compositions (given in delta notation, δD = 1000 × (D/Hsample/D/HSMOW−1), where D/Hsample is the measured D/H ratio measured in the sample and D/HSMOW is the D/H ratio (1.5576 × 10−4) of the standard mean ocean water) ranging from −990 ± 6‰ to 522 ± 440‰ and which are negatively correlated with the water abundances (Fig. 2). We also analysed water abundances and H isotope compositions along six transects in five glass beads (CE5#33,003, CE5#33,036, CE5#33,046, CE5#33,052 and CE5#33,076; Fig. 1), which show elevated water abundances at the rims that gradually decrease towards the cores while δD values decrease from cores to rims (Fig. 3 and Extended Data Figs. 4 and 5). Three glass beads display sharp decreases in water content from the rims towards the cores, while the other two glass beads display a smoother water abundance variation (Fig. 3 and Extended Data Figs. 4 and 5).

Fig. 1: Representative back-scattered electron images of CE5 impact glass beads after nanoSIMS analyses.
figure 1

a, CE5#33,002. b, CE5#33,003. c, CE5#33,036. d, CE5#33,046. e, CE5#33,052. f, CE5#33,076. These impact glass beads have spherical shapes and have homogeneous chemical compositions (Supplementary Table 2). Beads CE5#33,003 and CE5#33,036 are coated by finer agglutinates. Six profiles carried out on five impact glass beads in the second session are outlined in red. Blue squares indicate locations of NanoSIMS in the first analytical session. Each NanoSIMS analysis pit is 7 × 7 μm2.

Fig. 2: Water abundances and δD values measured in CE5 impact glass beads.
figure 2

All glass bead analyses define a negative correlation between the water abundances and δD values, which is consistent with a two-endmember mixing model (red lines) that indicates addition of solar wind-derived water with δD value of −990‰ into a low water abundance (H2O0 = 5–50 μg g−1) CE5 impact glass with δD value of +500‰ (ref. 50). De-gassing of H2 from CE5 impact glasses cannot be responsible for the large observed H2O–δD variations (brown line). The inset diagram is the zoom-in view of water abundance range of 0–100 μg g−1. For more details, see Methods. The literature data for melt inclusions, pyroclastic glasses and impact-produced agglutinates are also plotted for comparison, and the datasets are listed in Supplementary Table 4. All error bars correspond to 2 s.d. analytical uncertainties. SMOW, standard mean ocean water. H2O0 is the water abundance of impact glass beads prior to the inward diffusion of solar wind-derived water.

Fig. 3: Two profiles across CE5 impact glass bead CE5#33,036.
figure 3

a,b, Water abundance (a) and δD values (b) measured along Profile 1, as outlined in Fig. 1c. c,d, Water abundance (c) and δD values (d) measured along Profile 2, as outlined in Fig. 1c. The apparent distances along the profiles were measured from the rim on a side (Fig. 1c). A water-diffusion coefficient of 20.84 μm2 yr−1 at T = 360 K was used48. The maximum diffusion durations for Profile 1 and Profile 2 are 9 and 15 yr at T = 360 K, respectively. Profile 2 was further modelled by a post-diffusion de-gassing on the basis of the bell shape of water abundance on the right side. More modelling details can be found in Methods. All error bars correspond to 2 s.d. analytical uncertainties.

Solar wind origin of water

The rims of CE5 impact glass beads are generally characterized by higher water abundances (up to ~2,000 μg g−1) and lower δD values (~−990‰) compared with pyroclastic glasses, impact-produced agglutinates and melt inclusions reported in previous studies (Fig. 2 and Supplementary Table 4). The extremely low δD values measured at the rims of CE5 impact glasses (Fig. 3) are consistent with the H isotope composition of solar wind40,41. These low δD values are distinct from other potential lunar water reservoirs, such as indigenous water outgassed during lunar volcanism (δD ≈ −300 to 1,200‰) (refs. 22,36,42), water-rich carbonaceous chondrites (δD ≈ −200 to 800‰) (refs. 43,44) and comets (δD ≈0 to 2,500‰) (refs. 45,46), indicating that the water retained in CE5 impact glasses is of a solar wind origin.

There are two potential mechanisms, de-gassing and two-endmember mixing, that could explain the negative correlation between water abundances and δD values defined by CE5 impact glasses (Fig. 2). Because the initial D/H ratio is so low (D/H = 3.1 × 10−6 for the analysis with the highest H2O abundance), the H2 de-gassing modelling indicates that D/H fractionation due to H2 loss would yield an increase in δD values from about −980‰ up to only −926‰ for 99.9% loss, starting with an initial water abundance in CE5 impact glasses of ~2,000 μg g−1 (Fig. 2 and Methods). Therefore, H2 loss from CE5 impact glasses via de-gassing is not a viable mechanism to account for the observed δD variations from around −990‰ to +522‰ (Fig. 2). More important, the high water content at the rim but low water content in the core of the glass beads is not consistent with the de-gassing scenario.

In the other scenario, the CE5 impact glasses could have lost most of the original water inventory of their precursors during their formation at high temperature, following which solar wind-derived water started diffusing inwards into these glass beads on the Moon’s surface. This scenario is consistent with the notable lower abundances of volatile elements such as Na2O (most <0.05 wt%) and K2O (<0.06 wt%) in CE5 impact glasses compared with CE5 mare basalt fragments (Na2O ≈ 0.6 wt% and K2O ≈ 0.1–0.2 wt%, on average)39. This binary mixing model, between a water-poor but relatively high D/H endmember and a water-rich but extremely low D/H endmember, reproduces well the range of water abundances and δD values measured in CE5 impact glasses (Fig. 2). Modelling requires that the initial water abundances of CE5 impact glass beads are less than 50 μg g−1 before the addition of solar wind-derived water (Fig. 2). This requirement is consistent with the water abundances measured at the cores of CE5 impact glasses, which are comparable to the instrument H2O background (~10–30 μg g−1; Methods), indicating that any water potentially present in the precursor materials was lost during formation of the CE5 impact glass beads, before the addition of solar wind-derived water. Most important, the hydration profiles in CE5 impact glass beads show that solar wind-derived water has diffused inwards into the glasses (Fig. 3), unambiguously demonstrating that the addition of solar wind-derived water post-dated the formation of the impact glass beads.

Impact glasses as the driver for the lunar surface water cycle

The negative correlation between water abundances and δD values, coupled with the hydration profiles observed for CE5 impact glass beads, reveal that solar wind-derived water can diffuse into and be stored in lunar impact glasses on the Moon’s surface, probably via post-implantation diffusion19,47. To model the diffusion timescale, we used Raman spectroscopy to confirm that water hosted in CE5 impact glass beads is in the form of hydroxyl/molecular water (Extended Data Fig. 6). The diffusion coefficient of hydroxyl/water in glasses48 was used for modelling the duration timescale necessary for forming the observed hydration profiles of CE5 impact glass beads (Methods). Diffusion modelling results indicate that the diffusion time needed to produce the observed hydration profiles varies from 1 year to 15 years at a peak temperature (T) of 360 K, relevant for the landing site of CE532 (Fig. 3, Extended Data Figs. 4 and 5 and Methods). This short diffusion time indicates that the solar wind-derived water can be rapidly accumulated and stored in lunar impact glass beads. It is interesting to note that one profile on bead CE5#33,036 displays notable post-diffusion water loss at the very rim (Fig. 3, and Methods), suggesting that CE5 impact glass beads have the capacity of both incorporating solar wind-derived water and releasing it to the lunar exosphere if T increases. Therefore, we propose that impact glass beads in lunar soils are a prime water reservoir candidate able to drive the lunar surface water cycle (Fig. 4).

Fig. 4: Schematic diagram of the lunar surface water cycle associated with impact glass beads.
figure 4

a, Because of the high temperature of formation of impact glass beads, most water present in precursor materials would have been lost. b, After their deposition at the surface, solar wind-derived water would diffuse into impact glass beads after solar wind H+ implantation. c, Gardening would transfer the impact glass beads deeper into the soils, creating a water reservoir at the subsurface of the Moon. The impact glass beads at the surface would still be able to release water into the lunar exosphere, due to meteoroid impact, for example13. The amount of water held by impact glass beads in lunar soils is estimated to be 0.03–27.00 × 1013 kg, assuming a global depth of lunar soils of 3–12 m (refs. 56,57), a modal abundance of impact glass beads in lunar soils of 3–5 vol.% (ref. 51) and impact glass beads having a bulk H2O abundance of 132–1,570 μg g−1.

Linking a lunar surface water cycle with impact glass beads implies three major processes: (1) formation of impact glass beads, (2) solar wind-derived water implantation and diffusion into the beads and (3) water release to the lunar exosphere during T increase, possibly due to sunlight irradiation and/or meteoroid impacts (Fig. 4). Impact glass beads are produced during melting of lunar soil and regolith targets when impactors hit the Moon. Most of the water present in the precursor materials would have probably been lost during high-temperature melting that led to the formation of the impact glass beads. Solar wind hydrogen ions would then be continually implanted into lunar soils and impact glass beads and would combine with oxygen atoms to form the hydroxyl/water that diffused into the glass interiors and is stored in glass beads. This water stored in impact glass beads could then be released to the lunar exosphere due to sunlight irradiation19 and/or meteoroid impacts13.

Some of the CE5 impact glass rims contain higher water abundance (up to ~2,000 μg g−1) than that reported in Apollo agglutinates (<600 μg g−1) (Fig. 2); however, they may still not be water saturated compared with terrestrial basaltic glasses of similar chemical compositions that can contain up to 2.43 wt% H2O (ref. 49). This could be the result of dynamic diffusion and release of water in the impact glass beads controlled by the time-of-the-day temperature oscillations10. The dynamic ingress and egress of water in impact glass beads could have acted as a buffer to explain the global and daily variations of water abundance on the lunar surface and in the lunar exosphere7,8,10.

Inventory of water in lunar soils

Using the water abundance measured in CE5 impact glass beads and mare basalts50, we can estimate the relative contribution of solar wind implanted in impact glasses and basalts to the global water inventory in the CE5 regolith. Considering the zoning features of water in CE5 impact glass beads, we estimate that these beads have a bulk water abundance of 132–1,570 μg g−1, which is notably higher than that of the bulk CE5 mare basalts (7 ± 3 μg g−1) (ref. 50). The modal abundance of impact glass beads in lunar soils is in the region of 3–5 vol.% estimated from Apollo lunar soils51. Using this modal abundance yields a contribution by the CE5 impact glass beads to bulk soils of 4–78 μg g−1 equivalent water, which is consistent with bulk water abundance estimates derived from reflectance spectroscopy measurements carried out from the CE5 lander52. Importantly, the major water-bearing phase in CE5 basalts is apatite in which water occurs as hydroxyl occupying the X site in apatite crystal structure53,54. The daily temperature variation of 93 K to 423 K on the Moon’s surface is not high enough to drive hydroxyl out from the apatite structure32,55, indicating that impact glass beads are probably the dominant reservoir involved in the lunar surface water cycle, except for the space-weathered rims and possibly water ice in polar regions.

The global lunar surface is covered by lunar soils with a thickness ranging from about 3 m to 12 m (refs. 56,57). The amount of water in lunar soils contributed by impact glass beads varies from 3.0 × 1011 kg to 2.7 × 1014 kg, based on the bulk water abundance of 132–1,570 μg g−1 in impact glass beads and a modal abundance of ~3–5 vol.% (ref. 51) (Fig. 4). Importantly, this estimate provides a minimum quantity of water available at the lunar subsurface as it does not take into account the water contribution from permanently shadowed regions4,5,6,11, indigenous lunar rocks (for example, ref. 26), delivery by meteoritic impacts13,27 and solar wind-derived water preserved in the space-weathered rims of lunar soils34,35. These findings indicate that the lunar soils contain a much higher amount of solar wind-derived water than previously thought, which could be a water reservoir for in situ utilization in future lunar exploration. Indeed, this water entrapped in impact glass beads appears to be quite easy to extract as shown by the post-diffusion de-gassing profile for bead CE5#33,036 (Fig. 3 and Methods).

Implications for water on airless bodies

Asteroid and comet impacts are the major exogenous processes that reshape the surface morphologies of airless bodies, as evidenced by the widespread presence of impact craters on the Moon58, Mercury59 and asteroids60,61. These impacts probably create impact glasses and glass beads on any airless bodies, as evidenced by the studies of lunar soils62 and howardites63. This study demonstrates that impact glass beads have the capacity to store significant quantities of solar wind-derived water at the surface of airless bodies, in addition to the possible presence of water ice trapped in permanently shadowed areas in polar regions. The presence of water, stored in impact glass beads, is consistent with the remote detection of water at lower-latitude regions of the Moon5,6,7,8,9,10,11,12, 4 Vesta64 and Mercury65,66. Our findings indicate that the impact glasses on the surface of Solar System airless bodies are capable of storing solar wind-derived water and releasing it to space.

Methods

Sample preparation

The CE5 lunar soil allocated by the China National Space Administration and used in this study is CE5C0100YJFM00103, weighing ~1 g. The CE5 lunar soils were scooped by the robotic arm on board the CE5 lander. A total of 150 individual grains, varying in size from ~50 μm to ~1 mm in diameter, were hand picked under a binocular microscope in the ultraclean room at the Institute of Geology and Geophysics, Chinese Academy of Sciences (IGGCAS). The hand-picked grains were mounted in epoxy and prepared as a double polished section with a thickness of about 100 μm (Supplementary Fig. 1 and Supplementary Table 1). The surrounding epoxy was removed to improve the vacuum quality for the water abundance measurements using NanoSIMS67. The prepared section was cleaned using anhydrous ethanol before drying at 60 °C in a baking oven overnight.

SEM observation

Petrographic observations were carried out using field-emission SEM using FEI Nova NanoSEM 450 and Thermofisher Apreo instruments at the IGGCAS, using electron beam currents of 2.0–6.4 nA and an acceleration voltage of 15 kV. The prepared sections were initially coated with carbon for petrographic observations and analyses of mineral chemistry and then recoated with Au for in situ NanoSIMS measurement of water content and hydrogen isotopes. After NanoSIMS analyses, Au coating was removed, and the samples were recoated with carbon and observed by SEM to confirm the positions of the NanoSIMS spots.

Electron probe microanalysis

Before the NanoSIMS analyses, we used the CAMECA SXFive EPMA at IGGCAS to quantify the major and minor elemental abundances in CE5 impact glass beads. The operating accelerating voltage was 20 kV and the beam current was 10 nA, with a 10 μm beam diameter. The total counting time was 10 min for each analysis. One basalt glass (VG-2) and one komatiite glass (MPI-DING-GOR128) were used to monitor the analysis accuracy and precision68. The detection limits were (3σ) 0.02 wt% for Na, Ni, K, S and P, 0.04 wt% for Ca, 0.05 wt% for Al, Mg and Cr, 0.06 wt% for Si and Ti and 0.09 wt% for Mn and Fe. The EPMA data obtained for impact glass beads are listed in Supplementary Table 2.

In situ water abundance and hydrogen isotope analyses

Two analytical sessions were designed to measure the water abundances and hydrogen isotope compositions of CE5 impact glass beads using a CAMECA NanoSIMS 50L at the IGGCAS. In the first session, measurements at the cores and rims of CE5 impact glass beads were carried out to investigate the potential water zonation features in the glasses. In the second session, more-detailed profile measurements on five individual CE5 impact glasses were acquired, as well as several parallel measurements on other glass beads for cross checking. A liquid nitrogen cold trap was used to further improve the vacuum quality in the second session.

The sample and the standard prepared in different mounts were loaded in the same sample holder of NanoSIMS and were baked overnight at ~60 °C in the NanoSIMS airlock. The holders were then stored in the NanoSIMS sample chamber to improve the vacuum quality and minimize the H background49,67. The vacuum pressure in the analysis chamber was 1.5 × 10−10 to 1.9 × 10−10 Torr during analysis. Each 7 μm × 7 μm analysis area was pre-sputtered for 2 min with a Cs+ ion beam current of 2 nA to remove the surface coating and potential contamination. During analysis, the secondary anions 1H, 2D, 12C and 18O were simultaneously counted by electron multipliers (EMs) from the central 5 μm × 5 μm areas using the NanoSIMS electronic gate technique (49% blanking). A 44 ns dead time was corrected for all EMs, while the EM noise (<10−2 cps) was ignored. We used a primary ion beam current of ~0.5 nA for analysis, corresponding to a beam size of ~500 nm in diameter. The charging effect on the sample surfaces was compensated by an electron gun during analysis.

A chip of the nominally anhydrous San Carlos olivine reference, with a reported water content of 1.4 μg g−1 (ref. 69), was used for instrument H2O background corrections, following the relationship H/Obg = (Hcounts – Hbg)/Ocounts and D/Hmeasured = (1 – f) × D/Htrue + f × D/Hbg, where f is the proportion of H emitted from the instrumental background70. In the first session, D/Hbg was 2.83 (±1.98) × 10−4 and Hbg = 815 ± 457 cps (2 s.d., N = 16, corresponding to an instrument H2O background abundance of 29.8 ± 14.8 μg g−1 (2 s.d.)). In the second session, the measured H counts on the central areas of CE5 impact glass beads are even lower than that of San Carlos olivine. Therefore, the extremely low H/O ratio measurements on CE5 impact glass beads were treated as reference (Hbg = 215 ± 40 cps, corresponding to an instrument H2O background of 11.7 ± 2 μg g−1) for background correction.

After an instrument H2O background subtraction, the water abundances of CE5 impact glass beads were calculated from the background-subtracted H/O ratios multiplied by the slope of the calibration lines (Extended Data Fig. 2), which were determined by measuring two apatite standards, Durango apatite (H2O = 0.0478 wt% and δD = −120 ± 5‰) (refs. 26,71) and Kovdor apatite (H2O = 0.98 ± 0.07 wt% and δD = −66 ± 21‰) (ref. 72), the SWIFT mid-ocean-ridge basalt (MORB) glass (H2O = 0.258 wt% and δD = −73 ± 2‰) and two basaltic glasses, 519-4-1 (H2O = 0.17 wt%) (ref. 49) and 1833-11 (H2O = 1.2 wt%) (ref. 49) (Supplementary Table 3). Corrections for instrumental mass fractionation (IMF) on H isotopic compositions were carried out using repeated analyses of the SWIFT MORB glass standard and monitored by analysing both the Durango apatite and Kovdor apatite standards during the whole analytical session (Extended Data Fig. 3). The matrix effects on water abundance and IMF on hydrogen isotope composition are the same between apatite and silicate glass within analytical uncertainties67. Hydrogen isotopic compositions are given using the delta notation, δD = ((D/H)sample/(D/H)SMOW) − 1) × 1,000‰, where SMOW is the standard mean ocean water with a D/H ratio of 1.5576 × 10−4. More technical details can be found in ref. 67. All data are reported with their 2 s.d. analytical uncertainties, which include reproducibility of D/H measurements on the reference materials, uncertainty of H2O background subtraction and internal precision on each analysis (Supplementary Tables 3 and 4). The raw measured D/H ratios were corrected for the background, followed by correction for IMF.

Raman spectroscopy

The Raman spectra of impact glasses were collected using a WITec alpha 300 R confocal Raman system and carried out after the NanoSIMS measurements at the IGGCAS. An optical microscope with an objective of 0.9 NA and ×50 magnification was used to focus the Ar+ laser beam with a 484 nm excitation wavelength on the target phases. Raman spectra were acquired with a total integration time of 500 s, a 9.2 mW laser power and a beam size of 5 μm. Three CE5 impact glass beads CE5#33,076 with water abundances zoning from 1,909 μg g−1 at the rims to <10 μg g−1 in the cores were chosen for Raman analyses (Extended Data Fig. 6). Along with the CE5 impact glass beads, terrestrial MORB glasses, MRN-G1 (H2O = 1.5 wt%) (ref. 73) and EPR-G3 (H2O = 0.22 wt%)73, were also measured to demonstrate the variation of Raman peak intensity of water and/or hydroxyl at ~3,300 cm−1 to 3,700 cm−1 with their water abundance (Extended Data Fig. 6).

Correction of water abundances and D/H ratios for spallation effects

The measured D/H ratios have been corrected for the potential effects of spallation by cosmic ray, using a D production rate of 2.17 × 10−12 mol g–1 Myr–1 (ref. 74) for impact glasses. The correction errors induced by D spallation are around 5% on D/H ratios and negligible on water content36. Because no cosmic-ray exposure age is yet available for the CE5 impact glass beads, a cosmic-ray exposure age of 50 Ma estimated by Hu et al.50 was used for correction of spallation in this study (Supplementary Table 3). Another D production rate, of 9.2 × 10−13 mol g–1 Myr–1 (ref. 75), was also used for correction of spallation effects to account for the suggestion that the D production rate of 2.17 × 10−12 mol g–1 Myr–1 (ref. 74) may not be appropriate for impact glasses. As the D production rate in Ref. 75 is slower, corrections for spallation-produced D are even smaller, not making any notable difference overall for our analytical results (Supplementary Table 4).

De-gassing modelling

The hydrogen isotope fractionation during volatile loss into a vacuum is given by α2 = M1/M2, where M1 and M2 are the masses of the volatile phase isotopologues. The change of the isotopic composition of H during volatile loss by Rayleigh fractionation is given by R = R0 × f(α − 1), where R0 and R are the initial and final D/H ratios for a fraction f of remaining hydrogen76. De-gassing of H2 (M1 = 2 for H2 and M2 = 3 for HD) yields an α value of ~0.8165 (ref. 76) (Fig. 2). The de-gassing modelling assumes that the CE5 impact glass beads have an initial water abundance of 2,000 μg g−1, an initial D/H ratio of 3.1 × 10−6 (corresponding to that of the analysis with the highest H2O abundance) and de-gas in the form of H2, yielding a δD value increase from −980‰ at f = 1 to −926‰ at f = 0.001.

Diffusion modelling

The five sharp hydration profiles of CE5 impact glass beads were simulated by a simplified diffusion process that follows Fick’s second law as equation (1) (ref. 77).

$${{{{C}}}}\left( x \right) = \left( {{{{{C{\mathrm{s}}}}}} - {{{{C}}}}0} \right){{{\mathrm{erfc}}}}\left( {\frac{x}{{2\sqrt {Dt} }}} \right) + {{{{C}}}}0$$
(1)

In equation 1, D is the diffusion coefficient that depends on temperature, x is the distance to the rim, C0 is the initial water abundance of the CE5 impact glass beads, C(x) is the water abundance at distance x, Cs is the maximum water abundance at the rims, t is duration of diffusion and erfc is the complementary error function.

Considering the scooped samples were collected from less than 3 cm deep38, diffusion simulations were performed using a peak temperature of 360 K, relevant for a latitude of 45° N on the Moon32, because the peak temperature is the main controlling factor for the diffusion time35. The H diffusion coefficients of water/hydroxyl (DH) determined on silica glass48, DH = 20.84 μm2 yr−1 at 360 K, were used in our simulations. It is not possible to reconstruct the orientation of the exposed surface of the studied impact glass beads on the lunar surface. Thus, we carried out a simplistic diffusion modelling based on the measured profiles, assuming a constant Cs at the rims for each bead based on their diverse Pb–Pb ages39. The diffusion duration is roughly a maximum constraint for qualitatively estimating the efficiency of solar wind-derived OH/H2O recharged in impact glass beads. The water abundances at the rims (Cs) of the six CE5 impact glass bead profiles are extrapolated from the measured profiles using 2,000 μg g−1, 2,000 μg g−1, 1,600 μg g−1, 500 μg g−1, 2,600 μg g−1 and 3,000 μg g−1 and a C0 of 10 μg g−1 for Profiles 1–6, respectively. The diffusional durations vary from 3.0 to 9.0 years, 5.0 to 15.0 years, 1.0 to 4.5 years, 2.0 to 6.0 years, 0.5 to 4.0 years and 2.0 to 14.0 years at T = 360 K for Profiles 1–6, respectively, based on the apparent distance (Fig. 1) and a fixed Cs at the rims for each profile.

For the notable post-diffusion de-gassing feature of Profile 2, a two-stage modelling of water diffusion (Stage 1) and post-diffusion de-gassing (Stage 2) was carried out to fit the measured water abundance profile. The water-diffusion details in Stage 1 have been described in the preceding. A post-diffusion de-gassing of water (Stage 2) was modelled to fit the bell shape of water abundance on the right side (Fig. 3). The water abundance on the right side modelled in Stage 1 was taken as the initial condition for modelling of Stage 2. The de-gassing rate is assumed to be proportional to the surface water abundance along the profile, where transmission capacity of a boundary was formulated as an additional boundary condition for the diffusion equation (see ref. 35 for more details). The corresponding δD profiles are calculated by a binary mixing model (Two-endmember mixing modelling) considering a limited diffusional fractionation for H isotope78 compared with the large difference in δD values for the two endmembers (δD ≈ −980‰, D/H ≈ 3.0 × 10−6 for solar wind-derived water and δD ≈ 500‰, D/H ≈ 230 × 10−6 for the initial water in impact glass beads).

Two-endmember mixing modelling

In the two-endmember mixing model, the solar wind endmember is set to have a δD value of −980‰ (ref. 41) and a water abundance of 2,000 μg g−1, which corresponds to the maximum measured value in this work (Fig. 2). The initial water component in CE5 impact glass beads could have a minimum δD value of ~500‰ based on the measurements on CE5 apatite50, yielding a best fit of initial water abundance of 5–50 μg g−1 before the addition of solar wind-derived water (Fig. 2). On the basis of the following mass balance (equations (2) and (3)), we can model the relationship between the measured water abundances and δD values:

$${{{\mathrm{H}}}}_2{{{\mathrm{O}}}}_0 + {{{\mathrm{H}}}}_2{{{\mathrm{O}}}}_{{{{\mathrm{solar}}}}} = {{{\mathrm{H}}}}_2{{{\mathrm{O}}}}_{{{{\mathrm{measured}}}}}$$
(2)
$${{{\mathrm{H}}}}_2{{{\mathrm{O}}}}_0 \times {{{\mathrm{D}}}}/{{{\mathrm{H}}}}_0 + {{{\mathrm{H}}}}_2{{{\mathrm{O}}}}_{{{{\mathrm{solar}}}}} \times {{{\mathrm{D}}}}/{{{\mathrm{H}}}}_{{{{\mathrm{solar}}}}} = {{{\mathrm{H}}}}_2{{{\mathrm{O}}}}_{{{{\mathrm{measured}}}}} \times {{{\mathrm{D}}}}/{{{\mathrm{H}}}}_{{{{\mathrm{measured}}}}}$$
(3)

where H2O0 and D/H0 are the initial water abundance and D/H ratio of CE5 impact glass beads before the addition of solar wind-derived water, respectively, H2Osolar is the water abundance contributed by solar wind-derived water, D/Hsolar is the hydrogen isotopic composition of solar wind-derived water, and H2Omeasured and D/Hmeasured are the water abundance and D/H ratio measured from CE5 impact glass beads, respectively.

Petrography and chemistry of CE5 impact glass beads

Petrography

A total of 150 impact glass fragments were mounted in CE5#33 mount (Supplementary Fig. 1 and Supplementary Table 1). There are 117 glass fragments displaying ovoid to spherical shapes, which are mainly composed mainly of glass, 1 glass rod, 3 olivine and plagioclase monomict clasts, 1 basaltic fragment and 28 glass fragments showing irregular shapes, which are composed mainly of glass, mafic mineral fragments and large bubbles (Supplementary Fig. 1 and Supplementary Table 1). The ovoid to spherical fragments are referred to as glass beads in this work and further classified into two textural types, homogeneous and heterogeneous glass beads (Fig. 1, Supplementary Fig. 2 and Supplementary Table 1). The homogeneous glass beads display smooth surfaces and homogeneous contrast in back-scattered electron images while the heterogeneous glass beads usually contain some vesicles, pyroxene and plagioclase clasts and micrometre-sized sulfide spheres (Fig. 1 and Supplementary Fig. 2). The textural differences between homogeneous and heterogeneous glass beads could reflect a higher formation temperature for the homogeneous glass beads (Fig. 1 and Supplementary Fig. 2).

Chemistry

A total of 32 homogeneous impact glass beads were selected for EPMA analyses. The analytical results are listed in Supplementary Table 2. The analysed glass beads have a mare origin, except for two glass beads plotting in the highlands population and another two distributed in pyroclastic areas judging from their CaO/Al2O3 and MgO/Al2O3 ratios (Extended Data Fig. 1)79,80. These glasses plotting in non-mare fields are not discussed further in this work because their origin is currently uncertain. The homogeneous CE5 impact glass beads are composed mainly of SiO2 (31.2–44.5 wt%), TiO2 (3.95–6.54 wt%), Al2O3 (10.4–17.3 wt%), FeO (17.6–23.5 wt%), MgO (5.94–10.9 wt%) and CaO (10.1–15.4 wt%), with minor Cr2O3 (0.12–0.28 wt%), MnO (0.22–0.31 wt%), NiO (0.15–0.61 wt%), Na2O (<0.17 wt%, most <0.05 wt%) and K2O (<0.06 wt%) (Supplementary Table 2). Each single glass bead is homogeneous in terms of its major element chemical composition (Supplementary Table 2). Most CE5 homogeneous impact glass beads have compositions similar to the bulk CE5 basalt compositions reported by refs. 38,81, indicating that they could have been formed locally rather than ejected from other regions. This prediction is consistent with the similar geochemical characteristics between CE5 impact glass beads and CE5 bulk basalts82. The Pb–Pb dating of the CE5 impact glass beads indicates that they were formed in the range of 4–2,000 Ma by various local meteoritic impacts39.

Water abundances and hydrogen isotope compositions of the CE5 impact glass beads

A total of 32 CE5 impact glass beads with smooth exposed surface (Fig. 1 and Supplementary Fig. 2) and a mare origin (Extended Data Fig. 2) were selected for the NanoSIMS measurements of water abundances and hydrogen isotopic compositions. In the first analytical session, the objective was to investigate the homogeneity of water abundance within each glass bead; therefore, two analyses of the cores and the rims on relatively large glass beads were carried out to identify the hydration and or de-gassing candidate profiles. In the second session, six more-detailed profile analyses were carried out on five CE5 impact glass beads (Fig. 1), as well as several duplicate measurements on other glass beads (Supplementary Fig. 2). Overall, the analytical water abundances in the first session are in the range of the second session due to the heterogeneous distribution of water in the CE5 impact glass beads (Fig. 2 and Supplementary Table 4). The water abundances of the analysed CE5 impact glass beads vary from ~0 to 1,909 μg g−1 and are negatively correlated with their δD values (−990 ± 6‰ to +522 ± 440‰) (Fig. 2 and Supplementary Table 4). Water abundances measured in all CE5 impact glass beads display higher values at the rims and decrease towards the glass bead cores (Fig. 3 and Extended Data Figs. 4 and 5). The δD values display the reverse trends (Fig. 3 and Extended Data Figs. 4 and 5). Three profile measurements (Profile 3, Profile 4 and Profile 5) show a sharp decrease in water abundance from the rim towards the core (Fig. 3 and Extended Data Figs. 4 and 5), while the other three (Profile 1, Profile 2 and Profile 6) display a modest decrease of water abundance towards the cores (Fig. 3 and Extended Data Figs. 4 and 5). It is also noticeable that the latter profiles have higher water abundances (1,055–1,909 μg g−1) at the rim than do the former profiles (176–1,423 μg g−1) (Fig. 3, Extended Data Figs. 4 and 5 and Supplementary Table 4). All of the profile measurements acquired on the central areas of CE5 impact glass beads have water abundances comparable to the instrument H2O background (Supplementary Tables 3 and 4).