Generating precise and homogeneous orbits for Jason-1 and Jason-2

Driven by the GMES (Global Monitoring for Environment and Security) and GGOS (Global Geodetic Observing System) initiatives the user community has a strong demand for high-quality altimetry products. In order to derive such high-quality altimetry products, precise orbits for the altimetry satellites are a necessity. With the launch of the TOPEX/Poseidon mission in 1992 a still on-going time series of high-accuracy altimetry measurements of ocean topography started, continued by the altimetry missions Jason-1 in 2001 and Jason-2/OSTM in 2008. This paper contributes to the on-going orbit reprocessing carried out by several groups and presents the efforts of the Navigation Support Office at ESA/ESOC using its NAPEOS software for the generation of precise and homogeneous orbits referring to the same reference frame for the altimetry satellites Jason-1 and Jason-2. Data of all three tracking instruments on-board the satellites (beside the altimeter), i.e. GPS, DORIS, and SLR measurements, were used in a combined data analysis. About 7 years of Jason-1 data and more than 1 year of Jason-2 data were processed. Our processing strategy is close to the GDR-C standards. However, we estimated slightly different scaling factors for the solar radiation pressure model of 0.96 and 0.98 for Jason-1 and Jason-2, respectively. We used 30 s sampled GPS data and introduced 30 s satellite clocks stemming from ESOC’s reprocessing of the combined GPS/GLONASS IGS solution. We present the orbit determination results, focusing on the benefits of adding GPS data to the solution. The fully combined solution was found to give the best orbit results. We reach a post-fit RMS of the GPS phase observation residuals of 6 mm for Jason-1 and 7 mm for Jason-2. The DORIS post-fit residuals clearly benefit from using GPS data in addition, as the DORIS data editing improves. The DORIS observation RMS for the fully combined solution is with 3.5 mm and 3.4 mm, respectively, 0.3 mm better than for the DORIS-SLR solution. Our orbit solution agrees well with external solutions from other analysis centers, as CNES, LCA, and JPL. The orbit differences between our fully combined orbits and the CNES GDR-C orbits are of about 0.8 cm for Jason-1 and at 0.9 cm for Jason-2 in the radial direction. In the cross-track component we observe a clear improvement when adding GPS data to the POD process. The 3D-RMS of the orbit differences reveals a good orbit consistency at 2.7 cm and 2.9 cm for Jason-1 and Jason-2. Our resulting orbit series for both Jason satellites refer to the ITRF2005 reference frame and are provided in sp3 file format on our ftp server.     

Impact of Jason-2/T2L2 Ultra-Stable-Oscillator Frequency Model on DORIS stations coordinates and Earth Orientation Parameters

Doppler Orbitography and Radiopositioning Integrated by Satellite (DORIS) is a tracking technique based on a one-way ground to space Doppler link. For Low Earth Orbit (LEO) satellites, DORIS shows a robust capability in terms of data coverage and availability, due to a wide and well-distributed ground network, where data are made available by the International Doris Service (IDS). However, systematic errors remain in the DORIS data, such as instabilities of the on-board clock due to radiation encountered in space, which limit the accurate determination of station positions.The DORIS on-board clock frequency stability is degraded by the increased radiation found in the region of the South Atlantic Anomaly (SAA) and has been shown to degrade station position estimation. This paper introduces a new model correction to the DORIS data for the frequency of the Jason-2 Ultra Stable Oscillator (USO), derived from the Time Transfer by Laser Link (T2L2) experiment (Belli and Exertier, 2018). We show that a multi-satellite DORIS solution including this T2L2-corrected data applied to the frequency modelling for The DORIS data, improves the estimation of station coordinates. We show the tie residuals with respect to collocated GPS stations are improved by several millimeters. We also demonstrate that the 117-day (Jason-2) draconitic signal in the geophysical parameters is reduced, implying that the origin of this signal is not just solar radiation pressure mis-modeling, but also radiation-induced clock perturbations on the Jason-2 DORIS Ultra-Stable-Oscillator (USO). Finally we demonstrate through comparisons with the International Earth Rotations and Reference Systems Service (IERS) C04 series for Earth Orientation Parameters (EOP), that the estimation of EOP is improved in both a Jason-2 DORIS-only and a multi-satellite DORIS solution for EOP.     

