Spatial and temporal variations in Titan's surface temperatures from Cassini CIRS observations

We report a wide-ranging study of Titan's surface temperatures by analysis of the Moon's outgoing radiance through a spectral window in the thermal infrared at 19 μm (530 cm^?1) characterized by lower atmospheric opacity. We begin by modeling Cassini Composite Infrared Spectrometer (CIRS) far infrared spectra collected in the period 2004–2010, using a radiative transfer forward model combined with a non-linear optimal estimation inversion method. At low-latitudes, we agree with the HASI near-surface temperature of about 94 K at 10°S (Fulchignoni et al., 2005). We find a systematic decrease from the equator toward the poles, hemispherically asymmetric, of ～1 K at 60° south and ～3 K at 60° north, in general agreement with a previous analysis of CIRS data (Jennings et al., 2009), and with Voyager results from the previous northern winter. Subdividing the available database, corresponding to about one Titan season, into 3 consecutive periods, small seasonal changes of up to 2 K at 60°N became noticeable in the results. In addition, clear evidence of diurnal variations of the surface temperatures near the equator are observed for the first time: we find a trend of slowly increasing temperature from the morning to the early afternoon and a faster decrease during the night. The diurnal change is ～1.5 K, in agreement with model predictions for a surface with a thermal inertia between 300 and 600 J m^?2 s^?1/2 K^?1. These results provide important constraints on coupled surface–atmosphere models of Titan's meteorology and atmospheric dynamic.     

Radiative transfer analyses of Titan’s tropical atmosphere

Titan’s optical and near-IR spectra result primarily from the scattering of sunlight by haze and its absorption by methane. With a column abundance of 92 km amagat (11 times that of Earth), Titan’s atmosphere is optically thick and only ～10% of the incident solar radiation reaches the surface, compared to 57% on Earth. Such a formidable atmosphere obstructs investigations of the moon’s lower troposphere and surface, which are highly sensitive to the radiative transfer treatment of methane absorption and haze scattering. The absorption and scattering characteristics of Titan’s atmosphere have been constrained by the Huygens Probe Descent Imager/Spectral Radiometer (DISR) experiment for conditions at the probe landing site (Tomasko, M.G., Bézard, B., Doose, L., Engel, S., Karkoschka, E. [2008a]. Planet. Space Sci. 56, 624–247; Tomasko, M.G. et al. [2008b]. Planet. Space Sci. 56, 669–707). Cassini’s Visual and Infrared Mapping Spectrometer (VIMS) data indicate that the rest of the atmosphere (except for the polar regions) can be understood with small perturbations in the high haze structure determined at the landing site (Penteado, P.F., Griffith, C.A., Tomasko, M.G., Engel, S., See, C., Doose, L., Baines, K.H., Brown, R.H., Buratti, B.J., Clark, R., Nicholson, P., Sotin, C. [2010]. Icarus 206, 352–365). However the in situ measurements were analyzed with a doubling and adding radiative transfer calculation that differs considerably from the discrete ordinates codes used to interpret remote data from Cassini and ground-based measurements. In addition, the calibration of the VIMS data with respect to the DISR data has not yet been tested. Here, VIMS data of the probe landing site are analyzed with the DISR radiative transfer method and the faster discrete ordinates radiative transfer calculation; both models are consistent (to within 0.3%) and reproduce the scattering and absorption characteristics derived from in situ measurements. Constraints on the atmospheric opacity at wavelengths outside those measured by DISR, that is from 1.6 to 5.0 μm, are derived using clouds as diffuse reflectors in order to derive Titan’s surface albedo to within a few percent error and cloud altitudes to within 5 km error. VIMS spectra of Titan at 2.6–3.2 μm indicate not only spectral features due to CH₄ and CH₃D (Rannou, P., Cours, T., Le Mouélic, S., Rodriguez, S., Sotin, C., Drossart, P., Brown, R. [2010]. Icarus 208, 850–867), but also a fairly uniform absorption of unknown source, equivalent to the effects of a darkening of the haze to a single scattering albedo of 0.63 ± 0.05. Titan’s 4.8 μm spectrum point to a haze optical depth of 0.2 at that wavelength. Cloud spectra at 2 μm indicate that the far wings of the Voigt profile extend 460 cm^?1 from methane line centers. This paper releases the doubling and adding radiative transfer code developed by the DISR team, so that future studies of Titan’s atmosphere and surface are consistent with the findings by the Huygens Probe. We derive the surface albedo at eight spectral regions of the 8 × 12 km² area surrounding the Huygens landing site. Within the 0.4–1.6 μm spectral region our surface albedos match DISR measurements, indicating that DISR and VIMS measurements are consistently calibrated. These values together with albedos at longer 1.9–5.0 μm wavelengths, not sampled by DISR, resemble a dark version of the spectrum of Ganymede’s icy leading hemisphere. The eight surface albedos of the landing site are consistent with, but not deterministic of, exposed water ice with dark impurities.     

