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.     

The effect of geocenter motion on Jason-2 orbits and the mean sea level

We compute a series of Jason-2 GPS and SLR/DORIS-based orbits using ITRF2005 and the std0905 standards ( Lemoine et al., 2010). Our GPS and SLR/DORIS orbit data sets span a period of 2 years from cycle 3 (July 2008) to cycle 74 (July 2010). We extract the Jason-2 orbit frame translational parameters per cycle by the means of a Helmert transformation between a set of reference orbits and a set of test orbits. We compare the annual terms of these time-series to the annual terms of two different geocenter motion models where biases and trends have been removed. Subsequently, we include the annual terms of the modeled geocenter motion as a degree-1 loading displacement correction to the GPS and SLR/DORIS tracking network of the POD process. Although the annual geocenter motion correction would reflect a stationary signal in time, under ideal conditions, the whole geocenter motion is a non-stationary process that includes secular trends. Our results suggest that our GSFC Jason-2 GPS-based orbits are closely tied to the center of mass (CM) of the Earth consistent with our current force modeling, whereas GSFC’s SLR/DORIS-based orbits are tied to the origin of ITRF2005, which is the center of figure (CF) for sub-secular scales. We quantify the GPS and SLR/DORIS orbit centering and how this impacts the orbit radial error over the globe, which is assimilated into mean sea level (MSL) error, from the omission of the annual term of the geocenter correction. We find that for the SLR/DORIS std0905 orbits, currently used by the oceanographic community, only the negligence of the annual term of the geocenter motion correction results in a – 4.67 ± 3.40 mm error in the Z-component of the orbit frame which creates 1.06 ± 2.66 mm of systematic error in the MSL estimates, mainly due to the uneven distribution of the oceans between the North and South hemisphere.     

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.     

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.     

From TOPEX/Poseidon to Jason-2/OSTM in the Amazon basin

A major interest of radar altimetry over rivers is to monitor water resources and associated risk in basins where there is little or no conventional in situ data. The objective of the present study is to calibrate altimetry data in a place where conventional data are available, and use the results to estimate the potential error committed in the estimation of water levels in an ungauged or poorly gauged basin. The virtual stations extracted with Jason-2 in this study concern a very broad sample of river channel width and complexity. Minimum channel width has been estimated at 400 m. Unlike TOPEX/Poseidon (T/P), Jason-2 seems to have the capability to distinguish the river bed from its floodplain. The quality of the results obtained with Jason-2 is incomparably better than that obtained with T/P. Despite the fact that no absolute calibration has been assessed for river in this study, the bias calculated converge around 0, 35 m, which could be then the error estimated on the water stage derived from Jason-2 ranges, when no other validation is available. ICE3 algorithm seems to be performing as well as ICE1, and further research is needed to design retracking algorithm specifically for continental water.     

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.     

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.     

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.     

IDS contribution to ITRF2008

