Chemistry and evolution of Titan's atmosphere

The chemistry and evolution of Titan's atmosphere is reviewed in the light of the scientific findings from the Voyager mission. It is argued that the present N₂ atmosphere may be Titan's initial atmosphere rather than photochemically derived from an original NH₃ atmosphere. The escape rate of hydrogen from Titan is controlled by photochemical production from hydrocarbons. CH₄ is irreversibly converted to less hydrogen rich hydrocarbons, which over geologic time accumulate on the surface to a layer thickness of ～0.5 km. Magnetospheric electrons interacting with Titan's exosphere may dissociate enough N₂ into hot, escaping N atoms to remove ～0.2 of Titan's present atmosphere over geologic time. The energy dissipation of magnetospheric electrons exceeds solar e.u.v. energy deposition in Titan's atmosphere by an order of magnitude and is the principal driver of nitrogen photochemistry. The environmental conditions in Titan's upper atmosphere are favorable to building up complex molecules, particularly in the north polar cap region.     

Radiative equilibrium model of Titan's atmosphere

A simple global radiative equilibrium model is developed for Titan. It is restricted to the two-stream approximation, is vertically homogeneous in its scattering properties, and is spectrally divided into one thermal and two solar channels. A partially absorbing “violet” channel is responsible for heating in the stratosphere, while a conservatively scattering “red” channel permits heating at the surface. The optical thickness of the atmosphere in the red is 1 < τ₁^r < 3. Between 13 and 33% of the total incident solar radiation is absorbed at the planetary surface. The ratio of violet to thermal infrared absorption cross sections is between 30 and 60 in the stratosphere, leading to the large temperature inversion observed there. The observed and theoretically computed tropopause temperatures are 72 and 69°K, respectively, while their corresponding thermal optical depths are, respectively, ～0.1 and ～0.07. The spectrally integrated mass absorption coefficient at thermal wavelengths is approximately constant throughout the stratosphere and roughly linear with pressure in the troposphere. This in turn implies the presence of a uniformly mixed aerosol in the stratosphere, and suggests pressure-induced absorption by gaseous N₂CH₄H₂ in the troposphere. In addition there appear to be two regions of enhanced opacity near 30 and 500 mbar which may be due to C₂H₂C₂H₆C₃H₈ and CH₄ condensation clouds, respectively.     

Behavior of Titan's atmosphere during a total eclipse

The edge-on presentation of Saturn's rings and satellites system has provided a rare opportunity to observe total eclipses of Titan. During its emersion from the Saturnian shadow (1980, June 28), Titan has been observed simultaneously in the visible and the infrared ranges (6000–9000 Å, 11.8 μm and 20 μm). No change has been recorded in these three spectral ranges. Our observations tend to support the thick-atmosphere model, which has been shown to be valid by Voyager a few months later.     

Internal gravity waves in Titan's atmosphere observed by Voyager radio occultation

This paper discusses observational considerations in a search for extraterrestrial intelligence (SETI) program to detect extraterrestrial messenger probes in the solar system. Observable artifacts will most likely be found in a search space consisting of geocentric, selenocentric, Earth-Moon libration, and Earth-Moon halo orbits, which may be searched to a limiting artifact size of 1–10 m (p_v = 0.1) using existing or foreseeable instrumentation.     

Stability of Pluto's atmosphere

L. Trafton (1980, Icarus44, 53–61) has pointed out that a substantial methane atmosphere, observed on Pluto by U. Fink, B.A. Smith, D.C. Benner, J.R. Johnson, and H.J. Reitsema (1980, Icarus44, 62–71), appears to be unstable against blowoff. The difficulty is shown to disappear if the actual heat balance and thermal structure are considered, instead of the classic assumption that the upper atmosphere is isothermal. An energy-limited flux (referred to the surface area) of 3.9 × 10¹⁰ cm^?2 sec^?1 is found. The loss of methane ice over the age of the solar system is an acceptable 3 km.     

Pioneer 10 and 11 observations and the dynamics of Jupiter's atmosphere

