深空卫星重力测量计划研究综述

地球卫星重力测量计划CHAMP(CHAllenging Minisatellite Payload)、GRACE(Gravity Recovery and Climate Experiment)、GOCE(Gravity field and steady-state Ocean Circulation Explorer)和月球卫星重力测量计划(Gravity Recovery and Interior Laboratory,GRAIL)的成功实施,以及下一代地球重力卫星(GRACE Follow-On)的即将发射昭示着我们将迎来一个前所未有的高精度和高空间分辨的深空卫星重力探测时代。围绕深空卫星重力测量的研究背景、必要性、可行性、卫星重力反演软件平台构建、轨道摄动和未来研究方向开展了研究论证。研究表明:深空卫星重力测量作为新世纪重力探测技术,在精化量体重力场、提高惯性导航精度、天体动力学、天体物理学和军事技术的研究,以及促进国民经济发展和提高社会效益等方面具有广泛的应用前景。     

Improved kHz-SLR tracking techniques and orbit quality analysis for LEO missions

Satellite gravity field missions such as CHAMP, GRACE and GOCE are designed as low Earth orbiting spacecraft (LEO) with orbit heights of about 250–500 km. The challenging mission objectives require a very precise knowledge of the satellite orbit position in space. For these missions precise orbit information is typically provided by GPS satellite-to-satellite tracking (SST) observations supported by satellite laser ranging (SLR).     

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.     

Measuring atmospheric density using GPS–LEO tracking data

We present a method to estimate the total neutral atmospheric density from precise orbit determination of Low Earth Orbit (LEO) satellites. We derive the total atmospheric density by determining the drag force acting on the LEOs through centimeter-level reduced-dynamic precise orbit determination (POD) using onboard Global Positioning System (GPS) tracking data. The precision of the estimated drag accelerations is assessed using various metrics, including differences between estimated along-track accelerations from consecutive 30-h POD solutions which overlap by 6 h, comparison of the resulting accelerations with accelerometer measurements, and comparison against an existing atmospheric density model, DTM-2000. We apply the method to GPS tracking data from CHAMP, GRACE, SAC-C, Jason-2, TerraSAR-X and COSMIC satellites, spanning 12 years (2001–2012) and covering orbital heights from 400 km to 1300 km. Errors in the estimates, including those introduced by deficiencies in other modeled forces (such as solar radiation pressure and Earth radiation pressure), are evaluated and the signal and noise levels for each satellite are analyzed. The estimated density data from CHAMP, GRACE, SAC-C and TerraSAR-X are identified as having high signal and low noise levels. These data all have high correlations with anominal atmospheric density model and show common features in relative residuals with respect to the nominal model in related parameter space. On the contrary, the estimated density data from COSMIC and Jason-2 show errors larger than the actual signal at corresponding altitudes thus having little practical value for this study. The results demonstrate that this method is applicable to data from a variety of missions and can provide useful total neutral density measurements for atmospheric study up to altitude as high as 715 km, with precision and resolution between those derived from traditional special orbital perturbation analysis and those obtained from onboard accelerometers.     

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.     

AIUB-CHAMP02S: The influence of GNSS model changes on gravity field recovery using spaceborne GPS

The gravity field model AIUB-CHAMP02S, which is based on six years of CHAMP GPS data, is presented here. The gravity field parameters were derived using a two step procedure: In a first step a kinematic trajectory of a low Earth orbiting (LEO) satellite is computed using the GPS data from the on-board receiver. In this step the orbits and clock corrections of the GPS satellites as well as the Earth rotation parameters (ERPs) are introduced as known. In the second step this kinematic orbit is represented by a gravitational force model and orbit parameters.     

Finding the suitable drag-free acceleration noise level for future low-low satellite-to-satellite tracking geodesy missions

