Atmospheric Corrections for Altimetry Studies over Inland Water

A Dual-Threshold Algorithm for Ice-Covered Lake Water Level Retrieval Using Sentinel-3 SAR Altimetry Waveforms

Satellite altimetry has been proven to measure water levels in lakes and rivers effectively. The Sentinel-3A satellite is equipped with a dual-frequency synthetic aperture radar altimeter (SRAL), which allows for inland water levels to be measured with higher precision and improved spatial resolution. However, in regions at middle and high latitudes, where many lakes are covered by ice during the winter, the non-uniformity of the altimeter footprint can substantially impact the accuracy of water level estimates, resulting in abnormal readings when applying standard SRAL synthetic aperture radar (SAR) waveform retracking algorithms (retrackers). In this study, a modified method is proposed to determine the current surface type of lakes, analyzing changes in backscattering coefficients and brightness temperature. This method aligns with ground station observations and ensures consistent surface type classification. Additionally, a dual-threshold algorithm that addresses the limitations of the original bimodal algorithm by identifying multiple peaks without needing elevation correction is introduced. This innovative approach significantly enhances the precision of equivalent water level measurements for ice-covered lakes. The study retrieves and compares the water level data of nine North American lakes covered by ice from 2016-2019 using the dual-threshold and the SAMOSA-3 algorithm with in situ data. For Lake Athabasca, Cedar Lake, Great Slave Lake, Lake Winnipeg, and Lake Erie, the root mean square error (RMSE) of SAMOSA-3 is 39.58 cm, 46.18 cm, 45.75 cm, 42.64 cm, and 6.89 cm, respectively. However, the dual-threshold algorithm achieves an RMSE of 6.75 cm, 9.47 cm, 5.90 cm, 7.67 cm, and 5.01 cm, respectively, representing a decrease of 75%, 79%, 87%, 82%, and 27%, respectively, compared to SAMOSA-3. The dual-threshold algorithm can accurately estimate water levels in ice-covered lakes during winter. It offers a promising prospect for achieving long-term, continuous, and high-precision water level measurements for middle- and high-latitude lakes.

Synergy between Satellite Altimetry and Optical Water Quality Data towards Improved Estimation of Lakes Ecological Status

European countries are obligated to monitor and estimate ecological status of lakes under European Union Water Framework Directive (2000/60/EC) for sustainable lakes’ ecosystems in the future. In large and shallow lakes, physical, chemical, and biological water quality parameters are influenced by the high natural variability of water level, exceeding anthropogenic variability, and causing large uncertainty to the assessment of ecological status. Correction of metric values used for the assessment of ecological status for the effect of natural water level fluctuation reduces the signal-to-noise ratio in data and decreases the uncertainty of the status estimate. Here we have explored the potential to create synergy between optical and altimetry data for more accurate estimation of ecological status class of lakes. We have combined data from Sentinel-3 Synthetic Aperture Radar Altimeter and Cryosat-2 SAR Interferometric Radar Altimeter to derive water level estimations in order to apply corrections for chlorophyll a, phytoplankton biomass, and Secchi disc depth estimations from Sentinel-3 Ocean and Land Color Instrument data. Long-term in situ data was used to develop the methodology for the correction of water quality data for the effects of water level applicable on the satellite data. The study shows suitability and potential to combine optical and altimetry data to support in situ measurements and thereby support lake monitoring and management. Combination of two different types of satellite data from the continuous Copernicus program will advance the monitoring of lakes and improves the estimation of ecological status under European Union Water Framework Directive.

Monitoring Large-Scale Inland Water Dynamics by Fusing Sentinel-1 SAR and Sentinel-3 Altimetry Data and by Analyzing Causal Effects of Snowmelt

The warming climate is threatening to alter inland water resources on a global scale. Within all waterbody types, lake and river systems are vital not only for natural ecosystems but, also, for human society. Snowmelt phenology is also altered by global warming, and snowmelt is the primary water supply source for many river and lake systems around the globe. Hence, (1) monitoring snowmelt conditions, (2) tracking the dynamics of snowmelt-influenced river and lake systems, and (3) quantifying the causal effect of snowmelt conditions on these waterbodies are critical to understand the cryo-hydrosphere interactions under climate change. Previous studies utilized in-situ or multispectral sensors to track either the surface areas or water levels of waterbodies, which are constrained to small-scale regions and limited by cloud cover, respectively. On the contrary, in the present study, we employed the latest Sentinel-1 synthetic aperture radar (SAR) and Sentinel-3 altimetry data to grant a high-resolution, cloud-free, and illumination-independent comprehensive inland water dynamics monitoring strategy. Moreover, in contrast to previous studies utilizing in-house algorithms, we employed freely available cloud-based services to ensure a broad applicability with high efficiency. Based on altimetry and SAR data, the water level and the water-covered extent (WCE) (surface area of lakes and the flooded area of rivers) can be successfully measured. Furthermore, by fusing the water level and surface area information, for Lake Urmia, we can estimate the hypsometry and derive the water volume change. Additionally, for the Brahmaputra River, the variations of both the water level and the flooded area can be tracked. Last, but not least, together with the wet snow cover extent (WSCE) mapped with SAR imagery, we can analyze the influence of snowmelt conditions on water resource variations. The distributed lag model (DLM) initially developed in the econometrics discipline was employed, and the lagged causal effect of snowmelt conditions on inland water resources was eventually assessed.

