Titan solar occultation observed by Cassini/VIMS: Gas absorption and constraints on aerosol composition

A solar occultation by Titan's atmosphere has been observed through the solar port of the Cassini/VIMS instrument on January 15th, 2006. Transmission spectra acquired during solar egress probe the atmosphere in the altitude range 70 to 900 km at the latitude of 71° S. Several molecular absorption bands of CH₄ and CO are visible in these data. A line-by-line radiative transfer calculation in spherical geometry is used to model three methane bands (1.7, 2.3, 3.3 μm) and the CO 4.7 μm band. Above 200 km, the methane 2.3 μm band is well fit with constant mixing ratio between 1.4 and 1.7%, in agreement with in situ and other Cassini measurements. Under 200 km, there are discrepancies between models and observations that are yet fully understood. Under 480 km, the 3.3 μm CH₄ band is mixed with a large and deep additional absorption. It corresponds to the C-H stretching mode of aliphatic hydrocarbon chains attached to large organic molecules. The CO 4.7 μm band is observed in the lower stratosphere (altitudes below 150 km) and is well fit with a model with constant mixing ratio of 33±10 ppm. The continuum level of the observed transmission spectra provides new constraints on the aerosol content of the atmosphere. A model using fractal aggregates and optical properties of tholins produced by Khare et al. [Khare, B.N., Sagan, C., Arakawa, E.T., Suits, F., Callcott, T.A., Williams, M.W., 1984. Icarus 60, 127-137] is developed. Fractal aggregates with more than 1000 spheres of radius 0.05 μm are needed to fit the data. Clear differences in the chemical composition are revealed between tholins and actual haze particles. Extinction and density profiles are also retrieved using an inversion of the continuum values. An exponential increase of the haze number density is observed under 420 km with a typical scale height of 60 km.     

The mesosphere and lower thermosphere of Titan revealed by Cassini/UVIS stellar occultations

Stellar occultations observed by the Cassini/UVIS instrument provide unique data that probe the mesosphere and lower thermosphere of Titan at altitudes between 400 and 1400 km. This region is a site of complex photochemistry that forms hydrocarbon and nitrile species, and plays a crucial role in the formation of the organic hazes observed in the stratosphere, but has yet to be adequately characterized. We analyzed publicly available data obtained between flybys Tb in December 2004 and T58 in July 2009, with an emphasis on two stable occultations obtained during flybys T41 and T53. We derived detailed density profiles for CH₄, C₂H₂, C₂H₄, C₄H₂, HCN, HC₃N and C₆H₆ between ∼400 and 1200 km and extinction coefficients for aerosols between 400 and 900 km. Our analysis reveals the presence of extinction layers in the occultation data that are associated with large perturbations in the density profiles of the gaseous species and extinction profiles of the aerosols. These relatively stable features vary in appearance with location and change slowly over time. In particular, we identify a sharp extinction layer between 450 and 550 km that coincides with the detached haze layer. In line with recent images obtained by Cassini/ISS, the altitude of this layer changes rapidly around the equinox in 2009. Our results point to unexpected complexity that may have significant consequences for the dynamics and physical processes taking place in the upper atmosphere of Titan.     

An analysis of Pluto occultation light curves using an atmospheric radiative-conductive model

We use a radiative-conductive model to least-squares fit Pluto stellar occultation light curve data. This model predicts atmospheric temperature based on surface temperature, surface pressure, surface radius, and CH₄ and CO mixing ratios, from which model light curves are to be calculated. The model improves upon previous techniques for deriving Pluto’s atmospheric thermal structure from stellar occultation light curves by calculating temperature (as a function of height) caused by heating and cooling by species in Pluto’s atmosphere, instead of a general assumption that temperature follows a power law with height or some other idealized function. We are able to fit for model surface radius, surface pressure, and CH₄ mixing ratio with one of the 2006 datasets and for surface pressure and CH₄ mixing ratio for other datasets from the years 1988, 2002, 2006, and 2008. It was not possible to fit for CO mixing ratio and surface temperature because the light curves are not sensitive to these parameters. We determine that the model surface radius, under the assumption of a stratosphere only (i.e. no troposphere) model in radiative-conductive balance, is . The CH₄ mixing ratio results are more scattered with time and are in the range of 1.8-9.4 × 10⁻³. The surface pressure results show an increasing trend from 1988 to 2002, although it is not as dramatic as the factor of 2 from previous studies.     

