Mangroves are more and more recognized for nature-based mitigation of flood risks in low-lying coastal zones, which host a disproportionately large part of the global population. Aerial roots, branches, and canopy of mangrove trees exert friction on the water flow, thereby reducing the propagation of high water levels through the forest. Field measurements of high water level attenuation rates are limited so far to mangrove forests situated at higher-latitude subtropical settings, where aerial roots are much sparser and lower than in low-latitude tropical mangroves. Here, for the first time, we measured high water level attenuation in a tropical Rhizophora forest, where aerial roots are several meters high and water levels never exceeded the aerial root height. Our measurements reveal attenuation rates between 42 ± 9.8 and 46 ± 9.8 cm km−1, which are the highest attenuation rates ever recorded in a mangrove forest, but an exponential rate is more suited to quantify high water level attenuation. In contrast to observations inside the mangrove forest, our observations showed that the propagation of high water levels through a 20 km long tidal channel fringed by wide mangrove areas was amplified, but that high water level amplification was reduced for higher tides with deeper flooding of the fringing mangroves. Our results provide the first empirical assessment of flood protection by tropical Rhizophora mangroves. As Rhizophora is globally the most common genus among mangroves, we propose that our reported high attenuation rates should be incorporated in future assessments of nature-based flood risk mitigation by mangroves.

In recent decades, satellite radar altimetry has been widely used to assess volume changes over the Greenland Ice Sheet. In particular, melt events result in drastic changes in the volume scattering of firn, which induces a pronounced change in the parameters derived from radar altimeter data. Due to the recent and increasingly frequent melt events over Greenland, the impacts of these events on the firn condition, i.e. the formation of ice lenses and reduction in firn air content, need to be better understood. This study therefore exploits the ability of long-term CryoSat-2 data to indicate changes in firn volume scattering in order to assess the spatiotemporal firn condition variations in Greenland. More specifically, this study utilises the leading edge width (LeW) parameter derived from CryoSat-2 Low Resolution Mode data, which has been proven to be a parameter strongly sensitive to changes in volume scattering, and assesses its variation between September 2010 and September 2024. With a combined analysis of remote sensing observations, in situ observations, and outputs from regional climate models, our study demonstrates that the LeW drop induced by extreme melt events in the interior of Greenland experiences a gradual recovery, which can potentially be explained by new-snow deposition. However, in many high-elevation regions of Greenland where firn layers were originally dry, the recent recurrence of extensive melt has prevented a full recovery of the firn volume scattering to pre-2012 conditions, indicating a persistent increase in firn density under a changing climate. Finally, our study also confirms the utility of radar altimeter data for long-term monitoring of the impact of melt and refreezing events on the properties of the upper firn layer

Tidal models that incorporate satellite altimeter data have historically shown discrepancies in accuracy between shallow and deep marine environments. A recent study suggests that these differences may partly stem from neglecting the nonlinear tide-surge interactions in tidal analyses. In this study, we introduce a novel method for estimating tidal constituents from satellite altimeter data in shallow waters, leveraging a 2D hydrodynamic model that accounts for these nonlinear interactions. This approach substantially reduces the variance of unaccounted water level variability, thereby benefiting the estimation. A distinctive feature of our method is the treatment of prior model tidal constituents as stochastic, which helps manage the low temporal resolution of altimeter data by ensuring that unresolved tidal constituents are not updated. We tested our method in the data-rich northwest European continental shelf region, using the high-resolution 2D Dutch Continental Shelf Model version 7 (DCSM). Results show a substantial reduction in the standard deviations of residual water level time series in the shallow waters around Great Britain and in the German Bight, from 11 cm to 5 cm. In deep waters (>200 m), the median standard deviation decreased from 6.8 cm to 6.2 cm. When compared to state-of-the-art ocean tide and surge corrections from publicly available models, our method outperformed them in shallow waters (median standard deviation of 6.0 cm versus 7.5 cm), though the alternative products performed better in deep waters (median standard deviation of 5.5 cm versus 6.2 cm). An estimate of the accuracy at satellite crossovers resulted in an estimated total tidal error of about 1.5 cm (RSS VD). We acknowledge that comparisons in shallow waters are complicated, as alternative products do not account for nonlinear tide-surge interactions. Overall, the demonstration along-track tidal product developed in this study shows potential for improving the tidal representation in the DCSM model. In data-poor regions, the number of tidal constituents that can be reliably estimated using the method may be limited, and alternative strategies might be needed to evaluate the model’s uncertainty in representing tides.

