The generation of vertical thermospheric winds and gravity waves at auroral latitudes—II. Theory and numerical modelling of vertical winds

Recent observations of strong vertical thermospheric winds and the associated horizontal wind structures, using the 01(3P-1D)nm emission line, by ground-based Fabry-Perot interferometers in Northern Scandinavia have been described in an accompanying paper (Paper I). The high latitude thermosphere at a height of 200–300 km displays strong vertical winds (30–50m ms^?1)of a persistent nature in the vicinity of the auroral oval even during relatively quiet geomagnetic conditions. During an auroral substorm, the vertical (upward) wind in the active region, including that invaded by a Westward Travelling Surge, may briefly(10–30 min)exceed 150 m s^?1. Very large and rapid changes of horizontal wind structure (up to 500 m^?1 in 30 min) usually accompany such large impulsive vertical winds. Magnetospheric energy and momentum sources generate large vertical winds of both a quasi-steady nature and of a strongly time-dependent nature. The thermospheric effects of these sources can be evaluated using the UCL three-dimensional, time-dependent thermospheric model. The auroral oval is, under average geomagnetic conditions, a stationary source of significant vertical winds (10–40 m s^?1). In large convective events (directly driven by a strong momentum coupling from the solar wind) the magnitude may increase considerably. Auroral substorms and Westward Travelling Surges appear to be associated with total energy disposition rates of several tens to more than 100 erg cm^?2s^?1, over regions of a few hours local time, and typically 2–5° of geomagnetic latitude (approximately centred on magnetic midnight). Such deposition rates are needed to drive observed time-dependent vertical (upward) winds of the order of 100–200m s^?1.The response of the vertical winds to significant energy inputs is very rapid, and initially the vertical lifting of the atmosphere absorbs a large fraction (30% or more) of the total substorm input. Regions of strong upward winds tend to be accompanied in space (and time) by regions of rather lower downward winds, and the equatorward propagation of thermospheric waves launched by auroral substorms is extremely complex.     

The generation of vertical thermospheric winds and gravity waves at auroral latitudes—I. Observations of vertical winds

Observations of vertical and horizontal thermospheric winds, using the OI (3P-1D) 630 nm emission line, by ground-based Fabry-Perot interferometers in Northern Scandinavia and in Svalbard (Spitzbergen) have identified sources of strong vertical winds in the high latitude thermosphere. Observations from Svalbard (78.2N 15.6E) indicate a systematic diurnal pattern of strong downward winds in the period 06.00 U.T. to about 18.00 U.T., with strong upward winds between 20.00 U.T. and 05.00 U.T. Typical velocities of 30 m s^?1 downward and 50 m s^?1 upward occur, and there is day to day variability in the magnitude (30–80 m s^?1) and phase (+/- 3 h) in the basically diurnal variation. Strong and persistent downward winds may also occur for periods of several hours in the afternoon and evening parts of the auroral oval, associated with the eastward auroral electrojet (northward electric fields and westward ion drifts and winds), during periods of strong geomagnetic disturbances. Average downward values of 30–50 m s^?1 have been observed for periods of 4–6 h at times of large and long-lasting positive bay disturbances in this region. It would appear that the strong vertical winds of the polar cap and disturbed dusk auroral oval are not in the main associated with propagating wave-like features of the wind field. A further identified source is strongly time-dependent and generates very rapid upward vertical motions for periods of 15–30 min as a result of intense local heating in the magnetic midnight region of the auroral oval during the expansion phase of geomagnetic disturbances, and accompanying intense magnetic and auroral disturbances. In the last events, the height-integrated vertical wind (associated with a mean altitude of about 240 km) may exceed 100–150 m s^?1. These disturbances also invariably cause major time-dependent changes of the horizontal wind field with, for example, horizontal wind changes exceeding 500 m s^?1 within 30 min. The changes of vertical winds and the horizontal wind field are highly correlated, and respond directly to the local geomagnetic energy input. In contrast to the behaviour observed in the polar cap or in the disturbed afternoon auroral oval, the ‘expansion phase’ source, which corresponds to the classical ‘auroral substorm’, generates strong time-dependent wind features which may propagate globally. This source thus directly generates one class of thermospheric gravity waves. In this first paper we will consider the experimental evidence for vertical winds. In a second paper we will use a three-dimensional time-dependent model to identify the respective roles of geomagnetic energy and momentum in the creation of both classes of vertical wind sources, and consider their propagation and effects on global thermospheric dynamics.     

Magnetic ordering of the polar airglow