An investigation of Pluto’s troposphere using stellar occultation light curves and an atmospheric radiative-conductive-convective model

We use a radiative-conductive-convective model to assess the height of Pluto’s troposphere, as well as surface pressure and surface radius, from stellar occultation data from the years 1988, 2002, and 2006. The height of the troposphere, if it exists, is less than 1 km for all years analyzed. Pluto has at most a planetary boundary layer and not a troposphere. As in previous analyses of Pluto occultation light curves, we find that the surface pressure is increasing with time, assuming that latitude and longitude variations in Pluto’s atmosphere are negligible. The surface pressure is found to be slightly higher ( μbar in 1988,  μbar in 2002, and 18.5 ± 4.7 μbar in 2006) than in our previous analyses with the troposphere excluded. The surface radius is determined to be . Comparison of the minimum reduced chi-squared values between the best-fit radiative-conductive-convective (i.e., troposphere-included) model and best-fit radiative-conductive (i.e., troposphere-excluded) shows that the troposphere-included model is only a slightly better fit to the data for all 3 years. Uncertainties in the small-scale physical processes of Pluto’s lower atmosphere and consequently the functional form of the model troposphere lend more confidence to the troposphere-excluded results.     

New perspectives on Titan's upper atmosphere from a reanalysis of the Voyager 1 UVS solar occultations

We have reanalyzed the Voyager 1 UVS solar occultations by Titan to expand upon previous analyses and to resolve inconsistencies that have been noted in the scientific literature. To do so, we have developed a detailed model of the UVS detector and improved both the data reduction methods and retrieval techniques. In comparison to the values previously determined by Smith et al. (1982, J. Geophys. Res. 87, 1351-1359) we find N₂ densities that are 25-60% higher, CH₄ densities that are smaller by a factor of 3-7, and C₂H₂ densities that are roughly two orders of magnitude smaller. Our values for the thermospheric temperature are 153-158 K, which are approximately 20-40 K colder than previous estimates. We also report the first-ever determination from Voyager UVS data of density profile information for C₂H₄, HCN, and HC₃N. Finally, we present a simple engineering model that is consistent with our new results in the upper atmosphere and merges smoothly with the model of Yelle et al. (1997, in: HUYGENS Science, Payload and Mission, in: ESA SP, vol. 1177, pp. 243-256) in the lower atmosphere. Our results provide improved constraints for photochemical models and offer scientists a better understanding of Titan's upper atmosphere as we head into the Cassini era in the exploration of the saturnian system.     

The 10 October 1999 HIP 9369 occultation by the northern polar region of Jupiter: ingress and egress lightcurves analysis

The occultation of bright star HIP9369 by the northern polar region of Jupiter was observed from four locations in North and South America, providing four data sets for ingress and egress. The inversion of the eight occultation lightcurves provides temperature profiles at different latitudes ranging from 55°N to 73.2°N. We estimate the errors on the profiles due to the uncertainties of the inversion method and compare the value of the temperature at the deepest level probed (∼ 50 μbar) with previous observations. The shape of the temperature gradient profile is found similar to previous investigations of planetary atmospheres with propagating and breaking gravity waves. We analyze the small scale structures in both lightcurves and temperature profiles using the continuous wavelet transform. The calculated power spectra of localized fluctuations in the temperature profiles show slopes close to −3 for all eight profiles. We also isolate and reconstruct the three-dimensional geometry of a single wave mode with vertical and horizontal wavelengths of respectively 3 and 70 km. The identified wave is consistent with the gravity wave regime, with a horizontal phase speed nearly parallel to the planetary meridian. Nevertheless, the dissipation of the corresponding wave in Jupiter’s stratosphere should preclude its detection at the observed levels and an acoustic wave cannot be ruled out.     

Impact of seas/lakes on polar meteorology of Titan: Simulation by a coupled GCM-Sea model

