On the vertical wind shear of Saturn's Equatorial Jet at cloud level

In this work we analyze and compare the vertical cloud structure of Saturn's Equatorial Zone in two different epochs: the first one close to the Voyagers flybys (1979-1981) and the second one in 2004, when the Cassini spacecraft entered its orbit around the planet. Our goal is to retrieve the altitude of cloud features used as zonal wind tracers in both epochs. We reanalyze three different sets of photometrically calibrated published data: ground-based in 1979, Voyager 2 PPS and ISS observations in 1981, and we analyze a new set of Hubble Space Telescope images for 2004. For all situations we reproduced the observed reflectivity by means of a similar vertical model with three layers. The results indicate the presence of a changing tropospheric haze in 1979-1981 (Ptop∼100 mbar, τ∼10) and in 2004 (Ptop∼50 mbar, τ∼15) where the tracers are embedded. According to this model the Voyager 2 ISS images locate cloud tracers moving with zonal velocities of 455 to 465 (±2) m/s at a pressure level of 360 ± 140 mbar. For HST observations, our previous works had showed cloud tracers moving with zonal wind speeds of 280±10 m/s at a pressure level of about 50±10 mbar. All these values are calculated in the same region (3°±2° N). This speed difference, if interpreted as a vertical wind shear, requires a change of per scale height, two times greater than that estimated from temperature observations. We also perform an initial guess on Cassini ISS vertical sounding levels, retrieving values compatible with HST ones and Cassini CIRS derived vertical wind shear, but not with Voyager wind measurements. We conclude that the wind speed velocity differences measured between 1979-1981 and 2004 cannot be explained as a wind shear effect alone and demand dynamical processes.     

Numerical models of Saturn's long-lived anticyclones

New measurements of the dynamical properties of the long-lived Saturn's anticyclonic vortex known as “Brown Spot” (BS), discovered during the Voyager 1 and 2 flybys in 1980-1981 at latitude 43.1° N, and model simulations using the EPIC code, have allowed us to constrain the vertical wind shear and static stability in Saturn's atmosphere (vertically from pressure levels from 10 mbar to 10 bars) at this latitude. BS dynamical parameters from Voyager images include its size as derived from cloud albedo gradient (6100 km East-West times 4300 km North-South), mean tangential velocity ( at 2400 km from center) and mean vorticity (4.0±1.5×¹⁰⁻⁵ s⁻¹), lifetime >1 year, drift velocity relative to Voyager's System III rotation rate, mean meridional atmospheric wind profile at cloud level at its latitude and interactions with nearby vortices (pair orbiting and merging). An extensive set of numerical experiments have been performed to try to reproduce this single vortex properties and its observed mergers with smaller anticyclones by varying the vertical structure of the zonal wind and adjusting the static stability of the lower stratosphere and upper troposphere. Within the context of the EPIC model atmosphere, our simulations indicate that BS's drift velocity, longevity and merging behavior are very sensitive to these two atmospheric properties. The best results at the BS latitude occur for static stability conditions that use a Brunt-Väisäla frequency constant in the upper troposphere (from 0.5 to 10 bar) above 3.2×¹⁰⁻³ s⁻¹ and suggest that the wind speed slightly decays below the visible cloud deck from ∼0.5 to 10 bar at a rate per scale height. Changing the vortex latitude within the band domain introduces latitude oscillations in the vortex but not a significant meridional migration. Simulated mergers always showed orbiting movements with a typical merging time of about three days, very close to the time-span observed in the interaction of real vortices. Although these results are not unique in view of the unknowns of Saturn's deep atmosphere, they serve to constrain realistically its structure for ongoing Cassini observations.     

The three-dimensional structure of Saturn's equatorial jet at cloud level

