The martian mesosphere as revealed by CO2 cloud observations and General Circulation Modeling

Different missions have observed mesospheric clouds on Mars in the last years. The presence of these clouds implies, among other conditions, mesospheric temperatures below CO₂ condensation temperature. We use a General Circulation Model to study the mesospheric temperatures and compare the observed distribution of the mesospheric clouds and the predicted climatology of mesospheric temperatures. Although the model does not usually predict temperatures below condensation for daytime conditions, in some regions the predicted temperatures are close enough to condensation that perturbations caused by small scale processes could produce local excursions below condensation. The location and time of the lowest temperatures predicted by the GCM correspond to a first order with the two observed populations of mesospheric clouds: equatorial clouds observed before and after the Northern summer solstice, and mid-latitude clouds observed around the Northern winter solstice. For the equatorial clouds season, the model predicts temperatures close to condensation at the longitude, latitude, altitude and local time where they have been observed. We find that the diurnal migrating thermal tide and non-migrating tides are at the root of the spatial confinement of the equatorial clouds. For the mid-latitude clouds season, the temperatures predicted by the model at the location of the observed clouds is too high. Stereo observations by two different instruments allow for the determination of the zonal speed of these clouds producing a rare dataset of mesospheric winds. We compare the mesospheric zonal winds predicted by the model with these observations, finding a good agreement, although in some cases the observed variability exceeds that predicted by the model.     

Scale-dependent infrared radiative damping rates on Mars and their role in the deposition of gravity-wave momentum flux

Using a Curtis-matrix model of 15 μm CO₂ radiative cooling rates for the martian atmosphere, we have computed vertical scale-dependent IR radiative damping rates from 0 to 200 km altitude over a broad band of vertical wavenumbers ∣m∣ = 2π(1-500 km)⁻¹ for representative meteorological conditions at 40°N and average levels of solar activity and dust loading. In the middle atmosphere, infrared (IR) radiative damping rates increase with decreasing vertical scale and peak in excess of 30 days⁻¹ at ∼50-80 km altitude, before gradually transitioning to scale-independent rates above ∼100 km due to breakdown of local thermodynamic equilibrium. We incorporate these computed IR radiative damping rates into a linear anelastic gravity-wave model to assess the impact of IR radiative damping, relative to wave breaking and molecular viscosity, in the dissipation of gravity-wave momentum flux. The model results indicate that IR radiative damping is the dominant process in dissipating gravity-wave momentum fluxes at ∼0-50 km altitude, and is the dominant process at all altitudes for gravity waves with vertical wavelengths ?10-15 km. Wave breaking becomes dominant at higher altitudes only for “fast” waves of short horizontal and long vertical wavelengths. Molecular viscosity plays a negligible role in overall momentum flux deposition. Our results provide compelling evidence that IR radiative damping is a major, and often dominant physical process controlling the dissipation of gravity-wave momentum fluxes on Mars, and therefore should be incorporated into future parameterizations of gravity-wave drag within Mars GCMs. Lookup tables for doing so, based on the current computations, are provided.     

The density of the upper martian atmosphere measured by Lyman-α absorption with Mars Express SPICAM

Absorption of interplanetary Lyman-α emission by Mars’ nightside lower thermosphere was observed by Mars Express Spectrometer for Investigation of Characteristics of the Atmosphere of Mars (SPICAM), and is analyzed to derive the CO₂ density at 110 km during a martian year. The observed density seasonal variability is consistent with recent observations obtained by stellar occultations, proving that this method, though not as accurate as stellar occultations could be used complementary to them to characterize large variations of thermospheric density on Mars and provide a better spatial coverage by Lyman-α imagery.     

Titan's hydrodynamically escaping atmosphere: Escape rates and the structure of the exobase region

