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.     

On the information contents and regularisation of lunar gravity field solutions

Till the present day the recovery of the lunar gravity field from satellite tracking data depends in a crucial way on the level and method of regularisation. With Earth-based tracking only, the spatial data coverage is limited to only slightly more than 50% and the inverse problem remains severely ill-posed. The development of global gravity models suitable for precise orbit modeling as well as geophysical studies therefore requires a significant level of regularisation, limiting the solution power over the far-side where no gravity information is available. Unconstrained solutions, within the framework of global harmonic base functions, are only possible for very low degrees (< 10). Any significant change to this situation is only to be expected when global satellite-to-satellite tracking data of high quality becomes available early in the next decade. Yet, a rigorous analysis of the impact of the chosen method and level of regularisation is lacking. Most gravity models employ a Kaula-type signal smoothness constraint of 15 × 10⁻⁵ /l², which allows a good overall data fit as well as a smooth field over the far-side. Furthermore, a geographical type of constraint has been suggested, where surface accelerations have been introduced in areas of no data coverage. Modern numerical methods, on the other hand, offer direct tools and search mechanisms for the optimal level of regularisation. This paper presents a study of Tikhonov-type regularisation of lunar gravity solutions, with emphasis on the so-called L-curve and quasi-optimality methods for regularisation parameter estimation. Furthermore, new quality measures of lunar gravity solutions are presented, which account for the bias introduced by the regularisation.     

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.     

Chang’E-1 precision orbit determination and lunar gravity field solution

In this paper we present results assessing the role of Very Long Baseline Interferometry (VLBI) tracking data through precision orbit determination (POD) during the check-out phase for Chang’E-1, and the lunar gravity field solution CEGM-01 based on the orbital tracking data acquired during the nominal phase of the mission. The POD of Chang’E-1 is performed using S-band two-way Range and Range Rate (R&RR) data, together with VLBI delay and delay rate observations. The role of the VLBI data in the POD of Chang’E-1 is analyzed, and the resulting orbital accuracies are estimated for different solution strategies. The final orbital accuracies proved that the VLBI tracking data can improve the Chang’E-1 POD significantly. Consequently, CEGM-01 based on six-month tracking data during Chang’E-1 nominal mission phase is presented, and the accuracy of the model is assessed by means of the gravity field power spectrum, admittance and coherence between gravity and topography, lunar surface gravity anomaly and POD for both Chang’E-1 and Lunar Prospector (LP). Our analysis indicates that CEGM-01 has significant improvements over a prior model (i.e. GLGM-2), and shows the potential of Chang’E-1 tracking data in high resolution lunar gravity field model solution by combining with SELENE and LP tracking data.     

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.     

On the usefulness of atmospheric and oceanic angular momentum in recovering polar motion and gravity field variations in a unified process

Earth rotation parameters (ERPs) are excited by variations in the mass distribution on the Earth’s surface and the exchange of angular momentum between the atmosphere and oceans and the solid Earth. The same mass redistribution causes temporal changes in the gravity field coefficients with the second degree harmonics related to the rotational deformation and hence to changes in the Earth’s inertial tensor. If precise models of the atmospheric and oceanic angular momentum (AM) are available solution for polar motion and degree 2 Stokes harmonics can be unified. In this study we utilize SLR tracking of LAGEOS to compare (i) degree 2 harmonics from ERPs and gravitation, and (ii) LAGEOS excitation functions and geophysical data (mass + motion). Finally, we investigate to what extent a unified approach is possible with current models for AM data and gravity mass change estimated from ERP within orbit determinations.     

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.     

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.     

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.     

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.     

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.     

Autonomous navigation for lunar final approach based on gravity gradient measurements

Lunar final approach navigation is critical for pin-point lunar landing in future missions. This study investigates the use of lunar gravity gradient measurements for autonomous navigation of a lunar probe during the final approach phase. As the spacecraft approaches the Moon, the strength of gravity gradient signals improves. A spaceborne gravity gradiometer can precisely measure local gravity gradients, and the latest lunar gravity model GL1500E is used to provide reference values. The employed truncation degree and order of the gravity model are increased stepwise considering the decreasing altitude of the spacecraft in order to reach a compromise between computational costs and model accuracy. An iterative Kalman filter is developed for coupled orbit and attitude estimation using gravity gradient measurements and attitude quaternions obtained from star sensors. A simulated spacecraft with a gradiometer noise level of 0.01 E is considered. Simulation results show that the spacecraft’s position converges rapidly and achieves an accuracy of less than 100 m at the last epoch.     

Impact of far-side satellite tracking on gravity estimation in the SELENE project

Preliminary results of numerical simulation are presented to examine the gravity estimation capability in the Japanese lunar exploration project SELENE (SELenological and ENgineering Explorer), which will be launched in 2003. One of the new characteristics of the SELENE lunar gravimetry is 4-way satellite-to-satellite Doppler tracking of a low-altitude lunar orbiter by means of a high-altitude relay satellite. It is shown that planned satellites configuration will provide a good far-side data coverage of the lunar orbiter and will improve lunar gravity field as well as far-side selenoid.     

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.     

