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.     

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.     

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.     

Cross testing of ionosphere models IRI-2001 and IRI-2007, data from satellite altimeters (Topex/Poseidon and Jason-1) and global ionosphere maps

In this study, we use a great body of statistical data covering the entire 23rd solar cycle to cross test data of satellite altimeters, Global Ionosphere Maps and the International Reference Ionosphere models, IRI-2001 and IRI-2007. It is revealed that experimental TEC values of the satellite altimeters regularly exceed the model ones by ∼3 TECU (1 TECU = 10¹⁶ m⁻²). The best possible value of difference between TECs obtained from altimeter and GIM-map data significantly differs for different laboratories: the maximum for CODG data falls on 2.5 TECU, ESAG – 3 TECU, JPLG – 0 TECU, UPCG – 2 TECU. The dependence of experimental and model data root-mean-square deviation on the F10.7 index is shown to be nearly linear. IRI-2001 and IRI-2007 relative errors are characterized by considerable 11-year and annual variations. Given the geomagnetic planetary index K_p under 7, IRI-2001 and IRI-2007 reproduce TEC in the ionosphere with an accuracy of ∼30% relative to measurement data from satellite altimeters. The amplitude of absolute error variations resulting from the difference in ionization enhancement between the model and the real ionosphere during the morning solar terminator transit is ∼5 TECU.     

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 contribution of local gravimetric geoid models to the calibration of satellite altimetry data and an outlook of the latest GOCE GGM performance in Gavdos

The use of geoid heights has been one of the available methodologies utilized for the independent calibration/validation of altimeters on-board satellites. This methodology has been employed for long in the Gavdos dedicated cal/val facility (Crete, Greece), where calibration results for the Jason satellites have been estimated, both for ascending and descending passes. The present work gives a detailed overview of the methodology followed in order to estimate a high-resolution and accuracy gravimetric geoid model for the wider Gavdos area, in support of the on-going calibration work. To estimate the geoid model, the well-known remove-compute-restore method is used while residual geoid heights are estimated through least-squares collocation so that associated errors are determined as well. It is found that the estimated formal geoid errors from LSC along passes 018 and 109 of Jason satellites, used for the bias estimation, range between ±0.8–1.6 cm. The so-derived geoid heights are employed in the determination of the Jason-2 altimeter bias for all available cycles (cycles 1-114, spanning the period from July 2008 to August 2011) together with the RioMed DOT model. From the results acquired the Jason-2 bias has been estimated to be +196.1 ± 3.2 mm for pass 109 and +161.9 ± 5.1 mm for pass 018. Within the same frame, the GOCE/GRACE-based geopotential model GOCO02s has been used to estimate the mean dynamic ocean topography and the steady-state circulation in the area around Gavdos. The so-derived DOT model was used to estimate the Jason-2 bias in an effort to evaluate the performance of satellite-only geoid models and investigate whether their spatial resolution and accuracy provides some improvement w.r.t. traditional local gravimetric geoids. From the results acquired with geoid heights from GOCO02s, the estimated Jason-2 bias deviates significantly from that of the local gravimetric model, which can be attributed to a possible mean offset and the low resolution of GOCE-based GGMs. On the other hand, when the newly estimated GOCE-based DOT was employed with geoid heights from the local gravimetric geoid model, the Jason-2 bias has been estimated to be +185.1 ± 3.2 mm for pass 109 and +130.2 ± 5.1 mm for pass 018.     

The challenges in long-term altimetry calibration for addressing the problem of global sea level change

Long-term change of the global sea level resulting from climate change has become an issue of great societal interest. The advent of the technology of satellite altimetry has modernized the study of sea level on both global and regional scales. In combination with in situ observations of the ocean density and space observations of Earth’s gravity variations, satellite altimetry has become an essential component of a global observing system for monitoring and understanding sea level change. The challenge of making sea level measurements with sufficient accuracy to discern long-term trends and allow the patterns of natural variability to be distinguished from those linked to anthropogenic forcing rests largely on the long-term efforts of altimeter calibration and validation. The issues of long-term calibration for the various components of the altimeter measurement system are reviewed in the paper. The topics include radar altimetry, the effects of tropospheric water vapor, orbit determination, gravity field, tide gauges, and the terrestrial reference frame. The necessity for maintaining a complete calibration effort and the challenges of sustaining it into the future are discussed.     

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.     

R3D-B2 – Measurement of ionizing and solar radiation in open space in the BIOPAN 5 facility outside the FOTON M2 satellite

