This paper provides a comparative analysis of two different atmospheric drag models aimed at predicting the reentry of Resident Space Objects (RSO). We quantify the impact of considering a rotating atmosphere in the drag model, and analyze its effects on the Milankovitch elements, namely the angular momentum vector, and the Laplace-Runge-Lenz vector. The secular variation of the Milankovitch elements is expressed analytically through averaging. We evaluate the performance of the two formulations of the averaged equations of motion to provide accurate predictions of the orbital decay of RSO’s, by comparing simulated trajectories to those derived from the non-averaged dynamics, and to a documented spacecraft reentry event. Results show that accounting for atmospheric rotation in the averaged dynamics, provides simulated RSO trajectories closer to the non-averaged dynamics, and to the documented two-line element reentry data.