The high latitude circulation and temperature structure of the thermosphere near solstice

The neutral gas temperature and circulation of the thermosphere are calculated for December solstice conditions near solar cycle maximum using NCAR's thermospheric general circulation model (TGCM). High-latitude heat and momentum sources significantly alter the basic solar-driven circulation during solstice. At F-region heights, the increased ion density in the summer hemisphere results in a larger ion drag momentum source for the neutral gas than in the winter hemisphere. As a result there are larger wind velocities and a greater tendency for the neutral gas to follow the magnetospheric convection pattern in the summer hemisphere than in the winter hemisphere. There is about three times more Joule heating in the summer than the winter hemisphere for moderate levels of geomagnetic activity due to the greater electrical conductivity in the summer E-region ionosphere.

The results of several TGCM runs are used to show that at F-region heights it is possible to linearly combine the solar-driven and high-latitude driven solutions to obtain the total temperature structure and circulation to within 10–20%. In the lower thermosphere, however, non-linear terms cause significant departures and a linear superposition of fields is not valid.

The F-region winds at high latitudes calculated by the TGCM are also compared to the meridional wind derived from measurements by the Fabry-Perot Interferometer (FPI) and the zonal wind derived from measurements by the Wind and Temperature Spectrometer (WATS) instruments onboard the Dynamics Explorer (DE−2) satellite for a summer and a winter day. For both examples, the observed and modeled wind patterns are in qualitative agreement, indicating a dominant control of high latitude winds by ion drag. The magnitude of the calculated winds (400–500 m s⁻¹) for the assumed 60 kV cross-tail potential, however, is smaller than that of the measured winds (500–800 m s⁻¹). This suggests the need for an increased ion drag momentum source in the model calculations due to enhanced electron densities, higher ion drift velocities, or some combination that needs to be further denned from the DE−2 satellite measurements.     

Hydrogen atoms and ions in the thermosphere and exosphere

The distribution of atomic hydrogen in the thermosphere and exosphere is computed taking into account the upward flow which balances the escape flux. Because of the upward flow the number-density gradient is much steeper than it would be in a static atmosphere. Attention is drawn to the fact that the ratio of the amount of hydrogen above the 100 or 110km levels to the amount of hydrogen above the 200 or 300 km levels is a sensitive measure of the temperature of the exosphere. The evidence on the absolute abundance of atomic hydrogen is examined. It is concluded that the number density at the 120km level is probably about 5 × 10⁵/cm³. The Ly. absorption line at this level is beyond the linear part of the curve of growth.

Consideration is also given to the steady-state distributions of O⁺ and H⁺ ions. In the lower part of the exosphere the number density of O⁺ ions falls with increase in altitude (the associated scale height being twice that of the O atoms) and the number density of H⁺ ions rises at the same rate (as was first pointed out by Dungey). The altitude at which the number densities of O⁺ and H⁺ ions become equal is calculated on various assumptions regarding the temperature and hydrogen content of the exosphere. It is found to be about 1200 km when the temperature is 1250° K and the hydrogen content corresponds to the number density cited near the end of the preceding paragraph. The gradient of the predicted electrondensity distribution at several Earth radii is much less than that deduced from whistler studies.

The passage from charge transfer to diffusive equilibrium is discussed in an Appendix.     

Annual and semi-annual zonal wind components and corresponding temperature and density variations, 60–130 km

From rocket and radar-meteor wind observations, annual and semi-annual components of the zonal flow are derived for latitudes N at heights between 60 and 130 km. Height regions of maximum and minimum amplitude are described with reference to changes in phase. The annual components decrease with height throughout the mesosphere and, after a reversal of phase, enhance to 25 m/sec at 100 ± 5 km. The semi-annual components have maximum amplitudes of 25 m/sec over a wide range of latitude in two height regions at 90 and 120 km and in a limited range of latitude (near 50°) at 65 km.

