The variation of air density at 240 and 280 km from April to November 1967

Changes in the orbital periods of two satellites, 1962-βτ6 (Injun 3 rocket) and 1965-11D (Cosmos 54 rocket), have been used to deduce the air density at heights of 240 and 280 km during April–November 1967. At both heights the generally low density observed in July and the higher density in April and October were almost certainly part of a semi-annual variation similar to that observed at other heights in the thermosphere. The ratio of the maximum (October) to minimum (July) density was about 1·8 at 240 km and 2·2 at 280 km. Superimposed upon this variation were short-lived increases in density associated with magnetic storms, the largest being of 65 per cent at 280 km on 25 May, and a periodic variation with an amplitude of up to 25 per cent from the monthly mean density, related to the 10·7 cm solar radiation flux. A diurnal variation of density was also detected with a maximum density at 14 hr and a maximum to minimum ratio of 1·7 at 280 km.     

Upper-atmosphere features revealed by the orbit of 1980-43A

The satellite NOAA-B (1980-43A) was launched in May 1980 into an orbit with perigee height near 260 km and apogee height 1440 km, at an inclination of 92.2°.The lifetime was 11 months. The orbit has been determined at 40 epochs between October 1980 and May 1981 from about 3000 radar and optical observations. The average orbital accuracy, radial and cross-track, was about 100 m, with rather better accuracy in the final 14 days.The variation of orbital inclination has been analysed to determine two good values of atmospheric rotation rate, namely 1.10 ± 0.10 rev day^?1 at 300 km (average local time) and 1.15 ± 0.06 rev day^?1 at 225 km (evening).The effect of atmospheric rotation on the precession of the orbital plane of an actual satellite has never previously been detected; it is clearly apparent for 1980-43A in its last days and conforms to the expected theoretical change.The variation of perigee height has been analysed to determine ten values of atmospheric density scale height, for heights of 280–370 km. These values, accurate to about 3%, exceed by 15% the values indicated by the COSPAR International Reference Atmosphere. Solar activity was higher in the years 1980–1981 than at any time since early 1958 and it appears that the CIRA model underestimates the density and density scale height at high levels of solar activity.     

Determination and geophysical interpretation of the orbit of China 2 rocket (1971-18B)

The orbit of China 2 rocket, 1971-18B, has been determined at 114 epochs throughout its 5-yr life, using the RAE orbit refinement program PROP 6, with more than 7000 radar and optical observations from 83 stations.The rocket passed slowly enough through the resonances 14:1, 29:2, 15:1 and 31:2 to allow lumped geopotential harmonic coefficients to be calculated for each resonance, by least-squares fittings of theoretical curves to the perturbation-free values of inclination and eccentricity. These lumped coefficients can be combined with values from satellites at other inclinations, to obtain individual harmonic coefficients.The rotation rate of the upper atmosphere, at heights near 300 km, was estimated from the decrease in orbital inclination, and values of 1.15, 1.05, 1.10 and 1.05 rev/day were obtained between April 1971 and January 1976. From the variation in perigee height, 25 values of density scale height were calculated, from April 1971 to decay. Comparison with values from the COSPAR International Reference Atmosphere 1972 shows good agreement between April 1971 and October 1975, but the observational values are 10% lower, on average, than CIRA thereafter.A further 1400 observations, made during the final 15 days before decay, were used to determine 15 daily orbits. Analysis of these orbits reveals a very strong West-to-East wind, of 240 ± 40 ms^?1, at a mean height of 195 km under winter evening conditions, and gives daily values of density scale height in the last 7 days before decay.     

Features of the upper atmosphere revealed by analysis of the orbit of cosmos 1009 rocket, 1978-50B

COSMOS 1009 rocket was launched on 19 May 1978 into an orbit with initial perigee height 150 km and apogee 1100 km: its lifetime was only 17 days. The orbit has been determined daily during the final 14 days of its life, using the RAE orbit refinement program PROP6,with about 1100 observations supplied by NORAD. An average accuracy of about 60 m, radial and cross-track, was achieved.The orbits were analysed to reveal three features of the upper atmosphere at heights between 125 and 175 km. From the decrease in perigee height, five values of density scale height, accurate to ±4%, were obtained. The first three were within 10% of those from CIRA 1972; the fourth, after a magnetic storm, was higher than expected; the fifth gave evidence of the decrease in drag coefficient at heights below 130 km.Atmospheric oblateness produced a change of 4° in perigee position during the last four days of the life. Analysis showed that the ellipticity of the upper atmosphere was approximately equal to that of the Earth, f, for the first two of the four days, and about $\text{1}\text{2}\text{f}$ in the last two.The orbital inclination decreased during the 14 days by about 50 times its standard deviation, and the observed variation was analysed to determine zonal winds at heights of 150–160 km at latitudes near 47° north. The zonal wind was very weak (0±30 m/s) for 23–28 May at local times near 03h; and 90±30 m/s east-to-west for 29 May to 4 June at local times near 01 h.     

