Dynamic and non-dynamic ERS-1 radial orbit improvement from ERS-1/TOPEX dual-satellite altimetry

Dual-satellite altimeter crossover differences between ERS-1 and TOPEX/Poseidon have been included as supplementary tracking data in ERS-1 orbit computations from SLR and single-satellite crossover differences. It was found that including the dual-satellite crossover differences slightly improves the ERS-1 radial orbit accuracy of about 12 cm for orbits computed with the JGM-2 gravity field and also leads to a better ‘centering’ of the ERS-1 orbit in the terrestrial reference frame defined for TOPEX/Poseidon. In addition to this dynamic orbit improvement technique, a non-dynamic technique has been investigated that removes the larger part of the ERS-1 radial orbit error from the dual-satellite crossover difference residuals. For ERS-1 orbits computed with the GEM-T2 gravity field, it was found that the non-dynamic technique could improve the radial orbit accuracy from 140 cm to the same level of accuracy as the ERS-1 JGM-2 orbits.     

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.     

Earth gravity model improvement: An alternative method for Doppler-tracked satellites

A new method of Earth gravity model improvement based on an analytical formulation of Doppler residuals is presented here in prospect of future geodetic and altimetric missions (DORIS, TOPEX/POSEIDON, ERS1). After an intermediate step of orbit improvement, disturbing forces due to gravity field mismodeling are recovered above tracking stations at satellite altitude. Some significant simulation results for SEASAT and DORIS are presented.     

Precision orbit determination performance for CryoSat-2

In this paper we discuss our efforts to perform precision orbit determination (POD) of CryoSat-2 which depends on Doppler and satellite laser ranging tracking data. A dynamic orbit model is set-up and the residuals between the model and the tracking data is evaluated. The average r.m.s. of the 10?s averaged Doppler tracking pass residuals is approximately 0.39?mm/s; and the average of the laser tracking pass residuals becomes 1.42?cm. There are a number of other tests to verify the quality of the orbit solution, we compare our computed orbits against three independent external trajectories provided by the CNES. The CNES products are part of the CryoSat-2 products distributed by ESA. The radial differences of our solution relative to the CNES precision orbits shows an average r.m.s. of 1.25?cm between Jun-2010 and Apr-2017. The SIRAL altimeter crossover difference statistics demonstrate that the quality of our orbit solution is comparable to that of the POE solution computed by the CNES. In this paper we will discuss three important changes in our POD activities that have brought the orbit performance to this level. The improvements concern the way we implement temporal gravity accelerations observed by GRACE; the implementation of ITRF2014 coordinates and velocities for the DORIS beacons and the SLR tracking sites. We also discuss an adjustment of the SLR retroreflector position within the satellite reference frame. An unexpected result is that we find a systematic difference between the median of the 10 s Doppler tracking residuals which displays a statistically significant pattern in the South Atlantic Anomaly (SSA) area where the median of the velocity residuals varies in the range of ?0.15 to +0.15?mm/s.     

Gravity model development for precise orbit computations for satellite altimetry

Satellite altimeters have been developed which have the capability of measuring the surface topography of the ocean over entire basins for time periods of several years with a precision of a few cm. The combination of these data with accurate satellite position information, subsurface measurements and models of the density field can be used to determine the general circulation of the ocean and its variability. The Topex altimeter oceanographic experiment proposed for the early 1990's requires a radial orbit accuracy approaching 10 cm. Considerable progress has been made at the Goddard Space Flight Center in developing precision orbit determination to support satellite altimetry analysis. This progress is attributed to the development of more accurate models of the earth gravity field, earth and ocean tidal effects, a reference coordinate system as well as improved drag and solar radiation pressure computations. This computational system has been used for the Seasat data acquired during the 1978 time period and has produced an rms radial accuracy of about 50 cm. A first step towards the development of a new Interim gravity model is presented. Global orbit determination accuracies have been improved using these new PGS-T1 and T2 fields.     

SPOT-2 and TOPEX/Poseidon precise orbit determination from DORIS doppler tracking

