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.     

Semi-automatic determination of the Azores Current axis using satellite altimetry: Application to the study of the current variability during 1995–2006

Satellite altimetry has been widely used to study the variability of the ocean currents such as the Azores Current (AzC) in the North Atlantic. Most analyses are performed over the region that encloses the current, thus being somehow affected by other oceanographic signals, e.g., eddies. In this study, a new approach for extracting the axis of a zonal current solely based on satellite altimetry is presented. This is a semi-automatic procedure that searches for the maximum values of the gradient of absolute dynamic topography (ADT), using the geostrophic velocity as auxiliary information. The advantage of this approach is to allow the analyses to be performed over a buffer centered on the current axis instead of using a wider region. It is here applied to the AzC for the period June 1995–October 2006.     

Atmospheric correction of satellite altimetry observations and sea-level variability in the NE Atlantic

Satellite altimetry provides continuous and spatially regular measurements of the height of the sea surface. Sea level responds to density changes of the water, to mass changes, due to addition or reduction of water mass, and to changes in the atmosphere above it. The present study examines the influence of atmospheric effects on sea-level variability in the North-East Atlantic. The association between the height of the sea surface and the North Atlantic Oscillation (NAO) is investigated by considering different sets of altimetry measurements for which the atmospheric effects have been handled differently. Altimetry data not corrected for atmospheric effects are strongly anti-correlated with the state of the NAO, reflecting the hydrostatic response of sea-level to the NAO pressure dipole. The application of an atmospheric correction to satellite altimetry observations in the NE Atlantic decreases variability of the height time series by more than 70% and reduces the amplitude of the seasonal cycle by ∼5 cm. Altimetry data for which atmospheric effects are removed via an inverse barometer correction show a non-negligible correlation with the NAO index at some locations suggesting further indirect non-hydrostatic influences of the state of the NAO on sea level variability.     

Remote sensing using GNSS signals: Current status and future directions

The refracted, reflected and scattered signals of global navigation satellite systems (GNSS) have been successfully used to remotely sense the Earth’s surface and atmosphere. It has demonstrated its potential to sense the atmosphere and ionosphere, ocean, land surfaces (including soil moisture) and the cryosphere. These new measurements, although in need of refinement and further validation in many cases, can be used to complement existing techniques and sensors, e.g., radiosonde, ionosonde, radar altimetry and synthetic aperture radar (SAR). This paper presents the current status and new developments of remote sensing using GNSS signals as well as its future directions and applications. Some notable emerging applications include monitoring sea ice, dangerous sea states, ocean eddy and storm surges. With the further improvement of the next generation multi-frequency GNSS systems and receivers and new space-based instruments utilizing GNSS reflections and refractions, new scientific applications of GNSS are expected in various environment remote sensing fields in the near future.     

The X-TRACK/ALES multi-mission processing system: New advances in altimetry towards the coast

In the context of the ESA Climate Change Initiative project, a new coastal sea level altimetry product has been developed in order to support advances in coastal sea level variability studies. Measurements from Jason-1,2&3 missions have been retracked with the Adaptive Leading Edge Subwaveform (ALES) Retracker and then ingested in the X-TRACK software with the best possible set of altimetry corrections. These two coastal altimetry processing approaches, previously successfully validated and applied to coastal sea level research, are combined here for the first time in order to derive a 16-year-long (June 2002 to May 2018), high-resolution (20-Hz), along-track sea level dataset in six regions: Northeast Atlantic, Mediterranean Sea, West Africa, North Indian Ocean, Southeast Asia and Australia. The study demonstrates that this new coastal sea level product called X-TRACK/ALES is able to extend the spatial coverage of sea level altimetry data ~3.5 km in the land direction, when compared to the X-TRACK 1-Hz dataset. We also observe a large improvement in coastal sea level data availability from Jason-1 to Jason-3, with data at 3.6 km, 1.9 km and 0.9 km to the coast on average, for Jason-1, Jason-2 and Jason-3, respectively. When combining measurements from Jason-1 to Jason-3, we reach a distance of 1.2–4 km to the coast. When compared to tide gauge data, the accuracy of the new altimetry near-shore sea level estimations also improves. In terms of correlations with a large set of independent tide gauge observations selected in the six regions, we obtain an average value of 0.77. We also show that it is now possible to derive from the X-TRACK/ALES product an estimation of the ocean current variability up to 5 km to the coast. This new altimetry dataset, freely available, will provide a valuable contribution of altimetry in coastal marine research community.     