The visible airglow experiment on the Atmosphere Explorer-C satellite has gathered sufficient data over the Earth's polar regions to allow one to map the geographic distribution of particle precipitation using emissions at 3371 and 5200 Å. Both of these features exhibit large variations in space and time. The 3371 Å emission of N₂(C³π), excited by low energy electrons, indicates substantial energy inputs on the dayside in the vicinity of the polar cusp. More precipitation occurs in the morning than evening for the sample reported here, while the entire night sector between magnetic latitudes 65° and 77.5° is subjected to particle fluxes. Regions of enhanced 5200 Å emission from $\text{N}\text{(}{}^2\text{D)}$ are larger in horizontal extent than those at 3371 Å. This smearing effect is due to ionospheric motions induced by magnetospheric convection.     

A theoretical and empirical study of the response of the high latitude thermosphere to the sense of the “Y” component of the interplanetary magnetic field

The strength and direction of the Interplanetary Magnetic Field (IMF) controls the transfer of solar wind momentum and energy to the high latitude thermosphere in a direct fashion. The sense of “ Y” component of the IMF (BY) creates a significant asymmetry of the magnetospheric convection pattern as mapped onto the high latitude thermosphere and ionosphere. The resulting response of the polar thermospheric winds during periods when BY is either positive or negative is quite distinct, with pronounced changes in the relative strength of thermospheric winds in the dusk-dawn parts of the polar cap and in the dawn part of the auroral oval. In a study of four periods when there was a clear signature of BY, observed by the ISEE-3 satellite, with observations of polar winds and electric fields from the Dynamics Explorer-2 satellite and with wind observations by a ground-based Fabry-Perot interferometer located in Kiruna, Northern Sweden, it is possible to explain features of the high latitude thermospheric circulation using three dimensional global models including BY dependent, asymmetric, polar convection fields. Ground-based Fabry-Perot interferometers often observe anomalously low zonal wind velocities in the (Northern) dawn auroral oval during periods of extremely high geomagnetic activity when BY is positive. Conversely, for BY negative, there is an early transition from westward to southward and eastward winds in the evening auroral oval (excluding the effects of auroral substorms), and extremely large eastward (sunward) winds may be driven in the auroral oval after magnetic midnight. These observations are matched by the observation of strong anti-sunward polar-cap wind jets from the DE-2 satellite, on the dusk side with BY negative, and on the dawn side with BY positive.     

The westward thermospheric jet-stream of the evening auroral oval

One of the most consistent and often dramatic interactions between the high latitude ionosphere and the thermosphere occurs in the vicinity of the auroral oval in the afternoon and evening period. Ionospheric ions, convected sunward by the influence of the magnetospheric electric field, create a sunward jet-stream in the thermosphere, where wind speeds of up to 1 km s^?1 can occur. This jet-stream is nearly always present in the middle and upper thermosphere (above 200 km altitude), even during periods of very low geomagnetic activity. However, the magnitude of the winds in the jet-stream, as well as its location and range in latitude, each depend on geomagnetic activity. On two occasions, jet-streams of extreme magnitude have been studied using simultaneous ground-based and satellite observations, probing both the latitudinal structure and the local time dependence. The observations have then been evaluated with the aid of simulations using a global, three-dimensional, time-dependent model of thermospheric dynamics including the effects of magnetospheric convection and particle precipitation. The extreme events, where sunward winds of above 800 ms^?1 are generated at relatively low geomagnetic latitudes (60–70°) require a greatly expanded auroral oval and large cross-polar cap electric field ( ~ 150 kV). These in turn are generated by a persistent strong Interplanetary Magnetic Field, with a large southward component. Global indices such as K_p are a relatively poor indicator of the magnitude and extent of the jet-stream winds.     

Measurements of thermospheric response to auroral activities

The Joule heating produced by auroral electrojets and its thermospheric response can be studied by monitoring the thermospheric temperatures by means of optical methods; simultaneously investigating the concurrent auroral electrojet activities using geomagnetic records obtained from stations along a meridian close to the observation site of optical measurements. We report, in this paper, the measurements of thermospheric response to auroral activities which were made at Albany (42.68°N, 73.82°W), New York on 2 September 1978 (U.T.) when an isolated substorm occurred. The thermospheric temperatures were measured by using a high-resolution Fabry-Perot interferometer that determines the line profiles of the [OI] 6300 Å line emission. The intensities and latitudinal positions of auroral electrojets were obtained by the analysis of magnetograms from the IMS Fort Churchill meridian chain stations.     

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.     

Observations of the optical spectrum of the dayside magnetospheric cleft aurora

