Satellite altimetry has been designed for sea level monitoring over open ocean areas. However, for some years, this technology has also been used to retrieve water levels from reservoirs, wetlands and in general any inland water body, although the radar altimetry technique has been especially applied to rivers and lakes. In this paper, a new approach for the estimation of inland water level time series is described. It is used for the computation of time series of rivers and lakes available through the web service "Database for Hydrological Time Series over Inland Waters" (DAHITI). The new method is based on an extended outlier rejection and a Kalman filter approach incorporating cross-calibrated multi-mission altimeter data from Envisat, ERS-2, Jason-1, Jason-2, TOPEX/Poseidon, and SARAL/AltiKa, including their uncertainties. The paper presents water level time series for a variety of lakes and rivers in North and South America featuring different characteristics such as shape, lake extent, river width, and data coverage. A comprehensive validation is performed by comparisons with in situ gauge data and results from external inland altimeter databases. The new approach yields rms differences with respect to in situ data between 4 and 36 cm for lakes and 8 and 114 cm for rivers. For most study cases, more accurate height information than from other available altimeter databases can be achieved.

EOT20 is the latest in a series of empirical ocean tide (EOT) models derived using residual tidal analysis of multi-mission satellite altimetry at DGFI-TUM. The amplitudes and phases of 17 tidal constituents are provided on a global 0.125^∘ grid based on empirical analysis of seven satellite altimetry missions and four extended missions. The EOT20 model shows significant improvements compared to the previous iteration of the global model (EOT11a) throughout the ocean, particularly in the coastal and shelf regions, due to the inclusion of more recent satellite altimetry data as well as more missions, the use of the updated FES2014 tidal model as a reference to estimated residual signals, the inclusion of the ALES retracker and improved coastal representation. In the validation of EOT20 using tide gauges and ocean bottom pressure data, these improvements in the model compared to EOT11a are highlighted with the root sum square (RSS) of the eight major tidal constituents improving by ∼ 1.4 cm for the entire global ocean with the major improvement in RSS (∼ 2.2 cm) occurring in the coastal region. Concerning the other global ocean tidal models, EOT20 shows an improvement of ∼ 0.2 cm in RSS compared to the closest model (FES2014) in the global ocean. Variance reduction analysis was conducted comparing the results of EOT20 with FES2014 and EOT11a using the Jason-2, Jason-3 and SARAL satellite altimetry missions. From this analysis, EOT20 showed a variance reduction for all three satellite altimetry missions with the biggest improvement in variance occurring in the coastal region. These significant improvements, particularly in the coastal region, provide encouragement for the use of the EOT20 model as a tidal correction for satellite altimetry in sea-level research. All ocean and load tide data from the model can be freely accessed at https://doi.org/10.17882/79489 (Hart-Davis et al., 2021). The tide gauges from the TICON dataset used in the validation of the tide model, are available at https://doi.org/10.1594/PANGAEA.896587 (Piccioni et al., 2018 a).

Vertical land movements can cause regional relative sea-level changes to differ substantially from climate-driven absolute sea-level changes. Whereas absolute sea level has been accurately monitored by satellite altimetry since 1992, there are limited observations of vertical land motion. Vertical land motion is generally modelled as a linear process, despite some evidence of nonlinear motion associated with tectonic activity, changes in surface loading or groundwater extraction. As a result, the temporal evolution of vertical land motion, and its contribution to projected sea-level rise and its uncertainty, remains unresolved. Here we generate a probabilistic vertical land motion reconstruction from 1995 to 2020 to determine the impact of regional-scale and nonlinear vertical land motion on relative sea-level projections up to 2150. We show that regional variations in projected coastal sea-level changes are equally influenced by vertical land motion and climate-driven processes, with vertical land motion driving relative sea-level changes of up to 50 cm by 2150. Accounting for nonlinear vertical land motion increases the uncertainty in projections by up to 1 m on a regional scale. Our results highlight the uncertainty in future coastal impacts and demonstrate the importance of including nonlinear vertical land motions in sea-level change projections.

