Potassium map from Chang'E-2 constraints the impact of Crisium and Orientale basin on the Moon.

KREEP materials were thought to be last crystallized at the lunar crust and mantle boundary. Impact cratering and volcanism are mainly responsible for their distributions on the lunar surface. Therefore, observation of global KREEP materials and investigation of distributions in the areas of large basins are of critical importance to understand the geologic history of the Moon. Here we report the new global potassium distribution on the Moon detected by Chang'E-2 Gamma-ray Spectrometer. We found that our new measurements are in general agreement with previous observation. A new finding and an important difference is that relatively higher K abundances in the Mare Crisium and Mare Orientale than their surrounding rims were detected for the first time. In light of our observations in these two areas, we propose that Crisium and Orientale basin-forming impact events may have penetrated to the lower crust and excavate the deeper materials to the lunar surface.

In the Lunar Magma Ocean (LMO) hypotheses1, KREEP materials are proposed to lie between the mafic mantle and plagioclase crust in the process of compositional crystallization. Their presence on the lunar surface is the end product of a series of evolutionary processes that have changed the surface features of the Moon2. In light of LMO hypotheses and the general impact cratering theory3, impact events penetrating through the crust should excavate the KREEP materials on the lunar surface. Therefore, observation of the global KREEP material on the Moon and the investigation of KREEP material in the areas of large basins are of critical importance to understand the formation of such a magma ocean and constrain the depth of the huge impact events.

In the past few years, several spacecraft with improved compositional remote sensing instruments were sent into the orbits around the Moon and large quantities of valuable data have been obtained. Chang'E-2 spacecraft, the second lunar probe of China, was launched on 1 October 2010, and inserted into a circular orbit of 100 km altitude on 15 October 2010. It carried a suite of remote sensing and plasma instruments. Here we report the results from the gamma ray spectrometer during its 178-days measurements and their new implications on the impact depth of the Crisium and of Orientale basin, respectively.

CE-2 GRS is the first gamma ray spectrometer that used large LaBr3 crystal as the main detector to measure gamma rays emitted from the Moon within the energy region from 0.3-MeV to 9-MeV45. The energy resolution of the flight model is ~ 5% full width at half maximum at 662 keV, which is about two times better than LPGRS6. Since CE-2 GRS was designed to be mounted inside the spacecraft, a big well-shape CsI(Tl) crystal was employed in anticoincidence with the main crystal to suppress the Compton effects and the background gamma rays from the spacecraft materials generated by Galactic Cosmic Rays (GCRs). Gamma rays emitted from the lunar surface are mainly from natural radioactive elemental (e.g., U, Th, K) decay effects and major elemental (e.g., Si, O, Fe, Ti, Mg, Al, Ca) excitation from interactions with incident GCRs7. Gamma rays entering the main crystal were recorded by the CE-2 GRS with a cadence of 3 second in the 512-channel spectrum.

More than 25 full energy peaks could be identified in CE-2 GRS spectra corresponding to different elements on the lunar surface (Fig. 1). But owing to the limitation of intrinsic energy resolution of the scintillation detector, many spectral peaks contain contributions from two or more elements. The strongest peak around 1.46 MeV is a combination of lines resulting mainly from the K line at 1.46 MeV and LaBr3 self-activity because of the 138La decay (1.436 MeV + 32 keV)8 that contributes about 85% of the corresponding gamma-ray intensity.

Figure 1

Chang'E-2 gamma-ray spectrum plotted as counts per second (cps) of the whole Moon versus energy for its 178-days measurements with 100 km altitude.

The major peaks are labeled. The spectrum below 2.7 MeV are mainly contributed by the self-activities from LaBr3 crystal: natural 138La is primarily responsible for the spectrum below 1.7 MeV while 227Ac decay chain generates alpha emitters between 1.7 and 2.7 MeV, as label 3. K line (1.46 MeV) is partly contributed to the sum peak as Label 2 with the 138La decay to 138Ba (1.436 MeV + 32 keV) of the LaBr3 crystal; label 1 is the sum peak of elements from the Moon, and crystal beta decay of 138La (789 keV) to 138Ce. The inset figures are the fitting results of the energy region around 1.46 MeV. (A) The orbital spectrum and its fitting result. (B) The cruise spectrum and its fitting result.

