Abstract.
An algorithm for the determination of the spherical harmonic coefficients of the terrestrial gravitational field representation from the analysis of a kinematic orbit solution of a low earth orbiting GPS-tracked satellite is presented and examined. A gain in accuracy is expected since the kinematic orbit of a LEO satellite can nowadays be determined with very high precision, in the range of a few centimeters. In particular, advantage is taken of Newton's Law of Motion, which balances the acceleration vector with respect to an inertial frame of reference (IRF) and the gradient of the gravitational potential. By means of triple differences, and in particular higher-order differences (seven-point scheme, nine-point scheme), based upon Newton's interpolation formula, the local acceleration vector is estimated from relative GPS position time series. The gradient of the gravitational potential is conventionally given in a body-fixed frame of reference (BRF) where it is nearly time independent or stationary. Accordingly, the gradient of the gravitational potential has to be transformed from spherical BRF to Cartesian IRF. Such a transformation is possible by differentiating the gravitational potential, given as a spherical harmonics series expansion, with respect to Cartesian coordinates by means of the chain rule, and expressing zero- and first-order Ferrer's associated Legendre functions in terms of Cartesian coordinates. Subsequently, the BRF Cartesian coordinates are transformed into IRF Cartesian coordinates by means of the polar motion matrix, the precession–nutation matrices and the Greenwich sidereal time angle (GAST). In such a way a spherical harmonic representation of the terrestrial gravitational field intensity with respect to an IRF is achieved. Numerical tests of a resulting Gauss–Markov model document not only the quality and the high resolution of such a space gravity spectroscopy, but also the problems resulting from noise amplification in the acceleration determination process.
Similar content being viewed by others
Acknowledgments.
The authors are indebted to Dr. B. Richter from the DGFI (Deutsches Geodätisches Forschungsinstitut) in Munich for his advice concerning the implementation of the transformation between the qBRF and the qIRF. Furthermore, the authors thank the FESG (Forschungseinrichtung Satellitengeodäsie der TU München) and the IAPG (Institut für Astronomische und Physikalische Geodäsie) in Munich, in particular Prof. M. Rothacher and D. Svehla for providing an 11-day kinematic CHAMP orbit.
Author information
Authors and Affiliations
Corresponding author
Rights and permissions
About this article
Cite this article
Reubelt, T., Austen, G. & Grafarend, E. 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. Journal of Geodesy 77, 257–278 (2003). https://doi.org/10.1007/s00190-003-0322-9
Received:
Accepted:
Issue Date:
DOI: https://doi.org/10.1007/s00190-003-0322-9