Layer-Based Modelling of the Earth’s Gravitational Potential up to 10-km Scale in Spherical Harmonics in Spherical and Ellipsoidal Approximation

Global forward modelling of the Earth’s gravitational potential, a classical problem in geophysics and geodesy, is relevant for a range of applications such as gravity interpretation, isostatic hypothesis testing or combined gravity field modelling with high and ultra-high resolution. This study presents spectral forward modelling with volumetric mass layers to degree 2190 for the first time based on two different levels of approximation. In spherical approximation, the mass layers are referred to a sphere, yielding the spherical topographic potential. In ellipsoidal approximation where an ellipsoid of revolution provides the reference, the ellipsoidal topographic potential (ETP) is obtained. For both types of approximation, we derive a mass layer concept and study it with layered data from the Earth2014 topography model at 5-arc-min resolution. We show that the layer concept can be applied with either actual layer density or density contrasts w.r.t. a reference density, without discernible differences in the computed gravity functionals. To avoid aliasing and truncation errors, we carefully account for increased sampling requirements due to the exponentiation of the boundary functions and consider all numerically relevant terms of the involved binominal series expansions. The main outcome of our work is a set of new spectral models of the Earth’s topographic potential relying on mass layer modelling in spherical and in ellipsoidal approximation. We compare both levels of approximations geometrically, spectrally and numerically and quantify the benefits over the frequently used rock-equivalent topography (RET) method. We show that by using the ETP it is possible to avoid any displacement of masses and quantify also the benefit of mapping-free modelling. The layer-based forward modelling is corroborated by GOCE satellite gradiometry, by in-situ gravity observations from recently released Antarctic gravity anomaly grids and degree correlations with spectral models of the Earth’s observed geopotential. As the main conclusion of this work, the mass layer approach allows more accurate modelling of the topographic potential because it avoids 10–20-mGal approximation errors associated with RET techniques. The spherical approximation is suited for a range of geophysical applications, while the ellipsoidal approximation is preferable for applications requiring high accuracy or high resolution.