Analysis of the orbit of 1965-11D (COSMOS 54 rocket)

The satellite 1965-11D was the final-stage rocket used to launch Cosmos 54, 55 and 56 into orbit on 21 February 1965. The orbit of 1965-11D was inclined at 56° to the Equator, with an initial perigee height of 280 km; the lifetime was nearly 5 yr, with decay on 23 December 1969. The orbit has been determined at 75 epochs during the life, using the RAE orbit determination program PROP with over 4000 observations, photographic, visual and radar. Observations from the Hewitt camera at Malvern were available for 34 of the 75 orbits and typical accuracies for these orbits are 0.0005° in inclination and 100 m in perigee height.The variations in perigee height have been analyzed to determine reliable values of density scale height, at heights between 240 and 360 km. The analysis also revealed a rapid decrease of 5 km in perigee distance early in 1966, attributed to the escape of residual propellants.The variations in orbital inclination have been analyzed to determine upper-atmosphere zonal winds and 15th-order harmonics in the geopotential. The region of the upper atmosphere traversed by 1965-11D near its perigee is found to have had an average rotation rate of 1.10 ± 0.05 rev/day in 1966–1967, and 1.00 ± 0.03 rev/day between March 1968 and May 1969. In late 1969 there were probably wide variations in zonal winds, with east-to-west winds of order 100 m/s followed by west-to-east winds of order 200 m/s. The changes in inclination at the 15th-order resonance in July 1969 have been analyzed to give the first accurate values of lumped 15th-order harmonics obtained from a high-drag satellite. This success points the way towards similar analyses of the many other high-drag satellites that pass through 15th-order resonance, to evaluate individual geopotential coefficients of order 15 and even degree.     

Upper-atmosphere rotation rate,from analysis of the orbit of 1970-43B

The orbit of Cosmos 347 rocket (1970-43B) has been determined in the form of 23 sets of orbital elements at intervals during its 8-month life, with the aid of the RAE orbit improvement program PROP, using about 850 observations from 47 observing stations. The values of orbital inclination obtained, which had standard deviations between 0.7 and 10 sec of arc, were analysed to give a mean atmospheric rotation rate of 1.40 ± 0.05 rev/day at a mean height near 240 km, for dates between July and December 1970, and local times ranging from 1800 hr to midnight to 0900 hr. This value is higher than those obtained from other satellites at similar heights.     

Analysis of the orbit of 1970-65D,cosmos 359 rocket

Cosmos 359 rocket 1970-65D, was launched on 22 August 1970 into an orbit inclined at 51·2° to the Equator, with an initial perigee height of 209 km: it decayed on 6 October 1971 after a lifetime of 410 days. The orbit has been determined at 42 epochs during the lifetime, using the RAE orbit refinement program, PROP, with over 2600 observations. Observations from the Hewitt cameras at Malvern and Edinburgh were available for 10 of the 42 orbits.Ten values of density scale height, at heights between 185 and 261 km, have been determined from analysis of the variations in perigee height.Upper-atmosphere zonal winds and 15th-order harmonics in the geopotential have been evaluated from the changes in orbital inclination. The average atmospheric rotation rate, for heights near 220 km, is found to be 1·04 rev/day; but there are striking departures from the average, with well-established values of 1·30, 0·75, 1·35 and 0·95 over four successive 75-day intervals. The changes in inclination at the 15th-order resonance in November 1970 give values of lumped 15th-order harmonics, which will provide equations for evaluating coefficients of order 15 and even degree (16,18,…) and also show that useful results on the geopotential can be obtained from satellites with perigee as low as 200 km.     

Air density at heights near 400km from the orbit of 1963-27A