The mapping and forecasting of droughts and floods is an important potential field of application of global soil moisture and water storage products from satellites and models. Especially when extremes in near-surface soil moisture propagate into extremes in total water storage, agricultural production and water supply can be severely impacted. This study relates soil moisture from the WaterGAP Global Hydrology Model (WGHM) and the satellite sensors Advanced Microwave Scanning Radiometer—Earth Observing System (AMSR-E) and Advanced Scatterometer (ASCAT) to total water storage variations from the satellite gravity mission GRACE. A particular focus is on destructive hydrological extreme events, as listed in the International Disaster Database EM-DAT. Data sets are analyzed via correlation, time shift, and principal component analyses. The study area is the La Plata Basin in South America. The results indicate that most of the soil moisture anomalies are linked to periods of El Niño and La Niña and associated natural disasters. For the La Plata drought of 2008/2009 and the El Niño flooding of 2009/2010, soil moisture serves as an indicator for the later deficit or surplus in total water storage. These hydrological anomalies were strongest in the southern, central, and eastern parts of the basin, but more than one hundred thousand people were also affected in the northwestern part.

We provide a present-day surface-kinematics model for the Alpine region and surroundings based on a high-level data analysis of about 300 geodetic stations continuously operating over more than 12 years. This model includes a deformation model, a continuous surface-kinematic (velocity) field, and a strain field consistently assessed for the entire Alpine mountain belt. Special care is given to the use of the newest Global Navigation Satellite Systems (GNSS) processing standards to determine high-precision 3-D station coordinates. The coordinate solution refers to the reference frame IGb08, epoch 2010.0. The mean precision of the station positions at the reference epoch is ±1.1 mm in N and E and ±2.3 mm in height. The mean precision of the station velocities is ±0.2 mm a⁻¹ in N and E and ±0.4 mm a⁻¹ in height. The deformation model is derived from the point-wise station velocities using a geodetic least-squares collocation (LSC) approach with empirically determined covariance functions. According to our results, no significant horizontal deformation is detected in the Western Alps, while across the Southern and Eastern Alps the deformation vectors describe a progressive eastward rotation towards Pannonia. This kinematic pattern also makes evident an increasing magnitude of the deformation from 0.1 mm a⁻¹ in the western part of Switzerland up to about 1.3 mm a⁻¹ in the Austrian Alps. The largest shortening is observed along the southern front of the Eastern Alps (in the northern area of the Venetian-Friuli Basin) and in the northern part of the Apennine Peninsula, where rates reach 2 and 3 mm a⁻¹, respectively. The average accuracy of the horizontal deformation model is ±0.2 mm a⁻¹. Regarding the vertical kinematics, our results clearly show an ongoing average uplift rate of 1.8 mm a⁻¹ of the entire mountain chain, with the exception of the southern part of the Western Alps, where no significant uplift (less than 0.5 mm a⁻¹) is detected. The fastest uplift rates (more than 2 mm a⁻¹) occur in the central area of the Western Alps, in the Swiss Alps, and in the Southern Alps in the boundary region between Switzerland, Austria, and Italy. The general uplift observed across the Alpine mountain chain decreases towards the outer regions to stable values between 0.0 and 0.5 mm a⁻¹ and, in some cases, to subsidence like in the Liguro-Provençal and Vienna basins, where vertical rates of −0.8 and −0.3 mm a⁻¹ are observed, respectively. In the surrounding region, three regional subsidence regimes are identified: the Rhône-Bresse Graben with −0.8 mm a⁻¹, the Rhine Graben with −1.3 mm a⁻¹, and the Venetian-Friuli Basin with −1.5 mm a⁻¹. The estimated uncertainty of our vertical motion model across the Alpine mountain belt is about ±0.3 mm a⁻¹. The strain field inferred from the deformation model shows two main contrasting strain regimes: (i) shortening across the south-eastern front of the Alps and the northern part of the Dinarides and (ii) extension in the Apennines. The pattern of the principal strain axes indicates that the compression directions are more or less perpendicular to the thrust belt fronts, reaching maximum values of 20×10⁻⁹ a⁻¹ in the Venetian-Friuli and Po basins. Across the Alpine mountain belt, we observe a slight dilatation regime in the Western Alps, which smoothly changes to a contraction regime in western Austria and southern Germany, reaching maximum shortening values of 6×10⁻⁹ a⁻¹ in north-eastern Austria. The numerical results of this study are available at https://doi.pangaea.de/10.1594/PANGAEA.886889.