The detection of large hydrocarbon seas/lakes near the poles by the Cassini spacecraft raises the question as to whether and how polar seas affect the meteorology on Titan. The polar meteorology and methane hydrological cycle in the presence of seas are investigated by a three-dimensional atmospheric general circulation model coupled to a one-dimensional sea energy balance model considering the observed sea/lake geography. The sea composition has a large control on the seasonal evolution of seas, temperature and wind system in the polar region, particularly in the north where large seas are located. The surface of ethane-rich seas, which do not evaporate methane, undergo a large seasonal temperature variation and the sea surface is often warmer than the surrounding land surface. Land breeze in summer towards the seas causes a moisture convergence over the seas, which leads to enhanced summer precipitation in the sea area. On the other hand, methane-rich seas evaporate some methane and are therefore colder than the surroundings. This causes a sea breeze across the north pole in summer, which blows away the moisture from the polar region, so precipitation becomes scarce in the north polar region. The breeze can become stronger than the tidal wind. Sea evaporation peaks in winter, when the temperature and average methane mixing ratio in the planetary boundary layer become lowest. The sea level predominantly rises in summer by precipitation and retreats in winter by evaporation. The meteorology in the south polar region is less sensitive to the composition of the lakes because of the paucity and smallness of southern lakes. Lake-effect precipitation can occur either by moisture convergence by the breeze or humidity enhancement over the seas, but is more characteristic of warm seasons than of cold seasons.     

A re-analysis of the 1971 Beta Scorpii occultation by Jupiter: study of temperature fluctuations and detection of wave activity

Data from the 13 May 1971 β Scorpii occultation by the southern polar region of Jupiter (Vapillon et al., 1973, Astron. Astrophys. 29, 135-149) are re-analyzed with current methods. We correct the previous results for an inacurrate background estimation and calculate new temperature profiles, that are now consistent with the results of other observers of this occultation, as well as with the current knowledge of the jovian atmosphere. The characteristics of the profiles of temperature gradient and the spectral behavior of the temperature fluctuations are found to be similar to the results of previous investigations of planetary atmospheres and in agreement with the presence of atmospheric propagating gravity waves in the jovian atmosphere. We use a wavelet analysis of the temperature profiles to identify the dominant modes of wave activity and compare the reconstructed temperature fluctuations to model-generated gravity waves.     

The impact of methane thermodynamics on seasonal convection and circulation in a model Titan atmosphere

We identify mechanisms controlling the distribution of methane convection and large-scale circulation in a simplified, axisymmetric model atmosphere of Titan forced by gray radiation and moist (methane) convection. The large-scale overturning circulation, or Hadley cell, is global in latitudinal extent and provides fundamental control of precipitation and tropospheric winds. The precipitating, large-scale updraft regularly oscillates in latitude with seasons. The distance of greatest poleward excursion of the Hadley cell updraft is set by the mass of the convective layer of the atmosphere; convection efficiently communicates seasonal warming of the surface through the cold and dense lower atmosphere, increasing the heat capacity of the system. The presence of deep, precipitating convection introduces three effects relative to the case with no methane latent heating: (1) convection is narrowed and enhanced in the large-scale updraft of the Hadley cell; (2) the latitudinal amplitude of Hadley cell updraft oscillations is decreased; and (3) a time lag is introduced. These effects are observable in the location and timing of convective methane clouds in Titan’s atmosphere as a function of season. A comparison of simulations over a range of convective regimes with available observations suggest methane thermodynamic-dynamic feedback is important in the Titan climate.     

Meteorological assessment of the surface temperatures on Titan: constraints on the surface type

The latitudinal profile of near-surface air temperature on Titan retrieved by Voyager 1 has been difficult to understand and raised several speculations about possible exotic processes that might be occurring near Titan's surface, while the thermal properties of the surface itself are unknown. This study systematically investigates the seasonal and spatial variation of the surface temperature and air temperature in the lower troposphere by a 3-dimensional general circulation model for different putative surface types (porous icy regolith, rock-ice mixture, hydrocarbon lakes). For any viable surface type the surface temperature is unlikely to be constant through the year and should more or less vary seasonally and even diurnally, most likely by a few K. Recent observations of tropospheric clouds may be evidence of seasonal variation of the surface temperature and the model predicts in the case of solid surface the development of a convective layer with superadiabatic lapse rates near the surface exactly at those latitudes and seasons where clouds have been identified. The latitudinal profile of the surface temperature retrieved from Voyager 1 infrared spectra can be explained without invoking exotic effects, provided the thermal inertia of the surface is relatively small and/or the surface albedo is low. A dominance of water ice (high thermal inertia and high albedo) at the surface is unfavorable to reproduce the observation. The latitudinal gradient of the surface temperature is particularly large at the hydrocarbon lake surface due to low albedo and small surface drag. Local anomalies of the surface albedo or surface thermal inertia are likely to cause substantial inhomogeneities of the surface temperature. Quasi-permanent accumulation of stratospheric haze at both poles would create a perennial equator-to-pole contrast of the surface temperature, but also a substantially lower global-mean surface temperature due to an enhanced anti-greenhouse effect in summer. The air temperature in the lower troposphere exhibits a tiny latitudinal gradient and a pole-to-pole gradient due to the presence of a pole-to-pole Hadley circulation, indicating that the temperature within the planetary boundary layer may exhibit a vertical profile characteristic of season, location and scenario. There may be a shallow near-surface inversion layer in cold seasons and a shallow convective layer in warm seasons.     