Optical spectra of the cleft aurora in the region 5000–8500 Å were measured in December, 1977 at Cape Parry, N.W.T. A Michelson interferometer was used at a resolution of 10 cm^?1. The auroral features observed were OI (5577, 6300-64, 7774, 8446 Å), OII (7319-30 Å), NI (5200 Å), Hα, O₂ atm (1,1), some weak N₂1P bands and possibly some Meinel bands of N₂⁺. In addition, nightglow emissions of Na and OH were observed. Theoretical predictions of the OI and NI emission rates using the model of Link et al. (1980) fit the observed rates reasonably well if a 40 eV Maxwellian incident electron spectrum is assumed. The predicted rates for OII exceed the observed value by a factor of 4. It is suggested that the ionization cross-section may be over-estimated.     

A global circulation model of Saturn's thermosphere

We present the first 3-dimensional self-consistent calculations of the response of Saturn's global thermosphere to different sources of external heating, giving local time and latitudinal changes of temperatures, winds and composition at equinox and solstice. Our calculations confirm the well-known finding that solar EUV heating alone is insufficient to produce Saturn's observed low latitude thermospheric temperatures of 420 K. We therefore carry out a sensitivity study to investigate the thermosphere's response to two additional external sources of energy, (1) auroral Joule heating and (2) empirical wave heating in the lower thermosphere. Solar EUV heating alone produces horizontal temperature variations of below 20 K, which drive horizontal winds of less than 20 m/s and negligible horizontal changes in composition. In contrast, Joule heating produces a strong dynamical response with westward winds comparable to the sound speed on Saturn. Joule heating alone, at a total rate of 9.8 TW, raises polar temperatures to around 1200 K, but values equatorward of 30° latitude, where observations were made, remain below 200 K due to inefficient meridional energy transport in a fast rotating atmosphere. The primarily zonal wind flow driven by strong Coriolis forces implies that energy from high latitudes is transported equatorward mainly by vertical winds through adiabatic processes, and an additional 0.29-0.44 mW/m² thermal energy are needed at low latitudes to obtain the observed temperature values. Strong upwelling increases the H₂ abundances at high latitudes, which in turn affects the H⁺₃ densities. Downwelling at low latitudes helps increase atomic hydrogen abundances there.     

On the temperature dependence of the reaction O + NO → NO12

The published data on the temperature dependence of the radiative combination of atomic oxygen with nitric oxide at pressures near 1 torr is examined. Arguments are advanced to suggest that radiation near the cut-off wavelength (～ 3875Å) is coming from the unstabilized activated complex, No¹₂. At 4000Å a positive activation energy of 1 kcal mole^?1 is deduced. Application of this temperature dependence with the rate coefficient at 5200Å is made to airglow measurements in aurora. The deduced NO concentration is about 10⁹ cm^?3, in general agreement with that deduced from the measured NO⁺/O⁺₂ ratio as well as an auroral model prediction.     

Auroral heating and the composition of the neutral atmosphere

Heating of the neutral atmosphere by auroral particle fluxes and by orthogonal electric fields is responsible for large changes in the thermospheric composition that have been observed by satellite mass spectrometers. Vertical winds of a few meters per second are produced in the region subject to auroral heating; this vertical upwelling drives circulation cells that extend the effects of heating in the auroral region on a global scale. Our analysis focuses on the initial phase of a magnetic storm within the auroral region.     

Upper atmospheric temperatures from Doppler line widths—V. Auroral electron energy spectra and fluxes deduced from the 5577 and 6300 Å atomic oxygen emissions

This is a report upon further data obtained from the auroral OI 5577 Å emission with a Wide Angle Michelson Interferometer (WAMI), and upon our first observations made with it on the 6300 Å emission. The method used for converting emission intensities and temperatures to auroral electron fluxes and energy spectra is described. Data for the 5577 Å emission are presented for the (lack of) heating in auroral forms, vertical temperature profiles in aurora, electron flux and energy spectrum variations in pulsating aurora, and a ‘cold’ subvisual auroral arc. Data from the OI 6300 Å emission are presented for the diurnal variation of exospheric temperature and for the thermalization of $\text{O}\text{(}{}^1\text{D)}$ in the F-region.     

Diurnal variability of thermospheric N and NO

