Atmospheric expansion from Joule heating

An expression for the vertical velocity of the neutral atmosphere in the F-region is derived for Joule heating by the electric field that drives the auroral electrojet. When only vertical expansion is allowed, it is found that the vertical wind must always increase monotonically with altitude. The heating rate is proportional to the F-region ion density, so that appreciable heating, even during high electric fields, requires some production mechanism of ionization such as auroral secondary ionization or solar photoionization, in the lower F-region. Once started at night, when an ionizing source is present in the lower F-region, the expansion of the atmosphere transports ionization upward, thereby increasing the heating rate, and hence the expansion rate, i.e. positive feedback. Electric field strengths and F-region ion densities of 50 mV/m and 2 × 10¹¹e/m³, respectively, will produce vertal neutral wind speeds of several tens of m/sec in the 300–500 km altitude range. During periods of high magnetic activity, i.e. high electric field, Joule heating can produce large increases in the relative N₂ concentration in the upper F-region; computations made with a simple model suggest that tenfold increases can occur at 400 km altitude $\text{1}\text{2?1}\text{hr}$ after the onset of magnetic activity, a result in agreement with satellite observations. When the Joule heating theory is applied to incoherent scatter data taken during one period of high heating, the horizontal electric field in the F-region is found to decrease markedly, possibly approaching zero as the field penetrates a weak, discrete auroral arc; the decrease began 10–20 km from the arc.     

The thermospheric heating efficiency under electron precipitation conditions

A time-dependent model calculation of the interaction between energetic electrons and the thermosphere is performed to evaluate the heating efficiency. All energy channels which contribute to the neutral heating of the atmosphere are considered in detail. Neutral chemical reactions and deactivation of metastables are found to be the major heat sources. It is shown that, below 150km, the heating efficiency is on the order of 60% and nearly independent of altitude. At higher altitude, its value decreases in a manner depending on the atmosphere density and composition. It is concluded that, in the dayside polar cusp, particle heating may be an important source of thermospheric motion.     

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.     

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.     

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.     

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.     

Localized ionization patches in the nighttime ionosphere of Mars and their electrodynamic consequences

Using an electron transport model, we calculate the electron density of the electron impact-produced nighttime ionosphere of Mars and its spatial structure. As input we use Mars Global Surveyor electron measurements, including an interval when accelerated electrons were observed. Our calculations show that regions of enhanced ionization are localized and occur near magnetic cusps. Horizontal gradients in the calculated ionospheric electron density on the night side of Mars can exceed 10⁴ cm⁻³ over a distance of a few tens of km; the largest gradients produced by the model are over 600 cm⁻³ km⁻¹. Such large gradients in the plasma density have several important consequences. These large pressure gradients will lead to localized plasma transport perpendicular to the ambient magnetic field which will generate horizontal currents and electric fields. We calculate the magnitude of these currents to be up to 10 nA/m². Additionally, transport of ionospheric plasma by neutral winds, which vary in strength and direction as a function of local time and season, can generate large (up to 1000 nA/m²) and spatially structured horizontal currents where the ions are collisionally coupled to the neutral atmosphere while electrons are not. These currents may contribute to localized Joule heating. In addition, closure of the horizontal currents and electric fields may require the presence of vertical, field-aligned currents and fields which may play a role in high altitude acceleration processes.     

Estimated energy balance in the jovian upper atmosphere during an auroral heating event

We present an analysis of a series of observations of the auroral/polar regions of Jupiter, carried out between September 8 and 11, 1998, making use of the high-resolution spectrometer, CSHELL, on the NASA InfraRed Telescope Facility (IRTF), Mauna Kea, Hawaii; these observations spanned an “auroral heating event.” This analysis combines the measured line intensities and ion velocities with a one-dimensional model vertical profile of the jovian thermosphere/ionosphere. We compute the model line intensities both assuming local thermodynamic equilibrium (LTE) and, relaxing this condition (non-LTE), through detailed balance calculations, in order to compare with the observations. Taking the model parameters derived, we calculate the changes in heating rate required to account for the modelled temperature profiles that are consistent with the measured line intensities. We compute the electron precipitation rates required to give the modelled ion densities that are consistent with the measured line intensities, and derive the corresponding Pedersen conductivities. We compute the changes in heating due to Joule heating and ion drag derived from the measured ion velocities, and modelled conductivities, making use of ion-neutral coupling coefficients derived from a 3-D global circulation model. Finally, we compute the cooling due to the downward conduction of heat and the radiation-to-space from the molecular ion and hydrocarbons. Comparison of the various heating and cooling terms enables us to investigate the balance of energy inputs into the auroral/polar atmosphere. Increases in Joule heating and ion drag are sufficient to explain the observed heating of the atmosphere; increased particle precipitation makes only a minor heating contribution. But local cooling effects—predominantly radiation-to-space—are shown to be too inefficient to allow the atmosphere to relax back to pre-event thermal conditions. Thus we conclude that this event provides observational, i.e. empirical, evidence that heat must be transported away from the auroral/polar regions by thermally or mechanically driven winds.     