Calculated temperatures and log densities are discussed in terms of amplitude and phase as functions of height and latitude. Below 100 km a comparison is made with temperature amplitudes derived from independent temperature data. Above 100 km the annual temperature variation maximizes at 115 km and is particularly large at high latitudes (exceeding 50°K). On the other hand, the semi-annual component increases rapidly with height between 110 and 120 km at all latitudes maximizing at the 120 km level, where amplitudes exceed 25°K at high and low latitudes and 10°K at mid-latitudes. The annual component of log density, like the temperature variation, is largest at high latitudes up to 125 km. The semi-annual variation has a minimum at 110–115 km, above which amplitudes increase with height, reaching 5–12 per cent at 130 km according to latitude. The phases at and near 130 km for the annual and semi-annual density variations are very close to those found at greater heights from satellite orbits and amplitudes could be readily extrapolated to agree with those in the satellite region.     

Rocket measurements of the 6300 Å and 3914 Å dayglow features

Rocket results are presented on the OI 6300 Å line and on the N₂⁺ 3914 Å band in the dayglow. An altitude range of 78–335 km is covered. Theoretical interpretations are given, using results of simultaneous measurements of electron density and electron temperature. The apparent brightness of the 6300 Å line at the base of the emitting region is found to be 13 kR, of which 5.5 kR are ascribed to excitation through the Schumann-Runge dissociation of O₂ by the solar UV radiations, 0.55 kR to the dissociative recombination of O₂⁺ and NO⁺ ions, and 0.03 kR to the excitation of O by thermal electrons. An additional source of excitation above 280 km is suggested. The deactivation of O(¹D) by O₂(X³Σ_g⁻) is found to be appreciable below 200 km, and its rate coefficient is estimated to be 2 × 10⁻¹⁰ cm³/sec. The apparent brightness of the 3914 Å band at the base of the emitting region is found to be 6.5 kR, decreasing to 3.2 kR at 330 km. Assuming that fluorescent scattering of solar radiation is the mechanism involved the distribution of N₂⁺ ions is calculated. The rate coefficients for the loss of these ions are hence calculated.     

A study of F2 region night-time vertical ionization fluxes at millstone hill

Vertical fluxes of ionization in the F2 region have been measured by the incoherent scatter technique over Millstone Hill in 1969. The results obtained near midnight for the region above h_maxF2 have been examined to determine whether there is a significant flux of ionization from the magnetosphere to the ionosphere that serves to maintain the F-layer. It is found that H⁺ ions are a minor constituent over the altitude range in which useful measurements can be made, so that any conclusion must rest upon properly interpreting the observed O⁺ fluxes. By selecting periods when the layer did not appear to be decaying rapidly it was hoped to find cases where the O⁺ flux did not vary with altitude in the range 500 h 800 km (i.e. where losses are unimportant), since this would imply that the flux is of magnetospheric origin.

While three cases exhibited this behaviour, the majority exhibited a decrease in the O⁺ flux with height, indicating that the layer was descending. Attempts to correct for this were made, and the average flux from the magnetosphere was estimated as 3 × 10⁷ el/cm²/sec. This is in fair agreement with other recent estimates, and implies that at this latitude the ionosphere is not maintained solely by the magnetospheric flux. Moreover, large increases in flux that could give rise to nocturnal increases in the total content of the layer do not appear to have been seen.     

Some consequences of turbulent dissipation on diurnal thermal tide

It is shown that the phase difference between the horizontal components of the wind velocity induced by the first diurnal propagating tidal mode between 90 and 100 km altitude is modified by turbulent diffusion. Experimental data are compared with calculations made by two different methods. The calculated results agree well with the data for a coefficient of eddy diffusion K = 5 × 10⁶ cm² s⁻¹ which is a realistic one at this altitude. Incidentally, it is shown that some measurements which have been interpreted as a presence of evanescent modes are probably due to an increase in vertical wavelength of the first propagating diurnal mode due to turbulent diffusion.     

Atomic oxygen transport in the thermosphere

