Methane, ethane, and mixed clouds in Titan's atmosphere: Properties derived from microphysical modeling

Theoretical arguments point to and recent observations confirm the existence of clouds in Titan's atmosphere, yet we possess very little data on their particle size, composition and formation mechanism. A time-dependent microphysical model is used to study the evolution of ice clouds in Titan's atmosphere. The model simulates nucleation, condensational growth, evaporation, coagulation, and transport of particles in a column of atmosphere. A variety of cloud compositions are studied, including pure ethane clouds, pure methane clouds, and mixed methane-ethane clouds (all with tholin cores). The abundance of methane cloud particles may be limited by the number of ethane coated tholin nuclei rather than the number of tholins with hydrocarbon coatings. However, even the condensation of methane onto these relatively sparse ethane/tholin cloud particles is sufficient to keep the methane close to saturation. Typical methane supersaturations are of order 0.06 on the average. For simulations which take into account recent lab measurements indicating it is relatively easy for methane to nucleate onto tholin particles without an ethane-layer present, the three types of clouds (methane, ethane, and mixed) exist simultaneously. Pure methane clouds are the most abundant cloud type and serve to lower the supersaturation to about 0.04. Cloud production does not require a continuous surface source of methane. However, clouds produced by mean motions are not the visible methane clouds seen in recent Cassini and ground-based observations. Ethane clouds in the troposphere almost instantaneously nucleate methane to form mixed clouds. However, a thin ethane ‘haze’ remains just above the tropopause for some scenarios and the mixed clouds at the tropopause remain ?50% ethane by mass. Also, evaporation of methane from the mixed cloud particles near the surface leaves a thicker layer of ethane cloud particles at ∼10 km. Nevertheless, the precipitation rate of methane to Titan's surface is between 0.001 and 0.5 cm/terrestrial-year, depending on various initial conditions such as critical saturation, size and abundance of cloud condensation nuclei, surface sources and eddy diffusion.     

Evidence for layered methane clouds in Titan’s troposphere

Layered methane clouds in Titan’s troposphere with an upper methane ice cloud, a lower liquid methane-nitrogen cloud, and a gap in between were suggested from in situ measurements and ground-based observations. Here we report laboratory investigations under conditions that mimic Titan’s troposphere providing a detailed picture of the cloud layers. A solid methane cloud with a nitrogen content of less than 14% and a liquid methane-nitrogen cloud with a nitrogen content of ∼30% form above ∼19 km and below ∼16 km altitude, respectively. Contrary to previous assertions, long-lived supercooled liquid methane-nitrogen droplets can be sustained in the region in between. The results demonstrate that a cloud gap could only form in the presence of high amounts of other traces species (ethane nuclei, tholin particles, etc.).     

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

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

The structure of Titan’s atmosphere from Cassini radio occultations

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

Probing Titan’s lower atmosphere with acousto-optic tuning

Narrow-band images of Titan were obtained in November 1999 with the NASA/GSFC- built acousto-optic imaging spectrometer (AImS) camera. This instrument utilizes a tunable filter element that was used within the 500- to 1050-nm range, coupled to a CCD camera system. The images were taken with the Mount Wilson 2.54-m (100 in.) Hooker telescope, which is equipped with a natural guide star adaptive optics system. We observed Titan at 830 and 890 nm and at a series of wavelengths across the 940-nm window in Titan’s atmosphere where the methane opacity is relatively low. We determined the absolute reflectivity (I/F) of Titan and fit the values at 940 nm to a Minnaert function at Titan’s equator and at −30° latitude (closer to the subsolar point) and obtained average values for the Minnaert limb-darkening slope, k, of 0.661 ± 0.007 and 0.775 ± 0.018, respectively. Comparison with models suggests that the equatorial value of k is consistent with rain removal of haze in the lower atmosphere. The higher value of k at −30° is consistent with the southern hemisphere being brighter than the equator. However, the fits are not unique. The data and models at 890 are consistent with no limb brightening or darkening at this wavelength either at the equator or at −30°.     

Molecular hydrogen in Titan’s atmosphere: Implications of the measured tropospheric and thermospheric mole fractions