The French earth observation satellite SPOT-2 has served as a testbed for precise orbit determination from DORIS doppler tracking in anticipation of the TOPEX/Poseidon mission. Using the most up-to-data gravity field model, JGM-2, a radial orbit accuracy of about 2–9 cm was achieved, with an rms of fit of the tracking data of about 0.64 mm/s. Furthermore, it was found that the coordinates of the ground stations can be determined with an accuracy of the order of 2–5 cm after removal of common rotations, and translations.

Using a slightly different model for atmospheric drag, but the same gravity model, precise orbits of TOPEX/Poseidon from DORIS tracking data were determined with a radial orbit accuracy of the order of 4–5 cm, which is far within the 13 cm mission requirement. This conclusion is based on the analysis of 1-day overlap of successive 11-day orbits, and the comparisons with orbits computed from satellite laser tracking (SLR) and from the combination of SLR and DORIS tracking. Results indicate a consistency between the different orbits of 1–4 cm, 4–20 cm, and 6–13 cm in the radial, cross-track, and along-track directions, respectively. The residual rms is about 4–5 cm for SLR data and 0.56 mm/s for DORIS tracking. These numbers are roughly twice as large as the system noise levels, reflecting the fact that there are still some modeling errors left.     

Gravity field model adjustment from ERS-1 and TOPEX altimetry and SLR tracking using an analytical orbit perturbation theory

The JGM-2 gravity field model has been adjusted using 70 days of ERS-1 and 19 10-day repeat cycles of TOPEX/Poseidon SLR and single satellite altimeter crossover differences. In addition, dual satellite altimeter crossover differences between ERS-1 and TOPEX for the selected 70-day period of ERS-1 and TOPEX repeat cycle 18 have been used in the JGM-2 model adjustment. In the computation of the normal equations, use was made of the analytical Lagrange linear perturbation theory. The single satellite altimeter crossover difference rms was brought down from 15.2 to 14.1 cm for ERS-1 and from 10.6 to 10.5 cm for TOPEX. The dual satellite altimeter crossover difference rms was reduced from 17.7 to 16.9 cm. Furthermore, the weighted rms of fit of SLR measurements was brought down from 16.0 to 14.4 cm for ERS-1, and from 5.4 to 5.0 cm for TOPEX/Poseidon.     

Precision orbit determination for TOPEX/Poseidon using TDRSS doppler tracking data

Precision orbit determination on the TOPEX/Poseidon (T/P) altimeter satellite is now being routinely achieved with sub-5cm radial and sub-15 cm total positioning accuracy using state-of-the-art modeling with precision tracking provided by a combination of: (a) global Satellite Laser Ranging (SLR) and Doppler Orbitography and Radiopositioning Integrated by Satellite (DORIS), or (b) the Global Positioning System (GPS) Constellation which provides pseudo-range and carrier phase observations. The geostationary Tracking and Data Relay Satellite System (TDRSS) satellites are providing the operational tracking and communication support for this mission. The TDRSS Doppler data are of high precision (0.3 mm/s nominal noise levels). Unlike other satellite missions supported operationally by TDRSS, T/P has high quality independent tracking which enables absolute orbit accuracy assessments. In addition, the T/P satellite provides extensive geometry for positioning a satellite at geostationary altitude, and thus the TDRSS-T/P data provides an excellent means for determining the TDRS orbits. Arc lengths of 7 and 10 days with varying degrees of T/P spacecraft attitude complexity are studied. Sub-meter T/P total positioning error is achieved when using the TDRSS range-rate data, with radial orbit errors of 10.6 cm and 15.5 cm RMS for the two arcs studied. Current limitations in the TDRSS precision orbit determination capability include mismodeling of numerous TDRSS satellite-specific dynamic and electronic effects, and in the inadequate treatment of the propagation delay and bending arising from the wet troposphere and ionosphere.     

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.     

Impact of ITRS 2014 realizations on altimeter satellite precise orbit determination

