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.     

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

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.     

The International DORIS Service (IDS): Toward maturity

DORIS is one of the four space-geodetic techniques participating in the Global Geodetic Observing System (GGOS), particularly to maintain and disseminate the Terrestrial Reference Frame as determined by International Earth rotation and Reference frame Service (IERS). A few years ago, under the umbrella of the International Association of Geodesy, a DORIS International Service (IDS) was created in order to foster international cooperation and to provide new scientific products. This paper addresses the organizational aspects of the IDS and presents some recent DORIS scientific results. It is for the first time that, in preparation of the ITRF2008, seven Analysis Centers (AC’s) contributed to derive long-term time series of DORIS stations positions. These solutions were then combined into a homogeneous time series IDS-2 for which a precision of less than 10 mm was obtained. Orbit comparisons between the various AC’s showed an excellent agreement in the radial component, both for the SPOT satellites (e.g. 0.5–2.1 cm RMS for SPOT-2) and Envisat (0.9–2.1 cm RMS), using different software packages, models, corrections and analysis strategies. There is now a wide international participation within IDS that should lead to future improvements in DORIS analysis strategies and DORIS-derived geodetic products.     

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.     

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.     

Improving DORIS geocenter time series using an empirical rescaling of solar radiation pressure models

Even if Satellite Laser Ranging (SLR) remains the fundamental technique for geocenter monitoring, DORIS can also determine this geophysical parameter. Gobinddass et al. (2009) found that part of the systematic errors at 118 days and 1 year can be significantly reduced by rescaling the current solar radiation pressure models using satellite-dependent empirical models. Here we extend this study to all DORIS satellites and propose a complete set of empirical solar radiation parameter coefficients. A specific problem related to SPOT-5 solar panel realignment is also detected and explained. New DORIS geocenter solutions now show a much better agreement in amplitude with independent SLR solutions and with recent geophysical models. Finally, the impact of this refined DORIS data strategy is discussed in terms of Z-geocenter monitoring as well as for other geodetic products (altitude of high latitude station such as Thule in Greenland) and Precise Orbit Determination. After reprocessing the full 1993.0-2008.0 DORIS data set, we confirm that the proposed strategy allows a significant reduction of systematic errors at periods of 118 days and 1 year (up to 20 mm), especially for the most recent data after 2002.5, when more DORIS satellites are available for geodetic purposes.     

Impact of DORIS ground antennas environment on their radio signal quality

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

Refining DORIS atmospheric drag estimation in preparation of ITRF2008

In preparation of ITRF2008, all geodetic technique services (VLBI, SLR, GPS and DORIS) are generating new solutions based on combination of individual analysis centers solutions. These data reprocessing are based on a selection of models, parameterization and estimation strategy unique to each analysis center and to each technique. While a good agreement can be found for models between groups, thanks to the existence of the IERS conventions, a great diversity still exist for parameter estimation, allowing possible future improvements in this direction. The goal of this study is to focus on the atmospheric drag estimation used to generate the new DORIS/IGN ignwd08 time series prepared for ITRF2008. We develop here a method to inter-compare different processing strategies. In a first step, by analyzing single-satellite solutions for a few weeks of data but for a large number of possible analysis strategies, we demonstrate that estimating drag coefficient more frequently (typically every 1–2 h instead of previously every 4–8 h) for the lowest DORIS satellites (SPOTs and Envisat) provides better geodetic results for station coordinates and polar motion. This new processing strategy also solved earlier problem found when processing DORIS data during intense geomagnetic events, such as geomagnetic storms. Differences between drag estimation strategies can mostly be found during these few specific periods of extreme geomagnetic activity (few days per year). In such a case, when drag coefficient is only estimated every 6 h or less often for single-satellite solution, a significant degradation in station coordinate accuracy can be observed (120 mm vs. 20 mm) and significant biases arose in polar motion estimation (5 mas vs. 0.3 mas). In a second step, we reprocessed a full year of DORIS data (2003) in a standard multi-satellite mode. We were able to provide statistics on a more reliable data set and to strengthen these conclusions. Our proposed DORIS analysis is easy to implement in all software packages and is now already used by several analysis centers of the International DORIS Service (IDS) when submitting reprocessed solutions for ITRF2008.     

SPOT-2 and TOPEX/Poseidon precise orbit determination from DORIS doppler tracking

The French earth observation satellite SPOT-2 has served as a testbed for precise orbit determination from DORIS doppler tracking in anticipation of the TOPEX/Poseidon mission. Using the most up-to-data gravity field model, JGM-2, a radial orbit accuracy of about 2–9 cm was achieved, with an rms of fit of the tracking data of about 0.64 mm/s. Furthermore, it was found that the coordinates of the ground stations can be determined with an accuracy of the order of 2–5 cm after removal of common rotations, and translations.

Using a slightly different model for atmospheric drag, but the same gravity model, precise orbits of TOPEX/Poseidon from DORIS tracking data were determined with a radial orbit accuracy of the order of 4–5 cm, which is far within the 13 cm mission requirement. This conclusion is based on the analysis of 1-day overlap of successive 11-day orbits, and the comparisons with orbits computed from satellite laser tracking (SLR) and from the combination of SLR and DORIS tracking. Results indicate a consistency between the different orbits of 1–4 cm, 4–20 cm, and 6–13 cm in the radial, cross-track, and along-track directions, respectively. The residual rms is about 4–5 cm for SLR data and 0.56 mm/s for DORIS tracking. These numbers are roughly twice as large as the system noise levels, reflecting the fact that there are still some modeling errors left.     