Kinetics of methane clathrate formation and dissociation under Mars relevant conditions

Spectral observations have detected methane within the martian atmosphere (Formisano, V., Atreya, S., Encrenaz, T., Ignatiev, N., Giuranna, M. [2004]. Science 306, 1758–1761; Mumma, M.J. et al. [2009]. Science 323, 1041–1045), however, the origin of the methane has not been determined. Methane clathrate (also referred to as methane hydrate) has been suggested as a potential subsurface reservoir, storing and releasing biologic and/or abiogenic methane. In this study, rates of methane hydrate formation and dissociation were measured experimentally at 234–264 K and 1.4–4.7 MPa to test the clathrate reservoir hypothesis. Initial formation rates range from 4.3 × 10^?6 to 8.1 × 10^?5 mol m^?2 s^?1. Results show decreasing rates of formation over time in individual experiments, indicating initial rapid clathration, followed by diffusion-limited transport of methane into the ice through the previously formed hydrate. These experiments indicate increased pressure results in increased formation rates, likely the result of higher concentration gradients, enhancing the methane diffusion flux into the solid phase. Experiments conducted at elevated temperatures produced faster initial rates of formation, resulting from increased kinetic energy of methane molecules and/or thickening of the Quasi-Liquid Layer. Based on this temperature dependence, the activation energy for methane hydrate formation from ice was determined to be 35.9 kJ/mol. Hydrate dissociation experiments initiated by depressurization or warming at conditions between 222 K and 265 K and 0.1–2.0 MPa were conducted following each formation experiment, yielding methane hydrate dissociation rates from 3.01 × 10^?6 to 9.92 × 10^?5 mol m^?2 s^?1. While both hydrate dissociation and formation showed decreasing instantaneous rates over the course of each experiment, the transition between the initial rate of dissociation and the interpreted diffusion-limited period of continued dissociation was more abrupt than that observed in formation experiments, supporting an ice shielding effect. The initial concentration of methane in the solid phase had a significant effect on hydrate dissociation rates. Higher methane concentrations in the solid phase produce faster initial rates, likely due to increased concentration gradients, thus increasing the diffusion component of dissociation. Increased temperatures also produced faster dissociation rates, yielding an activation energy for dissociation of 32.7 kJ/mol. The rates determined within this study suggest that small near-surface methane hydrate reservoirs are a feasible source for recent methane plumes detected on Mars. Rates of methane release from gas hydrates also indicate that gas hydrate dissociation may have played a role in forming ancient chaos terrain and associated outflow channels.     

Optical constants from 370 nm to 900 nm of Titan tholins produced in a low pressure RF plasma discharge

Determining the optical constants of Titan aerosol analogues, or tholins, has been a major concern for the last three decades because they are essential to constrain the numerical models used to analyze Titan’s observational data (albedo, radiative transfer, haze vertical profile, surface contribution, etc.). Here we present the optical constant characterization of tholins produced with an RF plasma discharge in a (95%N₂–5%CH₄) gas mixture simulating Titan’s main atmospheric composition, and deposited as a thin film on an Al–SiO₂ substrate. The real and imaginary parts, n and k, of the tholin complex refractive index have been determined from 370 nm to 900 nm wavelength using spectroscopic ellipsometry. The values of n decrease from n = 1.64 (at 370 nm) to n = 1.57 (at 900 nm) as well as the values of k which feature two behaviors: an exponential decay from 370 nm to 500 nm, with k = 12.4 × e^?0.018λ (where λ is expressed in nm), followed by a plateau, with k = (1.8 ± 0.2) × 10^?3. The trends observed for the PAMPRE tholins optical constants are compared to those determined for other Titan tholins, as well as to the optical constants of Titan’s aerosols retrieved from observational data.     

Influence of methane concentration on the optical indices of Titan’s aerosols analogues

