Gravity field anomalies from recent GOCE satellite-based geopotential models and terrestrial gravity data: a comparative study over Saudi Arabia

Global geopotential models (GGMs) provide important information about the Earth gravity field in terms of precise gravity anomaly, gravity disturbance, and geoid height on global and regional scales. The GGMs based on the Gravity Recovery and steady-state Ocean Circulation Explorer (GOCE) satellite gravity mission provided by the International Center for Global Earth Models (ICGEM) differ in their accuracy and spatial resolution. In this paper, we evaluate the free-air gravity anomalies and geoid heights determined from several recent GOCE-based GGMs using the corresponding gravity functions (i.e., gravity anomalies and geoid heights) obtained from the Earth Gravitational Model 2008 (EGM2008) over the Kingdom of Saudi Arabia (KSA). Moreover, free-air gravity anomalies determined from GOCE-based GGMs have also been compared with the corresponding ones obtained from ground-based "terrestrial" gravity measurements over the study area. The latter comparison, using the available 3500 terrestrial gravity observations, applies the spectral enhanced method (SEM). We use the SEM to compensate the missed high-frequency component of geoid heights in GOCE-based GGMs using EGM2008 and a high-resolution digital terrain model based on the Shuttle Radar Topography Mission (SRTM). The statistical comparisons of the computed results in terms of free-air gravity anomalies and geoid heights obtained from GOCE-based GGMs are in general agreement. The estimated standard deviations of differences are 5 mGal (20 cm), 6-7 mGal (2224 cm), and 8 mGal (25.9 cm) at spherical harmonic (SH) degree/order (d/o) 200, 240, and 280, respectively. Free-air gravity anomalies computed from available terrestrial gravity data, using the SEM, have been compared with the corresponding ones estimated from GOCE-based GGMs, and the results indicate an improvement over EGM2008 in the medium wavelengths. Investigation of various releases of GOCE-based GGM, up to d/o 200, reveals that the fourth and fifth releases of GOCE-based GGMs developed with the use of time-wise solution strategies (TIM_R4 and TIM_R5) approximate the gravity field well over the KSA region, since they provide the least standard deviation differences at SH. Thus, the TIM_R5 and TIM_R4 are suggested as reference models for recovering the long wavelength up to SH d/o 200 and 240, respectively, when modeling the gravimetric quasi-geoid over the KSA area.