Improving the estimation of thermospheric neutral density via two-step assimilation of in situ neutral density into a numerical model

Earth, Planets and Space

Full Paper

Armin Corbin¹  Jürgen Kusche¹ 

[1] Institute of Geodesy and Geoinformation, University of Bonn, Nussallee 17, 53115, Bonn, Germany;

关键词: Thermosphere;    Neutral mass density;    Data assimilation;    Geomagnetic storm;   

DOI  :  10.1186/s40623-022-01733-z

 received in 2022-05-05, accepted in 2022-11-03,  发布年份 2022

来源: Springer

【 摘 要 】

Neutral thermospheric density is an essential quantity required for precise orbit determination of satellites, collision avoidance of satellites, re-entry prediction of satellites or space debris, and satellite lifetime assessments. Empirical models of the thermosphere fail to provide sufficient estimates of neutral thermospheric density along the orbits of satellites by reason of approximations, assumptions and a limited temporal resolution. At high solar activity these estimates can be off by 70% when comparing to observations at 12-hourly averages. In recent decades, neutral density is regularly observed with satellite accelerometers on board of low Earth orbiting satellites like CHAMP, GOCE, GRACE, GRACE-FO, or Swarm. When assimilating such along-track information into global models of thermosphere–ionosphere dynamics, it has been often observed that only a very local sub-domain of the model grid around the satellite’s position is updated. To extend the impact to the entire model domain we suggest a new two-step approach: we use accelerometer-derived neutral densities from the CHAMP mission in a first step to calibrate an empirical thermosphere density model (NRLMSIS 2.0). In a second step, we assimilate—for the first time—densities predicted for a regular three-dimensional grid into the TIE-GCM (Thermosphere Ionosphere Electrodynamics General Circulation Model). Data assimilation is performed using the Local Error-Subspace Transform Kalman Filter provided by the Parallel Data Assimilation Framework (PDAF). We test the new approach using a 2-week-long period containing the 5 April 2010 Geomagnetic storm. Accelerometer-derived neutral densities from the GRACE mission are used for additional evaluation. We demonstrate that the two-step approach globally improves the simulation of thermospheric density. We could significantly improve the density prediction for CHAMP and GRACE. In fact, the offset between the accelerometer-derived densities and the model prediction is reduced by 45% for CHAMP and 20% for GRACE when applying the two-step approach. The implication is that our approach allows one to much better ’transplant’ the precise CHAMP thermospheric density measurements to satellites flying at a similar altitude.Graphical Abstract

【 授权许可】

CC BY   
© The Author(s) 2022

【 图 表 】

Fig. 4

Fig. 3

Fig. 3

Fig. 5

Fig. 1

12888_2022_4451_Article_IEq2.gif

Fig. 5

12982_2022_119_Article_IEq19.gif

12864_2022_9026_Article_IEq187.gif

Fig. 1

Fig. 1

Fig.1

Fig. 6

12864_2022_9026_Article_IEq59.gif

Fig. 1

12864_2022_9026_Article_IEq57.gif

Fig. 1

12864_2022_9026_Article_IEq55.gif

Fig. 3

Fig. 1

41408_2022_764_Article_IEq12.gif

41408_2022_764_Article_IEq11.gif

41408_2022_764_Article_IEq10.gif

41408_2022_764_Article_IEq9.gif

41408_2022_764_Article_IEq8.gif

41408_2022_764_Article_IEq7.gif

41408_2022_764_Article_IEq6.gif

Fig. 2

41408_2022_764_Article_IEq3.gif

Fig. 1

41408_2022_764_Article_IEq1.gif

Fig. 2

Fig. 5

Fig. 1

Fig. 1

40249_2022_1045_Article_IEq27.gif

Fig. 5

40249_2022_1045_Article_IEq25.gif

40249_2022_1045_Article_IEq24.gif

40249_2022_1045_Article_IEq23.gif

40249_2022_1045_Article_IEq22.gif

40249_2022_1045_Article_IEq21.gif

40249_2022_1045_Article_IEq20.gif

40249_2022_1045_Article_IEq17.gif

40249_2022_1045_Article_IEq16.gif

Fig. 1

40249_2022_1045_Article_IEq15.gif

40249_2022_1045_Article_IEq14.gif

40249_2022_1045_Article_IEq13.gif

40249_2022_1045_Article_IEq12.gif

40249_2022_1045_Article_IEq11.gif

40249_2022_1045_Article_IEq10.gif

40249_2022_1045_Article_IEq9.gif

40249_2022_1045_Article_IEq8.gif

40249_2022_1045_Article_IEq7.gif

40249_2022_1045_Article_IEq6.gif

Fig. 1

40249_2022_1045_Article_IEq1.gif

Fig. 4

Fig. 1

Fig. 1

Fig. 3

Fig. 1

Fig. 2

Fig. 3

Fig. 2

13731_2022_257_Article_IEq4.gif

Fig. 2

Fig. 1

Fig. 2

40517_2022_243_Article_IEq12.gif

Fig. 1