Bring The Noise: Three Rules for Improving Thermospheric Density Retrieval From LEO POD Data

Thermospheric mass density knowledge remains a primary driver of operational uncertainty for Low Earth Orbit operations, yet current sources of density observations suffer from limited latency, accuracy, and/or licensing issues. The growing body of Precise Orbit Determination data from emerging Low Earth Orbit constellations offers a promising source of thermospheric mass density measurements. Realising this potential requires the application of methods which can convert Precise Orbit Determination data into reliable density observations. This work provides an assessment of the Energy Dissipation Rate technique. Using a simulated constellation of up to eighteen Starlink-like satellites, we quantify the relationship between the quality of retrieved mass densities and four major influencing factors: Precise Orbit Determination noise level, drag acceleration magnitude, estimation fit-span, and the number of co-orbital spacecraft whose independent density estimates are averaged. Densities retrieved via the Energy Dissipation Rate method are also compared with densities that have been assimilated into a background density model. We identify conditions under which non-assimilated retrievals can outperform their assimilated counterparts, demonstrating the benefits of averaging data from multiple satellites along the same orbital track. Our results reveal three rules of thumb that describe how retrieval error scales as a function of fit-span, satellite count, and drag level. Doubling the fit-span roughly halves the RMS percentage error; adding co-orbital satellites reduces error as ~ N-0.6 ∼ N - 0.6 (unassimilated) or ~ N-0.25 ∼ N - 0.25 (assimilated); and drag-acceleration sensitivity follows a logistic response with the steepest gains between 2× 10⁻⁸ 2 × 10 - 8 and 2× 10⁻⁷\, m\,s⁻² 2 × 10 - 7 m s - 2 . Together, these findings provide the first systematic quantification of how density-retrieval errors scale with Precise Orbit Determination noise, drag-acceleration level, fit-span, and satellite count. The resulting analytical approximation captures over 89% of the variance in simulated retrieval errors and is intended to support the design of next-generation assimilative thermospheric mass density models.