The simulation of lunar gravity field recovery from D-VLBI of Chang’E-1 and SELENE lunar orbiters

The lunar gravity field is a foundation to study the lunar interior structure, and to recover the evolution history of the Moon. It is still an open and key topic for lunar science. For above mentioned reasons, it becomes one of the important scientific objectives of recent lunar missions, such as KAGUYA (SELENE) the Japanese lunar mission and Chang’E-1, the Chinese lunar mission. The Chang’E-1 and the SELENE were successfully launched in 2007. It is estimated that these two missions can fly around the Moon longer than 6 months simultaneously. In these two missions, the Chinese new VLBI (Very Long Baseline Interferometry) network will be applied for precise orbit determination (POD) by using a differential VLBI (D-VLBI) method during the mission period. The same-beam D-VLBI technique will contribute to recover the lunar gravity field together with other conventional observables, i.e. R&RR (Range and Range Rate) and multi-way Doppler. Taking VLBI tracking conditions into consideration and using the GEODYNII/SOVLE software of GSFC/NASA/USA [8 and 10], we simulated the lunar gravity field recovering ability with and without D-VLBI between the Chang’E-1 and SELENE main satellite. The cases of overlapped flying and tracking period of 30 days, 60 days and 90 days have been analyzed, respectively. The results show that D-VLBI tracking between two lunar satellites can improve the gravity field recovery remarkably. The results and methods introduced in this paper will benefit the actual missions.     

Proposal of application of same beam VLBI measurements in precision orbit determination of lunar orbiter and return capsule and lunar gravity field simulation in Chang’E-3 mission

High accuracy differenced phase delay can be obtained by observing multiple point frequencies of two spacecraft using the same beam Very Long Baseline Interferometry (VLBI) technology. Its contribution in lunar spacecraft precision orbit determination has been performed during the Japanese lunar exploration mission SELENE. In consideration that there will be an orbiter and a return capsule flying around the moon during the Chinese lunar exploration future mission Chang’E-3, the contributions of the same beam VLBI in spacecraft precision orbit determination and lunar gravity field solution have been investigated. Our results show that the accuracy of precision orbit determination can be improved more than one order of magnitude after including the same beam VLBI measurements. There are significant improvements in accuracy of low and medium degree coefficients of lunar gravity field model obtained from combination of two way range and Doppler and the same beam VLBI measurements than the one that only uses two way range and Doppler data, and the accuracy of precision orbit determination can reach meter level.     

Comparison analyses on the 150 × 150 lunar gravity field models by gravity/topography admittance,correlation and precision orbit determination

We analyzed the 150 × 150 lunar gravity field models, LP150Q, GLGM-3 and SGM150, using the power spectrum on the lunar nearside and farside, the lunar global and localized gravity/topography admittance and correlation, and Chang’E-2 precision orbit determination to investigate which model is a more effective tool to estimate geophysical parameters and determine the lunar satellite precision orbit. Results indicate that all gravity field models can be used to estimate the lunar geophysical parameters of the nearside of the Moon. However, SGM150 is better in such computation of the farside. Additionally, SGM150 is shown to be the most useful model for determining the lunar satellite orbit.     

Simulation study for the determination of the lunar gravity field from PRARE-L tracking onboard the German LEO mission

A simulation study has been performed at GFZ Potsdam, which shows the anticipated improvement of the lunar gravity field model with respect to current (LP150Q model) or near-future (SELENE) knowledge in the framework of the planned German Lunar Explorations Orbiter (LEO) mission, based on PRARE-L (Precise Range And Range-rate Equipment – Lunar version) Satellite-to-Satellite (SST) and Satellite-Earth-Satellite (SEST) tracking observations. It is shown that the global mean error of the lunar gravity field can be reduced to less than 0.1 mGal at a spatial resolution of 50 km. In the spectral domain, this means a factor of 10 (long wavelengths) and some 100 (mid to short wavelengths) improvement as compared to predictions for SELENE or a factor of 1000 with respect to LP150Q. Furthermore, a higher spatial resolution of up to 28 km seems feasible and would correspond to a factor of 2–3 improvement of SELENE results. Moreover, PRARE-L is expected to derive the low-degree coefficients of the lunar gravity field with unprecedented accuracy. Considering long mission duration (at least 1 year is planned) this would allow for the first time a precise direct determination of the low-degree tidal Love numbers of the Moon and, in combination with high precision SEST, would provide an experimental basis to study relativistic effects such as the periselenium advance in the Earth–Moon system.     