The photodissociation of oxygen in the lower thermosphere is evaluated to obtain its global average value and the hemispheric imbalance. The observed concentrations of atomic oxygen do not reflect this imbalance in production due to the effect of seasonal wind patterns redistributing the atomic oxygen. The wind system necessary to compensate for the imbalance in solar thermal input into the lower thermosphere is found to transport an amount of atomic oxygen sufficient to compensate for the hemispheric imbalance in production. Ionospheric data indicate a winter enhancement in atomic oxygen concentration; to produce this, a higher degree of oxygen dissociation than that normally accepted (i.e. higher than an atomic to molecular oxygen ratio of unity at 120 km) is needed. The concept that the concentrations of atomic oxygen observed over the winter polar region are maintained by transport from lower latitudes requires that eddy diffusion coefficients derived from vertical transport at low latitudes (ignoring horizontal transport) be reduced by about 25 per cent.     

An auroral infrasonic substorm investigation using Chatanika radar and other geophysical sensors

Characteristics of the supersonic auroral arcs within the 0905 UT 2 April 1973 substorm were determined using data from (1) all-sky cameras; (2), surface magnetometers, (3) multispectral scanning photometers, (4) 30MHz riometers, (5) Chatanika incoherent-scatter radar, (6) Homer auroral radar, and (7) infrasonic microphone arrays at College and Stevens Village in Alaska. These data were analyzed to determine the properties of an auroral electrojet arc that generates auroral infrasonic waves (AIW).

An arc that was show to be the source of an AIW was found to have the following characteristics: (1) a velocity of 500 m/sec traveling from an azimuth of 350°; (2) an intensity in 4278 A of 26 Kr, (3) a maximum electron density of 2.8 × 10⁶ el/cm⁶ at 100km height, (4) an equivalent westward line current of 2.8 × 10⁶ A, (5) orientation of ΔH parallel to the AIW direction of travel and perpendicular to the arc's long axis, (6) a characteristic energy of the primary auroral electron spectrum of 3.0keV, and (7) an energy deposition rate for the auroral pdarticles of 100 erg/cm² sec.     

The consequences of high latitude particle precipitation on global thermospheric dynamics

A previous comparison of experimental measurements of thermospheric winds with simulations using a global self-consistent three-dimensional time-dependent model confirmed a necessity for a high latitude source of energy and momentum acting in addition to solar u.v. and e.u.v. heating. During quiet geomagnetic conditions, the convective electric field over the polar cap and auroral oval seemed able to provide adequate momentum input to explain the thermospheric wind distribution observed in these locations. However, it seems unable to provide adequate heating, by the Joule mechanism, to complete the energy budget of the thermosphere and, more importantly, unable to provide the high latitude input required to explain mean meridional winds at mid-latitudes. In this paper we examine the effects of low energy particle precipitation on thermospheric dynamics and energy budget. Modest fluxes over the polar cap and auroral oval, of the order of 0.4 erg cm ⁻²/s, are consistent with satellite observations of the particles themselves and with photometer observations of the OI and OII airglow emissions. Such particle fluxes, originating in the dayside magnetosheath cusp region and in the nightside central plasma sheet, heat the thermosphere and modify mean meridional winds at mid-latitudes without enhancing the OI 557.7 line, or the ionization of the lower thermosphere (and thus enhancing the auroral electrojets), neither of which would be consistent with observations during quiet geomagnetic conditions.     

Electron precipitation observations from a rocket flight through the dayside auroral oval

Energy spectra of electrons between 30 eV and 18 keV were obtained with a spectrometer on a Black Brant rocket launched from Cape Parry, N.W.T. (Λ = 75.2°) on December 6, 1974 to study the dayside magnetospheric cleft. The rocket flew to an apogee of 236 km and travelled poleward to 80° invariant latitude. The cleft was observed to extend from 76.9 to 78.4° invariant latitude. Equatorward of this electrons of a few keV energy were observed with a total energy flux of up to 2 erg/cm² sec ster. Variable fluxes of electrons with a spectrum fitted by a Maxwellian distribution of 150 eV characteristic energy were observed through most of the cleft. One inverted V structure was crossed. In that region, the electron energy increased to 650 eV and a total energy flux of 8 erg/cm² sec ster was measured. The event was a temporal one and only a few km in width, as deduced from optical data. Fluxes of about 10⁻² erg/cm² sec ster were recorded poleward of the cleft.     