The three-dimensional structure of Saturn's intense equatorial jet from latitudes 8° N to 20° S is revealed from detailed measurements of the motions and spectral reflectivity of clouds at visible wavelengths on high resolution images obtained by the Cassini Imaging Science Subsystem (ISS) in 2004 and early 2005. Cloud speeds at two altitude levels are measured in the near infrared filters CB2 and CB3 matching the continuum (effective wavelengths 750 and 939 nm) and in the MT2 and MT3 filters matching two methane absorption bands (effective wavelengths 727 and 889 nm). Radiative transfer models in selective filters covering an ample spectral range (250-950 nm) require the existence of two detached aerosol layers in the equator: an uppermost thin stratospheric haze extending between the pressure levels ∼20 and 40 mbar (tropopause level) and below it, a dense tropospheric haze-cloud layer extending between 50 mbar and the base of the ammonia cloud (between ∼1 and 1.4 bar). Individual cloud elements are detected and tracked in the tropospheric dense haze at 50 and 700 mbar (altitude levels separated by 142 km). Between latitudes 5° N and 12° S the winds increase their velocity with depth from 265 m s⁻¹ at the 50 mbar pressure level to 365 m s⁻¹ at 700 mbar. These values are below the high wind speeds of 475 m s⁻¹ measured at these latitudes during the Voyager era in 1980-1981, indicating that the equatorial jet has suffered a significant intensity change between that period and 1996-2005 or that the tracers of the flow used in the Voyager images were rooted at deeper levels than those in Cassini images.     

Saturn's south polar vortex compared to other large vortices in the Solar System

Observations made by the Imaging Science Subsystem (ISS), Visible and Infrared Mapping Spectrometer (VIMS) and the long-wavelength Composite Infrared Spectrometer (CIRS) aboard the Cassini spacecraft reveal that the large, long-lived cyclonic vortex at Saturn's south pole has a 4200-km-diameter cloud-free nearly circular region. This region has a 4 K warm core extending from the troposphere into the stratosphere, concentric cloud walls extending 20-70 km above the internal clouds, and numerous external clouds whose anticyclonic vorticity suggests a convective origin. The rotation speeds of the vortex reach . The Saturn polar vortex has features in common with terrestrial hurricanes and with the Venus polar vortex. Neptune and other giant planets may also have strong polar vortices.     

Gravitational tidal waves in Titan's upper atmosphere

Tidal waves driven by Titan's orbital eccentricity through the time-dependent component of Saturn's gravitational potential attain nonlinear, saturation amplitudes (|T^′|>10 K, , and ) in the upper atmosphere (?500 km) due to the approximate exponential growth as the inverse square root of pressure. The gravitational tides, with vertical wavelengths of ∼100-150 km above 500 km altitude, carry energy fluxes sufficient in magnitude to affect the energy balance of the upper atmosphere with heating rates in the altitude range of 500-900 km.     

On the HCN and CO2 abundance and distribution in Jupiter's stratosphere

