Measurements of methane absorption by the descent imager/spectral radiometer (DISR) during its descent through Titan's atmosphere

New low-temperature methane absorption coefficients pertinent to the Titan environment are presented as derived from the Huygens DISR spectral measurements combined with the in-situ measurements of the methane gas abundance profile measured by the Huygens Gas Chromatograph/Mass Spectrometer (GCMS). The visible and near-infrared spectrometers of the descent imager/spectral radiometer (DISR) instrument on the Huygens probe looked upward and downward covering wavelengths from 480 to 1620 nm at altitudes from 150 km to the surface during the descent to Titan's surface. The measurements at continuum wavelengths were used to determine the vertical distribution, single-scattering albedos, and phase functions of the aerosols. The gas chromatograph/mass spectrometer (GCMS) instrument on the probe measured the methane mixing ratio throughout the descent. The DISR measurements are the first direct measurements of the absorbing properties of methane gas made in the atmosphere of Titan at the pathlengths, pressures, and temperatures that occur there. Here we use the DISR spectral measurements to determine the relative methane absorptions at different wavelengths along the path from the probe to the sun throughout the descent. These transmissions as functions of methane path length are fit by exponential sums and used in a haze radiative transfer model to compare the results to the spectra measured by DISR. We also compare the recent laboratory measurements of methane absorption at low temperatures [Irwin et al., 2006. Improved near-infrared methane band models and k-distribution parameters from 2000 to 9500 cm⁻¹ and implications for interpretation of outer planet spectra. Icarus 181, 309-319] with the DISR measurements. We find that the strong bands formed at low pressures on Titan act as if they have roughly half the absorption predicted by the laboratory measurements, while the weak absorption regions absorb considerably more than suggested by some extrapolations of warm measurements to the cold Titan temperatures. We give factors as a function of wavelength that can be used with the published methane coefficients between 830 and 1620 nm to give agreement with the DISR measurements. We also give exponential sum coefficients for methane absorptions that fit the DISR observations. We find the DISR observations of the weaker methane bands shortward of 830 nm agree with the methane coefficients given by Karkoschka [1994. Spectrophotometry of the jovian planets and Titan at 300- to 1000-nm wavelength: the methane spectrum. Icarus 111, 174-192]. Finally, we discuss the implications of our results for computations of methane absorption in the atmospheres of the outer planets.     

Cassini RADAR constraint on Titan's winter polar precipitation

The Cassini RADAR instrument made a dedicated cloud backscatter observation near Titan's north pole, presently in winter darkness, to constrain the precipitation of material onto the surface. The detection limit is ∼5 orders of magnitude above that expected in methane rainstorms, and rules out ‘drizzle’ of more than , placing constraints on the winter accumulation of material on Titan's surface during polar winter.     

On the maintenance of thermal wind balance and equatorial superrotation in Titan's stratosphere

Titan's atmospheric winds, like those on Venus, exhibit superrotation at high altitudes. Titan general circulation models have yielded conflicting results on whether prograde winds in excess of 100 m s⁻¹ at the 1 mbar level are possible based on known physical processes that drive wind systems. A comprehensive two-dimensional (2D) model for Titan's stratosphere was constructed to systematically explore the physical mechanisms that produce and maintain stratospheric wind systems. To ensure conservation of angular momentum in the limit of no net exchange of atmospheric angular momentum with the solid satellite and no external sources and sinks, the zonal momentum equation was solved in flux form for total angular momentum. The relationships among thermal wind balance, meridional circulation, and zonal wind were examined with numerical experiments over a range of values for fundamental input parameters, including planetary rotation rate, radius, internal friction due to wave stresses, and net radiative drive. The magnitude of mid-latitude jets is most sensitive to a single parameter, the planetary rotation rate and results from the conversion of planetary angular momentum to relative angular momentum by the meridional circulation, whereas the strength of meridional circulation is mainly determined by the magnitude of the radiative drive. For Titan's slowly rotating atmosphere, the meridional temperature gradient is vanishingly small, even when the radiative drive is enhanced beyond reasonable magnitudes, and can be inferred from zonal winds in gradient/thermal wind balance. In our 2D model large equatorial superrotation in Titan's stratosphere can be only produced through internal drag forcing by eddy momentum fluxes, which redistribute angular momentum within the atmosphere, while still conserving the total angular momentum of the atmosphere with time. We cannot identify any waves, such as gravitational or thermal tides, that are sufficiently capable of generating the required eddy forcing of >50 m s⁻¹ Titan-day⁻¹ to maintain peak prograde winds in excess of 100 m s⁻¹ at the 1 mbar level.     