For the first time, the International DORIS Service (IDS) has produced a technique level combination based on the contributions of seven analysis centers (ACs), including the European Space Operations Center (ESOC), Geodetic Observatory Pecny (GOP), Geoscience Australia (GAU), the NASA Goddard Space Flight Center (GSFC), the Institut Géographique National (IGN), the Institute of Astronomy, Russian Academy of Sciences (INASAN, named as INA), and CNES/CLS (named as LCA). The ACs used five different software packages to process the DORIS data from 1992 to 2008, including NAPEOS (ESA), Bernese (GOP), GEODYN (GAU, GSC), GIPSY/OASIS (INA), and GINS (LCA). The data from seven DORIS satellites, TOPEX/Poseidon, SPOT-2, SPOT-3, SPOT-4, SPOT-5, Envisat and Jason-1 were processed and all the analysis centers produced weekly SINEX files in either variance–covariance or normal equation format. The processing by the analysis centers used the latest GRACE-derived gravity models, forward modelling of atmospheric gravity, updates to the radiation pressure modelling to improve the DORIS geocenter solutions, denser parameterization of empirically determined drag coefficients to improve station and EOP solutions, especially near the solar maximum in 2001–2002, updated troposphere mapping functions, and an ITRF2005-derived station set for orbit determination, DPOD2005. The CATREF software was used to process the weekly AC solutions, and produce three iterations of an IDS global weekly combination. Between the development of the initial solution IDS-1, and the final solution, IDS-3, the ACs improved their analysis strategies and submitted updated solutions to eliminate troposphere-derived biases in the solution scale, to reduce drag-related degradations in station positioning, and to refine the estimation strategy to improve the combination geocenter solution. An analysis of the frequency content of the individual AC geocenter and scale solutions was used as the basis to define the scale and geocenter of the IDS-3 combination. The final IDS-3 combination has an internal position consistency (WRMS) that is 15 to 20 mm before 2002 and 8 to 10 mm after 2002, when 4 or 5 satellites contribute to the weekly solutions. The final IDS-3 combination includes solutions for 130 DORIS stations on 67 different sites of which 35 have occupations over 16 years (1993.0–2009.0). The EOPs from the IDS-3 combination were compared with the IERS 05 C04 time series and the RMS agreement was 0.24 mas and 0.35 mas for the X and Y components of polar motion. The comparison to ITRF2005 in station position shows an agreement of 6 to 8 mm RMS in horizontal and 10.3 mm in height. The RMS comparison to ITRF2005 in station velocity is at 1.8 mm/year on the East component, to 1.2 mm/year in North component and 1.6 mm/year in height.     

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.     

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.     

Regional in situ validation of satellite altimeters: Calibration and cross-calibration results at the Corsican sites

The in situ validation of the satellite altimeter sea surface heights is generally performed either at a few local points directly flown over by the satellites or using the global tide gauge network. A regional in situ calibration method was developed by NOVELTIS in order to monitor the altimeter data quality in a perimeter of several hundred kilometres around a given in situ calibration site. The primary advantage of this technique is its applicability not only for missions flying over dedicated sites but also for missions on interleaved or non repetitive orbits. This article presents the altimeter bias estimates obtained with this method at the Corsican calibration site, for the Jason-1 mission on its nominal and interleaved orbits as well as for the Jason-2 and Envisat missions. The various regional bias estimates (8.2 cm and 7.4 cm for Jason-1 respectively on the nominal and interleaved orbits in Senetosa, 16.4 cm for Jason-2 in Senetosa and 47.0 cm for Envisat in Ajaccio, with an accuracy between 2.5 cm and 4 cm depending on the mission) are compared with the results obtained by the other in situ calibration teams. This comparison demonstrates the coherency at the centimetre level, the stability and the generic character of the method, which would also be of benefit to the new and future altimeter missions such as Cryosat-2, SARAL/AltiKa, Sentinel-3, Jason-3, Jason-CS.     

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.     

Estimation of differential code biases with Jason-2/3 onboard GPS observations

With the continuous deployment of Low Earth Orbit (LEO) satellites, the estimation of differential code biases (DCBs) based on GNSS observations from LEO has gained increasing attention. Previous studies on LEO-based DCB estimation are usually using the spherical symmetry ionosphere assumption (SSIA), in which a uniform electron density is assumed in a thick shell. In this study, we propose an approach (named the SHLEO method) to simultaneously estimate the satellite and LEO onboard receiver DCBs by modeling the distribution of the global plasmaspheric total electron content (PTEC) above the satellite orbit with a spherical harmonic (SH) function. Compared to the commonly used SSIA method, the SHLEO model improves the GPS satellite DCB estimation accuracy by 13.46% and the stability by 22.34%, respectively. Compared to the GPS satellite DCBs estimated based on the Jason-3-only observations, the accuracy and monthly stability of the satellite DCBs can be improved by 14.42% and 26.8% when both Jason-2 and Jason-3 onboard observations are jointly processed. Compared with the Jason-2 solutions, the GPS satellite DCB estimates based on the fusion of Jason-2 and Jason-3 observations have an improved consistency of better than 18.26% and 9.71% with the products provided by the Center for Orbit Determination in Europe (CODE) and Chinese Academy of Sciences (CAS). Taking the DCB products provided by the German Aerospace Center (DLR) as references, there is no improvement in accuracy of the GPS satellite DCB estimates based on the fusion of Jason-2 and Jason-3 observations than the Jason-2 solutions alone. A periodic variation is found in the time series of both the Jason-3 and Jason-2 onboard receiver DCB estimates. Preliminary analysis of the PTEC distribution based on the estimated SH coefficients are also presented.     

