Abstract:The numerical integration method has been routinely used by major institutions worldwide, for example, NASA Goddard Space Flight Center and German Research Center for Geosciences (GFZ), to produce global gravitational models from satellite tracking measurements of CHAMP and/or GRACE types. Such Earth’s gravitational products have found widest possible multidisciplinary applications in Earth Sciences. The method is essentially implemented by solving the differential equations of the partial derivatives of the orbit of a satellite with respect to the unknown harmonic coefficients under the conditions of zero initial values. From the mathematical and statistical point of view, satellite gravimetry from satellite tracking is essentially 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 mathematics and statistics, chemistry and physics, and satellite gravimetry, is groundless, mathematically and physically. Given the Newton’s nonlinear governing differential equations of satellite motion with unknown equation parameters and unknown initial conditions, we develop three methods to derive new local solutions around a nominal reference orbit, which are linked to measurements to estimate the unknown corrections to approximate values of the unknown parameters and the unknown initial conditions. Bearing in mind that satellite orbits can now be tracked almost continuously at unprecedented accuracy, we propose the measurement-based perturbation theory and derive global uniformly convergent solutions to the Newton’s nonlinear governing differential equations of satellite motion for the next generation of global gravitational models. Since the solutions are global uniformly convergent, theoretically speaking, they are able to extract smallest possible gravitational signals from modern and future satellite tracking measurements, leading to the production of global high-precision, high-resolution gravitational models. By directly turning the nonlinear differential equations of satellite motion into the nonlinear integral equations, and recognizing the fact that satellite orbits are measured with random errors, we further 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 unknown parameters, or equivalently, the weighted least squares estimation of unknown differential equation parameters with equality constraints, for the reconstruction of global high-precision, high-resolution gravitational models from modern (and future) satellite tracking measurements.

Effect of Gravitational Perturbations on Satellite Orbit Via Different Astronomy Standards

Study on Secular Effect Caused by Earth Non-Spherical and Compensation for Its Effect

Analytical Propagation on Relative Motion of Geostationary Earth Orbit Satellites with Orbital Perturbations

Satellite motion in a non-singular gravitational potential

Mathematical Modelling and a Numerical Solution for High Precision Satellite Ephemeris Determination

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

On Precise Inertial Force Modeling for Autonomous Orbit Propagation in Earth-Centered Fixed System for Earth Satellites

Extending higher-order model for non-conservative perturbing forces acting on Galileo satellites during eclipse periods

Orbit Propagation in Irregular and Uncertain Gravity Field Using Differential Algebra

A Parameterized Model to Achieve Orbit Extrapolation On-Board for Highly Elliptical Orbits

The determination of gravitational potential differences from satellite-to-satellite tracking

Analytical Solutions for Satellite Relative Motion under J2 Perturbation

Autonomous orbit determination using epoch-differenced gravity gradients and starlight refraction

Dynamic Models Of Satellite Relative Motion Around An Oblate Earth

Useful Relative Motion Description Method for Perturbations Analysis in Satellite Formation Flying

High-Precision, Practical Relative Orbit Propagation Using Least-Squares Method

A Method for Measuring Gravitational Potential of Satellite's Orbit Using Frequency Signal Transfer Technique between Satellites

Analysis of Orbit Accuracy for Non-cooperative Earth-Orbiting Objects

Propagation and Estimation of the Dynamical Behaviour of Gravitationally Interacting Rigid Bodies

Influence of the Low-Frequency Error of the Residual Orbit on Recovering Time-Variable Gravity Field from the Satellite-To-Satellite Tracking Mission

Numerical modeling of a Global Navigation Satellite System in a general relativistic framework