Improving Sea Level Anomaly Precision from Satellite Altimetry Using Parameter Correction in the Red Sea

An improved Fourier series model (FSM01) method is used in geophysical and environmental corrections to enhance the final product of the along-track Jason-2 sea level anomaly (SLA) data and extend it near the Red Sea borders. In this study, the ionospheric correction range, wet tropospheric correction range, sea state bias correction range, and dry tropospheric correction range are enhanced and improved using FSM01, which helped to retrieve three more tracks (106, 170, and 234) earlier neglected by the distribution centers and extend the tracks toward the coast. The FSM01 SLA is compared with Jason-2 SLA and Archiving Validation and Interpretation of Satellite Oceanographic (AVISO) SLA for the available five tracks, in which the FSM01 SLA shows a good agreement and higher correlation with the Jason-2 SLA compared with that of AVISO, in addition to filling the gaps in the times series of all tracks. The newly retrieved tracks are also compared with those retrieved by AVISO, and both data points show similar variability, with FSM01 SLA showing no gaps in the time series. The FSM01 SLA was also extended toward the coast and showed high correlation with the coastal tide measurements.

Modelling the Altitude Dependence of the Wet Path Delay for Coastal Altimetry Using 3-D Fields from ERA5

Wet path delay (WPD) for satellite altimetry has been provided from external sources, raising the need of converting this value between different altitudes. The only expression available for this purpose considers the same altitude reduction, irrespective of geographic location and time. The focus of this study is the modelling of the WPD altitude dependence, aiming at developing improved expressions. Using ERA5 pressure level fields (2010–2013), WPD vertical profiles were computed globally. At each location and for each vertical profile, an exponential function was fitted using least squares, determining the corresponding decay coefficient. The time evolution of these coefficients reveals regions where they are highly variable, making this modelling more difficult, and regions where an annual signal exists. The output of this modelling consists of a set of so-called University of Porto (UP) coefficients, dependent on geographic location and time. An assessment with ERA5 data (2014) shows that for the location where the Kouba coefficient results in a maximum Root Mean Square (RMS) error of 3.2 cm, using UP coefficients this value is 1.2 cm. Independent comparisons with WPD derived from Global Navigation Satellite Systems and radiosondes show that the use of UP coefficients instead of Kouba’s leads to a decrease in the RMS error larger than 1 cm.

Daily River Discharge Estimation Using Multi-Mission Radar Altimetry Data and Ensemble Learning Regression in the Lower Mekong River Basin

Estimating river discharge (Q) is critical for ecosystems and water resource management. Traditionally, estimating Q has depended on a single rating curve or the Manning equation. In contrast to the single rating curve, several rating curves at different locations have been linearly combined in an ensemble learning regression method to estimate Q (ELQ) at the Brazzaville gauge station in the central Congo River in a previous study. In this study, we further tested the proposed ELQ and apply it to the Lower Mekong River Basin (LMRB) with three locations: Stung Treng, Kratie, and Tan Chau. Two major advancements for estimating Q with ELQ are presented. First, ELQ successfully estimated Q at Tan Chau, downstream of Kratie, where hydrodynamic complexities exist. Since the hydrologic characteristics downstream of Kratie are extremely diverse and complex in time and space, most previous studies have estimated Q only upstream from Kratie with hydrologic models and statistical methods. Second, we estimated Q over the LMRB using ELQ with water levels (H) obtained from two radar altimetry missions, Envisat and Jason-2, which made it possible to estimate Q seamlessly from 2003 to 2016. Owing to ELQ with multi-mission radar altimetry data, we have overcome the problems of a single rating curve: Locations for estimating Q have to be close to virtual stations, e.g., a few tens of kilometers, because the performance of the single rating curve degrades as the distance between the location of Q estimation and a virtual station increases. Therefore, most previous studies had not used Jason-2 data whose cross-track interval is about 315 km at the equator. On the contrary, several H obtained from Jason-2 altimetry were used in this study regardless of distances from in-situ Q stations since the ELQ method compensates for degradation in the performance for Q estimation due to the poor rating curve with virtual stations away from in-situ Q stations. In general, the ELQ-estimated Q (QELQ) showed more accurate results compared to those obtained from a single rating curve. In the case of Tan Chau, the root mean square error (RMSE) of QELQ decreased by 1504/1338 m3/s using Envisat-derived H for the training/validation datasets. We successfully applied ELQ to the LMRB, which is one of the most complex basins to estimate Q with multi-mission radar altimetry data. Furthermore, our method can be used to obtain finer temporal resolution and enhance the performance of Q estimation with the current altimetry missions, such as Sentinel-3A/B and Jason-3.

Evaluation and Validation of CryoSat-2-Derived Water Levels Using In Situ Lake Data from China

