Efficient gravity field recovery from GOCE gravity gradient observations

An efficient algorithm is proposed for gravity field recovery from Gravity Field and Steady-State Ocean Circulation Explorer (GOCE) satellite gravity gradient observations. The mathematical model is formulated in the time domain, which allows the inclusion of realistic observational noise models. The algorithm combines the iterative solution of the normal equations, using a Richardson-type iteration scheme, with the fast computation of the right-hand side of the normal equations in each iteration step by a suitable approximation of the design matrix. The convergence of the iteration is investigated, error estimates are provided, and the unbiasedness of the method is proved. It is also shown that the method does not converge to the solution of the normal equations. The performance of the approach for white noise and coloured noise is demonstrated along a simulated GOCE orbit up to spherical harmonic degree and order 180. The results also indicate that the approximation error may be neglected.