The third most abundant species in Titan’s atmosphere is molecular hydrogen with a tropospheric/lower stratospheric mole fraction of 0.001 derived from Voyager and Cassini infrared measurements. The globally averaged thermospheric H₂ mole fraction profile from the Cassini Ion Neutral Mass Spectrometer (INMS) measurements implies a small positive gradient in the H₂ mixing ratio from the tropopause region to the lower thermosphere (∼950-1000 km), which drives a downward H₂ flux into Titan’s surface comparable to the H₂ escape flux out of the atmosphere (∼2 × 10¹⁰ cm⁻² s⁻¹ referenced to the surface) and requires larger photochemical production rates of H₂ than obtained by previous photochemical models. From detailed model calculations based on known photochemistry with eddy, molecular, and thermal diffusion, the tropospheric and thermospheric H₂ mole fractions are incompatible by a factor of ∼2. The measurements imply that the downward H₂ surface flux is in substantial excess of the speculative threshold value for methanogenic life consumption of H₂ (McKay, C.P., Smith, H.D. [2005], Icarus 178, 274-276. doi:10.1016/j.icarus.2005.05.018), but without the extreme reduction in the surface H₂ mixing ratio.     

Atomic and molecular hydrogen budget in Titan’s atmosphere

Using a one-dimensional model, we investigate the hydrogen budget and escape to space in Titan’s atmosphere. Our goal is to study in detail the distributions and fluxes of atomic and molecular hydrogen in the model, while identifying sources of qualitative and quantitative uncertainties. Our study confirms that the escape of atomic and molecular hydrogen to space is limited by the diffusion through the homopause level. The H distribution and flux inside the atmosphere are very sensitive to the eddy diffusion coefficient used above altitude 600 km. We chose a high value of this coefficient 1 × 10⁸ cm² s⁻¹ and a homopause level around altitude 900 km. We find that H flows down significantly from the production region above 500 km to the region [300-500] km, where it recombines into H₂. Production of both H and H₂ also occurs in the stratosphere, mostly from photodissociation of acetylene. The only available observational data to be compared are the escape rate of H deduced from Pioneer 11 and IUE observations of the H torus 1-3 × 10⁹ cm⁻² s⁻¹ and the latest retrieved value of the H₂ mole fraction in the stratosphere: (1.1 ± 0.1) × 10⁻³. Our results for both of these values are at least 50-100% higher, though the uncertainties within the chemical schemes and other aspects of the model are large. The chemical conversion from H to H₂ is essentially done through catalytic cycles using acetylene and diacetylene. We have studied the role of this diacetylene cycle, for which the associated reaction rates are poorly known. We find that it mostly affects C₄ species and benzene in the lower atmosphere, rather than the H profile and the hydrogen budget. We have introduced the heterogenous recombination of hydrogen on the surface of aerosol particles in the stratosphere, and this appears to be a significant process, comparable to the chemical processes. It has a major influence on the H distribution, and consequently on several other species, especially C₃H₄, C₄H₂ and C₆H₆. Therefore, this heterogenous process should be taken into account when trying to understand the stratospheric distribution of these hydrocarbons.     

The photochemical products of benzene in Titan’s upper atmosphere

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

A comprehensive numerical simulation of Io’s sublimation-driven atmosphere