Observations of Jupiter by Cassini/CIRS, acquired during the December 2000 flyby, provide the latitudinal distribution of HCN and CO₂ in Jupiter's stratosphere with unprecedented spatial resolution and coverage. Following up on a preliminary study by Kunde et al. [Kunde, V.G., and 41 colleagues, 2004. Science 305, 1582-1587], the analysis of these observations leads to two unexpected results (i) the total HCN mass in Jupiter's stratosphere in 2000 was (6.0±1.5)×¹⁰¹³ g, i.e., at least three times larger than measured immediately after the Shoemaker-Levy 9 (SL9) impacts in July 1994 and (ii) the latitudinal distributions of HCN and CO₂ are strikingly different: while HCN exhibits a maximum at 45° S and a sharp decrease towards high Southern latitudes, the CO₂ column densities peak over the South Pole. The total CO₂ mass is (2.9±1.2)×¹⁰¹³ g. A possible cause for the HCN mass increase is its production from the photolysis of NH₃, although a problem remains because, while millimeter-wave observations clearly indicate that HCN is currently restricted to submillibar (∼0.3 mbar) levels, immediate post-impact infrared observations have suggested that most of the ammonia was present in the lower stratosphere near 20 mbar. HCN appears to be a good atmospheric tracer, with negligible chemical losses. Based on 1-dimensional (latitude) transport models, the HCN distribution is best interpreted as resulting from the combination of a sharp decrease (over an order of magnitude in Kyy) of wave-induced eddy mixing poleward of 40° and an equatorward transport with velocity. The CO₂ distribution was investigated by coupling the transport model with an elementary chemical model, in which CO₂ is produced from the conversion of water originating either from SL9 or from auroral input. The auroral source does not appear adequate to reproduce the CO₂ peak over the South Pole, as required fluxes are unrealistically high and the shape of the CO₂ bulge is not properly matched. In contrast, the CO₂ distribution can be fit by invoking poleward transport with a velocity and vigorous eddy mixing (). While the vertical distribution of CO₂ is not measured, the combined HCN and CO₂ results imply that the two species reside at different stratospheric levels. Comparing with the circulation regimes predicted by earlier radiative-dynamical models of Jupiter's stratosphere, and with inferences from the ethane and acetylene stratospheric latitudinal distribution, we suggest that CO₂ lies in the middle stratosphere near or below the 5-mbar level.     

Saturn eddy momentum fluxes and convection: First estimates from Cassini images

We apply an automated cloud feature tracking algorithm to estimate eddy momentum fluxes in Saturn's southern hemisphere from Cassini Imaging Science Subsystem near-infrared continuum image sequences. Voyager Saturn manually tracked images had suggested no conversion of eddy to mean flow kinetic energy, but this was based on a small sample of <1000 wind vectors. The automated procedure we use for the Cassini data produces an order of magnitude more usable wind vectors with relatively unbiased sampling. Automated tracking is successful in and around the westward jet latitudes on Saturn but not in the vicinity of most eastward jets, where the linearity and non-discrete nature of cloud features produces ambiguous results. For the regions we are able to track, we find peak eddy fluxes and a clear positive correlation between eddy momentum fluxes and meridional shear of the mean zonal wind, implying that eddies supply momentum to eastward jets and remove momentum from westward jets at a rate . The behavior we observe is similar to that seen on Jupiter, though with smaller eddy-mean kinetic energy conversion rates per unit mass of atmosphere (). We also use the appearance and rapid evolution of small bright features at continuum wavelengths, in combination with evidence from weak methane band images where possible, to diagnose the occurrence of moist convective storms on Saturn. Areal expansion rates imply updraft speeds of over the convective anvil cloud area. As on Jupiter, convection preferentially occurs in cyclonic shear regions on Saturn, but unlike Jupiter, convection is also observed in eastward jet regions. With one possible exception, the large eddy fluxes seen in the cyclonic shear latitudes do not seem to be associated with convective events.     

Cassini RADAR constraint on Titan's winter polar precipitation

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

Dissipation of Titan's south polar clouds

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

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].     

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.     

First ENA observations at Mars: ENA emissions from the martian upper atmosphere

The neutral particle detector (NPD) on board Mars Express has observed energetic neutral atoms (ENAs) from a broad region on the dayside of the martian upper atmosphere. We show one such example for which the observation was conducted at an altitude of 570 km, just above the induced magnetosphere boundary (IMB). The time of flight spectra of these ENAs show that they had energies of 0.2-2 keV/amu, with an average energy of ∼1.1 keV/amu. Both the spatial distribution and the energy of these ENAs are consistent with the backscattered ENAs, produced by an ENA albedo process. This is the first observation of backscattered ENAs from the martian upper atmosphere. The origin of these ENAs is considered to be the solar wind ENAs that are scattered back by collision processes in the martian upper atmosphere. The particle flux and energy flux of the backscattered ENAs are and , respectively.     

Titan's hydrodynamically escaping atmosphere: Escape rates and the structure of the exobase region

