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.     

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.     

Martian atmospheric Lee waves

The Mariner 9 television pictures of Mars showed areas of extensive mountain lee ware phenomenon in the northern mid-latitudes during winter. In most cases the characteristic wavelength of the lee waves is readily observable and in a few cases the boundaries of the wave patterns, as well as the wavelength, are observed. The cloud patterns resulting from the waves generated by the flow across a mountain or crater are dependent upon the velocity profile of the air stream and the vertical stability of the atmosphere. Using the stability as inferred by the temperature structure obtained from the infrared spectrometer data, a two layer velocity model of the air stream is used in calculations based on the theory of mountain lee waves. The parameters that yield a pattern similar to that in a picture with a well-defined wave configuration are a lower 11 km deep air stream of 40 m/sec and an upper air stream of 85 m/sec. This calculation appears to be an upper limit of the wind speeds, with most of the pictures implying wind speeds of the lower layer to be less than 40 m/sec. These results yield magnitudes generally in agreement with circulation models, in particular, the Leovy and Mintz (1969, J. Atmos. Sci.26, 1167) two layer numerical model. Under not too different conditions they calculate winds of approximately 30 and 70 m/sec in the lower and upper layers, respectively.     

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.     

Midlatitude F-region neutral winds and temperatures during the priority regular world day 14 August 1980

Midlatitude F-region neutral winds and temperatures determined from Fabry-Perot interferometer measurements of the doppler shifts and widths of nightglow 630.0 nm line profiles are presented for the priority regular world day 14 August 1980. They exhibit, in many respects, the observed behavior for other summer, geomagnetically quiet nights at solar maximum. The neutral temperature decreases from 1500°K after sunset (21 h LT) to a minimum of ˜ 1200°K before dawn (05 h LT), except to the north of the observatory. The zonal winds are eastward at sunset at 50 m/sec, decrease to zero at 02 h LT and are westward just before dawn. The meridional winds are zero just after sunset and reach a maximum equatorward value of 50–70 m/sec at local midnight but do not decrease as predicted; instead, they remain at roughly these values towards dawn. The NCAR thermospheric general circulation model (TGCM) is used to predict the global upper atmospheric temperature and circulation patterns for this world day. The model predictions agree with the measured neutral temperatures and exhibit qualitative similarities to the measured neutral winds. It is concluded that inclusion in the model of ion drift at midlatitudes should improve the agreement with observations.     

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.     

An analysis of auroral E region neutral winds based on incoherent scatter radar observations at Chatanika

Auroral E region neutral winds determined from incoherent scatter radar observations at Chatanika, AK, during geomagnetic disturbances (15 May 1974) are compared with detailed theoretical calculations of neutral velocities for these conditions. The theoretical velocities are obtained by numerically solving the ion and neutral momentum equations in the ion drag approximation, including coriolis and viscous forces, using observed electric fields and electron densities. Large vertical gradients are found in the calculated velocities for altitudes below about 130 km. As a consequence of this structure and fluctuations in the electron density profiles, the data analysis procedure of Brekke et al. (1973) for obtaining neutral winds from radar data is found to underestimate the wind speed by up to 40%, but it determines the direction and temporal structure reasonably well. Comparison of observed neutral velocities with calculated values shows that ion drag alone cannot account for the observations. An equation is derived to estimate the pressure gradients required to resolve the discrepancy between calculated and observed neutral winds. Accelerations due to these pressure gradients are of the same order as those due to ion drag, but at least an order of magnitude larger than those due to solar heating. Directions of the horizontal pressure gradients are consistent with expected locations of auroral heating. During geomagnetic disturbances, ion drag and auroral heating both appear to play important roles in the generation and modification of neutral winds.     

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.     

Measurements of high latitude thermospheric winds by rocket and ground-based techniques and their interpretation using a three-dimensional,time-dependent dynamical model

A three-dimensional, time-dependent model of thermospheric dynamics has been used to interpret recent experimental measurements of high altitude winds by rocket-borne and ground-based techniques. The model is global and includes a self-consistent treatment of the non-linear, Coriolis and viscosity terms. The solar u.v. and e.u.v. energy input provides the major energy source for the thermosphere. Solar u.v. and e.u.v. heating appear to be inadequate to explain observed thermospheric temperatures if e.u.v. heating efficiency (ε) lies in the range 0.3 < ε < 0.35. If the recent solar e.u.v. data are correct, then a value of ε between 0.4 and 0.45 would bring fluxes and observed temperatures into agreement. The Heppner (1977) and Volland (1978) models of high-latitude electric field are used to provide sources of both momentum (via ion drag) and energy (via Joule heating). We find that the Heppner Model CO (equivalent to Volland Model 1) is most appropriate for very quiet geomagnetic conditions (K_p ? 2) while Model A (equivalent to Volland Model 2) provides the necessary enhancement at high latitudes for conditions of moderate activity (K_p ～ 4). Even with the addition of a polar electric field, there still appears to be a shortage of high-latitude energy input in that model winds tend to be 10 m s^?1 poleward of observed winds under quiet or average geomagnetic conditions. This extra energy cannot be provided by enhancing the polar electric fields since the extra momentum would cause disagreement with the observed high latitude winds. High latitude particulate sources of relatively low energies, ~100 eV, seem the most likely candidates depositing their energy above about 200km. Relatively modest amounts of energy are then required, < 10¹⁰W global, to bring the model into agreement with both high- and mid-latitude neutral wind results.     