Measurements of tides are relatively sparse in the Arctic. This paper studies GNSS buoy tracks to complement existing data. Existing methods to perform tidal harmonic analysis of the buoy data are inadequate in the Arctic region because these methods for tidal analysis combine data from multiple buoy tracks, which is often infeasible in the Arctic. Moreover, we find that there are significant spatial and temporal variations in amplitudes and phases in baroclinic zones. To address these complexities, we introduce a new approach–Model-derived Fitting Method–to estimate the tidal current constituents (TCC) from a single buoy trajectory. Our study assesses the proposed method by analyzing GNSS buoy data from three Arctic regions characterized by barotropic or baroclinic tidal currents. Through detailed case studies in the Barents Sea, Chukchi Sea, and Baffin Bay, our approach demonstrates accuracy, robustness, and operational capabilities. In the Barents Sea, TCC estimates from two buoys were compared at a common location within their trajectories and compared against model estimates. In the Chukchi Sea's barotropic dominant zone, our method's estimates were evaluated against nearby ADCP mooring data. In Baffin Bay, known for baroclinic currents, a synthetic evaluation confirmed the method's effectiveness. Our study also highlights that phase variations along buoy trajectories can lead to frequency shifts in the spectrum, similar to the Doppler shift effect, particularly notable in regions with baroclinic tides.

Arctic sea ice leads to a significant dissipation of tidal energy, necessitating its inclusion in global tidal models. However, most global tidal models either neglect or only partially incorporate the impact of sea ice on tides. This study proposes a method to model the dissipative forces exerted by sea ice on tides without directly coupling to a sea ice model, yet utilizing sea ice parameters such as thickness and concentration. Our approach involves (re)-categorizing the sea ice cover into regions dominated either by the velocity difference between sea ice and tides (Vertical Shear (VS)) or by the shear from drifting sea ice on tides (Horizontal Shear (HS)), which primarily govern the energy dissipation between tides and sea ice. The subdivision and resulting areas of these HS and VS regions are based on a nondimensional number referred to as the Friction number, which is the ratio of the internal stress of the sea ice field to the ice–water frictional stress, and directly depends on the thickness and concentration of the sea ice. The new parameterization is validated through a performance assessment comparing it to a commonly used approach of assuming all the sea ice to be stationary (landfast). The seasonal modulation of the M₂ tidal component, quantified as the March–September difference, serves as the performance metric, demonstrating that the new parameterization has better agreement with observations from altimeter- and tide gauge-derived seasonal modulation. The results indicate that the physics of ice–tide interaction is better represented with the new parameterization of sea ice-induced dissipation, making it suitable for investigating the effects of declining sea ice thickness on tides.

Traditionally, nadir-looking satellite radar altimeters provide water surface elevations of rivers only at intersections with the satellite’s ground track, called virtual stations. These observations have limited spatial coverage because such cross-overs are sparse, depending on the altimeter's orbit. In this work, we introduce the novel concept of Polygon-Informed Cross-Track Altimetry (PICTA), enabling precise measurement of water surface elevations at cross-track distances---for as long as the target's signal is recorded in the altimeter's range window. Using fully-focused SAR data from the Sentinel-6 altimetry mission, we demonstrate how the new approach can provide detailed river elevation profiles over a ground swath of about 14 km cross-track width and with an along-track resolution as fine as 10 m. On the one hand, this marks a drastic improvement in the number of available measurements when compared to the virtual station approach, on the other hand, for the first time, water surface slopes and level variations along the river, caused by rapids, dams, and sluices, can be directly observed using a nadir radar altimeter. The validation over two river segments in France reveals biases as low as [[EQUATION]]4 cm and random errors on the order of 3 to 8 cm at 30 m along-track resolution. The new PICTA concept can potentially be generalized to other targets, such as lakes or even coastlines.