Corrections of gain, solid angle, and cosmic ray variation have been carried out to CE-2 GRS measurements. The spatial resolution of CE-2 GRS is about 150 km in its 100 km working orbit, but in this work, all the corrected spectra are binned into equal-area pixels about 60 km × 60 km on the Moon. In our calculation we assumed gamma rays emitted from LaBr3 self-activity decay are constant in the whole mission, the spectra measured in the cruise period could therefore be used as the background. The net 40K gamma-rays emanating from the Moon in each pixel are obtained by calculating the peak fitting difference between binned spectra and the cruise spectrum (see inset of Fig. 1). Then all the net 40K count rates from 11306 pixels were compared to count rates in the detector derived from the calculated surfaced gamma-ray fluxes to determine the abundance7910.

Results

The K distribution from CE-2 GRS (Fig. 2) is in general agreement with Lunar Prospector91112, Kaguya10, and Chang'E-1 measurements13. The global average abundance is 620 ± 615 ppm with the upper limit of 3240 ppm, which is lower than the LP value of 755 ± 740 ppm with the upper limit of 4360 ppm9. The Procellarum KREEP Terrane (PKT)14 has an average value of 1768 ± 530 ppm, comparing to 2004 ± 708 ppm determined by LPGRS. In the South Pole-Aitken (SPA) basin, the abundance ranges from 500 ppm to 1350 ppm with an average value of 615 ± 235 ppm, which is also lower than the LP measurements of 765 ± 320 ppm. Different from the PKT and SPA basin, lunar highlands is depleted in the K with the abundance being less than 500 ppm according to the CE-2 GRS observations. This is consistent with LP measurements9.

Figure 2

The global K abundance (unit: ppm) of the Moon from CE-2 GRS 178-days measurements.

All the data are binned into quasi-equals 2° × 2° pixels in the equator and displayed in a cylindrical projection.

Mare basins on the nearside, such as Serenitatis, Tranquillitatis, Fecunditatis, and Crisium, have lower K abundance than mare regions. But after careful examination of the CE-2 GRS data, we found that the abundance within the Crisium impact basin is clearly higher than its surrounding highlands (Fig. 3a–b). This is inconsistent with previous LP results (Fig. 3c), which indicates that there is no distinct difference between the K abundance of Mare Crisium and that of its surrounding11.

Figure 3

Surface topography and K distributions of the 550-km-diameter Crisium Basin.

All the maps are displayed in a Lambert azimuthal equal-area projection superimposed on the shaded lunar surface relief. Two circles represent the proposed rims with diameter of 300 km and 700 km. NE - northeast highlands of Crisium basin; SW - southwest of Crisium basin that covered by ancient mare basalt. Picard - Picard crater with 30 km in diameter. (A) the topographic elevation of the Crisium basin from Lunar Orbiter Laser Altimeter32. (B) K abundance of the Crisium Basin measured by CE-2 GRS with data binned into quasi-equal 2° × 2° pixels. (C) K abundance of the Crisium basin measured by LPGRS9.

Discussion

The Crisium basin is a Nectarian-age multi-ring impact structure on the eastern nearside of the Moon. Its interior is covered by mare basalt from a protracted period of volcanism over 0.8 Ga15. The abundance in its inner mare basalt ranges from 430 ppm to 840 ppm with an average value of 600 ±100 ppm. On the other hand, the K abundance in the surrounding highlands is lower than 480 ppm with an average value of 290 ± 105 ppm. Because the southwest highlands is thought to be covered with ancient mare deposits16, not with ejecta from the Crisium basin, the abundance in its northeast highlands thus represents the true value of the ejecta for estimating the average value here (Fig. 3b). The highest K abundance in the Picard crater (14.6°N, 54.7°E) indicates the subjacent material has K abundance about 3 times (900 ppm) higher than highlands deposits surrounding the Crisium basin.

There is a long standing question concerning the origin of material excavated by the Crisium impact event1617. Some studies indicated that the Crisium impact event only excavated the upper or middle crustal material with a depth of about 40 km16 while others preferred that the impact event might have excavated mantle material17. According to recent observations of olivine sites around the Crisium basin17, the scenario that impact event only excavate upper or middle crustal material is probably excluded because olivine is a significant component associated with lunar mantle or Mg-rich plutons in lower crust18. Furthermore, if the Crisium impact event did penetrate through the crust and KREEP layer, exposing the upper mantle, KREEP (e.g., K here) and mafic materials should be found at its rim in light of the general impact cratering theory3, deep material will be spread around the near rim. However, higher K signal is not observed in Crisium rim and its basin terrain by CE-2 GRS, FeO is also found as lower as 2–4 wt.% in these areas19.