In Strobel [Strobel, D.F., 2008. Icarus, 193, 588-594] a mass loss rate from Titan's upper atmosphere, , was calculated for a single constituent, N₂ atmosphere by hydrodynamic escape as a high density, slow outward expansion driven principally by solar UV heating due to CH₄ absorption. It was estimated, but not proven, that the hydrodynamic mass loss is essentially CH₄ and H₂ escape. Here the individual conservation of momentum equations for the three major components of the upper atmosphere (N₂, CH₄, H₂) are solved in the low Mach number limit and compared with Cassini Ion Neutral Mass Spectrometer (INMS) measurements to demonstrate that light gases (CH₄, H₂) preferentially escape over the heavy gas (N₂). The lightest gas (H₂) escapes with a flux 99% of its limiting flux, whereas CH₄ is restricted to ?75% of its limiting flux because there is insufficient solar power to support escape at the limiting rate. The respective calculated H₂ and CH₄ escape rates are 9.2×¹⁰²⁷ and 1.7×¹⁰²⁷ s−1, for a total of . From the calculated densities, mean free paths of N₂, CH₄, H₂, and macroscopic length scales, an extended region above the classic exobase is inferred where frequent collisions are still occurring and thermal heat conduction can deliver power to lift the escaping gas out of the gravitational potential well. In this region rapid acceleration of CH₄ outflow occurs. With the thermal structure of Titan's thermosphere inferred from INMS data by Müller-Wodarg et al. [Müller-Wodarg, I.C.F., Yelle, R.V., Cui, J., Waite Jr., J.H., 2008. J. Geophys. Res. 113, doi:10.1029/2007JE003033. E10005], in combination with calculated temperature profiles that include sputter induced plasma heating at the exobase, it is concluded that on average that the integrated, globally average, orbit-averaged, plasma heating rate during the Cassini epoch does not exceed ().     

Saturn’s zonal wind profile in 2004-2009 from Cassini ISS images and its long-term variability

Five years of Cassini Imaging Science Subsystem images, from 2004 to 2009, are analyzed in this work to retrieve global zonal wind profiles of Saturn’s northern and southern hemispheres in the methane absorbing bands at 890 and 727 nm and in their respective adjacent continuum wavelengths of 939 and 752 nm. A complete view of Saturn’s global circulation, including the equator, at two pressure levels, in the tropopause (60 mbar to 250 mbar with the MT filters) and in the upper troposphere (from ∼350 mbar to ∼500 mbar with the CB filter set), is presented. Both zonal wind profiles (available at the Supplementary Material Section), show the same structure but with significant differences in the peak of the eastward jets and the equatorial region, including a region of positive vertical shear symmetrically located around the equator between the 10° < |φ_c| < 25° where zonal velocities close to the tropopause are higher than at 500 mbar. A comparison of previously published zonal wind sets obtained by Voyager 1 and 2 (1980-1981), Hubble Space Telescope, and ground-based telescopes (1990-2004) with the present Cassini profiles (2004-2009) covering a full Saturn year shows that the shape of the zonal wind profile and intensity of the jets has remained almost unchanged except at the equator, despite the seasonal insolation cycle and the variability of Saturn’s emitted power. The major wind changes occurred at equatorial latitudes, perhaps following the Great White Spot eruption in 1990. It is not evident from our study if the seasonal insolation cycle and its associated ring shadowing influence the equatorial circulation at cloud level.     

Interaction between eddies and mean flow in Jupiter's atmosphere: Analysis of Cassini imaging data