Jason-2 DORIS phase measurement processing

The DORIS instrument on Jason-2 is the first of a new generation. The satellite receivers have now seven simultaneous measurement channels, with synchronous dual frequency phase and pseudo-range measurements. These measurements are now described in a similar manner as GPS measurements and an extension of the RINEX 3.0 format has been defined for DORIS. Data are available to users with a shorter latency.     

DPOD2005: An extension of ITRF2005 for Precise Orbit Determination

For Precise Orbit Determination of altimetry missions, we have computed a data set of DORIS station coordinates defined for specific time intervals called DPOD2005. This terrestrial reference set is an extension of ITRF2005. However, it includes all new DORIS stations and is more reliable, as we disregard stations with large velocity formal errors as they could contaminate POD computations in the near future. About 1/4 of the station coordinates need to be defined as they do not appear in the original ITRF2005 realization. These results were verified with available DORIS and GPS results, as the integrity of DPOD2005 is almost as critical as its accuracy. Besides station coordinates and velocities, we also provide additional information such as periods for which DORIS data should be disregarded for specific DORIS stations, and epochs of coordinate and velocity discontinuities (related to either geophysical events, equipment problem or human intervention). The DPOD model was tested for orbit determination for TOPEX/Poseidon (T/P), Jason-1 and Jason-2. Test results show DPOD2005 offers improvement over the original ITRF2005, improvement that rapidly and significantly increases after 2005. Improvement is also significant for the early T/P cycles indicating improved station velocities in the DPOD2005 model and a more complete station set. Following 2005 the radial accuracy and centering of the ITRF2005-original orbits rapidly degrades due to station loss.     

Towards development of a consistent orbit series for TOPEX,Jason-1, and Jason-2

The TOPEX/Poseidon, Jason-1 and Jason-2 set of altimeter data now provide a time series of synoptic observations of the ocean that span nearly 17 years from the launch of TOPEX in 1992. The analysis of the altimeter data including the use of altimetry to monitor the global change in mean sea level requires a stable, accurate, and consistent orbit reference over the entire time span. In this paper, we describe the recomputation of a time series of orbits that rely on a consistent set of reference frames and geophysical models. The recomputed orbits adhere to the IERS 2003 standards for ocean and earth tides, use updates to the ITRF2005 reference frame for both the SLR and DORIS stations, apply GRACE-derived models for modeling of the static and time-variable gravity, implement the University College London (UCL) radiation pressure model for Jason-1, use improved troposphere modeling for the DORIS data, and apply the GOT4.7 ocean tide model for both dynamical ocean tide modeling and for ocean loading. The new TOPEX orbits have a mean SLR fit of 1.79 cm compared to 2.21 cm for the MGDR-B orbits. These new TOPEX orbits agree radially with independent SLR/crossover orbits at 0.70 cm RMS, and the orbit accuracy is estimated at 1.5–2.0 cm RMS over the entire TOPEX time series. The recomputed Jason-1 orbits agree radially with the Jason-1 GDR-C orbits at 1.08 cm RMS. The GSFC SLR/DORIS dynamic and reduced-dynamic orbits for Jason-2 agree radially with independent orbits from the CNES and JPL at 0.70–1.06 cm RMS. Applying these new orbits, and using the latest altimeter corrections for TOPEX, Jason-1, and Jason-2 from September 1992 to May 2009, we find a global rate in mean sea level of 3.0 ± 0.4 mm/yr.     

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.     

The X-TRACK/ALES multi-mission processing system: New advances in altimetry towards the coast

