Ice losses from the Greenland and Antarctic ice sheets have accelerated since the 1990s, accounting for a significant increase in the global mean sea level. Here, we present a new 29-year record of ice sheet mass balance from 1992 to 2020 from the Ice Sheet Mass Balance Inter-comparison Exercise (IMBIE). We compare and combine 50 independent estimates of ice sheet mass balance derived from satellite observations of temporal changes in ice sheet flow, in ice sheet volume, and in Earth's gravity field. Between 1992 and 2020, the ice sheets contributed 21.0±1.9 mm to global mean sea level, with the rate of mass loss rising from 105 Gt yr ⁻¹ between 1992 and 1996 to 372 Gt yr ⁻¹ between 2016 and 2020. In Greenland, the rate of mass loss is 169±9 Gt yr ⁻¹ between 1992 and 2020, but there are large inter-annual variations in mass balance, with mass loss ranging from 86 Gt yr ⁻¹ in 2017 to 444 Gt yr ⁻¹ in 2019 due to large variability in surface mass balance. In Antarctica, ice losses continue to be dominated by mass loss from West Antarctica (82±9 Gt yr ⁻¹ ) and, to a lesser extent, from the Antarctic Peninsula (13±5 Gt yr ⁻¹ ). East Antarctica remains close to a state of balance, with a small gain of 3±15 Gt yr ⁻¹ , but is the most uncertain component of Antarctica's mass balance. The dataset is publicly available at https://doi.org/10.5285/77B64C55-7166-4A06-9DEF-2E400398E452 (IMBIE Team, 2021).

The Gravity Recovery And Climate Experiment (GRACE) satellite mission recorded temporal variations in the Earth’s gravity field, which are then converted to Total Water Storage Change (TWSC) fields representing an anomaly in the water mass stored in all three physical states, on and below the surface of the Earth. GRACE provided a first global observational record of water mass redistribution at spatial scales greater than 63000 km². This limits their usability in regional hydrological applications. In this study, we implement a statistical downscaling approach that assimilates 0.5° × 0.5° water storage fields from the WaterGAP hydrology model (WGHM), precipitation fields from 3 models, evapotranspiration and runoff from 2 models, with GRACE data to obtain TWSC at a 0.5° × 0.5° grid. The downscaled product exploits dominant common statistical modes between all the hydrological datasets to improve the spatial resolution of GRACE. We also provide open access to scripts that researchers can use to produce downscaled TWSC fields with input observations and models of their own choice.

The water budget equation describes the exchange of water between the land, ocean, and atmosphere. Being able to adequately close the water budget gives confidence in our ability to model and/or observe the spatio-temporal variations in the water cycle and its components. Due to advances in observation techniques, satellite sensors, and modelling, a number of data products are available that represent the components of water budget in both space and time. Despite these advances, closure of the water budget at the global scale has been elusive. In this study, we attempt to close the global water budget using precipitation, evapotranspiration, and runoff data at the catchment scale. The large number of recent state-of-the-art datasets provides a new evaluation of well-used datasets. These estimates are compared to terrestrial water storage (TWS) changes as measured by the Gravity Recovery And Climate Experiment (GRACE) satellite mission. We investigated 189 river basins covering more than 90 % of the continental land area. TWS changes derived from the water balance equation were compared against GRACE data using two metrics: the Nash–Sutcliffe efficiency (NSE) and the cyclostationary NSE. These metrics were used to assess the performance of more than 1600 combinations of the various datasets considered. We found a positive NSE and cyclostationary NSE in 99 % and 62 % of the basins examined respectively. This means that TWS changes reconstructed from the water balance equation were more accurate than the long-term (NSE) and monthly (cyclostationary NSE) mean of GRACE time series in the corresponding basins. By analysing different combinations of the datasets that make up the water balance, we identified data products that performed well in certain regions based on, for example, climatic zone. We identified that some of the good results were obtained due to the cancellation of errors in poor estimates of water budget components. Therefore, we used coefficients of variation to determine the relative quality of a data product, which helped us to identify bad combinations giving us good results. In general, water budget components from ERA5-Land and the Catchment Land Surface Model (CLSM) performed better than other products for most climatic zones. Conversely, the latest version of CLSM, v2.2, performed poorly for evapotranspiration in snow-dominated catchments compared, for example, with its predecessor and other datasets available. Thus, the nature of the catchment dynamics and balance between components affects the optimum combination of datasets. For regional studies, the combination of datasets that provides the most realistic TWS for a basin will depend on its climatic conditions and factors that cannot be determined a priori. We believe that the results of this study provide a road map for studying the water budget at catchment scale.

The mass change information from the Gravity Recovery And Climate Experiment (grace) satellite mission is available in terms of noisy spherical harmonic coefficients truncated at a maximum degree (band-limited). Therefore, filtering is an inevitable step in post-processing of grace fields to extract meaningful information about mass redistribution in the Earth-system. It is well known from previous studies that a number can be allotted to the spatial resolution of a band-limited spherical harmonic spectrum and also to a filtered field. Furthermore, it is now a common practice to correct the filtered grace data for signal damage due to filtering (or convolution in the spatial domain). These correction methods resemble deconvolution, and, therefore, the spatial resolution of the corrected grace data have to be reconsidered. Therefore, the effective spatial resolution at which we can obtain mass changes from grace products is an area of debate. In this contribution, we assess the spatial resolution both theoretically and practically. We confirm that, theoretically, the smallest resolvable catchment is directly related to the band-limit of the spherical harmonic spectrum of the grace data. However, due to the approximate nature of the correction schemes and the noise present in grace data, practically, the complete band-limited signal cannot be retrieved. In this context, we perform a closed-loop simulation comparing four popular correction schemes over 255 catchments to demarcate the minimum size of the catchment whose signal can be efficiently recovered by the correction schemes. We show that the amount of closure error is inversely related to the size of the catchment area. We use this trade-off between the error and the catchment size for defining the potential spatial resolution of the grace product obtained from a correction method. The magnitude of the error and hence the spatial resolution are both dependent on the correction scheme. Currently, a catchment of the size ≈63,000 km ${}^2$ can be resolved at an error level of $2\mathrm{cm}$ in terms of equivalent water height.