The spectrum of electron content fluctuations in the ionosphere

Continuous records of the electron content of the ionosphere, from 1965 to 1970, are used to obtain power spectra covering periods from 30 sec to 2 yr at latitudes of 34°S and 42°S. At periods up to 5 min, amplitudes were less than 0.2 per cent of the total electron content. Variations produced by gravity waves were very common in the range 20–80 min, with no preferred periods. The r.m.s. amplitude per octave A₀ was about 10¹⁵ electrons/m², or 0.6 per cent of the mean electron content. The amplitude increased during the day, particularly in winter when periodic components predominated. The cut-off at about 17 min was sharply defined, giving a mean scale height for the neutral atmosphere (at 300 km) of about 43 km in summer, 47 km on winter days and 42 km on winter nights.

From 12 hr to 1 month A₀ was about 12 per cent of the mean electron content in both summer and winter at 34°S, and 10 per cent at 42°S. The 24 hr and 27 day peaks were largest just before sunspot maximum, and almost disappeared near sunspot minimum. Variations between 1 and 27 days reflect the random occurrence of ionospheric storms and show no consistent peaks. Day to day and night to night variations were both about 10 per cent of the background content for periods from 2 days to 2 yr, apart from a slight decrease between 1 and 6 months.     

Properties of energetic water-group ions in the extended pick-up region surrounding comet Giacobini-Zinner

The properties of energetic (65–95 keV) cometary water-group ions in the extended solar wind pick-up region surrounding comet Giacobini-Zinner are examined using data from the EPAS instrument on the ICE spacecraft. In the outer part of this region, extending from cometocentric distances of several hundred thousand to a few million kilometres (the limit of pick-up ion detectability), it is found that large modulations of the ion flux occur (with J_MAX/J_MIN 10²-10³) which are related to the direction of the magnetic field. It is also found that the ions stream in a direction which is intermediate between the directions of the solar wind flow and the E × B drift, and that ions are present at energies somewhat above the local pick-up energy. These properties indicate that the waves which are excited by the unstable “ring-beam” pick-up ion velocity distributions do result in significant scattering of the ions in this region, both in pitch angle and in energy, but that they have insufficient amplitude to scatter the ions into near isotropy in the solar wind frame. Closer to the comet (but still upstream from the bow shock), the ion flux modulations are considerably reduced in amplitude and the ions respond less to the E × B drift, indicating that the ions are scattered nearer to isotropy in this region. Inbound, this transition takes place relatively abruptly at a distance of 4 × 10⁵ km in association with an increase in the solar wind speed, after which the ion flux increases, and ceases to be modulated by the field direction, while the streaming direction is continuously antisolar and unmodulated by the direction of the E × B drift. Outbound, weak vestiges of the ring-beam ion anisotropy are present in the region immediately upstream from the bow shock (at −1 × 10⁵ km), but these become more marked at distances in excess of t4 × 10⁵ km, increasing gradually with increasing distance from the comet. It is shown that the evolution of the ion properties is qualitatively consistent with expectations based on quasi-linear diffusion of the ions by the magnetosonic waves observed during the encounter.     

An auroral F-region study using in situ measurements by the Atmosphere Explorer-C satellite

On 14 July 1974 the Atmosphere Explorer-C satellite flew through an aurora at F-region altitudes just after local midnight. The effects of the particle influx are clearly evident in the ion densities, the 6300 Å airglow, and the electron and ion temperatures. This event provided an opportunity to study the agreement between the observed ion densities and those calculated from photochemical theory using in situ measurements of such atmospheric parameters as the neutral densities and the differential electron energy spectra obtained along the satellite track. Good agreement is obtained for the ions O₂⁺, NO⁺ and N₂⁺ using photochemical theory and measured rate constants and electron impact cross sections. Atomic nitrogen densities are calculated from the observed [NO⁺]/[O₂⁺] ratio. In the region of most intense electron fluxes (20 erg cm⁻² sec⁻¹) at 280 km, the N density is found to be between 2 and 7 × 10⁷ cm⁻³. The resulting N densities are found to account for approx. 60% of the production of N⁺ through electron impact on N and the resonant charge exchange of O⁺(²P) with N(⁴S). This reaction also provides a significant source of O(¹S) in the aurora at F-region altitudes. In the region of intense fast electron influx, the reaction with atomic nitrogen is found to be the main loss of O⁺(²P).     