A model for diurnal variations of neutral and ionic nitrogen compounds in the thermosphere is reconstructed on the basis of a new photochemical aspect on N(²D), together with new observations of the NO density. The NO density so far measured must be reduced by a factor 2, due to a revision of the fluorescence coefficient for the NO γ-band airglow. Incorporating the quenching reaction of N(²D) with O in the model calculation results in a reduction of the NO density at heights as low as 100 km. These two effects are combined to lead to an evaluation that the N(²D) quantum yield for various possible reactions is as large as 0.9. A smaller rate coefficient for the quenching reaction than that measured in the laboratory, i.e. 1.0 × 10^?12cm³sec^?1 is favourable for the recent NO observation in the early morning, as well as the observed emission rates of the 5200 A airglow from N(²D) The present model predicts a significant day-to-night variation of N and NO densities at heights above 100 km. Below 100 km, the NO density is fairly stable because of its long chemical time constant. Since the rate coefficient for the conversion of N(⁴S) to NO is highly temperature dependent, the relative population of N(⁴S) and NO is very sensitive to the thermospheric temperature variation. Large variations of both N(⁴S) and NO densities due to the temperature change could occur especially at night. The model is in good agreement with the NO observations so far available in low and middle latitudes, as well as the N observation by the use of a rocket in the twilight.     

Atmospheric scattering effects on ground-based measurements of thermospheric winds

Inherent in observations of thermospheric winds from the ground with the Fabry-Perot interferometer is the assumption that the measured Doppler shift is a property of the source medium viewed by the instrumental line of sight. However, ground based airglow observations in regions of weak airglow emission near large intensity gradients may be contaminated by scattered light. Light from areas where the emission is strong can be scattered by the lower atmosphere into the field of view of the observations. Thermospheric winds deduced from the observed Doppler shifts will then show apparent convergence or divergence with respect to the site of observation. Examples of this effect are found in observations by the Michigan Airglow Observatory station located near the auroral zone at Calgary, Alberta. Simulation calculations based upon an experimental model for a significant scattering atmosphere also showed results with either convergence or divergence in the apparent neutral wind field observed by the station.     

Observation of vertical E-region neutral winds in two intense auroral arcs

Vertical winds have been observed by optical Doppler measurements of the 557.7 nm emission in the aurora, using a Fabry-Perot spectrometer. Both upward and downward winds were observed, of 15 m s^?1 magnitude. The upward winds were associated with westward overhead currents, and with low altitude aurora (～ 110 km) as determined by the auroral temperature, while a high altitude aurora (～ 135 km) and eastward currents were associated with the downward wind. The Lorentz force of these currents has the wrong direction to act as a direct forcing mechanism. It is concluded that Joule heating is directly responsible for the upward winds, while the divergence of horizontal winds is responsible for the downward winds.     

Electrical coupling of the E- and F-regions and its effect on F-region drifts and winds

Electric currents, generated by thermospheric winds, flow along the geomagnetic field lines linking the E-and F-regions. Their effects on the electric field distribution are investigated by solving the electrical and dynamical equations. The input data include appropriate models of the F-region tidal winds, the thermospheric pressure distribution and the E-and F-layer concentrations. At the magnetic equator, the calculated neutral air wind at 240 km height has a prevailling eastward component of 55 m sec^-1 and the west-east and vertical ion drifts agree in their general form with incoherent scatter data from Jicamarca     

The influence of thermospheric winds on exospheric hydrogen on Venus

Monte Carlo models of the distribution of atomic hydrogen in the exosphere of Venus were computed which simulate the effects of thermospheric winds and the production of a “hot” hydrogen component by charge exchange of H⁺ and H and O in the exosphere, as well as classic exospheric processes. A thermosphere wind system that is approximated by a retrograde rotating component with equatorial speed of 100 m/sec superimposed on a diurnal solar tide with cross-terminator day-to-night winds of 200 m/sec is shown to be compatible with the thermospheric hydrogen distribution deduced from Pioneer Venus orbiter measurements.     

Horizontal and vertical winds and temperatures in the equatorial thermosphere: Measurements from Natal,Brazil during August–September 1982

Fabry-Perot interferometer measurements of Doppler shifts and widths of the 630.0 nm nightglow line have been used to determine the neutral winds and temperatures in the equatorial thermosphere over Natal, Brazil during August–September 1982. During this period, in the early night (2130 U.T.) the average value of the horizontal wind vector was 95 m s^?1 at 100° azimuth, and the temperature varied from a low of 950 K during geomagnetically quiet conditions to a high of ~ 1400 K during a storm (6 September). The meridional winds were small, ?, 50 m s^?1, and the eastward zonal winds reached a maximum value 1–3 h after sunset, in qualitative agreement with TGCM predictions. On 26 August, an observed persistent convergence in the horizontal meridional flow was accompanied by a downward vertical velocity and an increase in the thermospheric temperature measured overhead. Oscillations with periods of 40–45 min in both the zonal and vertical wind velocities were observed during the geomagnetic storm of 6 September, suggesting gravity wave modulation of the equatorial thermospheric flow.