Traditionally, nadir-looking satellite radar altimeters provide water levels of rivers only at intersections with the satellite's ground track, called virtual stations. These observations have limited spatial coverage because such cross-overs are sparse, depending on the altimeter's orbit. In this work, we introduce the novel concept of Polygon-Informed Cross-Track Altimetry (PICTA), enabling accurate estimation of water levels at cross-track distances — for as long as the target's signal is recorded in the altimeter's range window. Using fully-focused SAR data from the Sentinel-6 altimetry mission, we demonstrate how the new approach can provide detailed river water level profiles over a ground swath of about 14 km cross-track width and with an along-track resolution as fine as 10 m. On the one hand, this marks a drastic improvement in the number of available measurements when compared to the virtual station approach, on the other hand, for the first time, water surface slopes and level variations along the river, caused by rapids, dams, and sluices, can be directly observed using a nadir radar altimeter. The validation over two river segments in France reveals biases as low as ±4 cm and random errors on the order of 3–8 cm at 30 m along-track resolution. The new PICTA concept can potentially be generalized to other targets, such as lakes or even coastlines.

Multi-mission satellite altimetry data have been used to study the spatial and temporal variability in global storm surge water levels. This was done by means of a time-dependent extreme value analysis applied to the monthly maximum detided water levels. To account for the limited temporal resolution of the satellite data, the data were first stacked on a 5^∘× 5^∘ grid. Moreover, additional scaling was applied to the extreme value analysis for which the scaling factors were determined by means of a resampling method using reanalysis data. In addition to the conventional analysis using data from tide gauges, this study provides an insight in the ocean-wide storm surge properties. Nonetheless, where possible, results were compared to similar information derived from tide gauge data. Except for secular changes, the satellite-derived results are comparable to the information derived from tide gauges (correlation > 0.5), although the tide gauges show more local variability. Where limited correlation was observed for the secular change, it was suggested that the satellites may not be able to fully capture the temporal variability in the short-lived, tropical storms, as opposed to extra-tropical storms.

Earlier work has empirically demonstrated some advantages of an increased posting rate of Synthetic Aperture Radar (SAR) altimeters beyond the expected ground resolution of about 320 m in Delay-Doppler (unfocused SAR, UFSAR) processing, corresponding to ∼20 Hz sampling. Higher posting rates of 40–80 Hz were shown to prevent spectral aliasing of the signal, enable to measure swell wave related signal distortions and may lead to a reduced root mean square error of 1 Hz estimates of Sea Surface Heights (SSH), radar cross section (sigma0) and Significant Wave Heights (SWH) from current SAR altimeters. These improvements were explained by the narrow noise autocovariance function of the waveform signal's power speckle noise in along-track direction on one hand, and frequency doubling by power detection (squaring of the signal) on the other. It has not been explained, however, why the power speckle noise decorrelates faster than anticipated by the predicted Doppler resolution, and whether this decorrelation depends on the altimeter and processing configuration. Also, it has not been shown explicitly that the estimates of SSH, SWH and sigma0 decorrelate in the same way. Describing the noise autocovariance function – or equivalently the noise power spectral density via the Wiener-Kintchin theorem – is necessary on two counts: Knowing the noise autocovariance allows to apply optimal filtering strategies that maximize precision on one hand, while the noise power spectral density predicts the frequencies contained in the noise (and signals), which in turn determines the required sampling frequency according to the Nyquist theorem. Using a newly derived analytic noise autocovariance model for UFSAR-processed altimeter data, we show that the swift signal decorrelation is mainly due to the observation geometry. Furthermore, our results demonstrate that the noise autocovariance functions of power speckle, SSH, SWH and sigma0 estimates in along-track direction are different and depend on the sea state. On top of that, the noise autocovariance functions are strongly dependent on the number of Doppler beams used for multilooking, the used retracker, and the processing choices such as antenna gain pattern compensation and windowing within the UFSAR processing (Level-1b). We validated our noise autocovariance model with segments of 42 Sentinel-3B overpasses. Our findings are in accordance to all earlier work, but indicate that the reported precision improvements with respect to 20 Hz may have been too optimistic and that the SSH, SWH and sigma0 generally decorrelate slower than the power speckle noise. We found that the required posting rate is always higher or equal to 40 Hz. Our results will potentially enable improved spectral analysis and optimal filtering of any UFSAR altimetry data. More importantly, our results can be used to trade off different aspects for determining an optimal posting rate in UFSAR altimeter processing in different sea states and with changing processing parameters, which is necessary in view of strict precision requirements of existing and future SAR altimetry missions.

