Description: Over the last 15 years, the Gravity Recovery and Climate Experiment (GRACE) mission has provided measurements of temporal changes in mass redistribution at and within the Earth that affect polar motion. The newest generation of GRACE temporal models, are evaluated by conversion into the equatorial components of hydrological polar motion excitation and compared with the residuals of observed polar motion excitation derived from geodetic measurements of the pole coordinates. We analyze temporal variations of hydrological excitation series and decompose them into linear trends and seasonal and non-seasonal changes, with a particular focus on the spectral bands with periods of 1000–3000, 450–1000, 100–450, and 60–100 days. Hydrological and reduced geodetic excitation series are also analyzed in four separated time periods which are characterized by different accuracy of GRACE measurements. The level of agreement between hydrological and reduced geodetic excitation depends on the frequency band considered and is highest for interannual changes with periods of 1000–3000 days. We find that the CSR RL06, ITSG 2018 and CNES RL04 GRACE solutions provide the best agreement with reduced geodetic excitation for most of the oscillations investigated.

Description: Different Earth orientation parameter (EOP) time series are publicly available that typically arise from the combination of individual space geodetic technique solutions. The applied processing strategies and choices lead to systematically differing signal and noise characteristics particularly at the shortest periods between 2 and 8 days. We investigate the consequences of typical choices by introducing new experimental EOP solutions obtained from combinations at either normal equation level processed by Deutsches Geodätisches Forschungsinstitut at the Technical University of Munich (DGFI‐TUM) and Federal Agency for Cartography and Geodesy (BKG), or observation level processed by European Space Agency (ESA). All those experiments contribute to an effort initiated by ESA to develop an independent capacity for routine EOP processing and prediction in Europe. Results are benchmarked against geophysical model‐based effective angular momentum functions processed by Helmholtz Centre Potsdam GFZ German Research Centre for Geosciences (ESMGFZ). We find, that a multitechnique combination at normal equation level that explicitly aligns a priori station coordinates to the ITRF2014 frequently outperforms the current International Earth Rotation and Reference Systems Service (IERS) standard solution 14C04. A multi‐Global Navigation Satellite System (GNSS)‐only solution already provides very competitive accuracies for the equatorial components. Quite similar results are also obtained from a short combination at observation level experiment using multi‐GNSS solutions and SLR from Sentinel‐3A and Sentinel‐3B to realize space links. For ΔUT1, however, very long baseline interferometry (VLBI) information is known to be critically important so that experiments combining only GNSS and possibly SLR at observation level perform worse than combinations of all techniques at normal equation level. The low noise floor and smooth spectra obtained from the multi‐GNSS solution nevertheless illustrates the potential of this most rigorous combination approach so that further efforts to include in particular VLBI are strongly recommended.

Description: Based on the latest GFZ release 06 of monthly gravity fields from GRACE satellite mission, area-averaged barystatic sea-level is found to rise by 2.02 mm/a during the period April 2002 until August 2016 in the open ocean with a 1000 km coastal buffer zone when low degree coefficients are properly augmented with information from satellite laser ranging. Alternative spherical harmonics solutions from CSR, JPL and TU Graz reveal rates between 1.94 and 2.08 mm/a, thereby demonstrating that systematic differences among the centers are much reduced in the latest release. The results from the direct integration in the open ocean can be aligned to associated solutions of the sea-level equation when fractional leakage derived from two differently filtered global gravity fields is explicitly considered, leading to a global mean sea-level rise of 1.72 mm/a. This result implies that estimates obtained from a 1000 km coastal buffer zone are biased 0.3 mm/a high due the systematic omission of regions with below-average barystatic sea-level rise in regions close to substantial coastal mass losses induced by the reduced gravitational attraction of the remaining continental ice and water masses.

Description: The German Research Centre for Geosciences GFZ provides non-tidal atmospheric and oceanic loading (NTAOL) displacements for tracking stations located on the Earth’s surface based on ECMWF numerical weather models (NWM) and compatible with GRACE atmosphere and ocean de-aliasing products (AOD1B). We apply both, the dynamic (AOD1B) and the geometric (NTAOL) models, to Precise Orbit Determination (POD) of the altimetry missions ENVISAT, Jason-1, and Jason-2 with Satellite Laser Ranging (SLR) and Doppler Orbitography and Radiopositioning Integrated by Satellite (DORIS) observations. Alternatively to the GFZ displacements, those provided by the International Mass Loading Service (IMLS) based on the MERRA NWM are applied also. The differences in POD are qualified in terms of orbital fits, orbit differences and altimeter cross-over differences. It turns out that the major effect comes from AOD1B, that from NTAOL is considerably smaller, but visible and of benefit for POD. Long-term systematics introduced by the loading models can not be detected. Geographicaly correlated systematics important for altimetry show up, however well below the millimeter. The effect of loading on altimeter cross-over results is even smaller and close to the point where sensitivity is lost.