Solar and space radiation have been monitored using the R3D-B2 radiation risks radiometer-dosimeter on board a recent space flight on the Russian satellite Foton M2 within the ESA Biopan 5 facility mounted on the outside of the satellite exposed to space conditions. The solar radiation has been assayed in four wavelength bands (UV-C, 170–280 nm, UV-B, 280–315 nm), UV-A (315–400 nm) and PAR (photosynthetic active radiation, 400–700 nm). The data show an increasing tumbling rotation of the satellite during the mission. The photodiodes do not show a cosine response to the incident light which has been corrected. After calibration of the signals using the extraterrestrial spectrum, doses have been calculated for each orbit, for each day and for the total mission as basic data for the biological material which has been exposed in parallel in the Biopan facility. Cosmic ionizing radiation has been monitored and separated in 256 deposited energy spectra, which were further used for determination of the absorbed dose rate and flux. Basic data tables were prepared to be used by other Biopan 5 experiments. The paper summarizes the results for the Earth radiation environment at the altitude (262–304 km) of the Foton M2 spacecraft. Comparisons with the predictions of NASA Earth radiation environment experimental models AE-8 and AP-8, and the PSB97 model are also presented, which calculate the fluxes of ionizing radiation from a simulation. AP-8 is a model for trapped radiation.     

Influence of different atmospheric models and seeing effects on the nonlinear calibration

This paper describes a technique intended to improve calibration of vector magnetograms achieved by the spectral line Fei 5324.19 Å which was first described in detail by Hagyard and Kineke [Hagyard, M.J., Kineke, J.I., Improved method for calibrating filter vector magnetographs, Solar Phys., 158, 11–28, 1995.]. The technique attempts to compensate for nonlinear behavior of the measured linear and circular polarizations of the line as functions of inclination angle and field strength. This method depends on different atmospheric models. To the three different atmospheric models, the VAL-C quiet photosphere model, the Ding and Fang penumbra model and Allen umbral model, the circular polarization signal is directly linearly proportional to the longitudinal field strength in the range of 800 G, 1000 G and 1300 G, respectively. For the transverse field, however, this linear relationship holds over much smaller ranges which are 200 G, 300 G and 200 G. The improved method improves the accuracy in calibrating the polarization signals over a wide range of field strengths and inclination angles. With the seeing effects decreasing in the magnetic observation, the longitudinal fields and the transverse ones are underestimated. The differences in the zero-iterative longitudinal and transverse field are relatively minor, but differences in the multi-iterative longitudinal and transverse field can be as high as 260 G and 500 G. The improved calibration method is more sensitive to the seeing effects variation than the former calibration method.     

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.     

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.     

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

Single photon laser altimeter simulator and statistical signal processing

Spaceborne altimeters are common instruments onboard the deep space rendezvous spacecrafts. They provide range and topographic measurements critical in spacecraft navigation. Simultaneously, the receiver part may be utilized for Earth-to-satellite link, one way time transfer, and precise optical radiometry. The main advantage of single photon counting approach is the ability of processing signals with very low signal-to-noise ratio eliminating the need of large telescopes and high power laser source. Extremely small, rugged and compact microchip lasers can be employed. The major limiting factor, on the other hand, is the acquisition time needed to gather sufficient volume of data in repetitive measurements in order to process and evaluate the data appropriately. Statistical signal processing is adopted to detect signals with average strength much lower than one photon per measurement. A comprehensive simulator design and range signal processing algorithm are presented to identify a mission specific altimeter configuration. Typical mission scenarios (celestial body surface landing and topographical mapping) are simulated and evaluated. The high interest and promising single photon altimeter applications are low-orbit (∼10 km) and low-radial velocity (several m/s) topographical mapping (asteroids, Phobos and Deimos) and landing altimetry (∼10 km) where range evaluation repetition rates of ∼100 Hz and 0.1 m precision may be achieved. Moon landing and asteroid Itokawa topographical mapping scenario simulations are discussed in more detail.     

The impact of orbital errors on the estimation of satellite clock errors and PPP

