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.     

Modeling and analysis of Earth’s gravity field measurement performance by inner-formation flying system

The Earth’s gravity field can be measured with high precision by constructing the purely gravitational orbit of the inner-satellite in Inner-formation Flying System (IFS), which is independently proposed by Chinese scholars and offers a new way to carry out gravity field measurement by satellite without accelerometers. In IFS, for the purpose of quickly evaluating the highest degree of recovered gravity field model and geoid error as well as analyzing the influence of system parameters on gravity field measurement, an analytical formula was established by spectral analysis method. The formula can reflect the analytical relationship between gravity field measurement performance and system parameters such as orbit altitude, the inner-satellite orbit determination error, the inner-satellite residual disturbances, data sampling interval and total measurement time. This analytical formula was then corrected by four factors introduced from numerical simulation of IFS gravity field measurement. By comparing computation results from corrected analytical formula and the actual gravity field measurement performance by CHAMP, the correctness and rationality of this analytical formula were verified. Based on this analytical formula, the influences of system parameters on IFS gravity field measurement were analyzed. It is known that gravity field measurement performance is a monotone decreasing function of orbit altitude, the inner-satellite orbit determination error, the inner-satellite residual disturbances, data sampling interval and the reciprocal of total measurement time. There is a match relationship between the inner-satellite orbit determination error and residual disturbances, in other words, the change rate of gravity field measurement performance with one of them is seriously restricted by their relative size. The analytical formula can be used to quantitatively evaluate gravity field measurement performance fast and design IFS parameters optimally. It is noted that the analytical formula and corresponding conclusions are applied to any gravity satellite which measures gravity field by satellite perturbation orbit.     

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.     

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.     

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.     

“嫦娥4号”中继星任务分析与系统设计

作为"嫦娥4号"任务的重要组成部分,中继星将为着陆器和巡视器提供中继通信支持。不同于其它月球探测器,中继星首次选择了绕地月L2平动点运行的晕(Halo)轨道以保证对月球背面的着陆器和巡视器提供连续的中继通信服务,面临诸多技术挑战。在对中继星任务特点进行分析的基础上,梳理了研制中的技术难题,包括使命轨道的选择、使命轨道的到达和长期维持、中继通信体制选择等,并提出了解决方案。中继星的总体设计方案概述也在文中给出。     

An empirical model of electron density in low latitude at 600 km obtained by Hinotori satellite

An empirical model of electron density (Ne) was constructed by using the data obtained with an impedance probe on board Japanese Hinotori satellite. The satellite was in circular orbit of the height of 600 km with the inclination of 31 degrees from February 1981 to June 1982. The constructed model gives Ne at any local time with the time resolution of 90 min and between −25 and 25 degrees in magnetic latitude with its resolution of 5 degrees in the range of F_10.7 from 150 to 250 under the condition of K_p < 4. Spline interpolations are applied to the functions of day of year, geomagnetic latitude and solar local time, and linear interpolation is applied to the function of F_10.7. Longitude dependence of Ne is not taken into account. Our density model can reproduce solar local time variation of electron density at 600 km altitude better than current International Reference Ionosphere (IRI2001) model which overestimates Ne in night time and underestimates Ne in day time. Our density model together with electron temperature model which has been constructed before will enable more understanding of upper ionospheric phenomenon in the equatorial region.     

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.     

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 present status of the Japanese Penetrator Mission: LUNAR-A

The scientific objective of the LUNAR-A Japanese Penetrator Mission is to explore the lunar interior by seismic and heat-flow experiments. Two penetrators containing two-component seismometer and heat-flow probes will be deployed from a spacecraft onto the lunar surface, one on the nearside and the other on the farside of the moon. The data obtained by the penetrators will be transmitted to the ground station by way of the LUNAR-A mother spacecraft orbiting at an altitude of about 200 km. The seismic observations are expected to provide key data on the size of the lunar core, as well as data on the deep mantle structure. The heat-flow measurements at two different sites will also provide important data on the thermal structure and bulk concentrations of heat-generating elements in the Moon. These data will provide much stronger geophysical constraints on the origin and evolution of the Moon than has ever been obtained. The LUNAR-A mission was supposed to be launched in 2004. However, a malfunction of spacecraft subsystem and technical issues for penetrator system occurred during the course of the qualification level test. Therefore, further improvements and some modifications were considered to be required for reliability and robustness. The development of the mother spacecraft was temporarily suspended, while we have put a three-year program into effect to solve the penetrator technology issues.     

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.     

“嫦娥4号”月球背面软着陆任务设计

