Determination of ellipsoidal surface mass change from GRACE time-variable gravity data

Abstract

SUMMARY The problem of determining mass redistribution within the Earth system from time-variable gravity (TVG) data is non-unique. Over seasonal and decadal time-scales, mass redistribution likely takes place on the Earth’s surface. By approximating the Earth’s surface by a sphere, surface mass variation can be uniquely determined from TVG data. Recently, using the improved GRACE TVG data, Li et al. and Ditmar found that such spherical approximation is no longer tenable and suggested practical approaches to accommodate the elliptical shape of the Earth. In this study, we develop a rigorous method of determining surface mass change on the Earth’s reference ellipsoid. We derive a unique one-to-one relationship between ellipsoidal spectra of surface mass and gravitational potential for the ellipsoidal geometry. In conjunction with our ellipsoidal formulation, the linear transformation between spherical and ellipsoidal harmonic coefficients of the geopotential field enables us to determine mass redistribution on the ellipsoid from GRACE TVG data. Using the Release 6 of GRACE TVG data to degree 60, we show that the ellipsoidal approach reconciles surface mass change rate significantly better than the spherical computation by 3–4 cm yr−1, equivalent to 10–15  per cent increase of total signal, in Greenland and West Antarctica. We quantify the spherical approximation error over the polar regions using GRACE Level-2 TVG data as well as mascon solutions, and demonstrate that the systematic error increases linearly with the maximum degree used for the synthesis. The terrestrial water storage computation is less affected by the spherical approximation because of geographic location of major river basins (lower latitude) and signal characteristics. The improvement of TVG data from GRACE and its Follow-On necessitates the ellipsoidal computation, particularly for quantifying mass change in polar regions.