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.     

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.     

Precise orbit determination using the batch filter based on particle filtering with genetic resampling approach

In this study, genetic resampling (GRS) approach is utilized for precise orbit determination (POD) using the batch filter based on particle filtering (PF). Two genetic operations, which are arithmetic crossover and residual mutation, are used for GRS of the batch filter based on PF (PF batch filter). For POD, Laser-ranging Precise Orbit Determination System (LPODS) and satellite laser ranging (SLR) observations of the CHAMP satellite are used. Monte Carlo trials for POD are performed by one hundred times. The characteristics of the POD results by PF batch filter with GRS are compared with those of a PF batch filter with minimum residual resampling (MRRS). The post-fit residual, 3D error by external orbit comparison, and POD repeatability are analyzed for orbit quality assessments. The POD results are externally checked by NASA JPL’s orbits using totally different software, measurements, and techniques. For post-fit residuals and 3D errors, both MRRS and GRS give accurate estimation results whose mean root mean square (RMS) values are at a level of 5 cm and 10–13 cm, respectively. The mean radial orbit errors of both methods are at a level of 5 cm. For POD repeatability represented as the standard deviations of post-fit residuals and 3D errors by repetitive PODs, however, GRS yields 25% and 13% more robust estimation results than MRRS for post-fit residual and 3D error, respectively. This study shows that PF batch filter with GRS approach using genetic operations is superior to PF batch filter with MRRS in terms of robustness in POD with SLR observations.     

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.     

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.     

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.     

10 Years of LAGEOS-1 and 15 years of LAGEOS-2 spin period determination from SLR data

Satellite Laser Ranging (SLR) stations measure distance to the satellites equipped with Corner Cube Reflectors (CCRs). These range measurements contain information about spin parameters of the spacecraft. In this paper we present results of spin period determination of two passive satellites from SLR data only: 10 years of LAGEOS-1 (10426 values), and 15 years of LAGEOS-2 (15580 values). The measurements have been made by standard 10 Hz SLR systems and the first 2 kHz SLR system from Graz (Austria). The obtained data allowed calculation of the initial spin period of the satellites: 0.61 s for LAGEOS-1 and 0.906 s for LAGEOS-2. Long time series of the spin period values show that the satellite’s slowing down rate is not constant but is oscillating with a period of 846 days for LAGEOS-1 and 578 days for LAGEOS-2. The results presented here definitely prove that the SLR is a very efficient technique able to measure spin period of the geodetic satellites.     

凸轮测量采点间距选取方法

对凸轮测量中如何选取采点间距的问题进行了分析 ,指出了传统的等角度间隔采点的不合理性 ,提出按被测凸轮测量误差的要求选取采点间距的方法 ,并推导出采点间距的计算公式。     

参数估计法在规则图形三维形貌测量中的应用

引入了参数估计法对规则图形三维形貌测量中的关键问题提高差分干涉系统横向定位准确度的方法进行了理论分析。仿真计算表明，在采用普通激光束径和工作台的条件下，系统横向定位准确度优于0.1μm，同时系统的复杂性和造价大大降低。     

Assessing attitude error of FORMOSAT-3/COSMIC satellites and its impact on orbit determination

An attitude determination and control system (ADCS) is critical to satellite attitude maneuvers and to the coordinate transformation from the inertial frame to the spacecraft frame. This paper shows specific sensors in the ADCS of the satellite mission FORMOSAT-3/COSMIC (F3/C) and the impact of the ADCS quality on orbit accuracy. The selection of main POD antenna depends on the beta angles of the different F3/C satellites (for FM2 and FM4) during the inflight phase. In particular, under the eclipse, alternative attitude sensors are activated to replace the Sun sensors, and such a sensor change leads to anomalous GPS phase residuals and a degraded orbit accuracy. Since the nominal attitude serves as a reference for ADCS, the 3-dimensional attitude-induced errors in reduced dynamic orbits over selected days in 2010 show 9.35, 10.78, 4.97, 5.48, 7.18, and 6.89 cm for FM1–FM6. Besides, the 3-dimensional velocity errors induced by the attitude effect are 0.10, 0.10, 0.07, 0.08, 0.09, and 0.10 for FM1–FM6. We analyze the quality of the observed attitude transformation matrix of F3/C and its impact on kinematic orbit determination. With 249 days of GPS in 2008, the analysis leads to the following averaged 3-dimensional attitude-induced orbit errors: 2.72, 2.62, 2.37, 1.90, 1.70, and 1.99 cm for satellites FM1–FM6. Critical suggestions of geodetic payloads for the follow-on mission of F3/C are presented based on the current result.     