New results and questions of lunar exploration from SELENE,Chang’E-1, Chandrayaan-1 and LRO/LCROSS

The moon has longstanding questions such as lunar environments, origin, formation and evolution, magnetization of crustal rocks, internal structure and possible life. The recent lunar missions, e.g., SELenological and ENgineering Explorer “KAGUYA” (SELENE), Chang’E-1, Chandrayaan-1, and Lunar Reconnaissance Orbiter/Lunar CRater Observation and Sensing Satellite (LRO/LCROSS), have provided new opportunities to explore and understand these issues. In this paper, we reviewed and presented the results and findings in the fields of lunar gravity, magnetic field, atmosphere, surface geomorphology and compositional variations, volcano, craters, internal structure, water and life science from new lunar exploration missions. In addition, the new objectives and scientific questions on lunar explorations in near future are presented and discussed.     

Selenodetic experiments of SELENE: Relay subsatellite,differential VLBI,and laser altimeter

Since 1960s, the gravitational potential of the Moon has been extensively studied from Doppler tracking data between a ground station and spacecraft orbiting in front of the Moon (e. g., Lorell and Sjogren, 1968; Bills and Ferrari, 1980; Konopliv et al., 1993; Lemoine et al., 1997). Because direct radio communication is interrupted while spacecraft is orbiting behind the Moon, however, the coverage of tracking data has been limited mostly to the nearside of the Moon so far. In order to compensate for such lack of tracking data, we propose satellite-to-satellite Doppler measurement by using a relay subsatellite in Japanese mission to the Moon in 2003. A complete coverage of Doppler tracking from an orbiter at sufficiently low altitude will significantly improve lunar gravity model and will contribute for future geophysical study of interior and tectonics on the Moon. Further, we propose differential VLBI experiment between the subsatellite and a propulsion module landed on the surface of the Moon. The differential VLBI is about 10 times more accurate than conventional Doppler measurement for long-wavelength gravity field. Besides, differential VLBI is sensitive to the displacement perpendicular to the line of sight. Thus the VLBI experiment provides precise estimates of the lunar gravity potential at low degree. The last proposal for selenodetic experiments is a laser altimeter. Global topography model has been already developed from the analysis of Clementine LIDAR data (Zuber et al., 1994), but it is suggested that the model includes appreciable anisotropy between NS and E-W directions due to highly eccentric orbit of Clementine spacecraft (Bills and Lemoine, 1995). The laser altimeter experiment from an orbiter in nearly circular orbit will provide a new reference for the isotropic lunar topography model.     

A simulation study for anticipated accuracy of lunar gravity field model by SELENE tracking data

Results of numerical simulations are presented to examine the global gravity field recovery capability of the Japanese lunar exploration project SELENE (SELenological and ENgineering Explorer) which will be launched in 2007. New characteristics of the SELENE lunar gravimetry include 4-way satellite-to-satellite Doppler tracking of main orbiter and differential VLBI tracking of two small free-flier satellites. It is shown that the proposed satellite constellation will provide the first truly global satellite tracking data coverage. The expected results from these data are; (1) drastic reduction in far-side gravity error, (2) estimation of many gravity coefficients by the observation, not by a priori information, and (3) one order of magnitude improvement over existing gravity models for low-degree field.     

Gamma-ray detector on board lunar mission Chang’E-1

