Venus mesosphere and thermosphere temperature structure: II. Day-night variations

A two-dimensional nonlinear hydrodynamic model has been developed for studying the global scale winds, temperature, and compositional structure of the mesosphere and thermosphere of Venus. The model is driven by absorption of solar radiation. Ultraviolet radiation produces both heating and photodissociation. Infrared solar heating and thermal cooling are also included with an accurate NLTE treatment. The most crucial uncertainty in determining the solar drive is the efficiency by which $\text{λ < 1080}\text{A}\text{?}$ solar radiation is converted to heat. This question was analyzed in Part I, where it was concluded that essentially all hot atom and O(¹D) energy may be transferred to vibrational-rotational energy of CO₂ molecules. If this is so, the minimum possible euv heating occurs and is determined by the quenching of the resulting excess rotational energy. The hydrodynamic model is integrated with this minimum heating and neglecting any small-scale vertical eddy mixing. The results are compared with predictions of another model with the same physics except that it assumes that 30% of $\text{λ < 1080}\text{A}\text{?}$ radiation goes into heat and that the heating from longer-wavelength radiation includes the O(¹D) energy. For the low-efficiency model, exospheric temperatures are ?300°K on the dayside and drop to < 180°K at the antisolar point. For the higher-efficiency model, the day-to-night temperature variation is from ?600°K to ?250°K. Both versions of the model predict a wind of several hundred meters per second blowing across the terminator and abruptly weakening to small values on the nightside with the mass flow consequently going into a strong tongue of downward motion on the nightside of the terminator. The presence of this circulation could be tested observationally by seeing if its signature can be found in temperature measurements. Both versions of the model indicate that a self-consistent large-scale circulation would maintain oxygen concentrations with ?5% mixing ratios near the dayside F-1 ionospheric peak but ?40% at the antisolar point at the same pressure level.     

On the heat conduction in a high-temperature plasma in solar flares

We have developed three types of mathematical models to describe the mechanisms of plasma heating in the corona by intense heat fluxes from a super-hot (T _e ? 10⁸ K) reconnecting current layer in connection with the problem of energy transport in solar flares. We show that the heat fluxes calculated within the framework of self-similar solutions using Fourier’s classical law exceed considerably the real energy fluxes known from present-day multi-wavelength observations of flares. This is because the conditions for the applicability of ordinary heat conduction due to Coulomb collisions of thermal plasma electrons are violated. Introducing anomalous heat conduction due to the interaction of thermal runaway electrons with ion-acoustic turbulence does not give a simple solution of the problem, because it produces unstable temperature profiles. Themodels incorporating the effect of collisional heat flux relaxation describe better the heat transport in flares than Fourier’s law and anomalous heat conduction.     

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.     

On the structure of Mars' atmosphere at 120–220 km

Altitude dependences of [CO₂] and [CO₂⁺] are deduced from Mariner 6 and 7 CO₂⁺ airglow measurements. CO₂ densities are also obtained from n_e radio occultation measurements. Both [CO₂] profiles are similar and correspond to the model atmosphere of Barth et al. (1972) at 120 km, but at higher altitudes they diverge and at 200–220 km the obtained [CO₂] values are three times less the model. Both the airglow and radio occultation observations show that a correction factor of 2.5 should be included into the values for solar ionization flux given by Hinteregger (1970). The ratio of [CO₂⁺]/n_e is 0.15–0.2 and, hence, [O]/[CO₂] is ～3% at 135 km. An atmospheric and ionospheric model is developed for 120–220 km. The calculated temperature profile is characterized by a value of T ≈ 370°K at h ? 220 km, a steep gradient (～2°/km) at 200-160 km, a bend in the profile at 160 km, a small gradient (～0.7°/km) below and a value of T ≈ 250°K at 120 km. The upper point agrees well with the results of the Lyman-α measurements; the steep gradient may be explained by molecular viscosity dissipation of gravity and acoustical waves (the corresponding energy flux is 4 × 10^?2 erg cm^?2sec^?1 at 180 km). The bend at 160 km may be caused by a sharp decrease of the eddy diffusion coefficient and defines K ≈ 2 × 10⁸cm²sec^?1; and the low gradient gives an estimate of the efficiency of the atmosphere heating by the solar radiation as ? ≈ 0.1.     