The influence of neutral winds on He+ distributions in the equatorial ionosphere

The effects of F-region neutral winds on the distribution of He⁺ in the equatorial ionosphere have been examined using a theoretical model and an observational data set. It is shown by the model that components of neutral wind in the magnetic meridian up to only 50 m s^? can produce He⁺ gradients in the northern and southern sectors of a flux tube that differ by more than 80%. This is associated with interhemisphere transport velocities of He⁺ as large as 15 m s^?1 at 800 km. A substantial latitude gradient in the He⁺ distribution across the dip equator also results from the redistribution of He⁺ The changes in the He⁺ concentration at the dip equator and the latitude distribution of He⁺ in response to different neutral wind components is determined from the model and used to construct longitude distributions of He⁺ to compare with observations made at equinox. Good agreement between the calculations and observations is obtained both at the geographic and geomagnetic equators using the relationship between neutral winds, interhemispheric transport velocity and He⁺ concentration derived from the model. If these relationships can be extrapolated to accommodate the different conditions expected during solstice, we can also discuss the He⁺ distributions expected during this season.     

Ionospheric winds in the F-region and their effects on the limiting periods of gravity waves

Spectrum analyses of ionospheric electron density and content fluctuations show periods with a lower limit near 5 min. Interpretation of this cut off in terms of gravity waves in a windless atmosphere leads to unacceptably low thermospheric temperatures near 180°K. It is concluded that neutral winds reduce the apparent cut-off period in the ionosphere. The maximum horizontal wind speed obtained from cut-off data is about 100 m/sec.     

The effect of displaced geomagnetic and geographic poles on the thermospheric neutral winds

A three-dimensional, semi-empirical dynamic model of the neutral thermosphere is used to examine the effect of the displaced geomagnetic and geographic poles on the daily variation of neutral gas motion. The global-scale pressure distribution to drive the neutral gas motion is derived from the empirical model of Jacchia (1965). The ionization distribution is obtained from the Pennsylvania State M.K 1 model ionosphere using the first few longitudinal Fourier coefficients of the ionization distribution. The calculations were made at various latitudes at equinox and solstice and for various values of solar activity. The results show that the calculated neutral winds for the case where the geomagnetic and geographic poles are coincident differ at most only a few per cent from the winds calculated assuming the poles displaced. With the poles coincident, longitude and local time are interchangeable, and one dimension in any dynamic model of the thermosphere may be eliminated.     

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 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.     

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.     

On the reality of the Venus winds

An examination of the effect of assumptions in the interpretation of the Venera wind data is made as a rebuttal to the suggestion by A.T. Young that the 140 m/sec Venera 8 horizontal wind at 45 km may be either spurious or anomalous. The Venera measurements of wind speed along with the Mariner measurements of a lower region of strong turbulence are evidence for a wide band of variable high-speed retrograde horizontal winds which girdle Venus at the equator. In the prevalent interpretation of the Mariner 10 uv photographs, the region of the top of the visible cloud is characterized by variable high-speed retrograde horizontal winds which orbit Venus with an average period of 4 Earth days, and by many features indicating vertical convection. This interpretation, together with the possibility of atmospheric corotation due to frictional coupling, suggests that the Venera-Mariner band of winds at 45 km extends well beyond the top of the visible cloud, and that the upper region of strong turbulence detected by the Mariners may result in part from vertical convection currents carried along by high-speed horizontal winds. In an alternate interpretation of the Mariner 10 uv photographs Young suggests that the predominant motions may be traveling wavelike disturbances with a 4-day period rather than bulk motion of the atmosphere. For this case the upper region of strong turbulence is interpreted as due mostly to vertical wind shear resulting from a rapid decrease in wind speed within a relatively short distance above the Venera-Mariner band of high-speed winds.     

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.     

Atmospheric tides on Venus: III. The planetary boundary layer

Diurnal solar heating of Venus' surface produces variable temperatures, winds, and pressure gradients within a shallow layer at the botton of the atmosphere. The corresponding asymmetric mass distribution experiences a tidal torque tending to maintain Venus' slow retrograde rotation. It is shown that including viscosity in the boundary layer does not materially affect the balance of torques. On the other hand, friction between the air and ground can reduce the predicted wind speeds from ∽5 to ∽1 m/sec in the lower atmosphere, more consistent with the observations from Venus landers and descent probes. Implications for aeolian activity on Venus' surface and for future missions are discussed.