Dissipation of Titan's south polar clouds

Nearly all adaptive optics images of Titan taken between December 2001 and November 2004 showed tropospheric clouds located within 30° of the south pole. We report here on a dissipation of Titan's south polar clouds observed in twenty-nine Keck and Gemini images taken between December 2004 and April 2005. The near complete lack of south polar cloud activity during this time, and subsequent resurgence months later at generally higher latitudes, may be the beginning of seasonal change in Titan's weather. The ∼5 month decrease in cloud activity may also have been caused by methane rainout from a large cloud event in October 2004. Understanding the seasonal evolution of Titan's clouds, and of any precipitation associated with them, is essential for interpreting the geological observations of fluid flow features observed over a wide range of Titan latitudes with the Cassini/Huygens spacecraft.     

A large cloud outburst at Titan's south pole

Images of Titan acquired over five nights in October 2004 using the adaptive optics system at the Keck Observatory show dramatic increases in tropospheric cloud activity at the south pole compared with all other images of Titan clouds to date. During this time, Titan's south polar clouds brightened to more than 18 times their typical values. The Cassini Ta flyby of Titan occurred as this storm was rapidly dissipating. We find that the brightness of this cloud outburst is consistent with the dramatic transient brightening of Titan observed in atmospheric windows on two nights in 1995 by Griffith et al. [Griffith, C.A., Owen, T., Miller, G.A., Geballe, T., 1998. Nature 395 (6702) 575-578] if we scale the brightness of the cloud by projecting it onto the equator. While apparently infrequent, the fact that large cloud events have been observed in different seasons of Titan's year indicates that these large storms might be a year-round phenomenon on Titan. We propose possible mechanisms to explain these occasional short-term increases in Titan's cloud activity.     

Shading under Titan's sky

During the descent of the Huygens probe in January 2005, its Descent Imager/Spectral Radiometer (DISR) will take the first close up images of Titan's surface. The shading imposed by the illumination of a planetary surface contains information on its topography. For planetary bodies without an optically thick atmosphere, the light can be assumed to stem from a point source. In this case, methods are available in order to estimate the shape of surface features from shading. The situation is quite different for Titan, as its atmosphere is optically thick at optical wavelengths. The sun is visible from the surface, but the illumination is dominated by diffuse radiance. In order to investigate the characteristics of shading under Titan's sky and to assess methods to retrieve the shape, different digital terrain models (DTMs) are used to simulate images according to different types of illumination. For an idealized DTM, the shape is retrieved from the shading in the simulated images. Deriving the shape from shading under Titan's sky using existing methods is only possible if the topography is relatively flat, i.e. in the absence of steep slopes.     

New laboratory measurements of CH4 in Titan's conditions and a reanalysis of the DISR near-surface spectra at the Huygens landing site

Laboratory spectra of methane-nitrogen mixtures have been recorded in the near-infrared range (1.0-1.65 μm) in conditions similar to Titan's near surface, to facilitate the interpretation of the DISR/DLIS (DISR—Descent Imager/Spectral Radiometer) spectra taken during the last phase of the descent of the Huygens Probe, when the surface was illuminated by a surface-science lamp. We used a 0.03 cm⁻¹ spectral resolution, adequate to resolve the lines at high pressure (p_N2∼1.5 bar). By comparing the laboratory spectra with synthetic calculations in the well-studied ν₂+2ν₃ band (7515-7620 cm⁻¹), we determine a methane absorption column density of 178±20 cm atm and a temperature of 118±10 K in our experiment. From this, we derive the methane absorption coefficients over 1.0-1.65 μm with a 0.03 cm⁻¹ sampling, allowing for the extrapolation of the results to any other methane column density under the relevant pressure and temperature conditions. We then revisit the calibration and analysis of the Titan “lamp-on” DLIS spectra. We infer a 5.1±0.8% methane-mixing ratio in the first 25 m of Titan's atmosphere. The CH₄ mixing ratio measured 90 s after landing from a distance of 45 cm is found to be 0.92±0.25 times this value, thus showing no post-landing outgassing of methane in excess of ∼20%. Finally, we determine the surface reflectivity as seen between 25 m and 45 cm and find that the 1500 nm absorption band is deeper in the post-landing spectrum as compared to pre-landing.     