Three new results of the Pioneer 10 and 11 mission are discussed. The first is that effective temperature is the same at the poles and equator in spite of the large differences in solar energy deposition. This is consistent with theories of convection which suggest that an extremely small equator-to-pole temperature difference at the level of infrared emission could suppress the internal heat flux at the equator relative to the pole by an amount sufficient to balance the difference in solar energy deposition.The second result is that the effective temperature of belts is 3 to 4K greater than that of zones, which is almost exactly accounted for by the lower albedo of belts. This result cannot be interpreted uniquely, but is consistent with a model in which the internal heat flux is the same under belts and zones, and the horizontal atmospheric heat flux is zero.The third observation provides evidence of instability along the south edges of zones in the northern hemisphere. These are the latitudes of minimum prograde velocity, where instability is most likely to occur in a barotropic fluid, as pointed out by Ingersoll and Cuzzi (1969). A more realistic baroclinic stability analysis suggests instability at these same latitudes.     

The composition of liquid methane-nitrogen aerosols in Titan’s lower atmosphere from Monte Carlo simulations

Molecular level Monte Carlo simulations have been performed with various model potentials for the CH₄-N₂ vapor-liquid equilibrium at conditions prevalent in the atmosphere of Saturn’s moon Titan. With a single potential parameter adjustment to reproduce the vapor-liquid equilibrium at a higher temperature, Monte Carlo simulations are in excellent agreement with available laboratory measurements. The results demonstrate the ability of simple pair potential models to describe phase equilibria with the requisite accuracy for atmospheric modeling, while keeping the number of adjustable parameters at a minimum. This allows for stable extrapolation beyond the range of available laboratory measurements into the supercooled region of the phase diagram, so that Monte Carlo simulations can serve as a reference to validate phenomenological models commonly used in atmospheric modeling. This is most important when the relevant region of the phase diagram lies outside the range of laboratory measurements as in the case of Titan. The present Monte Carlo simulations confirm the validity of phenomenological thermodynamic equations of state specifically designed for application to Titan. The validity extends well into the supercooled region of the phase diagram. The possible range of saturation levels of Titan’s troposphere above altitudes of 7 km is found to be completely determined by the remaining uncertainty of the most recent revision of the Cassini-Huygens data, yielding a saturation of 100 ± 6% with respect to CH₄-N₂ condensation up to an altitude of about 20 km.     

Remarks concerning the diffusion of ion clouds in the earth's upper atmosphere

The diffusive motion of initially ellipsoidal plasma irregularities or ion clouds in the Earth's upper atmosphere is studied theoretically using a model similar to that described by Pickering (1972) for an initially spherical cloud. The work presented here concerns irregularities with major to minor axis ratio between 10:1 and 200:1 at each of the altitudes 97.5 km, 102 km and 114 km (where the ionization could be produced by meteors) and between approximately 200:1 and 1000:1 for altitudes 210 km and 300 km. In particular the effect of the space-charge electric field on the nature of the diffusion process is discussed. The possible effects of ionospheric electric fields and possible relevance to artificial Ba⁺ clouds released in the upper atmosphere are discussed in the second section.     

Jupiter's outer atmosphere

The current state of the theory of Jupiter's outer atmosphere is briefly reviewed. The similarities and dissimilarities between the terrestrial and Jovian upper atmospheres are discussed, including the interaction of the solar wind with the planetary magnetic fields. Estimates of Jovian parameters are given, including magnetosphere and auroral zone sizes, ionospheric conductivity, energy inputs, and solar wind parameters at Jupiter. The influence of the large centrifugal force on the cold plasma distribution is considered. The Jovian Van Alien belt is attributed to solar wind particles diffused in towards the planet by dynamo electric fields from ionospheric neutral winds and consequences of this theory are given.     

The rotation period of Neptune's upper atmosphere