介绍了"嫦娥4号"月球背面软着陆任务设计方案。着陆区初步选定为月球背面南极–艾特肯(South PoleAitken,SPA)盆地内的冯·卡门(Von Kármán)撞击坑内。采用中继星实现着陆器和巡视器的对地通信,并选择环绕地月拉格朗日L2点的halo轨道作为其使命轨道。采用CZ-4C火箭和CZ-3B火箭,分别完成中继星和着陆器–巡视器组合体的发射。两器一星上共配置了6台国内研制科学载荷和3台国际合作科学载荷,开展以低频射电天文观测、巡视区形貌、矿物组份及浅层结构为主的科学探测。此外,还搭载了2颗月球轨道编队飞行微卫星、月面微型生态圈和大孔径激光角反射镜,分别开展超长波天文干涉测量试验、月面生态系统试验和超过地月距离的激光测距试验。通过创新设计顶层任务,充分继承成熟技术和产品,增加中继通信功能模块,开放资源引入高性能载荷和搭载项目,将实现一次低成本、短周期、大开放、高效益的月球探测任务。     

Spin axis orientation of Ajisai determined from Graz 2 kHz SLR data

The Graz 2 kHz Satellite Laser Ranging (SLR) measurements allow determination of the spin axis orientation of the geodetic satellite Ajisai. The high repetition rate of the laser makes it possible to determine the epoch time when the laser is pointing directly between two corner cube reflector (CCR) rings of the satellite. Identification of many such events during a few (up to 3) consecutive passes allows to state the satellite orientation in the celestial coordinate system. Six years of 2 kHz SLR data (October 2003–October 2009) delivered 331 orientation values which clearly show precession of the axis along a cone centered at 14^h56^m2.8^s in right ascension and 88.512° in declination (J2000.0 celestial reference frame) and with an half-aperture angle θ of 1.405°. The spin axis precesses with a period of 117 days, which is equal to the period of the right ascension of the ascending node of Ajisai’s orbit. We present a model of the axis precession which allows prediction of the satellite orientation – necessary for the envisaged laser time transfer via Ajisai mirrors.     

月球的构造格架及其演化差异

根据以GRAIL重力场数据和LOLA地形数据计算的月壳厚度,将月球构造格架初步划分为三个构造单元:主要覆盖月球正面风暴洋区域的月海构造单元、主要覆盖月球背面高地的月陆构造单元以及主要位于南半球背面的南极艾肯盆地构造单元。结合最新的研究成果,对各构造单元上的重大地质事件包括岩浆事件、火山事件和撞击事件分别进行了简单的梳理,结果表明月球不同构造单元的演化事件具有明显的差异。     

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.     

GOCE orbit analysis: Long-wavelength gravity field determination using the acceleration approach

The restricted sensitivity of the Gravity field and steady-state Ocean Circulation Explorer (GOCE) gradiometer instrument requires satellite gravity gradiometry to be supplemented by orbit analysis in order to resolve long-wavelength features of the geopotential. For the hitherto published releases of the GOCE time-wise (TIM) and GOCE space-wise gravity field series—two of the official ESA products—the energy conservation method has been adopted to exploit GPS-based satellite-to-satellite tracking information. On the other hand, gravity field recovery from data collected by the CHAllenging Mini-satellite Payload (CHAMP) satellite showed the energy conservation principle to be a sub-optimal choice. For this reason, we propose to estimate the low-frequency part of the gravity field by the point-wise solution of Newton’s equation of motion, also known as the acceleration approach. This approach balances the gravitational vector with satellite accelerations, and hence is characterized by (second-order) numerical differentiation of the kinematic orbit. In order to apply the method to GOCE, we present tailored processing strategies with regard to low-pass filtering, variance–covariance information handling, and robust parameter estimation. By comparison of our GIWF solutions (initials GI for “Geodätisches Institut” and IWF for “Institut für WeltraumForschung”) and the GOCE-TIM estimates with a state-of-the-art gravity field solution derived from GRACE (Gravity Recovery And Climate Experiment), we conclude that the acceleration approach is better suited for GOCE-only gravity field determination as opposed to the energy conservation method.     

Spin parameters of nanosatellite BLITS determined from Graz 2 kHz SLR data

The nanosatellite BLITS (Ball Lens In The Space) is the first object designed as a passive, spherical retroreflector of the Luneburg type, dedicated for Satellite Laser Ranging (SLR). The 2 kHz SLR station Graz measures spin parameters of this satellite, providing information about the rotational dynamics of the body. The measurements obtained during the period from September 26, 2009 to November 24, 2010 show a significant change of the spin configuration. The spin axis was dynamically precessing since the launch and currently is sinus-like behaving between coordinates RA 120°…150°, Dec 30°…60° (J2000 inertial reference frame). The angle between the symmetry axis and the spin axis of BLITS is not constant, but is decreasing since the launch, while its spin period is rather stable with a mean value of 5.613 s (clockwise rotation). The satellite was dynamically changing its attitude during the first three months after deployment; after this time the spin parameters are relatively stable.     

月球中继通信卫星系统发展综述与展望

中继通信是一些月球探测任务中必须解决的关键问题,特别是月球背面和两极地区的着陆和巡视探测任务以及载人登月任务。对国内外月球中继通信卫星的研究和发展情况进行了综述,在对月球中继通信任务需求分析的基础上,从中继通信体制选择、轨道选择等方面对月球中继通信任务的实现途径进行了分析,并对月球中继通信技术的未来发展进行了展望。     

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.