This work deals with the optical constant characterization of Titan aerosol analogues or “tholins” produced with the PAMPRE experimental setup and deposited as thin films onto a silicon substrate. Tholins were produced in different N₂–CH₄ gaseous mixtures to study the effect of the initial methane concentration on their optical constants. The real (n) and imaginary (k) parts of the complex refractive index were determined using the spectroscopic ellipsometry technique in the 370–1000 nm wavelength range. We found that optical constants depend strongly on the methane concentrations of the gas phase in which tholins are produced: imaginary optical index (k) decreases with initial CH₄ concentration from 2.3 × 10^?2 down to 2.7 × 10^?3 at 1000 nm wavelength, while the real optical index (n) increases from 1.48 up to 1.58 at 1000 nm wavelength. The larger absorption in the visible range of tholins produced at lower methane percentage is explained by an increase of the secondary and primary amines signature in the mid-IR absorption. Comparison with results of other tholins and data from Titan observations are presented. We found an agreement between our values obtained with 10% methane concentration, and Imanaka et al. (Imanaka, H., Khare, B.N., Elsila, J.E., Bakes, E.L.O., McKay, C.P., Cruikshank, D.P., Sugita, S., Matsui, T., Zare, R.N. [2004]. Icarus, 168, 344–366) values, in spite of the difference in the analytical method. This confirms a reliability of the optical properties of tholins prepared with various setups but with similar plasma conditions. Our comparison with Titan’s observations also raises a possible inconsistency between the mid-IR aerosol signature by VIMS and CIRS Cassini instruments and the visible Huygens-DISR derived data. The mid-IR VIMS and CIRS signatures are in agreement with an aerosol dominated by an aliphatic carbon content, whereas the important visible absorption derived from the DISR measurement seems to be incompatible with such an important aliphatic content, but more compatible with an amine-rich aerosol.     

Line broadening in the neutral and ionized mercury spectra

The neutral, singly, doubly and triply ionized mercury (Hg I–IV, respectively) spectral line shapes and line center positions have been investigated in the laboratory helium plasma at electron densities ranging between 9.3 × 10²² m^?3 and 1.93 × 10²³ m^?3 and electron temperatures around 19,500 K, both interesting for astrophysics. The mercury (natural isotope composition) atoms were sputtered from the cylindrical amalgamated gold plates located in the homogenous part of the pulsed helium discharge operating at a pressure of 665 Pa in a flowing regime. The mercury spectral line profiles were recorded using the McPherson model 209 spectrograph and the Andor ICCD camera as the detection system. This research presents Stark broadening parameters, the width (W) and the shift (d), of one Hg I, 19 Hg II, 6 Hg III and 4 Hg IV lines, not investigated so far. Our experimental W values were compared with the data calculated applying various approaches. The shape and intensity of astrophysically important 398.4 nm Hg II spectral line was discussed taking into account the isotope shift, hyperfine structure and Penning effects. At the mentioned plasma parameters the Stark broadening is found to be a main line broadening mechanism of the lines (λ > 200 nm) in the Hg I–IV spectra.     

Constraints on Saturn’s tropospheric general circulation from Cassini ISS images

An automated cloud tracking algorithm is applied to Cassini Imaging Science Subsystem high-resolution apoapsis images of Saturn from 2005 and 2007 and moderate resolution images from 2011 and 2012 to define the near-global distribution of zonal winds and eddy momentum fluxes at the middle troposphere cloud level and in the upper troposphere haze. Improvements in the tracking algorithm combined with the greater feature contrast in the northern hemisphere during the approach to spring equinox allow for better rejection of erroneous wind vectors, a more objective assessment at any latitude of the quality of the mean zonal wind, and a population of winds comparable in size to that available for the much higher contrast atmosphere of Jupiter. Zonal winds at cloud level changed little between 2005 and 2007 at all latitudes sampled. Upper troposphere zonal winds derived from methane band images are ～10 m s^?1 weaker than cloud level winds in the cores of eastward jets and ～5 m s^?1 stronger on either side of the jet core, i.e., eastward jets appear to broaden with increasing altitude. In westward jet regions winds are approximately the same at both altitudes. Lateral eddy momentum fluxes are directed into eastward jet cores, including the strong equatorial jet, and away from westward jet cores and weaken with increasing altitude on the flanks of the eastward jets, consistent with the upward broadening of these jets. The conversion rate of eddy to mean zonal kinetic energy at the visible cloud level is larger in eastward jet regions (5.2 × 10^?5 m² s^?3) and smaller in westward jet regions (1.6 × 10^?5 m² s^?3) than the global mean value (4.1 × 10^?5 m² s^?3). Overall the results are consistent with theories that suggest that the jets and the overturning meridional circulation at cloud level on Saturn are maintained at least in part by eddies due to instabilities of the large-scale flow near and/or below the cloud level.     