In the context of the ESA Climate Change Initiative project, a new coastal sea level altimetry product has been developed in order to support advances in coastal sea level variability studies. Measurements from Jason-1,2&3 missions have been retracked with the Adaptive Leading Edge Subwaveform (ALES) Retracker and then ingested in the X-TRACK software with the best possible set of altimetry corrections. These two coastal altimetry processing approaches, previously successfully validated and applied to coastal sea level research, are combined here for the first time in order to derive a 16-year-long (June 2002 to May 2018), high-resolution (20-Hz), along-track sea level dataset in six regions: Northeast Atlantic, Mediterranean Sea, West Africa, North Indian Ocean, Southeast Asia and Australia. The study demonstrates that this new coastal sea level product called X-TRACK/ALES is able to extend the spatial coverage of sea level altimetry data ~3.5 km in the land direction, when compared to the X-TRACK 1-Hz dataset. We also observe a large improvement in coastal sea level data availability from Jason-1 to Jason-3, with data at 3.6 km, 1.9 km and 0.9 km to the coast on average, for Jason-1, Jason-2 and Jason-3, respectively. When combining measurements from Jason-1 to Jason-3, we reach a distance of 1.2–4 km to the coast. When compared to tide gauge data, the accuracy of the new altimetry near-shore sea level estimations also improves. In terms of correlations with a large set of independent tide gauge observations selected in the six regions, we obtain an average value of 0.77. We also show that it is now possible to derive from the X-TRACK/ALES product an estimation of the ocean current variability up to 5 km to the coast. This new altimetry dataset, freely available, will provide a valuable contribution of altimetry in coastal marine research community.     

Impact of DORIS ground antennas environment on their radio signal quality

Among the factors which may disrupt the DORIS measurements quality, the ground antennas environment is of high importance. For a set of 15 selected DORIS beacon, the differences between the effective and theoretical power received on-board the satellites (SPOT-5 and Envisat) have been analyzed in terms of spatial direction around the antenna. Such antenna maps have also been established regarding the Doppler residuals of the least-square precise orbit adjustment. Thanks to 360° views from the antennas and aerial views of the sites, the impact of the signal obstructions (trees, roofs, antennas …) on power attenuation and Doppler residuals is discussed. Depending on the nature of the obstructed object, the attenuation level can reach more than 5 dB, and the residual RMS of the orbit adjustment may be doubled from the nominal value, reaching 1 mm/s locally. The nature of the ground at the foot of the antennas has been correlated to DORIS signal quality at high elevation: reflections on flat surfaces (e.g. roofs) affect the signal more significantly than reflections on natural ground (e.g. soil). In particular, a modeling of the multipath phenomenon affecting Fairbanks site has been established and fits remarkably with the observations. Finally, an evaluation of the direct impact of obstructing objects on the orbit has also been performed. The example of a scaffolding at Kauai site displays a few millimeters error in the along-track position of the satellite.     

Status of the T2L2/Jason2 Experiment

The T2L2 (Time Transfer by Laser Link) project, developed by CNES and OCA will permit the synchronization of remote ultra stable clocks and the determination of their performances over intercontinental distances. The principle of the experiment derives from Satellite Laser Ranging (SLR) technology with dedicated space equipment. T2L2 was accepted in 2005 to be on board the Jason2 altimetry satellite. The payload consists of both event timer and photo detection modules. The system uses the ultra-stable quartz oscillator of DORIS as on-board reference clock on one hand, and the Laser Reflector Array, making T2L2 a real two-way time transfer system on the other hand. The expected time stability of the T2L2 instrument (detection and timing), referenced by the DORIS oscillator and including all internal error sources should be at the level of 10–12 ps at 1 s and <1 ps at 1000 s. The metrological specifications of T2L2 should permit to maintain a precision of 1 to a few ps when measuring the phase of a clock during around 1000 seconds.     

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.     

嫦娥一号卫星飞行控制模拟器的设计与研制