Gamma-ray spectrometer (GRS) is included in the payload of Chinese first lunar mission Chang’E-1 that will be launched in 2007. Specific objectives of the GRS are to map abundance of O, Si, Fe, Ti, U, Th, K, and perhaps, Mg, Al, and Ca to depths of about 20 cm. There are remarkable advantages for GRS application to remote sensing elemental materials over the entire lunar surface: large effective area and good ability for background rejection. We will describe the design of GRS and present its performance in this paper. Moreover, the GRS calibration will also be introduced.     

Three-dimensional tracking of a lunar satellite with differential Very-Long-Baseline-Interferometry

Lunar gravimetry mission in the Japanese lunar exploration project SELENE (Selenological and Engineering Explorer) is characterized by inter-satellite tracking by means of a relay satellite in a high eccentric orbit, combined with differential Very-Long-Baseline-Interferometry (ΔVLBI) and conventional 2-way Doppler tracking. ΔVLBI provides information on the satellite position and velocity complementary to conventional range and range rate measurement, and allows us to measure lunar gravitational accelerations in all the three components. In this article, ΔVLBI and 2-way Doppler numerical simulation results are compared to those obtained from 2-way Doppler observations only, so that we can evaluate the contribution of ΔVLBI to the SELENE lunar gravimetry mission.     

GPS-only gravity field recovery with GOCE,CHAMP, and GRACE

Gravity missions such as the Gravity field and steady-state Ocean Circulation Explorer (GOCE) are equipped with onboard Global Positioning System (GPS) receivers for precise orbit determination (POD), instrument time-tagging, and the extraction of the long wavelength part of the Earth’s gravity field. The very low orbital altitude of the GOCE satellite and the availability of dense 1 s GPS tracking data are ideal characteristics to exploit the contribution of GPS high-low Satellite-to-Satellite Tracking (hl-SST) to gravity field determination. We present gravity field solutions based on about 8 months of GOCE GPS hl-SST data from 2009 and compare the results with those obtained from the CHAllenging Minisatellite Payload (CHAMP) and Gravity Recovery And Climate Experiment (GRACE) missions. The very low orbital altitude of GOCE significantly improves gravity field recovery from GPS hl-SST data above degree 20, but not for the degrees below 20, where the quality of the spherical harmonic coefficients remains essentially unchanged. Despite the limited time span of GOCE data used, the gravity field of the Earth can be resolved up to about degree 115 using GPS data only. Empirically determined phase center variations (PCVs) of the GOCE onboard GPS helix antenna are, however, mandatory to achieve this performance.     

Lunar international science coordination/calibration targets (L-ISCT)

Eight lunar areas, each ∼200 km in diameter, are identified as targets for coordinated science and instrument calibration for the orbital missions soon to be flown. Instrument teams from SELENE, Chang’E, Chandrayaan-1, and LRO are encouraged to participate in a coordinated activity of early-release data that will improve calibration and validation of data across independent and diverse instruments. The targets are representative of important lunar terrains and geologic processes and thus will also provide a broad introduction to lunar science for new investigators. We briefly identify additional cross-calibration issues for instruments that produce time series data rather than maps.     

Design and analysis of an extended mission of CE-2: From lunar orbit to Sun–Earth L2 region

Chang’E-2 (CE-2) has firstly successfully achieved the exploring mission from lunar orbit to Sun–Earth L2 region. In this paper, we discuss the design problem of transfer trajectory and at the same time analyze the visible segment of Tracking, Telemetry & Control (TT&C) system for this mission. Firstly, the four-body problem of Sun–Earth–Moon and Spacecraft can be decoupled in two different three-body problems (Sun–Earth + Moon Restricted Three-Body Problems (RTBPs) and Earth–Moon ephemeris model). Then, the transfer trajectory segments in different model are computed, respectively, and patched by Poincaré sections. The full-flight trajectory including transfer trajectory from lunar orbit to Sun–Earth L2 region and target Lissajous orbit is obtained by the differential correction method. Finally, the visibility of TT&C system at the key time is analyzed. Actual execution of CE-2 extended mission shows that the trajectory design of CE-2 mission is feasible.     