Impact of ITRS 2014 realizations on altimeter satellite precise orbit determination

This paper evaluates orbit accuracy and systematic error for altimeter satellite precise orbit determination on TOPEX, Jason-1, Jason-2 and Jason-3 by comparing the use of four SLR/DORIS station complements from the International Terrestrial Reference System (ITRS) 2014 realizations with those based on ITRF2008. The new Terrestrial Reference Frame 2014 (TRF2014) station complements include ITRS realizations from the Institut National de l’Information Géographique et Forestière (IGN) ITRF2014, the Jet Propulsion Laboratory (JPL) JTRF2014, the Deutsche Geodätisches Forschungsinstitut (DGFI) DTRF2014, and the DORIS extension to ITRF2014 for Precise Orbit Determination, DPOD2014. The largest source of error stems from ITRF2008 station position extrapolation past the 2009 solution end time. The TRF2014 SLR/DORIS complement impact on the ITRF2008 orbit is only 1–2 mm RMS radial difference between 1992–2009, and increases after 2009, up to 5 mm RMS radial difference in 2016. Residual analysis shows that station position extrapolation error past the solution span becomes evident even after two years, and will contribute to about 3–4 mm radial orbit error after seven years. Crossover data show the DTRF2014 orbits are the most accurate for the TOPEX and Jason-2 test periods, and the JTRF2014 orbits for the Jason-1 period. However for the 2016 Jason-3 test period only the DPOD2014-based orbits show a strong and statistically significant margin of improvement. The positive results with DTRF2014 suggest the new approach to correct station positions or normal equations for non-tidal loading before combination is beneficial. We did not find any compelling POD advantage in using non-linear over linear station velocity models in our SLR & DORIS orbit tests on the Jason satellites. The JTRF2014 proof-of-concept ITRS realization demonstrates the need for improved SLR+DORIS orbit centering when compared to the Ries (2013) CM annual model. Orbit centering error is seen as an annual radial signal of 0.4 mm amplitude with the CM model. The unmodeled CM signals show roughly a 1.8 mm peak-to-peak annual variation in the orbit radial component. We find the TRF network stability pertinent to POD can be defined only by examination of the orbit-specific tracking network time series. Drift stability between the ITRF2008 and the other TRF2014-based orbits is very high, the relative mean radial drift error over water is no larger than 0.04 mm/year over 1993–2015. Analyses also show TRF induced orbit error meets current altimeter rate accuracy goals for global and regional sea level estimation.     

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.     

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

Laser measurements to space debris from Graz SLR station

In order to test laser ranging possibilities to space debris objects, the Satellite Laser Ranging (SLR) Station Graz installed a frequency doubled Nd:YAG pulse laser with a 1 kHz repetition rate, a pulse width of 10 ns, and a pulse energy of 25 mJ at 532 nm (on loan from German Aerospace Center Stuttgart – DLR). We developed and built low-noise single-photon detection units to enable laser ranging to targets with inaccurate orbit predictions, and adapted our standard SLR software to include a few hundred space debris targets. With this configuration, we successfully tracked – within 13 early-evening sessions of each about 1.5 h – 85 passes of 43 different space debris targets, in distances between 600 km and up to more than 2500 km, with radar cross sections from >15 m² down to <0.3 m², and measured their distances with an average precision of about 0.7 m RMS.     

Status of the T2L2/Jason2 Experiment

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

Contribution of satellite laser ranging to combined gravity field models

In the framework of satellite-only gravity field modeling, satellite laser ranging (SLR) data is typically exploited to recover long-wavelength features. This contribution provides a detailed discussion of the SLR component of GOCO02S, the latest release of combined models within the GOCO series. Over a period of five years (January 2006 to December 2010), observations to LAGEOS-1, LAGEOS-2, Ajisai, Stella, and Starlette were analyzed. We conducted a series of closed-loop simulations and found that estimating monthly sets of spherical harmonic coefficients beyond degree five leads to exceedingly ill-posed normal equation systems. Therefore, we adopted degree five as the spectral resolution for real data analysis. We compared our monthly coefficient estimates of degree two with SLR and Gravity Recovery and Climate Experiment (GRACE) time series provided by the Center for Space Research (CSR) at Austin, Texas. Significant deviations in C₂₀ were noted between SLR and GRACE; the agreement is better for the non-zonal coefficients. Fitting sinusoids together with a linear trend to our C₂₀ time series yielded a rate of (−1.75 ± 0.6) × 10⁻¹¹/yr; this drift is equivalent to a geoid change from pole to equator of 0.35 ± 0.12 mm/yr or an apparent Greenland mass loss of 178.5 ± 61.2 km³/yr. The mean of all monthly solutions, averaged over the five-year period, served as input for the satellite-only model GOCO02S. The contribution of SLR to the combined gravity field model is highest for C₂₀, and hence is essential for the determination of the Earth’s oblateness.