Significant variations in the near-infrared brightness of Neptune during July and August 1980 were observed. These observations show a well-defined, large-amplitude variation in Neptune's J-K color, with a period of 17.73 ± 0.1 hr and are interpreted as diurnal variations resulting from the 17.73-hr rotation period of the upper atmosphere of Neptune in the presence of inhomogeneous weather. These results qualitatively corroborate those of D. P. Cruikshank (1978, Astrophys. J.220, L57-L59) in an earlier study using similar techniques. In addition, variations were observed in the 5-μm spectral region which are in phase with the variations seen at shorter wavelengths. A new 5-μm measurement of Uranus is also reported.     

The dual origin of the terrestrial atmosphere

The origin of the terrestrial atmosphere is one of the most puzzling enigmas in the planetary sciences. It is suggested here that two sources contributed to its formation, fractionated nebular gases and accreted cometary volatiles. During terrestrial growth, a transient gas envelope was fractionated from nebular composition. This transient atmosphere was mixed with cometary material. The fractionation stage resulted in a high Xe/Kr ratio, with xenon being more isotopically fractionated than krypton. Comets delivered volatiles having low Xe/Kr ratios and solar isotopic compositions. The resulting atmosphere had a near-solar Xe/Kr ratio, almost unfractionated krypton delivered by comets, and fractionated xenon inherited from the fractionation episode. The dual origin therefore provides an elegant solution to the long-standing “missing xenon” paradox. It is demonstrated that such a model could explain the isotopic and elemental abundances of Ne, Ar, Kr, and Xe in the terrestrial atmosphere.     

Constraints on bulk composition,seasonal variation,and global dynamics of Pluto's atmosphere

Cosmic abundance, vapor pressure, and molecular weight considerations restrict the likely gas candidates for Pluto's atmosphere to Ne, N₂, CO, O₂, and Ar, in addition to the already detected CH₄. The vapor pressures and cosmic abundances of these gases indicate that all except Ne should be saturated in Pluto's atmosphere. The vapor pressure of Ne is so high that the existence of solid or liquid Ne on Pluto's surface is very unlikely; cosmic abundance arguments imply that Ne cannot attain saturation in Pluto's atmosphere. At both perihelion, N₂ should dominate the saturated gases. CO₂ should have the next highest mixing ratio, followed by O₂ and Ar. CH₄ should have the smallest mixing ratio. Because vapor pressures of these gases vary with temperature at diverse rates, the bulk and constituent mixing ratios of Pluto's atmosphere should vary with season. Between perihelion and aphelion, the column abundance of CH₄ may change by a factor of 260 while that of N₂ changes by only a factor of 52. The potential seasonal variation of Pluto's atmosphere was investigated by considering the behavior of these gases when individually mixed with CH₄. The effects of diurnal and latitudinal variation of insolation and eclipses on the atmosphere also were investigated. Seasonal effects are shown to dominate. It was shown that the atmospheric bulk may not be a minimum near aphelion but rather at intermediate distances from the Sun during summer/winter inadequate ice deposits may allow the atmosphere to collapse by freezing out over winter latitudes. If the atmosphere does not collapse, its weight is sufficient to keep it distributed uniformly around Pluto's surface. In this case, the atmosphere tends to regulate the surface temperature to a seasonally dependent value which is uniform over the globe.Finally, the likely global circulation regimes for each model atmosphere as a function of temperature were investigated and it was concluded that if CH₄, O₂, or CO dominates the atmosphere, Pluto will exhibit cyclic variations between an axially symmetric circulation system at perihelion and a baroclinic wave regime at aphelion. However, if N₂ dominates, as is likely, the wave regime should hold continuously. If the atmosphere collapses to a thin halo during summer/winter seasons, only a weak, symmetric circulation should occur.     

The structure of Titan’s atmosphere from Cassini radio occultations