Heat balance in Titan's atmosphere

The recent measurements of the vertical distribution and optical properties of haze aerosols as well as of the absorption coefficients for methane at long paths and cold temperatures by the Huygens entry probe of Titan permit the computation of the solar heating rate on Titan with greater certainty than heretofore. We use the haze model derived from the Descent Imager/Spectral Radiometer (DISR) instrument on the Huygens probe [Tomasko, M.G., Doose, L., Engel, S., Dafoe, L.E., West, R., Lemmon, M., Karkoschka, E., See, C., 2008a. A model of Titan's aerosols based on measurements made inside the atmosphere. Planet. Space Sci., this issue, doi:10.1016/j.pss.2007.11.019] to evaluate the variation in solar heating rate with altitude and solar zenith angle in Titan's atmosphere. We find the disk-averaged solar energy deposition profile to be in remarkably good agreement with earlier estimates using very different aerosol distributions and optical properties. We also evaluated the radiative cooling rate using measurements of the thermal emission spectrum by the Cassini Composite Infrared Spectrometer (CIRS) around the latitude of the Huygens site. The thermal flux was calculated as a function of altitude using temperature, gas, and haze profiles derived from Huygens and Cassini/CIRS data. We find that the cooling rate profile is in good agreement with the solar heating profile averaged over the planet if the haze structure is assumed the same at all latitudes. We also computed the solar energy deposition profile at the 10°S latitude of the probe-landing site averaged over one Titan day. We find that some 80% of the sunlight that strikes the top of the atmosphere at this latitude is absorbed in all, with 60% of the incident solar energy absorbed below 150 km, 40% below 80 km, and 11% at the surface at the time of the Huygens landing near the beginning of summer in the southern hemisphere. We compare the radiative cooling rate with the solar heating rate near the Huygens landing site averaging over all longitudes. At this location, we find that the solar heating rate exceeds the radiative cooling rate by a maximum of 0.5 K/Titan day near 120 km altitude and decreases strongly above and below this altitude. Since there is no evidence that the temperature structure at this latitude is changing, the general circulation must redistribute this heat to higher latitudes.     

Titan's 5-micron lightcurve

We report on repeated mid-resolution (R∼2000) spectroscopic observations of Titan, acquired between November 2002 and January 2003 with ISAAC at the ESO/VLT and covering the 4.84-5.05 μm range. These observations, which sample four different positions of Titan around Saturn, clearly indicate a variability of the 5-μm continuum albedo, with Titan's geometric albedo decreasing by ∼40% from Titan's leading side to the trailing side. This Titan 5-μm “lightcurve” appears to be in phase with the other near-infrared lightcurves. This can be understood in terms of a surface model in which water ice coexists in minor and variable proportions (13-25%, if pure) with a second, dark, component.     

The detached haze layer in Titan's mesosphere

By comparing observations from the Cassini imaging system, UV spectrometer, and Huygens atmospheric structure instrument, we determine an apparent radius of ∼40 nm, an imaginary index <0.3 at 187.5 nm and a number density of ∼30 particles cm⁻³ for the detached haze layer at 520 km in Titan's mesosphere. We point out that the detached haze layer is coincident with a local maximum in the measured temperature profile and show that the temperature maximum is caused by absorption of sunlight in the detached haze layer. This rules out condensation as the source of the layer. The derived particle size is in good agreement with that estimated for the size of the monomers in the aggregate particles that make up the main haze layer. Calculations of the sedimentation velocity of the haze particles coupled with the derived number density imply a mass flux , which is approximately equal to the mass flux required to explain the main haze layer. Because the aerosol size and mass flux derived for the detached layer agree with those determined for the main layer, we suggest that the main haze layer in Titan's stratosphere is formed primarily by sedimentation and coagulation of particles in the detached layer. This implies that high-energy radical and ion chemistry in the thermosphere is the main source of haze on Titan.     

Titan's atmospheric engine: an overview

Heating occurs in Titan's stratosphere from the absorption of incident solar radiation by methane and aerosols. About 10% of the incident sunlight reaches Titan's surface and causes heating there. Thermal radiation redistributes heat within the atmosphere and cools to space. The resulting vertical temperature profile is stable against convection and a state of radiative equilibrium is established. Equating theoretical and observed temperature profiles enables an empirical determination of the vertical distribution of thermal opacity. A uniformly mixed aerosol is responsible for most of the opacity in the stratosphere, whereas collision-induced absorption of gases is the main contributor in the troposphere. Occasional clouds are observed in the troposphere in spite of the large degrees of methane supersaturation found there. Photochemistry converts CH₄ and N₂ into more complex hydrocarbons and nitriles in the stratosphere and above. Thin ice clouds of trace organics are formed in the winter and early spring polar regions of the lower stratosphere. Precipitating ice particles serve as condensation sites for supersaturated methane vapor in the troposphere below, resulting in lowered methane degrees of supersaturation in the polar regions. Latitudinal variations of stratospheric temperature are seasonal, and lag instantaneous response to solar irradiation by about one season for two reasons: (1) an actual instantaneous thermal response to a latitudinal distribution of absorbing gases, themselves out of phase with the sun by about one season, and (2) a sluggish dynamical response of the stratosphere to the latitudinal transport of angular momentum, induced by radiative heating and cooling. Mean vertical abundances of stratospheric organics and aerosols are determined primarily by atmospheric chemistry and condensation, whereas latitudinal distributions are more influenced by meridional circulations. In addition to preferential scavenging by precipitating ice particles from above, the polar depletion of supersaturated methane results from periodic scavenging by short-lived tropospheric clouds, coupled with the steady poleward march of the continuously drying atmosphere due to meridional transport.     

Compositional evidence for Titan's stratospheric tilt

Five years of Cassini CIRS infrared spectra have been used to determine the tilt of Titan's stratospheric symmetry axis with respect to the solid body rotation axis. Measurements of HCN abundance centred around 5 mbar (125 km altitude) at equatorial latitudes show the symmetry axis is tilted by 4.0±1.5^° in a direction 70±40^°W of the sub-solar point. This value is consistent with tilts determined from temperature and haze measurements by Achterberg et al. (2008a) and Roman et al. (2009). The consistency of results from three independent methods suggests that Titan's entire stratosphere is tilted and provides a powerful constraint on the underlying atmospheric dynamics.     

A spherical model for computing polarized radiation in Titan's atmosphere

The Huygens descent through Titan's atmosphere in January 2005 will provide invaluable information about Titan's atmospheric composition and aerosol properties. The Descent Imager/Spectral Radiometer (DISR) will perform upward and downward looking radiation observations at various spectral ranges and spatial resolutions. To prepare the DISR data interpretation we have developed a new model for radiation transfer in Titan's atmosphere. The model solves for the full three-dimensional polarized radiation field in spherical geometry. However, the atmosphere itself is assumed to be spherically symmetric. The model is initialized with a fast-to-compute plane–parallel solution based on the doubling and adding algorithm that incorporates a spherical correction for the incoming direct solar beam. The full three-dimensional problem is then solved using the characteristics method combined with the Picard iterative approximation as described in Rozanov et al. (J. Quant. Spectrosc. Radiat. Transfer 69 (2001) 491). Aerosol scattering properties are calculated with a new microphysical model. In this formulation, aerosols are assumed to be fractal aggregates and include methane gas absorption embedded into the extinction coefficient. The resulting radiance of the model atmosphere's internal field is presented for two prescribed DISR wavelengths.     

Laboratory studies of methane and ethane adsorption and nucleation onto organic particles: Application to Titan's clouds

Titan, Saturn's largest moon, has a thick nitrogen/methane atmosphere. The temperature and pressure conditions in Titan's atmosphere are such that the methane vapor should condense near the tropopause to form clouds. Several ground-based measurements have observed sparse cloud-like features in Titan's atmosphere, while the Cassini mission to Saturn has provided large scale images of the clouds. However, Titan's cloud formation conditions remain poorly constrained. Heterogeneous nucleation (from the vapor phase onto a solid or liquid aerosol surface) greatly enhances cloud formation relative to homogeneous nucleation. In order to elucidate the cloud formation mechanism near the tropopause, we have performed laboratory measurements of the adsorption of methane and ethane onto solid organic particles (tholins) representative of Titan's photochemical haze. We find that monolayers of methane adsorb onto tholin particles at saturation ratios less than unity. We also find that solid methane nucleates onto the adsorbed methane at a saturation ratio of S=1.07±0.008. This implies that Titan's methane clouds should form easily. This is consistent with recent measurements of the column of methane ruling out excessive methane supersaturation. In addition, we find ethane adsorbs onto tholin particles in a metastable phase prior to nucleation. However, ethane nucleation onto the adsorbed ethane occurs at a relatively high saturation ratio of S=1.36±0.08. These findings are consistent with the recent report of polar ethane clouds in Titan's lower stratosphere.     

Westward rotation of the atmospheric angular momentum vector of Titan by thermal tides

The orientation of the atmospheric angular momentum vector of Titan and its temporal variation predicted by a general circulation model are analysed and interpreted. The atmospheric angular momentum vector is tilted by a few degrees from the polar axis and the vector rotates (precesses) westward with a constant period of 1 Titan day. The fast westward rotation is likely to be caused by migrating diurnal thermal tides. The tilt is almost cancelled out in the troposphere by the wavenumber 2 pattern of Saturn's gravitational tide, but is more pronounced in the stratosphere, where thermal tides are significant. The predicted tilt angle and the equatorial angular momentum vary with season and maximize when the hemispheric asymmetry of the axial angular momentum or superrotation attains its peak.     

Condensate clouds in Titan's north polar stratosphere

Analysis of the 250-560 cm⁻¹ spectral continuum of Titan's north polar hood just after spring equinox reveals, in addition to the ubiquitous aerosol, a tenuous but relatively uniform cloud of small particles permeating the lower stratosphere at altitudes between 58 and 90 km. Voyager 1 IRIS data suggest the particles are highly scattering, almost certainly condensed organics, with radii between 1 and 5 μm. Mole fractions for the condensed material range between 4×¹⁰⁻⁸ and 4×¹⁰⁻⁶, depending upon particle size. Vapor pressure arguments imply condensed nitriles near 90 km, the most likely being HCN, with condensed hydrocarbons such as C₂H₆ restricted to regions considerably nearer the tropopause. No direct chemical identification is possible. Negligible methane supersaturation in the troposphere at 67.4° N latitude, when compared with degrees of supersaturation at other latitudes, hints at precipitation fluxes of north polar stratospheric condensates during the previous northern winter that were perhaps three orders of magnitude greater than those at low latitudes during that time. A scale height of 1.5 times the density scale height above 160 km is reaffirmed for the photochemical aerosol of the north polar hood. There appears to be a depletion of aerosol somewhere below 160 km. An aerosol mole fraction ∼8×¹⁰⁻⁸ at 160 km is inferred, about 33% greater than the value derived in a previous study. The Cassini CIRS instrument, with its expanded spectral range and higher spectral resolution, should be able to provide highly complementary information for the time period covering most of the northern winter season.     

A coupled dynamics-microphysics model of Titan's atmosphere

We have developed a coupled general circulation model of Titan's atmosphere in which the aerosol haze is treated with a microphysical model and is advected by the winds. The radiative transfer accounts for the non uniform haze distribution and, in turn, drives the dynamics. We analyze the GCM results, especially focusing on the difference between a uniform haze layer and a haze layer coupled to the dynamics. In the coupled simulation the aerosols tend to accumulate at the poles, at latitudes higher than ±60°. During winter, aerosols strongly radiate at thermal infrared wavelengths enhancing the cooling rate near the pole. Since this tends to increase the latitudinal gradients of temperature the direct effect of this cooling excess, in contrast to the uncoupled haze case, is to increase the strength of the meridional cells as well as the strength of the zonal winds and profile. This is a positive feedback of the haze on dynamics. The coupled model reproduces observations about the state of the atmosphere better than the uniform haze model, and in addition, the northern polar hood and the detached haze are qualitatively reproduced.     

Optical properties of aerosols in Titan's atmosphere

In the frame of fractal modeling of tholin aggregates we made a systematic analysis of their optical properties. Ballistic particle-cluster aggregation (BPCA) and diffusion-limited aggregation (DLA) of spherical primary particles (monomers) identical in material composition were considered. Aggregates composed of identical particles (monodisperse cluster), as well as of size-distributed particles (polydisperse cluster), were simulated. To calculate the light-scattering models, the code based on the superposition T-matrix method is used. Orientationally averaged properties of light scattering by model particles were extracted, and the normalized phase function and the degree of linear polarization were calculated as functions of scattering angle. We concluded that: (a) aggregation mechanism as well as specific internal structure of the clusters play only a minor role, and for the future it is not necessary to investigate aggregates of different types; (b) the intensity is very sensitive both to the size parameter of forming particles x and to the size parameter of the aggregates X; (c) characterization of the aggregates by the number of monomers is insufficient to retrieve physical properties of aggregates from optical measurement; and (d) it is very desirable to include into the analysis polarization data calculated for the different clusters.