Studies of equatorial 630.0 nm airglow enhancements produced by a chemical release in the F-region

A sky-mapping filter photometer has been used to determine the 630.0 nm airglow enhancement produced by explosive release of 3 × 10²⁶ CO₂ molecules into the F-region at 320 km altitude on 8 September 1982 as part of project BIME. The enhancement is produced when COg molecules engage in atom transfer with the F-region O⁺ ions to form O₂⁺ ions, which subsequently dissociatively recombine with the ambient electrons to produce O(¹D) atoms to yield the 630.0 nm radiation. The morphology of the enhanced airglow region has been traced in a series of 630.0 nm intensity contour maps as a function of time, the enhancement reaching a central brightness of approximately 400 R about 2 min after release and a diameter of 250 km some 3 min after release. The measurements of central intensity and enhanced region radius as a function of time are compared with model calculations by Mendillo and Herniter of diffusive expansion of CO₂ molecules from either a point release or from an initial, extended volume. While peak intensities are reasonably reproduced, the measured decay of the 630.0 nm intensity and the growth in size of the enhanced region are rather different from the model predictions. The measured 200 m/s drift southeastward of the enhanced region is consistent with the motion of the neutral thermosphere determined from optical doppler shifts less than an hour earlier.     

Response of mesopheric ozone to particle precipitation

In the mesosphere, water vapor photolysis is the major source of odd hydrogen (H, OH and HO₂) under normal conditions. The odd hydrogen produced may then be converted to H₂ by the reaction H + HO₂→ H₂ + O₂. This process is responsible for the calculated decrease in the H₂O mixing ratio and accompanying increase in the H₂ mixing ratio with altitude in the upper mesosphere and lower thermosphere. Charged particle precipitation events are calculated to produce the same effect, particularly in the 70–85 km region, thus temporarily resulting in enhanced conversion of H₂O to H₂ following such an event. Since odd hydrogen is produced predominantly by water vapor photolysis at these altitudes, decreased odd hydrogen concentrations are also anticipated. Odd hydrogen processes dominate ozone destruction in this region, and so an increase in ozone may occur if odd hydrogen concentrations decrease. We have examined the calculated time behavior of these processes in a numerical model using the August 1972 solar proton event as an example, and we present calculations indicating what might be observed in future events.     

The composition,structure, temperature and dynamics of the upper thermosphere in the polar regions during October to December 1981

During the period October to December 1981, the Dynamics Explorer-2 (DE-2) spacecraft successively observed the South polar and the North polar regions, and recorded the temperature, composition and dynamical structure of the upper thermosphere. In October 1981, perigee was about 310 km altitude, in the vicinity of the South Pole, with the satellite orbit in the 09.00–21.00 L.T. plane. During late November and December, the perigee had precessed to the region of the North Pole, with the spacecraft sampling the upper thermosphere in the 06.00 18.00 L.T. plane. DE-2 observed the meridional wind with a Fabry-Perot interferometer (FPI), the zonal wind with the wind and temperature spectrometer (WATS), the neutral temperature with the FPI, and the neutral atmosphere composition and density with the neutral atmosphere composition spectrometer (NACS). A comparison between the South (summer) Pole and the North (winter) Pole data shows considerable seasonal differences in all neutral atmosphere parameters. The region of the summer pole, under similar geomagnetic and solar activity conditions, and at a level of about 300 km, is about 300 K warmer than that of the winter pole, and the density of atomic oxygen is strongly depleted (and nitrogen enhanced) around the summer pole (compared with the winter pole). Only part of the differences in temperature and composition structure can be related to the seasonal variation of solar insolation, however, and both polar regions display structural variations (with latitude and Universal Time) which are unmistakeable characteristics of strong magnetospheric forcing. The magnitude of the neutral atmosphere perturbations in winds, temperature, density and composition within both summer and winter polar regions all increase with increasing levels of geomagnetic activity.The UCL 3-dimensional time dependent global model has been used to simulate the diurnal, seasonal and geomagnetic response of the neutral thermosphere, attempting to follow the major features of the solar and geomagnetic inputs to the thermosphere which were present during the late 1981 period.In the UCL model, geomagnetic forcing is characterized by semi-empirical models of the polar electric field which show a dependence on the Y component of the Interplanetary Magnetic Field, due to Heppner and Maynard (1983), It is possible to obtain an overall agreement, in both summer and winter hemispheres, with the thermospheric wind structure at high latitudes, and to explain the geomagnetic control of the combined thermal and compositional structure both qualitatively and quantitatively. To obtain such agreement, however, it is essential to enhance the polar ionosphere as a consequence of magnetospheric particle precipitation, reflecting both widespread auroral (kilovolt) electrons, and “soft” cusp and polar cap sources. Geomagnetic forcing of the high latitude thermosphere cannot be explained purely by a polar convective electric field, and the thermal as well as ionising properties of these polar and auroral electron sources are crucial components of the total geomagnetic input.     

