An inter-comparison of zenith tropospheric delays derived from DORIS and GPS data

Doppler Orbitography Radiopositioning Integrated by Satellite (DORIS) and Global Positioning System (GPS) techniques are similarly affected by propagation delays in the neutral atmosphere (troposphere) and hence make use of similar data processing strategies for reducing this effect. We compare Zenith Tropospheric Delays (ZTDs) estimated from 52 DORIS and GPS station pairs co-located at 35 sites over the 2005–2008 period. We find an overall systematic negative mean bias of −4 mm and a median bias of −2 mm, with a large site-to-site scatter and especially stronger biases over South America, potentially linked to remaining problems related to the South Atlantic Anomaly (SAA) in the current DORIS data processing. The standard deviation of ZTD differences is in the range 4–12 mm over the globe (8 mm on average), with larger values located in the southern hemisphere. The spatial variability of differences is consistent with previous work but remains largely unexplained. DORIS is shown to be much less sensitive to instrumental changes than GPS (only the switch from Alcatel to Starec antenna at Toulouse is detected as an offset of −4 mm in the ZTD time series). On the opposite, discontinuities and spurious annual signals are found in the GPS ZTD solutions. A discontinuity of +5 mm is found on 5 November 2006, linked to the switch from relative to absolute GPS antenna models used in the data processing. The use of modified GPS antennas (e.g. at GODE) or improved antenna models is shown to reduce the spurious annual signal (e.g. from 5 mm to 2 mm at METS). Overall, the agreement between both techniques is good, though DORIS shows a significantly larger random scatter. The high stability and good spatial and temporal coverage make DORIS a potential candidate technique for meteorology and climate studies as long as reasonable time averaging can be applied (e.g. differences are reduced from 8.6 to 2.4 mm with 5-day averages) and no real-time application is considered. This technique could be considered as a potential contributor to Global Geodetic Observing System (GGOS) for climatology.     

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.     

DORIS system: The new age

The boarding of the first DGXX DORIS instrument on Jason-2 mission gives us the opportunity to present the improvements that have been implemented on the DORIS system. The goal of this paper is to present information about the new capacities of the DORIS system and to give the current status of its components. An overview of the DORIS system, the International DORIS Service and the Jason-2 satellite mission are first presented. Then the new characteristics of the on-board instrument are detailed. The capacity to track up to seven ground beacons simultaneously dramatically increases the number of measurements performed: a factor of three increase over Jason-1 is observed at the altitude of 1330 km. It also increases the diversity of directions of observation and allows low elevation measurements from 0°. The new phase measurements capability allows now phase processing. The instability of the Jason-1 USOs (Ultra-Stable Oven-controlled quartz oscillator) while crossing the South Atlantic Anomaly has been solved by decreasing the sensitivity to radiation by a factor of 10. New features of the on-board software enhance the coastal and inland water altimetry and increase the robustness of the data. The new software also improves the real time orbit accuracy for operational altimetry. The improvements introduced concurrently on the ground segment have also significantly enhanced capability. The new RINEX exchange formats provide simultaneous phase and pseudo-range measurements. The maintenance of the DORIS Beacons Network and the work done by the DORIS Signal Integrity monitoring team lead to an increased availability of the Network from 75% to 90% and so to a more homogenous orbit coverage.     

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.     

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.     

GSFC DORIS contribution to ITRF2008

The NASA GSFC DORIS analysis center has provided weekly DORIS solutions from November 1992 to January 2009 (839 SINEX files) of station positions and Earth Orientation Parameters for inclusion in the DORIS contribution to ITRF2008. The NASA GSFC GEODYN orbit determination software was used to process the orbits and produce the normal equations. The weekly SINEX gscwd10 submissions included DORIS data from Envisat, TOPEX/Poseidon, SPOT-2, SPOT-3, SPOT-4, SPOT-5. The orbits were mostly seven days in length (except for weeks with data gaps or maneuvers). The processing used the GRACE-derived EIGEN-GL04S1 gravity model, updated modeling for time-variable gravity, the GOT4.7 ocean tide model and tuned satellite-specific macromodels for SPOT-2, SPOT-3, SPOT-4, SPOT-5 and TOPEX/Poseidon. The University College London (UCL) radiation pressure model for Envisat improves nonconservative force modeling for this satellite, reducing the median residual empirical daily along-track accelerations from 3.75 × 10⁻⁹ m/s² with the a priori macromodel to 0.99 × 10⁻⁹ m/s² with the UCL model. For the SPOT and Envisat DORIS satellite orbits from 2003 to 2008, we obtain average RMS overlaps of 0.8–0.9 cm in the radial direction, 2.1–3.4 cm cross-track, and 1.7–2.3 cm along-track. The RMS orbit differences between Envisat DORIS-only and SLR & DORIS orbits are 1.1 cm radially, 6.4 cm along-track and 3.7 cm cross-track and are characterized by systematic along-track mean offsets due to the Envisat DORIS system time bias of ±5–10 μs. We obtain a good agreement between the geometrically-determined geocenter parameters and geocenter parameters determined dynamically from analysis of the degree one terms of the geopotential. The intrinsic RMS weekly position repeatability with respect to the IDS-3 combination ranges from 2.5 to 3.0 cm in 1993–1994 to 1.5 cm in 2007–2008.     