Explaining the redox imbalance between the H and O escape fluxes at Mars by the oxidation of methane

From a comparison between the different observations of Martian methane existing today, including the new TES methane maps (Fonti and Marzo, 2010), we show that all sets of data are globally consistent with each other, and that a well definite seasonal cycle of methane has been at work for at least 10 yr. With a simple model of the balance between the loss fluxes of H and O, using up-to-date values of the escape fluxes, we show that the long-standing enigma of the imbalance between H and O escape fluxes may be solved by assuming that the missing sink of oxygen is the oxidation of methane. If no H₂ is released together with CH₄, a good agreement is found between the present CH₄ flux and the value imposed by the balance between H and O escape fluxes, an average over the last ≈10³ yr. If H₂ is released together with CH₄, as expected if CH₄ originates in serpentinization, the average level of CH₄ during the last 10³ yr should have been at least ten times lower than the present one. The lack of present H₂ release could suggest a long-term storage of methane in the subsurface under the form of clathrates, whereas H₂ has been lost to the atmosphere shortly after being produced. We suggest that the thin layer of CO₂ ice covering the permanent southern polar cap could result from the release of methane since the end of the last obliquity transition (time scale: 1 Myr), at an average rate of 0.1 Mt yr^?1, consistent with the values derived from: (i) the present observations of methane (time scale: 10 yr), (ii) the estimate from the observed imbalance between the H and O escape fluxes (time scale: 1 kyr). If so, the present release of methane from subsurface clathrates would have acted at a similar rate since at least 3 Myr.     

Titan’s internal structure and the evolutionary consequences

Titan’s moment of inertia (MoI), estimated from the quadrupole gravity field measured by the Cassini spacecraft, is 0.342, which has been interpreted as evidence of a partially differentiated internal mass distribution. It is shown here that the observed MoI is equally consistent with a fully differentiated internal structure comprising a shell of water ice overlying a low-density silicate core; depending on the chemistry of Titan’s subsurface ocean, the core radius is between 1980 and 2120 km, and its uncompressed density is 2570–2460 kg m^?3, suggestive of a hydrated CI carbonaceous chondrite mineralogy. Both the partially differentiated and fully differentiated hydrated core models constrain the deep interior to be several hundred degrees cooler than previously thought. I propose that Titan has a warm wet core below, or buffered at, the high-pressure dehydration temperature of its hydrous constituents, and that many of the gases evolved by thermochemical and radiogenic processes in the core (such as CH₄ and ⁴⁰Ar, respectively) diffuse into the icy mantle to form clathrate hydrates, which in turn may provide a comparatively impermeable barrier to further diffusion. Hence we should not necessarily expect to see a strong isotopic signature of serpentinization in Titan’s atmosphere.     

Regional differences in gully occurrence on Mars: A comparison between the Hale and Bond craters

The observation of gullies on Mars raised questions about the presence of liquid water in the recent past. In some regions like Hale and Bond crater, gullies occur in one crater (Hale) but do not in another crater nearby (Bond). These regional differences have been interpreted as an argument for a formation of the gullies related to groundwater. The formation of gullies on Earth depends on rainfall and/or melting of snow as well as on several parameters such as the presence of steep slopes and sufficient amounts of fines and debris. We investigated the Hale/Bond region for differences in crater wall morphology and texture, slopes, and thermal properties to determine whether the gully formation is dependent on factors such as steep slope angles and availability of fine-grained material. Morphologically there exist two kinds of gullies in the Hale crater: Gullies on the south- and east-facing crater slopes have a pristine appearance with deep channels eroded into the talus material and well-preserved aprons. Gully-like features on the north- and west-facing slopes are degraded and superposed by craters, indicating that they are old in comparison to the pristine ones. However, their formation process is unclear and might be due to debris flows, surface runoff or dry mass wasting processes or a combination of these processes. The crater walls of Bond do not show gullies. Their morphology is most likely consistent with a degraded mantle deposit. Slope measurements reveal that the gullies in Hale crater occur on slopes between ～20° and ～30° in contrast to the slopes without gullies in Bond that are between ～10° and ～20° steep. Mean thermal inertia values on slopes with younger gullies are ～175 J m^?2 K^?1 s^?1/2 corresponding to higher amounts of fine-grained material. At slopes with older gully-like features mean thermal inertia values are ～315 J m^?2 K^?1 s^?1/2 corresponding to higher amounts of bedrock or possibly indurated grain sizes. Mean thermal inertia values of the Bond crater walls are ～230 J m^?2 K^?1 s^?1/2 indicating more consolidated terrain possibly due to the cementation of the dissected mantle material. From our investigation we conclude that the occurrence of gullies in the Hale/Bond region most likely depends on the distribution of unconsolidated material and steep slopes. The regional and local gully distribution on Mars likely varies due to differences in topography and surface material properties. Their proposed clustered distribution on Mars is not an argument for a groundwater formation mechanism of the gullies.     