飞行控制模拟器,简称飞控模拟器,是飞控中心用来验证飞行程序和地面测控方案,训练操作人员的模拟仿真系统,为系统演练构造一个尽量逼真的任务环境.介绍了嫦娥一号飞控模拟器的功能、结构和使用模式.提出了带轨道根数跳时的新方法,使模拟器具备加速跳时至任意轨道段的能力,提高了任务演练的灵活性.建立了高精度轨道动力学模型,并首次在模拟器中直接嵌入了部分星载软件,大大提高了模拟器的保真度.同时专门设计了图形化的集中监控软件,方便了用户的操作和管理.该模拟器在嫦娥一号卫星任务准备中得到了充分的应用,为嫦娥一号卫星任务的圆满完成奠定了坚实的基础.根据该模拟器的开发经验,对未来模拟器的研制提出了几点建议.     

基于VxWorks的小天体撞击任务的星载GNC软件设计

针对小天体撞击任务,应用VxWorks嵌入式实时操作系统,设计小天体高速撞击器的星载GNC软件部分．对星载GNC系统结构进行简要描述;在此基础上,综合考虑小天体撞击任务的实时性要求、不同飞行模式的耦合关系、轨道确定的数学运算量以及对不同敏感器数据采集的周期性控制等多方面因素,对小天体撞击任务进行模块化分解,提出各个任务模块间的同步方式与通信手段;在PC-104嵌入式计算机与dSPACE实时仿真平台的联合环境下,对所设计的星载GNC软件进行仿真验证,结果表明,基于VxWorks嵌入式实时操作系统所设计的小天体撞击GNC软件完全可以满足小天体撞击任务的实时性要求,为撞击任务的顺利进行提供有效的保证．     

Calibration of Envisat radar altimeter over Lake Issykkul

This study presents the results of calibration/validation (C/V) of Envisat satellite radar altimeter over Lake Issykkul located in Kyrgyzstan, which was chosen as a dedicated radar altimetry C/V site in 2004. The objectives are to estimate the absolute altimeter bias of Envisat and its orbit based on cross-over analysis with TOPEX/Poseidon (T/P), Jason-1 and Jason-2 over the ocean. We have used a new method of GPS data processing in a kinematic mode, developed at the Groupe de Recherche de Geodesie Spatiale (GRGS), which allows us to calculate the position of the GPS antenna without needing a GPS reference station. The C/V is conducted using various equipments: a local GPS network, a moving GPS antenna along the satellites tracks over Lake Issykkul, In Situ level gauges and weather stations. The absolute bias obtained for Envisat from field campaigns conducted in 2009 and 2010 is between 62.1 and 63.4 ± 3.7 cm, using the Ice-1 retracking algorithm, and between 46.9 and 51.2 cm with the ocean retracking algorithm. These results differ by about 10 cm from previous studies, principally due to improvement of the C/V procedure. Apart from the new algorithm for GPS data processing and the orbit error reduction, more attention has been paid to the GPS antenna height calculation, and we have reduced the errors induced by seiche over Lake Issykkul. This has been assured using cruise data along the Envisat satellite track at the exact date of the pass of the satellite for the two campaigns. The calculation of the Envisat radar altimeter bias with respect to the GPS levelling is essential to allow the continuity of multi-mission data on the same orbit, with the expected launch of SARAL/Altika mission in 2012. Implications for hydrology in particular, will be to produce long term homogeneous and reliable time series of lake levels worldwide.     

GNSS-R concept extended by a fine orbit tuning

Small changes in semimajor axis of the orbits selected for the GNSS-R [R as Reflectometry] satellites, so-called fine orbit tuning, known from the ESA’s Gravity and steady-state Ocean Circulation Explorer mission, can dramatically increase the number of nadir and off-nadir reflecting points and, in turn, can enhance the capability of the concept of bistatic altimetry (GNSS Reflectometry) without additional costs. The application of our suggestion is feasible for a satellite which will be equipped by thrusters for the orbit keeping. During the mission lifetime several orbit tunings are feasible, just to transfer from one to another orbit. Then we can study short-periodic or longer-periodic features, according to scientific goals defined for the mission. The shortest cycles (few days), corresponding to the required revisit time (defined by ESA), may be subcycles of much longer cycles (repeat periods).     