Low-latitude ionospheric scintillations and total electron content obtained with the CITRIS instrument on STPSat1 using radio transmissions from DORIS ground beacons

The primary objective of the Scintillation and Tomography Receiver in Space (CITRIS) is to detect ionospheric irregularities from space at low latitude. For this purpose, the satellite receiver uses the UHF and S-Band transmissions of the ground network of Doppler Orbitography and Radiopositioning Integrated by Satellite (DORIS) beacons. CITRIS, developed at the Naval Research Laboratory, differs from the normal DORIS receiver by being able to capture and store the complex amplitude of the 401.25 and 2036.25 MHz transmissions at 200 Hz sample rate. Ground processing of the CITRIS data yields total electron content (TEC) and both phase and amplitude scintillations. With CITRIS flying on the US Space Test Program (STP) satellite STPSat1, 2 years of data were collected and processed to determine the fluctuations in ionospheric TEC and radio scintillations associated with equatorial irregularities. CITRIS flights over DORIS transmitters yield direct measurements of the horizontal plasma density fluctuations associated with equatorial plasma bubbles. Future flights of CITRIS can provide valuable complements to other satellite instruments such as GPS occultation receivers used to estimate vertical electron density profiles in the ionosphere.     

The GLONASS-M satellite yaw-attitude model

The proper modeling of the satellites’ yaw-attitude is a prerequisite for high-precision Global Navigation Satellite System (GNSS) positioning and poses a particular challenge during periods when the satellite orbital planes are partially eclipsed. Whereas a lot of effort has been put in to examine the yaw-attitude control of GPS satellites that are in eclipsing orbits, hardly anything is known about the yaw-attitude behavior of eclipsing GLONASS-M satellites. However, systematic variations of the carrier phase observation residuals in the vicinity of the orbit’s noon and midnight points of up to ±27 cm indicate significant attitude-related modeling issues. In order to explore the GLONASS-M attitude laws during eclipse seasons, we have studied the evolution of the horizontal satellite antenna offset estimates during orbit noon and orbit midnight using a technique that we refer to as “reverse kinematic precise point positioning”. In this approach, we keep all relevant global geodetic parameters fixed and estimate the satellite clock and antenna phase center positions epoch-by-epoch using 30-second observation and clock data from a global multi-GNSS ground station network. The estimated horizontal antenna phase center offsets implicitly provide the spacecraft’s yaw-attitude. The insights gained from studying the yaw angle behavior have led to the development of the very first yaw-attitude model for eclipsing GLONASS-M satellites. The derived yaw-attitude model proves to be much better than the nominal yaw-attitude model commonly being used by today’s GLONASS-capable GNSS software packages as it reduces the observation residuals of eclipsing satellites down to the normal level of non-eclipsing satellites and thereby prevents a multitude of measurements from being incorrectly identified as outliers. It facilitates continuous satellite clock estimation during eclipse and improves in particular the results of kinematic precise point positioning of ground-based receivers.     

DORIS processing at the European Space Operations Centre

This paper gives an overview of the DORIS related activities at the Navigation Support Office of the European Space Operations Centre. The DORIS activities were started in 2002 because of the launch of the Envisat satellite where ESOC is responsible for the validation of the Envisat Precise Orbits and a brief overview of the key Envisat activities at ESOC is given. Typical orbit comparison RMS values between the CNES POE (GDR-C) and the ESOC POD solution is 6.5, 18.8 and 23.1 mm in radial-, along- and cross-track direction. In the framework of the generation of the ITRF2008 ESOC participated in the reprocessing of all three space geodetic techniques; DORIS, SLR, and GPS. Here the main results of our DORIS reprocessing, in the framework of the International DORIS Service (IDS), are given. The WRMS of the weekly ESOC solution (esawd03) for the 2004–2009 period compared to the IDS-1 combined solution is of the order of 12 mm. Based on the long time series of homogeneously processed data a closer look is taken at the estimated solar radiation pressure parameters of the different satellites used in this DORIS analysis. The main aim being the stabilization of the Z-component of the geocentre estimates. We conclude that the ESOC participation to the IDS ITRF2008 contribution has been beneficial for both ESOC and the IDS. ESOC has profited significantly from the very open and direct communications and comparisons that took place within the IDS during the reprocessing campaign.     

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.     