Description: Temporal variations in the total ocean mass representing the barystatic part of present-day global-mean sea-level rise can be directly inferred from time-series of global gravity fields as provided by the GRACE and GRACE-FO missions. A spatial integration over all ocean regions, however, largely underestimates present-day rates as long as the effects of spatial leakage along the coasts of in particular Antarctica, Greenland, and the various islands of the Canadian Archipelago are not properly considered. Based on the latest release 06 of monthly gravity fields processed at GFZ, we quantify (and subsequently correct) the contribution of spatial leakage to the post-processed mass anomalies of continental water storage and ocean bottom pressure. We find that by utilizing the sea level equation to predict spatially variable ocean mass trends out of the (leakage-corrected) terrrestial mass distributions from GRACE and GRACE-FO consistent results are obtained also from spatial integrations over ocean masks with different coastal buffer zones ranging from 400 to 1000 km. However, the results are critically dependent on coefficients of degree 1, 2 and 3, that are not precisely determined from GRACE data alone and need to be augemented by information from satellite laser ranging. We will particularly discuss the impact of those low-degree harmonics on the secular rates in global barystatic sea-level.

Description: GRACE/GRACE-FO Level-3 product based on COST-G RL01 Level-2B products (Dahle & Murböck, 2020) representing Ocean Bottom Pressure (OBP) variations provided at 1° latitude-longitude grids as defined over ocean areas. The OBP grids are provided in NetCDF format divided into yearly batches. The files each contain seven different variables: 1) 'barslv': gravity-based barystatic sea-level pressure 2) 'std_barslv': gravity-based barystatic sea-level pressure uncertainties 3) 'resobp': gravity-based residual ocean circulation pressure resobp 4) 'std_resobp': gravity-based residual ocean circulation pressure uncertainties 5) 'leakage': apparent gravity-based bottom pressure due to continental leakage 6) 'model_ocean': background-model ocean circulation pressure 7) 'model_atmosphere': background-model atmospheric surface pressure These Level-3 products are visualized at GFZ's web portal GravIS (http://gravis.gfz-potsdam.de). Link to data products: ftp://isdcftp.gfz-potsdam.de/grace/GravIS/COST-G/Level-3/OBP

Description: GRACE/GRACE-FO Level-3 product based on COST-G RL01 Level-2B products (Dahle & Murböck, 2020) representing Terrestrial Water Storage (TWS) anomalies provided at 1° latitude-longitude grids as defined over all continental regions except Greenland and Antarctica. The TWS anomaly grids are provided in NetCDF format divided into yearly batches. The files each contain four different variables: 1) 'tws': gravity-based TWS 2) 'std_tws': gravity-based TWS uncertainties 3) 'leakage': spatial leakage contained in TWS 4) 'model_atmosphere': background model atmospheric mass These Level-3 products are visualized at GFZ's web portal GravIS (http://gravis.gfz-potsdam.de). Link to data products: ftp://isdcftp.gfz-potsdam.de/grace/GravIS/COST-G/Level-3/TWS

Description: Temporal variations in the total ocean mass representing the barystatic part of present-day global mean sea-level rise can be unambiguously inferred from time-series of global gravity fields as provided by the GRACE and GRACE-FO missions. A spatial integration over all ocean regions, however, largely underestimates present-day rates as long as the effects of spatial leakage along the coasts of in particular Antarctica, Greenland, and the various islands of the Canadian Archipelago are not properly considered. Based on the recent release 06 of monthly gravity fields processed at GFZ, we quantify (and subsequently correct) the contribution of spatial leakage to the post-processed mass anomalies of continental water storage and ocean bottom pressure. Utilising the sea level equation allows to predict spatially variable ocean mass trends out of the (leakage-corrected) terrestrial mass distributions from GRACE and GRACE-FO. Consistent results for the global mean barystatic sea-level rise are obtained also from spatial integrations over ocean masks with different coastal buffer zones ranging from 400 to 1000 km, thereby confirming the robustness of our method. Residual month-to-month variations in ocean bottom pressure are indicative for errors in the monthly-mean estimates of the applied de-aliasing model AOD1B RL06 and will be thus contrasted against very recent MPIOM experiments considered for AOD1B RL07. The in this way improved leakage correction will be implemented in future GravIS versions (http://gravis.gfz-potsdam.de).

Description: Quantifying and monitoring terrestrial water storage (TWS) is an essential task for understanding the Earth's hydrosphere cycle, its susceptibility to climate change, and concurrent impacts for ecosystems, agriculture, and water management. Changes in TWS manifest as anomalies in the Earth's gravity field, which are routinely observed from space. However, the complex underlying distribution of water masses in rivers, lakes, or groundwater basins remains elusive. We combine machine learning, numerical modeling, and satellite altimetry to build a downscaling neural network that recovers simulated TWS from synthetic space‐borne gravity observations. A novel constrained training is introduced, allowing the neural network to validate its training progress with independent satellite altimetry records. We show that the neural network can accurately derive the TWS in 2019 after being trained over the years 2003 to 2018. Further, we demonstrate that the constrained neural network can outperform the numerical model in validated regions.