Venus' superrotation,mixing length theory and eddy diffusion: A parametric study

The Hadley mechanism is adopted to describe the axisymmetric four day superrotation in the Venus atmosphere, with solar driven meridional winds redistributing energy and momentum, and eddy diffusion describing the actions of three dimensional transient eddies. We address the question how the eddy diffusion coefficients are related to the properties of the circulation. For the atmosphere of a slowly rotating planet such as Venus, we show that a form of the non-linear closure is suggested by the mixing length hypothesis, which constrains the magnitude of the eddy diffusion coefficients. Combining this constraint with the concept of the Rossby radius of deformation yields zonal velocities on the order of 100 m sec^–1. A steady state, non-linear, one-layer spectral model is used for a parametric study to find a relationship between heat source, meridional circulation and eddy diffusion coefficients, which yields the large zonal velocities observed. This analysis leads to the following conclusions: (1) Proportional changes in the heat source and eddy diffusion coefficients do not significantly change the zonal velocities. (2) The meridional velocity is virtually constant for large eddy diffusion coefficients. (3) Below a threshold in the diffusion rate, the meridional velocity decreases, commensurate with the mixing length hypothesis. Eddy heat conduction becomes important and shares with the Hadley cell advection in balancing the solar heating. The zonal velocities then reach large values near 100 m sec^–1. (4) For large eddy diffusion and small heating rates, the zonal velocities decrease with decreasing planetary rotation rates. However, under condition (3), the zonal velocities are independent of the planetary rotation rate. Ramifications are discussed for related parameterizations in GCMs, emphasizing that eddy diffusion coefficients are governed by solar forcing and cannot be chosen independently.     

The structure of cometary atmospheres

The temperature distributions in cometary atmospheres at various heliocentric distances for comets of Bennett and Encke types have been calculated by taking into account heating due to the absorption of solar ultraviolet radiation, cooling by H₂O far infrared emission, and various dynamical processes (expansion, advection, and thermal conduction). The agreement of the results with the observations is in general satisfactory. The conversion of CH₄ and NH₃ to CO and N₂ through thermochemical reaction with H₂O is concluded to be impossible, since the temperature is too low at a heliocentric distance 1.5 AU where CO⁺ ions begin to be observable.     

On the dispersal of artificially-injected gases in the night-time atmosphere

We have studied the extent to which various transport processes affect the dispersal of a gas artificially injected into the night-time atmosphere at F-region altitudes. In addition to diffusion, we have found that nonlinear acceleration, viscous stress, and thermospheric winds affect the dispersal of the injected gas. The magnitude of the effect depends on the atmospheric density, which is a function of solar activity. For an injected H₂ gas, non-linear acceleration and viscous stress rapidly become more important than diffusion above about 300 km for low solar activity (T_∞ = 750K), 340 km for moderate solar activity (T_∞ = 1000K), and 400 km for high solar activity (T_∞ = 1500K). For an injected H₂O gas, the corresponding altitudes are 350, 400, and 470 km for low, moderate and high solar activity, respectively. The effect of nonlinear acceleration and viscous stress is to retard the expansion of the injected gas. Thermospheric winds of 150–400 m s^?1 are important at altitudes near and below the F-region peak electron density. These winds act to transport the injected gas in the wind direction and this affects the shape and temporal development of the subsequent ionospheric hole. Because the H₂O diffusion coefficient is smaller than the H₂ diffusion coefficient, winds are more important for H₂O than for H₂.     

Some problems in Venus' aeronomy