Atmospheric acoustics of Titan, Mars, Venus, and Earth

Planetary atmospheres are complex dynamical systems whose structure, composition, and dynamics intimately affect the propagation of sound. Thus, acoustic waves, being coupled directly to the medium, can effectively probe planetary environments. Here we show how the acoustic absorption and speed of sound in the atmospheres of Venus, Mars, Titan, and Earth (as predicted by a recent molecular acoustics model) mirror the different environments. Starting at the surface, where the sound speed ranges from ∼200 m/s for Titan to ∼410 m/s for Venus, the vertical sound speed profiles reveal differences in the atmospheres' thermal layering and composition. The absorption profiles are relatively smooth for Mars, Titan, and Earth while Venus stands out with a noticeable attenuation dip occurring between 40 and 100 km. We also simulate a descent module sampling the sound field produced by a low-frequency “event” near the surface noting the occurrence of acoustic quiet zones.     

New Horizons Alice ultraviolet observations of a stellar occultation by Jupiter’s atmosphere

The Alice ultraviolet spectrograph onboard the New Horizons spacecraft observed two occultations of the bright star χ Ophiucus by Jupiter’s atmosphere on February 22 and 23, 2007 during the approach phase of the Jupiter flyby. The ingress occultation probed the atmosphere at 32°N latitude near the dawn terminator, while egress probed 18°N latitude near the dusk terminator. A detailed analysis of both the ingress and egress occultations, including the effects of molecular hydrogen, methane, acetylene, ethylene, and ethane absorptions in the far ultraviolet (FUV), constrains the eddy diffusion coefficient at the homopause level to be  cm² s⁻¹, consistent with Voyager measurements and other analyses (Festou, M.C., Atreya, S.K., Donahue, T.M., Sandel, B.R., Shemansky, D.E., Broadfoot, A.L. [1981]. J. Geophys. Res. 86, 5717-5725; Vervack Jr., R.J., Sandel, B.R., Gladstone, G.R., McConnell, J.C., Parkinson, C.D. [1995]. Icarus 114, 163-173; Yelle, R.V., Young, L.A., Vervack Jr., R.J., Young, R., Pfister, L., Sandel, B.R. [1996]. J. Geophys. Res. 101 (E1), 2149-2162). However, the actual derived pressure level of the methane homopause for both occultations differs from that derived by [Festou et al., 1981] and [Yelle et al., 1996] from the Voyager ultraviolet occultations, suggesting possible changes in the strength of atmospheric mixing with time. We find that at 32°N latitude, the methane concentration is  cm⁻³ at 70,397 km, the methane concentration is  cm⁻³ at 70,383 km, the acetylene concentration is  cm⁻³ at 70,364 km, and the ethane concentration is  cm⁻³ at 70,360 km. At 18°N latitude, the methane concentration is  cm⁻³ at 71,345 km, the methane concentration is  cm⁻³ at 71,332 km, the acetylene concentration is cm⁻³ at 71,318 km, and the ethane concentration is  cm⁻³ at 71,315 km. We also find that the H₂ occultation light curve is best reproduced if the atmosphere remains cold in the microbar region such that the base of the thermosphere is located at a lower pressure level than that determined by in situ instruments aboard the Galileo probe (Seiff, A., Kirk, D.B., Knight, T.C.D., Young, R.E., Mihalov, J.D., Young, L.A., Milos, F.S., Schubert, G., Blanchard, R.C., Atkinson, D. [1998]. J. Geophys. Res. 103 (E10), 22857-22889) - the Sieff et al. temperature profile leads to too much absorption from H₂ at high altitudes. However, this result is highly model dependent and non-unique. The observations and analysis help constrain photochemical models of Jupiter’s atmosphere.     

Analysis of Cassini/CIRS limb spectra of Titan acquired during the nominal mission II: Aerosol extinction profiles in the 600–1420 cm spectral range

