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.     

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.     

Global distribution of helium in the upper atmosphere during solar minimum

Helium in the Earth's thermosphere traces the dynamical systems that redistribute energy and mass. Measurements of the global helium distribution in the thermosphere, using Atmosphere Explorer satellite C. (AE-C), show a gradual seasonal change in the number density of helium for all latitudes. The enhancement in helium over the winter pole (the helium bulge) changes in magnitude slowly as seasons progress. The bulge builds and recedes following the progression of winter North to South and back again. This progression of the winter helium enhancement is presented in this paper using latitudinal profiles of helium number density for each month during the year. The absolute magnitude of the winter helium enhancement in the auroral regions is affected by auroral heating at low altitudes. The reduction in the winter helium bulge at low altitudes shown in AE-C data can be traced to this localized heating. The gradual variation in helium concentration measured at many latitudes for all seasons of the year implies that global thermospheric wind systems change gradually with the seasons.     

Eddy mixing and circulation at ionospheric levels

The rates of heat input into the mesosphere and lower thermosphere are calculated and compared with the heat losses. The worldwide average eddy diffusion coefficient required to maintain continuity in the heat budget is calculated and found to vary from about 10⁷ cm²/sec at 120 km down to about 10⁵ cm²/sec at 60 km. From the global asymmetry in heating at the solstice, it is concluded that a systematic pattern of vertical velocities prevails ranging from less than 1 cm/sec in the mesosphere up to 10 cm/sec near 120 km, upward over the summer polar region and downward over the winter polar region. This can be balanced by a wind system towards the winter polar region with velocities near 1 m/sec at 60 km increasing to 30 m/sec at 120 km. Such a wind system provides an explanation for the helium bulge in the upper thermosphere over the winter polar region.     

Geomagnetic storm effects on the thermosphere and the ionosphere revealed by in situ measurements from OGO 6

The temporal response of ion and neutral densities to a geomagnetic storm has been investigated on a global scale with data from consecutive orbits of OGO-6 (>400km) for 4 days covering both magnetically quiet and disturbed conditions. The first response of the neutral atmosphere to the storm takes place in the H and He densities which start to decrease near the time of the storm sudden commencement. The maximum decreases in H and He were more than 40% of the normal density at high latitudes. A subsequent increase in O and N₂ densities occurs about 8 hours later than the change in H and He densities, while the relative O and N₂ density changes indicate a depletion of atomic oxygen in the lower thermosphere by more than a factor of two. The overall features of the change in the neutral atmosphere, especially the patterns of change for individual species, strongly support the physical picture that energy is deposited primarily at high latitudes during the storm, and the thermosphere structure changes through (1) heating of the lower thermosphere and (2) generation of large scale circulation in the atmosphere with upwelling at high latitudes and subsidence at the equator. The storm-time response of H⁺ occurs in two distinct regions separated by the low latitude boundary of the light ion trough. While on the poleward side of the boundary the H⁺ density decreases in a similar manner to the decrease in H density, on the equatorward side of the boundary the H⁺ decrease occurs about half a day later. It is shown that the decrease of H⁺ density is principally caused by the decrease in H density for both regions. The difference in H⁺ response between the two regions is interpreted as the difference in H⁺ dynamics outside and inside the plasmasphere. The O⁺ density shows an increase, the pattern of which is rather similar to that for O. Two possibilities for explaining the observed change in O⁺ density are suggested. One attributes the observed increase in O⁺ density to an increase in the plasma temperature during the storm. The other possibility is that the increase in the production rate of O⁺ due to an increase in O density exceeds the increase in the loss rate of O⁺ due to an increase in N₂ density, especially around the time of sunrise. Hence the change in O⁺ density in the F-region may actually be controlled by the change in O density.     

Comments on stubbe's deviations from the barometric law

An approximation for the height distribution of nitrogen and atomic oxygen in the thermosphere between 120 and 500 km is obtained by the assumption that these constituents are in diffusive equilibrium. The horizontal wind system in the thermosphere may cause deviations from the diffusive equilibrium distribution due to the non-zero divergence of the horizontal flow and the vertical mass motion associated with it. Stubbe has examined this problem and deduced considerable deviations from diffusive equilibrium due to this effect. The problem is re-examined here, and it is shown that there exist arguments that the deviations from diffusive equilibrium due to horizontal wind fields are not as large as suggested by Stubbe.     

Latitude variations of exospheric hydrogen and the polar wind