We present results from the two radio occultations of the Cassini spacecraft by Titan in 2006, which probed mid-southern latitudes. Three of the ingress and egress soundings occurred within a narrow latitude range, 31-34°S near the surface, and the fourth at 52.8°S. Temperature-altitude profiles for all four occultation soundings are presented, and compared with the results of the Voyager 1 radio occultation (Lindal, G.F., Wood, G.E., Hotz, H.B., Sweetnam, D.N., Eshleman, V.R., Tyler, G.L. [1983]. Icarus 53, 348-363), the HASI instrument on the Huygens descent probe (Fulchignoni, M. et al. [2005]. Nature 438, 785-791), and Cassini CIRS results (Flasar, F.M. et al. [2005]. Science 308, 975-978; Achterberg, R.K., Conrath, B.J., Gierasch, P.J., Flasar, F.M., Nixon, C.A. [2008b]. Icarus 194, 263-277). Sources of error in the retrieved temperature-altitude profiles are also discussed, and a major contribution is from spacecraft velocity errors in the reconstructed ephemeris. These can be reduced by using CIRS data at 300 km to make along-track adjustments of the spacecraft timing. The occultation soundings indicate that the temperatures just above the surface at 31-34°S are about 93 K, while that at 53°S is about 1 K colder. At the tropopause, the temperatures at the lower latitudes are all about 70 K, while the 53°S profile is again 1 K colder. The temperature lapse rate in the lowest 2 km for the two ingress (dawn) profiles at 31 and 33°S lie along a dry adiabat except within ∼200 m of the surface, where a small stable inversion occurs. This could be explained by turbulent mixing with low viscosity near the surface. The egress profile near 34°S shows a more complex structure in the lowest 2 km, while the egress profile at 53°S is more stable.     

The microwave spectrum of ammonia in Jupiter's atmosphere

High-frequency-resolution observations of the microwave inversions lines of ammonia in Jupiter have been compared with the models of the temperature inversion in the stratosphere of the planet to deduce that much of the ammonia must be frozen out in the cloud layer, leaving a smaller mixing ration above.     

The photochemical products of benzene in Titan’s upper atmosphere

The Cassini spacecraft detected benzene high in Titan’s atmosphere as well as the presence of large mass positive and negative ions. Previous work has suggested that these large mass ions could be composed of fused-ring polycyclic aromatic hydrocarbon compounds. These fused-ring PAHs, such as naphthalene and anthracene, are usually the result of high temperature processes that may not occur in Titan’s thin, cold, upper thermosphere. Here we suggest that a different class of aromatic compounds, polyphenyls, may be a better explanation of the data. Polyphenyls can grow to be large polymeric structures and could condense to form the aerosols seen in Titan’s cloud and hazes. They have similar properties to fused-ring PAHs (for example, electron affinity, ionization potential) and could be the negative ion species seen in the CAPS instrument data from the Cassini spacecraft.     

NN detectability in Pluto’s atmosphere

