Reevaluation of thermosphere heating by auroral electrons

This paper presents a new calculation of neutral gas heating by precipitating auroral electrons. It is found that the heating rate of the neutral gas is significantly lower than previous determinations below 200 km altitude. The neutral gas heating arises from the many exothermic chemical reactions that take place from the ions and excited species created by the energetic electrons. The calculations show that less than half the energy initially deposited ends up heating the neutral gases. The rest is radiated or lost in the dissociation of O₂ because the O atoms do not recombine in the thermosphere. This paper also presents a new way of calculating the heating rate per ionization that can be used for efficient determination of the overall neutral gas heating for global thermosphere models. The heating rates are relatively insensitive to the neutral atmosphere when plotted against pressure rather than altitude coordinates. At high altitudes, the heating rates are sensitive to the thermal electron density and long-lived species. The calculations were performed with the Field Line Interhemispheric Plasma (FLIP) model using a 2-stream auroral electron precipitation model. The heating rate calculations in this paper differ from previous heating rate calculations in the treatment of backscattered electrons to produce better agreement with observed flux spectra. This paper shows that more realistic model auroral electron spectra can be obtained by reflecting the up going flux back to the ionosphere at the upper boundary of the model. In this case, the neutral gas heating rates are 20%–25% higher than when the backscattered flux escapes from the ionosphere.     

Large-scale circulation of atomic oxygen in the upper mesosphere and lower thermosphere

The variability of the atomic oxygen green line airglow at 557.7 nm, originating from the O(¹S) level, has a long history of observation. Only recently, global observations of the oxygen airglow, interpreted with the help of global circulation models have provided some understanding. Satellite observations of winds and temperatures clearly demonstrate the dynamical influence of tides on the daily variations. Both annual and semi-annual components of the variation have been identified as occurring over the course of the year. From the large-scale circulation of the atmosphere one expects downwelling in the winter at high latitudes, causing enhanced atomic oxygen and enhanced airglow while the opposite is expected in summer. In the present investigation a search is made for this large-scale signature using data from the WIND Imaging Interferometer on the Upper Atmosphere Research Satellite, model results from the Thermosphere–Ionosphere–Mesosphere-Electrodynamics-Global Circulation Model, and earlier published results. A well-defined semi-annual variation of emission rate is found in the tropics, apparently the result of the semi-annual variation of the diurnal tide. Annual and semi-annual patterns are found at mid-latitudes, in satellite and model data, while ground-based observations detect only an annual variation, with a maximum in the autumn. At still higher latitudes the fall peak persists, but with a deep depletion of atomic oxygen in the springtime; this appears to be the signature of the large-scale circulation.     

Properties of the mesosphere and lower thermosphere

The variability and systematic variations of the properties of the upper mesosphere and lower thermosphere are probably the least well known aspects of the terrestrial atmosphere. Satellite measurements of this region are very limited and rocket and remote sounding techniques do not provide comprehensive coverage. Progress is being made in theoretical studies of this region, primarily with regard to tidal effects, and some progress is being made in analyzing the relatively sparse experimental data that are available. Turbulence dynamics of the region has been studied by analyzing structure measurements at Kwajalein, wind data from Natal and systematic variations of the turbopause altitude determined from measurements of the diffusive separation of argon. One question that is being raised at this time, and it is appropriate at a time near solar maximum, is the extent of solar activity control of the properties of this region of the atmosphere. The occurrence rates and magnitudes of the turbulent diffusivity in the 70 to 90 km altitude region appear to correlate with solar activity with a time lag, as do also the incidence of aurora and the atomic oxygen green line intensity. Solar cycle dependence has been identified in mean zonal wind speeds in the 65 to 110 km altitude region above Saskatoon and in lower thermosphere temperatures measured at Heiss Island and at St. Santin. Millstone Hill data show that the mean meridional wind changes during a solar cycle. Solar cycle variations have also been detected in the stratosphere and troposphere.     

A comparison of semi-empirical models in the lower thermosphere

