Exploring gravity field determination from orbit perturbations of the European Gravity Mission GOCE

Abstract

A comparison was made between two methods for gravity field recovery from orbit perturbations that can be derived from global positioning system satellite-to-satellite tracking observations of the future European gravity field mission GOCE (Gravity Field and Steady-State Ocean Circulation Explorer). The first method is based on the analytical linear orbit perturbation theory that leads under certain conditions to a block-diagonal normal matrix for the gravity unknowns, significantly reducing the required computation time. The second method makes use of numerical integration to derive the observation equations, leading to a full set of normal equations requiring powerful computer facilities. Simulations were carried out for gravity field recovery experiments up to spherical harmonic degree and order 80 from 10 days of observation. It was found that the first method leads to large approximation errors as soon as the maximum degree surpasses the first resonance orders and great care has to be taken with modeling resonance orbit perturbations, thereby loosing the block-diagonal structure. The second method proved to be successful, provided a proper division of the data period into orbital arcs that are not too long.