Atmospheric density is a critical component for studies of meteor flight, artificial space debris reentry and space object (satellite and debris) lifetime.
Atmospheric density varies with altitude, geographical latitude and longitude, time, season and atmospheric dynamics. In the upper atmosphere (>100 km) it also depends strongly on the solar soft x-ray and extreme ultraviolet (EUV) output. It also varies according to the geomagnetic activity. In the lower atmosphere a mean density as a function of altitude is normally sufficient for models of meteors and space debris reentry. In the upper atmosphere for lifetime studies which consider atmospheric drag the density must be specified as a function of altitude, solar ten centimetre flux and the geomagnetic activity index Ap.
The following sections will consider a number of simple atmospheric density models that are useful for the above applications.
AN ISOTHERMAL ATMOSPHERE
An isothermal atmosphere is one that has a constant temperature. In such an atmosphere the density (ρ) can be written as
Below 100km the variation of gravity with altitude is insignificant. Thus the only factor that varies significantly is the temperature.
THE REAL ATMOSPHERE
|The graph at the left shows the temperature variation with height that is used in the US Standard Atmosphere model. It is a mean of mostly mid-latitude conditions in the northern hemisphere. Some is derived from extensive data (eg below 20 km). Some is derived from quite limited data (eg between 50 and 80 km).
From the surface the temperature intially decreases at what is called the standard adiabatic lapse rate of about 7 degrees per kilometre. It then plateaus at -55 Celcius for almost 10 km, before rising again to another plateau just below 0 C. It decreases again in the Mesosphere before finally rising when it reaches the Thermosphere.
The changing temperature causes the scale height to change. From the surface to 100 km altitude H varies between about 5.5 and 8.0 km , rising and falling along with the temperature.
Despite the changing scale height, a simple isothermal model with a fixed scale height can represent the atmosphere up to 100 km with an error less than 50%, and up to 130 km with an error no larger than a factor of 2. These errors are quite comparable to the variation that the real atmosphere shows from day to day from the mean 'standard' model. This model is quite suitable for demonstrating the physics of meteor flight and space debris reentries. High velocity meteoroids (~60 km/s) start ablating around 130 km altitude, but space debris reentries, with the much lower velocity of ~ 8 km/s only start to show significant atmospheric interaction below 90 km.
The simple isothermal model is specified by the formula:
THE CIRA MODEL
The CIRA model is the COSPAR International Reference Atmosphere. COSPAR is the United Nations Committee for Space Research. Graphs of CIRA atmospheric temperature and density to 500 km altitude are shown below:
The changing temperature is the predominant parameter that changes the scale height and hence the density. However, this is reinforced by the decrease in gravity (g) with increasing altitude and also the reduction in mean molecular mass above 100 km altitude. Below 100km the Sun has only a minor influence on the density. Above 200 km it becomes quite significant. This is because the solar EUV and soft X-rays which affect the upper atmospheric density are absorbed around the 100 km altitude level.
MODELS TO 180 KM
We can use the mean CIRA model to 180 km by specifying the density values every 5 or 10 km. Because the density varies over many orders of magnitude over this range it is more useful to use the logarithm (base 10) of these values. Linear interpolation can be used to calculate the values intermediate between the specified points at the 10 km intervals. The actual density can then be calculated by raising the interpolated log values as a
power of ten. That is, ρ = 10logval.
An alternative is to fit a polynomial to the log points such that:
= ((((((a7*h + a6)*h + a5)*h + a4)*h + a3)*h + a2)*h + a1)*h + a0
[The lower form is identical but computationally more efficient]
A MODEL FROM 180 to 500 KM
Above 180 km we can no longer use a mean reference model. We need a model that varies according to the space weather conditions at the time, specifically the solar EUV flux and the geomagnetic activity. We use the solar radio ten centimetre flux (F10) as a proxy for the solar EUV output, and the geomagnetic Ap index as a proxy for the geomagnetic activity. The model we have developed is suitable for investigating orbital decay and satellite and debris lifetimes in the altitude range from 180 to 500 km.
The algorithm to calculate the atmospheric density in this range is given by:
|T = 900 + 2.5 (F10 – 70) + 1.5 Ap||(Kelvin)|
|μ = 27 – 0.012 (h – 200)||180 < h (km) < 500|
|H = T / μ||(km)|
|ρ = 6x10-10 exp( - (h –175) / H )||( kg m-3)|
Values of the model atmospheric density are graphed below for the cases of minimal (F10=70) and maximal (F10=300) solar EUV and soft X-ray output:
Note: it has been found that an average F10 value is the more appropriate parameter to use rather than the daily value. Typically the last 90 days are used to form the average.
Above 500 km altitude the lifetime of an orbital body is larger than the solar period of ~11 years, and so once again an average or mean atmospheric model may be appropriate.
Australian Space Academy