藏匿在冯·卡门坑南部下方的月球质量瘤

月球着陆探测的候选着陆区域的质量异常分布,预示着对应区域发生过复杂的动力学演化过程,存在着异常丰富的物质组成。对月球的质量异常分布或重力场异常的探测是基于绕月卫星重力探测方法获取的,结合“嫦娥1号”月球探测获得的月球重力场异常模型与地形探测模型的数据,分析揭示存在于月球背面的质量瘤区域,研究结果表明：使用多次任务的数据首次证认了一个位于中尺度撞击区域,冯·卡门（Von Kármán）撞击坑的南半部和边缘的正下部的质量瘤,暂定编号CEFC04,冯·卡门撞击坑是该区域形成后,又一次大型撞击产生的。这类显著的被藏匿的中尺度月球质量瘤的发现,对探测研究月球的演化有重要意义。     

Influence of lunar topography on simulated surface temperature

The surface temperature of the Moon is one of the essential parameters for the lunar exploration, especially to evaluate the Moon thermophysical features. The distribution of the temperature is heavily influenced by the Moon topography, which, however, is rarely studied in the state-of-art surface temperature models. Therefore, this paper takes the Moon topography into account to improve the surface temperature model, Racca model. The main parameters, such as slopes along the longitude and latitude directions, are estimated with the topography data from Chang’E-1 satellite and the Horn algorithm. Then the effective solar illumination model is then constructed with the slopes and the relative position to the subsolar point. Finally, the temperature distribution over the Moon surface is obtained with the effective illumination model and the improved Racca model. The results indicate that the distribution of the temperature is very sensitive to the fluctuation of the Moon surface. The change of the surface temperature is up to 150 K in some places compared to the result without considering the topography. In addition, the variation of the surface temperature increases with the distance from the subsolar point and the elevation, along both latitude and longitude directions. Furthermore, the simulated surface temperature coincides well with the brightness temperature in 37 GHz observed by the microwave sounder onboard Chang’E-2 satellite. The corresponded emissivity map not only eliminates the influence of the topography, but also hints the inherent properties of the lunar regolith just below the surface. Last but not the least, the distribution of the permanently shadowed regions (PSRs) in the lunar pole area is also evaluated with the simulated surface temperature result.     

Lunar gravity field modeling critical analysis and challenges

This paper summarizes and provides a critical analysis of the historical developments of lunar gravitational models from the earliest use of ground based tracking systems of the Lunar Orbiter to the Lunar Prospector mission. This encompasses a comprehensive and critical analysis of the various methods used in the estimation of the gravity coefficients and the processing of large batches of diverse measurements and data types. It has been shown that weakness exists in the current models of the lunar gravity field, which is primarily due to the lack of far side lunar tracking data information, which makes the lunar potential modeling difficult but expected to be overcome as data from SELENE satellite-to-satellite tracking becomes available. Comparisons of various lunar models reveal an agreement in the low order coefficients of the spherical harmonics. However, substantial differences in the models exist in the higher-order harmonics. A numerical comparison has been presented showing the performance of all the contemporary lunar gravitational models used within the astrodynamics community and available in public domain. Improvements to the current models are part of a continuing process and the recent model improvements and future possibilities in lunar gravity modeling are discussed. A brief review of the recent missions has been presented. It is hoped that this critical review will benefit the researchers by presenting the historical as well as state of the art in this field.     

Long-wavelength lunar gravity field recovery from simulated orbit and inter-satellite tracking data

As has been demonstrated recently, inter-satellite Ka-band tracking data collected by the GRAIL (Gravity Recovery And Interior Laboratory) spacecraft have the potential to improve the resolution and accuracy of the lunar gravity field by several orders of magnitude compared to previous models. By means of a series of simulation studies, here we investigate the contribution of inter-satellite ranging for the recovery of the Moon’s gravitational features; the evaluation of results is made against findings from ground-based Doppler tracking. For this purpose we make use of classical dynamic orbit determination, supported by the analysis of satellite-to-satellite tracking observations. This study sheds particularly light on the influence of the angular distance between the two satellites, solar radiation modeling and the co-estimation of the lunar Love number k₂. The quality of the obtained results is assessed by gravity field power spectra, gravity anomalies and precision orbit determination. We expect our simulation results to be supportive for the processing of real GRAIL data.     

