Validation of GOCE densities and evaluation of thermosphere models

Atmospheric densities from ESA’s GOCE satellite at a mean altitude of 270 km are validated by comparison with predictions from the near real time model HASDM along the GOCE orbit in the time frame 1 November 2009 through 31 May 2012. Except for a scale factor of 1.29, which is due to different aerodynamic models being used in HASDM and GOCE, the agreement is at the 3% (standard deviation) level when comparing daily averages. The models NRLMSISE-00, JB2008 and DTM2012 are compared with the GOCE data. They match at the 10% level, but significant latitude-dependent errors as well as errors with semiannual periodicity are detected. Using the 0.1 Hz sampled data leads to much larger differences locally, and this dataset can be used presently to analyze variations down to scales as small as 150 km.     

Use of two-line element data for thermosphere neutral density model calibration

Traditional empirical thermospheric density models are widely used in orbit determination and prediction of low-Earth satellites. Unfortunately, these models often exhibit large density errors of up to around 30% RMS. Density errors translate into orbit errors, adversely affecting applications such as re-entry operations, manoeuvre planning, collision avoidance and precise orbit determination for geodetic missions. The extensive database of two-line element (TLE) orbit data contains a wealth of information on satellite drag, at a sufficiently high spatial and temporal resolution to allow a calibration of existing neutral density models with a latency of one to two days. In our calibration software, new TLE data for selected objects is converted to satellite drag data on a daily basis. The resulting drag data is then used in a daily adjustment of density model calibration parameters, which modify the output of an existing empirical density model with the aim of increasing its accuracy. Two different calibration schemes have been tested using TLE data for about 50 objects during the year 2000. The schemes involve either height-dependent scale factors to the density or corrections to CIRA-72 model temperatures, which affect the density output based on a physical model. Both schemes have been applied with different spherical harmonic expansions of the parameters in latitude and local solar time. Five TLE objects, varying in perigee altitude between 280 and 530 km, were deliberately not used during calibration, in order to provide independent validation. Even with a single daily parameter, the RMS density model error along their tracks can already be reduced from the 30% to the 15% level. Adding additional parameters results in RMS errors lower than 12%.     

Horizontal and vertical thermospheric cross-wind from GOCE linear and angular accelerations

Thermospheric wind measurements obtained from linear non-gravitational accelerations of the Gravity field and steady-state Ocean Circulation Explorer (GOCE) satellite show discrepancies when compared to ground-based measurements. In this paper the cross-wind is derived from both the linear and the angular accelerations using a newly developed iterative algorithm. The two resulting data sets are compared to test the validity of wind derived from angular accelerations and quantify the uncertainty in accelerometer-derived wind data. In general the difference is found to be less than 50?m/s vertically after high-pass filtering, and 100?m/s horizontally. A sensitivity analysis reveals that continuous thrusting is a major source of uncertainty in the torque-derived wind, as are the magnetic properties of the satellite. The energy accommodation coefficient is identified as a particularly promising parameter for improving the consistency of thermospheric cross-wind data sets in the future. The algorithm may be applied to obtain density and cross-wind from other satellite missions that lack accelerometer data, provided the attitude and orbit are known with sufficient accuracy.     

High-fidelity geometry models for improving the consistency of CHAMP,GRACE, GOCE and Swarm thermospheric density data sets

During the last two decades, accelerometers on board of the CHAMP, GRACE, GOCE and Swarm satellites have provided high-resolution thermosphere density data to improve our knowledge on atmospheric dynamics and coupling processes in the thermosphere-ionosphere region. Most users of the data have focused on relative density variations. Scale differences between datasets and models have been largely neglected or removed using ad hoc scale factors. The origin of these scale differences arises from errors in the aerodynamic modelling, specifically in the modelling of the satellite outer surface geometry and of the gas-surface interactions. Therefore, the first step to remove the scale differences is to enhance the geometry modelling. This work forms the foundation for the future improvement of characterization of satellite aerodynamics and gas-surface interactions models at TU Delft, as well as for extending the use of sideways and angular accelerations in the aerodynamic analysis of accelerations and derivation of thermosphere datasets. Although work to improve geometry and aerodynamic force models by other authors has focused on CHAMP and GRACE, this paper includes the GOCE and Swarm satellites as well. In addition, it uses a density determination algorithm that is valid for arbitrary attitude orientations, enabling a validation making use of attitude manoeuvres. The results show an improvement in the consistency of density data between these four missions, and of data obtained before, during and after attitude manoeuvres of CHAMP and Swarm. The new models result in larger densities, compared to the previously used panel method. The largest average rescaling of density, by switching to the new geometry models is reached for Swarm at 32%, the smallest for GRACE at 5%. For CHAMP and GOCE, mean differences of 11% and 9% are obtained respectively. In this paper, an overview of the improvements and comparisons of data sets is provided together with an introduction to the next research phase on the gas-surface interactions.     

CHAMP and GRACE accelerometer calibration by GPS-based orbit determination

Current and planned Earth observation missions are equipped with highly sensitive accelerometers. Before using the data, the instrument has to be calibrated by determining scale and bias parameters for each axis. Here, the accelerometer measurements are used in a GPS-based reduced-dynamic orbit determination approach, replacing the non-gravitational force models, and nominally daily calibration parameters are estimated. Additional empirical accelerations are estimated to account for deficiencies in the applied force models. This method is applied to 5 years of CHAMP and GRACE data, resulting in an orbit precision at the level of a few centimeters. In along-track direction the calibration parameters can be estimated freely, scale factors of 0.96 ± 0.014 and 0.95 ± 0.015 are obtained for GRACE A and B, and 0.85 ± 0.024 for CHAMP. A constant scale factor results in the smoothest bias series, with clear trends and occasional jumps. In radial and cross-track direction tight constraints to a priori biases have to be applied. Furthermore, the determined orbits are analyzed with respect to reference trajectories, and SLR, phase and KBR residuals are presented.