Models are presented for the height distribution of various photochemically active gases in Venus' upper atmosphere. Attention is directed to the chemistry and vertical transport of odd hydrogen (H, OH, HO₂, H₂O₂), odd oxygen (O, O₃), free chlorine (Cl, ClO, ClOO, Cl₂), CO, O₂, H₂ and H₂O. Supply of O₂ may play a limiting role in the formation of a possible H₂SO₄ cloud on Venus. The supply rate is influenced by both chemical and dynamical processes in the stratosphere, and an analysis of recent spectroscopic data for O₂ implies a lower limit to the appropriate eddy coefficient of about 3 × 10⁵ cm²/sec. The abundances of thermospheric O and CO are determined largely by vertical mixing, and an analysis of Mariner 10 measurements of Venus' Lyman α airglow suggests that the eddy coefficient in the lower thermosphere may be as large as 5 × 10⁷cm²sec. The corresponding values for the mixing ratios of O and CO at the ionospheric peak are approximately 1 per cent. The Lyman α data could be reconciled with larger values for thermospheric O, and smaller values for the vertical eddy coefficient, if non-thermal loss processes were to play a dominant role in hydrogen escape, and if the corresponding flux were to exceed 10⁷ atoms/cm²/sec. A sink of this magnitude would imply major depletion of Venus' atmospheric water over geologic time, and would appear to require mixing ratios of H₂O in the lower atmosphere in excess of 4 × 10^?4. The extensive component to the Lyman α emission measured by Mariner 5 may be due to resonance scattering of sunlight by hot atoms formed by charge transfer with O⁺. The H scale height, therefore, may reflect the temperature of positive ions in Venus' topside ionosphere.     

Atmospheric chemistry on Venus, Earth, and Mars: Main features and comparison

This paper deals with two common problems and then considers major aspects of chemistry in the atmospheres of Mars and Venus. (1) The atmospheres of the terrestrial planets have similar origins but different evolutionary pathways because of the different masses and distances to the Sun. Venus lost its water by hydrodynamic escape, Earth lost CO₂ that formed carbonates and is strongly affected by life, Mars lost water in the reaction with iron and then most of the atmosphere by the intense meteorite impacts. (2) In spite of the higher solar radiation on Venus, its thermospheric temperatures are similar to those on Mars because of the greater gravity acceleration and the higher production of O by photolysis of CO₂. O stimulates cooling by the emission at 15 μm in the collisions with CO₂. (3) There is a great progress in the observations of photochemical tracers and minor constituents on Mars in the current decade. This progress is supported by progress in photochemical modeling, especially by photochemical GCMs. Main results in these areas are briefly discussed. The problem of methane presents the controversial aspects of its variations and origin. The reported variations of methane cannot be explained by the existing data on gas-phase and heterogeneous chemistry. The lack of current volcanism, SO₂, and warm spots on Mars favor the biological origin of methane. (4) Venus’ chemistry is rich and covers a wide range of temperatures and pressures and many species. Photochemical models for the middle atmosphere (58-112 km), for the nighttime atmosphere and night airglow at 80-130 km, and the kinetic model for the lower atmosphere are briefly discussed.     

Temperature effect on the photodissociation rates in the atmospheres of Mars and Venus

The available solar flux at a given altitude in the atmospheres of Mars and Venus is attenuated mainly by CO₂ (molecular absorption and Rayleigh scattering) with an extra contribution due to SO₂ on Venus. The dissociation cross section of CO₂ depends on temperature. At temperatures appropriate for these atmospheres (～250°K), the cross sections are about 15% lower than those at room conditions (Y.L. Yung and W.B. De More, 1982, Icarus, 51, 199). It is shown that this temperature effect cannot be neglected in the evaluation of photolysis rates. Calculations of the photodissociation coefficients of CO₂, SO₂, HCl, and H₂O are presented. For example, at the surface of Mars, the coefficient of H₂O is nearly multiplied by a factor of 10!     

The diurnal heat budget of the thermosphere

