Low Earth orbit (LEO) satellites are leveraged to support new position, navigation, and timing (PNT) service alternatives to GNSS. These alternatives require accurate propagation of satellite position and velocity with a realistic quantification of uncertainty. It is commonly assumed that the propagated uncertainty distribution is Gaussian; however, the validity of this assumption can be quickly compromised by the mismodeling of atmospheric drag. We develop a machine learning approach that corrects error growth in the argument of latitude for a diverse set of LEO satellites. The improved orbit propagation accuracy extends the applicability of the Gaussian assumption and modeling of the errors with a corrected mean and covariance. We compare the performance of a time-conditioned neural network and a Gaussian Process on datasets computed with an open source orbit propagator and publicly available Vector Covariance Message (VCM) ephemerides. The learned models predict the argument of latitude error as a Gaussian distribution given parameters from a single VCM epoch and reverse propagation errors. We show that this one-dimensional model captures the effect of mismodeled drag, which can be mapped to the Cartesian state space. The correction method only updates information along the dimensions of dominant error growth, while maintaining the physics-based propagation of VCM covariance in the remaining dimensions. We therefore extend the utility of VCM ephemerides to longer time horizons without modifying the functionality of the existing propagator.
Figures & tables
Fig. 1: Machine learning model predicts errors and variance in the argument of latitude to correct both the propagated orbit and covariance for a variety of LEO orbit and satellite types.
Fig. 2: Orbit propagation correction pipeline using VCMs, open source SP propagator, and machine learning model.
Parameter
Description
Epoch time (UTC)
Orbital estimate epoch time, expressed in Coordinated Universal Time.
J2K Pos/Vel X-Y-Z ( km , km/s )
The position and velocity vector components of the satellite expressed in Earth-Centered Inertial reference frame.
Ballistic Coef. ( m2/kg )
Estimated ballistic coefficient at epoch.
Solar Rad Press Coeff ( m2/kg )
Estimated solar radiation pressure coefficient at epoch.
Geo Pot
Degree and order of geopotential model used.
Drag
Atmospheric density model used.
TABLE I: Summary of VCM parameters used in modeling and dataset generation
Fig. 3: Accelerations due to individual force models in Orekit orbit propagator over various altitudes for an example satellite in an initially circular orbit with a ballistic coefficient of 0.13 m2/kg .
Perturbation Force Models
Solar Radiation Pressure
70 × 70 Gravity (EGM-96)
Solid Earth Tides
Sun Third Body
Moon Third Body
Atmospheric Drag (NRLMSISE-00)
TABLE II: Orekit Force Models
Fig. 4: Example of error growth for propagation of a single VCM (NORAD ID 47). True AOL error maps to predominantly along-track position and radial velocity dimensions. Errors in the remaining dimensions are caused by errors in other orbital elements and therefore are not captured by argument of latitude.
Fig. 5: Over 95% of the norm position error is captured by the argument of latitude after 5 day propagation (NORAD ID 47).
Fig. 6: AOL errors used to train the ML models consisting of 50 VCMs propagated 7 days for 800 satellites.
σRr km
σRs km
σRw km
σVr m/s
σVs m/s
σvw m/s
0.23
12
0.013
13
0.043
0.017
TABLE III: Validation dataset standard deviations in RSW frame after seven days of propagation
Fig. 7: TCNN - Corrected and Uncorrected Error Distributions Over Time.
Fig. 8: TCNN - Error reduction in all dimensions of ECI Cartesian frame after AOL correction.
Dataset
σRs km
σVr m/s
σ∥R∥ km
σ∥V∥ m/s
% DMu,C
% DM[R,V],C
Uncorrected
12
13
11
12
-
7.7
TCNN
6.2
6.8
5.5
6.1
95
81
HGP
6.8
7.5
6.1
6.7
97
81
TABLE IV: Correction performance for HGP and TCNN for 7 day propagation
Fig. 9: Covariance realism improvement after AOL correction in ECI Cartesian frame.
Dataset
σRs km
σVr m/s
σ∥R∥ km
σ∥V∥ m/s
% DMu,C
Uncorrected
12
13
11
12
-
TCNN
11
12
11
12
64
HGP
9.9
11
9.3
10
96
TABLE V: Correction performance for HGP and TCNN for 7 day propagation with time generalization
Uncertainty propagation studies how measurement, model, and data uncertainty flows through computations and affects downstream predictions and decisions.