We have analyzed the continuum emission of limb spectra acquired by the Cassini/CIRS infrared spectrometer in order to derive information on haze extinction in the 3–0.02 mbar range (∼150–350 km). We focused on the 600–1420 cm⁻¹ spectral range and studied nine different limb observations acquired during the Cassini nominal mission at 55°S, 20°S, 5°N, 30°N, 40°N, 45°N, 55°N, 70°N and 80°N. By means of an inversion algorithm solving the radiative transfer equation, we derived the vertical profiles of haze extinction coefficients from 17 spectral ranges of 20-cm⁻¹ wide at each of the nine latitudes. At a given latitude, all extinction vertical profiles retrieved from various spectral intervals between 600 and 1120 cm⁻¹ display similar vertical slopes implying similar spectral characteristics of the material at all altitudes. We calculated a mean vertical extinction profile for each latitude and derived the ratio of the haze scale height (H_haze) to the pressure scale height (H_gas) as a function of altitude. We inferred H_haze/H_gas values varying from 0.8 to 2.4. The aerosol scale height varies with altitude and also with latitude. Overall, the haze extinction does not show strong latitudinal variations but, at 1 mbar, an increase by a factor of 1.5 is observed at the north pole compared to high southern latitudes. The vertical optical depths at 0.5 and 1.7 mbar increase from 55°S to 5°N, remain constant between 5°N and 30°N and display little variation at higher latitudes, except the presence of a slight local maximum at 45°N. The spectral dependence of the haze vertical optical depth is uniform with latitude and displays three main spectral features centered at 630 cm⁻¹, 745 cm⁻¹ and 1390 cm⁻¹, the latter showing a wide tail extending down to ∼1000 cm⁻¹. From 600 to 750 cm⁻¹, the optical depth increases by a factor of 3 in contrast with the absorbance of laboratory tholins, which is generally constant. We derived the mass mixing ratio profiles of haze at the nine latitudes. Below the 0.4-mbar level all mass mixing ratio profiles increase with height. Above this pressure level, the profiles at 40°N, 45°N, 55°N, at the edge of the polar vortex, display a decrease-with-height whereas the other profiles increase. The global increase with height of the haze mass mixing ratio suggest a source at high altitudes and a sink at low altitudes. An enrichment of haze is observed at 0.1 mbar around the equator, which could be due to a more efficient photochemistry because of the strongest insolation there or an accumulation of haze due to a balance between sedimentation and upward vertical drag.     

Latitudinal transport by barotropic waves in Titan's stratosphere.: II. Results from a coupled dynamics-microphysics-photochemistry GCM

We present a 2D general circulation model of Titan's atmosphere, coupling axisymmetric dynamics with haze microphysics, a simplified photochemistry and eddy mixing. We develop a parameterization of latitudinal eddy mixing by barotropic waves based on a shallow-water, longitude-latitude model. The parameterization acts locally and in real time both on passive tracers and momentum. The mixing coefficient varies exponentially with a measure of the barotropic instability of the mean zonal flow. The coupled GCM approximately reproduces the Voyager temperature measurements and the latitudinal contrasts in the distributions of HCN and C₂H₂, as well as the main features of the zonal wind retrieved from the 1989 stellar occultation. Wind velocities are consistent with the observed reversal time of the North-South albedo asymmetry of 5 terrestrial years. Model results support the hypothesis of a non-uniform distribution of infrared opacity as the cause of the Voyager temperature asymmetry. Transport by the mean meridional circulation, combined with polar vortex isolation may be at the origin of the latitudinal contrasts of trace species, with eddy mixing remaining restricted to low latitudes most of the Titan year. We interpret the contrasts as a signature of non-axisymmetric motions.     

Products of meteoric metal ion chemistry within planetary atmospheres. 1. Mg at Titan

We report results of quantum chemical calculations of Mg⁺/ligand bond dissociation energies involving ligands identified as major constituents of Titan's upper atmosphere. Trends identified in these results allow elucidation of the important bimolecular and termolecular reactions of Mg⁺, and of simple molecular ions containing Mg⁺, arising from meteoric infall into Titan's atmosphere. Our study highlights, and includes calculated rate coefficients for, crucial ligand-switching and ligand-stripping reactions which ensure that a dynamic equilibrium exists between atomic and molecular ions of Mg⁺. Neutralization of ionized meteoric Mg is expected to produce the radical MgNC in high yield. The highly polar MgNC radical should provide an excellent nucleation site for condensation of polar (e.g., HCN, CH₃CN, and HC₃N) and highly unsaturated (e.g., C₂H₂, C₄H₂, and C₂N₂) neutrals at comparatively high altitude, leading to precipitation of Mg-doped tholin-like material. The implications for Titan's prebiotic chemical evolution, of the surface deposition of such material (which may feasibly contain magnesium porphyrins, or other bioactive Mg-containing complexes) remain to be assessed.     