Semi-empirical models are derived predominantly from satellite-borne observations. The nature of these observations restricts the applicability of the models mainly to the atmospheric regions sampled, i.e. the upper thermosphere. Current models are only capable of reproducing a zero-order approximation of the structure of the lower thermosphere. Based on selected examples, the progress in atmospheric research since CIRA-72 as well as the continuing deficiencies are demonstrated.     

Planetary and gravity waves in the mesosphere and lower thermosphere

Planetary and gravity waves contribute significantly to the variability of atmospheric parameters in the middle atmosphere. In the mesosphere and lower thermosphere the wave fluctuations are sufficiently large to often mask the prevailing or mean state of the atmosphere. This review summarizes current knowledge about the motion, temperature and density fields associated with both large and small scale waves and stresses improved understanding that has come from recent ground based, rocket and satellite investigations.     

Semi-empirical model of middle atmosphere wind from the ground to the lower thermosphere

During recent years, special attention has been paid to understanding the background circulation of the middle atmosphere. Particularly in the mesosphere/lower thermosphere (MLT) region, this has involved including data from a range of new radar measurements. It has also involved the comparison of existing empirical middle atmosphere wind models, such as CIRA-86 and HWM-93 to the new data. This has led to the construction of empirical models of MLT winds such as the Global Empirical Wind Model (GEWM). Further investigations are aimed at the construction of new empirical and semi-empirical wind models of the entire middle atmosphere including these new experimental results. The results of a new wind climatology (0–100 km) are presented here, based upon the GEWM, a reanalysis of stratospheric data, and a numerical model which is used to fill the gap between data from the stratospheric and MLT regions.     

Excitations of nonmigrating diurnal tides in the mesosphere and lower thermosphere simulated by the Kyushu-GCM

Excitation mechanisms of nonmigrating diurnal tides in the MLT region simulated by the Kyushu-GCM are examined. It is shown that the westward propagating diurnal tide with zonal wavenumber s = 2 is mainly excited by nonlinear interactions between the migrating diurnal tide and the stationary planetary wave with zonal wavenumber s = 1, while the nonlinear excitation of the standing diurnal tide with zonal wavenumber s = 0 is less important than the excitation by tropospheric heating. Nonlinear interactions between the migrating diurnal tide and the stationary planetary wave with zonal wavenumber s = 2 are not dominant to excite the westward propagating diurnal tide with zonal wavenumber s = 3, and it is shown that the excitation by tropospheric heating is comparable to the nonlinear excitation. It is also shown that other nonmigrating diurnal tides are excited by tropospheric heating.     

Detection of long-term trends in the mesosphere/lower thermosphere from ground-based radio propagation measurements

Mesospheric temperature trends can be derived from LF phase-height observations in mid-latitudes supported by ionospheric absorption and ionosonde observations. Analysing the full observation period from 1959 until 2003, a mean yearly temperature trend has been derived with −0.25 K/yr for the height interval from 48 to 82 km. Subdividing the whole observation interval in two parts before and after 1979, the trend is markedly stronger in the second period with −0.38 K/yr compared with −0.20 K/yr in the first part before 1979. These differences can at least partly be explained by a steeper CO₂ increase and ozone decrease in the second interval. The differences in the mesospheric temperature trends are most evidently expressed during winter months and are markedly smaller during summer season. The reason of this seasonal difference is not quite clear; it may be related with detected ozone trends which are clearly stronger during winter months on both hemispheres.     

A global model of circulation and temperature for the upper mesosphere - lower thermosphere region

The zonal mean model of zonal and meridional wind for the Northern and Southern Hemisphere based the analysis of meteor radar wind and partial reflection drift data for the 70–110-km height interval is constructed. The height-latitudinal cross-sections of vertical wind are calculated from data on the latitudinal structure of a meridional wind field using the continuity equation. The temperature field cross-sections from the zonal wind model using the thermal wind equation are derived.     

On the structure and dynamics of the thermosphere