Io’s sublimation-driven atmosphere is modeled using the direct simulation Monte Carlo (DSMC) method. These rarefied gas dynamics simulations improve upon earlier models by using a three-dimensional domain encompassing the entire planet computed in parallel. The effects of plasma heating, planetary rotation, inhomogeneous surface frost, molecular residence time of SO₂ on the exposed (non-volatile) rocky surface, and surface temperature distribution are investigated. Circumplanetary flow is predicted to develop from the warm dayside toward the cooler nightside. Io’s rotation leads to a highly asymmetric frost surface temperature distribution (due to the frost’s high thermal inertia) which results in circumplanetary flow that is not axi-symmetric about the subsolar point. The non-equilibrium thermal structure of the atmosphere, specifically vibrational and rotational temperatures, is also examined. Plasma heating is found to significantly inflate the atmosphere on both the dayside and nightside. The plasma energy flux causes high temperatures at high altitudes but plasma energy depletion through the dense gas column above the warmest frost permits gas temperatures cooler than the surface at low altitudes. A frost map (Douté, S., Schmitt, B., Lopes-Gautier, R., Carlson, R., Soderblom, L., Shirley, J., and the Galileo NIMS Team [2001]. Icarus 149, 107-132) is used to control the sublimated flux of SO₂ which can result in inhomogeneous column densities that vary by nearly a factor of four for the same surface temperature. A short residence time for SO₂ molecules on the “rock” component is found to smooth lateral atmospheric inhomogeneities caused by variations in the surface frost distribution, creating an atmosphere that looks nearly identical to one with uniform frost coverage. A longer residence time is found to agree better with mid-infrared observations (Spencer, J.R., Lellouch, E., Richter, M.J., López-Valverde, M.A., Jessup, K.L, Greathouse, T.K., Flaud, J. [2005]. Icarus 176, 283-304) and reproduce the observed anti-jovian/sub-jovian column density asymmetry. The computed peak dayside column density for Io assuming a surface frost temperature of 115 K agrees with those suggested by Lyman-α observations (Feaga, L.M., McGrath, M., Feldman, P.D. [2009]. Icarus 201, 570-584). On the other hand, the peak dayside column density at 120 K is a factor of five larger and is higher than the upper range of observations (Jessup, K.L., Spencer, J.R., Ballester, G.E., Howell, R.R., Roesler, F., Vigel, M., Yelle, R. [2004]. Icarus 169, 197-215; Spencer et al., 2005).     

Power spectral analysis of Jupiter’s clouds and kinetic energy from Cassini

We present suggestive evidence for an inverse energy cascade within Jupiter’s atmosphere through a calculation of the power spectrum of its kinetic energy and its cloud patterns. Using Cassini observations, we composed full-longitudinal mosaics of Jupiter’s atmosphere at several wavelengths. We also utilized image pairs derived from these observations to generate full-longitudinal maps of wind vectors and atmospheric kinetic energy within Jupiter’s troposphere. We computed power spectra of the image mosaics and kinetic energy maps using spherical harmonic analysis. Power spectra of Jupiter’s cloud patterns imaged at certain wavelengths resemble theoretical spectra of two-dimensional turbulence, with power-law slopes near −5/3 and −3 at low and high wavenumbers, respectively. The slopes of the kinetic energy power spectrum are also near −5/3 at low wavenumbers. At high wavenumbers, however, the spectral slopes are relatively flatter than the theoretical prediction of −3. In addition, the image mosaic and kinetic energy power spectra differ with respect to the location of the transition in slopes. The transition in slope is near planetary wavenumber 70 for the kinetic energy spectra, but is typically above 200 for the image mosaic spectra. Our results also show the importance of calculating spectral slopes from full 2D velocity maps rather than 1D zonal mean velocity profiles, since at large wavenumbers the spectra differ significantly, though at low wavenumbers, the 1D zonal and full 2D kinetic energy spectra are practically indistinguishable. Furthermore, the difference between the image and kinetic energy spectra suggests some caution in the interpretation of power spectrum results solely from image mosaics and its significance for the underlying dynamics. Finally, we also report prominent variations in kinetic energy within the equatorial jet stream that appear to be associated with the 5 μm hotspots. Other eddies are present within the flow collar of the Great Red Spot, suggesting caution when interpreting snapshots of the flow inside these features as representative of a time-averaged state.     

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.     

Temporal variations of Titan’s middle-atmospheric temperatures from 2004 to 2009 observed by Cassini/CIRS

