Effects of thermosphere total density perturbations on LEO orbits during severe geomagnetic conditions (Oct–Nov 2003) using DORIS and SLR data

An exceptional solar activity event occurred at the end of October 2003. On October 29th, seven groups of sunspots were visible on the Sun’s surface and the geomagnetic index Kp reached the extreme value of 9, leading to beautiful auroras, among other effects. Similar events also occurred in November 2003. These events have been an exceptional opportunity to estimate effects of the thermosphere total density perturbations in extreme conditions on a Low Earth Orbit (LEO orbit corresponding in this study at altitudes ranging from 800 km up to 1400 km). Specifically, we study the best way to get reliable geodetic products even during such solar events, and how well the thermosphere models (DTM-78, DTM-94, DTM-2000, and MSIS-86 models) allow us to accomplish this and to predict the observed perturbations on the orbital elements. Thus our analysis is twofold. First, we have computed orbits of satellites equipped with the on-board DORIS tracking system, at an altitude ranging from 800 to 1330 km and for periods of time including these exceptional events (October 29–31 and November 20, 2003). We have computed 30-hour arc orbits, estimating the drag coefficient parameter very frequently (every minute) using a tight random-walk constraint or using a white noise reset in the Gipsy/Oasis software. We show that significant errors are obtained for the considered thermosphere models, but can be greatly improved using a more sophisticated data processing. We also investigated how these proposed processing strategies affect the quality of the DORIS geodetic products. Significant improvements were found for weekly stations coordinates estimations as well as for polar motion determination. In addition, we have investigated the global quality of the modeling of the variations of the mean orbital elements in using thermosphere models over long periods of time including such severe geomagnetic events. This original approach is based on averaging techniques implemented in the CODIOR software. We have analyzed the continuous long-term evolution of the semi-major axis of the geodetic satellites Starlette, Stella and Ajisai tracked by the Satellite Laser Ranging (SLR) network. In this second part, we give: (i) a single global empirical coefficient per satellite, to validate the models over the whole duration of the orbital arc, including the exceptional solar activity events, and (ii) one coefficient per a period of about 2 months to quantify the seasonal differences between the models and the observations. As a result of both investigations, we compare the capability of recent models of thermosphere to allow us to get reliable geodetic products and mean orbital elements variations over short and long periods of time as well as for the recent exceptional geomagnetic events. Different altitudes are considered: around 800–900 km for the Stella and Starlette satellites tracked by laser ranging systems and SPOT tracked by DORIS; around also 1300–1400 km for the Ajisai satellite tracked by laser techniques and TOPEX and Jason also tracked by DORIS. It appears in our study that laser and DORIS data are complementary to probe the thermosphere, and to test the quality of thermosphere models in specific conditions for getting accurate geodetic products.

[1]  C. Lathuillère,et al.  From the Sun’s atmosphere to the Earth’s atmosphere: an overview of scientific models available for space weather developments , 2002 .

[2]  R. Walterscheid,et al.  Solar cycle effects on the upper atmosphere - Implications for satellite drag , 1989 .

[3]  D. Gambis,et al.  Monitoring Earth orientation using space-geodetic techniques: state-of-the-art and prospective , 2004 .

[4]  Y. Bar-Sever,et al.  Estimating horizontal gradients of tropospheric path delay with a single GPS receiver , 1998 .

[5]  M. Watkins,et al.  The gravity recovery and climate experiment: Mission overview and early results , 2004 .

[6]  V. Obridko,et al.  Cyclic variation of the global magnetic field indices , 1992 .

[7]  Y. Bar-Sever,et al.  One-Centimeter Orbit Determination for Jason-1: New GPS-Based Strategies , 2004 .

[8]  Pascal Willis,et al.  First assessment of GPS-based reduced dynamic orbit determination on TOPEX/Poseidon , 1994 .

[9]  F. Barlier,et al.  A thermospheric model based on satellite drag data , 2017 .

[10]  R. Pirjola,et al.  Space weather research elucidates risks to technological infrastructure , 2004 .

[11]  L. Jacchia Variations in the Earth's Upper Atmosphere as Revealed by Satellite Drag , 1963 .

[12]  Pierre Exertier,et al.  Semi-analytical theory of the mean orbital motion. , 1995 .

[13]  Ronald J. Muellerschoen,et al.  Topex/Jason combined GPS/DORIS orbit determination in the tandem phase , 2003 .

[14]  Pierre Exertier,et al.  Why the 18.6 year tide cannot explain the change of sign observed in j 2 , 2003 .

[15]  Gérard Thuillier,et al.  The DTM-2000 empirical thermosphere model with new data assimilation and constraints at lower boundary: accuracy and properties , 2003 .

[16]  W. A. Radasky,et al.  Advanced geomagnetic storm forecasting: a risk management tool for electric power system operations , 2000 .

[17]  F. Meyer Mathematical modelling of the sunspot cycle , 1981 .

[18]  Zuheir Altamimi,et al.  ITRF2000: A new release of the International Terrestrial Reference Frame for earth science applications , 2002 .

[19]  J. P. Berthias,et al.  Comportement de l'oscillateur DORIS/Jason au passage de l'anomalie sud-atlantique , 2004 .

[20]  An Analytical Solution of the Lagrange Equations Valid also for Very Low Eccentricities: Influence of a Central Potential , 2006 .

[21]  F. Barlier,et al.  Response of the equatorial thermosphere to magnetic activity analysed with accelerometer total density data. Asymmetrical structure , 1981 .

[22]  John C. Ries,et al.  Current status of the doris pilot experiment and the future international doris service , 2002 .

[23]  R. Biancale,et al.  An assessment of new satellite total density data for improving upper atmosphere models , 1999 .

[24]  P. Willis,et al.  DORIS as a potential part of a global geodetic observing system , 2005 .

[25]  Stephen H. Knowles,et al.  The Effect of Atmospheric Drag on Satellite Orbits During the Bastille Day Event , 2001 .

[26]  R. Biancale,et al.  Improvement of the empirical thermospheric model DTM: DTM94 – a comparative review of various temporal variations and prospects in space geodesy applications , 1998 .

[27]  A. Hedin MSIS‐86 Thermospheric Model , 1987 .

[28]  J. Lemoine,et al.  A high‐quality global gravity field model from CHAMP GPS tracking data and accelerometry (EIGEN‐1S) , 2002 .

[29]  Markus Rothacher,et al.  The International GPS Service (IGS): An interdisciplinary service in support of Earth sciences , 1999 .

[30]  C. Valorge,et al.  Precise Centre National d'Etudes Spatiales orbits for TOPEX/POSEIDON: Is reaching 2 cm still a challenge? , 1994 .

[31]  M. Menvielle,et al.  The K-derived planetary indices - Description and availability , 1991 .

[32]  Pascal Willis,et al.  External validation of the GRACE GGM01C gravity field using GPS and DORIS positioning results , 2004 .

[33]  S. Bruinsma,et al.  Semi-analytical theory of mean orbital motion: A new tool for computing ephemerides , 1997 .

[34]  Ch. Jayles,et al.  DORIS-DIODE: two-years results of the first European Navigator , 2002 .

[35]  F. LeMoine,et al.  The 1-Centimeter Orbit: Jason-1 Precision Orbit Determination Using GPS, SLR, DORIS, and Altimeter Data Special Issue: Jason-1 Calibration/Validation , 2003 .