A Monte-Carlo technique for weight estimation in satellite geodesy

Abstract. Precise orbit determination, and satellite-geodetic applications such as gravity field modelling or satellite altimetry, rely on different observation types and groups that have to be processed in a common parameter estimation scheme. Naturally, the choice of the relative weights for these data sets as well as for added prior information is of importance for obtaining reliable estimates of the unknown parameters and their associated covariance matrices. If the observations are predominantly affected by random errors and systematic errors play a minor role, variance component models can be applied. However, most of the methods proposed so far for variance component estimation involve repeated inversion of large matrices, resulting in intensive computations and large storage requirements if more than a few hundred unknowns are to be determined. In addition, these matrices are not necessarily provided as standard output from common geodetic least-squares estimation software. Therefore, a method is proposed which is based on Monte-Carlo estimation of the redundancy contributions of disjunctive observation groups. The method can handle unknown variance components without the need for repeated inversion of matrices. It is computationally simple, numerically stable and easy to implement. Its application is demonstrated in an experiment concerning low-medium-degree gravity field recovery from simulated orbit perturbations of the GOCE mission, and compared in performance with Lerch's method of subset solutions.

[1]  J. R. Lucas,et al.  MINQUE for block diagonal bordered systems such as those encountered in VLBI data analysis , 1998 .

[2]  Francis J. Lerch,et al.  Optimum data weighting and error calibration for estimation of gravitational parameters , 1989 .

[3]  Robert James Purser,et al.  A Bayesian technique for estimating continuously varying statistical parameters of a variational assimilation , 2003 .

[4]  R. Koop,et al.  Exploring gravity field determination from orbit perturbations of the European Gravity Mission GOCE , 2001 .

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

[6]  George H. Ling,et al.  Die Ausgleichungsrechnung nach der Methode der Kleinsten Quadrate , 1907 .

[7]  Karl-Rudolf Koch,et al.  Bayesian Inference with Geodetic Applications , 1990 .

[8]  P. Visser,et al.  GPS-based precise orbit determination of the very low Earth-orbiting gravity mission GOCE , 2000 .

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

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

[11]  M. Şahin,et al.  Variance component estimation applied to satellite laser ranging , 1992 .

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

[13]  Muhammed Sahin,et al.  SLR precision analysis for LAGEOS I and II , 2000 .

[14]  W. N. Venables,et al.  Estimation of Variance Components and Applications. , 1990 .

[15]  N. Crocetto,et al.  Simplified formulae for the BIQUE estimation of variance components in disjunctive observation groups , 2000 .

[16]  Cheng Huang,et al.  Application of variance component estimation to precise orbit determination for ERS-2 , 2001 .

[17]  A. Girard A fast ‘Monte-Carlo cross-validation’ procedure for large least squares problems with noisy data , 1989 .

[18]  M. Hutchinson A stochastic estimator of the trace of the influence matrix for laplacian smoothing splines , 1989 .

[19]  E. C. Pavlis,et al.  A new gravitational model for the earth from satellite tracking data - GEM-T1 , 1988 .

[20]  Karl-Rudolf Koch,et al.  Maximum likelihood estimate of variance components , 1986 .

[21]  N. L. Johnson,et al.  Linear Statistical Inference and Its Applications , 1966 .

[22]  Calyampudi R. Rao,et al.  Linear Statistical Inference and Its Applications. , 1975 .

[23]  Ronald Bremer,et al.  Estimation of variance components and applications , 1988 .

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