Low-degree gravity field coefficients based on inverse GNSS method: insights into hydrological and ice mass change studies

The relative displacements of stations from a global network of Global Navigation Satellite System (GNSS) sites provide information on global mass transport. In this study, we use 19 years of global GNSS station displacements from the 3rd International GNSS Service reprocessing campaign to estimate the coefficients of the spherical harmonics of the Earth's gravity field up to degree and order 8 using the inverse GNSS method based on elastic loading theory. The results indicate that the C30 coefficient can be derived based on GNSS station displacements as an alternative to solutions provided by Satellite Laser Ranging (SLR) and Gravity Recovery and Climate Experiment (GRACE). GNSS may support GRACE solutions that face the problems of deriving C30, which has fundamental meaning in estimating ice mass changes in polar regions. The recovery of Antarctic ice sheet mass change from January 2007 to December 2020 based on coefficients replaced by GNSS estimates results in a linear trend of - 152 +/- 4 Gt/year, in comparison to - 149 +/- 2 Gt/year for the replacement based on SLR from GRACE Technical Note #14. The results indicate that the spatial and seasonal patterns of terrestrial water storage changes derived from GNSS are consistent with those estimated using GRACE/GRACE Follow-On and SLR at a few-millimeter level in the Amazon and Brahmaputra River basin regions.