In Strobel [Strobel, D.F., 2008. Icarus, 193, 588-594] a mass loss rate from Titan's upper atmosphere, , was calculated for a single constituent, N₂ atmosphere by hydrodynamic escape as a high density, slow outward expansion driven principally by solar UV heating due to CH₄ absorption. It was estimated, but not proven, that the hydrodynamic mass loss is essentially CH₄ and H₂ escape. Here the individual conservation of momentum equations for the three major components of the upper atmosphere (N₂, CH₄, H₂) are solved in the low Mach number limit and compared with Cassini Ion Neutral Mass Spectrometer (INMS) measurements to demonstrate that light gases (CH₄, H₂) preferentially escape over the heavy gas (N₂). The lightest gas (H₂) escapes with a flux 99% of its limiting flux, whereas CH₄ is restricted to ?75% of its limiting flux because there is insufficient solar power to support escape at the limiting rate. The respective calculated H₂ and CH₄ escape rates are 9.2×¹⁰²⁷ and 1.7×¹⁰²⁷ s−1, for a total of . From the calculated densities, mean free paths of N₂, CH₄, H₂, and macroscopic length scales, an extended region above the classic exobase is inferred where frequent collisions are still occurring and thermal heat conduction can deliver power to lift the escaping gas out of the gravitational potential well. In this region rapid acceleration of CH₄ outflow occurs. With the thermal structure of Titan's thermosphere inferred from INMS data by Müller-Wodarg et al. [Müller-Wodarg, I.C.F., Yelle, R.V., Cui, J., Waite Jr., J.H., 2008. J. Geophys. Res. 113, doi:10.1029/2007JE003033. E10005], in combination with calculated temperature profiles that include sputter induced plasma heating at the exobase, it is concluded that on average that the integrated, globally average, orbit-averaged, plasma heating rate during the Cassini epoch does not exceed ().     

N2 escape rates from Pluto's atmosphere

Hydrodynamic escape of N₂ molecules from Pluto's atmosphere is calculated under the assumption of a high density, slow outflow expansion driven by solar EUV heating by N₂ absorption, near-IR and UV heating by CH₄ absorption, and CO cooling by rotational line emission as a function of solar activity. At 30 AU, the N₂ escape rate varies from in the absence of heating, but driven by an upward thermal heat conduction flux from the stratosphere, for lower boundary temperatures varying from 70-100 K. With solar heating varying from solar minimum to solar maximum conditions and a calculated lower boundary temperature, 88.2 K, the N₂ escape rate range is , respectively. LTE rotational line emission by CO reduces the net solar heat input by at most 35% and plays a minor role in lowering the calculated escape rates, but ensures that the lower boundary temperature can be calculated by radiative equilibrium with near-IR CH₄ heating. While an upward thermal conduction heat flux at the lower boundary plays a fundamental role in the absence of heating, with solar heating it is downward at solar minimum, and is, at most, 13% of the integrated net heating rate over the range of solar activity. For the arrival of the New Horizons spacecraft at Pluto in July 2015, predictions are lower boundary temperature, T₀∼81 K, and N₂ escape rate , and peak thermospheric temperature ∼103 K at 1890 km, based on expected solar medium conditions.     

The Venus nightglow: Ground-based observations and chemical mechanisms

In 1999, observations of the Venus nightglow with the Keck I telescope showed that the 5577 Å oxygen green line was a significant feature, comparable in intensity to the terrestrial green line. Subsequent measurements have been carried out at the Apache Point Observatory (APO) and again at Keck I, confirming the presence of the line with substantially varying intensity. The Herzberg II emission intensity, from the O₂(c-X) transition, was found to have an intensity near 3 kR in one APO run, comparable to the value found on all previous measurements. Thus, of the three oxygen features seen at Venus—the green line, the Herzberg II emission system, and the 1.27-μ 0-0 band of the IR atmospheric system—the first is quite variable, the second is relatively constant, while the third also shows large variations. The reaction between O₂(, v=0) and CO is considered as a possible mechanism to explain green line production and its variability, as well as the variability of the 1.27-μ emission and the stability of the CO₂ atmosphere. This reaction may catalyze CO₂ recombination some five orders of magnitude faster than the slow three-body O + CO reaction.     