Detailed numerical calculations of thermospheric heat sources and sinks are presented and their relative importance is discussed in reference to the energy balance phenomena of the neutral atmosphere. It is shown that the thermal energy available from the absorption in the Schumann-Runge continuum leading to photo-dissociation of O₂ is by far the largest energy source in the lower thermosphere. Other sources of varying importance in different altitude ranges are: (1) energy from photoelectrons; (2) energy exchange from thermal plasma; (3) chemical reaction (ion-electron dissociative recombination) energy gain; (4) kinetic and dissipative energy associated with the neutral wind. The energy sinks of importance are (1) thermal conduction at the lower boundary (120km); and (2) radiative cooling of atomic oxygen.It is shown that the combined energy from processes 2–4 constitutes only a small fraction of the total energy available from photoelectrons and is in phase with the latter. These secondary sources (processes 2–4), therefore, do not constitute a significant energy source and their contribution can be simply incorporated into photoelectron energy (process 1) by defining an effective photo-ionization heating efficiency. The heating efficiencies for photo-ionization (including processes 2–4) and photo-dissociation are estimated to be 0.5 and 0.3, respectively.As the important heat input (photo-dissociation) and loss (conduction and radiation) rates are basically governed by the O₂ and O densities, any diurnal or seasonal variation in these constituents at the lower boundary would have profound effects on the thermal structure of the overlying atmosphere. For this and other reasons, it is suggested that a choice of lower boundary much below 120km, e.g. near the mesopause level (90 km), should be more appropriate for general thermospheric studies.     

Modelling the atmospheric CO2 10-μm non-thermal emission in Mars and Venus at high spectral resolution

A study of the CO₂ atmospheric emissions at 10-μm in the upper atmospheres of Mars and Venus is performed in order to explain a number of ground-based measurements of these emissions recently taken at very high spectral resolution in both planets. The measurements are normally used to derive atmospheric temperatures and winds, but uncertainties on the actual emission layers were so far a serious drawback for their correct interpretation. The non-LTE models used for Mars and Venus in the present analysis are entirely similar in order to perform consistent comparisons between the two planets. In particular, the same scheme of CO₂ states and ro-vibrational bands are used, with similar assumptions on collisional routes and rate coef?cients, and also the same radiative transfer approximations. The emissions at 10-μm are produced in both atmospheres by the same excitation mechanism: radiative pumping of the CO₂(00⁰1) vibrational state by direct solar absorption(at 4.3 μm) and indirect absorption (at 2.7 μm, followed by collisional quenching). The computed radiances are specially strong in the upper mesosphere and lower thermosphere of the two planets during maximum solar illumination, producing a population inversion in such conditions with the lower states of the bands, the CO₂ (10⁰0) and CO₂(02⁰0). We obtained that other population inversions are also possible, involving higher energy CO₂ states. The larger solar ?ux available on Venus is found to produce larger vibrational populations and stronger emissions than equivalent atmospheric layers on Mars, in agreement with the observations. A number of perturbation studies were used to determine the exact emission altitudes, or weighting function peaks, for usual nadir sounding. The sensitivity of the emission to non-LTE model uncertainties and to atmospheric variations in temperature and CO₂ density is also presented. The dependence with the solar zenith angle and with the emission angle, as obtained with this model, could also be useful for guiding future observations.     

On the diurnal variations of total mass density,number density and temperature in the upper thermosphere

Theoretical work on the time of occurrence of the diurnal maxima of total mass density (?), number densities (n) and temperature (T) in the upper thermospheric region is discussed. It is suggested that the observed ρ-T phase difference is due to the neutral air wind acting as a heating or cooling agent and to variations in n and T at the lower boundary of the region. The observed n(O)-n(N₂)-T phase differences may be caused by vertical diffusion of atomic oxygen through molecular nitrogen and by variations in n and T at the lower boundary.     

Ionospheric photoelectrons: Comparing Venus, Earth, Mars and Titan

The sunlit portion of planetary ionospheres is sustained by photoionization. This was first confirmed using measurements and modelling at Earth, but recently the Mars Express, Venus Express and Cassini-Huygens missions have revealed the importance of this process at Mars, Venus and Titan, respectively. The primary neutral atmospheric constituents involved (O and CO₂ in the case of Venus and Mars, O and N₂ in the case of Earth and N₂ in the case of Titan) are ionized at each object by EUV solar photons. This process produces photoelectrons with particular spectral characteristics. The electron spectrometers on Venus Express and Mars Express (part of ASPERA-3 and 4, respectively) were designed with excellent energy resolution (ΔE/E=8%) specifically in order to examine the photoelectron spectrum. In addition, the Cassini CAPS electron spectrometer at Saturn also has adequate resolution (ΔE/E=16.7%) to study this population at Titan. At Earth, photoelectrons are well established by in situ measurements, and are even seen in the magnetosphere at up to 7R_E. At Mars, photoelectrons are seen in situ in the ionosphere, but also in the tail at distances out to the Mars Express apoapsis (∼3R_M). At both Venus and Titan, photoelectrons are seen in situ in the ionosphere and in the tail (at up to 1.45R_V and 6.8R_T, respectively). Here, we compare photoelectron measurements at Earth, Venus, Mars and Titan, and in particular show examples of their observation at remote locations from their production point in the dayside ionosphere. This process is found to be common between magnetized and unmagnetized objects. We discuss the role of photoelectrons as tracers of the magnetic connection to the dayside ionosphere, and their possible role in enhancing ion escape.     