Thermospheric temperature, composition and wind measurements from the Dynamics Explorer satellite (DE-2) are interpreted using a three dimensional, multiconstituent spectral model. The analysis accounts for tides driven by the absorbed solar radiation as well as energy and momentum coupling involving the magnetosphere and lower atmosphere. We discuss phenomena associated with the annual tide, polar circulation, magnetic storms and substorms.     

Properties of the neutral density and composition in the thermosphere

Measurements of the density and composition of the thermosphere between 150 and 500 km, which were obtained by the S3-1 satellite, have been compared with the Jacchia and MSIS models. The measurements of the densities of O, N₂, N and Ar show some differences from the current models which should be considered during the preparation of the next CIRA model. The Ar measurements are particularly useful in examining the response of the neutral atmosphere to geomagnetic heating. These results are useful in establishing the appropriate lower boundary conditions for modeling of the thermosphere.     

Solar EUV energy budget of the thermosphere

Using the accumulation of experimental data and theoretical studies conducted on the terrestrial thermosphere since the mid seventies, we have re-evaluated the channels by which solar UV energy is transferred to the atmosphere. As an outcome of this evaluation we have redetermined the solar EUV heating efficiency for the thermosphere and find this to be considerably different from that established in earlier studies. The heating efficiency has strong altitude, solar cycle and diurnal dependencies. The values of this parameter vary from less than 10% to greater than 100%, with peak midday values of 50–55% In recent papers we have presented the results of this new UV heating efficiency determination using a steady state solution of the ionospheric model. In this paper we present the results obtained solving a time dependent model over a diurnal cycle. The time dependent effects are found to be significant, with certain longlived species acting as temporary reservoirs of latent heat that is released to the neutral atmosphere at later times.     

On the formation of auroral arcs

A new mechanism for auroral arc formation is presented. The characteristic linear shape of auroral arcs is determined by magnetically connected plasma clouds in the distant equatorial magnetosphere. These clouds originate as high speed plasma beams in the magnetotail and in the solar wind. It is found that the free energy for driving an auroral arc can be provided by the difference of pressure between the cloud and the ambient plasma.     

A three-dimensional modelling study of the processes leading to mid latitude nitric oxide increases in the lower thermosphere following periods of high geomagnetic activity

The processes leading to enhancements in mid latitude nitric oxide (NO) densities following geomagnetic storms have been investigated using the University College London (UCL) Coupled Middle Atmosphere and Thermosphere (CMAT) general circulation model. A comparison of calculated storm time and quiet time NO densities at 110 km altitude reveals the presence of aurorally produced NO at both high and mid latitudes for several days after subsidence of activity. At 150 km, the NO enhancements are shorter lived and remain for up to approximately 2 days after the storm. By separating the contribution of chemical production and loss, horizontal and vertical advection, and molecular and eddy diffusion in the calculation of NO densities, we show that at 150 km altitude, horizontal transport must be taken into consideration if post-storm mid latitude enhancements are to be reproduced. Chemical production of NO at high latitudes continues for up to 2 days after subsidence of a storm at altitudes of around 150 km. We show that equatorward winds at this altitude are sufficiently strong to transport the aurorally produced NO to mid latitudes. Vertical diffusion transports NO from altitudes of 150 km and above, to lower altitudes where it is longer lived. At 110 km altitude, chemical, diffusive and advective terms must all be included in the calculation of NO density in order to simulate realistic mid latitude enhancements. We propose that it is the combined effects of increased chemical production, downward diffusion from altitudes of 150 km and above, and transport by winds that lead to increases in mid latitude NO density at altitudes of around 110 km. This is the first detailed study of the causes of post-storm mid latitude NO enhancements to use a three-dimensional general circulation model.     

Empirical models of the thermosphere and requirements for improvements