Validity assessment of SAMOSA retracking for fully-focused SAR altimeter waveforms

We demonstrate in this work how we can take advantage of known unfocused SAR (UF-SAR) retracking methods (e.g. the physical SAMOSA model) for retracking of fully-focused SAR (FF-SAR) waveforms. Our insights are an important step towards consistent observations of sea surface height, significant wave height and backscatter coefficient (wind speed) with both UF-SAR and FF-SAR. This is of particular interest for SAR altimetry in the coastal zone, since coastal clutter may be filtered out more efficiently in the high-resolution FF-SAR waveform data, which has the potential to improve data quality. We implemented a multi-mission FF-SAR altimetry processor for Sentinel-3 (S3) and Sentinel-6 Michael Freilich (S6), using a back-projection algorithm, and analysed ocean waveform statistics compared to multilooked UF-SAR. We find for Sentinel-3 that the averaged power waveforms of UF-SAR and FF-SAR over ocean are virtually identical, while for Sentinel-6 the FF-SAR power waveforms better resemble the UF-SAR zero-Doppler beam. We can explain and model the similarities and differences in the data via theoretical considerations of the waveform integrals. These findings suggest to use the existing UF-SAR SAMOSA model for retracking S3 FF-SAR waveforms but the SAMOSA zero-Doppler beam model for S6 FF-SAR waveforms, instead. Testing the outlined approach over short track segments, we obtain range biases between UF-SAR and FF-SAR lower than 2 mm and significant wave height biases lower than 5 cm.     

Analysis of coastal altimetry in the Mexican Caribbean

A reprocessing of sea-level anomalies (SLA) resulting from X-TRACK coastal altimetry was carried out for the ENVISAT (2002–2010) and TOPEX/POSEIDON-Jason (1992–2019) satellite missions in the coastal area of the Mexican Caribbean. This consisted of applying a tidal correction to coastal altimetry sea level observations. Harmonic analysis of five coastal tide gauge records was performed to estimate the most important tidal components of the area, resulting on M2, N2, O1, S2, K1, MF, and MM. The tidal signal was reconstructed with the seven tidal components using the TPXO9 model. The SLA signals corrected with the seven tidal components were validated with in situ data from coastal tide gauges. The validation showed that the TPXO9 tidal barotropic model (1/30° grid) used to reconstruct the tidal signal with the seven representative tidal components performed better than the FES2012 global model (1/16° grid) that uses 33 tidal components. The reprocessed SLAs showed clear seasonality with significant signals at 4, 6, and 12 months, with the annual signal being the dominant one. In the Mexican Caribbean coastal zone, oceanographic processes with different scales (from coastal to mesoscale) converge, showing their complexity in the different SLA signals observed. The aim of this work is to contribute to the analysis of coastal altimetry data and understanding the sea level variations in the Mexican Caribbean. This work is the first step in the implementation of methodologies that take advantage of coastal satellite altimetry in the Caribbean Sea.     

基于QZSS L1信号的岸基GNSS-R码延迟海面测高

为了探究GNSS L1波段中信号不同调制方式的GNSS-R码相位延迟海面的测高性能，开展了双天线岸基GNSS-R海面测高实验，收集了采样率为40 MHz的原始中频数据。利用自主开发的GNSS-R测高软件接收机和反演软件对实验数据进行了处理分析，同时获得了基于QZSS L1C/A码和L1C码的GNSS-R海面测高结果，并分别与岸基同步观测的雷达高度计的测高值进行对比，以对反演结果进行精度评定。实验结果表明:基于QZSS L1C/A码和L1C码的海面测高精度最优可达0.63 m和0.4 m，L1C码的延迟GNSS-R测高精度明显高于L1C/A码。此外，GNSS-R测高精度会随着卫星高度角的增加而有所增加。      

Coastal sea level changes in Europe from GPS,tide gauge,satellite altimetry and GRACE, 1993–2011

