Frequency-dependent data weighting in global gravity field modeling from satellite data contaminated by non-stationary noise

Satellite data that are used to model the global gravity field of the Earth are typically corrupted by correlated noise, which can be related to a frequency dependence of the data accuracy. We show an opportunity to take such noise into account by using a proper noise covariance matrix in the estimation procedure. If the dependence of noise on frequency is not known a priori, it can be estimated on the basis of a posteriori residuals. The methodology can be applied to data with gaps. Non-stationarity of noise can also be dealt with, provided that the necessary a priori information exists. The proposed methodology is illustrated with CHAllenging Mini-satellite Payload (CHAMP) data processing. It is shown, in particular, that the usage of a proper noise model can make the measurements of non-gravitational satellite accelerations unnecessarily. This opens the door for high-quality modeling of the Earth’s gravity field on the basis of observed orbits of non-dedicated satellites (i.e., satellites without an on-board accelerometer). Furthermore, the processing of data from dedicated satellite missions – GRACE (Gravity Recovery and Climate Experiment) and GOCE (Gravity field and steady-state Ocean Circulation Explorer) – may also benefit from the proposed methodology.

[1]  J. Kusche,et al.  Regularization of geopotential determination from satellite data by variance components , 2002 .

[2]  Pavel Ditmar,et al.  A technique for modeling the Earth’s gravity field on the basis of satellite accelerations , 2004 .

[3]  Jürgen Kusche,et al.  Stochastic model validation of satellite gravity data: A test with CHAMP pseudo-observations , 2005 .

[4]  Pavel Ditmar,et al.  M ODELING THE EARTH ' S GRAVITY FIELD FROM PRECISE SATELLITE ORBIT DATA : THE ACCELERATION APPROACH WORKS ! , 2004 .

[5]  N. K. Pavlis,et al.  The Development of the Joint NASA GSFC and the National Imagery and Mapping Agency (NIMA) Geopotential Model EGM96 , 1998 .

[6]  Victor Zlotnicki,et al.  Time‐variable gravity from GRACE: First results , 2004 .

[7]  Battista Benciolini,et al.  Hotine-Marussi Symposium on Mathematical Geodesy , 1990 .

[8]  M. Rothacher,et al.  Kinematic Precise Orbit Determination for Gravity Field Determination , 2005 .

[9]  J. Bouman.  Quality assessment of satellite-based global gravity field models , 2000 .

[10]  J. Lemoine,et al.  Earth Gravity Field and Seasonal Variability from CHAMP , 2005 .

[11]  P. Teunissen.  Adjustment Theory: an introduction , 2000 .

[12]  Fernando Sansò,et al.  A Window on the Future of Geodesy , 2005 .

[13]  A. N. Tikhonov,et al.  Solutions of ill-posed problems , 1977 .

[14]  S.M. Kay,et al.  Spectrum analysis—A modern perspective , 1981, Proceedings of the IEEE.

[15]  M Krügel.  Dynamic Planet 2005 , 2005 .

[16]  Christopher Jekeli,et al.  Gravity, Geoid and Space Missions , 2008 .

[17]  R. Koop,et al.  Validation of fast pre-mission error analysis of the GOCE gradiometry mission by a full gravity field recovery simulation , 2002 .

[18]  Baltimore County.  The Development of the Joint NASA GSFC and the National Imagery and Mapping Agency (NIMA) Geopotential Model EGM96 , 1998 .

[19]  C. Reigber,et al.  Gravity field recovery from satellite tracking data , 1989 .

[20]  J. Kusche,et al.  Regularization of gravity field estimation from satellite gravity gradients , 2002 .

[21]  Roland Klees,et al.  ‘DEOS_CHAMP-01C_70’: a model of the Earth’s gravity field computed from accelerations of the CHAMP satellite , 2006 .

[22]  M. Watkins,et al.  GRACE Measurements of Mass Variability in the Earth System , 2004, Science.

[23]  R. Klees,et al.  How to handle colored noise in large least-squares problems in the presence of data gaps? , 2004 .

[24]  U. Meyer,et al.  An Earth gravity field model complete to degree and order 150 from GRACE: EIGEN-GRACE02S , 2005 .

[25]  W. Schuh,et al.  The Processing of Band-Limited Measurements; Filtering Techniques in the Least Squares Context and in the Presence of Data Gaps , 2003 .

[26]  Roland Klees,et al.  Computation of spherical harmonic coefficients from gravity gradiometry data to be acquired by the GOCE satellite: regularization issues , 2003 .

[27]  Fernando Sansò,et al.  Space-wise approach to satellite gravity field determination in the presence of coloured noise , 2004 .

[28]  M. Hestenes,et al.  Methods of conjugate gradients for solving linear systems , 1952 .

[29]  P. Broersen,et al.  How to handle colored observation noise in large least-squares problems , 2003 .

[30]  G. Blewitt,et al.  A New Global Mode of Earth Deformation: Seasonal Cycle Detected , 2001, Science.

[31]  A. Eicker,et al.  ITG-CHAMP01: a CHAMP gravity field model from short kinematic arcs over a one-year observation period , 2005 .

[32]  M. Cheng,et al.  GGM02 – An improved Earth gravity field model from GRACE , 2005 .

[33]  C. Readings,et al.  Gravity field and steady-state ocean circulation mission , 1996 .

[34]  Fernando Sansò,et al.  International Summer School of Theoretical Geodesy on “Theory of Satellite Geodesy and Gravity Field Determination” , 1990 .

[35]  E. Grafarend,et al.  Harmonic analysis of the Earth's gravitational field by means of semi-continuous ephemerides of a low Earth orbiting GPS-tracked satellite. Case study: CHAMP , 2003 .