Identification of the Gravitational Constant from Satellite Trajectories Using Neural Ordinary Differential Equations

Published in SPAICE 2026, 2026

Recommended citation: Katarina Dyreby, Claudia Soares, "Identification of the Gravitational Constant from Satellite Trajectories Using Neural Ordinary Differential Equations." SPAICE 2026, 2026. https://openreview.net/forum?id=Gf8ij9SHjP

Katarina Dyreby, Claudia Soares

Abstract

The rapid growth of satellite constellations in low Earth orbit has increased the probability of collisions involving satellites and debris, making accurate orbit propagation increasingly important. However, accurate orbit prediction remains challenging. Propagators depend on spacecraft properties, environmental parameters, and non-conservative forces that are often uncertain, causing propagation errors to accumulate over time. Estimating these values directly from orbital data offers a solution that could improve propagation accuracy. This work investigates whether a global physical constant can be identified directly from real orbital data using Neural Ordinary Differential Equations (Neural ODEs). A hybrid Neural ODE model is developed in which the equations of motion of satellites are split into two parts: known gravitational dynamics, with the gravitational constant treated as a trainable parameter, and a Neural Network that learns the residual accelerations representing the unmodelled forces. Using different initial biases, the model successfully recovers the gravitational constant for relative errors of up to 150%. Increasing the initialization error primarily increases the number of training epochs required for convergence, as the final relative errors achieved are below 0.02% for all different initial biases used. At higher initial errors, the Neural Network instead compensates for the incorrect gravitational term by learning residual accelerations several orders of magnitude larger than physically expected, preventing recovery of the gravitational constant. These results show that Neural ODEs can recover physical parameters from real orbital data, while also revealing the limits of parameter identifiability when unknown dynamics are learned simultaneously.

Access paper here

BibTeX

@inproceedings{spaice2026Gf8ij9SHjP,
    author = "Dyreby, Katarina and Soares, Claudia",
    title = "Identification of the Gravitational Constant from Satellite Trajectories Using Neural Ordinary Differential Equations",
    booktitle = "SPAICE 2026",
    year = "2026",
    month = "oct",
    day = "21",
    url = "https://openreview.net/forum?id=Gf8ij9SHjP",
    abstract = "The rapid growth of satellite constellations in low Earth orbit has increased the probability of collisions involving satellites and debris, making accurate orbit propagation increasingly important. However, accurate orbit prediction remains challenging. Propagators depend on spacecraft properties, environmental parameters, and non-conservative forces that are often uncertain, causing propagation errors to accumulate over time. Estimating these values directly from orbital data offers a solution that could improve propagation accuracy. This work investigates whether a global physical constant can be identified directly from real orbital data using Neural Ordinary Differential Equations (Neural ODEs). A hybrid Neural ODE model is developed in which the equations of motion of satellites are split into two parts: known gravitational dynamics, with the gravitational constant treated as a trainable parameter, and a Neural Network that learns the residual accelerations representing the unmodelled forces. Using different initial biases, the model successfully recovers the gravitational constant for relative errors of up to 150\\%. Increasing the initialization error primarily increases the number of training epochs required for convergence, as the final relative errors achieved are below 0.02\\% for all different initial biases used. At higher initial errors, the Neural Network instead compensates for the incorrect gravitational term by learning residual accelerations several orders of magnitude larger than physically expected, preventing recovery of the gravitational constant. These results show that Neural ODEs can recover physical parameters from real orbital data, while also revealing the limits of parameter identifiability when unknown dynamics are learned simultaneously.",
    note = "Accepted for oral presentation"
}