Monte Carlo computations of the escape of atomic nitrogen from Mars

Monte Carlo simulations were carried out to compute the escape flux of atomic nitrogen for the low and high solar activity martian thermospheres. The total escape of atomic nitrogen at low and high solar activities was found to be 3.03×¹⁰⁵ and , respectively. The escape flux of atomic nitrogen at low and high solar activities from photodissociation of N₂ was found to be 2.75×¹⁰⁵ and , respectively. The remainder of the contribution is from dissociative recombination, which is only important at high solar activity were it comprises about 25% of the total escape. The relative contributions to the total N escape flux from thermal motion of the background atmosphere, winds and co-rotation, and photoionization and subsequent solar wind pickup are also considered here. We find that the total predicted escape fluxes are observed to increase by 20 and 25% at low and high solar activities owing to thermal motion of the background atmosphere. At low and high solar activities, we find that the co-rotation and wind velocities combined translate to a maximum transferable energy of ∼0.0103 and 0.0181 eV, respectively, and that the total escape flux contribution from winds and co-rotation is negligible. Photoionization was found to be a minor process only impacting those source atoms produced with energies close to the escape energy, between 1.5 and 2 eV. The contributions to the total escape fluxes at low and high solar activities from photoionization and subsequent solar wind pickup are found to be about 8 and 13%, respectively.     

Production of hot nitrogen atoms in the martian thermosphere

The energy and altitude distributions for nitrogen atoms produced from photodissociation of N₂ and dissociative recombination of N⁺₂ have been computed for the low and high solar activity martian thermospheres. We find that the fraction of nitrogen atoms with E>Eesc at the exobase is ∼2.5% for photodissociation as compared to ∼22.6% for dissociative recombination. However, the peak magnitudes of the production rate profiles for dissociative recombination are a factor of ∼1000 and ∼200 times smaller than those of photodissociation at low and high solar activities, respectively. Thus, our findings suggest that photodissociation of N₂ is the dominant escape mechanism for atomic nitrogen for both the low and high solar activity martian thermospheres. The general consensus in the field prior to these research results was that dissociative recombination of N⁺₂ was the main escape mechanism. We find that it is the dissociation of N₂ along repulsive states that results in the production of very energetic nitrogen atoms.     

Modeling the atmospheric limb emission of CO2 at 4.3 μm in the terrestrial planets

The MIPAS instrument on board Envisat, in Earth orbit, the PFS and OMEGA instruments on Mars Express, and VIRTIS on board Venus Express are currently providing a dataset of limb measurements of the CO₂ atmospheric fluorescence emission at 4.3‐μm from the upper atmosphere of the three planets. These measurements represent an excellent dataset to perform comparative studies between the terrestrial planets’ upper atmospheres, and also to test our theoretical understanding of these emissions. In order to exploit these datasets, we apply a set of non-local thermodynamic equilibrium (non-LTE) models developed at the IAA/CSIC, in Granada, Spain, to a selection of data. In general, the models can explain the main spectral features of the measurements, and also the altitude and solar zenith angle variations. However, the simulations for Mars and Venus give an incorrect ratio of the emissions at two wavelengths, 4.4 and . In order to explain this deficiency, a revision of the most uncertain non-LTE energy transfer parameters has been performed. The quenching rate of ν₃ quanta of high-energy CO₂ states by CO₂ itself could reduce the model-data discrepancy if increased by a factor 2-4, still within its current uncertainty range. This factor, however, is subject to the uncertainty in the thermal structure. A number of simulations with the non-LTE models were also used to study and compare the role of radiative transfer in this spectral region in the three terrestrial planets. Sensitivity studies of density and temperature are also presented, and they permit an analysis of how the differences between the planets and between the three instruments affect their sounding capabilities.