This paper investigates the full variance–covariance (VC) matrix of two high-resolution regional quasi-geoid models, utilizing a spherical radial basis function parameterization. Model parameters were estimated using weighted least-squares techniques and variance component estimation (VCE) for data weighting. The first model, known as the “RCR model,” is computed through the remove–compute–restore method, incorporating various local gravity and radar altimeter datasets. The second model, the “combined model,” includes the GOCO05s satellite-only global geopotential model as an additional dataset with a full-noise VC matrix. Validation of the noise VC matrix scaling for each quasi-geoid model is achieved by comparing observed and formal noise standard deviations of differences between geometric and gravimetric height anomalies at GPS height markers in the Netherlands. Analysis of the noise VC matrix of height anomalies at grid nodes reveals significantly smaller formal noise standard deviations for the RCR model compared to the combined model. This difference is attributed to VCE assigning larger weights to the GOCO05s dataset, which exhibits greater noise standard deviations for the specific spatial scales used. Additionally, the formal noise standard deviations of height anomaly differences, relevant for GNSS-heighting, favor the RCR model. However, the disparity between the two models is smaller than implied by the height anomaly noise standard deviations. This is due to the combined model’s noise autocorrelation function displaying a longer correlation length (67 km) in contrast to the RCR model’s (17 km). Consequently, the combined model exhibits a greater reduction in noise variance for height anomaly differences relative to white noise compared to the RCR model.

All realizations of the European Vertical Reference System (EVRS) computed so far are solely based on geopotential differences obtained by spirit leveling/gravimetry. As such, there are no direct connections between height benchmarks separated by large water bodies. In this study, such connections are added by means of model-based hydrodynamic leveling resulting in a new, yet unofficial realization of the EVRS. The model-derived mean water levels used in computing the hydrodynamic leveling connections were obtained from the Nemo-Nordic (Baltic Sea) and 3D DCSM-FM (northwest European continental shelf) hydrodynamic models. The impact of model-based hydrodynamic leveling on the European Vertical Reference Frame is significant, especially for France and Great Britain. Compared to a solution which only uses spirit leveling/gravimetry, the differences in these countries reach tens to hundreds of kgalmm . We also observed an improved agreement with normal heights obtained by differencing GNSS and the European gravimetric quasi-geoid 2015 (EGG2015) heights. In Great Britain, the south-north slope of 48 mm deg ^{- 1} present in the solution which uses only spirit leveling/gravimetry data reduced to 2.2 mm deg ^{- 1} . In France, the improvement is confined to the southwest. The choice of the period over which water levels are averaged has an impact on the results as it determines, among others, the set of tide gauges available to establish the hydrodynamic leveling connections. When using an averaging period that can be considered as the least preferred choice based on three established criteria, the positive impact for France has gone. For Great Britain, the estimated south-north slope became 12.6 mm deg ^{- 1} . This is larger than the slope obtained using the most preferred averaging period but still substantially lower compared to the slope associated with a solution that uses only spirit leveling/gravimetry.

Intertidal wetlands, such as mangroves in the tropics, are increasingly recognized for their role in nature-based mitigation of coastal flood risks. Yet it is still poorly understood how effective they are at attenuating the propagation of extreme sea levels through large (order of 100 km2) estuarine or deltaic systems, with complex geometry formed by networks of branching channels intertwined with mangrove and intertidal flat areas. Here, we present a delta-scale hydrodynamic modelling study, aiming to explicitly account for these complex landforms, for the case of the Guayas delta (Ecuador), the largest estuarine system on the Pacific coast of Latin America. Despite coping with data scarcity, our model accurately reproduces the observed propagation of high water levels during a spring tide. Further, based on a model sensitivity analysis, we show that high water levels are most sensitive to the mangrove platform elevation and degree of channelization but to a much lesser extent to vegetation-induced friction. Mangroves with a lower surface elevation, lower vegetation density, and higher degree of channelization all favour a more efficient flooding of the mangroves and therefore more effectively attenuate the high water levels in the deltaic channels. Our findings indicate that vast areas of channelized mangrove forests, rather than densely vegetated forests, are most effective for nature-based flood risk mitigation in a river delta.

