The Gravity Recovery and Climate Experiment (GRACE) satellite mission, which was in operation from March 2002 to June 2017, was the first remote sensing mission to provide temporal variations of Terrestrial Water Storage (TWS), which is the sum of the water masses that were contained in the soil column (i.e., snow, surface water, soil moisture, and groundwater), at a spatial resolution of a few hundred kilometers. As in situ level measurements are generally not sufficiently available for monitoring groundwater changes at the regional-scale, this unique dataset, combined with external information, is widely used to quantify the interannual variations of groundwater storage in the world’s major aquifers. GRACE-based groundwater changes revealed significant aquifer depletion over large regions, such as the Middle East, the northwest India aquifer, the North China Plain aquifer, the Murray-Darling Basin in Australia, the High Plains, and the California Central Valley aquifers in the United States of America (USA), but were also used to estimate groundwater-related parameters such as the specific yield, which relates groundwater level to storage, or to define the indices of groundwater depletion and stress. In this review, the approaches used for estimating groundwater storage variations are presented along with the main applications of GRACE data for groundwater monitoring. Issues that were related to the use of GRACE-based TWS are also addressed.

Monitoring changes in coastal sea levels is necessary given the impacts of climate change. Information on the sea level and its changes are important parameters in connection to climate change processes. In this study, radar altimetry data from successive satellite missions, European Remote Sensing-2 (ERS-2), Jason-1, Envisat, Jason-2, and Satellite with ARgos and ALtiKa (SARAL), were used to measure sea surface heights (SSH). Altimetry-derived SSH was validated for the southern Bay of Biscay, using records from seven tide gauges located along the French Atlantic coast. More detailed comparisons were performed at La Rochelle, as this was the only tide gauge whose records covered the entire observation period for the different radar altimetry missions. The results of the comparison between the altimetry-based and in-situ SSH, recorded from zero to five kilometers away from the coast, had root mean square errors (RMSE) ranging from 0.08 m to 0.21 m, 0.17 m to 0.34 m, 0.1 m to 0.29 m, 0.18 m to 0.9 m, and 0.22 m to 0.89 m for SARAL, Jason-2, Jason-1, ENVISAT, and ERS-2, respectively. Comparing the missions on the same orbit, ENVISAT had better results than ERS-2, which can be accounted for by the improvements in the sensor mode of operation, whereas the better results obtained using SARAL are related to the first-time use of the Ka-band for an altimetry sensor. For Jason-1 and Jason-2, improvements were found in the ocean retracking algorithm (MLE-4 against MLE-3), and also in the bi-frequency ionosphere and radiometer wet troposphere corrections. Close to the shore, the use of model-based ionosphere (GIM) and wet troposphere (ECMWF) corrections, as applied to land surfaces, reduced the error on the SSH estimates.

In this study, three months of records (January–March 2010) that were acquired by a geodetic Global Navigation Satellite Systems (GNSS) station from the permanent network of RGP (Réseau GNSS Permanent), which was deployed by the French Geographic Institute (IGNF), located in Socoa, in the south of the Bay of Biscay, were used to determine the tide components and identify the signature of storms on the signal to noise ratio (SNR) during winter 2010. The Xynthia storm hit the French Atlantic coast on the 28th of February 2010, causing large floods and damages from the Gironde to the Loire estuaries. Blind separation of the tide components and of the storm signature was achieved while using both a singular spectrum analysis (SSA) and a continuous wavelet transform (CWT). A correlation of 0.98/0.97 and root mean square error (RMSE) of 0.21/0.28 m between the tide gauge records of Socoa and our estimates of the sea surface height (SSH) using the SSA and the CWT, respectively, were found. Correlations of 0.76 and 0.7 were also obtained between one of the modes from the SSA and atmospheric pressure from a meteorological station and a mode of the SSA. Particularly, a correlation reaches to 0.76 when using both the tide residual that is associated to surges and atmospheric pressure variation.

GRACE spherical harmonics are well-adapted for representation of hydrological signals in river drainage basins of large size such as the Amazon or Mississippi basins. However, when one needs to study smaller drainage basins, one comes up against the low spatial resolution of the solutions in spherical harmonics. To overcome this limitation, we propose a new approach based on Slepian functions which can reduce the energy loss by integrating information in the spatial, spectral and time domains. Another advantage of these regionally-defined functions is the reduction of the problem dimensions compared to the spherical harmonic parameters. This also induces a drastic reduction of the computational time. These Slepian functions are used to invert the GRACE satellite data to restore the water mass fluxes of different hydro-climatologic environments in Africa. We apply them to two African drainage basins chosen for their size of medium scale and their geometric specificities: the Congo river basin with a quasi-isotropic shape and the Nile river basin with an anisotropic and more complex shape. Time series of Slepian coefficients have been estimated from real along-track GRACE geopotential differences for about ten years, and these coefficients are in agreement with both the spherical harmonic solutions provided by the official centers CSR, GFZ, JPL and the GLDAS model used for validation. The Slepian function analysis highlights the water mass variations at sub-basin scales in both basins.

We demonstrate a new approach to recover water mass changes from GRACE satellite data at a daily temporal resolution. Such a product can be beneficial in monitoring extreme weather events that last a few days and are missing by conventional monthly GRACE data. The determination of the distribution of these water mass sources over networks of juxtaposed triangular tiles was made using Kalman Filtering (KF) of daily GRACE geopotential difference observations that were reduced for isolating the continental hydrology contribution of the measured gravity field. Geopotential differences were obtained from the along-track K-Band Range Rate (KBRR) measurements according to the method of energy integral. The recovery approach was validated by inverting synthetic GRACE geopotential differences simulated using GLDAS/WGHM global hydrology model outputs. Series of daily regional and global KF solutions were estimated from real GRACE KBRR data for the period 2003–2012. They provide a realistic description of hydrological fluxes at monthly time scales, which are consistent with classical spherical harmonics and mascons solutions provided by the GRACE official centers but also give an intra-month/daily continuity of these variations.

In this study, we propose to estimate the steric sea-level variations over a < 2-year period (April 2002 through December 2003) by combining global mean sea level (GMSL) based on Topex/Poseidon (T/P) altimetry with time-variable geoid averaged over the oceans, as observed by the GRACE (Gravity Recovery and Climate Experiment) satellite. In effect, altimetry-derived GMSL changes results from two contributions: Steric (thermal plus salinity) effects due to sea water density change and ocean mass change due to water exchange with atmosphere and continents. On the other hand, GRACE data over the oceans provide the ocean mass change component only. The paper first discusses the corrections to apply to the GRACE data. Then the steric contribution to the GMSL is estimated using GRACE and T/P data. Comparison with available thermal expansion based on in situ hydrographic data is performed.