Several satellite experiments have measured the solar Lyman-α line, either in scattering from upper atmospheric atomic hydrogen (the Lyman-α airglow) or directly at line center (which determines the hydrogen column density along the line of sight). Recent analyses of data from the above experiments consistently reveal the presence of an atomic hydrogen depletion at high latitudes. In situ determinations of hydrogen at lower altitude show no evidence of such behaviour. This has led us to postulate two mechanisms which may be more effective in reducing the high-latitude density at the high altitudes of the exospheric measurements (500–2000 km). The first is the polar wind loss of protons, which depletes atomic hydrogen through a charge exchange reaction. The second is a high-latitude magnetospheric heating of protons, followed by charge exchange. Opposing the above loss mechanisms are the influences of ballistic lateral flow and mean meriodional winds. We have shown by means of a three-dimensional exospheric transport model that none of the above mechanisms can reconcile the disparate results in the two altitude regimes, nor can they provide the large outward hydrogen fluxes and the correct seasonal variations observed at high latitudes.     

Supernovae effects on the terrestrial atmosphere

The first effects of a nearby (～ 10 parsec) supernova on the Earth's atmosphere will be caused by ultraviolet radiation dissociating molecular oxygen. The event will be of about one month's duration. Several months later nuclear gamma radiation may arrive, causing a decrease in atmospheric ozone. Cosmic radiation from the supernova remnant will not intercept the Earth for at least 1000 years at which time ozone will be seriously depleted.Supernova ultraviolet radiation increases column ozone and atomic oxygen. Atmospheric thermal structure is modified with a large temperature increase in the mesosphere and lower thermosphere and a decrease at higher altitudes caused by enhanced heat loss due to atomic oxygen radiation and conduction.     

Thermal equilibrium and hydrostatic equilibrium for excited states in planetary atmospheres

Thermal equilibrium and hydrostatic equilibrium are mutually exclusive for any particular quantum state of an atmospheric constituent in a non-isothermal atmosphere. As a result, there is a flux of rotationally, vibrationally, and electronically excited atoms and molecules down the temperature gradient, balanced by an up-gradient transport of ground-state atoms and molecules, resulting in a net transport of excitation energy, but with no net mass transport. The energy flux is first formulated as a molecular process and applied to vibrationally excited molecular nitrogen and rotationally excited atomic oxygen in the Earth's lower thermosphere, then reformulated as a bulk process and applied to the Venusian atmosphere, where it is shown that the CO₂ vibrational flux is a significant contribution to the total eddy energy flux in the 0–60 km region.     

Impact of CIR Storms on Thermosphere Density Variability during the Solar Minimum of 2008

The solar minimum of 2008 was exceptionally quiet, with sunspot numbers at their lowest in 75 years. During this unique solar-minimum epoch, however, solar-wind high-speed streams emanating from near-equatorial coronal holes occurred frequently and were the primary contributor to the recurrent geomagnetic activity at Earth. These conditions enabled the isolation of forcing by geomagnetic activity on the preconditioned solar minimum state of the upper atmosphere caused by Corotating Interaction Regions (CIRs). Thermosphere density observations around 400 km from the CHAMP satellite are used to study the thermosphere density response to solar-wind high-speed streams/CIRs. Superposed epoch results show that the thermosphere density responds to high-speed streams globally, and the density at 400 km changes by 75% on average. The relative changes of neutral density are comparable at different latitudes, although its variability is largest at high latitudes. In addition, the response of thermosphere density to high-speed streams is larger at night than in daytime, indicating the preconditioning effect of the thermosphere response to storms. Finally, the thermosphere density variations at the periods of 9 and 13.5 days associated with CIRs are linked to the spatial distribution of low?–?middle latitude coronal holes on the basis of the EUVI observations from STEREO.     

Stratospheric superrotation in the TitanWRF model

TitanWRF general circulation model simulations performed without sub-grid-scale horizontal diffusion of momentum produce roughly the observed amount of superrotation in Titan’s stratosphere. We compare these results to Cassini-Huygens measurements of Titan’s winds and temperatures, and predict temperature and winds at future seasons. We use angular momentum and transformed Eulerian mean diagnostics to show that equatorial superrotation is generated during episodic angular momentum ‘transfer events’ during model spin-up, and maintained by similar (yet shorter) events once the model has reached steady state. We then use wave and barotropic instability analysis to suggest that these transfer events are produced by barotropic waves, generated at low latitudes then propagating poleward through a critical layer, thus accelerating low latitudes while decelerating the mid-to-high latitude jet in the late fall through early spring hemisphere. Finally, we identify the dominant waves responsible for the transfers of angular momentum close to northern winter solstice during spin-up and at steady state. Problems with our simulations include peak latitudinal temperature gradients and zonal winds occurring ∼60 km lower than observed by Cassini CIRS, and no reduction in zonal wind speed around 80 km, as was observed by Huygens. While the latter may have been due to transient effects (e.g. gravity waves), the former suggests that our low (∼420 km) model top is adversely affecting the circulation near the jet peak, and/or that we require active haze transport in order to correctly model heating rates and thus the circulation. Future work will include running the model with a higher top, and including advection of a haze particle size distribution.     