Sea level changes are threatening the human living environments, particularly along the European Coasts with highly dense population. In this paper, coastal sea level changes in western and southern Europe are investigated for the period 1993–2011 using Global Positioning System (GPS), Tide Gauge (TG), Satellite Altimetry (SA), Gravity Recovery and Climate Experiment (GRACE) and geophysical models. The mean secular trend is 2.26 ± 0.52 mm/y from satellite altimetry, 2.43 ± 0.61 mm/y from TG+GPS and 1.99 ± 0.67 mm/y from GRACE mass plus steric components, which have a remarkably good agreement. For the seasonal variations, annual amplitudes of satellite altimetry and TG+GPS results are almost similar, while GRACE Mass+Steric results are a little smaller. The annual phases agree remarkably well for three independent techniques. The annual cycle is mainly driven by the steric contributions, while the annual phases of non-steric (mass component) sea level changes are almost a half year later than the steric sea level changes.     

First quality assessment of the Cryosat-2 altimetric system over ocean

This paper presents the results of a Calval analysis performed for the Cryosat-2 mission over ocean. The data set used in this analysis consists of products generated by the Cryosat-2 Processor Prototype developed by CNES (Center National d’Etudes Spatiales). This data set has been analysed focusing on LRM (Low Resolution Mode) mode only. One major objective of this paper is to illustrate the potential of Cryosat-2 data over ocean, mainly for waves and Sea Level Anomaly applications. All the results indicate very good performances of the SIRAL (SAR/Interferometric Altimeter) altimeter over ocean. Crossover standard deviation is close to 6.5 cm over the analysed period (3 months) which is close to the Jason-2 and ENVISAT performance. All these results confirm that Cryosat’s altimeter can provide data almost as valuable as other flying altimetric missions, and that it has the potential to contribute to oceanography (e.g. multi-mission climate record, mesoscale monitoring in near real time) and to geodesy (e.g. mean sea surface, bathymetry).     

High-performance altimeter Doppler processing for measuring sea level height under varying sea state conditions

This paper describes an innovative method for processing nadir altimeter data acquired in Synthetic Aperture Radar (SAR) mode, enhancing the system performances over open ocean. Similarly to the current SAR data processing scheme, the so-called LR-RMC (Low Resolution with Range Migration Correction) method, originally designed by Phalippou and Demeester (2011), includes Doppler beam forming, Doppler shift correction and range correction. In LR-RMC, however, an alternative and less complex averaging (stacking) operation is used so that all the Doppler beams produced in a radar cycle (4 bursts of 64 beams for the open-burst Sentinel-3-mode altimeter) are incoherently combined to form a multi-beam echo. In that manner, contrarily to the narrow-band SAR technique, the LR-RMC processing enlarges the effective footprint to average out the effects of surface waves and particularly those from small sub-mesoscale structures (<1 km) that are known to impact SAR-mode performances. On the other hand, the number of averaged beams is as high as in current SAR-mode processing, thus providing a noise reduction at least equally good. The LR-RMC method has the added benefit of reducing the incoherent integration time with respect to the SAR-mode processing (50 ms compared to 2.5 s) limiting possible surface movement effects. By processing one year of Sentinel-3A SRAL SAR-mode data using the LR-RMC method, it is shown that the swell impact on the SAR altimeter performances is totally removed and that an improvement of 10–50% is obtained in the measurement noise of the sea surface height and significant wave height with respect to SAR mode. Additionally, observational capabilities over the middle scales are enhanced potentially allowing the ocean mesoscale features to be retrieved and observations assimilated more usefully in ocean models.     

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.     

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.     

Preliminary gravity recovery from CryoSat-2 data in the Baffin Bay

A number of geophysical phenomenons in the open ocean are still unresolved by conventional altimetry, but could be resolved through the potential improvements offered by Synthetic Aperture Radar (SAR), also called Delay-Doppler, altimetry. The SAR altimeter offers the following benefits with respect to conventional satellite altimetry: factor of 20 improvements in the along-track resolution, the along-track footprint length which does not vary with wave height (sea state), and improved precision in sea surface height measurements or sea surface slope measurements.     

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.     

Interannual variability of the zonal sea surface slope in the equatorial Atlantic during the 1990s

Altimetry data have proven themselves essential for the early detection, analysis and monitoring of large scale tropical anomalies associated with El Niño in the Pacific. Warm events in the Atlantic are much weaker than in the Pacific and are partially masked by the strong seasonal cycle. Satellite altimetric data permits one to estimate the zonal sea surface slope variations at the equator in the Atlantic with sufficient accuracy for resolving interannual sea surface slope variations. The altimetry-derived slope is here shown to detect Atlantic warm events. For all warm events, anomalies in sea surface slope tend to lead SST. In the mid-1990s’, the equatorial interannual variability is dominated by 17-month period events which exhibit the structure observed in local coupled ocean–atmosphere warm events (zonal wind stress weakening and zonal surface slope relaxation, warm SST, excess precipitation). The frequency of occurrence of these Atlantic warm events is seen to have increased during the mid-1990s.     