Charting thermal emission variability at Pele,Janus Patera and Kanehekili Fluctus with the Galileo NIMS Io Thermal Emission Database (NITED)

Using the NIMS Io Thermal Emission Database (NITED), a collection of over 1000 measurements of radiant flux from Io’s volcanoes (Davies, A.G. et al. [2012]. Geophys. Res. Lett. 39, L01201. doi:10.1029/2011GL049999), we have examined the variability of thermal emission from three of Io’s volcanoes: Pele, Janus Patera and Kanehekili Fluctus. At Pele, the 5-μm thermal emission as derived from 28 night time observations is remarkably steady at 37 ± 10 GW μm^?1, re-affirming previous analyses that suggested that Pele an active, rapidly overturning silicate lava lake. Janus Patera also exhibits relatively steady 5-μm thermal emission (≈20 ± 3 GW μm^?1) in the four observations where Janus is resolved from nearby Kanehekili Fluctus. Janus Patera might contain a Pele-like lava lake with an effusion rate (Q_F) of ≈40–70 m³ s^?1. It should be a prime target for a future mission to Io in order to obtain data to determine lava eruption temperature. Kanehekili Fluctus has a thermal emission spectrum that is indicative of the emplacement of lava flows with insulated crusts. Effusion rate at Kanehekili Fluctus dropped by an order of magnitude from ≈95 m³ s^?1 in mid-1997 to ≈4 m³ s^?1 in late 2001.     

Restitution coefficients and sticking velocities of a chondrule analogue colliding on a porous silica layer at impact velocities between 0.1 and 80 m s−1

To study the accretional growth of rimmed chondrules and their agglomerates in the solar nebula, we measured the restitution coefficients, ε, and the sticking velocities to a porous silica layer, v_c, by impacting the silica layer with a glass ball at velocities from 0.1 to 80 m s^?1. We used a porous silica layer covering a basalt block with thicknesses ranging from 1/5 of the glass ball radius to equal to the glass ball radius as a rimmed chondrule analogue, and the porosity of the silica layer was set to be 70%, 80%, 85%, and 90%. Collisional experiments were conducted by means of the free fall method or by the use of a spring gun or a gas gun, allowing us to vary the impact velocity. We used a laser displacement meter to estimate the impact and rebound velocities as well as the acceleration during the collision at impact velocities below 1 m s^?1. As a result, the sticking velocity, v_c, of 90%- and 85%- porosity layers with a thickness equal to 1/2 of the glass ball diameter was 0.44 and 2.4 m s^?1, respectively. On the other hand, we found a distinct barrier to sticking for smaller-porosity layers: the silicate layer with a porosity smaller than 80% never exhibited sticking at any impact velocity below 1 m s^?1. Instead, we observed a rebound effect with restitution coefficients larger than 0.2. In the case of a silica layer with a porosity smaller than 80%, we observed the sub-sticking condition defined by ε < 0.1 at velocities extending from 5 m s^?1 to 70 m s^?1.     

The growth of wind-waves in Titan’s hydrocarbon seas

