Low Earth Orbit (LEO) satellites are significantly affected by variable air drag forces, which are mainly altered by atmospheric expansion/contraction driven by solar and geomagnetic activity. Air drag reduces the orbital velocity of a satellite, its nominal altitude, and shortens its lifespan. The effect of air drag pressure on the position of a satellite orbiting at an altitude of around 450 km may drag around 3 m per revolution in the along-track axis, limiting the satellite's lifespan to approximately 5-10 years. In applications, such as remote sensing or satellite altimetry and gravity, the orbital trajectory and velocity (ephemeris) of satellites must be known to an accuracy of a few millimeters. Moreover, the exponential increase in presence of space debris (consider the recent destructive events of Fengyun-1C, Iridium, and Mission Shakti) has highlighted the importance of orbital tracking and prediction of potential collisions. The dynamic Precise Orbit Determination (POD) method tracks and predicts the orbital ephemeris by calculating an orbital trajectory through a double integration and linearization of Newton-Euler's equation of motion (Montenbruck & Gill, 2013). In the POD method, by combining force models with empirical observations, used for example in laser-ranging, Doppler, accelerometer, or Global Navigation Satellite System (GNSS) measurements, the position, and velocity of a satellite can be stochastically estimated with significant accuracy (Jin & Su, 2020;Tapley et al., 2004).Due to variable air drag force being so important, in the last decade, thermospheric mass density (TMD) variations driven by solar and geomagnetic activity have been investigated using satellite technology to a great extent (e.g.