DPOD2014: A new DORIS extension of ITRF2014 for precise orbit determination

To support precise orbit determination of the altimetry missions, the International DORIS Service (IDS) regularly estimates the DPOD (DORIS terrestrial reference frame for Precise Orbit Determination) solution which includes mean positions and velocities of all the DORIS stations. This solution is aligned to the current realization of the International Terrestrial Reference Frame (ITRF) and so, can be seen as a DORIS extension of the ITRF. In 2016, moving to the IDS Combination Center, the DPOD construction scheme changed. The new DPOD solution is produced from a DORIS cumulative position and velocity solution. We present the new methodology used to compute DPOD2014 and its validation procedure. In order to present geophysical applications and interpretations of these results, we show two examples: (1) the Gorkha earthquake (M7.8 – April 2015) generates a 3-D mis-positioning of nearly 55?mm of the EVEB DORIS station at the Everest base camp 90?km from the epicenter. (2) Applying the results the DPOD2014 realization, we show that the most recent vertical velocity of Thule, Greenland is similar to that observed between 2006 and 2010, indicating further ongoing ice mass loss in the Thule region of northwest Greenland.     

Derniers resultats de calcul d'orbite SPOT et TOPEX/POSEIDON obtenus avec DORIS

DORIS is a tracking system developed by CNES to support precise orbit computation of the US-French TOPEX/POSEIDON project. Moreover DORIS data are currently processed to compute SPOT2 and SPOT3 orbits. Although the SPOT satellites are at 800 km altitudes, their orbits reach the decimetric level, to be compared to the 3–4 cm RMS on the T/P radial component.For each type of orbit, there is an adaptive period of a few months which is used for improving the precision. The paper describes what has been done to reduce some items in the total error budget. The latest results will be presented as well as the criteria which are settled to characterize the improvements.Finally, the future developments of DORIS and the potential projects flowing DORIS will be outlined such as ENVISAT and TOPEX/POSEIDON Follow On.     

The Huygens Probe and mission design

The mission of the Huygens Probe is to enter into the Titan atmosphere and to descend in a controlled manner to the surface such that scientific measurements can be taken and transmitted to earth via the Orbiter. The entry is the most challenging phase of the mission. The deceleration is achieved by an aeroshell. Thermal protection against entry heat fluxes is achieved by an ablator rejecting heat mainly by radiation. The descent profile in terms of descent velocity, spin rate and attitude stability is designed to allow imaging, sampling and measurements at the various atmospheric layers within the overall time and energy constraints. The Huygens system has to operate in an autonomous and fault tolerant manner, as no ground control and failure recovery during the descent phase is feasible. This autonomy and fault tolerance is the main driver for the chosen architecture of the on-board system.     

DORIS and GPS monitoring of the Gavdos calibration site in Crete

Due to its specific geographical location as well as its geodetic equipment (DORIS, GNSS, microwave transponder and tide gauges), the Gavdos station in Crete, Greece is one of the very few sites around the world used for satellite altimetry calibration. To investigate the quality of the Gavdos geodetic coordinates and velocities, we analyzed and compared here DORIS and GPS-derived results obtained during several years of observations. The DORIS solution is the latest ignwd11 solution at IGN, expressed in ITRF2008, while the GPS solution was obtained using the GAMIT software package. Current results show that 1–2 mm/yr agreement can be obtained for 3-D velocity, showing a good agreement with current geophysical models. In particular, the agreement obtained for the vertical velocity is around 0.3–0.4 mm/yr, depending on the terrestrial reference frame. As a by-product of these geodetic GPS and DORIS results, Zenith Tropospheric Delays (ZTDs) estimations were also compared in 2010 between these two techniques, and compared to ECMWF values, showing a 6.6 mm agreement in dispersion without any significant difference between GPS and DORIS (with a 97.6% correlation), but with a 13–14 mm agreement in dispersion when comparing to ECMWF model (with only about 90% correlation for both techniques). These tropospheric delay estimations could also provide an external calibration of the tropospheric correction used for the geophysical data of satellite altimetry missions.