Influence of gravity wave activity on lower thermospheric photochemistry and composition

The wind and temperature oscillations of internal gravity waves can cause horizontal variations of a factor of two in minor gas number densities in the lower thermosphere over length scales of several hundred kilometers. The variations are due both to vertical transport of constituents whose lifetimes are long compared to the wave period and to chemical activity driven by temperature dependent reaction rate coefficients. The nightglow emission of the hydroxyl radical provides a remote sensor of wave activity between 80 and 90 km. Theoretical calculations show that the horizontal variations in the atomic hydrogen distribution are the largest single contributor to wave structure in the nightglow followed by the effects of temperature fluctuations on the rate coefficient of the reaction H + O₃ → O₂ + OHv′ > 0).     

Venus zonal wind above the cloud layer

The altitude variation of the zonal wind velocity in the Venus atmosphere above the cloud layer is deduced from the structure of the wavenumber 2 solar tide. Results show that the amplitude of the zonal wind increases with respect to altitude near the equator, but decreases for latitudes greater than 30°. Thus, the zonal wind becomes concentrated at lower latitudes by 100 km altitude.     

The atmospheric circulation and dust activity in different orbital epochs on Mars

A general circulation model is used to evaluate changes to the circulation and dust transport in the martian atmosphere for a range of past orbital conditions. A dust transport scheme, including parameterized dust lifting, is incorporated within the model to enable passive or radiatively active dust transport. The focus is on changes which relate to surface features, as these may potentially be verified by observations. Obliquity variations have the largest impact, as they affect the latitudinal distribution of solar heating. At low obliquities permanent CO₂ ice caps form at both poles, lowering mean surface pressures. At higher obliquities, solar insolation peaks at higher summer latitudes near solstice, producing a stronger, broader meridional circulation and a larger seasonal CO₂ ice cap in winter. Near-surface winds associated with the main meridional circulation intensify and extend polewards, with changes in cap edge position also affecting the flow. Hence the model predicts significant changes in surface wind directions as well as magnitudes. Dust lifting by wind stress increases with obliquity as the meridional circulation and associated near-surface winds strengthen. If active dust transport is used, then lifting rates increase further in response to the larger atmospheric dust opacities (hence circulation) produced. Dust lifting by dust devils increases more gradually with obliquity, having a weaker link to the meridional circulation. The primary effect of varying eccentricity is to change the impact of varying the areocentric longitude of perihelion, l, which determines when the solar forcing is strongest. The atmospheric circulation is stronger when l aligns with solstice rather than equinox, and there is also a bias from the martian topography, resulting in the strongest circulations when perihelion is at northern winter solstice. Net dust accumulation depends on both lifting and deposition. Dust which has been well mixed within the atmosphere is deposited preferentially over high topography. For wind stress lifting, the combination produces peak net removal within western boundary currents and southern midlatitude bands, and net accumulation concentrated in Arabia and Tharsis. In active dust transport experiments, dust is also scoured from northern midlatitudes during winter, further confining peak accumulation to equatorial regions. As obliquity increases, polar accumulation rates increase for wind stress lifting and are largest for high eccentricities when perihelion occurs during northern winter. For dust devil lifting, polar accumulation rates increase (though less rapidly) with obliquity above o=25°, but increase with decreasing obliquity below this, thus polar dust accumulation at low obliquities may be increasingly due to dust lifted by dust devils. For all cases discussed, the pole receiving most dust shifts from north to south as obliquity is increased.     

Modelling the Global Solar Corona II: Coronal Evolution and Filament Chirality Comparison