Analysis on navigation accuracy of high orbit spacecraft based on BDS GEO and IGSO satellites

To make up for the insufficiency of earth-based TT&C systems, the use of GNSS technology for high-orbit spacecraft navigation and orbit determination has become a new technology. It is of great value to applying Geosynchronous Earth Orbit (GEO) and Inclined GeoStationary Orbit (IGSO) navigation satellites for supporting the navigation of high-orbit spacecraft since there are three different types of navigation satellites in BeiDou Navigation Satellite System (BDS): Medium Earth Orbit (MEO), GEO and IGSO. This paper conducts simulation experiments based on Two-Line Orbital Element (TLE) data to analyze and demonstrate the role of these satellites in the navigation of high-orbit spacecraft. Firstly, the spacecraft in GEO was used as the target satellite to conduct navigation experiments. Experiments show that for the spacecraft on the GEO orbit, after adding GEO and IGSO respectively on the basis of receiving MEO navigation satellite signals, the accuracies were improved by 7.22 % and 6.06 % respectively. When adding both GEO and IGSO navigation satellites at the same time, the accuracy can reach 16 m. In the second place, navigation and positioning experiments were carried out on three high elliptical orbit (HEO) satellites with different semimajor axis (32037.2 km, 42385.9 km, 67509.6 km). The experiments show that the number of visible satellites has been improved significantly after adding GEO and IGSO navigation satellites at the same time. The visible satellites in these three orbits were improved by 32.84 %, 41.12 % and 37.68 %, respectively compared with only observing MEO satellites.The RMS values of the navigation positioning errors of these three orbits are 25.59 m, 87.58 m and 712.48 m, respectively.     

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.     

Performance assessment of real-time orbit determination for the Haiyang-2D using onboard BDS-3/GPS observations

The Global Navigation Satellite System (GNSS) receivers equipped on the Haiyang-2D (HY-2D) satellite is capable of tracking the signals of both the third generation of BeiDou satellite navigation System (BDS-3) and the Global Positioning System (GPS), which make it feasible to assess the performance of real-time orbit determination (RTOD) for the HY-2D using onboard GNSS observations. In this study, the achievable accuracy and convergence time of RTOD for the HY-2D using onboard BDS-3 and GPS observations are analyzed. Benefiting from the binary-offset-carrier (BOC) modulation, the BDS-3 C1X signal includes less noise than the GPS C1C signal, which has the same signal frequency and chipping rate. The root mean squares (RMS) of the noises of C1X and C1C code measurements are 0.579 m and 1.636 m, respectively. Thanks to a ten-times higher chipping rate, the code measurements of BDS-3 C5P, GPS C1W and C2W are less noisy. The RMS of code noises of BDS-3 C5P, GPS C1W, and C2W are 0.044 m, 0.386 m, and 0.272 m, respectively. For the HY-2D orbit, the three-dimensional (3D) and radial accuracies can reach 31.8 cm and 7.5 cm with only BDS-3 observations, around 50 % better than the corresponding accuracies with GPS. Better performance of the BDS-3 in RTOD for the HY-2D is attributed to the high quality of its broadcast ephemeris. When random parameters are used to absorb ephemeris errors, substantial improvement is seen in the accuracy of HY-2D orbit with either BDS-3 or GPS. The 3D RMS of HY-2D orbit errors with BDS-3 and GPS are enhanced to 23.1 cm and 33.6 cm, and the RMS of the radial components are improved to 6.1 cm and 13.3 cm, respectively. The convergence time is 41.6 and 75.5 min for the RTOD with BDS-3 and GPS, while it is reduced to 39.2 and 27.4 min after the broadcast ephemeris errors are absorbed by random parameters. Overall, the achievable accuracy of RTOD with BDS-3 reaches decimeter level, which is even better than that with GPS, making real-time navigation using onboard BDS-3 observations a feasible choice for future remote sensing missions.     

Numerical orbit integration based on Lie series with use of parallel computing techniques