Based on the vapor pressure behavior of Pluto’s surface ices, Pluto’s atmosphere is expected to be predominantly composed of N₂ gas. Measurement of the N₂ isotopologue ¹⁵N/¹⁴N ratio within Pluto’s atmosphere would provide important clues to the evolution of Pluto’s atmosphere from the time of formation to its present state. The most straightforward way of determining the N₂ isotopologue ¹⁵N/¹⁴N ratio in Pluto’s atmosphere is via spectroscopic observation of the ¹⁴N¹⁵N gas species. Recent calculations of the 80–100 nm absorption behavior of the ¹⁴N₂ and ¹⁴N¹⁵N isotopologues by Heays et al. (Heays, A.N. et al. [2011]. J. Chem. Phys. 135, 244301), Lewis et al. (Lewis, B.R., Heays, A.N., Gibson, S.T., Lefebvre-Brion, H., Lefebvre, R. [2008]. J. Chem. Phys. 129, 164306); Lewis et al. (Lewis, B.R., Gibson, S.T., Zhang, W., Lefebvre-Brion, H., Robbe, J.-M. [2005]. J. Chem. Phys. 122, 144302), and Haverd et al. (Haverd, V.E., Lewis, B.R., Gibson, S.T., Stark, G. [2005]. J. Chem. Phys. 123, 214304) show that the peak magnitudes of the ¹⁴N₂ and ¹⁴N¹⁵N absorption bandhead cross-sections are similar, but the locations of the bandhead peaks are offset in wavelength by ∼0.05–0.1 nm. These offsets make the segregation of the ¹⁴N₂ and ¹⁴N¹⁵N absorption signatures possible. We use the most recent N₂ isotopologue absorption cross-section calculations and the atmospheric density profiles resulting from photochemical models developed by Krasnopolsky and Cruickshank (Krasnopolsky, V.A., Cruickshank, D.P. [1999]. J. Geophys. Res. 104, 21979–21996) to predict the level of solar light that will be transmitted through Pluto’s atmosphere as a function of altitude during a Pluto solar occultation. We characterize the detectability of the isotopic absorption signature per altitude assuming ¹⁴N¹⁵N concentrations ranging from 0.1% to 2% of the ¹⁴N₂ density and instrumental spectral resolutions ranging from 0.01 to 0.3 nm. Our simulations indicate that optical depth of unity is attained in the key ¹⁴N¹⁵N absorption bands located between 85 and 90 nm at altitudes ∼1100–1600 km above Pluto’s surface. Additionally, an ¹⁴N¹⁵N isotope absorption depth ∼4–15% is predicted for observations obtained at these altitudes at a spectral resolution of ∼0.2–0.3 nm, if the N₂ isotopologue ¹⁵N/¹⁴N percent ratio is comparable to the 0.37–0.6% ratio observed at Earth, Titan and Mars. If we presume that the predicted absorption depth must be at least 25% greater than the expected observational uncertainty, then it follows that a statistically significant detection of these signatures and constraint of the N₂ isotopologue ¹⁴N/¹⁵N ratio within Pluto’s atmosphere will be possible if the attainable observational signal-to noise (S/N) ratio is ?9. The New Horizons (NH) Mission will be able to obtain high S/N, 0.27–0.35 nm full-width half-max 80–100 nm spectral observations of Pluto using the Alice spectrograph. Based on the NH/Alice specifications we have simulated 0.3 nm spectral resolution solar occultation spectra for the 1100–1600 km altitude range, assuming 30 s integration times. These simulations indicate that NH/Alice will obtain spectral observations within this altitude range with a S/N ratio ∼25–50, and should be able to reliably detect the ¹⁴N¹⁵N gas absorption signature between 85 and 90 nm if the ¹⁴N¹⁵N concentration is ∼0.3% or greater. This, additionally, implies that the non-detection of the ¹⁴N¹⁵N species in the 1100–1600 km range by NH/Alice may be used to reliably establish an upper limit to the N₂ isotopologue ¹⁵N/¹⁴N ratio within Pluto’s atmosphere. Similar results may be derived from 0.2 to 0.3 nm spectral resolution observations of any other N₂-rich Solar System or exoplanet atmosphere, provided the observations are attained with similar S/N levels.     

The tilt effect for Saturn's rings

Multiple-scattering computations are carried out to explain the variation of the observed brightness of the A and B rings of Saturn with declination of the Earth and Sun. These computations are performed by a doubling scheme for a homogeneous plane-parallel scattering medium. We test a range of choices for the phase function, albedo for single scattering, and optical depth of both the rings. Isotropic scattering and several other simple phase functions are ruled out, and we find that the phase function must be moderately peaked in both the forward and backward directions. The tilt effect can be explained by multiple scattering in a homogeneous layer, but, for ring B, this requires a single-scattering albedo in excess of 0.8. The brightest part of ring B must have an optical depth greater than 0.9. We find that the tilt effect for ring A can be reproduced by particles having the same properties as those in ring B with the optical depth for the A ring in the range 0.4 to 0.6.     

Physical processes in Jupiter's ring: Clues to its origin by Jove!