DORIS geodesy: A dynamic determination of geocentre location

Geoscience Australia contributed a multi-satellite, multi-year weekly time series to the International DORIS Service combined submission for the construction of International Terrestrial Reference Frame 2008 (ITRF2008). This contributing solution was extended to a study of the capability of DORIS to dynamically estimate the variation in the geocentre location. Two solutions, comprising different constraint configurations of the tracking network, were undertaken. The respective DORIS satellite orbit solutions (SPOT-2, SPOT-4, SPOT-5 and Envisat) were verified and validated by comparison with those produced at the Goddard Space Flight Center (GSFC), DORIS Analysis Centre, for computational consistency and standards. In addition, in the case of Envisat, the trajectories from the GA determined SLR and DORIS orbits were compared. The results for weekly dynamic geocentre estimates from the two constraint configurations were benchmarked against the geometric geocentre estimates from the IDS-2 combined solution. This established that DORIS is capable of determining the dynamic geocentre variation by estimating the degree one spherical harmonic coefficients of the Earth’s gravity potential. It was established that constrained configurations produced similar results for the geocentre location and consequently similar annual amplitudes. For the minimally constrained configuration Greenbelt–Kitab, the mean of the uncertainties of the geocentre location were 2.3, 2.3 and 7.6 mm and RMS of the mean uncertainties were 1.9, 1.2 and 3.5 mm for the X, Y and Z components, respectively. For GA_IDS-2_Datum constrained configuration, the mean of the uncertainties of the geocentre location were 1.7, 1.7 and 6.2 mm and RMS of the mean uncertainties were 0.9, 0.7 and 2.9 mm for the X, Y and Z components, respectively. The mean of the differences of the two DORIS dynamic geocentre solutions with respect to the IDS-2 combination were 1.6, 4.0 and 5.1 mm with an RMS of the mean 21.2, 14.0 and 31.5 mm for the Greenbelt–Kitab configuration and 4.1, 3.9 and 4.3 mm with an RMS 8.1, 9.0 and 28.6 mm for the GA_IDS-2_Datum constraint configuration. The annual amplitudes for each component were estimated to be 5.3, 10.8 and 11.0 mm for the Greenbelt–Kitab configuration and 5.3, 9.3 and 9.4 mm for the GA_IDS-2_Datum constraint configuration. The two DORIS determined dynamic geocentre solutions were compared to the SLR determined dynamic solution (which was determined from the same process of the GA contribution to the ITRF2008 ILRS combination) gave mean differences of 3.3, −4.7 and 2.5 mm with an RMS of 20.7, 17.5 and 28.0 mm for the X, Y and Z components, respectively for the Greenbelt–Kitab configuration and 1.1, −5.4 and 4.4 mm with an RMS of 9.7, 13.3 and 24.9 mm for the GA_IDS-2_Datum configuration. The larger variability is reflected in the respective amplitudes. As a comparison, the annual amplitudes of the SLR determined dynamic geocentre are 0.9, 1.0 and 6.8 mm in the X, Y and Z components. The results from this study indicate that there is potential to achieve precise dynamically determined geocentre from DORIS.     

Demonstration of the use of the Doppler Orbitography and Radio positioning Integrated by Satellite (DORIS) measurements to validate GPS ionospheric imaging

There is a lack of independent ionospheric data that can be used to validate GPS imaging results at mid latitudes over severe storm times. Doppler Orbitography and Radio positioning Integrated by Satellite (DORIS), a global network of dual-frequency ground to satellite observations, provides this missing data and here is employed as verification to show the accuracy of the ionospheric GPS images in terms of the total electron content (TEC). In this paper, the large-scale ionospheric structures that appeared during the strong geomagnetic storm of 20 November 2003 are reconstructed with a GPS tomographic algorithm, known as MIDAS, and validated with DORIS TEC measurements. The main trough shown in an extreme equatorward position in the ionospheric imaging over mainland Europe is confirmed by DORIS satellite measurements. Throughout the disturbed day, the variations of relative slant TECs between DORIS data and MIDAS results agree quite well, with the average of the mean differences about 2 TECu. We conclude that as a valuable supplement to GPS data, DORIS ionospheric measurements can be used to analyse TEC variations with a relatively high resolution, ∼10 s in time and tens of kilometres in space. This will be very helpful for identification of some highly dynamic structures in the ionosphere found at mid-latitudes, such as the main trough, TID (Travelling Ionospheric Disturbances) and SED (Storm Enhanced Density), and could be used as a valuable auxiliary data source in ionospheric imaging.     

