To improve the lunar DEM accuracy derived from CE-1 altimeter data, CE-1 laser altimeter data are calibrated in this paper. Orbit accuracy and ranging accuracy are the two most important factors to affect the application of altimeter data in the lunar topography. An empirical method is proposed to calibrate CE-1 altimeter data, using gridded LOLA DEM to correct systematic errors of CE-1 altimeter data, and the systematic bias is about -139.52 m. A new lunar DEM grid model based on calibrated CE-1 altimeter data with the spatial resolution of 0.0625° × 0.0625° is obtained as well as a spherical harmonic model at 1400th order. Furthermore, the DEM accuracy is assessed through the comparison with the nearside landmarks of the Moon, and the results show that the DEM accuracy is improved from 127.3 m to 48.7 m after the calibration of laser altimeter data.
National Natural Science Foundation of China1147305611403076Science and Technology Commission of Shanghai12DZ2273300Chinese Academy of SciencesLunar Exploration Project of ChinaKey Laboratory of Space Object Measurement1. Introduction
Lunar topography plays an important role in lunar exploration [1]. Benefiting from several lunar explorers launched in recent 20 years, humans have reached a high level to understand the lunar topography. After Clementine was launched in 1994, Japan, China, and India launched lunar explorers SELENE (Selenological and Engineering Explorer), CE-1 (ChangE-1), and Chandrayaan-1, respectively, between the year of 2007 and 2008. In 2009 the US launched LRO (Lunar Reconnaissance Orbiter), which provided accurate laser altimeter data and CCD image data for scientists to study the lunar topography [2–7].
With the increasing topography measurements, the study of lunar topography no longer depends on one single probe or one single type of data. Data from different devices of the same probe can be used to verify each other and improve the data quality. For example, the combined adjustment of CE-1 laser altimeter data and CCD image data can improve the spatial resolution of a topography model effectively [8, 9]. Barker solved lunar topography by combining data from LRO Lunar Orbiter Laser Altimeter (LOLA) and TC (Terrain Camera) deployed on SELENE’s main satellite Kaguya, and the results showed that TC data can make up for gaps in the spatial distribution of LOLA and improve the resolution of the topography model [10].
LAM (laser altimeter) deployed on CE-1 obtained about total 9.12 million measurements covering the whole lunar surface. Ping et al. developed a lunar global digital elevation model (DEM) CLTM-s01 with 360th order based on the 3 million effective laser altimetry observations acquired during the first official flight of CE-1. The space resolution was better than 0.5°, and the elevation accuracy of the model was about 31 m [11]. Li et al. developed a DEM with 3 km resolution. The model had a planar positioning accuracy of 445 m and an elevation accuracy of 60 m [12]. Cai et al. proposed a hierarchical many-knot spline method and developed a global lunar DEM with a 0.0625° resolution based on more than 8.2 million data obtained from CE-1 [13]. Hu reprocessed LAM data by using the crossover adjustment method. The lunar surface was divided into 32 regions for the least square adjustment, and the newly obtained topography model was improved significantly on the regional scale [14].
The orbit accuracy and ranging accuracy are the two most important factors to affect the application of altimeter data in the lunar topography. The altimeter measurements include systematic errors and random errors, and the systematic errors should be calibrated carefully before the lunar topography calculation. The mean radius of the Moon derived from Chang’E-1 is 1737013±2 m [11], which has an obvious bias of nearly 140 m when compared with those derived from LRO (1737153±10 m) and LALT (1737150±10 m) [2, 15], and indicates that CE-1 LAM data has a systematic bias and should be removed before DEM calculation. Ping et al. also estimated the possible systematic bias of reference frequency in CE-1 LAM oscillator, and the result is about 145 m [16]. In this paper, a different method was used to calibrate CE-1 LAM data systematic errors, using LOLA data as the benchmark to calibrate and remove systematic errors of the data from CE-1 LAM data, and a new lunar DEM was generated using calibrated LAM data.
2. Orbit Improvement and Calibration of LAM Data
CE-1 was tracked by Unified S-Band (USB) ranging and Very Long Baseline Interferometry (VLBI) technique. To improve the orbit accuracy of CE-1, the altimeter crossover data were used in the orbit determination. We reprocessed the CE-1 orbit by combining the ground tracking data as well as the altimetry crossover data, and the orbit accuracy was increased from 147 m to 105 m, which was about ~30% improvement compared with the result using ground tracking data only [17].
It is important to eliminate altimeter data errors by cross-correcting. The altimeter data contain systematic errors consisted of orbit errors, positioning errors of laser footprints on lunar surface, time tag, and instrument errors. In this paper, an empirical data calibration method was used to remove the systematic errors remaining in LAM data.
As Table 1 shows, LOLA ranging error was less than 1 m, which was better than LAM’s 5 m accuracy and LALT’s 4 m accuracy, and the orbit accuracy of LRO was also much better than that of CE-1 and SELENE. Therefore, the LOLA data were used as the benchmark to calibrate possible systematic errors in LAM data before the DEM calculation.
Specifications of altimetry instrument on CE-1/SELENE/LRO.
CE-1
SELENE
LRO
Orbit altitude
200 km
100 km
50 km
Instrument
Laser AltiMeter
Laser ALTtimeter
Lunar Orbiter Laser Altimeter
Ranging error
5 m
4.1 m
<1 m
Sampling rate
1 Hz
1 Hz
28 Hz
Radial orbit error
15 m
1 m
1 m
3D orbit error
~200 m
~50 m
~20 m
Footprint
<200 m
~40 m
5 m/beam
In this paper, the CE-1 LAM observations, namely, LAM 2B-SHAO data, were about 7 million, in about 5 months (2007.11.28–2008.02.06 and 2008.05.15–2008.07.07), which were obtained based on the improved orbit [17]. The LOLA data used here were about 1 year during the LRO mapping orbit from 2009.9.26 to 2010.8.20 [15], which were the initial observations from LRO, and LOLA data acquisition continued up to 2015 [18].
The LOLA data were firstly used to develop a reliable topography model, which would be used to calibrate LAM data. The following model was used to calculate systematic errors in LAM data.(1)Δht=ht-hM=b+h˙τ+∑k=1lSksinωkt-t0+Ckcosωkt-t0
In formula (1), h(t) was the elevation derived from LAM data at the time t. hM was the elevation of the same point as h(t) calculated from the LOLA topography model. And Δh(t) was the elevation residual. h˙ was the rate of height change. We modeled the elevation residual Δh(t): b was the constant bias, τ was the time tag deviation, ωk was the frequency, Sk and Ck were used to characterize the error relating to the orbit, and t0 was the initial epoch time. b, τ, Sk, and Ck were the parameters to be estimated using the least square method. The value of ωk could be obtained by doing Fourier transform on Δh(t).
Figure 1 showed the time series of the residual Δh(t) of LAM data. In the calculation of hM, nearest neighbor interpolation was used. In order to determine the characteristic frequency, we needed to carry out Fourier spectrum analysis on Δh(t), and the results were shown in Figure 2. Judging from the spectral analysis in Figure 2, the main periodic terms in the elevation residuals were ~129 min (about 1 orbit period) for LAM data. Because the orbit error was dominated by the error at a frequency of once per orbital revolution (1 cpr), only 1 cpr orbital error was estimated (i.e., l=1 in (1)).
Residuals of LAM data compared with LOLA topography model.
Spectral analysis of LAM residuals.
The calibration results were shown in Table 2. For comparison, the calibration result of SELENE LALT data was also shown in Table 2. The result indicated that errors in LALT data were smaller than those in LAM data. The constant bias of LAM was -139.52m, which was the same level as Ping’s analysis of the overall deviation of ground crystal oscillator [16]. The time tag deviations in LAM and LALT were very small. The orbital error based on LAM data was about 5m, and the orbital error based on LALT data was less than 1m, both of which were smaller than the real radial orbit errors of CE-1 and SELENE, respectively, indicating that the method cannot remove orbital errors totally. We used the above estimated results to correct LAM data, and the corrected data were adopted in the DEM calculation. For the sake of clarity, the modified LAM 2B-SHAO data were named as LAM 2B-SHAO-cali.
Results of LAM and LALT data calibration.
b/m
τ/s
Sk/m
Ck/m
LAM
-139.52
-0.0073
5.25
-4.21
LALT
3.03
-0.0780
0.26
-0.24
3. The Lunar DEM Based on the Calibrated CE-1 LAM Data
After LAM data systematic errors correction, we developed a lunar DEM based on the calibrated CE-1 altimeter data with a 0.0625° resolution. In this paper, we used a quadratic fitting method to interpolate grid points. The formulae were as follows:(2)hλi,φi=Ai·ξ(3)Ai=1λi-λ0φi-φ0λi-λ0φi-φ0λi-λ02φi-φ02(4)ξ=a0a1a2a3a4a5T
where (λi,φi) and h were the lunar longitudes, latitudes, and elevation of the laser altimetry footprints in the grid, respectively. (λ0,φ0) were the longitude and latitude of the grid point to be interpolated, and a0~a5 were 6 parameters to be estimated. The fitting result a0 was the elevation value of the grid points to be interpolated.
Because CE-1 was a polar-orbiting satellite, the altimeter data were mainly distributed along the longitude. In order to avoid the abnormal solution as a result of few or poor distribution of the altimeter data, when the correlation coefficient of the fitting result was too large (greater than 105) or the measurement in the grid was less than 50, we applied the inverse distance weighted interpolation method whose distance weight function followed Gaussian attenuation to calculate the elevation of the grid point:(5)h^=∑j=1nwjhj∑j=1nwj(6)wj=exp-djD2cosφjwhere dj was the distance of hj to the grid center. D was the value of grid width at 0.0625° in this paper.
Figure 3 showed the results of the DEM, in which the elevation value referred to the mean lunar radius of 1737.4 km. Figure 3 offered a global Mollweide projection, and the central longitude was 90°W.
Lunar topography model based on the calibrated LAM data.
Three topography models based on the original LAM data, LAM data with orbit improvement, and LAM data with calibration, respectively, were compared to the topography model derived from LOLA data. Figure 4 showed the differences. In Figure 4(a), stripped differences could be seen clearly, which were mainly caused by the orbit error. With the new orbit results, no more significant stripped differences could be found in Figure 4(b). There was still an obvious bias in Figures 4(a) and 4(b) from the color bar, which was weakened a lot in Figure 4(c) based on the calibrated LAM data after orbit improvement.
Differences of DEMs compared to LOLA DEM. (a) model from original LAM data, (b) model from LAM data with orbit improvement, and (c) model from LAM data with calibration.
Original LAM data
LAM data with orbit improvement
LAM data with calibration
Table 3 exhibited the highest and lowest points of lunar topography obtained from different DEMs. Among them, the result of CE-1 (Li Chunlai) was obtained by Li Chunlai et al. using all of about 9.12 million altimeter data [12]. The result of LOLA was obtained by Smith using LOLA measurements over seven years up to 2015 [18]. The positions and elevations of the highest points in all models were almost the same, which were located on the uplifts adjacent to the Korolev and Dirichlet-Jackson basins on the far side of the Moon. The lowest point of DEM obtained using calibrated LAM data in this paper was quite different from the other two models. According to Li Chunlai’s analysis, the lowest point from LAM data was obtained at Beijing time 2008.07.15-13: 52: 29.778, but the LAM data used in this paper ended in 2008-07-07 for the limit of precise orbit. Therefore, the position and elevation of the lowest point in the DEM obtained by LAM-2B-SHAO-cali data were quite different from those in other DEMs. This result was similar to the CLTM-s01 model obtained by the first 3 months of CE-1 LAM data [11].
The highest and lowest points from different lunar DEMs.
Data source
Highest point of lunar surface
Lowest point of lunar surface
Max elevation difference /km
Longitude
Latitude
Elevation /km
Longitude
Latitude
Elevation /km
CE-1 (Li Chunlai)
158.656°W
5.441°N
10.629
172.413°W
70.368°S
-9.178
19.807
LAM-2B-SHAO -cali
158.594°W
5.469°N
10.723
148.625°W
61.375°S
-8.531
19.254
LOLA
158.63°W
5.341°N
10.792
172.48°W
70.36°S
-9.129
19.921
The grid model of topography could be expanded to a spherical harmonic form according to(7)Hλ,φ=∑l=0N∑m=0lP¯lmsinφC¯lmcosmλ+S¯lmsinmλC¯lm=ΔλΔθ4π∑i=1A∑j=1BHλi,φjP¯lmsinφjcosmλicosφiS¯lm=ΔλΔθ4π∑i=1A∑j=1BHλi,φjP¯lmsinφjsinmλicosφi
In formula (7), λ and φ referred to the longitude and latitude, respectively, P-lm referred to the normalized Legendre function, C-lm and S-lm referred to the normalized spherical harmonic coefficients, the summation index N referred to the max order, and A and B referred to the number of grid points in the latitudinal and longitudinal directions, respectively. We expanded the spherical harmonic coefficients to the 1400th order according to the above formula.
Based on the results of spherical harmonic expansion, the mean radius of the Moon was 1737.152 km, the polar radius was 1735.658 km, the equatorial radius was 1737.613 km, and the lunar flattening rate 1/f was 1/887.70. The lunar COF (Center of Figure) and COM (Center of Mass) had a deviation of 0.241 km from the north of the lunar rotation axis, and a deviation of -1.782 km and -0.728 km, respectively, in the x and y directions of the Moon-fixed coordinate system. The deviation was mainly along the direction of the Earth-Moon. It was estimated that the lunar crust on the far side was 10 km thicker than on the near side, which might be a major cause of the deviation of the lunar COM and COF towards the Earth’s direction [18].
4. DEM Accuracy Assessment
In this section, we compared the accuracy of different DEMs using the positions of known control points on the lunar surface. The spectral analysis was also adopted to evaluate the DEM accuracy.
4.1. Accuracy Assessment with LLRR and ALSEP
As shown in Figure 5, LLRR (Lunar Laser Ranging Reflector) was a series of laser reflectors placed by Lunakhod, Apollo 11, Apollo 14, and Apollo 15 on the near side of the Moon. ALSEP (Apollo Lunar Surface Experiments Package) was an experimental setup equipped with radio transponders placed during Apollo missions 12, 14, 15, 16, and 17. The coordinate accuracy of LLRR and ALSEP was less than 3 m and 30 m, respectively [19, 20]. We could interpolate the topography model to calculate the elevation of the coordinates of LLRR and ALSEP landmarks and compared them with the elevation values of LLRR and ALSEP themselves to measure the accuracy of the lunar DEM. The results were listed in Table 4.
LLRR and ALSEP elevation error for different DEMs.
Reflector
Original LAM
LAM-2B-SHAO
LAM-2B-SHAO _cali
LOLA (initial phase)
Lunakhod 1 (LRRR)
-146.14
-141.66
-3.27
-3.66
Lunakhod 2 (LRRR)
-11.40
-48.37
65.62
-17.19
Apollo 11 (LRRR)
-140.17
-99.99
87.75
-0.4
Apollo 14 (LRRR)
-145.59
-85.45
-1.04
6.3
Apollo 15 (LRRR)
-180.92
-185.89
-49.72
-20.18
Apollo 12 (ALSEP)
-158.49
-83.7
-14.52
-23.15
Apollo 14 (ALSEP)
-160.59
-100.45
-16.04
-8.7
Apollo 15 (ALSEP)
-180.92
-185.89
-49.72
-20.18
Apollo 16 (ALSEP)
-187.24
-184.35
-16.72
-36.29
Apollo 17 (ALSEP)
-66.29
-46.61
77.48
-11.81
RMS/m
147.70
127.31
48.70
17.98
Distribution of LLRR and ALSEP on the lunar surface.
In Table 4, the difference between original LAM and LAM-2B-SHAO was from the orbit improvement, which had an improvement from 147.7 m to 127.3 m. The accuracy improvement of the LAM-2B-SHAO-cali DEM when compared with LAM-2B-SHAO mainly came from the correction of ~139 m bias. The improvement affected elevation directly, and the elevation accuracy of LLRR and ALSEP sites was improved obviously from 127.3 m to 48.7 m.
For comparison, the result of LOLA was also showed in Table 4, and we could see that the accuracy of LOLA derived DEM is better than that of LAM data, not only for its higher precision of altimeter data but also for its higher space resolution.
4.2. Spectral Analysis
The variance σl of various orders of the spherical harmonic model could reflect the intensity of the spatial spectrum signal of the lunar topographic model:(8)σl=∑m=0lC¯lm2+S¯lm22l+1
It could be seen from Figure 6 that CE-1 DEM and LRO DEM were highly consistent at a lower order. After the 200th order, the power spectrum of the LOLA model was higher than the LAM models. This meant that the LOLA DEM had a higher spatial resolution at small scales and could provide more topography information.
Power spectrum of CE-1 DEM and LRO DEM.
5. Conclusion
In this paper, an empirical calibration method was proposed to calibrate the systematic errors in CE-1 LAM data. A new lunar DEM was derived using the calibrated LAM data. The DEM accuracy was evaluated by the comparison with LLRR and ALSEP reference coordinates. The results showed that the DEM accuracy was improved after the correction of the systematic errors in LAM data.
The result suggested that the constant bias of LAM was about -139.52 m. A lunar DEM grid model with the spatial resolution of 0.0625° × 0.0625° was obtained based on the calibrated LAM data. After calibration, the DEM accuracy was improved from about 127.3 m to about 48.7 m, compared with the coordinates of LLRR and ALSEP landmarks on the lunar surface.
Data Availability
The data used to support the findings of this study are available from the corresponding author upon request.
Conflicts of Interest
The authors declare that they have no conflicts of interest.
Acknowledgments
This work was supported by the National Natural Science Foundation of China (Grants nos. 11473056 and 11403076), the Science and Technology Commission of Shanghai (Grant no. 12DZ2273300), the Planetary Sciences Laboratory of Chinese Academy of Sciences, the Lunar Exploration Project of China, and the Key Laboratory of Space Object Measurement. The authors also thank C.K. Shum (The Ohio State University) and H.S. Fok (Wuhan University) for their support.
PottsN. J.GulliksonA. L.CurranN. M.DhaliwalJ. K.LeaderM. K.RegeR. N.KlausK. K.KringD. A.Robotic traverse and sample return strategies for a lunar farside mission to the Schrödinger basin2015554124112542-s2.0-8492159744310.1016/j.asr.2014.11.028ArakiH.TazawaS.NodaH.IshiharaY.GoossensS.SasakiS.KawanoN.KamiyaI.OtakeH.OberstJ.ShumC.Lunar global shape and polar topography derived from Kaguya-LALT laser altimetry200932359168979002-s2.0-6014909135010.1126/science.1164146OuyangZ. Y.Science Result of Chang’e-1 Lunar Orbiter and Mission Goals of Chang’e-2201019516YanJ.PingJ.LiF.Chang'E-1 precision orbit determination and lunar gravity field solution2010461505710.1016/j.asr.2010.03.0022-s2.0-79957876522YanJ.GoossensS.MatsumotoK.PingJ.HaradaY.IwataT.NamikiN.LiF.TangG.CaoJ.HanadaH.KawanoN.CEGM02: An improved lunar gravity model using Chang'E-1 orbital tracking data2012621192-s2.0-8485697493710.1016/j.pss.2011.11.010Nath GoswamiJ.Chandrayaan-1: India's first planetary science mission2008964VondrakR.KellerJ.ChinG.GarvinJ.Lunar reconnaissance orbiter (LRO): Observations for lunar exploration and science20101501-47222-s2.0-7795018881610.1007/s11214-010-9631-5DiK.HuW.LiuY.PengM.Co-registration of Chang'E-1 stereo images and laser altimeter data with crossover adjustment and image sensor model refinement20125012161516282-s2.0-8486766332710.1016/j.asr.2012.06.037ZhaoS.RanX.FuJ.GuoQ.Combined adjustment of CE-1 stereo camera image and laser altimeter data20144312122412292-s2.0-84920523015BarkerM. K.MazaricoE.NeumannG. A.ZuberM. T.HaruyamaJ.SmithD. E.A new lunar digital elevation model from the Lunar Orbiter Laser Altimeter and SELENE Terrain Camera20162733463552-s2.0-8493919563010.1016/j.icarus.2015.07.039Ping SJ.HuangQ.Yan GJ.Lunar topographic model CLTM-s01 from ChangE-1 laser altimeter2008381116011612Li LC.RenX.LiuJ.Laser Altimetry Data of Chang’E-1 and the Global Lunar DEM Model20103281293Cai CZ.LiangY.LiJ.Digital Elevation Model of the Moon from the ChangE-1 Laser Altimeter201025411531160HuW.DiK.LiuZ.PingJ.A new lunar global DEM derived from chang'E-1 laser altimeter data based on crossover adjustment with local topographic constraint2013871731822-s2.0-8488837386810.1016/j.pss.2013.08.004SmithD. E.ZuberM. T.NeumannG. A.Initial observations from the Lunar Orbiter Laser Altimeter201037, article L18204PingJ.SuX.LiuJ.ShuR.YanJ.External calibration for laser altimetry mission in Chang’E-1201343111438144710.1360/132013-309ChangS.HuangY.LiP.HuX.FanM.The use of laser altimetry data in Chang'E-1 precision orbit determination2016169839010.1088/1674-4527/16/9/142Smith ED.Zuber TM.Neumann AG.Summary of the results from the lunar orbiter laser altimeter after seven years in lunar orbit20162837091DaviesM. E.ColvinT. R.Lunar coordinates in the regions of the Apollo landers2000105820277202802-s2.0-003376580410.1029/1999JE001165FokH. S.ShumC. K.YiY.ArakiH.PingJ.WilliamsJ. G.FotopoulosG.NodaH.GoossensS.HuangQ.IshiharaY.MatsumotoK.OberstJ.SasakiS.Accuracy assessment of lunar topography models201163115232-s2.0-8492256344010.5047/eps.2010.08.005