The particles making up the Jovian ring may be debris which has been excavated by micrometeoroids from the surfaces of many unseen (R ? 1 km) parent bodies (or “mooms” as we will occasionally call them) residing in the ring. A distribution of particle sizes exists: large objects are sources for the small visible ring particles and also account for the absorption of charged particles noted by Pioneer; the small grains are generated by micrometeoroid impacts, by jostling collisions among different-sized particles, and by self-fracturing due to electrostatic stresses. The latter are most effective in removing surface asperities to thereby produce smooth and crudely equidimensional grains. The presence of intermediate-sized (radius of several to several hundred microns) objects is also expected; these particles will have a total area comparable to the area of the visible ring particles. The nominal size (?2 μm) of the visible particles derived from their forward-scattering characteristics is caused, at least in part, by a selection effect but may also reflect a fundamental grain size or the preferential generation of certain sizes along with the destruction of others. The tiny ring particles have short lifetimes (?10²?10³ years) limited by erosion due to sputtering and meteoroid impacts. Plasma drag significantly modifies orbits in ～10² years but Poynting-Robertson drag is not effective (T_PR ～ 10⁵ years) in removing debris. The ring width is influenced by the distribution of source satellites, by the initial ejection velocity off them, by electromagnetic scattering, and by solar radiation forces. In the absence of electromagnetic forces, debris will reimpact a mother satellite or collide with another particle in about 10 years. A relative drift between different-sized particles, caused by a lessened effective gravity due to the Lorentz force, will substantially shorten these times to less than a month. The ring thickness is determined by a balance between initial conditions (abetted perhaps by electromagnetic scattering) and collisional damping; existence of the “halo” over the diffuse disk compared to its relative absence over the bright ring indicates the presence of mooms in the bright ring but not in the faint disk. Small satellites (R ? 1 km) will not reaccumulate colliding dust grains whereas satellites having the size of J14 or J16 may be able to do so, depending upon their precise shape, size, density, and location. Visible ring structure could indicate separate source satellites. The particles in the faint inner disk are delivered from the bright ring by orbital evolution principally under plasma drag. The halo is comprised of small particles (～0.1 μm) partially drawn out of the faint disk by interactions with the tilted Jovian magnetic field.     

Saturn's rings: II. Condensations of light and optical thickness of Cassini's division

“Condensations” of light have been observed when Saturn's rings are seen almost edge on, and the Sun and the Earth are on opposite sides of the ring plane. These condensations are associated with ring C and Cassini's division. If the relative brightness between the two condensations and the optical thickness of ring C are known, we can calculate the optical thickness of Cassini's division, τ_CASS. Using Barnard's and Sekiguchi's measurements, we have obtained 0.01 ? τ_CASS ? 0.05. A brightness profile of the condensations which agrees well with visual observations is also presented.We are able to set an upper limit of about 0.01 for the optical thickness of any hypothetical outer ring. This rules out a ring observed by C. Cragg in 1954, but does not eliminate the D′ ring observed by Feibelman in 1967.It is known that the outer edge of ring B is almost at the position of the 1/2 resonance with Mimas. Franklin, Colombo, and Cook explained this fact in 1971, postulating a total mass of ring B of 10^?6M_SATURN. We have derived a formula for the mass of the rings, which is a linear function of the mean particle size. We find that 10^?6M_SATURN implies large particles (～70m). If the particles are small (～10cm), as currently believed, the total mass of ring B is not enough to shift the outer edge. We conclude that the above explanation and current size estimates are inconsistent.     

Photochemical aerosols in Titan's atmosphere

The optical properties of polymers, produced photolytically from ethylene, which was detected in Titan's atmosphere and from acetylene or hydrogen cyanide which may be present there, were studied experimentally. It is shown that an aerosol consisting of polyethylene provides an excellent fit to the variation of Titan's albedo with wavelength, while polymers of acetylene or hydrogen cyanide do not. This fit seems to remove the requirement of nitrogen-bearing polymers, which was proposed earlier to account for Titan's red coloration. Therefore, Titan's coloration does not necessarily imply the presence of nitrogen in its atmosphere. It is also proposed that above the layer of larger aerosol particles, whose scattering determines the phase function, there are smaller particles of the same material, which act as an absorbing haze to darken and slightly redden the underlying aerosol. This high-altitude haze also causes the observed strong limb-darkening.