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.     

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.     

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

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.     

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.     

Integrating Non-Tidal Sea Level data from altimetry and tide gauges for coastal sea level prediction

The main objective of this paper is to integrate Non-Tidal Sea Level (NSL) from the joint TOPEX, Jason-1 and Jason-2 satellite altimetry with tide gauge data at the west and north coast of the United Kingdom for coastal sea level prediction. The temporal correlation coefficient between altimetric NSLs and tide gauge data reaches a maximum higher than 90% for each gauge. The results show that the multivariate regression approach can efficiently integrate the two types of data in the coastal waters of the area. The Multivariate Regression Model is established by integrating the along-track NSL from the joint TOPEX/Jason-1/Jason-2 altimeters with that from eleven tide gauges. The model results give a maximum hindcast skill of 0.95, which means maximum 95% of NSL variance can be explained by the model. The minimum Root Mean Square Error (RMSe) between altimetric observations and model predictions is 4.99 cm in the area. The validation of the model using Envisat satellite altimetric data gives a maximum temporal correlation coefficient of 0.96 and a minimum RMSe of 4.39 cm between altimetric observations and model predictions, respectively. The model is furthermore used to predict high frequency NSL variation (i.e., every 15 min) during a storm surge event at an independent tide gauge station at the Northeast of the UK (Aberdeen).     

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.     

GPS-based sea level measurements to help the characterization of land contamination in coastal areas

The Corsica site has been established in 1996 to perform altimeter calibration on TOPEX/Poseidon and then on its successors Jason-1 and Jason-2. The first chosen location was under the #85 ground track that overflight the Senetosa Cape. In 2005, it was decided to develop another location close to Ajaccio, to be able to perform the calibration of Envisat and in a next future of SARAL/AltiKa that will flight over the same ground tracks. Equipped with various instruments (tide gauges, permanent GPS, GPS buoy, weather station…) the Corsica calibration site is able to quantify the altimeter Sea Surface Height bias but also to give an input on the origin of this bias (range, corrections, orbits, …). Due to the size of Corsica (not a tiny island), the altimeter measurement system (range and corrections) can be contaminated by land. The aim of this paper is to evaluate this land contamination by using GPS measurements from a fixed receiver on land and from another receiver onboard a life buoy. Concerning the altimeter land contamination, we have quantify that this effect can reach 8 mm/km and then affects the Sea Surface Height bias values already published in the framework of the Corsica calibration site by 5–8 mm for TOPEX and Jason missions. On the other hand, the radiometer measurements (wet troposphere correction) are also sensitive to land and we have been able to quantify the level of improvement of a dedicated coastal algorithm that reconciles our results with those coming from other calibration sites. Finally, we have also shown that the standard deviation of the GPS buoy sea level measurements is highly correlated (∼87%) with the Significant Wave Height derived from the altimeters and can be used to validate such parameter.     

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.     

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.     

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

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

GPS monitoring of vertical seafloor motion at Platform Harvest

We describe results from two decades of monitoring vertical seafloor motion at the Harvest oil platform, NASA’s prime verification site for the TOPEX/Poseidon and Jason series of reference altimeter missions. Using continuous GPS observations, we refine estimates of the platform subsidence—due most likely to fluid withdrawal linked to oil production—and describe the impact on estimates of stability for the altimeter measurement systems. The cumulative seafloor subsidence over 20 yrs is approximately 10 cm, but the rate does not appear constant. The apparent non-linear nature of the vertical motion, coupled with long-period GPS errors, implies that the quality of the seafloor motion estimates is not uniform over the 20-yr period. For the Jason-1 era (2002–2009), competing estimates for the subsidence show agreement to better than 1 mm yr⁻¹. Longer durations of data are needed before the seafloor motion estimates for the Jason-2 era (2008–present) can approach this level of accuracy.     

Dayside ionospheric uplift during strong geomagnetic storms as detected by the CHAMP,SAC-C,TOPEX and Jason-1 satellites

Ionosphere response to severe geomagnetic storms that occurred in 2001–2003 was analyzed using data of global ionosphere maps (GIM), altimeter data from the Jason-1 and TOPEX satellites, and data of GPS receivers on-board CHAMP and SAC-C satellites. This allowed us to study in detail ionosphere redistribution due to geomagnetic storms, dayside ionospheric uplift and overall dayside TEC increase. It is shown that after the interplanetary magnetic field turns southward and intensifies, the crests of the equatorial ionization anomaly (EIA) travel poleward and the TEC value within the EIA area increases significantly (up to ∼50%). GPS data from the SAC-C satellite show that during the main phase of geomagnetic storms TEC values above the altitude of 715 km are 2–3 times higher than during undisturbed conditions. These effects of dayside ionospheric uplift occur owing to the “super-fountain effect” and last few hours while the enhanced interplanetary electric field impinged on the magnetopause.     

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

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

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.     

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

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

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.     

Coastal sea level measurements using a single geodetic GPS receiver

This paper presents a method to derive local sea level variations using data from a single geodetic-quality Global Navigation Satellite System (GNSS) receiver using GPS (Global Positioning System) signals. This method is based on multipath theory for specular reflections and the use of Signal-to-Noise Ratio (SNR) data. The technique could be valuable for altimeter calibration and validation. Data from two test sites, a dedicated GPS tide gauge at the Onsala Space Observatory (OSO) in Sweden and the Friday Harbor GPS site of the EarthScope Plate Boundary Observatory (PBO) in USA, are analyzed. The sea level results are compared to independently observed sea level data from nearby and in situ tide gauges. For OSO, the Root-Mean-Square (RMS) agreement is better than 5 cm, while it is in the order of 10 cm for Friday Harbor. The correlation coefficients are better than 0.97 for both sites. For OSO, the SNR-based results are also compared with results from a geodetic analysis of GPS data of a two receivers/antennae tide gauge installation. The SNR-based analysis results in a slightly worse RMS agreement with respect to the independent tide gauge data than the geodetic analysis (4.8 cm and 4.0 cm, respectively). However, it provides results even for rough sea surface conditions when the two receivers/antennae installation no longer records the necessary data for a geodetic analysis.