This paper evaluates the impact of residual acceleration noise on the estimation of the Earth’s time-varying gravity field for future low-low satellite-to-satellite tracking missions. The goal is to determine the maximum level of residual acceleration noise that does not adversely affect the estimation error. The Gravity Recovery And Climate Experiment (GRACE) has provided monthly average gravity field solutions in spherical harmonic coefficients for more than a decade. It provides information about land and ocean mass variations with a spatial resolution of ～350?km and with an accuracy within 2?cm throughout the entire Earth. GRACE Follow-on was launched in May 2018 to advance the work of GRACE and to test a new laser ranging interferometer, which measures the range between the two satellites with higher precision than the K-Band ranging system used in GRACE. Moreover, there have been simulation studies that show, an additional pair of satellites in an inclined orbit increases the sampling frequency and reduces temporal aliasing errors. Given the fact that future missions will likely continue to use the low-low satellite-to-satellite tracking formation with laser ranging interferometry, it is expected that the residual acceleration noise will become one of the largest error contributor for the time-variable gravity field solution. We evaluate three different levels of residual acceleration noise based on demonstrated drag-free systems to find a suitable drag-free performance target for upcoming geodesy missions. We analyze both a single collinear polar pair and the optimal double collinear pair of drag-free satellites and assume the use of a laser ranging interferometer. A partitioned best linear unbiased estimator that was developed, incorporating several novel features from the ground up is used to compute the solutions in terms of spherical harmonics. It was found that the suitable residual acceleration noise level is around 2?×?10^?12?ms^?2?Hz^?1/2. Decreasing the acceleration noise below this level did not result in more accurate gravity field solutions for the chosen mission architecture.     

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.     

Characteristics and accuracies of the GRACE inter-satellite pointing

For almost 10 years, the Gravity Recovery and Climate Experiment (GRACE) has provided information about the Earth gravity field with unprecedented accuracy. Efforts are ongoing to approach the GRACE baseline accuracy as there still remains an order of magnitude between the present error level of the gravity field solutions and the GRACE baseline. At the current level of accuracy, thorough investigation of sensor related effects is necessary as they are one of the potential contributors to the error budget. In the science mode operations, the twin satellites are kept precisely pointed with their KBR antennas towards each other. It is the task of the onboard attitude and orbit control system (AOCS) to keep the satellites in the required formation. We analyzed long time series of the inter-satellite pointing variations as they reflect the AOCS performance and characteristics. We present significant systematic effects in the inter-satellite pointing and discuss their possible sources. Prominent features are especially related to the magnetic torquer characteristics, star cameras’ performance and KBR antenna calibration parameters. The relation between the magnetic torquer attitude control and the Earth magnetic field, impact of the different performance of the two star camera heads on the attitude control and the features due to uncertainties in the calibration parameters relating the star camera frame to K-frame are discussed in detail. Proper understanding of these effects will help to reduce their impact on the science data and subsequently increase the accuracy of the gravity field solutions. Moreover, understanding the complexity of the onboard system is essential not only for increasing the accuracy of the GRACE data but also for the development of the future gravity field satellite missions.     

Recovery of the Earth’s gravity field from formation-flying satellites: Temporal aliasing issues

Temporal and mean gravity field models derived from the twin-satellite, leader–follower mission GRACE have provided a new type of information for Earth sciences. In this contribution, we study the potential of various alternative satellite formations for gravity field determination in the post-GRACE era in a simulation environment. In particular, the effects of spherical harmonic truncation and of temporal aliasing in the processing of gravity products from such future formations are investigated.     

重力卫星星间高精度测距技术研究

卫卫跟踪(SST)技术是当前地球重力场测量最有价值和应用前景的方法之一.高精度星间测距系统是低低卫卫跟踪(SST-ll)重力卫星的关键有效载荷.GRACE卫星携带的K波段测距系统(KBR K Band Ranging System)是一微米量级的测距系统,通过处理高精度的星间距离和距离变化率数据,可以恢复出地球重力场.GRACE后续计划又提出了一种更高精度的激光干涉测距系统.在研究KBR及激光干涉测距系统测量原理的基础上,提出了一种KBR系统的基本结构,详细分析了两种测距系统的关键技术及国内目前的研究水平,提出了我国开展星间测距系统研究的一些建议.      

Magnetic torquer induced disturbing signals within GRACE accelerometer data

The GRACE (Gravity Recovery And Climate Experiment) gravity field satellite mission was launched in 2002. Although many investigations have been carried out, not all disturbances and perturbations upon satellite instruments and sensors are resolved yet. In this work the issue of acceleration disturbances onboard of GRACE due to magnetic torquers is investigated and discussed. Each of the GRACE satellites is equipped with a three-axes capacitive accelerometer to measure non-gravitational forces acting on the spacecraft. We used 10 Hz Level 1a raw accelerometer data in order to determine the impact of electric current changes on the accelerometer. After reducing signals which are induced by highly dominating processes in the low frequency range, such as thermospheric drag and solar radiation pressure, which can easily be done by applying a high-pass filter, disturbing signals from onboard instruments such as thruster firing events or heater switch events need to be removed from the previously filtered data. Afterwards the spikes which are induced by the torquers can be very well observed. Spikes vary in amplitude with respect to an increasing or decreasing current used for magnetic torquers, and can be as large as 20 nm/s². Furthermore, we were able to set up a model for the spikes of each scenario with which we were able to compute model spike time series. With these time series the spikes can successfully be removed from the 10 Hz raw accelerometer data. Spectral analysis of the time series reveal that an influence onto gravity field determination due to these effects is very unlikely, but can theoretically not be excluded.     