Studies of the geocenter motion using 16-years DORIS data

An accuracy of geocenter motion estimation is strongly dependent on the geodetic network size and stations distribution over the Earth’s surface. From this point of view DORIS system has an advantage, as its ground network of beacons consists of more than 50 sites, equally distributed over the Earth’s surface. Aiming to study variations of the geocenter movements, the results of DORIS data analysis for the time span 1993.0–2009.0 (inawd06.snx series), performed at the Analysis Centre of the Institute of astronomy of the Russian Academy of Sciences, have been used. DORIS data processing was made with GIPSY/OASIS II software, developed by Jet Propulsion Laboratory and modified for DORIS data processing by Institute Géographique National. Standard deviations of stations coordinates are estimated at the level 0.5–4.0 cm (internal consistency), depending on the number of satellites used in the solution. RMS of estimated components of the DORIS satellites orbits, compared with the solutions of other IDS analysis centres, do not exceed 1–2 cm. Weekly solutions for coordinates have been transformed from free network solutions (inawd06.snx series) to a well defined terrestrial reference frame ITRF2005 with the use of seven parameters of Helmert transformation, which were examined with a view to study variations of the geocenter movements (ina05wd01.geoc time series). In order to estimate linear trend, amplitudes, periods and phases of geocenter variation a method of linear regression was applied. The evaluated amplitudes of annual variations are of the order of 5–7 mm for X and Y components and 27–29 mm for Z component. Semi-annual amplitudes are also noticeable in all components (1–34 mm for X, Y and Z components). Secular trends in the DORIS geocenter coordinates are: −1.2, −0.1 and −0.3 mm/year for X, Y and Z directions respectively.     

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.     

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.     

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.     

Contributions of the French Institut Géographique National (IGN) to the International DORIS Service

DORIS is one of the four geodetic space techniques participating to the Global Geodetic Observing System (GGOS). Since the early development of this system, the Institut Géographique National played a specific and active role. Within, the International DORIS Service (IDS), IGN is in a particular position. While it is responsible for the installation and the maintenance of the DORIS ground tracking stations, it also handles one of the two IDS data center delivering DORIS data and products and has been an IDS Analysis Center for years, providing all possible IDS products, in particular the latest ignwd08 time series in preparation of ITRF2008. This paper explains the different aspects of the IGN contribution to IDS from an historical point of view, presents current activities and scientific results and provides a perspective for future activities. Recent DORIS results show a 10 mm precision or better when more than four DORIS satellites are available. Comparisons between recent DORIS solutions (ign07d02 and ign09d02) and past ITRF realizations show that errors are shared between the DORIS and the ITRF realizations. Some problems related to DORIS data processing are also discussed and possible ways to solve them in the future are discussed. In particular, we can now reject the tropospheric origin of the problem detected in the Envisat data after the software upload of October 12, 2004. A few applications in geodesy (terrestrial reference frame, Earth’s polar motion) and geophysics are also discussed as a natural extension of these service-type activities.     

Quality Assessment of the IDS Contribution to ITRF2008

Doppler Orbitography Radiopositionning Integrated by Satellite (DORIS) is one of the four fundamental techniques contributing to the ITRF. The optimal coverage over the globe of the DORIS observing sites and sites co-located with GPS, allow a strong embedding of DORIS within the ITRF network. DORIS contributes to the access to ITRF through precise orbit determination of altimetric satellites with onboard DORIS receivers. The DORIS contribution to the ITRF2008 is enhanced by the fact that the solutions of seven analysis centers were included in the submitted combined time series of weekly station positions and daily polar motion. We evaluate the quality of the DORIS combined solution in terms of its agreement with the other techniques (VLBI, SLR, GPS) contributing to the ITRF2008 combination. We show in particular that the precisions of the current IDS products range between 1.5 to 2.6 mm for station positions (at the epochs of minimum variances); better than 1 mm/yr in velocities and between 170 and 260 micro-arc-seconds for polar motion, a significant improvement by a factor of three to five, compared to past data used in the ITRF2005 combination. This improvement is certainly due to improved analysis strategies employed by the seven IDS analysis centers that contributed to the combined weekly submitted solutions of station positions and polar motion. A spectral analysis of DORIS station height time series indicates that annual and semi-annual signals are dominant. However, TOPEX draconitic period of about 118 days is still detected in about 20% of the station position power spectra. DORIS height annual signals correlate well with GPS annual signal estimated at some co-located stations, which show that DORIS technique is able to detect loading signals.