Using altimetry and seafloor pressure data to estimate vertical deformation offshore: Vanuatu case study

Measuring ground deformation underwater is essential for understanding Earth processes at many scales. One important example is subduction zones, which can generate devastating earthquakes and tsunamis, and where the most important deformation signal related to plate locking is usually offshore. We present an improved method for making offshore vertical deformation measurements, that involve combining tide gauge and altimetry data. We present data from two offshore sites located on either side of the plate interface at the New Hebrides subduction zone, where the Australian plate subducts beneath the North Fiji basin. These two sites have been equipped with pressure gauges since 1999, to extend an on-land GPS network across the plate interface. The pressure series measured at both sites show that Wusi Bank, located on the over-riding plate, subsides by 11 ± 4 mm/yr with respect to Sabine Bank, which is located on the down-going plate. By combining water depths derived from the on-bottom pressure data with sea surface heights derived from altimetry data, we determine variations of seafloor heights in a global reference frame. Using altimetry data from TOPEX/Poseidon, Jason-1, Jason-2 and Envisat missions, we find that the vertical motion at Sabine Bank is close to zero and that Wusi Bank subsides by at least 3 mm/yr and probably at most 11 mm/yr.     

Influence of time variable geopotential models on precise orbits of altimetry satellites,global and regional mean sea level trends

During the last decade a significant progress has been reached in the investigation of the gravity field of the Earth. Besides static, also time variable geopotential models have been recently created. In this paper we investigate the impact of the recent time variable geopotential models on altimetry satellite orbits and such altimetry products based on these orbits, as global and regional mean sea level trends. We show that the modeling of time variable gravity improves the orbit solutions, at least for the GRACE period where time variable gravity is sufficiently accurately observed by this mission. Our analysis includes six geopotential models jointly developed by GFZ German Research Centre for Geosciences and Space Geodesy Research Group (CNES/GRGS) Toulouse: the stationary model EIGEN-GL04S, a stationary version of EIGEN-6S (EIGEN-6S_stat), a corrected version of EIGEN-6S and three enhanced versions of EIGEN-6S called EIGEN-6S2, EIGEN-6S2A and EIGEN-6S2B. By “stationary” we mean “containing periodic parameters such as annual and semi-annual variations, but no secular (drift) terms”. We computed precise orbits for the radar altimetry satellites ERS-1, ERS-2, TOPEX/Poseidon, and Envisat over 20 years between 1991 and 2011. The orbit, single-mission and multi-mission altimetry crossover analyses show that the time variable models EIGEN-6S_corrected, EIGEN-6S2 and its two precursors EIGEN-6S2A/B perform notably better than the stationary models for the GRACE period from 2003 onwards. Thus, using EIGEN-6S2 and EIGEN-6S2A/B we have got 3.6% smaller root mean square fits of satellite laser ranging observations for Envisat, as when using EIGEN-GL04S. However, for the pre-GRACE period 1991–2003, the stationary geopotential models EIGEN-GL04S and EIGEN-6S_stat as well as EIGEN-6S2 having no drift terms for degree 3–50 at this time interval perform superior compared to EIGEN-6S_correct and EIGEN-6S2A/B which contain drifts for this period. We found, that the time variable geopotential models have a low (0.1–0.2 mm/yr) impact on our results for the global mean sea level trend. However, we found strong East/West differences up to 3 mm/yr in the regional mean sea level trends when using orbits of all four satellites based on time variable and stationary geopotential models. We show that these differences are related to the relative drifts of the centers-of-origin between the orbit solutions based on the time variable and stationary geopotential models. From the results of our detailed study, we conclude that the final version of the time variable gravity field model EIGEN-6S2 performs best for the four satellites tested. This model provides the most reliable and mission-consistent sea level estimates for the whole time period from 1992 to 2010. This model is of maximum spherical harmonic degree and order 260 and contains time series for drifts as well as annual and semiannual variations of the spherical harmonic coefficients for degree 2–50.     

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.     

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.     

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.