The different types of variation in the thermosphere are briefly examined and the solar-activity effect is singled out for special attention. To this day, empirical models have made use of the decimetric solar flux F_10.7 as an index of the variable XUV radiation from the sun. To account for the change in the relative intensity of the different types of emissions in the course of the solar cycle, F_10.7 is made to perform double duty: The daily values are used to represent the day-to-day and “27-day” variations, while its averages over several solar rotations are used to represent the variations with the 11-year cycle. The availability of direct solar XUV data should eventually eliminate the need for such a make-shift procedure. Accuracy and continuity requirements of XUV intensity measurements are discussed and a strategy is outlined for sorting out the relevant features from the observational material and putting them to practical use in thermospheric modeling. It is suggested that future models of the diurnal and the geomagnetic variation use as a guide theoretical models which have achieved considerable success in qualitatively representing the observed phenomena.     

An investigation of the solar cycle impact on the lower thermosphere O(S) nightglow emission as observed by WINDII/UARS

The yearly variation of the integrated emission rate of the O(¹S) nightglow in the lower thermosphere is studied and the solar cycle impact is examined from the observations of the Wind Imaging Interferometer (WINDII) operated on the Upper Atmosphere Research Satellite (UARS). More than 300,000 volume emission rate profiles of the O(¹S) nightglow observed by WINDII for 40°S–40°N latitudes during November 1991–August 1997 over half of a solar cycle are utilized. These profiles are vertically integrated for the altitude range of 80–100 km and the equivalent column integrated emission rates are then zonally averaged for bins with 10° latitude and 3 month intervals. It is found that for each latitude the O(¹S) nightglow emission rate appears to increase with increasing solar F10.7 cm flux, following a linear relationship. This characterizes the solar cycle impact on the O(¹S) nightglow, while the solar influence is modulated by a seasonal variation. Based on these variations, an empirical formula is derived for predicting the three-month averages of the O(¹S) nightglow integrated emission rate. The standard error of the estimated values from the formula is smaller than 30 Rayleigh.     

Advances in auroral imaging from space

Space imaging of the aurora provides unique quantitative information on the instantaneous global distribution and pattern of the aurora - a direct result of particle acceleration and precipitation processes occurring throughout the magnetosphere. Techniques for imaging from space have usually involved systematically scanning a detector across the earth using some combination of satellite motion and spin, scanning mirrors, and/or electronic scanning. The recent advent of CCD image sensor arrays opens up possibilities for vastly improved sensitivity and resolution in satellite imagers, although scanning systems retain significant advantages in situations where scattered light is a problem and where narrow band optical interference filters are needed.A new type of imaging instrument being built in Canada for the Swedish Viking satellite will employ a two-dimensional CCD array detector which is fibre-optically coupled to a micro-channel plate intensifier. It will produce images of low light level auroral emissions at the rate of one image per spin of the satellite (20 seconds). The instrument operates in the vacuum ultraviolet and uses a fast inverse Cassegrain optical system. Charge is shifted in the CCD imaging array in synchronism with the moving image so that exposure times of approximately one second will be achieved.     

Models of the geomagnetic effect in the earth's thermosphere

A study was made of the variations in exospheric temperature in the disturbed thermosphere as a function of geomagnetic latitude and local time and of the level of disturbance. Temperatures were derived from the N₂ densities measured by the gas analyzer aboard the polar-orbiting satellite ESRO 4. Results indicated two maxima in the temperature response in high geomagnetic latitudes that are probably associated with particle precipitation. These maxima occur near the poleward limit of the auroral belts, at about 9 and 0 hours of local magnetic time. There is also evidence of the effect of Joule heating due to the auroral electrojets in the morning and evening auroral belts. Enhancement of the temperature response in middle latitudes is observed throughout the night side and in a region centered in the late afternoon.     

Perspectives of auroral X-ray measurements in the 1980s

The measurement of auroral X-rays with balloon-borne instruments is an efficient means to study the behaviour of electrons with energies above about 30 keV in the magnetosphere during disturbed periods. Possibilities will be discussed to continue such measurements in the 1980's. It will be pointed out, what kind of investigations may be performed. Recently developed payloads will be described that can be used as a basis for further technical developments. Satellite projects scheduled for the 1980's will be presented that are suited for coordinated X-ray measurements.