Ocean wave growth on Titan is considered. The classic Sverdrup–Munk theory for terrestrial wave growth is applied to Titan, and is compared with a simple energy balance model that exposes the effect of Titan’s environmental parameters (air density, gravity, and fluid density). These approaches are compared with the only previously-published (semi-empirical) model (Ghafoor, N.A.-L., Zarnecki, J.C., Challenor, P., Srokosz, M.A. [2000] J. Geophys. Res. 105, 12,077–12,091, hereafter G2k), and allow the impact of various parameters such as atmospheric density to be transparently explored.Our model, like G2k, suggests fully-developed significant wave heights on Titan H_s = 0.2 U², where U is the windspeed (SI units): in dimensionless terms this is rather close to H_s = 0.2 U²/g, a rule of thumb previously noted for terrestrial waves (we find various datasets where the prefactor varies by ～2). It is noted that liquid and air densities affect the growth rate of waves, but not their fully-developed height: for 1 m/s winds wave amplitude reaches 0.15 m (75% of fully-developed) with a fetch of only 1 km, rather faster than predicted by G2k. Liquid viscosity has no major effect on gravity wave growth, but does influence the threshold windspeed at which gravity–capillary waves form in the first place.The model is used to develop predicted ranges for wave height to guide the design of the Titan Mare Explorer (TiME), a proposed Discovery-class mission to float a capsule on Ligeia Mare in 2023. For the expected maximum 1 m/s winds, a significant wave height of 0.2 m and wavelength of ～4 m can be expected. Assuming that wave heights follow Rayleigh statistics as they do on Earth, then given the wave period of ～4 s, individual waves of ～0.6 m might be encountered over a 3 month period.For predicted Titan winds at Kraken Mare, significant wave heights may reach ～0.6 m in the peak of summer but do not exceed the tidal amplitude at its northern end, consistent with the area around Mayda Insula being a tidal flat, while elsewhere on Kraken and Ligeia and at Ontario Lacus, shorelines may be wave- or tidally-dominated, depending on the specific location.     

Search for methane and upper limits to ethane and SO2 on Mars

The IRTF/CSHELL observations in February 2006 at L_S = 10° and 63–93°W show ～10 ppb of methane at 45°S to 7°N and ～3 ppb outside this region that covers the deepest canyon Valles Marineris. Observations in December 2009 at L_S = 20° and 0–30°W included spectra of the Moon at a similar airmass as a telluric calibrator. A technique for extraction of the martian methane line from a combination of the Mars and Moon spectra has been developed. The observations reveal no methane with an upper limit of 8 ppb. The results of both sessions agree with the observations by Mumma et al. (Mumma, M.J. et al. [2009]. Science 323, 1041–1045) at the same season in February 2006 and are smaller than those in the PFS and TES maps. Production and removal of the biological methane on Mars do not significantly change the redox state of the atmosphere and the balance of hydrogen. A search for ethane at 2977 cm^?1 results in an upper limit of 0.2 ppb. However, this limit does not help to establish the origin of methane on Mars. Reanalysis of our search for SO₂ using TEXES confirms the recently established upper limit of 0.3 ppb. Along with the lack of hot spots and gas vents with endogenic heat sources in the THEMIS observations, the very low upper limit to SO₂ on Mars does not favor geological methane that is less abundant than SO₂ in the outgassing from the terrestrial planets.     

Titan’s aerosol and stratospheric ice opacities between 18 and 500 μm: Vertical and spectral characteristics from Cassini CIRS