The roles of particle precipitation and Joule heating in the energy balance of the Jovian thermosphere

Order of magnitude calculations have been carried out to compare particle precipitation and Joule heating with solar radiation as sources of energy in the Jovian thermosphere. Calculations based on a detailed atomic cross section approach to energy deposition show that the efficiency of conversion of energetic particle precipitation energy into thermal energy is 0.33, larger than on Earth. This emphasizes the role of particle precipitation heating, which may serve as a source for gravity waves. In contrast to the terrestrial case, Joule heating is found to be of only minor significance in the Jovian atmosphere.     

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.     

High altitude haze: Influence of monomer particles on Titan’s temperature profile

The present study investigates the role of high altitude monomer particles in the energy balance of Titan’s upper atmosphere above an assumed low and high aggregate formation altitude of 385 km and 535 km. A ‘single particle approach’ was applied, where the starting point is the energy balance of an individual aerosol. In our analysis 0.01–0.06 μm radius aerosol particles were chosen for the proposed monomer formation regions. These particles absorb solar radiation, emit in the infrared, and are energetically linked to the surrounding gas by thermal conduction. To compute the monomer particle heating effect, the aerosols are assumed to radiate directly to space. We found that high altitude monomers may affect the profile of Titan’s thermosphere from 2 to 20 K depending on the formation altitude of fluffy non-spherical aggregates, the monomer size and distribution. The actual Titan temperature profile in this altitude range including all heating effects will be measured by the HASI instrument during the descent of the Huygens probe.     

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.     

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.     

On the chemistry of H2O H2 and meteoritic ions in the mesosphere and lower thermosphere

In the Earth's lower thermosphere and mesosphere, water vapor is photodissociated by absorption of Lyman alpha radiation. The hydrogen containing free radicals produced by this process lead to the formation of molecular hydrogen. Therefore, very small water vapor mixing ratios are expected at high altitudes, particularly in summer, when photolysis is especially rapid. We present one and two-dimensional model calculations regarding the distribution of H₂O and H₂ in the upper atmosphere.The ion chemistry of meteor ions in the lower thermosphere is also examined and it is shown that silicon ion densities can be used to infer water vapor concentrations near 100 km. The water vapor mixing ratios obtained are generally well below one part per million and are in good agreement with the model calculations.     

Heating of Jupiter’s thermosphere by the dissipation of upward propagating acoustic waves

Thunderstorms in Jupiter’s atmosphere are likely to be prodigious generators of acoustic waves, as are thunderstorms in Earth’s atmosphere. Accordingly, we have used a numerical model to study the dissipation in Jupiter’s thermosphere of upward propagating acoustic waves. Model simulations are performed for a range of wave periods and horizontal wavelengths believed to characterize these acoustic waves. The possibility that the thermospheric waves observed by the Galileo Probe might be acoustic waves is also investigated. Whereas dissipating gravity waves can cool the upper thermosphere through the effects of sensible heat flux divergence, it is found that acoustic waves mainly heat the Jovian thermosphere through effects of molecular dissipation, sensible heat flux divergence, and Eulerian drift work. Only wave-induced pressure gradient work cools the atmosphere, an effect that operates at all altitudes. The sum of all effects is acoustic wave heating at all heights. Acoustic waves and gravity waves heat and cool the atmosphere in fundamentally different ways. Though the amplitudes and mechanical energy fluxes of acoustic waves are poorly constrained in Jupiter’s atmosphere, the calculations suggest that dissipating acoustic waves can locally heat the thermosphere at a significant rate, tens to a hundred Kelvins per day, and thereby account for the high temperatures of Jupiter’s upper atmosphere. It is unlikely that the waves detected by the Galileo Probe were acoustic waves; if they were, they would have heated Jupiter’s thermosphere at enormous rates.     