Our awareness of ice caps' and mountain glaciers' sensitivity to climate change has driven major advances in the application of remote sensing techniques during the past decade. Regarding ESA's SARIn altimeter CryoSat-2, processing the full waveform to generate swaths of elevation estimates has become standard practice in regions of complex topographies. This technique provides information on areas where we would be blind otherwise. In this article, we discuss systematic errors and analyze their impact on surface elevation measurements and change rates of two test areas. In particular, we focus on periodically occurring errors in elevation swaths, caused by the superposition of coherent signals from range-ambiguous surfaces. They can lead to measurement errors in excess of 10 m, affect most measurements in mountainous regions, are difficult to exclude with established post-processing techniques, and occur repeatedly for satellite revisits introducing a 369-day periodicity-difficult to distinguish from the annual cycle. We show a correlation between derived elevation swaths and the sensor view angle and explore the influence of common data exclusion choices on higher level products. Our results indicate that these systematic errors hold a substantial share of the error budget and that the choice of thresholds impacts higher level products. We conclude that error correlations need to be considered to characterize the data accuracy. With the established data editing strategies, systematic errors prevent resolving seasonal mass changes of single mountain glacier basins and impact aggregates over larger areas or longer periods.

Both empirical and assimilative global ocean tidal models are significantly more accurate in the deep ocean than in shelf and coastal waters. In this study, we answered whether this is due to the quality of the models used to reduce tide and surge or the general approach to treat tide and surge as two separate components of the water level obtained from stand-alone models, which ignores the nonlinear tide–surge interaction. In doing so, we used tide gauge observations as partially synthetic altimeter time series, tide–surge water-level time series obtained with the 2D Dutch Continental Shelf Model–Flexible Mesh (DCSM), and tide and surge water-level time series obtained using the DCSM, FES2014 (FES) and the Dynamic Atmospheric Correction (DAC) product. Expressed in the root-sum-square (RSS) of the eight main tidal constituents, we obtained a reduction (Formula presented.) % when removing the DCSM tide–surge water levels compared to when we removed the sum of the DCSM tide and DCSM surge water levels. The RSS obtained in the latter case was only 3.3% lower than with FES and DAC. We conclude that the lower tidal estimates accuracy in shelf-coastal waters derives from the missing nonlinear tide–surge interactions.

Estimating the three geophysical variables significant wave height (SWH), sea surface height, and wind speed from satellite altimetry continues to be challenging in the coastal zone because the received radar echoes exhibit significant interference from strongly reflective targets such as sandbanks, sheltered bays, ships etc. Fully focused SAR (FF-SAR) processing exhibits a theoretical along-track resolution of up to less than half a metre. This suggests that the application of FF-SAR altimetry might give potential gains over unfocused SAR (UF-SAR) altimetry to resolve and mitigate small-scale interferers in the along-track direction to improve the accuracy and precision of the geophysical estimates. The objective of this study is to assess the applicability of FF-SAR-processed Sentinel-6 Michael Freilich (S6-MF) coastal altimetry data to obtain SWH estimates as close as possible to the coast. We have developed a multi-mission FF-SAR processor and applied the coastal retracking algorithm CORALv2 to estimate SWH. We assess different FF-SAR and UF-SAR processing configurations, as well as the baseline Level-2 product from EUMETSAT, by comparison with the coastal, high-resolution SWAN-Kuststrook wave model from the Deltares RWsOS North Sea operational forecasting system. This includes the evaluation of the correlation, the median offset, and the percentage of cycles with high correlation as a function of distance to the nearest coastline. Moreover, we analyse the number of valid records and the L2 noise of the records. The case study comprises five coastal crossings of S6-MF that are located along the Dutch coast and the German coast along the East Frisian Islands in the North Sea. We observe that accurate and precise SWH records can be estimated in the nearshore zone within 1–3 km from the coast using satellite SAR altimetry. We find that the FF-SAR-processed dataset with a Level-1b posting rate of 140 Hz shows the greatest similarity with the wave model. We achieve a correlation of ∼0.8 at 80% of valid records and a gain in precision of up to 29% of FF-SAR vs UF-SAR for 1–3 km from the coast. FF-SAR shows, for all cycles, a high correlation of greater than or equal to 0.8 for 1–3 km from the coast. We estimate the decay of SWH from offshore at 30 km to up to 1 km from the coast to amount to 26.4% ± 3.1%.

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.