This paper considers the hemispheric pattern of solar filaments using newly developed simulations of the real photospheric and 3D coronal magnetic fields over a six-month period, on a global scale. The magnetic field direction in the simulation is compared directly with the chirality of observed filaments, at their observed locations. In our model the coronal field evolves through a continuous sequence of nonlinear force-free equilibria, in response to the changing photospheric boundary conditions and the emergence of new magnetic flux. In total 119 magnetic bipoles with properties matching observed active regions are inserted. These bipoles emerge twisted and inject magnetic helicity into the solar atmosphere. When we choose the sign of this active-region helicity to match that observed in each hemisphere, the model produces the correct chirality for up to 96% of filaments, including exceptions to the hemispheric pattern. If the emerging bipoles have zero helicity, or helicity of the opposite sign, then this percentage is much reduced. In addition, the simulation produces a higher proportion of filaments with the correct chirality after longer times. This indicates that a key element in the evolution of the coronal field is its long-term memory, and the build-up and transport of helicity from low to high latitudes over many months. It highlights the importance of continuous evolution of the coronal field, rather than independent extrapolations at different times. This has significant consequences for future modelling such as that related to the origin and development of coronal mass ejections.     

Titan global climate model: A new 3-dimensional version of the IPSL Titan GCM

We have developed a new 3-dimensional climate model for Titan’s atmosphere, using the physics of the IPSL Titan 2-dimensional climate model with the current version of the LMDZ General Circulation Model dynamical core. Microphysics and photochemistry are still computed as zonal averages. This GCM covers altitudes from surface to 500 km altitude, with barotropic waves now being resolved and the diurnal cycle included. The boundary layer scheme has been changed, yielding a strong improvement in the tropospheric zonal wind profile modeled at Huygens descent position and season. The potential temperature profile is fairly consistent with Huygens observations in the lowest 10 km. The latitudinal profile of the near-surface temperature is close to observed values. The minimum of zonal wind observed by the Huygens probe just above the tropopause is also present in these simulations, and its origin is discussed by comparing solar heating and dynamical transport of energy. The stratospheric temperature and wind fields are consistent with our previous works. Compared to observations, the zonal wind peak is too weak (around 120 m/s) and too low (around 200 km). The temperature structures appear to be compressed in altitude, and depart strongly from observations in the upper stratosphere. These discrepancies are correlated, and most probably related to the altitude of the haze production. The model produces a detached haze layer located more than 150 km lower than observed by the Cassini instruments. This low production altitude is due to the current position of the GCM upper boundary. However, the temporal behaviour of the detached haze layer in the model may explain the seasonal differences observed between Cassini and Voyager 1. The waves present in the GCM are analyzed, together with their respective roles in the angular momentum budget. Though the role of the mean meridional circulation in momentum transport is similar to previous work, and the transport by barotropic waves is clearly seen in the stratosphere, a significant part of the transport at high latitudes is done all year long through low-frequency tropospheric waves that may be baroclinic waves.     

Asymmetry of solar active prominences separately at low and high latitudes from 1957 to 1998

The paper presents the results of a study of the asymmetry of the solar active prominences (SAP) at low (≤40°) and high (≥50°) latitudes, respectively, from 1957 through 1998 (solar cycles 19–22). A quantitative analysis of the hemispheric distribution of the SAP is given. We found that the annual hemispheric asymmetry indeed exists at low latitudes, but strangely, a similar asymmetry does not seem to occur for SAPs at high latitudes. We found that the north–south (N–S) asymmetry of the solar active prominences at high latitudes is always north dominated during solar cycles 19–22 while the N–S asymmetry of the SAPs at low latitudes is shifted to a dominance in the southern hemisphere for solar cycle 21 and remains south dominated even in cycle 22. Thus, the hemispheric asymmetry of the solar active prominences at high latitudes in a cycle appears to have little connection with the asymmetry of the solar activity at low latitudes.     

The density and temperature structure near the exobase of Saturn from Cassini UVIS solar occultations

We analyzed 15 solar occultations observed by the Cassini UVIS instrument to constrain the density and temperature structure near the exobase of Saturn. We retrieved the density of H₂ and thus the temperature at altitudes higher than 1900 km above the 1 bar level by analyzing the ionization continuum of H₂ at wavelengths shorter than 804 Å. We find that the exospheric temperature ranges from 370 K to 540 K, with a typical uncertainty of less than 20 K. According to our data the temperature increases with latitude from the equator to the poles by 100–150 K. At similar latitudes, the temperature varies by 20–50 K at different times with no evidence for any systematic diurnal trend so far. Based on our data, the exobase of Saturn is 2700–3000 km above the 1 bar level and the thermal escape parameter near the exobase ranges from 260 to 340, implying that thermal escape from Saturn is firmly in the Jeans regime. The mixing ratio of H₂ is close to unity at all altitudes below the exobase. We find that the pressure levels in the thermosphere deviate significantly from a simple spheroid predicted by potential theory. This is consistent with significant meridional temperature variations in the lower thermosphere. A global analysis of the temperature structure at different depths in the atmosphere is required to constrain both the shape and the deposition and redistribution of energy in the upper atmosphere further.     

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.