Spiral eddies in the Baltic,Black and Caspian seas as seen by satellite radar data

Despite spiral eddies were first seen on the sea surface more than 40 years ago, there is still a lot of uncertainty concerning these eddies. The present paper is aimed to provide the comprehensive results on the occurrence and statistics of small-scale eddies in the three inner seas (the Baltic, Black and Caspian) using satellite synthetic aperture radar (SAR) images. The dataset used includes over 2000 medium resolution Envisat ASAR and ERS-2 SAR images obtained in 2009–2010 in the different parts of the seas mentioned. As a result of the analysis performed ∼14,000 vortical structures were detected. 71% of them were visualized due to surfactant films (“black” eddies), while 29% due to wave/current interactions (“white” eddies). Practically all the eddies detected were cyclonically rotating. Their diameter was within 1–20 km. Characteristic size of the “black” eddies in all the basins was discovered to be less than that of the “white” eddies. Characteristic eddy size for the Baltic, Black and Caspian seas proved to be strictly proportional to the values of the baroclinic Rossby radius of deformation typical for these basins. The “black” eddies did not demonstrate a significant connection with the basin- and meso-scale surface circulation of the seas. Most of the “white” eddies detected were attributed to the zones with the most intense drift currents, i.e. those along the western boundaries and (in the Baltic Sea only) in the elongated parts of the basin.     

Vertical TEC in the equatorial anomaly region from satellite altimetry data

The high-precision demands imposed by the ocean altimetry community of the late 1980 resulted in the TOPEX/Poseidon mission. This mission was the first to carry as its main instrument a dual-frequency sea-altimeter on board a satellite. This instrument together with other state-of-the-art technologies involved in the mission, led to sea-height determinations with precision better than 2 cm. As a by-product, the TOPEX/Poseidon mission provided vertical TEC determinations that since they became available, have demonstrated to be a powerful tool for ionospheric studies.     

Mapping mean lake surface from satellite altimetry and GPS kinematic surveys

Lake water height is a key variable in water cycle and climate change studies, which is achievable using satellite altimetry constellation. A method based on data processing of altimetry from several satellites has been developed to interpolate mean lake surface (MLS) over a set of 22 big lakes distributed on the Earth. It has been applied on nadir radar altimeters in Low Resolution Mode (LRM: Jason-3, Saral/AltiKa, CryoSat-2) in Synthetic Aperture Radar (SAR) mode (Sentinel-3A), and in SAR interferometric (SARin) mode (CryoSat-2), and on laser altimetry (ICESat). Validation of the method has been performed using a set of kinematic GPS height profiles from 18 field campaigns over the lake Issykkul, by comparison of altimetry’s height at crossover points for the other lakes and using the laser altimetry on ICESat-2 mission. The precision reached ranges from 3 to 7 cm RMS (Root Mean Square) depending on the lakes. Currently, lake water level inferred from satellite altimetry is provided with respect to an ellipsoid. Ellipsoidal heights are converted into orthométric heights using geoid models interpolated along the satellite tracks. These global geoid models were inferred from geodetic satellite missions coupled with absolute and regional anomaly gravity data sets spread over the Earth. However, the spatial resolution of the current geoid models does not allow capturing short wavelength undulations that may reach decimeters in mountaineering regions or for rift lakes (Baikal, Issykkul, Malawi, Tanganika). We interpolate in this work the geoid height anomalies with three recent geoid models, the EGM2008, XGM2016 and EIGEN-6C4d, and compare them with the Mean Surface of 22 lakes calculated using satellite altimetry. Assuming that MLS mimics the local undulations of the geoid, our study shows that over a large set of lakes (in East Africa, Andean mountain and Central Asia), short wavelength undulations of the geoid in poorly sampled areas can be derived using satellite altimetry. The models used in this study present very similar geographical patterns when compared to MLS. The precision of the models largely depends on the location of the lakes and is about 18 cm, in average over the Earth. MLS can serve as a validation dataset for any future geoid model. It will also be useful for validation of the future mission SWOT (Surface Water and Ocean Topography) which will measure and map water heights over the lakes with a high horizontal resolution of 250 by 250 m.