Performance of GPS-based accelerometry: A simulation experiment

Recent studies have shown that with the availability of high-quality CHAMP and GRACE gravity field models, it is feasible to determine accurate non-gravitational accelerations for low Earth orbiting satellites indirectly from precise GPS satellite-to-satellite observations. Possible applications of this so-called GPS-based accelerometry approach consist of accelerometer calibration and atmospheric density and wind computations. With the growing number of high-quality space-borne GPS receivers, this method could be applied to a large range of satellites. In this paper an extensive simulation study has been carried out, based on real accelerometer data from the GRACE mission, in order to determine the optimal processing strategy and the resulting accuracy of the estimated non-gravitational accelerations. It is shown that the optimal processing strategy consists of a piecewise linear parameterization of the estimated empirical accelerations, together with short 6-h orbit arcs. The GPS-based accelerometry approach makes use of triple-differenced GPS observations and the impact of considering the correlated observation noise was found to be marginal in the presence of other error sources such as GPS ephemeris errors. Using a priori non-gravitational force models improves the recovery of low temporal resolution accelerations, except during huge geomagnetic storms. With this strategy, non-gravitational accelerations can be recovered during high solar activity with an accuracy of better than 10% of the total signal in along-track direction and around 25–40% in cross-track direction, at time resolutions of around 8–20 min. During solar minimum conditions, the relative recovery error will increase to approximately 50% in along-track direction and around 60–70% in cross-track direction, due to the reduced atmospheric drag signal. Unfortunately, GPS-based accelerometry is hardly sensitive in the radial direction.     

NRLMSISE-00大气模型与GRACE和CHAMP卫星大气密度数据的对比分析

利用GRACE(Gravity Recovery And Climate Experiment)和CHAMP(Challenging Mini-Satellite Payload)卫星2002-2008年的大气密度数据与NRLMSISE-00大气模型密度结果进行比较,分析了模型密度误差及其特点.结果显示,NRLMSISE-00大气模型计算的密度值普遍偏大,其相对误差随经纬度变化,在高纬度相对较小;相对误差随地方时变化,在02:00LT和15:00LT左右较大,10:00LT和20:00LT左右较小.通过模型密度相对误差与太阳F_10.7指数的对比分析发现,在太阳活动低年模型相对误差最大,而在太阳活动高年相对误差较小;将模型结果分别与GRACEA/B双星和CHAMP卫星的密度数据进行比较,发现对于轨道高度更高的GRACE卫星轨道,模型相对误差更大;在地磁平静期,相对误差与地磁ap指数(当前3h)相关性不强,但是在大磁暴发生时,误差急剧增大.      

A gravity wave analysis near to the Andes Range from GPS radio occultation data and mesoscale numerical simulations: Two case studies

Global maps of potential wave energy per unit mass, recently performed with the Global Positioning System (GPS) Radio Occultation (RO) technique and different satellite missions (CHAMP and SAC-C since 2001, GRACE and COSMIC since 2006) revealed in Argentina, at the eastern side of the highest Andes Mountains, a considerable wave activity (WA) in comparison with other extra-tropical regions. The main gravity wave (GW) sources in this natural laboratory are deep convection (mainly during late Spring and Summer), topographic forcing and geostrophic adjustment.     

Impact of covariance information of kinematic positions on orbit reconstruction and gravity field recovery

Gravity missions are equipped with onboard Global Positioning System (GPS) receivers for precise orbit determination (POD) and for the extraction of the long wavelength part of the Earth’s gravity field. As positions of low Earth orbiters (LEOs) may be determined from GPS measurements at each observation epoch by geometric means only, it is attractive to derive such kinematic positions in a first step and to use them in a second step as pseudo-observations for gravity field determination. The drawback of not directly using the original GPS measurements is, however, that kinematic positions are correlated due to the ambiguities in the GPS carrier phase observations, which in principle requires covariance information be taken into account. We use GRACE data to show that dynamic or reduced-dynamic orbit parameters are not optimally reconstructed from kinematic positions when only taking epoch-wise covariance information into account, but that essentially the same orbit quality can be achieved as when directly using the GPS measurements, if correlations in time are taken into account over sufficiently long intervals. For orbit reconstruction covariances have to be considered up to one revolution period to avoid ambiguity-induced variations of kinematic positions being erroneously interpreted as orbital variations. For gravity field recovery the advantage is, however, not very pronounced.     