The main objective of this study is to develop and analyze an empirical noise model for model-derived coastal summer mean water levels (SMWLs) and use that to obtain a more realistic quality impact of combining hydrodynamic leveling and Unified European Leveling Network (UELN) data in realizing the European Vertical Reference System (EVRS). We considered three state-of-the-art hydrodynamic models for the Northeast Atlantic Ocean, including the North Sea and Wadden Sea; AMM7, DCSMv6-ZUNOv4, and 3D DCSM-FM. Moreover, we assess the spatiotemporal performance of these three models in representing coastal SMWLs. The empirical noise models are determined from the differences between observation- and model-derived SMWLs at coastal tide gauges. All three noise models show that the model noise is indeed correlated over sea distances up to hundreds of kilometers. At the same time, they all show a relatively large discontinuity at the origin (i.e., nugget effect); between 12.1 cm² (3D DCSM-FM) and 16.3 cm² (DCSMv6-ZUNOv4). The variance (i.e., covariance at zero sea distance) for these two models is 15.3 cm² and 21.7 cm², respectively. Averaging the water levels over three summers, lowered the variance and nugget effect for 3D DCSM-FM to 12.7 cm² and 10.0 cm², respectively. Our analysis also showed that between 30 and 50% of the variance has to be attributed to errors in the vertical referencing of the tide gauges. We lacked the information to assess what proportion of the observed noise covariances should be attributed to these errors. The performance assessments revealed significant variations over both space and time as well as among the three hydrodynamic models. The results suggest that there is still room for model improvement. In the final experiments, we used the noise model of the best overall performing model (i.e., 3D DCSM-FM) to reassess the quality impact of combining hydrodynamic leveling and UELN data in realizing the EVRS. The results suggest that not including the noise covariances leads to an overestimation of the total quality impact by 7 % and 8 % , when we average the water levels over one and three summer periods, respectively.

We present exact, closed-form expressions for the complete RTM correction and the harmonic correction to disturbing potential, gravity disturbance, gravity anomaly, and height anomaly. They need to be applied in quasi-geoid modelling whenever data points are buried inside the masses after residual terrain model (RTM) reduction and analytically downward-continued functionals of the disturbing potential at the original locations of the data points are required. Compared to recent work of the authors published in this journal, no Taylor series enter the expressions and numerical instabilities of the harmonic downward continuation from the RTM surface to the Earth’s surface are avoided as are inaccuracies in the free-air upward continuation from the Earth’s surface to the RTM surface caused by a lack of precise information about higher-order derivatives of the disturbing potential. The new expressions can easily be implemented in any existing RTM software package and do not require additional computational resources. For two test areas located in western Norway and the Auvergne in France, we compute the complete RTM correction and the harmonic correction to the afore-mentioned functionals of the disturbing potential. Overall, all harmonic corrections are non-negative with maximum values of 1.54 m ²/ s ², 263.0 μ Gal, 263.9 μ Gal, and 15.7 m (Norway) and 1.55 m ²/ s ², 263.3 μ Gal, 263.3 μ Gal, and 15.8 cm (Auvergne) for disturbing potential, gravity disturbance, gravity anomaly, and height anomaly, respectively. The medians are 0.02 m ²/ s ², 33.6 μ Gal, 33.7 μ Gal, and 0.3 cm (Norway) and 0.01 m ²/ s ², 19.2 μ Gal, 19.2 μ Gal, and 0.1 cm (Auvergne). We also show that the nth Taylor polynomials used in the recent work of the authors published in this journal may have large remainders depending on the topography in the vicinity of the evaluation point no matter how n is chosen. Finally, we show that the commonly used expression for the harmonic correction to gravity anomaly introduced in 1981 is almost exact, though it was derived along a completely different line of reasoning. The errors do not exceed 49 μ Gal in both test areas. Moreover, the errors have a negligible impact on the computed height anomalies in one-centimetre quasi-geoid modelling, as the mean error does not exceed a few μ Gal in both test areas.