Observations of the 10-μm natural laser emission from the mesospheres of Mars and Venus

Nonthermal emission occurs in the cores of the 9.4- and 10.4-μm CO₂ bands on Mars, and has been recently identified as a natural atmospheric laser. This paper presents observations of the total flux and center-to-limb dependence of this emission for Mars and Venus. The emission is believed to be excited by absorption of solar flux in the near-ir CO₂ bands, followed by collisional transfer to the 00°1 state of CO₂. A comparison is made between the observations and a detailed theoretical model based on this mechanism. It is found that the theoretical model successfully reproduces the observed center-to-limb dependence of this emission, to within the limits imposed by the spatial resolution of the observations. A comparison is also made between the observed fluxes and the predictions of the theoretical models. The observed flux from Mars agrees closely with the prediction of the model; the flux observed from Venus is 74% of the flux predicted by the model. This emission is utilized to obtain the kinetic temperatures of the Martian and Venusian mesospheres. For Mars near 70 km altitude, a rotational temperature analysis using five lines gives T = 135 ± 20°K. The frequency width of the emission is also analyzed to derive a temperature of 126 ± 6°K. In the case of the Venusian mesosphere near 109 km, the frequency width of the emission gives T = 204 ± 10°K.     

Solar activity variations of thermospheric temperatures on Mars and a problem of CO in the lower atmosphere

Long-term MGS drag density observations at 390 km reveal variations of the density with season L_S (by a factor of 2) and solar activity index F_10.7 (by a factor of 3 for F_10.7 = 40-100). According to Forbes et al. (Forbes, J.M., Lemoine, F.G., Bruinsma, S.L., Smith, M.D., Zhang, X. [2008]. Geophys. Res. Lett. 35, L01201, doi:10.1029/2007GL031904), the variation with F_10.7 reflects variations of the exospheric temperature from 192 to 284 K. However, the derived temperature range corresponds to variation of the density at 390 km by a factor of 8, far above the observed factor of 3. The recent thermospheric GCMs agree with the derived temperatures but do not prove their adequacy to the MGS densities at 390 km. A model used by Forbes et al. neglects effects of eddy diffusion, chemistry and escape on species densities above 138 km. We have made a 1D-model of neutral and ion composition at 80-400 km that treats selfconsistently chemistry and transport of species with F_10.7, T_∞, and [CO₂]_80 km as input parameters. Applying this model to the MGS densities at 390 km, we find variation of T_∞ from 240 to 280 K for F_10.7 = 40 and 100, respectively. The results are compared with other observations and models. Temperatures from some observations and the latest models disagree with the MGS densities at low and mean solar activity. Linear fits to the exospheric temperatures are T_∞ = 122 + 2.17F_10.7 for the observations, T_∞ = 131 + 1.46F_10.7 for the latest models, and T_∞ = 233 + 0.54F_10.7 for the MGS densities at 390 km. Maybe the observed MGS densities are overestimated near solar minimum when they are low and difficult to measure. Seasonal variations of Mars’ thermosphere corrected for the varying heliocentric distance are mostly due to the density variations in the lower and middle atmosphere and weakly affect thermospheric temperature. Nonthermal escape processes for H, D, H₂, HD, and He are calculated for the solar minimum and maximum conditions.Another problem considered here refers to Mars global photochemistry in the lower and middle atmosphere. The models gave too low abundances of CO, smaller by an order of magnitude than those observed. Our current work shows that modifications in the boundary conditions proposed by Zahnle et al. (Zahnle, K., Haberle, R.M., Catling, D.C., Kasting, J.F. [2008]. J. Geophys. Res. 113, E11004, doi:10.1029/2008JE003160) are reasonable but do not help to solve the problem.     