On the basis of our observation and previous results161719, we prefer the scenario that the Crisium impact event only penetrated to the plutonic lithosphere resulting from the magmatic intrusion into the lower crust. The magmatic intrusion, which is enriched in the PKT, is the hybrid source with low-density Mg-rich olivine at the base crust resulting from the overturn of the early mantle cumulate, KREEP, and crustal anorthosite20. While the Crisium impact event penetrated through the crust, to the upper boundary of the magmatic intrusion, the deepest materials would be ejected to the surrounding rim area. The redistribution of excavated material resulted in the lower K and FeO in the highlands19 bounding the basin and is in agreement with the observations of olivine exposure on the basin rim17. The constraint of impact to the upper boundary of the magmatic intrusion is consistent with the higher K abundance in subjacent highlands material covered by mare basalt as detected in the Picard crater, and is also consistent with the non-zero crustal thickness in the Crisium basin21. The big basin formation impact excavated tremendous volume material and formed a thinner crust with fracturing underlying lithosphere22 in the Crisium area. Owing to the density difference between lower crust and partly melted magma, the buoyancy contrast23 enabled the eruption of basalts at the Crisium basin. The partly melted magma is from the Moon's deep interior and produced by slowly heating up from the decay of radioactive elements. This supports the observations of higher K abundance (about 600 ppm) and very higher FeO contents19 (up to 20 wt.% for areas of the central areas) in the mare basalts covered in the Crisium basin. Located near the PKT region and consistent K abundance with mare basalts in the PKT boundary, high-heat production by radioactive decay is suggested to be responsible to the long duration of volcanic extrusions at the Crisium basin16.

Similar characteristics are also found in the Orientale basin (Fig. 4a–d), the youngest multi-ring basin on the Moon24. CE-2 GRS shows that the central areas (5–6 pixels in diameter with each pixel of ~ 60 km) of the Orientale basin corresponding to the Mare Orientale, the largest mare deposit (~ 320 km in diameter) in the basin, has obviously high K abundance with the average of 360 ± 115 ppm (Fig. 4b). This value is about a factor of two higher than the surrounding areas, including the Montes Rook Formation (MRF) occurring mostly between the edge of the Inner Depression (ID) and Cordillera Ring (CR), and Hevelius Formation (HF) lying outside the CR. The MRF, consisting of smooth plains and corrugated deposits, has K abundance of 115 ± 95 ppm (Fig. 4b), which is only one-third of that in the Mare Orientale. Most areas of HF, which represent the radial ejecta deposits from the Orientale basin, have an average abundance a little higher (about 80 ppm) than the deposits in MRF, except two areas located at the southwest (K-1) and northwest (K-2) portion of the basin (Fig. 4d). These two areas, related to the well-preserved ejecta units in MRF (part of K-1) and HF25, have abundance of 270 ± 80 ppm (K-1) and 255 ± 90 ppm (K-2), respectively. Because turbulent impact modification after Orientale impact event has less influence on these two areas26, the K abundances in these two areas therefore represent the original value of excavated ejecta from the basin-forming impact.

Figure 4

Orientale Basin and the K distribution from CE-2 GRS and LPGRS in cylindrical equidistant projection superimposed on the shaded lunar surface relief.

(A) the topographic elevation of the Orientale basin from Lunar Orbiter Laser Altimeter32 with the black circles represent the rings of the Orientale basin (CR, ORR), but the Inner Rook Ring is not shown here; CR - Cordillera Ring; ORR - Out Rook Ring; ID - Inner Depression; MRF - Montes Rook Formation; HF - Hevelius Formation; M-R - Mendel-Rydberg basin. (B) K abundance in the Orientale Basin measured by CE-2 GRS with data binned into quasi-equal 2° × 2° pixels. (C) K distribution in the Orientale basin measured by LPGRS9. (D) the smoothed K distribution from CE-2 GRS measurements; K-1/2 represent two areas with relatively higher K Abundances outside the CR. The yellow lines represent the contours of K abundance in K-1/2 area.