The Agena B upper-stage rocket 1963-27A was launched into a near-circular orbit, inclined at 82.3° to the Equator, on 29 June 1963. Its orbit is determined at 52 epochs over the 16 month interval prior to its decay on 26 October 1969. The resulting orbital elements are used to obtain 95 atmospheric density values, at heights near 400km. Corrected to fixed heights, and normalised to a common exospheric temperature, these values reveal the semi-annual variation in density. A comparison between the observed variation and that of a recent model atmosphere is made. Although agreement between the two is generally good, their principal differences are discussed.     

The helium drag effect on the orbit of Explorer 19

The orbit of Explorer 19 (1963-53A) has been determined at 60 epochs between February 1976 and October 1976 from over 3000 observations. Using values of the orbital decay rate corrected for the effects of solar radiation pressure, 58 values of air density at a height of 900 km have been evaluated. After correcting for solar and geomagnetic activity and seasonal-latitudinal and diurnal variations in the exospheric temperature, the residual variation exhibited modulations associated with the ‘winter helium bulge’.An examination of three different models of the helium variation has indicated a procedure, which combines distinct features of the CIRA (1972) and Jacchia (1977) model atmospheres, for determining the atmospheric drag effect on Explorer 19. It is proposed that this technique may be equally applicable to any satellite in near-polar orbit at an equivalent height.     

Analysis of the orbital inclination of 1963-27A

The Agena B upper-stage rocket 1963-27A was launched into near-circular orbit, inclined at 82.3° to the Equator, on 29 June 1963. Its orbital elements were available at 52 epochs over the 16 month interval prior to its decay on 26 October 1969. During this period the satellite passed through 31:2 resonance and the variation in orbital inclination near this event was analysed to obtain two lumped geopotential harmonics of order 31. Since the resonance was a weak feature in the data, the resulting values are poorly defined.Either side of the resonance period, the inclination was used to estimate the mean atmospheric rotation rate Λ rev day^?1. The values obtained were Λ = 0.85 ± 0.18 at a height of 440 km for the period June 1968 to February 1969 and Λ = 1.13 ± 0.10 at 338 km for the period June to October 1969.     

Analysis of the orbit of 1970-97b (cosmos 378 rocket)

Cosmos 378 rocket, 1970-97B, entered orbit on 17 November 1970, with orbital inclination 74.0°, period 105 min and perigee height 230 km, and decayed on 30 September 1972 after 683 days in orbit. The RAE computer program PROP was used, with more than 1900 observations from 64 stations, to determine the orbit at 39 epochs between February 1971 and September 1972.The main aim of the analysis was to determine the atmospheric rotation rate from the decrease in orbital inclination, which was determined with a mean standard deviation of 0.0010° and a best standard deviation of 0.0003°. After removal of relevant perturbations, analysis of the variation in inclination between July 1971 and April 1972 yields the surprisingly low average atmospheric rotation rate of 0.75 ± 0.05 rev/day, at a mean height of 250 km. The local time at perigee is however strongly biassed towards daytime values (07–16 hr), so the results lend support to the picture of east-to-west winds by day and west-to-east winds by night.Values of scale height are obtained by analysis of the change in perigee height.     

The orbit of 1972-05b in its final phase,with geophysical inferences

The polar orbit of HEOS 2 second-stage rocket, 1972-05B, has been determined on each of the final 16 days before its decay in September 1978, using the RAE orbit refinement program, PROP 6, with about 1360 observations. An accuracy of 30–70 m, both radial and across track, was achieved.Eleven values of density scale height have been determined from the decrease in perigee height, with a 2% error; seven of these values are within 6% of the CIRA 1972 reference-atmosphere values, the rms value being 4% higher than CIRA.The rotation rate of the upper atmosphere, A, was determined from the decrease in orbital inclination as Λ = 1.40 ± 0.05 rev day^?1; i.e. a strong west-to-east zonal wind of 160 ± 20m s^?1, at a mean height of about 240 km. The local time was 01–02 h; solar activity was high; and the latitude of perigee moved steadily from 10°N to 67°S.     

Some aspects of a global thermospheric density model deduced from the analysis of the orbit of Intercosmos 13 rocket (1975-22B)

An empirical density formula is explored as a practical model for atmospheric variations and satellite drag analyses. Expanding neutral air density as a series of spherical harmonics and normalizing to a fixed height, an analytical expression for the rate of change of the mean motion is developed for an oblate atmosphere with density scale height varying linearly with altitude. A subset of the coefficients in the density expansion is determined by least-squares adjustment to the observed orbital decay rate of Intercosmos 13 rocket (1975-22B) for the period May 1975–December 1979. Comparisons against four thermospheric models are undertaken for the solar activity effect and the diurnal and semi-annual variations. Given the even spread of data and the increase in solar activity from low to moderate, the air density variation with solar activity is particularly well determined. The results support the “J77” model revealing a greater increase in density with the daily solar index than either the “MSIS” or “DTM” thermospheric models near the solar minimum. Analyses of the diurnal and semi-annual variations are less exact.     