We use five and one-half years of limb- and nadir-viewing temperature mapping observations by the Composite Infrared Radiometer-Spectrometer (CIRS) on the Cassini Saturn orbiter, taken between July 2004 and December 2009 (L_S from 293° to 4°; northern mid-winter to just after northern spring equinox), to monitor temperature changes in the upper stratosphere and lower mesosphere of Titan. The largest changes are in the northern (winter) polar stratopause, which has declined in temperature by over 20 K between 2005 and 2009. Throughout the rest of the mid to upper stratosphere and lower mesosphere, temperature changes are less than 5 K. In the southern hemisphere, temperatures in the middle stratosphere near 1 mbar increased by 1-2 K from 2004 through early 2007, then declined by 2-4 K throughout 2008 and 2009, with the changes being larger at more polar latitudes. Middle stratospheric temperatures at mid-northern latitudes show a small 1-2 K increase from 2005 through 2009. At north polar latitudes within the polar vortex, temperatures in the middle stratosphere show a ∼4 K increase during 2007, followed by a comparable decrease in temperatures in 2008 and into early 2009. The observed temperature changes in the north polar region are consistent with a weakening of the subsidence within the descending branch of the middle atmosphere meridional circulation.     

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.     

Small-scale composition and haze layering in Titan’s polar vortex

Fine scale layering of haze and composition in Titan’s stratosphere and mesosphere was investigated using visible/UV images from Cassini’s Imaging Science Sub-system (ISS) and IR spectra from Cassini’s Composite Infra-Red Spectrometer (CIRS). Both ISS and CIRS independently show fine layered structures in haze and composition, respectively, in the 150-450 km altitude range with a preferred vertical wavelength of around 50 km. Layers are most pronounced around the north polar winter vortex, although some weaker layers do exist at more southerly latitudes. The amplitude of composition layers in each trace gas profile is proportional to the relative enrichment of that species in the winter polar vortex compared to equatorial latitudes. As enrichment is caused by polar subsidence, this suggests a dynamical origin. We propose that the polar layers are caused by cross-latitude advection across the vortex boundary. This is analogous to processes that lead to ozone laminae formation around Earth’s polar vortices.     

Impact-induced N2 production from ammonium sulfate: Implications for the origin and evolution of N2 in Titan’s atmosphere

Chemical reactions and volatile supply through hypervelocity impacts may have played a key role for the origin and evolution of both planetary and satellite atmospheres. In this study, we evaluate the role of impact-induced N₂ production from reduced nitrogen-bearing solids proposed to be contained in Titan’s crust, ammonium sulfate ((NH₄)₂SO₄), for the replenishment of N₂ to the atmosphere in Titan’s history. To investigate the conversion of (NH₄)₂SO₄ into N₂ by hypervelocity impacts, we measured gases released from (NH₄)₂SO₄ that was exposed to hypervelocity impacts created by a laser gun. The sensitivity and accuracy of the measurements were enhanced by using an isotope labeling technique for the target. We obtained the efficiency of N₂ production from (NH₄)₂SO₄ as a function of peak shock pressure ranging from ∼8 to ∼45 GPa. Our results indicate that the initial and complete shock pressures for N₂ degassing from (NH₄)₂SO₄ are ∼10 and ∼25 GPa, respectively. These results suggest that cometary impacts on Titan (i.e., impact velocity v_i > ∼8 km/s) produce N₂ efficiently; whereas satellitesimal impacts during the accretion (i.e., v_i < 4 km/s) produce N₂ only inefficiently. Even when using the proposed small amount of (NH₄)₂SO₄ content in the crust (∼4 wt.%) (Fortes, A.D. et al., 2007. Icarus 188, 139-153), the total amount of N₂ provided through cometary impacts over 4.5 Ga reaches ∼2-6 times the present atmospheric N₂ (i.e., ∼7 × 10²⁰-2 × 10²¹ [mol]) based on the measured production efficiency and results of a hydrodynamic simulation of cometary impacts onto Titan. This implies that cometary impacts onto Titan’s crust have the potential to account for a large part of the present N₂ through the atmospheric replenishment after the accretion.     

Simulation of Io’s auroral emission: Constraints on the atmosphere in eclipse