CryoSat-2 altimetry has become a valuable tool for monitoring the water level of lakes. In this study, a concentrated probability density function (PDF) method was proposed for preprocessing CryoSat-2 Geophysical Data Record (GDR) data. CryoSat-2 altimetry water levels were preprocessed and evaluated by in situ gauge data from 12 lakes in China. Results showed that the accuracy of the raw GDR data was limited due to outliers in most of the along-track segments. The outliers were generally significantly lower than the in situ values by several meters, and some by more than 30 m. Outlier detection, therefore, improves upon the accuracy of CryoSat-2 measurements. The concentrated PDF method was able to greatly improve the accuracy of CryoSat-2 measurements. The preprocessed CryoSat-2 measurements were able to observe lake levels with a high accuracy at nine of the twelve lakes, with an absolute mean difference of 0.09 m, an absolute standard deviation difference of 0.04 m, a mean root mean square error of 0.27 m, and a mean correlation coefficient of 0.84. Overall, the accuracy of CryoSat-2-derived lake levels was validated in China. In addition, the accuracy of Database for Hydrological Time Series of Inland Waters (DAHITI) and HYDROWEB water level products was also validated by in situ gauge data.

Analysis of Retrackers’ Performances and Water Level Retrieval over the Ebro River Basin Using Sentinel-3

Satellite altimeters have been used to monitor river and reservoir water levels, from which water storage estimates can be derived. Inland water altimetry can, therefore, play an important role in continental water resource management. Traditionally, satellite altimeters were designed to monitor homogeneous surfaces such as oceans or ice sheets, resulting in poor performance over small inland water bodies due to the contribution from land contamination in the returned waveforms. The advent of synthetic aperture radar (SAR) altimetry (with its improved along-track spatial resolution) has enabled the measurement of inland water levels with a better accuracy and an increased spatial resolution. This study aimed to retrieve water levels from Level-1B Sentinel-3 data with focus on the minimization of the land contamination over small- to middle-sized water bodies (130 m to 4.5 km), where continuous clean waveforms rarely exist. Three specialized algorithms or retrackers, together with a new waveform portion selection method, were evaluated to minimize land contamination in the waveforms and to select the nadir return associated with the water body being overflown. The waveform portion selection method, with consideration of the Digital Elevation Model (DEM), was used to fit the multipeak waveforms that arise when overflying the continental water bodies, exploiting a subwaveform-based approach to pick up the one corresponding to the nadir. The performances of the proposed waveform portion selection method with three retrackers, namely, the threshold retracker, Offset Center of Gravity (OCOG) retracker and two-step SAR physical-based retracker, were compared. No significant difference was found in the results of the three retrackers. However, waveform portion selection using DEM information great improved the results. Time series of water levels were retrieved for water bodies in the Ebro River basin (Spain). The results show good agreement with in situ measurements from the Ebro Reservoir (approximately 1.8 km wide) and Ribarroja Reservoir (approximately 400 m wide), with unbiased root-mean-square errors (RMSEs) down to 0.28 m and 0.16 m, respectively, depending on the retracker.

Absolute Calibration or Validation of the Altimeters on the Sentinel-3A and the Jason-3 over Lake Issykkul (Kyrgyzstan)

Calibration/Validation (C/V) studies using sites in the oceans have a long history and protocols are well established. Over lakes, C/V allows addressing problems such as the performance of the various retracking algorithms and evaluating the accuracy of the geophysical corrections for continental waters. This is achievable when measurements of specific and numerous field campaigns and a ground permanent network of level gauges and weather stations are processed. C/V consists of installation of permanent sites (weather stations, limnigraphs, and GPS reference points) and the organization of regular field campaigns. The lake Issykkul serves as permanent site of C/V, for a multi-mission purpose. The objective of this paper is to calculate the altimeter biases of Jason-3 and Sentinel-3A, both belonging to an operational satellite system which is used for the long-term monitoring of lake level variations. We have also determined the accuracy of the altimeters of these two satellites, through a comparison analysis with in situ data. In 2016 and 2017, three campaigns have been organized over this lake in order to estimate the absolute bias of the nadir altimeter onboard the Jason-3 and Sentinel-3A. The fieldwork consisted of measuring water height using a GPS system, carried on a boat, along the track of the altimeter satellite across the lake. It was performed at the time of the pass of the altimeter. Absolute altimeter biases were calculated by averaging the water height differences along the pass of the satellite (GPS from the boat system versus altimetry). Jason-3 operates in a Low Resolution Mode (LRM), while the Sentinel-3A operates in Synthetic Aperture Radar (SAR) mode. In this study we found that the absolute biases measured for Jason-3 were −28 ± 40 mm with the Ocean retracker and 206 ± 30 mm with the Ice-1 retracker. The biases for Sentinel-3A were −14 ± 20 mm with the Samosa (Ocean like) retracker and 285 ± 20 mm with the OCOG (Ice-1-like) retracker. We have also evaluated the accuracy of these two altimeters over Lake Issykkul which reached to 3 cm, for both the instruments, using the Ocean retracker.