High-fidelity geometry models for improving the consistency of CHAMP,GRACE, GOCE and Swarm thermospheric density data sets

During the last two decades, accelerometers on board of the CHAMP, GRACE, GOCE and Swarm satellites have provided high-resolution thermosphere density data to improve our knowledge on atmospheric dynamics and coupling processes in the thermosphere-ionosphere region. Most users of the data have focused on relative density variations. Scale differences between datasets and models have been largely neglected or removed using ad hoc scale factors. The origin of these scale differences arises from errors in the aerodynamic modelling, specifically in the modelling of the satellite outer surface geometry and of the gas-surface interactions. Therefore, the first step to remove the scale differences is to enhance the geometry modelling. This work forms the foundation for the future improvement of characterization of satellite aerodynamics and gas-surface interactions models at TU Delft, as well as for extending the use of sideways and angular accelerations in the aerodynamic analysis of accelerations and derivation of thermosphere datasets. Although work to improve geometry and aerodynamic force models by other authors has focused on CHAMP and GRACE, this paper includes the GOCE and Swarm satellites as well. In addition, it uses a density determination algorithm that is valid for arbitrary attitude orientations, enabling a validation making use of attitude manoeuvres. The results show an improvement in the consistency of density data between these four missions, and of data obtained before, during and after attitude manoeuvres of CHAMP and Swarm. The new models result in larger densities, compared to the previously used panel method. The largest average rescaling of density, by switching to the new geometry models is reached for Swarm at 32%, the smallest for GRACE at 5%. For CHAMP and GOCE, mean differences of 11% and 9% are obtained respectively. In this paper, an overview of the improvements and comparisons of data sets is provided together with an introduction to the next research phase on the gas-surface interactions.     

GPS single-frequency orbit determination for low Earth orbiting satellites

The determination of high-precision orbits for Low Earth Orbiting (LEO) satellites (e.g., CHAMP, GRACE, MetOp-A) is based on dual-frequency tracking data from on-board GPS (Global Positioning System) receivers. The two frequencies allow it to eliminate the first order ionosphere effects. Data screening and precise orbit determination (POD) procedures are optimized under the assumption of the availability of two frequencies.     

ERS-1 radial position refinement by dual crossover analysis with TOPEX/Poseidon

Achieving orbital accuracies in the radial direction for ERS-1 commensurate with those for TOPEX/Poseidon is of utmost importance for the integration of the two altimeter data sets. This paper outlines a procedure whereby the radial orbit error for ERS-1 is recovered as a time series expansion in the form of a finite Fourier series with additional terms for atmospheric drag, solar radiation pressure, and initial state vector mismodelling. Using a least squares collocation method with constraints derived from the JGM2 gravity field co-variance matrix, the radial error is recovered using both dual crossovers and ERS-1 single satellite crossovers. Aggregate arcs are then used to derive the ERS-1 orbit error over the repeat period of 35 days. The results are presented in the improvement of fit in the dual crossover, ERS-1 crossover and altimetry data sets as well as the recovery of an altimeter bias for the two satellites.     

重力卫星高精度测距系统观测量的降速率滤波算法

在卫星时频传递、矢量测量技术领域，有效观测数据的高精度滤波处理十分重要，其基本要求是在保留有效观测信息的同时能够降低噪声干扰，实现原始信息的高效、高精度传递。以中国地球重力场测量卫星高精度测距系统的工程研制为背景，为实现高频噪声抑制、保留低频重力场信息，基于对FIR滤波器结构的优化和窗函数的改进，提出了一种对10 Hz采样率的原始测距数据降速率、降噪滤波的算法。相比其他类似算法，所提方法同时具备降速滤波和差分运算两种功能，噪声抑制抑制度优于44%，且频率截止特性陡峭，第一旁瓣抑制度达到–94 dB。该算法在提高测距精度、充分保留重力场信息的同时，实现了地面数据后处理流程的简化，为微波测量系统、激光测量系统的数据处理提供了一种优选方案。