This paper evaluates orbit accuracy and systematic error for altimeter satellite precise orbit determination on TOPEX, Jason-1, Jason-2 and Jason-3 by comparing the use of four SLR/DORIS station complements from the International Terrestrial Reference System (ITRS) 2014 realizations with those based on ITRF2008. The new Terrestrial Reference Frame 2014 (TRF2014) station complements include ITRS realizations from the Institut National de l’Information Géographique et Forestière (IGN) ITRF2014, the Jet Propulsion Laboratory (JPL) JTRF2014, the Deutsche Geodätisches Forschungsinstitut (DGFI) DTRF2014, and the DORIS extension to ITRF2014 for Precise Orbit Determination, DPOD2014. The largest source of error stems from ITRF2008 station position extrapolation past the 2009 solution end time. The TRF2014 SLR/DORIS complement impact on the ITRF2008 orbit is only 1–2 mm RMS radial difference between 1992–2009, and increases after 2009, up to 5 mm RMS radial difference in 2016. Residual analysis shows that station position extrapolation error past the solution span becomes evident even after two years, and will contribute to about 3–4 mm radial orbit error after seven years. Crossover data show the DTRF2014 orbits are the most accurate for the TOPEX and Jason-2 test periods, and the JTRF2014 orbits for the Jason-1 period. However for the 2016 Jason-3 test period only the DPOD2014-based orbits show a strong and statistically significant margin of improvement. The positive results with DTRF2014 suggest the new approach to correct station positions or normal equations for non-tidal loading before combination is beneficial. We did not find any compelling POD advantage in using non-linear over linear station velocity models in our SLR & DORIS orbit tests on the Jason satellites. The JTRF2014 proof-of-concept ITRS realization demonstrates the need for improved SLR+DORIS orbit centering when compared to the Ries (2013) CM annual model. Orbit centering error is seen as an annual radial signal of 0.4 mm amplitude with the CM model. The unmodeled CM signals show roughly a 1.8 mm peak-to-peak annual variation in the orbit radial component. We find the TRF network stability pertinent to POD can be defined only by examination of the orbit-specific tracking network time series. Drift stability between the ITRF2008 and the other TRF2014-based orbits is very high, the relative mean radial drift error over water is no larger than 0.04 mm/year over 1993–2015. Analyses also show TRF induced orbit error meets current altimeter rate accuracy goals for global and regional sea level estimation.     

ERS-1 and ERS-2 operational and precise orbit determination

This paper presents the European Space Operations Centre's orbit determination and prediction systems for the ERS-1 mission. The routine operational orbit determination and prediction subsystem is discussed briefly, and statistics of the accuracy compared to the requirements are given. The precise orbit determination subsystem is then described, and the accuracy of its results are compared to those of the operational orbit system and to the D-PAF preliminary orbit solutions. Some geophysical results from the altimeter data, processed in these orbit determinations, are also presented. The ESOC/OAD ‘ERS-1 Orbit Report’ is introduced as a document providing this information on a monthly basis. Finally, this paper describes how the experience gained with the precise orbit determination will be exploited to further improve the accuracy of the routine system that will be used for ERS-2, and provides an estimate of this accuracy.     

New improved orbit solutions for the ERS-1 and ERS-2 satellites

Improved orbit solutions of the European Remote Sensing Satellites ERS-1 and ERS-2 have been computed in the ITRF2005 terrestrial reference frame using the recent models based mainly on IERS Conventions 2003. These solutions cover the periods 3 August 1991 to 8 July 1996 for ERS-1, and 3 May 1995 to 4 July 2003 for ERS-2. For each satellite, the final orbit solution is based on a combination of three separate orbit solutions independently computed at the Delft Institute of Earth Observation and Space Systems (DEOS) of the Delft University of Technology (The Netherlands), the Navigation Support Office of the European Space Operations Centre (ESOC, Germany) and the Helmholtz Centre Potsdam GFZ German Research Centre for Geosciences (Germany) using three different software packages for precise orbit determination, but using the same models in the same terrestrial reference frame within the European Space Agency (ESA) project ‘Reprocessing of Altimeter Products for ERS (REAPER)’. Validation using radar altimeter data indicates that the new combined orbits of ERS-1 and ERS-2 computed by us are significantly more accurate, approaching the 2–3 cm level in radial direction, than previously available orbit solutions.     

Geostationary orbit determination using SATRE