Analysis of Titan's neutral upper atmosphere from Cassini Ion Neutral Mass Spectrometer measurements

In this paper we present an in-depth study of the distributions of various neutral species in Titan's upper atmosphere, between 950 and 1500 km for abundant species (N₂, CH₄, H₂) and between 950 and 1200 km for other minor species. Our analysis is based on a large sample of Cassini/INMS (Ion Neutral Mass Spectrometer) measurements in the CSN (Closed Source Neutral) mode, obtained during 15 close flybys of Titan. To untangle the overlapping cracking patterns, we adopt Singular Value Decomposition (SVD) to determine simultaneously the densities of different species. Except for N₂, CH₄, H₂ and ⁴⁰Ar (as well as their isotopes), all species present density enhancements measured during the outbound legs. This can be interpreted as a result of wall effects, which could be either adsorption/desorption of these molecules or heterogeneous surface chemistry of the associated radicals on the chamber walls. In this paper, we provide both direct inbound measurements assuming ram pressure enhancement only and abundances corrected for wall adsorption/desorption based on a simple model to reproduce the observed time behavior. Among all minor species of photochemical interest, we have firm detections of C₂H₂, C₂H₄, C₂H₆, CH₃C₂H, C₄H₂, C₆H₆, CH₃CN, HC₃N, C₂N₂ and NH₃ in Titan's upper atmosphere. Upper limits are given for other minor species.The globally averaged distributions of N₂, CH₄ and H₂ are each modeled with the diffusion approximation. The N₂ profile suggests an average thermospheric temperature of 151 K. The CH₄ and H₂ profiles constrain their fluxes to be and , referred to Titan's surface. Both fluxes are significantly higher than the Jeans escape values. The INMS data also suggest horizontal/diurnal variations of temperature and neutral gas distribution in Titan's thermosphere. The equatorial region, the ramside, as well as the nightside hemisphere of Titan appear to be warmer and present some evidence for the depletion of light species such as CH₄. Meridional variations of some heavy species are also observed, with a trend of depletion toward the north pole. Though some of the above variations might be interpreted by either the solar-driven models or auroral-driven models, a physical scenario that reconciles all the observed horizontal/diurnal variations in a consistent way is still missing. With a careful evaluation of the effect of restricted sampling, some of the features shown in the INMS data are more likely to be observational biases.     

Clouds on Titan during the Cassini prime mission: A complete analysis of the VIMS data

We use data from the VIMS instrument on board the Cassini spacecraft to construct high sensitivity and high spatial-resolution maps of the locations of tropospheric clouds on Titan in the late northern winter season during which the Cassini prime mission took place. These observations show that, in this season, clouds on Titan are strongly hemispherically asymmetric. Mid-latitude clouds, in particular, occur only in the southern hemisphere and have not ever been observed in the north. Such an asymmetry is in general agreement with circulation models where sub-solar surface heating controls the locations of clouds and appears in conflict with models where perennial polar hazes prevent significant summertime polar heating from affecting the circulation. The southern mid-latitude clouds appear to be distributed uniformly in longitude, in contrast to some previous observations. Southern high-latitude clouds exhibit a significant concentration, however, between about 180° and 270°E longitude. A spatially and temporally uniform cloud always appears northward of ∼50°N latitude. This cloud appears unchanged over the course of the observations, consistent with the interpretation that it is caused by continuous ethane condensation as air subsides and radiatively cools through the tropopause. The location of this cloud likely provides a direct tracer of elements of north polar atmospheric circulation, potentially allowing continuous monitoring of circulation changes as Titan passes through equinox into north polar spring and summer. We show that a similar analysis of this dataset by Rodriguez et al. (2009) contains substantial errors and should not be used.     

Latitudinal transport by barotropic waves in Titan's stratosphere.: I. General properties from a horizontal shallow-water model