The Solar Wind interactions with Lunar Magnetic Anomalies: A case study of the Chang’E-2 plasma data near the Serenitatis antipode

In this paper, we present the first and preliminary results on the near-Moon plasma environment, based on the spectrogram data obtained with the Solar Wind Ion Detector (SWID) onboard Chang’E-2 from four lunar orbits on 10–11 Oct 2010. These orbits, at a constant altitude of ∼100 km, approach gradually the Lunar Magnetic Anomaly (LMA) at the Serenitatis antipode. The data reveal tentatively a region with decrement in proton density and enhancement in temperature. The near coincidence of this region with the Serenitatis antipode probably suggests the presence of a minimagnetosphere associated with the LMA, which effectively shields and heats the incident Solar Wind (SW) protons.     

VLBI for better gravimetry in SELENE

The Japanese lunar explorer SELENE (SElenological and Engineering Explorer), to be launched in 2007, will for the first time utilize VLBI observations in lunar gravimetry investigations. This will particularly improve the accuracy to which the low degree gravitational harmonics and the gravity field near the limb can be measured, and when combined with Doppler measurements will enable three-dimensional information to be extracted. Differential VLBI Radio sources called VRAD experiment involves two on-board sub-satellites, Rstar and Vstar. These will be observed using differential VLBI to measure the trajectories of the satellites with the Japanese network named VERA (VLBI Exploration of Radio Astrometry) and an international VLBI network.     

Testing Modified Newtonian Dynamics with LISA Pathfinder

We suggest that LISA Pathfinder, a technology demonstrator for the future gravitational wave observatory LISA, could be used to carry out a direct experimental test of Modified Newtonian Dynamics (MOND). The LISA Pathfinder spacecraft is currently being built and the launch date is just a few years away. No modifications of the spacecraft are required, nor any interference with its nominal mission. The basic concept is to fly LISA Pathfinder through the region around the Sun-Earth saddle point, in an extended mission phase, once the original mission goals are achieved. We examine various strategies to reach the saddle point, and find that the preferred strategy, yielding relatively short transfer times of just over 1 year, probably involves a lunar fly-by. LISA Pathfinder will be able to probe the intermediate MOND regime, i.e. the transition between deep MOND and Newtonian gravity. We present robust estimates of the anomalous gravity gradients that LISA Pathfinder should be exposed to, based on MONDian effects as derived from the Tensor-Vector-Scalar (TeVeS) theory. The spacecraft speed and spatial scale of the MOND signal combine in a way that the spectral signature of the signal falls precisely into LISA Pathfinder’s measurement bandwidth. We find that if the gravity gradiometer on-board the spacecraft achieves its currently predicted sensitivity, these anomalous gradients could not just be detected, but measured in some detail.     

从日地系统L2出发借力月球飞越近地小行星

对于停留在日地系统L2的“嫦娥2号”探测器,其后续飞行方案有多个选项,例如主动撞月或重返月球轨道、返回地球轨道或再入大气、飞往地月系统L1/L2或日地系统L1、进入深空飞越近地小行星(最终,“嫦娥2号”于2012年12月13日成功地实现了对Toutatis小行星的近距离飞越)。探讨上述的飞行方案需要对飞行轨道进行初步设计,总的速度脉冲限制在100 m/s以内并且需要考虑探测器同时受到太阳、地球、月球的引力作用。本研究设计了探测器从日地系统L2出发借力月球实现Toutatis小行星飞越的飞行方案,与直接飞越方案相比,借力月球可以进一步节省探测器的燃料消耗,其等效速度脉冲设计值为58.47 m/s。