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.     

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.     

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.     

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.     

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.     

基于引力场不对称性的三体系统轨道自主导航

以月球背面的中继通信为背景,提出了基于三体系统引力场不对称特性的星–星测距自主定轨方案。该方案以环月极轨卫星和地–月L2点Halo轨道卫星组成中继通信网,以实现对月球两极和背面的覆盖。通过采集极轨卫星与Halo轨道卫星的测距信息,结合卡尔曼滤波在日–地–月动力学模型下获得两颗卫星的绝对轨道。数值仿真结果表明:本文方法能将导航的位置精度和速度精度分别提高到百米和厘米/秒量级。该自主导航方法还可以扩展到不规则引力场小天体附近星群运动的自主导航。     

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.     

联合北斗导航与星间链路的大椭圆卫星定轨方法

传统的地面测控和GNSS均无法实现HEO卫星全弧段的跟踪观测.在分析北斗导航信号及其星间链路信号对典型HEO的观测几何及覆盖特性的基础上,利用北斗导航及其星间链路对HEO测控支持形成互补的特点,提出了一种卫星导航与星间链路相结合的自主导航方法.对HEO定轨进行分段划分并基于EKF设计了卫星导航与星间链路数据融合定轨的自主导航算法.分析结果表明,本文提出的方法能够从全弧段上改善HEO的观测几何,定轨精度比仅使用卫星导航提高了2个数量级,并且仅需较少的星间链路资源.      

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.     

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.     

Orbit determination of BDS-3 satellite based on regional ground tracking station and inter-satellite link observations

The BeiDou global navigation satellite system (BDS-3) has established the Ka-band inter-satellite link (ISL) to realize a two-way ranging function between satellites, which provides a new observation technology for the orbit determination of BDS-3 satellites. Therefore, this study presents a BDS satellite orbit determination model based on ground tracking station (GTS) observations and ISL ranging observations firstly to analyze the impact of the ISL ranging observations on the orbit determination of BDS-3 satellites. Subsequently, considering the data fusion processing, the variance component estimation (VCE) algorithm is applied to the parameter estimation process of the satellite orbit determination. Finally, using the measured data from China’s regional GTS observations and BDS-3 ISL ranging observations, the effects of ISL ranging observations on the orbit determination accuracy of BDS-3 satellites are analyzed. Moreover, the impact of the VCE algorithm on the fusion data processing is evaluated from the aspects of orbit determination accuracy, Ka-band hardware delay parameter stability, and ISL ranging observation residuals. The results show that for China’s regional GTSs, the addition of BDS-3 ISL ranging observations can significantly improve the orbit determination accuracy of BDS-3 satellites. The observed orbit determination accuracy of satellite radial component is improved from 48 cm to 4.1 cm. In addition, when the initial weight ratio between GTS observations and ISL ranging observations is not appropriate, the various indicators which include orbit determination accuracy, ISL hardware delay, and ISL observation residuals were observed to have improved after the adjustment of the VCE algorithm. These results validate the effectiveness of the VCE algorithm for the fusion data processing of the GTS observations and ISL ranging observations.     

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.     

Formation control using μ-synthesis for Inner-Formation Gravity Measurement Satellite System

Inner-Formation Gravity Measurement Satellite System (IFGMSS) is used to map the gravity field of Earth. The IFGMSS consists of two satellites in which one is called “inner satellite” and the other one is named as “outer satellite”. To measure the pure Earth gravity, the inner satellite is located in the cavity of the outer satellite. Because of the shield effect of the cavity, the inner satellite is affected only by the gravitational force, so it can sense Earth gravity precisely. To avoid the collision between the inner satellite and the outer satellite, it is best to perform a real-time control on the outer satellite. In orbit, the mass of the outer satellite decreases with the consumption of its propellant. The orbit angular rate of the inner satellite varies with time due to various disturbing forces. These two parameters’ uncertainties make the C–W function be not so accurate to describe the formation behavior of these two satellites. Furthermore, the thrusters also have some uncertainties due to the unmodelled dynamics. To cancel the effects caused by the above uncertainties, we have studied the robust control method based on the μ-synthesis. This μ-synthesis eliminates the conservativeness and improves the control efficiency comparing with the H_∞ method. Finally, to test the control method, we simulate an IFGMSS mission in which the satellite runs in a sun synchronous circular orbit with an altitude of 300 km. The simulation results show the effectiveness of the robust control method. The performances of the closed-loop system with the μ-controller are tested by the μ-analysis. It has found that the nominal performance, the robust stability and the robust performance are all achieved. The transient simulation results further prove the control response is fast and the accuracy of the relative position meets the demand of the gravity measurement.     

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.     

北斗导航星座星间通信速率控制方法

北斗导航星座可以通过星间测距和传输链路实现自主定轨和性能增强。导航卫星间进行数据传输时,卫星相对位置时变,传输信道特性也随之不断发生变化。针对导航卫星间传输链路时变特性,提出了一种基于星历的星间通信速率控制方法。在满足传输服务质量的需求下,根据导航卫星自有的高精度星历资源定量计算星间最优通信速率,通过速率的动态调整提高星座的传输效能。仿真结果表明,采用所提方法,北斗导航星座星间传输效能可以提高1.92倍,验证了方法的有效性。      

卫星编队飞行的相对轨道自主确定研究

为了研究卫星编队飞行相对轨道的自主确定,基于相对轨道根数建立编队卫星间的相对运动方程,利用测量所得到的星间距离和方位信息作为观测量。不同于目前广泛采用的扩展卡尔曼滤波算法,设计Unscented Kalman Filter(UKF)算法实现卫星编队飞行的相对轨道自主确定。仿真结果表明这种相对轨道自主确定方案能获得较高的定轨精度。     

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.     

The influence of orbital maneuver on autonomous orbit determination of an extended satellite navigation constellation

The right ascension of the ascending node is unobservable if only the inter-satellite ranging is used for autonomous orbit determination (AOD) of an Earth navigation constellation. However, if an Earth-Moon libration point satellite is added to the Earth navigation constellation to construct an extended navigation constellation, all the orbital elements can be determined with only the inter-satellite ranging. Furthermore, the extended navigation constellation can provide navigation information for interplanetary probes. For such an extended navigation constellation, orbital control needs to be considered due to the instability of the libration-point satellite orbit. This study concerns the influence of satellite orbital maneuver on the AOD of the extended navigation constellation. An AOD method under orbital maneuver is proposed. A low thrust controller is designed to achieve libration point satellite autonomous orbit maintenance by using AOD results. A navigation constellation consisting of three GPS satellites and one libration point satellite are designed for simulation. The simulation results show that libration point satellites can achieve autonomous navigation and autonomous orbit maintenance by only using inter-satellite ranging information. The rotation drift error of the Earth navigation constellation is also suppressed.