Vertical distributions and spectral characteristics of Titan’s photochemical aerosol and stratospheric ices are determined between 20 and 560 cm^?1 (500–18 μm) from the Cassini Composite Infrared Spectrometer (CIRS). Results are obtained for latitudes of 15°N, 15°S, and 58°S, where accurate temperature profiles can be independently determined.In addition, estimates of aerosol and ice abundances at 62°N relative to those at 15°S are derived. Aerosol abundances are comparable at the two latitudes, but stratospheric ices are ～3 times more abundant at 62°N than at 15°S. Generally, nitrile ice clouds (probably HCN and HC₃N), as inferred from a composite emission feature at ～160 cm^?1, appear to be located over a narrow altitude range in the stratosphere centered at ～90 km. Although most abundant at high northern latitudes, these nitrile ice clouds extend down through low latitudes and into mid southern latitudes, at least as far as 58°S.There is some evidence of a second ice cloud layer at ～60 km altitude at 58°S associated with an emission feature at ～80 cm^?1. We speculate that the identify of this cloud may be due to C₂H₆ ice, which in the vapor phase is the most abundant hydrocarbon (next to CH₄) in the stratosphere of Titan.Unlike the highly restricted range of altitudes (50–100 km) associated with organic condensate clouds, Titan’s photochemical aerosol appears to be well-mixed from the surface to the top of the stratosphere near an altitude of 300 km, and the spectral shape does not appear to change between 15°N and 58°S latitude. The ratio of aerosol-to-gas scale heights range from 1.3–2.4 at about 160 km to 1.1–1.4 at 300 km, although there is considerable variability with latitude. The aerosol exhibits a very broad emission feature peaking at ～140 cm^?1. Due to its extreme breadth and low wavenumber, we speculate that this feature may be caused by low-energy vibrations of two-dimensional lattice structures of large molecules. Examples of such molecules include polycyclic aromatic hydrocarbons (PAHs) and nitrogenated aromatics.Finally, volume extinction coefficients Nχ_E derived from 15°S CIRS data at a wavelength of λ = 62.5 μm are compared with those derived from the 10°S Huygens Descent Imager/Spectral Radiometer (DISR) data at 1.583 μm. This comparison yields volume extinction coefficient ratios Nχ_E(1.583 μm)/Nχ_E(62.5 μm) of roughly 70 and 20, respectively, for Titan’s aerosol and stratospheric ices. The inferred particle cross-section ratios χ_E(1.583 μm)/χ_E(62.5 μm) appear to be consistent with sub-micron size aerosol particles, and effective radii of only a few microns for stratospheric ice cloud particles.     

The heat flow during the formation of ribbon terrains on Venus

Ribbons are regularly spaced, between 2 and 6 km, troughs that exist on venusian tesserae, which are mainly located in, and characterize to, venusian crustal plateaus. Independent of the geological or temporal relations with other features, regularly and similarly spaced ribbons on several tesserae strongly suggest a thermal control on the thickness of the deformed layer. This can be used to constraint the heat flow at the time of ribbon formation, which holds important implications for the viability of the hypotheses that address the origin and evolution of crustal plateaus. For a brittle–ductile transition ∼1–3 km deep (as proposed from ribbon spacing), realistic strain rates, and a present-day surface temperature of 740 K, the implied heat flow is very high, 130–780 mW m⁻². If Venus has experienced higher surface temperatures due to climate forcing by massive volcanism, then the heat flow could be greatly reduced. For surface temperatures of 850–900 K the heat flow is 190–560, 60–230 and 20–130 mW m⁻² for brittle–ductile transition depths of 1, 2 and 3 km, respectively. Heat flow values around 80–100 mW m⁻² are reasonable for venusian hotspots, based on terrestrial analogs, but hardly consistent with coldspot settings. High surface temperatures are also required to maintain the crustal solidus deeper than a few kilometers during the formation of ribbon terrains. For the obtained heat flows, a solidus deeper than ∼30 km (the likely mean value for the crustal thickness) is difficult to achieve. This suggests that a substantial proportion of the crust beneath crustal plateaus was emplaced subsequently to the time when ribbon terrains were formed. Alternatively, at that time a magma reservoir inside the crust could have existed.     

Assessing the long-term variability of Venus winds at cloud level from VIRTIS–Venus Express

The Venus Express (VEX) mission has been in orbit to Venus for more than 4 years now. The Visible and Infrared Thermal Imaging Spectrometer (VIRTIS) instrument onboard VEX observes Venus in two channels (visible and infrared) obtaining spectra and multi-wavelength images of the planet that can be used to sample the atmosphere at different altitudes. Day-side images in the ultraviolet range (380 nm) are used to study the dynamics of the upper cloud at 66–72 km while night-side images in the near infrared (1.74 μm) map the opacity of the lower cloud deck at 44–48 km. Here we present a long-term analysis of the global atmospheric dynamics at these levels using a large selection of orbits from the VIRTIS-M dataset covering 860 Earth days that extends our previous work (Sánchez-Lavega, A. et al. [2008]. Geophys. Res. Lett. 35, L13204) and allows studying the variability of the global circulation at the two altitude levels. The atmospheric superrotation is evident with equatorial to mid-latitudes westward velocities of 100 and 60 m s^?1 in the upper and lower cloud layers. These zonal velocities are almost constant in latitude from the equator to 50°S. From 50°S to 90°S the zonal winds at both cloud layers decrease steadily to zero at the pole. Individual cloud tracked winds have errors of 3–10 m s^?1 with a mean of 5 m s^?1 and the standard deviations for a given latitude of our zonal and meridional winds are 9 m s^?1. The zonal winds in the upper cloud change with the local time in a way that can be interpreted in terms of a solar tide. The zonal winds in the lower cloud are stable at mid-latitudes to the tropics and present variability at subpolar latitudes apparently linked to the activity of the South polar vortex. While the upper cloud presents a net meridional motion consistent with the upper branch of a Hadley cell with peak velocity v = 10 m s^?1 at 50°S, the lower cloud meridional motions are less organized with some cloud features moving with intense northwards and southwards motions up to v = ±15 m s^?1 but, on average, with almost null global meridional motions at all latitudes. We also examine the long-term behavior of the winds at these two vertical layers by comparing our extended wind tracked data with results from previous missions.     