模拟月壤铺粉过程DEM数值仿真

空基激光选区熔化（SLM）技术与原位资源利用（ISRU）概念结合，有望解决地外大规模基地建设的工程难题。SLM铺粉过程对成形件性能和质量有重要影响。基于非球形粒子叠加球模型方法，建立模拟月壤颗粒几何模型；基于线性弹簧-阻尼接触作用模型、Hamaker理论及牛顿运动定律，建立颗粒动力学模型；采用三维离散单元方法（DEM）及软球模型，进行不同工况下模拟月壤在铺粉过程中的流变特性研究。结果显示:所提模型和方法能开展指定工况和环境参数的模拟月壤颗粒系统流动性和堆积行为数值仿真研究；月面低重力环境导致粉床表面粗糙度变大、堆积密度和平均配位数变小；通过降低铺粉速度和优化刮刀型面，可以有效改善月基铺粉的粉床质量，获得更密实和均匀的粉床结构。      

Lunar laser altimetry in the SELENE project

Lunar laser altimeter (LALT) has been developed for the Japanese lunar exploration SELENE (SELenological and ENgineering Explorer) in 2003. Lunar laser altimetry and two other selenodetic missions are put together and called RISE (Researches In SElenodesy) project that has been developed mainly by National Astronomical Observatory of Japan. Main objective of this measurements is to construct a more global and accurate topographic model of the moon than Clementine's GLTM-2. The newly developed topography of the moon will make a great contribution to many problems of lunar geology, geophysics, and the reduction of lunar occultation data. In order to expand the coverage in latitude, LALT has a function of slant ranging in the direction about 40 degrees aside from nadir. LALT is designed to have a capability of detecting most of returned pulses from the lunar surface even for the slant ranging with proper threshold level.     

An alternative computation of a gravity field model from GOCE

GOCE is the first satellite with a gravitational gradiometer (SGG). This allows to determine a gravity field model with high spatial resolution and high accuracy. Four of the six independent components of the gravitational gradient tensors (GGT) are measured with high accuracy in the so-called measurement band (MB) from 5 to 100 mHz by the GOCE gradiometer. Based on more than 1 year of GOCE measurements, two gravity field models have been derived. Here, we introduce a strategy for spherical harmonic analysis (SHA) from GOCE measurements, with a bandpass filter applied to the SGG data, combined with orbit analysis based on the integral equation approach, and additional constraints (or stabilization) in the polar areas where no observation is available due to the orbit geometry. In addition, we combined the GOCE SGG part with a set of GRACE normal equations. This improves the accuracy of the gravity field in the long-wavelength parts, due to the complementarity of GOCE and GRACE. Comparison with other models and with external data shows that our results are rather close to the GPS-levelling data in well-selected test regions, with an uncertainty of 4–7 cm, for truncation at degree 200.     

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.     

Analysis of data of “Clementine” and “KAGUYA” missions and “ULCN” and “KSC-1162” catalogues

In this paper an analysis of data coordinate systems from selenographic catalogues and space missions was carried out. The lunar macrorelief models were made on basis of the software package ASNI USTU using method of the spherical harmonic expansion. These models accurately describe the global features of the lunar figure. To construct these models the following sources of topographic information were used: “Clementine” and “KАGUYА” (Selena, Japan mission) missions, “KSC-1162” (Kazan selenocentric catalogue), “Kiev” (selenodesic catalogue), “SAI” (Chuikova (1975)), “Bills, Ferrari”, “ULCN” (The Unified Lunar Control Network 2005). Direct comparison hypsometric information “KSС-1162” catalogue data with “Clementine” mission was carried out. These researches confirmed a good agreement of the hypsometric information of compared systems. The normalized coefficients were obtained on basis of the hypsometric information expansion for eight sources. The displacement of the lunar center of mass (LCM) relatively to the lunar center of figure (LCF) was obtained by using topographic data selenodetical catalogues and space missions.     

Balloon-borne,high-altitude gravimetry

Gravity measurements from a high-altitude balloon can verify global and upward-continued gravity models. A gravimeter suspended beneath a balloon is in a dynamic, and largely unpredictable, environment sensing accelerations due to gravity and balloon motions. Independent measurements of balloon motions using inertial navigation data combined with ground tracking data will allow for separation of balloon-induced accelerations from gravitational accelerations. Analysis of these data must estimate: 1) vertical gravimeter accelerations due to motion and gravity, 2) horizontal velocity to estimate the Eötvös effect, and 3) gravimeter position for comparison with gravity models. The first engineering test flight occurred on 11 October 1983, during the seasonal wind reversal and was very successful. Flight duration was approximately seven hours, with two hours of data collected at each of 30 km and 26 km altitudes. The results include gravity estimates, design criteria for future flights and feasibility analysis for vertical gravity profiles during ascent and descent.