Beebe et al. [Beebe, R.F., et al., 1980. Geophys. Res. Lett. 17, 1-4] and Ingersoll et al. [Ingersoll, A.P., et al., 1981. J. Geophys. Res. 86, 8733-8743] used images from Voyagers 1 and 2 to analyze the interaction between zonal winds and eddies in Jupiter's atmosphere. They reported a high positive correlation between Jupiter's eddy momentum flux, , and the variation of zonal velocity with latitude, . This correlation implied a surprisingly high rate of conversion of energy from eddies to zonal flow: , a value more than 10% of Jupiter's thermal flux emission. However, Sromovsky et al. [Sromovsky, L.A., et al., 1982. J. Atmos. Sci. 39, 1413-1432] argued that possible biases in the analysis could have caused an artificially high correlation. In addition, significant differences in the derived eddy flux between datasets put into question the robustness of any one result. We return to this long-standing puzzle using images of Jupiter from the Cassini flyby of December 2000. Our method is similar to previous analyses, but utilizes an automatic feature tracker instead of the human eye. The number of velocity vectors used in this analysis is over 200,000, compared to the 14,000 vectors used by Ingersoll et al. We also find a positive correlation between and and derive a global average power per unit mass, , ranging from . Utilizing Ingersoll et al.'s estimate of the mass per unit area involved in the transport, this would imply a rate of energy conversion of . We discuss the implications of this result and employ several tests to demonstrate its robustness.     

Diagnostics of Titan's stratospheric dynamics using Cassini/CIRS data and the 2-dimensional IPSL circulation model

The dynamics of Titan's stratosphere is discussed in this study, based on a comparison between observations by the CIRS instrument on board the Cassini spacecraft, and results of the 2-dimensional circulation model developed at the Institute Pierre-Simon Laplace, available at http://www.lmd.jussieu.fr/titanDbase [Rannou, P., Lebonnois, S., Hourdin, F., Luz, D., 2005. Adv. Space Res. 36, 2194-2198]. The comparison aims at both evaluating the model's capabilities and interpreting the observations concerning: (1) dynamical and thermal structure using temperature retrievals from Cassini/CIRS and the vertical profile of zonal wind at the Huygens landing site obtained by Huygens/DWE; and (2) vertical and latitudinal profiles of stratospheric gases deduced from Cassini/CIRS data. The modeled thermal structure is similar to that inferred from observations (Cassini/CIRS and Earth-based observations). However, the upper stratosphere (above 0.05 mbar) is systematically too hot in the 2D-CM, and therefore the stratopause region is not well represented. This bias may be related to the haze structure and to misrepresented radiative effects in this region, such as the cooling effect of hydrogen cyanide (HCN). The 2D-CM produces a strong atmospheric superrotation, with zonal winds reaching 200 m s⁻¹ at high winter latitudes between 200 and 300 km altitude (0.1-1 mbar). The modeled zonal winds are in good agreement with retrieved wind fields from occultation observations, Cassini/CIRS and Huygens/DWE. Changes to the thermal structure are coupled to changes in the meridional circulation and polar vortex extension, and therefore affect chemical distributions, especially in winter polar regions. When a higher altitude haze production source is used, the resulting modeled meridional circulation is weaker and the vertical and horizontal mixing due to the polar vortex is less extended in latitude. There is an overall good agreement between modeled chemical distributions and observations in equatorial regions. The difference in observed vertical gradients of C₂H₂ and HCN may be an indicator of the relative strength of circulation and chemical loss of HCN. The negative vertical gradient of ethylene in the low stratosphere at 15° S, cannot be modeled with simple 1-dimensional models, where a strong photochemical sink in the middle stratosphere would be necessary. It is explained here by dynamical advection from the winter pole towards the equator in the low stratosphere and by the fact that ethylene does not condense. Near the winter pole (80° N), some compounds (C₄H₂, C₃H₄) exhibit an (interior) minimum in the observed abundance vertical profiles, whereas 2D-CM profiles are well mixed all along the atmospheric column. This minimum can be a diagnostic of the strength of the meridional circulation, and of the spatial extension of the winter polar vortex where strong descending motions are present. In the summer hemisphere, observed stratospheric abundances are uniform in latitude, whereas the model maintains a residual enrichment over the summer pole from the spring cell due to a secondary meridional overturning between 1 and 50 mbar, at latitudes south of 40-50° S. The strength, as well as spatial and temporal extensions of this structure are a difficulty, that may be linked to possible misrepresentation of horizontally mixing processes, due to the restricted 2-dimensional nature of the model. This restriction should also be kept in mind as a possible source of other discrepancies.     

