Analysis of the orbital decay of spherical satellites using different solar flux proxies and atmospheric density models

Abstract In order to highlight the performances of the JR-71 and MSISE-90 atmospheric density models and also to investigate whether the use of the E 10.7 index, instead of the F 10.7 cm solar radio flux, improves the trajectory modeling of low Earth objects, the orbital decay of 11 spherical satellites was analyzed. The data were distributed over a full solar activity cycle and covered altitudes in between 150 and 1500 km. The observed evolution of the mean semi-major axis was obtained by processing the historical two-line elements record. For each satellite, environmental condition, atmospheric model and solar flux proxy, a drag coefficient was computed by fitting the observed semi-major axis decay, comparing the results with the theoretical estimates. Below 400 km, both models overestimated the atmosphere density, independently on the solar flux proxy and environmental conditions. MSISE-90, with F 10.7 , proved to be the best model to compute the air density, in low solar activity conditions. JR-71, with F 10.7 , was generally more accurate over the maximum of the 23rd solar cycle, while MSISE-90 revealed to be the best model in the increasing phase of the solar flux, when using E 10.7 . Between 750 and 850 km, the average accuracy of the models was generally better than 15% and the choice of the atmospheric density model was relatively inconsequential, both in conditions of low and high solar activity. A quite large uncertainty was instead observed in moderate conditions, when using E 10.7 in both models. At 1500 km, both models underestimated the air density by more than 35%, in low solar activity conditions. For high solar fluxes, the density was well represented by JR-71, both using F 10.7 and E 10.7 , while it resulted to be underestimated, by about 25%, when using MSISE-90, both with E 10.7 and F 10.7 . In general, more accurate orbital decay fits were obtained with E 10.7 , which seems able to better describe the atmosphere density changes induced by a varying solar activity. In conclusion, none of the models analyzed was able to accurately describe the atmosphere density in all orbital regimes and environmental conditions, independently on the solar flux proxy.

[1]  W. Tobiska,et al.  New developments in SOLAR2000 for space research and operations , 2006 .

[2]  Gary J. Rottman,et al.  The SOLAR2000 empirical solar irradiance model and forecast tool , 2000 .

[3]  Edward M Gaposchkin Calculation of Satellite Drag Coefficients , 1994 .

[4]  Carmen Pardini,et al.  Comparison and accuracy assessment of semi-empirical atmosphere models through the orbital decay of spherical satellites , 2001 .

[5]  Graham G. Swinerd,et al.  A free molecule aerodynamic investigation using multiple satellite analysis , 1996 .

[6]  J. Cappellari,et al.  Mathematical theory of the Goddard trajectory determination system , 1976 .

[7]  W. Tobiska Validating the solar EUV Proxy, E10.7 , 2001 .

[8]  Requirements for accurate near-real time atmospheric density correction , 2000 .

[9]  D. King-hele,et al.  Satellite orbits in an atmosphere : theory and applications , 1987 .

[10]  G. E. Cook Drag coefficients of spherical satellites , 1965 .

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

[12]  A. Hedin Extension of the MSIS Thermosphere Model into the middle and lower atmosphere , 1991 .

[13]  F. Mignard,et al.  Reassessment of the charge and neutral drag of LAGEOS and its geophysical implications , 1985 .

[14]  C. Chao,et al.  An Evaluation of Jacchia 71 and MSIS90 Atmosphere Models with NASA ODERACS Decay Data , 1997 .

[15]  J. Lean,et al.  Global change in the thermosphere: Compelling evidence of a secular decrease in density , 2004 .