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.     

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.     

Fault rupture model of the 2008 Dangxiong (Tibet,China) Mw 6.3 earthquake from Envisat and ALOS data

On October 6, 2008, an Mw 6.3 earthquake occurred in Dangxiong county, southern Tibetan Plateau. In this study, Synthetic Aperture Radar (SAR) images from Envisat ASAR C-band descending Track 176 and ALOS PALSAR L-band ascending Track 500 are processed to generate the coseismic deformation caused by the earthquake. To estimate the source model, a downhill simplex non-linear inversion method is used to determine the fault rupture geometry, and an automatic fault discretization technique is employed to divide the fault plane to construct the optimal slip model, in which the uncertainties of the fault parameters are assessed by a Monte Carlo method. The inversion results show that the earthquake strikes almost south–north and has a normal faulting focal mechanism with rake angle and slip of −111.7° and 1.33 m, respectively. Peak slip of 2.15 m is located at a depth of 7.5 km. The estimated geodetic moment is 4.06 × 10¹⁸ N m (Mw 6.37), 71.2% of which is released in the depth range 4.5–11 km. The slip model suggests that coseismic slip also takes place at some fault patches near the earth’s surface and postseismic afterslip occurs below the coseismic rupture area after the earthquake.     

An estimation of Envisat’s rotational state accounting for the precession of its rotational axis caused by gravity-gradient torque

The rotational state of Envisat is re-estimated using the specular glint times in optical observation data obtained from 2013 to 2015. The model is simplified to a uniaxial symmetric model with the first order variation of its angular momentum subject to a gravity-gradient torque causing precession around the normal of the orbital plane. The sense of Envisat’s rotation can be derived from observational data, and is found to be opposite to the sense of its orbital motion. The rotational period is estimated to be $(120.674\pm 0.068)·\exp \left((4.5095\pm 0.0096)\times {10}^{-4}·t\right)\text{s}$, where t is measured in days from the beginning of 2013. The standard deviation is 0.760?s, making this the best fit obtained for Envisat in the literature to date. The results demonstrate that the angle between the angular momentum vector and the negative normal of the orbital plane librates around a mean value of $8.53°\pm 0.42°$ with an amplitude from about $0.7°$ (in 2013) to $0.5°$ (in 2015), with the libration period equal to the precession period of the angular momentum, from about 4.8?days (in 2013) to 3.4?days (in 2015). The ratio of the minimum to maximum principal moments of inertia is estimated to be $0.0818\pm 0.0011$, and the initial longitude of the angular momentum in the orbital coordinate system is $40.5°\pm 9.3°$. The direction of the rotation axis derived from our results at September 23, 2013, UTC 20:57 is similar to the results obtained from satellite laser ranging data but about $20°$ closer to the negative normal of the orbital plane.     

Progress in tropospheric ammonia retrieval from the MIPAS satellite instrument

The Michelson Interferometer for Passive Atmospheric Sounding is a polar orbiting high resolution mid-infrared emission limb sounder with a nominal vertical resolution of 3 km. Work to extend the list of routinely retrieved species led to the examination of ammonia, NH₃. In this paper, we explore the feasibility of upper tropospheric retrievals around 6-km altitude, discussing appropriate methods. We demonstrate a weak north to south gradient outside the error bars, using co-addition to improve signal-to-noise and compare these first results with a model.     

Transponder calibration of the Envisat RA-2 altimeter Ku band sigma naught

Although the history of spaceborne altimeters goes back to the early seventies, the absolute calibration of the backscattering coefficient has never been deeply investigated. This information has been primarily used to infer the wind speed via an empirical model, and the intercalibration among different satellite altimeters has revealed to be suitable for this purpose, being the wind retrieval based on an empirical relationship. As far as Ku band system is concerned, the sigma naught absolute calibration of the Envisat altimeter (RA-2) has been performed using an active reference target provided by a transponder. This has been exploited during the 6-month Commissioning phase to generate early calibration results. In order to monitor the RA-2 backscatter calibration during the Envisat lifetime, a continuous calibration effort has been carried out by operating the transponder as much as possible. This paper aims to review the entire effort for calibrating the RA-2 sigma naught measurements, which lasted for almost seven years. It presents in detail the adopted methodology and the final outcome of the activity, providing the users with the correction (bias) to get the calibrated sigma naught and analyzing its stability during almost the entire Envisat lifetime. Specifically, it is concluded that the RA-2 backscatter measurements were quite stable, even if a bias of about 1 dB should be considered with respect to the actually released product. Some small changes in the bias as function of time can be identified during most of the Envisat lifetime, consisting in a slight increase in the first two years, followed by a more stable period and a final drop observed at the end of 2009, until the conclusion of the calibration activity (corresponding to the change in Envisat orbit).     

Detection of Envisat RA2/ICE-1 retracked radar altimetry bias over the Amazon basin rivers using GPS

Altimetry is now routinely used to monitor stage variations over rivers, including in the Amazon basin. It is desirable for hydrologic studies to be able to combine altimetry from different satellite missions with other hydrogeodesy datasets such as leveled gauges and watershed topography. One requirement is to accurately determine altimetry bias, which could be different for river studies from the altimetry calibrated for deep ocean or lake applications. In this study, we estimate the bias in the Envisat ranges derived from the ICE-1 waveform retracking, which are nowadays widely used in hydrologic applications. As a reference, we use an extensive dataset of altitudes of gauge zeros measured by GPS collocated at the gauges. The thirty-nine gauges are spread along the major tributaries of the Amazon basin. The methodology consists in jointly modeling the vertical bias and spatial and temporal slope variations between altimetry series located upstream and downstream of each gauge. The resulting bias of the Envisat ICE-1 retracked altimetry over rivers is 1.044 ± 0.212 m, revealing a significant departure from other Envisat calibrations or from the Jason-2 ICE-1 calibration.     

Spin motion determination of the Envisat satellite through laser ranging measurements from a single pass measured by a single station

The Satellite Laser Ranging (SLR) technology is used to accurately determine the position of space objects equipped with so-called retro-reflectors or retro-reflector arrays (RRA). This type of measurement allows to measure the range to the spacecraft with high precision, which leads to determination of very accurate orbits for these targets. Non-active spacecraft, which are not attitude controlled any longer, tend to start to spin or tumble under influence of the external and internal torques and forces.If the return signal is measured for a non-spherical non-active rotating object, the signal in the range residuals with respect to the reference orbit is more complex. For rotating objects the return signal shows an oscillating pattern or patterns caused by the RRA moving around the satellite’s centre of mass. This behaviour is projected onto the radial component measured by the SLR.In our work, we demonstrate how the SLR ranging technique from one sensor to a satellite equipped with a RRA can be used to precisely determine its spin motion during one passage. Multiple SLR measurements of one target over time allow to accurately monitor spin motion changes which can be further used for attitude predictions. We show our solutions of the spin motion determined for the non-active ESA satellite Envisat obtained from measurements acquired during years 2013–2015 by the Zimmerwald SLR station, Switzerland. All the necessary parameters are defined for our own so-called point-like model which describes the motion of a point in space around the satellite centre of mass.     

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.