Extension of atmospheric dust loading to high altitudes during the 2001 Mars dust storm: MGS TES limb observations

Mars Global Surveyor (MGS) visible (solarband bolometer) and thermal infrared (IR) spectral limb observations from the Thermal Emission Spectrometer (TES) support quantitative profile retrievals for dust opacity and particle sizes during the 2001 global dust event on Mars. The current analysis considers the behavior of dust lifted to altitudes above 30 km during the course of this storm; in terms of dust vertical mixing, particle sizes, and global distribution. TES global maps of visible (solarband) limb brightness at 60 km altitude indicate a global-scale, seasonally evolving (over 190-240° solar longitudes, L_S) longitudinal corridor of vertically extended dust loading (which may be associated with a retrograde propagating, wavenumber 1 Rossby wave). Spherical radiative transfer analysis of selected limb profiles for TES visible and thermal IR radiances provide quantitative vertical profiles of dust opacity, indicating regional conditions of altitude-increasing dust mixing ratios. Observed infrared spectral dependences and visible-to-infrared opacity ratios of dust scattering over 30-60 km altitudes indicate particle sizes characteristic of lower altitudes (cross-section weighted effective radius, ), during conditions of significant dust transport to these altitudes. Conditions of reduced dust loading at 30-60 km altitudes present smaller dust particle sizes . These observations suggest rapid meridional transport at 30-80 km altitudes, with substantial longitudinal variation, of dust lifted to these altitudes over southern hemisphere atmospheric regions characterized by extraordinary (m/s) vertical advection velocities. By L_S=230° dust loading above 50 km altitudes decreased markedly at southern latitudes, with a high altitude (60-80 km) haze of fine (likely) water ice particles appearing over 10°S-40°N latitudes.     

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.     

N2 escape rates from Pluto's atmosphere

Hydrodynamic escape of N₂ molecules from Pluto's atmosphere is calculated under the assumption of a high density, slow outflow expansion driven by solar EUV heating by N₂ absorption, near-IR and UV heating by CH₄ absorption, and CO cooling by rotational line emission as a function of solar activity. At 30 AU, the N₂ escape rate varies from in the absence of heating, but driven by an upward thermal heat conduction flux from the stratosphere, for lower boundary temperatures varying from 70-100 K. With solar heating varying from solar minimum to solar maximum conditions and a calculated lower boundary temperature, 88.2 K, the N₂ escape rate range is , respectively. LTE rotational line emission by CO reduces the net solar heat input by at most 35% and plays a minor role in lowering the calculated escape rates, but ensures that the lower boundary temperature can be calculated by radiative equilibrium with near-IR CH₄ heating. While an upward thermal conduction heat flux at the lower boundary plays a fundamental role in the absence of heating, with solar heating it is downward at solar minimum, and is, at most, 13% of the integrated net heating rate over the range of solar activity. For the arrival of the New Horizons spacecraft at Pluto in July 2015, predictions are lower boundary temperature, T₀∼81 K, and N₂ escape rate , and peak thermospheric temperature ∼103 K at 1890 km, based on expected solar medium conditions.     

Detection of CH3D on Titan

We report the detection of ¹³CH₃D in Titan's stratosphere from Cassini/CIRS infrared spectra near 8.7 μm. Fitting simultaneously the ν₆ bands of both ¹³CH₃D and ¹²CH₃D and the ν₄ band of CH₄, we derive a D/H ratio equal to and a ¹²C/¹³C ratio in deuterated methane of , consistent with that measured in normal methane.