This article outlines necessary steps to perform numerical orbit integrations based on a Lie series approach. Its implementation requires an efficient evaluation of resulting series coefficients. As an example we treat the classical main problem in satellite orbit calculation (_J2$J_2$ only) and the case of a 4×4$4\times 4$-gravity field. All calculations were performed in very high precision with up to 100 significant digits. In comparison to independent third party computations this approach led to superior results referring to the verifiable constancy of various integrals of motion. To achieve a performance similar to classical numerical integrations in terms of acceptable computing time, at least for non-Keplerian motion problems, we exploited parallel computing capabilities. For our examples, run times were improved by several orders of magnitude, depending on the actual chosen precision level (up to a factor of 50,000 in case of double precision). Here we present the mathematical framework of the proposed orbital integration scheme as well as the work flow for its application in a multi-core, parallel computing environment.     

A model-free method for attitude estimation and inertial parameter identification of a noncooperative target

Determining the attitude and inertial parameters of a noncooperative target is essential in an on-orbit servicing mission. Various methods based on machine vision have been proposed, but most of them require the 3D model of the target. This paper proposes a model-free method through sequentially registering point clouds captured by a depth camera. Our main contributions are the avoidance of the ambiguity in registration, and the combination of the multiplicative extended Kalman filter and the pose graph optimization to reduce the effect of measurement noise and drift error. A hardware experiment was performed to generate the sequence of point clouds of a three-axis free-floating target and validate our method. The result shows that the proposed method outperforms existing methods and effectively identifies the inertial parameters, including the normalized principal moments of inertia and the orientation of principal axes.     

Pose measurement and motion estimation of non-cooperative satellite based on spatial circle feature

Malfunctioned satellites have seriously threatened orbital safety, and the capture of these satellites is of great significance. The pose measurement and the motion estimation of the tumbling satellite is the premise of capture. In this paper, the docking ring of the satellite is identified, which is equivalent to a spatial circle. Combined with the nozzle feature, the pose duality of the spatial circle can be eliminated. And the measurement accuracy is improved by minimizing the reprojection error of the docking ring and the nozzle. Due to the symmetry of the docking ring, the measured pose has only five degrees of freedom, losing the degree of freedom of rotation around the normal vector. In the motion estimation algorithm, the observability of the tumbling motion is firstly analyzed, then an error-state Kalman filter with inertia ratio constraints is designed. To improve the convergence speed and stability of the filter, a rough estimation algorithm of filter initial value based on linear term extraction and particle swarm optimization is proposed. The effectiveness of the pose measurement and motion estimation method is verified by simulations.     

Early analysis of precise orbit and clock offset determination for the satellites of the global BeiDou-3 system

Eight new-generation BeiDou satellites (BeiDou-3) have been launched into Medium Earth Orbit (MEO), allowing for global coverage since March 2018, and they are equipped with new hydrogen atomic clocks and updated rubidium clocks. Firstly, we analyzed the signals for the carrier-to-noise-density ratio (C/N0) and pseudorange multipath (MP) by using international GNSS (Global Navigation Satellite System) Monitoring and Assessment System (iGMAS) station data, and found that B1C has a lower C/N0, and B2a has the same level of C/N0 as the B1I and B3I signals. For pseudorange multipath, compared with the BeiDou-2 satellites, the obvious systematic variation of MP scatters related to the elevation angle is greatly improved for the BeiDou-3 and BeiDou-3e satellites signals. For the signals of the BeiDou-3 satellites, the order of the Root Mean Square (RMS) values of multipath and noise is B3I?<?B1I?<?B2a?<?B1C. Then, the comparison of the precise orbit determination and clock offset determination for the BeiDou-2, BeiDou-3, and BeiDou-3 experimental (BeiDou-3e) satellites was done by using 10 stations from iGMAS. The 3D precision of the 24?h orbit overlap is 24.55, 25.61, and 23.35?cm for the BeiDou-3, BeiDou-3e, and BeiDou-2 satellites, respectively. BeiDou-3 satellite has a comparable precision to that of the BeiDou-2 satellite. For the precision of clock offset estimation, the Standard Deviation (STD) of the BeiDou-3 MEO satellite is 0.350?ns, which is an improvement of 0.042?ns over that of the BeiDou-2 MEO satellite. The stabilities of the BeiDou-3 and BeiDou-3e onboard clocks are better than those of BeiDou-2 by factors of 2.84 and 1.61 at an averaging time of 1000 and 10,000?s, respectively.     

