A model of Titan's aerosols based on measurements made inside the atmosphere

The descent imager/spectral radiometer (DISR) instrument aboard the Huygens probe into the atmosphere of Titan measured the brightness of sunlight using a complement of spectrometers, photometers, and cameras that covered the spectral range from 350 to 1600 nm, looked both upward and downward, and made measurements at altitudes from 150 km to the surface. Measurements from the upward-looking visible and infrared spectrometers are described in Tomasko et al. [2008a. Measurements of methane absorption by the descent imager/spectral radiometer (DISR) during its descent through Titan's atmosphere. Planet. Space Sci., this volume]. Here, we very briefly review the measurements by the violet photometers, the downward-looking visible and infrared spectrometers, and the upward-looking solar aureole (SA) camera. Taken together, the DISR measurements constrain the vertical distribution and wavelength dependence of opacity, single-scattering albedo, and phase function of the aerosols in Titan's atmosphere.Comparison of the inferred aerosol properties with computations of scattering from fractal aggregate particles indicates the size and shape of the aerosols. We find that the aggregates require monomers of radius 0.05 μm or smaller and that the number of monomers in the loose aggregates is roughly 3000 above 60 km. The single-scattering albedo of the aerosols above 140 km altitude is similar to that predicted for some tholins measured in laboratory experiments, although we find that the single-scattering albedo of the aerosols increases with depth into the atmosphere between 140 and 80 km altitude, possibly due to condensation of other gases on the haze particles. The number density of aerosols is about 5/cm³ at 80 km altitude, and decreases with a scale height of 65 km to higher altitudes. The aerosol opacity above 80 km varies as the wavelength to the −2.34 power between 350 and 1600 nm.Between 80 and 30 km the cumulative aerosol opacity increases linearly with increasing depth in the atmosphere. The total aerosol opacity in this altitude range varies as the wavelength to the −1.41 power. The single-scattering phase function of the aerosols in this region is also consistent with the fractal particles found above 60 km.In the lower 30 km of the atmosphere, the wavelength dependence of the aerosol opacity varies as the wavelength to the −0.97 power, much less than at higher altitudes. This suggests that the aerosols here grow to still larger sizes, possibly by incorporation of methane into the aerosols. Here the cumulative opacity also increases linearly with depth, but at some wavelengths the rate is slightly different than above 30 km altitude.For purely fractal particles in the lowest few km, the intensity looking upward opposite to the azimuth of the sun decreases with increasing zenith angle faster than the observations in red light if the single-scattering albedo is assumed constant with altitude at these low altitudes. This discrepancy can be decreased if the single-scattering albedo decreases with altitude in this region. A possible explanation is that the brightest aerosols near 30 km altitude contain significant amounts of methane, and that the decreasing albedo at lower altitudes may reflect the evaporation of some of the methane as the aerosols fall into dryer layers of the atmosphere. An alternative explanation is that there may be spherical particles in the bottom few kilometers of the atmosphere.     

Simulations of Titan's brightness by a two-dimensional haze model

We have used a 2-D microphysics model to study the effects of atmospheric motions on the albedo of Titan's thick haze layer. We compare our results to the observed variations of Titan's brightness with season and latitude. We use two wind fields; the first is a simple pole-to-pole Hadley cell that reverses twice a year. The second is based on the results of a preliminary Titan GCM. Seasonally varying wind fields, with horizontal velocities of about 1 cm sec-1 at optical depth unity, are capable of producing the observed change in geometric albedo of about 10% over the Titan year. Neither of the two wind fields can adequately reproduce the latitudinal distribution of reflectivity seen by Voyager. At visible wavelengths, where only haze opacity is important, upwelling produces darkening by increasing the particle size at optical depth unity. This is due to the suspension of larger particles as well as the lateral removal of smaller particles from the top of the atmosphere. At UV wavelengths and at 0.89 micrometers the albedo is determined by the competing effects of the gas the haze material. Gas is bright in the UV and dark at 0.89 micrometers. Haze transport at high altitudes controls the UV albedo and transport at low altitude controls the 0.89 micrometers albedo. Comparisons between the hemispheric contrast at UV, visible, and IR wavelengths can be diagnostic of the vertical structure of the wind field on Titan.     

Influence of high abundances of aerosols on the electrical conductivity of the Titan atmosphere

Observations of optical depth and scattering by instrumentation onboard the Huygens probe have been used by Tomasko et al. [Tomasko et al., 2005. Rain, winds and haze during Huygens probe's descent to Titan's surface. Nature 438 (8), 765-778] to deduce that the size and abundance of Titan aerosols could be nearly independent of altitude. Here we show that by assuming a constant mass flux with altitude and using the measured optical depth as a constraint, we obtain more realistic size and abundance distributions. In particular, the calculated abundance decreases from 3.5×10⁷ m⁻³ at 100 km to 8×10⁶ m⁻³ near the surface while the particle radius varies from 0.25 μm at 150 km to 1.1 μm at the surface. These distributions are consistent with the reported measurements for these quantities. Our results are then employed to compute electron and ion densities and conductivities for various solar UV photoelectron emission thresholds. Our model shows that to get agreement with the published (preliminary) conductivity measurements, photoemission cannot be an important source of electrons and ions. To get agreement with the electron and ion conductivity observations, both an additional population of aerosol embryos above 50 km and a trace amount of an electrophillic molecular species below 50 km are needed.     

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.     

Limits on the size of aerosols from measurements of linear polarization in Titan’s atmosphere

The Descent Imager/Spectral Radiometer (DISR) instrument on the Huygens probe into the atmosphere of Titan yielded information on the size, shape, optical properties, and vertical distribution of haze aerosols in the atmosphere of Titan [Tomasko, M.G., Doose, L., Engel, S., Dafoe, L.E., West, R., Lemmon, M., Karkoschka, E., 2008. Planet. Space Sci. 56, 669-707] from photometric and spectroscopic measurements of sunlight in Titan’s atmosphere. This instrument also made measurements of the degree of linear polarization of sunlight in two spectral bands centered at 491 and 934 nm. Here we present the calibration and reduction of the polarization measurements and compare the polarization observations to models using fractal aggregate particles which have different sizes for the small dimension (monomer size) of which the aggregates are composed. We find that the Titan aerosols produce very large polarizations perpendicular to the scattering plane for scattering near 90° scattering angle. The size of the monomers is tightly constrained by the measurements to a radius of 0.04 ± 0.01 μm at altitudes from 150 km to the surface. The decrease in polarization with decreasing altitude observed in red and blue light is as expected by increasing dilution due to multiple scattering at decreasing altitudes. There is no indication of particles that produce small amounts of linear polarization at low altitudes.     

Long slit spectropolarimetry of Jupiter and Saturn

We present ground-based limb polarization measurements of Jupiter and Saturn consisting of full disk imaging polarimetry for the wavelength 7300 Å and spatially resolved (long-slit) spectropolarimetry covering the wavelength range 5200-9350 Å.For the polar region of Jupiter we find for λ = 6000 Å a very strong radial (perpendicular to the limb) fractional polarization with a seeing corrected maximum of about +11.5% in the South and +10.0% in the North. This indicates that the polarizing haze layer is thicker at the South pole. The polar haze layers extend down to 58° in latitude. The derived polarization values are much higher than reported in previous studies because of the better spatial resolution of our data and an appropriate consideration of the atmospheric seeing. Model calculations demonstrate that the high limb polarization can be explained by strongly polarizing (p ≈ 1.0), high albedo (ω ≈ 0.98) haze particles with a scattering asymmetry parameter of g ≈ 0.6 as expected for aggregate particles of the type described by West and Smith (West, R.A., Smith, P.H. [1991]. Icarus 90, 330-333). The deduced particle parameters are distinctively different when compared to lower latitude regions.The spectropolarimetry of Jupiter shows a decrease in the polar limb polarization towards longer wavelengths and a significantly enhanced polarization in strong methane bands when compared to the adjacent continuum. This is a natural outcome for a highly polarizing haze layer above an atmosphere where multiple scatterings are suppressed in absorption bands. For lower latitudes the fractional polarization is small, negative, and it depends only little on wavelength except for the strong CH₄-band at 8870 Å.The South pole of Saturn shows a lower polarization (p ≈ 1.0-1.5%) than the poles of Jupiter. The spectropolarimetric signal for Saturn decrease rapidly with wavelength and shows no significant enhancements in the fractional polarization in the absorption bands. These properties can be explained by a vertically extended stratospheric haze region composed of small particles <100 nm as suggested previously by Karkoschka and Tomasko (Karkoschka, E., Tomasko, M. [2005]. Icarus 179, 195-221).In addition we find in the V- and R-band a previously not observed strong polarization feature (p = 1.5-2.0%) near the equator of Saturn. The origin of this polarization signal is unclear but it could be related to a seasonal effect.Finally we discuss the potential of ground-based limb polarization measurements for the investigation of the scattering particles in the atmospheres of Jupiter and Saturn.     

Vertical distribution of scattering hazes in Titan's upper atmosphere

Radial intensity scars of a Voyager 2 high phase angle image of Titan have been inverted to yield vertical extinction profiles at 1° intervals around the limb. A detached haze layer with peak particle number densities ～0.2 cm^?3 exists at all latitudes south of ～45°N, and at an altitude of 300–350 km. The optical depth 0.01 level lies at a radius of 2932 ± 5 km at the equator and at a radius of 2915 ± 10 km over the poles (altitudes of 357 ± 5 and 340 ± 10 km, respectively). In addition to the haze layer at 300–350 km, there is a small enhancement in the extinction at ～450 km which exists at all latitudes between 75°S and ～60°N.     

Two formation regions for Titan's hazes: indirect clues and possible synthesis mechanisms

It is suggested that aerosol particles forming the detached and main haze layers of Titan's atmosphere do not originate in the same atmospheric levels. Particles present above approximately 350 km could be formed of polyacetylenes synthetized in the 500-800 km altitude range through successive insertion reactions involving the C2H radical under the action of solar ultraviolet photons (Yung et al., Astrophys. J. Suppl. 55, 465, 1984). They might contain C-N oligomers in comparable amounts, as well as C-H-N oligomers synthetized at high altitude (900-1000 km) by the action of suprathermal Saturn plasma electrons. Physically, they are expected to consist of fluffy aggregates of density approximately 0.01-0.1 g cm-3. Their mass production rate is small (10(-15)-10(-14) kg m-2 s-1), that is typically 10% or less of the main haze production rate. Due to their low fall velocity, they are very sensitive to large scale horizontal motions and one substantial part of them may be swept away by meridional circulation at the detached haze level. The altitude range where these aerosols are created is well above the range proposed by Cabane et al. (Planet. Space Sci. 41, 257, 1993) for aerosols of the main haze layer, on the basis of a new fractal microphysical modeling of Titan's aggregates, that is approximately 350-400 km. A natural outcome of this apparent discrepancy is to suppose that there is a second formation region, below approximately 400 km altitude, giving rise to the main haze layer. The aim of the present paper is to review the different possible formation mechanisms of this main haze layer and assess their ability to account for the observed characteristics of the haze. Several conditions are established. The first one, called "condition A", concerns the formation altitude range imposed by fractal modeling. Possible chemical and energy sources are examined. Two additional constraints, relative to the minimum gas mass ("condition B") and input energy ("condition C") required for efficient conversion of gas into aerosols, are defined. By comparing the production rates of the haze, as derived from microphysical models, and of gaseous chemical species, as derived from photochemical models, five possible source constituents are identified: N2, CH4, C2H2, C2H6 and HCN. Polymerization of C2H2 into (C2H2)n through action of solar ultraviolet photons is shown to be rather improbable (condition A is hardly satisfied). From both our current knowledge of the gaseous phase photochemistry, through modeling and laboratory experiments, and existing models of the interaction between Saturn magnetosphere and Titan atmosphere, the formation of C-H-N polymers through action of Saturn magnetospheric energetic particles (E approximately 100 keV), is proposed as the basic polymerization mechanism in the lower formation region (conditions A,B and C are jointly satisfied).     

Photometry and polarimetry of Jupiter at large phase angles: I. Analysis of imaging data of a prominent belt and a zone from pioneer 10

Limb-darkening curves are derived from Pioneer 10 imaging data for Jupiter's STrZ (?18 to ?21° latitude) and SEBn (?5 to ?8° latitude) in red and blue light at phase angles of 12, 23, 34, 109, 120, 127, and 150°. Inhomogeneous scattering models are computed and compared with the data to constrain the vertical structure and the single-scattering phase functions of the belt and the zone in each color. The very high brightness observed at a 150° phase angle seems to require the presence of at lleast a thin layer of reasonably bright and strongly forward-scattering haze particles at pressure levelsof about 100 mbar or less above both belts and zones. Marginally successful models have been constructed in which a moderate optical thickness (τ ≥ 0.5) of haze particles was uniformly distributed in the upper 25 km-amagats of H₂. Excellent fits to the data were obtained with models having a thin (optical depths of a few tenths) haze conentraated above most of the gas. Following recent spectrospcopicanalyses, we have placed the main “cloud” layer or layers beneath about 25 km-amagats of H₂, although successful fits to our continuum data probably could be achieved also if the clouds were permitted to extend all the way up to the thin haze layer. Similarly, below the haze level our data cannot distinguish between models having two clouds separated by a clear space as suggested by R. E. Danielson and M. G. Tomasko and models with a single extensive diffuse cloud having an H₂ abundance of a few kilometer-amagats per scattering mean free path as described by W. D. Cochran. In either case, the relative brightness of the planet at each phase angle primarily serves to constrain the single-scattering phase functions of the Jovian clouds at the corresponding scattering angles. The clouds in these models are characterized by single-scattering phase functions having strong forward peaks and modest backward-scattering peaks, indicating cloud particles with dimensions larger than about 0.6 μm. In our models, a lower single-scattering albedo of the cloud particles in the belt relative to the zone accounts for the contrast between these regions. If an increased abundance of absorbing dust above uniformly bright clouds is used to explain the contrast between belts and zones at visible wavelengths, the limb darkening is steeper than that observed for the SEBn in blue light at small phase angles. The phase integral for the planet calculated for either the belt or the zone model in either color lies in the range 1.2 to 1.3. If a value of 1.25 is used with D.J. Taylor's bolometric geometric albedo of 0.28, the planet emits 2.25 or 1.7 times the energy it absorbs from the Sun if it effective temperature is 134 or 125°K, respectively—roughly as expected from current theories of the cooling of Jupiter's interior.     

Venus: Mesospheric hazes of ice,dust, and acid aerosols

We speculate on the origin and physical properties of haze in the upper atmosphere of Venus. It is argued that at least four distinct types of particles may be present. The densest and lowest haze, normally seen by spacecraft, probably consists of a submicron sulfuric acid aerosol which extends above the cloud tops (at ～70 km) up to ～80 km; this haze represents an extension of the upper cloud deck. Measurements of the temperature structure between 70 and 120 km indicate that two independent water ice layers may occasionally appear. The lower one can form between 80 and 100 km and is probably the detached haze layer seen in high-contrast limb photography. This ice layer is likely to be nucleated on sulfuric acid aerosols, and is analogous to the nacreous (stratospheric) clouds on Earth. At the Venus “mesopause” near 120 km, temperatures are frequently cold enough to allow ice nucleation on meteoric dust or ambient ions. The resulting haze (which is analogous to noctilucent clouds on Earth) is expected to be extremely tenous, and optically invisible. On both Earth and Venus, meteoric dust is present throughout the upper atmosphere and probably has similar properties.     

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.     

The structure and microphysical properties of the Venus clouds: Venera 9, 10, and 11 data

Data processing and interpretation of the nephelometer measurements made in the Venus atmosphere aboard the Venera 9, 10 and 11 landers in the sunlit hemisphere near the equator are discussed. These results were used to obtain the aerosol distribution and its microphysical properties from 62 km to the surface. The main aerosol content is found in the altitude range between 62 km (where measurements began) and 48 km, the location of the cloud region. Three prominent layers labeled as I (between 62 and 57 km), II (between 57 and 51 km) and III (between 51 and 48 km), each with different particle characteristics are discovered within the clouds. The measured light-scattering patterns can be intrepreted as having been produced by particles with effective radii from 1 to 2 μm depending on height and indices of refractivity from 1.45 in layer I to 1.42 in layer III. These values do not contradict the idea that the droplets are made of sulfuric acid. In layers II and III the particle size distribution is at least bimodal rather than uni-modal. The index of refraction is found to decrease to 1.33 in the lower part of layer II, suggesting a predominant abundance of larger particles of different chemical origin, and chlorine compounds are assumed to be relevant to this effect. In the entire heightrange of the Venera 9–11 craft descents, the clouds are rather rarefied and are characterized by a mean volume scattering coefficient σ ～ 2 × 10^?5 cm^?1 that corresponds to the mean meteorological range of visibility of about 2 km. The average mass content of condensate is estimated to be equal to 4 × 10^?9 g/cm³, and the total optical depth of clouds to τ ～ 35. Near the bottom of layer III clouds are strongly variable. In the subcloud atmosphere a haze was observed between 48 and 32 km; that haze is mainly made of submicron particles, r_eff ～ 0.1μm. The atmosphere below that is totally transparent but separate (sometimes possibly disappearing) layers may be present up to a height of 8 km above the surface. A model of this region with a very low particle density (N ? 2–3 cm^?3) strongly refractive large particles (r_eff ? 2.5 μm; 1.7 < n < 2.0) provided satisfactory agreement. The optical depth of aerosol in the atmosphere below the subcloud haze does not exceed 2.5.     

Methane on Uranus: The case for a compact CH4 cloud layer at low latitudes and a severe CH4 depletion at high-latitudes based on re-analysis of Voyager occultation measurements and STIS spectroscopy

Lindal et al. (Lindal, G.F., Lyons, J.R., Sweetnam, D.N., Eshleman, V.R., Hinson, D.P. [1987]. J. Geophys. Res. 92 (11), 14987-15001) presented a range of temperature and methane profiles for Uranus that were consistent with 1986 Voyager radio occultation measurements of refractivity versus altitude. A localized refractivity slope variation near 1.2 bars was interpreted to be the result of a condensed methane cloud layer. However, models fit to near-IR spectra found particle concentrations much deeper in the atmosphere, in the 1.5-3 bar range (Sromovsky, L.A., Irwin, P.G.J., Fry, P.M. [2006]. Icarus 182, 577-593; Sromovsky, L.A., Fry, P.M. [2010]. Icarus 210, 211-229; Irwin, P.G.J., Teanby, N.A., Davis, G.R. [2010]. Icarus 208, 913-926), and a recent analysis of STIS spectra argued for a model in which aerosol particles formed diffusely distributed hazes, with no compact condensation layer (Karkoschka, E., Tomasko, M. [2009]. Icarus 202, 287-309). To try to reconcile these results, we reanalyzed the occultation observations with the He volume mixing ratio reduced from 0.15 to 0.116, which is near the edge of the 0.033 uncertainty range given by Conrath et al. (Conrath, B., Hanel, R., Gautier, D., Marten, A., Lindal, G. [1987]. J. Geophys. Res. 92 (11), 15003-15010). This allowed us to obtain saturated mixing ratios within the putative cloud layer and to reach above-cloud and deep methane mixing ratios compatible with STIS spectral constraints. Using a 5-layer vertical aerosol model with two compact cloud layers in the 1-3 bar region, we find that the best fit pressure for the upper compact layer is virtually identical to the pressure range inferred from the occultation analysis for a methane mixing ratio near 4% at 5°S. This strongly argues that Uranus does indeed have a compact methane cloud layer. In addition, our cloud model can fit the latitudinal variations in spectra between 30°S and 20°N, using the same profiles of temperature and methane mixing ratio. But closer to the pole, the model fails to provide accurate fits without introducing an increasingly strong upper tropospheric depletion of methane at increased latitudes, in rough agreement with the trend identified by Karkoschka and Tomasko (Karkoschka, E., Tomasko, M. [2009]. Icarus 202, 287-309).     

Latitudinal variations in Titan’s methane and haze from Cassini VIMS observations

We analyze observations taken with Cassini’s Visual and Infrared Mapping Spectrometer (VIMS), to determine the current methane and haze latitudinal distribution between 60°S and 40°N. The methane variation was measured primarily from its absorption band at 0.61 μm, which is optically thin enough to be sensitive to the methane abundance at 20-50 km altitude. Haze characteristics were determined from Titan’s 0.4-1.6 μm spectra, which sample Titan’s atmosphere from the surface to 200 km altitude. Radiative transfer models based on the haze properties and methane absorption profiles at the Huygens site reproduced the observed VIMS spectra and allowed us to retrieve latitude variations in the methane abundance and haze. We find the haze variations can be reproduced by varying only the density and single scattering albedo above 80 km altitude. There is an ambiguity between methane abundance and haze optical depth, because higher haze optical depth causes shallower methane bands; thus a family of solutions is allowed by the data. We find that haze variations alone, with a constant methane abundance, can reproduce the spatial variation in the methane bands if the haze density increases by 60% between 20°S and 10°S (roughly the sub-solar latitude) and single scattering absorption increases by 20% between 60°S and 40°N. On the other hand, a higher abundance of methane between 20 and 50 km in the summer hemisphere, as much as two times that of the winter hemisphere, is also possible, if the haze variations are minimized. The range of possible methane variations between 27°S and 19°N is consistent with condensation as a result of temperature variations of 0-1.5 K at 20-30 km. Our analysis indicates that the latitudinal variations in Titan’s visible to near-IR albedo, the north/south asymmetry (NSA), result primarily from variations in the thickness of the darker haze layer, detected by Huygens DISR, above 80 km altitude. If we assume little to no latitudinal methane variations we can reproduce the NSA wavelength signatures with the derived haze characteristics. We calculate the solar heating rate as a function of latitude and derive variations of ∼10-15% near the sub-solar latitude resulting from the NSA. Most of the latitudinal variations in the heating rate stem from changes in solar zenith angle rather than compositional variations.     

Titan global climate model: A new 3-dimensional version of the IPSL Titan GCM

We have developed a new 3-dimensional climate model for Titan’s atmosphere, using the physics of the IPSL Titan 2-dimensional climate model with the current version of the LMDZ General Circulation Model dynamical core. Microphysics and photochemistry are still computed as zonal averages. This GCM covers altitudes from surface to 500 km altitude, with barotropic waves now being resolved and the diurnal cycle included. The boundary layer scheme has been changed, yielding a strong improvement in the tropospheric zonal wind profile modeled at Huygens descent position and season. The potential temperature profile is fairly consistent with Huygens observations in the lowest 10 km. The latitudinal profile of the near-surface temperature is close to observed values. The minimum of zonal wind observed by the Huygens probe just above the tropopause is also present in these simulations, and its origin is discussed by comparing solar heating and dynamical transport of energy. The stratospheric temperature and wind fields are consistent with our previous works. Compared to observations, the zonal wind peak is too weak (around 120 m/s) and too low (around 200 km). The temperature structures appear to be compressed in altitude, and depart strongly from observations in the upper stratosphere. These discrepancies are correlated, and most probably related to the altitude of the haze production. The model produces a detached haze layer located more than 150 km lower than observed by the Cassini instruments. This low production altitude is due to the current position of the GCM upper boundary. However, the temporal behaviour of the detached haze layer in the model may explain the seasonal differences observed between Cassini and Voyager 1. The waves present in the GCM are analyzed, together with their respective roles in the angular momentum budget. Though the role of the mean meridional circulation in momentum transport is similar to previous work, and the transport by barotropic waves is clearly seen in the stratosphere, a significant part of the transport at high latitudes is done all year long through low-frequency tropospheric waves that may be baroclinic waves.     

Models for Polar Haze Formation in Jupiter's Stratosphere

We present coupled chemical-microphysical models for the formation, growth, and physical properties of the jovian polar haze based on a gas-phase photochemical model for the auroral regions developed by A. S. Wong et al. (2000, Astrophys. J.534, L215-217). In this model, auroral particle precipitation provides an important energy source for enhanced decomposition of methane and production of benzene and polycyclic aromatic hydrocarbons (PAHs). We find that at high altitude, A₄ (pyrene, a hydrocarbon consisting of four fused aromatic rings) should homogeneously nucleate to form tiny primary particles. At lower altitudes, A₃ (phenanthrene) and A₂ (naphthalene) heterogeneously nucleate on the A₄ nuclei. These particles subsequently grow by additional condensation of A₂ on the nucleated particles and by coagulation and eventually sediment out to the troposphere. We run different cases of the aerosol microphysical model for different assumptions regarding the fractal dimension of aggregate particles formed by the coagulation process. If coagulation is assumed to produce spherical particles (of dimensionality 3), then their mean radius at altitudes below the 20-mbar pressure level is computed to be approximately 0.1 μm. If coagulation produces fractal aggregates of dimension 2.1, then their equivalent mean radius below the 20-mbar level is much larger, of order 0.7 μm. Aggregates with fractal dimensions between 2.1 and 3 form with equivalent mean radii between 0.1 and 0.7 μm. In every case, mean particle radius is found to decrease with increasing altitude, as expected for a system approximately in sedimentation-coagulation equilibrium. The predicted range of altitudes where aerosol formation occurs and the mean size to which particles grow are found to be generally consistent with observations. However, our calculations cannot presently account for the large amount of total aerosol loading inferred by M. G. Tomasko et al. (1986, Icarus65, 218-243). We suggest that the primarily neutral chemical pathway to heavy hydrocarbon and PAH formation proposed by Wong et al. (2000) may proceed too slowly to produce a sufficient amount of condensible material. Inclusion of ion and ion-neutral reactions in the chemical scheme could potentially lead to the prediction of higher PAH production rates, higher nucleation rates, and greater aerosol loading, producing better agreement with the observations.     

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.     

A model of Titan's haze of fractal aerosols constrained by multiple observations

We use Titan's geometric albedo to constrain the vertical distribution of the haze. Microphysical models incorporating fractal aggregates do not readily fit the methane features at 0.62 μm band and the dark 0.88 μm of the albedo spectrum simultaneously. We take advantage of this apparent discrepancy to constrain the haze vertical profile.We used the geometric albedo and several results and constraints from other works to better constrain the vertical haze extinction profile, especially in the low stratosphere. The objective of this model is to give a solution that simultaneously fits the main constraints known to apply to the haze.We find that the haze extinction increases with decreasing altitude with a scale height about equal to the atmospheric scale height down to 100 km. Below this altitude, extinction must decrease down to 30 km. This is necessary in order to have enough haze to sustain a relatively high albedo (0.076) in the dark 0.88 μm methane band and to show the 0.62 μm band in the haze continuum. We set the haze production rate around 7×10⁻¹⁴ kgm⁻² s⁻¹, and the aerosols production altitude around 400 km (or at pressure 1.5 Pa).The physical processes which generate such a profile are not clear. However, purely one-dimensional effects such as condensation, sedimentation, and rainout can be ruled out, and we believe that this relative clearing in Titan's troposphere and lower stratosphere is due to particle horizontal transport by the mean circulation.     

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.