Variations in air density,satellite drag coefficient and atmospheric rotation rate from analysis of the orbit of 1966-92D

Variations in air density, the satellite drag coefficient, and the atmospheric rotation rate at 60°S lat and 120–130 km height during the period September 1968–June 1969 have been determined from analysis of the high-eccentricity orbit of the 4th Molyniya 1 upper-stage rocket body, 1966-92D. The results show good correlation between density increases and strong geomagnetic activity, although solar flares of equal geomagnetic index value do not consistently produce density changes of equal magnitude. A 30 per cent semi-annual variation was observed, but there was no indication of the 50 per cent lower thermosphere seasonal-latitudinal variation that was predicted from the CIRA 1972 atmosphere. The satellite drag coefficient was observed to begin decreasing with height at an altitude where the molecular mean free path, λ, was twice the satellite's length. The coefficient decreased to a value approaching 1.0 as the satellite's perigee height fell below the altitude where λ was one-half the length. A mean atmospheric rotation rate of 1.1 ± 0.1 Earth rot/day was obtained for the last 20 days of decay. However, variations were observed with west-to-east wind speeds of ?100 m/sec measured for a local time of 13 hr.     

Analysis of the orbit of 1970-114F in its last 20 days

The satellite 1970-114F, the final-stage rocket of the Molniya 1S communications satellite, decayed in the atmosphere on 3 March 1973. During the last 20 days of its life the orbit suffered exceptionally rapid decay, with the apogee height decreasing from 7000 to 1000 km while the perigee height remained near 110 km. About 650 observations, made by visual observers in Britain and by U.S. Navy sensors, have been used with the PROP6 orbit refinement program to determine orbits at 14 epochs. Although the decay rate was more than ten times greater than in any previous orbit determination with PROP, good orbits were obtained, the standard deviation in inclination being less than 0.002° on eight orbits.The combination of high drag and good accuracy allows three techniques in orbital analysis to be successfully applied for the first time. Since zonal winds have little effect on the orbit, the changes in inclination are analysed to determine meridional winds near perigee, at heights of 110–120 km, latitudes of 63–65°S, and 6–12 hr LT. The changes in right ascension of the node are also successfully analysed for the same purpose. The two methods agree in indicating a south-to-north wind of 40 ± 30 m/sec from 11 to 21 February, a geomagnetically quiet period, and a south-to-north wind averaging 150 ± 30 m/sec from 22 February to 3 March, a geomagnetically disturbed period. Thirdly, the changes in the argument of perigee are analysed to determine atmospheric oblateness, which is found to be equal to the Earth's oblateness, to within ±20%. Lastly, the drag coefficient in transition flow is evaluated and found to be 0.85 ± 0.20.     

Analysis of thermospheric density variations neglected in modern atmospheric models using accelerometer data

The accuracy of the atmospheric density and its spatiotemporal variations given by the NRLMSISE-00 atmosphere model at the solar minimum are estimated using density measurements of the CACTUS microaccelerometer at heights of 270–600 km. The model errors are found to be noticeably (by a factor of 2–3) higher than the errors in atmospheric densities obtained from satellite drag data at solar minima. Microaccelerometer density data are used to study short-period (during one orbit) spatiotemporal density variations. The analysis of density variations over one orbit reveals orographic and continental effects. The amplitudes of the continental and orographic effects are estimated at 10–15% at a height of 270 km and 40% at a height of 600 km.     

Air density at heights near 300 km,from analysis of the orbit of China 2 rocket (1971-18b)

China 2 rocket, 1971-18B, was launched on 3rd March 1971 into an orbit inclined at 69.9° to the Equator, with an initial perigee height of 265 km. Analysis of its orbit has yielded values of air density at average intervals of 6 days between July 1971 and January 1972. When corrected to a fixed height, the density exhibits a correlation with the geomagnetic index A_p and the solar 10.7-cm radiation. With values of density extending over seven months it is possible to examine a complete cycle of the semi-annual variation at a height near 300 km. The values of density, corrected for the day-to-night variation and for solar and geomagnetic activity, reveal minima in mid-August and late January; at the intervening maximum, in early November, the density is almost 40% higher than at the minima.     

