A photochemical model of Titanʹs atmosphere and ionosphere

A global-mean model of coupled neutral and ion chemistry on Titan has been developed. Unlike the previous coupled models, the model involves ambipolar diffusion and escape of ions, hydrodynamic escape of light species, and calculates the H2 and CO densities near the surface that were assigned in some previous models. We tried to reduce the numbers of species and reactions in the model and remove all species and reactions that weakly affect the observed species. Hydrocarbon chemistry is extended to C12H10 for neutrals and C10H+11 for ions but does not include PAHs. The model involves 415 reactions of 83 neutrals and 33 ions, effects of magnetospheric electrons, protons, and cosmic rays. UV absorption by Titanʹs haze was calculated using the Huygens observations and a code for the aggregate particles. Hydrocarbon, nitrile, and ion chemistries are strongly coupled on Titan, and attempt to calculate them separately (e.g., in models of ionospheric composition) may result in significant error. The model densities of various species are typically in good agreement with the observations except vertical profiles in the stratosphere that are steeper than the CIRS limb data. (A model with eddy diffusion that facilitates fitting to the CIRS limb data is considered as well.) The CO densities are supported by the O+ flux from Saturnʹs magnetosphere. The ionosphere includes a peak at 80 km formed by the cosmic rays, steplike layers at 500–700 and 700–900 km and a peak at 1060 km (SZA = 60°). Nighttime densities of major ions agree with the INMS data. Ion chemistry dominates in the production of bicyclic aromatic hydrocarbons above 600 km. The model estimates of heavy positive and negative ions are in reasonable agreement with the Cassini results. The major haze production is in the reactions C6H + C4H2, C3N + C4H2, and condensation of hydrocarbons below 100 km. Overall, precipitation rate of the photochemical products is equal to 4–7 kg cm−2 Byr−1 (50–90 m Byr−1 while the global-mean depth of the organic sediments is ∼3 m). Escape rates of methane and hydrogen are 2.9 and 1.4 kg cm−2 Byr−1, respectively. The model does not support the low C/N ratio observed by the Huygens ACP in Titanʹs haze.