We present a numerical study of barotropic waves in Titan's stratosphere based on a shallow-water model. The forcing of the zonal flow by the mean meridional circulation is represented by a relaxation towards a barotropically unstable wind profile. The relaxation profile is consistent with observations and with previous results from a 3D general circulation model. The time constant of the forcing that best matches the northward eddy-transport of zonal momentum from the 3D model is τ∼5 Titan days. The eddy wind field is a zonal wavenumber-2 wave with a peak amplitude about 10% of the mean wind speed. The latitudinal transport of angular momentum by the wave tends to keep the flow close to marginal stability by carrying momentum upgradient, from the core of the jets into the low latitudes. Although the strongest eddy motions occur at the latitudes of the wind maxima, the strongest mixing takes place at the barotropically unstable regions, close to ±30° and spanning about 30° in latitude. An eddy-mixing time constant of the order of 1 Titan day is inferred within these regions, and of a few tens of days within regions of stable flow. Horizontal gradients in transient tracer fields are less than 10% of the latitudinal gradient of the meridional tracer profile. Cassini's detection of such waves could provide a direct observation of wind speeds at stratospheric levels.     

The near-IR spectrum of Titan modeled with an improved methane line list

We have obtained spatially resolved spectra of Titan in the near-infrared J, H and K bands at a resolving power of ∼5000 using the near-infrared integral field spectrometer (NIFS) on the Gemini North 8 m telescope. Using recent data from the Cassini/Huygens mission on the atmospheric composition and surface and aerosol properties, we develop a multiple-scattering radiative transfer model for the Titan atmosphere. The Titan spectrum at these wavelengths is dominated by absorption due to methane with a series of strong absorption band systems separated by window regions where the surface of Titan can be seen. We use a line-by-line approach to derive the methane absorption coefficients. The methane spectrum is only accurately represented in standard line lists down to ∼2.1 μm. However, by making use of recent laboratory data and modeling of the methane spectrum we are able to construct a new line list that can be used down to 1.3 μm. The new line list allows us to generate spectra that are a good match to the observations at all wavelengths longer than 1.3 μm and allow us to model regions, such as the 1.55 μm window that could not be studied usefully with previous line lists such as HITRAN 2008. We point out the importance of the far-wing line shape of strong methane lines in determining the shape of the methane windows. Line shapes with Lorentzian, and sub-Lorentzian regions are needed to match the shape of the windows, but different shape parameters are needed for the 1.55 μm and 2 μm windows. After the methane lines are modeled our observations are sensitive to additional absorptions, and we use the data in the 1.55 μm region to determine a D/H ratio of 1.77 ± 0.20 × 10⁻⁴, and a CO mixing ratio of 50 ± 11 ppmv. In the 2 μm window we detect absorption features that can be identified with the ν₅ + 3ν₆ and 2ν₃ + 2ν₆ bands of CH₃D.     

The role of organic haze in Titan's atmospheric chemistry: I. Laboratory investigation on heterogeneous reaction of atomic hydrogen with Titan tholin

In Titan's atmosphere consisting of N₂ and CH₄, large amounts of atomic hydrogen are produced by photochemical reactions during the formation of complex organics. This atomic hydrogen may undergo heterogeneous reactions with organic aerosol in the stratosphere and mesosphere of Titan. In order to investigate both the mechanisms and kinetics of the heterogeneous reactions, atomic deuterium is irradiated onto Titan tholin formed from N₂ and CH₄ gas mixtures at various surface-temperatures of the tholin ranging from 160 to 310 K. The combined analyses of the gas species and the exposed tholin indicate that the interaction mechanisms of atomic deuterium with the tholin are composed of three reactions; (a) abstraction of hydrogen from tholin resulting in gaseous HD formation (HD recombination), (b) addition of D atom into tholin (hydrogenation), and (c) removal of carbon and/or nitrogen (chemical erosion). The reaction probabilities of HD recombination and hydrogenation are obtained as ηabst=1.9(±0.6)×¹⁰⁻³×exp(−300/T) and ηhydro=2.08(±0.64)×exp(−1000/T), respectively. The chemical erosion process is very inefficient under the conditions of temperature range of Titan's stratosphere and mesosphere. Under Titan conditions, the rates of hydrogenation > HD recombination ? chemical erosion. Our measured HD recombination rate is about 10 times (with an uncertainty of a factor of 3-5) the prediction of previous theoretical model. These results imply that organic aerosol can remove atomic hydrogen efficiently from Titan's atmosphere through the heterogeneous reactions and that the presence of aerosol may affect the subsequent organic chemistry.