Uncertainties in radiation pressure modelling play a significant role in the thermospheric density and crosswind observations derived from the GRACE-FO accelerometer, especially during low solar activity. Under such conditions, the radiation pressure acceleration matches the magnitude of the aerodynamic acceleration along the track and exceeds it in the cross-track direction. The GRACE-FO mission has been operating for several years at such high altitudes during both low and rising solar activity, providing a perfect opportunity to study the effects of radiation pressure. This research uses ray tracing based on a high-fidelity satellite geometry model to calculate the radiation pressure acceleration. We numerically fine-tuned the coefficients describing the thermo-optical surface properties to obtain more accurate radiation pressure accelerations than those specified in the GRACE-FO mission manual. We also used in situ temperature measurements from thermistors on the solar arrays to model the satellite's thermal emission. These temperature measurements allowed a realistic setup of the thermal model, extended by the parameter describing the efficiency of the solar cells, and reproduced the acceleration of the thermal emission with an accuracy of RMS 0.148 nms−2. The combination of the updated thermal model and the fine-tuning of the surface coefficients improved the accuracy of the crosswind acceleration to an RMS of 0.55 nms−2, compared to an RMS of 4.22 nms−2 when using panel models and instantaneous thermal radiation. We compared the observed crosswind with two models: HWM14 and TIE-GCM. While both models capture most of the salient features of the observed crosswind, HWM14 shows particularly good agreement at high latitudes. Compared to the previously employed radiation pressure model, the crosswind observations have been improved in low and mid-latitudes, especially during periods of higher solar activity. Since the effect of radiation pressure is most significant in the crosswind direction, the effect on density was small compared to previously published datasets.@en

We present new neutral mass density and crosswind observations for the CHAMP, GRACE, and GRACE-FO missions, filling the last gaps in our database of accelerometer-derived thermosphere observations. For consistency, we processed the data over the entire lifetime of these missions, noting that the results for GRACE in 2011- 2017 and GRACE-FO are entirely new. All accelerometer data are newly calibrated. We modeled the temperature-induced bias variations for the GRACE accelerometer data to counter the detrimental effects of the accelerometer thermal control deactivation in April 2011. Further, we developed a new radiation pressure model, which uses ray tracing to account for shadowing and multiple reflections and calculates the satellitea's thermal emissions based on the illumination history. The advances in calibration and radiation pressure modeling are essential when the radiation pressure acceleration is significant compared to the aerodynamic one above 450 km altitude during low solar activity, where the GRACE and GRACE-FO satellites spent a considerable fraction of their mission lifetime. The mean of the new density observations changes only marginally, but their standard deviation shows a substantial reduction compared to thermosphere models, up to 15% for GRACE in 2009. The mean and standard deviation of the new GRACE-FO density observations are in good agreement with the GRACE observations. The GRACE and CHAMP crosswind observations agree well with the physics-based TIE-GCM winds, particularly the polar wind patterns. The mean observed crosswind is a few tens of m·s-1 larger than the model one, which we attribute primarily to the crosswind errors being positive by the definition of the retrieval algorithm. The correlation between observed and model crosswind is about 60%, except for GRACE in 2004- 2011 when the signal was too small to retrieve crosswinds reliably.@en