The aries auroral modelling campaign: characterization and modelling of an evening auroral arc observed from a rocket and a ground-based line of meridian scanners

An auroral arc system excited by soft electrons was studied with a combination of in situ rocket measurements and optical tomographic techniques, using data from a photometer on a horizontal, spinning rocket and a line of three meridian scanning photometers. The ground-based scanner data at 4709, 5577, 8446 and 6300 Å were successfully inverted to provide a set of volume emission rate distributions in the plane of the rocket trajectory, with a basic time resolution of 24 s. Volume emission rate profiles, derived from these distributions peaked at about 150 km for 5577 and 4709 Å, while the 8446 Å emission peaked at about 170 km with a more extended height distribution. The rocket photometer gave comparable volume emission rate distributions for the 3914 Å emission as reported in a separate paper by McDade et al. (1991, Planet. Space Sci. 39, 895). Instruments on the rocket measured the primary electron flux during the flight and, in particular, the flux precipitating into the auroral arc overflown at apogee (McEwen et al., 1991; in preparation). The local electron density and temperature were measured by probes on the rocket (Margot and McNamara (1991; Can. J. Phys. 69, 950). The electron density measurements on the downleg were modelled using ion production rate data derived from the optical results. Model calculations of the emission height profile based on the measured electron flux agree with the observed profiles. The height distribution of the N₂⁺ emission in the equatorward band, through which the rocket passed during the descent, was measured by both the rocket and the ground-based tomographic techniques and the results are in good agreement. Comparison of these profiles with model profiles indicates that the exciting primary spectrum may be represented by an accelerated Maxwellian or a Gaussian distribution centered at about 3 keV. This distribution is close to what would be obtained if the electron flux exciting the poleward form were accelerated by a 1–2 kV upward potential drop. The relative height profiles for the volume emission rate of the 5577 Å OI emission and the 4709 Å N₂⁺ emission were almost indistinguishable from each other for both the forms measured, with ratios in the range 38–50; this is equivalent to I(5577)/I(4278) ratios of 8–10. The auroral intensities and intensity ratios measured in the magnetic zenith from the ground during the period before and during the rocket flight are consistent with the primary electron fluxes and height distributions measured from the rocket. Values of I(5577)/I(4278) in the range 8–10 were also measured directly by the zenith ground photometers over which the arc system passed. These values are slightly higher than those reported by Gattinger and Vallance-Jones (1972) and this may possibly indicate an enhancement of the atomic oxygen concentration at the time of the flight. Such an enhancement would be consistent with our result, that the observed values of I(5577) and I(8446) are also significantly higher than those modelled on the basis of the electron flux spectrum measured at apogee.     

Air density at heights near 200km from the orbit of 1970-65D

Cosmos 359 rocket, 1970-65D, entered orbit on 22 August 1970, with an initial perigee height of 209 km and inclination 51·2°, and decayed on 6 October 1971. Using the values of perigee height from RAE orbits and decay rates from USAF Spacetrack bulletins, 146 values of air density have been calculated between August 1970 and September 1971, mainly at heights between 180 and 230 km.On ten occasions in 1971 when there were substantial geomagnetic disturbances there were sudden increases in density, the largest being about 32 per cent.When the density was corrected to a fixed height and allowance was made for the day-tonight variation and the effects of solar activity, the dominant feature was a semi-annual variation, with maxima in density centred at 6 November 1970 and 7 April 1971, and minima centred at 5 January and 28 July 1971. The maxima in density are nearly equal and exceed the minima by about 50 per cent.     

Air density at heights near 435 km from the orbit of skylab 1 (1973-27A)

On February 8, 1974, Skylab 1 was manoeuvred into a near circular orbit of inclination 50.04° and perigee near 420km. Orbital parameters have been computed at forty-six epochs thereafter using all available observations. Using these orbital elements, supplemented by orbital decay rates derived from NORAD bulletins, 193 values of air density were determined between 23 February, 1974 and 11, August, 1976. Corrected to a fixed height and normalised with respect to exospheric temperature these values reveal the semi-annual variation, exhibiting maxima in March–April and October–November, and minima in January–February and July–August. For 1974–1976 the July minima are more pronounced than the January minima whilst the April and October maxima appear equal. Overall the variation is greater than that indicated by CIRA 1972.