High resolution investigation of the 7 μm region of the ethane spectrum

Building upon previous studies, we re-investigated the ethane spectrum between 1330 and 1610 cm^?1 by combining unapodized spectra obtained at room temperature with a Bruker Fourier transform spectrometer (FTS) in Brussels and at 131 K with a Bruker FTS in Pasadena. The maximum optical path differences (MOPD) of the two datasets were 450 and 323.7 cm, corresponding to spectral resolutions of 0.0020 and 0.0028 cm^?1, respectively. Of the 15,000 lines observed, over 4592 transitions were assigned to the ν₆ (at 1379 cm^?1), ν₈ (at 1472 cm^?1), ν₄+ν₁₂ (at 1481 cm^?1) and 2ν₄+ν₉ (at 1388 cm^?1) bands, and another 1044 transitions were located for the ν₄+ν₈?ν₄ hot band (at 1472 cm^?1). Our new analysis included an improved implementation of the Hamiltonian calculation needed to interpret the complex spectral structures caused by numerous interactions affecting these four modes of vibration. From these results, we created the first line-by-line database containing the molecular parameters for over 20,000 ¹²C₂H₆ transitions at 7 μm.     

Io’s atmosphere: Constraints on sublimation support from density variations on seasonal timescales using NASA IRTF/TEXES observations from 2001 to 2010

We present an analysis of 19 μm spectra of Io’s SO₂ atmosphere from the TEXES mid-infrared high spectral resolution spectrograph on NASA’s Infrared Telescope Facility, incorporating new data taken between January 2005 and June 2010 and a re-analysis of earlier data taken from November 2001 to January 2004. This is the longest set of contiguous observations of Io’s atmosphere using the same instrument and technique thus far. We have fitted all 16 detected blended absorption lines of the ν₂ SO₂ vibrational band to retrieve the subsolar values of SO₂ column abundance and the gas kinetic temperature. By incorporating an existing model of Io’s surface temperatures and atmosphere, we retrieve sub-solar column densities from the disk-integrated data. Spectra from all years are best fit by atmospheric temperatures <150 K. Best-fit gas kinetic temperatures on the anti-Jupiter hemisphere, where SO₂ gas abundance is highest, are low and stable, with a mean of 108 (±18) K. The sub-solar SO₂ column density between longitudes of 90–220° varies from a low of 0.61 (±0.145) × 10^?17 cm^?2, near aphelion in 2004, to a high of 1.51 (±0.215) × 10¹⁷ cm^?2 in 2010 when Jupiter was approaching its early 2011 perihelion. No correlation in the gas temperature was seen with the increasing SO₂ column densities outside the errors.Assuming that any volcanic component of the atmosphere is constant with time, the correlation of increasing SO₂ abundance with decreasing heliocentric distance provides good evidence that the atmosphere is at least partially supported by frost sublimation. The SO₂ frost thermal inertias and albedos that fit the variation in atmospheric density best are between 150–1250 W m^?2 s^?1/2 K^?1 and 0.613–0.425 respectively. Photometric evidence favors albedos near the upper end of this range, corresponding to thermal inertias near the lower end. This relatively low frost thermal inertia produces larger amplitude seasonal variations than are observed, which in turn implies a substantial additional volcanic atmospheric component to moderate the amplitude of the seasonal variations of the total atmosphere on the anti-Jupiter hemisphere. The seasonal thermal inertia we measure is unique both because it refers exclusively to the SO₂ frost surface component, and also because it refers to relatively deep subsurface layers (few meters) due to the timescales of many years, while previous studies have determined thermal inertias at shallower levels (few centimeters), relevant for timescales of ～2 h (eclipse) or ～2 days (diurnal curves).