Variations in ion and neutral composition at Venus: Evidence of solar control of the formation of the predawn bulges in H+ and He1

A comparison of ion and neutral composition measurements at Venus for periods of greatly different solar activity provides qualitative evidence of solar control of the day-to-night transport of light ion and neutral species. Concentrations of H⁺ and He in the predawn bulge near solar maximum in November, 1979, exhibit a depletion signature correlated with a pronounced modulation in the solar F_10.7 and EUV fluxes. This perturbation, not observed in the predawn region during an earlier period of relative quiet solar conditions, is interpreted as resulting from pronounced changes in solar heating and photoionization on the dayside, which in turn modulate the transport of ions and neutrals into the bulge region.     

Buffering of diurnal temperature variations on Venus

Previous calculations of the role of degassing of CO₂ from calcite in buffering the surface temperature of Venus ignored the efficacy with which energy could be conducted into the subsurface. Consuction into the subsurface plays a minor role in the energy balance, however, such that less than 10% of the insolation will conduct into the subsurface and go into degassing CO₂, with the remainder heating up the surface and atmosphere; negligible buffering of the surface temperature will occur.     

Liquid water on Mars: An energy balance climate model for CO2/H2O atmospheres

Geologic evidence of the prior existence of liquid water on Mars suggests surface temperatures T_s were once considerably warmer than at present; and that such a condition may have arisen from a larger atmospheric greenhouse. Here we develop a simple climate model for a CO₂/H₂O Mars atmosphere including water vapor-longwave opacity feedback in the atmosphere and temperature-albedo feedback at surface icecaps, under the assumption that once the Martian surface pressure was p_s ≥ 1 atm CO₂. Longwave flux to space is computed as a function of T_s and p_s using band-absorption models for the effect of the 15-μm fundamental, and the 10- and 15-μm hot bands, of the CO₂ molecule; as well as the pure rotation bands and e continuum of H₂O. The derived global radiative balance predicts a global mean surface temperature of 283°K at 1 atm CO₂. When the emission model is coupled to a latitudinally resolved energy balance climate model, including the effect of poleward heat transfer by atmospheric baroclinic eddies, the solutions vary, depending on p_s. We considered two cases: (1) the present Mars (p_s ? 0.007 atm) with pressure-buffering by solid CO₂ icecaps, and limited poleward heat flux by the atmosphere; and (2) a hypothetical “hot Mars” (p_s ? 1.0 atm), whose much higher CO₂ amount augmented by H₂O evaporative feedback yields a theoretical T_s distribution with latitude admitting liquid water over 95% of the surface, water icecaps at the poles, and a diminished equator-to-pole temperature gradient relative to the present.     

Modeling of the 10-μm natural laser emission from the mesospheres of Mars and Venus

Nonthermal emission which occurs in the cores of the 9.4- and 10.4-μm CO₂ bands on Mars has been recently identified as a natural atmospheric laser. The emission is believed to be excited by absorption of near-ir solar flux, followed by collisional transfer to the 00°1 state of CO₂. A numerical model based on this mechanism is developed which includes the solar pumping contributed by ～2 × 10⁴ near-ir CO₂ lines. The non-Local Thermodynamic Equilibrium (NLTE) radiative transfer problem is solved to obtain the 00°1 vibrational state population. This model successfully reproduces the existing center-to-limb observations, although higher spatial resolution observations are needed for a definitive test. The model also predicts total fluxes which are close to the observed values. The strength of the emission is predicted to be closely related to the instantaneous near-ir solar heating rate. Calculation of the CO₂ level populations in this model supports the identification of this emission as a natural atmospheric laser. Both Mars and Venus are predicted to exhibit level inversions; at large zenith distances several percent of the emergent intensity will be due to stimulated emission. While the stimulated emission present in these lines is not large enough to be of importance to mesospheric radiative equilibrium, it has other interesting consequences. The construction of large-volume radiation-pumped lasers, which utilize CO₂ planetary mesospheres as a gain medium, is theoretically possible.