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.

Resolution Estimation of Earth Gravity Field Recovery Through the Low-Low Satellite to Satellite Technology by Numerical Simulation

Precision Real-Time Orbit Determination Algorithm and Hardware-in-the-Loop Simulation for ZDPS-2 Satellite

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

Analysis for Satellite Gravity Field Measurement Based on Relative Weight of Gravitational Coefficients

Gravity field modeling using space frequency signal transfer technique between satellites

Gravity field recovery of inter-satellite links between Beidou Navigation Satellite System (BDS) and LEO based on Geodesy and Time Reference in Space (GETRIS)

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

Accurate Establishment of Error Models for the Satellite Gravity Gradiometry Recovery and Requirements Analysis for the Future GOCE Follow-On Mission

The Analytical Methods for Gravity Field Measurement Performance of Satellite Formation

Recovered GRACE time-variable gravity field based on dynamic approach with the non-linear corrections

Improved recovery of temporal variations of the Earth’s gravity field from satellite kinematic orbits using an epoch-difference scheme

A hybrid approach for recovering high-resolution temporal gravity fields from satellite laser ranging

Performance Analysis and Error Budget for Gravity-Measuring Satellite

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

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

High-precision ground simulator for laser tracking of gravity satellite

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

Orbit design for a future geodetic satellite and gravity field recovery

Expected Precision of Gravity Gradient Recovered from Ka-Band Radar Interferometer Observations and Impact of Instrument Errors

Modeling and Analysis of Earth'S Gravity Field Measurement Performance by Inner-Formation Flying System

A New GNSS-R Altimetry Algorithm Based on Machine Learning Fusion Model and Feature Optimization to Improve the Precision of Sea Surface Height Retrieval