A new strategy of precise orbit determination (POD) for GEO (Geostationary Earth Orbit) satellite using SATRE (SAtellite Time and Ranging Equipment) is presented. Two observation modes are proposed and different channels of the same instruments are used to construct different observation modes, one mode receiving time signals from their own station and the other mode receiving time signals from each other for two stations called pairs of combined observations. Using data from such a tracking network in China, the results for both modes are compared. The precise orbit determination for the Sino-1 satellite using the data from 6 June 2005 to 13 June 2005 has been carried out in this work. The RMS (Root-Mean-Square) of observing residuals for 3-day solutions with the former mode is better than 9.1 cm. The RMS of observing residuals for 3-day solutions with the latter mode is better than 4.8 cm, much better than the former mode. Orbital overlapping (3-day orbit solution with 1-day orbit overlap) tests show that the RMS of the orbit difference for the former mode is 0.16 m in the radial direction, 0.53 m in the along-track direction, 0.97 m in the cross-track direction and 1.12 m in the 3-dimension position and the RMS of the orbit difference for the latter mode is 0.36 m in the radial direction, 0.89 m in the along-track direction, 1.18 m in the cross-track direction and 1.52 m in the 3-dimension position, almost the same as the former mode. All the experiments indicate that a meter-level accuracy of orbit determination for geostationary satellite is achievable.     

基于UKF的雷达高度计自主定轨

探讨了利用推广卡尔曼滤波估计非线性系统状态时存在的问题,从而介绍了目前广泛使用的分步逼近的卡尔曼滤波(UKF,Unscented Kalman Filter).为了提高导航的可靠性和准确性,在星敏感器导航系统中引入雷达高度计作为一个新的测量设备,提出了一种基于星上雷达测高仪及星敏感器联合进行卫星自主定轨的算法.建立了比较复杂的地球海平面模型,并考虑了其中风生重力波的影响. 利用雷达测高仪的测量结果和地球形状模型,计算地心矢量在卫星本体中坐标系的方向.利用UKF滤波定轨算法,明显提高了自主定轨的精度.数值仿真结果表明,UKF定轨精度要远优于推广卡尔曼滤波.      

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.     

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).     

基于光学测角数据的风云四号同步轨道卫星精密定轨

针对风云四号同步卫星的精密定轨和精度评估需求,首先利用地面光学测角数据对FY-4A卫星进行精密定轨,定轨后方位角和高度角的残差rms分别为0.25＂和0.45＂。与基于测距数据的轨道相比,位置精度在有测角数据的弧段内小于50m。进一步联合测角数据和测距数据对FY-4A卫星进行联合定轨,定轨后轨道重叠精度优于15m。利用联合定轨结果评估了基于测距数据的实时轨道产品精度,可以明显发现轨道精度随着测距数据的积累而逐步提高。     

基于机动力模型的GEO卫星恢复期间定轨

北斗卫星导航系统(BDS)中GEO卫星频繁的轨道机动对高精度、实时不间断的导 航服务需求提出了更高要求, 如何在短弧跟踪条件下提高GEO卫星轨道快速 恢复能力, 是提升导航系统服务精度的关键因素. 针对该问题, 本文提出了基 于机动力模型的动力学定轨方法, 尝试利用高精度的C波段转发式测距数据, 辅 以机动期间的遥测遥控信息建立机动力模型, 联合轨控前后的观测数据进行动 力学长弧定轨. 利用BDS中GEO卫星实测数据进行了定轨试验与分析, 结果表明, 恢复期间需要采用解算机动推力的定轨方法, 联合机动前、机动期间和机 动后4h数据定轨的轨道位置精度在20m量级, 径向精度优于2.5m. 该方 法克服了短弧跟踪条件下动力学法定轨和单点定位中的诸多问题, 提供了解决 GEO卫星机动后轨道快速恢复问题的技术方法.      

Arctic geodynamics: The opening of the Arctic Ocean

A study of the SEASAT satellite altimetry data from the Arctic has revealed gravity anomalies which may be related to the initial events of the opening of the Arctic Ocean and formation of the margins. Because of the density of the satellite altimetry data available accurate 5′ × 5′ mean gravity anomalies could be formed and analysed.An interpretation of the gravity anomalies in terms of extensional processes during margin formation has been made and it appears that the data indicates the existence of large boundinage structures in the margin. These structures are discussed in terms of their significance to the interpretation of the original rifting processes prior to and during the opening of the Arctic seafloor.