Solar and thermal radiation pressure modelling for improving the GOCE horizontal wind dataset

The Gravity Field and Steady-State Ocean Circulation Explorer (GOCE) satellite, which operated at an altitude of ∼250km, provided neutral thermosphere mass density and crosswind observations in the dawn-dusk sectors throughout most of its operational lifetime (2009–2013). As a result of its Sun-synchronous orbit, GOCE’s large solar panels remained at a near-perpendicular angle to the incoming solar radiation, leading to a significant radiation pressure acceleration. In this research, we focused on revisiting and reprocessing GOCE thermosphere mass density and crosswind data. We selected the coefficients describing the thermo-optical surface properties and employed a high-fidelity satellite geometry in a ray-racing simulation. Additionally, we distinguished between the solar flux in the visible and infrared bands and introduced a model for the satellite’s thermal emission. The availability of the in-situ thermistor measurements allowed for the validation of the thermal model. Moreover, we replaced the Level-1b ion thruster data with raw telemetry, filling multiple data gaps. We analysed how incremental improvements in the radiation pressure modelling affected the observed crosswind speed. By replacing the panel model with the high-fidelity satellite geometry, the crosswind speed decreased up to 5ms -1 . The biggest difference reduction of 40ms -1 resulted from introducing the thermal model. Splitting the solar flux further decreases the observed crosswind speed by up to 8ms -1 . The reduction in crosswind speed was most prominent during the first years of the mission when the solar activity was low. We compared the newly processed GOCE zonal wind data with respect to the most recent previous release. We observed a median absolute deviation decrease of 10ms -1 around the south magnetic pole in the dawn sector. The yearly consistency of low-latitude zonal winds did not change significantly. The main obstacle in quantifying the improvement compared to the previous crosswind dataset stemmed from the fact that the previous and new datasets were generated with different crosswind estimation algorithms. The difference in thermosphere density compared to previously published datasets is minor since the effect of radiation pressure is most prominent in the cross-track direction. Finally, we verified the assumption about the energy accommodation coefficient of 0.82 and concluded that it remains valid after implementing the radiation pressure modelling improvements.