Measurement-based perturbation theory and differential equation parameter estimation with applications to satellite gravimetry

TL;DR

This paper introduces a measurement-based perturbation theory for satellite gravimetry, providing globally convergent solutions to nonlinear differential equations, improving the accuracy of gravitational models from satellite data.

Contribution

It challenges the traditional zero initial condition approach, proposing a new method that yields globally convergent solutions for high-precision gravitational modeling.

Findings

Proves zero initial derivatives are mathematically and physically invalid.

Derives new local solutions to satellite motion equations.

Demonstrates global convergence enables extraction of minute gravitational signals.

Abstract

The numerical integration method has been routinely used to produce global standard gravitational models from satellite tracking measurements of CHAMP/GRACE types. It is implemented by solving the differential equations of the partial derivatives of a satellite orbit with respect to the unknown harmonic coefficients under the conditions of zero initial values. From the mathematical point of view, satellite gravimetry from satellite tracking is the problem of estimating unknown parameters in the Newton's nonlinear differential equations from satellite tracking measurements. We prove that zero initial values for the partial derivatives are incorrect mathematically and not permitted physically. The numerical integration method, as currently implemented and used in satellite gravimetry and statistics, is groundless. We use three different methods to derive new local solutions to the Newton's nonlinear governing differential equations of motion with a nominal reference orbit. Bearing in mind that satellite orbits can now be tracked almost continuously at unprecedented high accuracy, we propose the measurement-based perturbation theory and derive global uniformly convergent solutions to the Newton's nonlinear governing differential equations of motion. Since the solutions are global uniformly convergent, they are able to extract smallest possible gravitational signals from modern and future satellite tracking measurements for global high-precision, high-resolution gravitational models. By directly turning the nonlinear differential equations of satellite motion into the nonlinear integral equations, we reformulate the links between satellite tracking measurements and the global uniformly convergent solutions to the Newton's governing differential equations as a condition adjustment model with equality constraints.