Assuming the Out Rook Ring is the closest approximation to the basin excavation cavity, CR is believed to have been formed by the inward collapse at the edge of structural uplift along the basin of the displaced zone25. MRF therefore represents the modification of rim ejecta and the mixture with the upper crust, which could be consistent with the similar K abundance in HF and MRF. The relative higher K abundance (about 80 ppm) in HF implies that the deeper material excavated by the impact event has relatively higher K abundance. The LMO scenario proposes that the KREEP layer lies between the anorthositic lunar crust and mafic cumulated mantle in the long-time compositional differentiation1. Therefore, in contrast to previous works27 which suggested that Orientale impact event excavated only to the upper half of the lunar crust in its target region with the depth of 40–50 km, we prefer that the Orientale impact event penetrated into a source of radioactive-rich material at the lower crust. This scenario supports the lower FeO (2–4 wt.%)19 in the surrounding of the Orientale basin and is also consistent with the estimated crustal thickness in this area21. This scenario is also consistent with the enrichment of low-Ca pyroxene in the central peaks of the Maunder crater, a post-Orientale crater with 55 km diameter in the interior depression of the Orientale basin, which is thought to be sampling the noritic lower crust26.

The internal accumulation of radiogenic heat sources causing partial melting of the mantle is the basic postulation for the formation of mare basalt28. The details of the generation, ascent and eruption of the mare basalts in the Orientale basin need further investigations, but the similar K abundance in the Mare Orientale to the nearside mare basalts, especially to the boundary of PKT, suggests that mare deposits within Orientale basin are similar in composition to nearside mare basalts27. Furthermore, two areas with higher K abundance in the Mare Orientale definitely correspond to the eastern and northwestern region of the inner depression covered by mare basalts, which have almost the identical age, ~ 3.2–3.3 Ga2930. The same concentration of FeO, TiO219, and K in these two areas supports the interpretation that the covered basalts in the center of Orientale basin are from single volcanic events2931.

Methods

All raw data received from the spacecraft were converted to time series gamma ray spectra that were used to derive elemental maps. However, in order to ensure the correctness of spectral analysis, mapping, and abundance calculations, data selection was carried out at first. We eliminated the data which have bad sync words, bad check sums, and off-limit values. The measurements during solar energetic particle event were also excluded. In this work, the energetic proton data taken from the Earth orbiting GOES satellites33 were used to identify the solar energetic particle events.

Three major corrections to the time series data were made for the variations of gain, galactic cosmic ray, and solid angle before analyzing the spectra. To monitor the gain changes, we calculated the peak positions of two strong peaks (0.511 MeV, 7.60 MeV) for every 15-min and every 3-hour interval. The time series data were then corrected according to the gain and offset values for each 3-hour and 15-min interval34. The time-varying galactic cosmic ray flux was monitored by measuring the counting rate of oxygen gamma ray line (6.13 MeV) for all latitudes poleward of + 85°. In order to apply the galactic cosmic ray correction, the 6.13 MeV counting rate were smoothed and normalized to the counting rate seen at the beginning of the orbital mission. All spectra were then divided by this normalized flux for the correction of the galactic cosmic ray variations34. In addition, all spectra were normalized to a spacecraft height of 100 km to normalize the variations of solid angle seen by the detector due to height variation from 85 to 125 km.

All corrected spectra were binned into approximately equal-area pixels about 60 km × 60 km by using their latitude and longitude positions at the times they were taken. Gamma-ray spectra were accumulated in each pixel to present the spectrum measured from such area. To calculate the K counting rate in each pixel, the integrated area beneath K peak should be determined. In this work, peak fitting method (Gaussian shaped, linear background line) was used to determine K peak area in each pixel. Since the peak around 1.46 MeV is a combination of lines resulting mainly from the K line at 1.46 MeV and LaBr3 self-activity because of the 138La decay (1.436 MeV + 32 keV)8, the contribution of LaBr3 self-activity should be eliminated from each calculation. We assumed that gamma rays emitted from LaBr3 self-activity decay were constant in the whole mission. The peak area estimated from the accumulated spectrum in the cruise period was used to eliminate the contribution of LaBr3 self-activity. The net 40K gamma-rays emanating from the Moon in each pixel were obtained by calculating the area difference between binned spectra and the cruise spectrum.

Author Contributions

M.Z. drafted the manuscript. J.C. is the principle investigator of the GRS, and contributed to the data process and analysis with M.Z. W.I. and W.F. contributed scientific interpretations with M.Z. J.C., T.M., J.W., Y.G., M.C. and Y.H. were responsible for development and operation of the GRS system. All of the authors, including A.X. and Z.T. discussed and provided significant comments on the results and manuscript.