We study the morphology of Io’s aurora by comparing simulation results of a three-dimensional (3D) two-fluid plasma model to observations by the high-resolution Long-Range Reconnaissance Imager (LORRI) on-board the New Horizons spacecraft and by the Hubble Space Telescope Advanced Camera for Surveys (HST/ACS). In 2007, Io’s auroral emission in eclipse has been observed simultaneously by LORRI and ACS and the observations revealed detailed features of the aurora, such as a huge glowing plume at the Tvashtar paterae close to the North pole. The auroral radiation is generated in Io’s atmosphere by collisions between impinging magnetospheric electrons and various neutral gas components. We calculate the interaction of the magnetospheric plasma with Io’s atmosphere-ionosphere and simulate the auroral emission. Our aurora model takes into account not only the direct influence of the atmospheric distribution on the morphology and intensity of the emission, but also the indirect influence of the atmosphere on the plasma environment and thus on the exciting electrons. We find that the observed morphology in eclipse can be explained by a smooth (non-patchy) equatorial atmosphere with a vertical column density that corresponds to ∼10% of the column density of the sunlit atmosphere. The atmosphere is asymmetric with two times higher density and extension on the downstream hemisphere. The auroral emission from the Tvashtar volcano enables us to constrain the plume gas content for the first time. According to our model, the observed intensity of the Tvashtar plume implies a mean column density of ∼5 × 10¹⁵ cm⁻² for the plume region.     

Determining a tilt in Titan’s north-south albedo asymmetry from Cassini images

Analysis of Titan’s hemispheric brightness asymmetry from mapped Cassini images reveals an axis of symmetry that is tilted with respect to the rotational axis of the solid body. Twenty images taken from 2004 through 2007 show a mean axial offset of 3.8 ± 0.9° relative to the solid body’s pole, directed 79 ± 24° to the west of the sub-solar longitude. These values are consistent with recent measurements of an implied atmospheric spin axis determined from isothermal mapping by [Achterberg, R.K., Conrath, B.J., Gierasch, P.J., Flasar, F.M., Nixon, C.A., 2008. Icarus 197, 549-555].     

Titan's hydrodynamically escaping atmosphere

The upper atmosphere of Titan is currently losing mass at a rate , by hydrodynamic escape as a high density, slow outward expansion driven principally by solar UV heating by CH₄ absorption. The hydrodynamic mass loss is essentially CH₄ and H₂ escape. Their combined escape rates are restricted by power limitations from attaining their limiting rates (and limiting fluxes). Hence they must exhibit gravitational diffusive separation in the upper atmosphere with increasing mixing ratios to eventually become major constituents in the exosphere. A theoretical model with solar EUV heating by N₂ absorption balanced by HCN rotational line cooling in the upper thermosphere yields densities and temperatures consistent with the Huygens Atmospheric Science Investigation (HASI) data [Fulchignoni, M., and 42 colleagues, 2005. Nature 438, 785-791], with a peak temperature of ∼185-190 K between 3500-3550 km. This model implies hydrodynamic escape rates of and , or some other combination with a higher H₂ escape flux, much closer to its limiting value, at the expense of a slightly lower CH₄ escape rate. Nonthermal escape processes are not required to account for the loss rates of CH₄ and H₂, inferred by the Cassini Ion Neutral Mass Spectrometer (INMS) measurements [Yelle, R.V., Borggren, N., de la Haye, V., Kasprzak, W.T., Niemann, H.B., Müller-Wodarg, I., Waite Jr., J.H., 2006. Icarus 182, 567-576].     

Convective cloud heights as a diagnostic for methane environment on Titan

The appearance of convective clouds in Titan’s troposphere has been documented from ground-based observation for more than a decade. Cloud tops have been reported between 14 and 25 km. Higher resolution Cassini data have shown smaller portions of the cloud system can reach up to 42 km. We use the Titan Regional Atmospheric Modeling System (TRAMS) to explore environments which allow convective clouds to reach the tropopause. In general, cloud tops remain below 30 km, but for environments where the surface humidity of methane is greater than 50%, a small portion at the center of the cloud rises briefly to higher altitudes; for ?65% humidity, the cloud top reaches nearly to the tropopause (∼40 km). A number of other parameters also have noticeable affects on cloud top such as nucleation critical saturation, haze abundance, and collisional growth of cloud particles.