Dynamical properties of the Venus mesosphere from the radio-occultation experiment VeRa onboard Venus Express

The dynamics of Venus’ mesosphere (60–100 km altitude) was investigated using data acquired by the radio-occultation experiment VeRa on board Venus Express. VeRa provides vertical profiles of density, temperature and pressure between 40 and 90 km of altitude with a vertical resolution of few hundred meters of both the Northern and Southern hemisphere. Pressure and temperature vertical profiles were used to derive zonal winds by applying an approximation of the Navier–Stokes equation, the cyclostrophic balance, which applies well on slowly rotating planets with fast zonal winds, like Venus and Titan. The main features of the retrieved winds are a midlatitude jet with a maximum speed up to 140 ± 15 m s^?1 which extends between 20°S and 50°S latitude at 70 km altitude and a decrease of wind speed with increasing height above the jet. Cyclostrophic winds show satisfactory agreement with the cloud-tracked winds derived from the Venus Monitoring Camera (VMC/VEx) UV images, although a disagreement is observed at the equator and near the pole due to the breakdown of the cyclostrophic approximation. Knowledge of both temperature and wind fields allowed us to study the stability of the atmosphere with respect to convection and turbulence. The Richardson number Ri was evaluated from zonal field of measured temperatures and thermal winds. The atmosphere is characterised by a low value of Richardson number from ～45 km up to ～60 km altitude at all latitudes that corresponds to the lower and middle cloud layer indicating an almost adiabatic atmosphere. A high value of Richardson number was found in the region of the midlatitude jet indicating a highly stable atmosphere. The necessary condition for barotropic instability was verified: it is satisfied on the poleward side of the midlatitude jet, indicating the possible presence of wave instability.     

Mars surface and atmospheric temperature during the 2001 global dust storm

The dramatic growth and evolution of the 2001 martian global dust storm were captured using the Submillimeter Wave Astronomy Satellite (SWAS). While the lower and middle atmosphere (pressures greater than 50 μbar, up to ∼45 km altitude) showed rapid heating of up to 40 K, the average surface brightness temperature plummeted by ∼20 K at the peak of the storm. The storm appears to have had little impact on the global temperature structure at altitudes above ∼ 10 μbar (∼ 60 km).     

Propagation of tropospheric gravity waves into the upper atmosphere of Mars

We study the propagation of gravity waves in the martian atmosphere using a linearized one-dimensional full-wave model. Calculations are carried out for atmospheric parameters characteristic of Mars Orbiter Laser Altimeter (on Mars Global Surveyor MGS) observations of apparent gravity waves in high latitude clouds and MGS radio occultation measurements of temperature variations with height suggestive of gravity wave activity. Waves that reach the thermosphere produce fluctuations in density comparable in amplitude with the density variations detected in Mars Odyssey aerobraking data. Gravity waves of modest amplitude are found to deposit momentum and generate significant heating and cooling in the martian atmosphere. The largest heating and cooling effects occur in the thermosphere, at altitudes between about 130 and 150 km, with heating occurring at the lower altitudes and cooling taking place above.     

Mass density of the upper atmosphere derived from Starlette’s Precise Orbit Determination with Satellite Laser Ranging

The atmospheric mass density of the upper atmosphere from the spherical Starlette satellite’s Precise Orbit Determination is first derived with Satellite Laser Ranging measurements at 815 to 1115 km during strong solar and geomagnetic activities. Starlette’s orbit is determined using the improved orbit determination techniques combining optimum parameters with a precise empirical drag application to a gravity field. MSIS-86 and NRLMSISE-00 atmospheric density models are compared with the Starlette drag-derived atmospheric density of the upper atmosphere. It is found that the variation in the Starlette’s drag coefficient above 800 km corresponds well with the level of geomagnetic activity. This represents that the satellite orbit is mainly perturbed by the Joule heating from geomagnetic activity at the upper atmosphere. This result concludes that MSIS empirical models strongly underestimate the mass density of the upper atmosphere as compared to the Starlette drag-derived atmospheric density during the geomagnetic storms. We suggest that the atmospheric density models should be analyzed with higher altitude acceleration data for a better understanding of long-term solar and geomagnetic effects.     

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.