Precise satellite orbit and clocks are essential for providing high accuracy real-time PPP (Precise Point Positioning) service. However, by treating the predicted orbits as fixed, the orbital errors may be partially assimilated by the estimated satellite clock and hence impact the positioning solutions. This paper presents the impact analysis of errors in radial and tangential orbital components on the estimation of satellite clocks and PPP through theoretical study and experimental evaluation. The relationship between the compensation of the orbital errors by the satellite clocks and the satellite-station geometry is discussed in details. Based on the satellite clocks estimated with regional station networks of different sizes (∼100, ∼300, ∼500 and ∼700 km in radius), results indicated that the orbital errors compensated by the satellite clock estimates reduce as the size of the network increases. An interesting regional PPP mode based on the broadcast ephemeris and the corresponding estimated satellite clocks is proposed and evaluated through the numerical study. The impact of orbital errors in the broadcast ephemeris has shown to be negligible for PPP users in a regional network of a radius of ∼300 km, with positioning RMS of about 1.4, 1.4 and 3.7 cm for east, north and up component in the post-mission kinematic mode, comparable with 1.3, 1.3 and 3.6 cm using the precise orbits and the corresponding estimated clocks. Compared with the DGPS and RTK positioning, only the estimated satellite clocks are needed to be disseminated to PPP users for this approach. It can significantly alleviate the communication burdens and therefore can be beneficial to the real time applications.     

DORIS-based point mascons for the long term stability of precise orbit solutions

In recent years non-tidal Time Varying Gravity (TVG) has emerged as the most important contributor in the error budget of Precision Orbit Determination (POD) solutions for altimeter satellites’ orbits. The Gravity Recovery And Climate Experiment (GRACE) mission has provided POD analysts with static and time-varying gravity models that are very accurate over the 2002–2012 time interval, but whose linear rates cannot be safely extrapolated before and after the GRACE lifespan. One such model based on a combination of data from GRACE and Lageos from 2002–2010, is used in the dynamic POD solutions developed for the Geophysical Data Records (GDRs) of the Jason series of altimeter missions and the equivalent products from lower altitude missions such as Envisat, Cryosat-2, and HY-2A. In order to accommodate long-term time-variable gravity variations not included in the background geopotential model, we assess the feasibility of using DORIS data to observe local mass variations using point mascons. In particular, we show that the point-mascon approach can stabilize the geographically correlated orbit errors which are of fundamental interest for the analysis of regional Mean Sea Level trends based on altimeter data, and can therefore provide an interim solution in the event of GRACE data loss. The time series of point-mass solutions for Greenland and Antarctica show good agreement with independent series derived from GRACE data, indicating a mass loss at rate of 210 Gt/year and 110 Gt/year respectively.     

Orbit determination and prediction of GEO satellite of BeiDou during repositioning maneuver

In order to establish a continuous GEO satellite orbit during repositioning maneuvers, a suitable maneuver force model has been established associated with an optimal orbit determination method and strategy. A continuous increasing acceleration is established by constructing a constant force that is equivalent to the pulse force, with the mass of the satellite decreasing throughout maneuver. This acceleration can be added to other accelerations, such as solar radiation, to obtain the continuous acceleration of the satellite. The orbit determination method and strategy are illuminated, with subsequent assessment of the orbit being determined and predicted accordingly. The orbit of the GEO satellite during repositioning maneuver can be determined and predicted by using C-Band pseudo-range observations of the BeiDou GEO satellite with COSPAR ID 2010-001A in 2011 and 2012. The results indicate that observations before maneuver do affect orbit determination and prediction, and should therefore be selected appropriately. A more precise orbit and prediction can be obtained compared to common short arc methods when observations starting 1 day prior the maneuver and 2 h after the maneuver are adopted in POD (Precise Orbit Determination). The achieved URE (User Range Error) under non-consideration of satellite clock errors is better than 2 m within the first 2 h after maneuver, and less than 3 m for further 2 h of orbit prediction.     

In orbit performance of butane propulsion system

The present work focuses on the determination of the in orbit performance of the Alsat-1 microsatellite propulsion system. The satellite mass is 90 kg, of which 6.2 kg is the propulsion system dry mass. The system is a butane propulsion system using low power resistojet thruster with 2.3 kg of propellant. The liquefied butane gas was selected due to its higher storage density and safety compared to the other propellants used for microsatellites. The purpose of this paper is the analysis of the firings performed after the launch of the satellite and to evaluate the system specific impulse and thrust level during the system lifetime. A total of 273 firings were performed on the Alsat-1 propulsion system in the period between the end of 2002 and mid 2009, the cumulated firing time is more than 12 h 49 min. The analysis of all the propulsion telemetry data shows that the system provides a total mission delta V of 25.3 m/s which is more than the 10 m/s specified for this mission. Furthermore, the mission average specific impulse and thrust are respectively 99.9 s and 48.8 mN.     

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.