Near real-time estimation of tropospheric water vapour content from ground based GNSS data and its potential contribution to weather now-casting in Austria

The importance of high resolution meteorological analysis of the atmosphere increased over the past years. A detailed analysis of the humidity field is an important precondition for a better monitoring of local and regional extreme precipitation events and for forecasts with improved spatial resolution. For this reason, the Austrian Meteorological Agency (ZAMG) is operating the spatial and temporal high resolution INCA system (Integrated Now-casting through Comprehensive Analysis) since begin of 2005. Errors in this analysis occur mainly in the areas of rapidly changing and hard to predict weather conditions or rugged topography with extreme differences in height such as the alpine area of Austria. The aim of this work is to provide GNSS based measurements of the tropospheric water vapour content with a temporal resolution of 1 h and a temporal delay of less than 1 h to assimilate these estimates into the INCA system. Additional requirement is an accuracy of better than 1 mm of the precipitable water (PW) estimates.     

Comparison of RO tropopause height based on different tropopause determination methods

The tropopause is an important boundary in the Earth’s atmosphere, and has been the subject of close attention from atmosphere and climate researchers. To monitor the global tropopause using radio occultation (RO) data, there are two primary methods, one is the widely used temperature lapse rate method, and the other is the bending angle covariance transform method which is unique to RO data. We use FengYun3-C (FY3C) and Meteorological Operational Satellite Program (MetOp) RO data and European Centre for Medium-Range Weather Forecasts (ECMWF) operational analysis data to determine differences in RO tropopause height calculated by these two methods. We compute biases of the RO lapse rate tropopause height (LRTH) and the RO bending angle tropopause height (BATH) relative to the ECMWF LRTH. The dependences of the tropopause height biases on tropopause height (TPH) retrieval method, latitude, season and RO mission are investigated. The results indicate that BATH show a consistent 0.8–1.2 km positive bias over the tropics and high latitude regions compared with LRTH, however, over 25° to 40° latitude in both hemisphere, BATH results are less stable. Furthermore, the mean bias between BATH and LRTH displayed a different symmetrical characteristic from 2017.12 to 2018.2 (DJF) compared to 2018.6–2018.8 (JJA). However, except for some bias over Antarctica, the mean value of both LRTH and BATH show a similar tropopause variation, indicating the consistency of both methods.     

Ionospheric behavior of the F2 peak parameters foF2 and hmF2 at Hainan and comparisons with IRI model predictions

Monthly median values of foF2, hmF2 and M(3000)F2 parameters, with hourly time interval resolution for the diurnal variation, obtained with DPS-4 digisonde observations at Hainan (19.4°N, 109.0°E) are used to study the low latitude ionospheric variation behavior. The observational results are compared with the International Reference Ionospheric Model (IRI) predictions. The time period coverage of the data used for the present study is from March 2002 to February 2005. Our present study showed that: (1) In general, IRI predictions using CCIR and URSI coefficients follow well the diurnal and seasonal variation patterns of the experimental values of foF2. However, CCIR foF2 and URSI foF2 IRI predictions systematically underestimate the observed results during most time period of the day, with the percentage difference ΔfoF2 (%) values changing between about −5% and −25%, whereas for a few hours around pre-sunrise, the IRI predictions generally overestimate the observational ones with ΔfoF2 (%) sometimes reaching as large as ∼30%. The agreement between the IRI results and the observational ones is better for the year 2002 than for the other years. The best agreement between the IRI results and the observational ones is obtained in summer when using URSI coefficients, with the seasonal average values of ΔfoF2 (%) being within the limits of ±10%. (2) In general, the IRI predicted hmF2 values using CCIR M(3000)F2 option shows a poor agreement with the observational results. However, when using the measured M(3000)F2 as input, the diurnal variation pattern of hmF2 given by IRI2001 has a much better agreement with the observational one with the detailed fine structures including the pre-sunrise and post-sunset peaks reproduced reasonably well. The agreement between the IRI predicted hmF2 values using CCIR M(30,000)F2 option and the observational ones is worst for the afternoon to post-midnight hours for the high solar activity year 2002. During daytime hours the agreement between the hmF2 values obtained with CCIR M(30,000)F2 option and the observational ones is best for summer season. The discrepancy between the observational hmF2 and that obtained with CCIR M(30,000)F2 option stem from the CCIR M(3000)F2 model, which does not produce the small scale structures observed in the measured M(3000)F2.