An accurate estimation of upper atmospheric densities is crucial for precise orbit determination (POD), prediction low Earth satellites, and scientific studies the Earth's atmosphere. But estimated using satellite tracking data are always uncertain up to drag-coefficient assumed in inversion method. This work develops a new framework simultaneously estimate density satellites with time-varying ...