In this paper, we derive improved expressions for the harmonic correction to gravity and, for the first time, expressions for the harmonic correction to potential and height anomaly. They need to be applied at stations buried inside the masses to transform internal values into harmonically downward continued values, which are then input to local quasi-geoid modelling using least-squares collocation or least-squares techniques in combination with the remove-compute-restore approach. Harmonic corrections to potential and height anomaly were assumed to be negligible so far resulting in yet unknown quasi-geoid model errors. The improved expressions for the harmonic correction to gravity, and the new expressions for the harmonic correction to potential and height anomaly are used to quantify the approximation errors of the commonly used harmonic correction to gravity and to quantify the magnitude of the harmonic correction to potential and height anomaly. This is done for two test areas with different topographic regimes. One comprises parts of Norway and the North Atlantic where the presence of deep, long, and narrow fjords suggest extreme values for the harmonic correction to potential and height anomaly and corresponding large errors of the commonly used approximation of the harmonic correction to gravity. The other one is located in the Auvergne test area with a moderate topography comprising both flat and hilly areas and therefore may be representative for many areas around the world. For both test areas, two RTM surfaces with different smoothness are computed simulating the use of a medium-resolution and an ultra-high-resolution reference gravity field, respectively. We show that the errors of the commonly used harmonic correction to gravity may be as large as the harmonic correction itself and attain peak values in areas of strong topographic variations of about 100 mGal. Moreover, we show that this correction may introduce long-wavelength biases in the computed quasi-geoid model. Furthermore, we show that the harmonic correction to height anomaly can attain values on the order of a decimetre at some points. Overall, however, the harmonic correction to height anomaly needs to be applied only in areas of strong topographic variations. In flat or hilly areas, it is mostly smaller than one centimetre. Finally, we show that the harmonic corrections increase with increasing smoothness of the RTM surface, which suggests to use a RTM surface with a spatial resolution comparable to the finest scales which can be resolved by the data rather than depending on the resolution of the global geopotential model used to reduce the data.

Satellite radar altimetry has been an important tool for cryospheric applications such as measuring ice-sheet height or assessing anomalies in snow and ice properties (e.g. the extensive melt in Greenland in 2012). Although accurate height measurements are key for such applications, slope-induced errors due to undulating topography within the kilometre-wide beam-limited footprint can cause multi-metre errors. Two main correction methods that have been developed (referred to as the slope- and point-based methods) neglect either the actual topography or the actual footprint that can be estimated by a combination of the leading edge and topography. Therefore, a leading edge point-based (LEPTA) method is presented that corrects for the slope-induced error by including the leading edge information of the radar waveform to determine the impact point. The principle of the method is that only the points on the ground that are within the range determined by the beginning and end of the leading edge are used to determine the impact point. Benchmarking of the LEPTA method against the slope- and point-based methods based on CryoSat-2 Low Resolution Mode (LRM) acquisitions over Greenland in 2019 shows that, when compared to ICESat-2 observations, the LEPTA method has a stable performance both in the flat, interior regions of Greenland and in regions with more complex topography. The median difference between the slope-corrected CryoSat-2 heights using LEPTA and the ICESat-2 heights is at the millimetre level, whereas the slope and point-based methods can have a 0.21 and 0.48 m difference, respectively, and the Level-2I (L2I) data provided by ESA have a 0.01 m difference. The median absolute deviation of height differences between CryoSat-2 and ICESat-2, which we use as an indicator of the variation in errors, is also the lowest for LEPTA (0.09 m) in comparison to the aforementioned methods (0.19 m for slope method and 0.10 m for point-based method) and ESA Level-2 data (0.14 m). Although ESA Level-2 products and the point-based method have good performance in either the median or the median absolute deviation, LEPTA shows a good performance in both metrics. Based on that, we recommend considering LEPTA for obtaining accurate height measurements with radar altimetry data, especially towards the margins of the LRM coverage where the surface slopes increase.