A global vertical datum defined by the conventional geoid potential and the Earth ellipsoid parameters

Abstract

Abstract The geoid, according to the classical Gauss–Listing definition, is, among infinite equipotential surfaces of the Earth’s gravity field, the equipotential surface that in a least squares sense best fits the undisturbed mean sea level. This equipotential surface, except for its zero-degree harmonic, can be characterized using the Earth’s global gravity models (GGM). Although, nowadays, satellite altimetry technique provides the absolute geoid height over oceans that can be used to calibrate the unknown zero-degree harmonic of the gravimetric geoid models, this technique cannot be utilized to estimate the geometric parameters of the mean Earth ellipsoid (MEE). The main objective of this study is to perform a joint estimation of W0, which defines the zero datum of vertical coordinates, and the MEE parameters relying on a new approach and on the newest gravity field, mean sea surface and mean dynamic topography models. As our approach utilizes both satellite altimetry observations and a GGM model, we consider different aspects of the input data to evaluate the sensitivity of our estimations to the input data. Unlike previous studies, our results show that it is not sufficient to use only the satellite-component of a quasi-stationary GGM to estimate W0. In addition, our results confirm a high sensitivity of the applied approach to the altimetry-based geoid heights, i.e., mean sea surface and mean dynamic topography models. Moreover, as W0 should be considered a quasi-stationary parameter, we quantify the effect of time-dependent Earth’s gravity field changes as well as the time-dependent sea level changes on the estimation of W0. Our computations resulted in the geoid potential W0 = 62636848.102 ± 0.004 m2 s−2 and the semi-major and minor axes of the MEE, a = 6378137.678 ± 0.0003 m and b = 6356752.964 ± 0.0005 m, which are 0.678 and 0.650 m larger than those axes of GRS80 reference ellipsoid, respectively. Moreover, a new estimation for the geocentric gravitational constant was obtained as GM = (398600460.55 ± 0.03) × 106 m3 s−2.