Measurements of attenuation and absorption lengths with the KASCADE experiment

The attenuation of the electron shower size beyond the shower maximum is studied with the KASCADE extensive air shower (EAS) experiment in the primary energy range of about 10¹⁴–10¹⁶ eV. Attenuation and absorption lengths are determined by applying different approaches, including the method of constant intensity, the decrease of the flux of EASs with increasing zenith angle, and its variation with ground pressure. We observe a significant dependence of the results on the applied method. The determined values of the attenuation length ranges from 175 to 196 g/cm² and of the absorption length from 100 to 120 g/cm². The origin of these differences is discussed emphasizing the influence of intrinsic shower fluctuations.     

Analysis of the variation in inclination of Intercosmos 13 Rocket, 1975-22B

The orbit of Intercosmos 13 rocket (1975-22B) has been determined at 103 epochs between 30 April 1975 and 10 April 1980 from almost 7000 observations. One hundred and three values of inclination have been determined and corrections incoporated for the effects due to zonal harmonic, lunisolar and tesseral harmonic perturbations, precession, and solid Earth tides. The modified data have been analysed to yield values of the atmospheric rotation rate, Λ rev day⁻¹, viz. Λ = 0.94 ± 0.10 at an average height of 322 ± 6 km and Λ = 1.27 ± 0.02 at 288 km. Analysis of the inclination near 14th-order resonance has indicated lumped harmonic values 10^{9 1.0}_1.4 = − 76.13 ± 12.47, 10^{9 1,0}₁₄ = − 29.89 ± 32.64, 10^{9 −1.2}₁₄ = − 63.11 ± 15.44 10^{9 −1.2}₁₄ = − 32.52 ± 26.96, for inclination 82.952°.     

Mesospheric O3, H and H2O at high latitudes: A theoretical model

The high latitude thermosphere is characterized by a large heat input which produces a strong wind field. In a previous work, the modification of the vertical transport mechanisms produced by high latitude horizontal winds was studied and the resulting concentration profiles are used here to study the influence at mesospheric levels. We obtain an improved agreement with experimental measurements. High solar zenithal angles could explain other differences between the high and middle latitude's mesosphere, especially below 80 km, approximately.     

A zonal-averaged dynamical model for the middle atmosphere including gravity wave mean flow interaction: Solstice conditions

Using a set of transformed Eulerian equations the zonal-averaged circulation of the middle atmosphere (10–110 km) is calculated on a global scale for solstice conditions. The emphasis lies on an improved modelling of the zonal momentum balance of the mesophere and lower thermosphere. For this purpose an internal gravity wave mean flow interaction model suggested by T. Matsuno has been incorporated in a slightly modified version. With this model the observed reversal of the zonal wind with height in the summer upper mesosphere and lower thermosphere can be reproduced. The coefficient of eddy momentum diffusion and the Rayleigh friction coefficient used in this model have been made temperature dependent by describing them as a function of the local static stability parameter.