Neptune’s cloud and haze variations 1994-2008 from 500 HST-WFPC2 images

The analysis of all suitable images taken of Neptune with the Wide Field Planetary Camera 2 on the Hubble Space Telescope between 1994 and 2008 revealed the following results. The activity of discrete cloud features located near Neptune’s tropopause remained roughly constant within each year but changed significantly on the time scale of ∼5 years. Discrete clouds covered 1% of the disk on average, but more than 2% in 2002. The other ∼99% of the disk probed Neptune’s hazes at lower altitudes. At red and near-infrared wavelengths, two dark bands around −70° and 10° latitude were perfectly steady and originated in the upper two scale heights of the troposphere, either by decreased haze opacity or by an increased methane relative humidity. At blue wavelengths, a dark band between −60° and −30° latitude was most obvious during the early years, caused by dark aerosols below the 3-bar level with single scattering albedos reduced by ∼0.04, and this contrast was constant between 410 and 630 nm wavelength. The dark band decayed exponentially with a time constant of 5 ± 1 years, which can be explained by settling of the dark aerosols at a rate of 1 bar pressure difference per year. The other latitudes brightened with the same time constant but lower amplitudes. The only exception was a darkening event in the 15-30° latitude region between 1994 and 1996, which coincides with two dark spots observed in the same region during the same time period, the only dark spots seen since Voyager. The dark aerosols had a similar latitudinal distribution as the discrete clouds near the tropopause, although both were separated by four scale heights. Photometric analysis revealed a phase coefficient of 0.0028 ± 0.0010 mag/deg for the 0-2° phase-angle range observable from Earth. Neptune’s sub-Earth latitude varied by less than 3° throughout the observation period providing a data set with almost constant viewing geometry. The trends observed up to 2008 continued into 2010 based on images taken with the Wide Field Camera 3.     

Vertical and meridional distribution of ethane, acetylene and propane in Saturn’s stratosphere from CIRS/Cassini limb observations

Measuring the spatial distribution of chemical compounds in Saturn’s stratosphere is critical to better understand the planet’s photochemistry and dynamics. Here we present an analysis of infrared spectra in the range 600-1400 cm⁻¹ acquired in limb geometry by the Cassini spacecraft between March 2005 and January 2008. We first determine the vertical temperature profiles from 3 to 0.01 hPa, at latitudes ranging from 70°N to 80°S. We infer a similar meridional temperature gradient at 1-2 hPa as in recent previous studies [Fletcher, L.N., Irwin, P.G.J., Teanby, N.A., Orton, G.S., Parrish, P.D., de Kok, R., Howett, C., Calcutt, S.B., Bowles, N., Taylor, F.W., 2007. Icarus 189, 457-478; Howett, C.J.A., Irwin, P.G.J., Teanby, N.A., Simon-Miller, A., Calcutt, S.B., Fletcher, L.N., de Kok, R., 2007. Icarus 190, 556-572]. We then retrieve the vertical profiles of C₂H₆ and C₂H₂ from 3 to 0.01 hPa and of C₃H₈ around 1 hPa. At 1 hPa, the meridional variation of C₂H₂ is found to follow the yearly averaged solar insolation, except for a strong equatorial mole fraction of 8×10^-7, nearly two times higher than expected. This enhancement in abundance can be explained by the descent of hydrocarbon-rich air, with a vertical wind speed at the equator of 0.25±0.1 mm/s at 1 hPa and 0.4±0.15 mm/s at 0.1 hPa. The ethane distribution is relatively uniform at 1 hPa, with only a moderate 25% increase from 35°S to 80°S. Propane is found to increase from north to south by a factor of 1.9, suggesting that its lifetime may be shorter than Saturn’s year at 1 hPa. At high altitudes (1 Pa), C₂H₂ and C₂H₆ abundances depart significantly from the photochemical model predictions of Moses and Greathouse [Moses, J.I., Greathouse, T.K., 2005. J. Geophys. Res. 110, 9007], except at high southern latitudes (62, 70 and 80°S) and near the equator. The observed abundances are found strongly depleted in the 20-40°S region and enhanced in the 20-30°N region, the latter coinciding with the ring’s shadow. We favor a dynamical explanation for these anomalies.     

Saturn’s tropospheric composition and clouds from Cassini/VIMS 4.6-5.1 μm nightside spectroscopy

The latitudinal variation of Saturn’s tropospheric composition (NH₃, PH₃ and AsH₃) and aerosol properties (cloud altitudes and opacities) are derived from Cassini/VIMS 4.6-5.1 μm thermal emission spectroscopy on the planet’s nightside (April 22, 2006). The gaseous and aerosol distributions are used to trace atmospheric circulation and chemistry within and below Saturn’s cloud decks (in the 1- to 4-bar region). Extensive testing of VIMS spectral models is used to assess and minimise the effects of degeneracies between retrieved variables and sensitivity to the choice of aerosol properties. Best fits indicate cloud opacity in two regimes: (a) a compact cloud deck centred in the 2.5-2.8 bar region, symmetric between the northern and southern hemispheres, with small-scale opacity variations responsible for numerous narrow light/dark axisymmetric lanes; and (b) a hemispherically asymmetric population of aerosols at pressures less than 1.4 bar (whose exact altitude and vertical structure is not constrained by nightside spectra) which is 1.5-2.0× more opaque in the summer hemisphere than in the north and shows an equatorial maximum between ±10° (planetocentric).Saturn’s NH₃ spatial variability shows significant enhancement by vertical advection within ±5° of the equator and in axisymmetric bands at 23-25°S and 42-47°N. The latter is consistent with extratropical upwelling in a dark band on the poleward side of the prograde jet at 41°N (planetocentric). PH₃ dominates the morphology of the VIMS spectrum, and high-altitude PH₃ at p < 1.3 bar has an equatorial maximum and a mid-latitude asymmetry (elevated in the summer hemisphere), whereas deep PH₃ is latitudinally-uniform with off-equatorial maxima near ±10°. The spatial distribution of AsH₃ shows similar off-equatorial maxima at ±7° with a global abundance of 2-3 ppb. VIMS appears to be sensitive to both (i) an upper tropospheric circulation (sensed by NH₃ and upper-tropospheric PH₃ and hazes) and (ii) a lower tropospheric circulation (sensed by deep PH₃, AsH₃ and the lower cloud deck).     

Improved measurement of Asteroid (4) Vesta’s rotational axis orientation

We report an improved measurement of the rotational axis orientation of Asteroid (4) Vesta. By analyzing and combining all previous measurements using a limb-fitting technique from ground/HST data collected from 1983 to 2006, we derive a pole solution of (RA = 304.5°, Dec = 41.5°). Images of Vesta acquired with the Wide Field Camera 3 onboard the Hubble Space Telescope (HST) in February 2010 are combined with images from the Wide Field Planetary Camera 2 on HST obtained in 1994, 1996, and 2007 at similar spatial resolution and wavelengths to perform new measurements. Control point stereogrammetry returns a pole solution of (305.1°, 43.4°). An alternate method tracks surface features and fits their projected paths with ellipses to determine a great circle containing the pole for each HST observation. Combined, the four great circles yield a pole solution of (309.3°, 41.9°). These three solutions obtained with almost independent methods are within 3.5° of each other, suggesting a robust solution. Combining the results from all three techniques, we propose an improved value of the rotational axis of Vesta as RA = 305.8° ± 3.1°, Dec = 41.4° ± 1.5° (1-σ error). This new solution changes from (301°, 41°) reported by Thomas et al. (Thomas, P.C., Binzel, R.P., Gaffey, M.J., Zellner, B.H., Storrs, A.D., Wells, E. [1997a]. Icarus 128, 88-94) by 3.6°, and from (306°, 38°) reported by Drummond and Christou (Drummond, J.D., Christou, J. [2008]. Icarus 197, 480-496) by 3.4°. It changes the obliquity of Vesta by up to ∼3°, but increases the Sun-centered RA of Vesta at equinox by ∼8°, and postpones the date of equinox by ∼35 days. The change of the pole position is less than the resolution of all previous images of Vesta, and should not change the main science conclusions of previous research about Vesta.     

Ground-based measurements of the methane distribution on Titan

Using spectra taken with NIRSPEC (Near Infrared Spectrometer) and adaptive optics on the Keck II telescope, we resolved the latitudinal variation of the 3ν₂ band of CH₃D at 1.56 μm. As CH₃D is less abundant than CH₄ by a factor of 50±10×10^-5, these CH₃D lines do not saturate in Titan’s atmosphere, and are well characterized by laboratory measurements. Thus they do not suffer from the large uncertainties of the CH₄ lines that are weak enough to be unsaturated in Titan. Our measurements of the methane abundance are confined to the latitude range of 32°S-18°N and longitudes sampled by a 0.04″ slit centered at ∼195°W. The methane abundance below 10 km is constant to within 20% in the tropical atmosphere sampled by our observations, consistent with the low surface insolation and lack of surface methane [Griffith, C.A., McKay, C.P., Ferri, F., 2008. Astrophys. J. 687, L41-L44].     

Jovian temperature and cloud variability during the 2009-2010 fade of the South Equatorial Belt

Mid-infrared 7-20 μm imaging of Jupiter from ESO’s Very Large Telescope (VLT/VISIR) demonstrate that the increased albedo of Jupiter’s South Equatorial Belt (SEB) during the ‘fade’ (whitening) event of 2009-2010 was correlated with changes to atmospheric temperature and aerosol opacity. The opacity of the tropospheric condensation cloud deck at pressures less than 800 mbar increased by 80% between May 2008 and July 2010, making the SEB (7-17°S) as opaque in the thermal infrared as the adjacent equatorial zone. After the cessation of discrete convective activity within the SEB in May 2009, a cool band of high aerosol opacity (the SEB zone at 11-15°S) was observed separating the cloud-free northern and southern SEB components. The cooling of the SEBZ (with peak-to-peak contrasts of 1.0 ± 0.5 K), as well as the increased aerosol opacity at 4.8 and 8.6 μm, preceded the visible whitening of the belt by several months. A chain of five warm, cloud-free ‘brown barges’ (subsiding airmasses) were observed regularly in the SEB between June 2009 and June 2010, by which time they too had been obscured by the enhanced aerosol opacity of the SEB, although the underlying warm circulation was still present in July 2010. Upper tropospheric temperatures (150-300 mbar) remained largely unchanged during the fade, but the cool SEBZ formation was detected at deeper levels (p > 300 mbar) within the convectively-unstable region of the troposphere. The SEBZ formation caused the meridional temperature gradient of the SEB to decrease between 2008 and 2010, reducing the vertical thermal windshear on the zonal jets bounding the SEB. The southern SEB had fully faded by July 2010 and was characterised by short-wave undulations at 19-20°S. The northern SEB persisted as a narrow grey lane of cloud-free conditions throughout the fade process.The cool temperatures and enhanced aerosol opacity of the SEBZ after July 2009 are consistent with an upward flux of volatiles (e.g., ammonia-laden air) and enhanced condensation, obscuring the blue-absorbing chromophore and whitening the SEB by April 2010. These changes occurred within cloud decks in the convective troposphere, and not in the radiatively-controlled upper troposphere. NH₃ ice coatings on aerosols at p < 800 mbar are plausible sources of the suppressed 4.8 and 8.6-μm emission, although differences in the spatial distribution of opacity at these two wavelengths suggest that enhanced attenuation by a deeper cloud (p > 800 mbar) also occurred during the fade. Revival of the dark SEB coloration in the coming months will ultimately require sublimation of these ices by subsidence and warming of volatile-depleted air.     

Changes in Jupiter’s zonal velocity between 1979 and 2008

We show that the peak velocity of Jupiter’s visible-cloud-level zonal winds near 24°N (planetographic) increased from 2000 to 2008. This increase was the only change in the zonal velocity from 2000 to 2008 for latitudes between ±70° that was statistically significant and not obviously associated with visible weather. We present the first automated retrieval of fast (∼130 m s⁻¹) zonal velocities at 8°N planetographic latitude, and show that some previous retrievals incorrectly found slower zonal winds because the eastward drift of the dark projections (associated with 5-μm hot spots) “fooled” the retrieval algorithms.We determined the zonal velocity in 2000 from Cassini images from NASA’s Planetary Data System using a global method similar to previous longitude-shifting correlation methods used by others, and a new local method based on the longitudinal average of the two-dimensional velocity field. We obtained global velocities from images acquired in May 2008 with the Wide Field Planetary Camera 2 (WFPC2) on the Hubble Space Telescope (HST). Longer-term variability of the zonal winds is based on comparisons with published velocities based on 1979 Voyager 2 and 1995-1998 HST images. Fluctuations in the zonal wind speeds on the order of 10 m s⁻¹ on timescales ranging from weeks to months were found in the 1979 Voyager 2 and the 1995-1998 HST velocities. In data separated by 10 h, we find that the east-west velocity uncertainty due to longitudinal fluctuations are nearly 10 m s⁻¹, so velocity fluctuations of 10 m s⁻¹ may occur on timescales that are even smaller than 10 h. Fluctuations across such a wide range of timescales limit the accuracy of zonal wind measurements. The concept of an average zonal velocity may be ill-posed, and defining a “temporal mean” zonal velocity as the average of several zonal velocity fields spanning months or years may not be physically meaningful.At 8°N, we use our global method to find peak zonal velocities of ∼110 m s⁻¹ in 2000 and ∼130 m s⁻¹ in 2008. Zonal velocities from 2000 Cassini data produced by our local and global methods agree everywhere, except in the vicinity of 8°N. There, the local algorithm shows that the east-west velocity has large variations in longitude; vast regions exceed ∼140 m s⁻¹. Our global algorithm, and all of the velocity-extraction algorithms used in previously-published studies, found the east-west drift velocities of the visible dark projections, rather than the true zonal velocity at the visible-cloud level. Therefore, the apparent increase in zonal winds between 2000 and 2008 at 8°N is not a true change in zonal velocity.At 7.3°N, the Galileo probe found zonal velocities of 170 m s⁻¹ at the 3-bar level. If the true zonal velocity at the visible-cloud level at this latitude is ∼140 m s⁻¹ rather than ∼105 m s⁻¹, then the vertical zonal wind shear is much less than the currently accepted value.     

Neptune’s rotational period suggested by the extraordinary stability of two features

The interior rotation and motions in giant planets have generally been probed only at radio wavelengths from spacecraft near the planet, except for Jupiter’s radio emission detectable from Earth. Here I suggest that Neptune’s interior can be indirectly probed at visible wavelength by tracking 10 features that are connected with a stationary latitudinal speed pattern of 7 m/s amplitude. All 10 features remained aligned at the same longitude throughout the Voyager observation period in 1989. Two of them, the South Polar Wave and South Polar Feature, have been observed from Earth for ∼20 years, but their extraordinary rotational stability was never recognized. They probably pinpoint Neptune’s rotational period (15.9663 ± 0.0002 h), one of the largest improvements in 346 years of measuring the giant planets’ rotations. The previous best estimate of Neptune’s rotational period (16.108 ± 0.006 h) was based on Voyager 2 radio data (Lecacheux, A., Zarka, P., Desch, M.D., Evans, D.R. [1993]. Geophys. Res. Lett. 20, 2711-2714). The new result suggests an upward revision of the mass of Neptune’s core. This finding may also question the accepted value of Uranus’ rotational period. The first reliable wind measurements within 15° of Neptune’s South Pole, based on tracking four features in Voyager images, show a 300 m/s eastward jet peaking near 76° South, while the area within 4° of the South Pole seems to be rotationally locked to the interior. These new observations of the stationary features and winds could address the long-standing question about the depth of the atmospheric circulation and may allow some constraints on convection currents in Neptune’s interior.     

HCN and CN in Comet 2P/Encke: Models of the non-isotropic, rotation-modulated coma and CN parent life time

Axisymmetric models of the outgassing of a cometary nucleus have been constructed. Such models can be used to describe a nucleus with a single active region. The models may include a solar zenith angle dependence of the outgassing. They retrieve the outgassing flux at distances from the nucleus where collisions between molecules are unimportant, as function of the angle with respect to the outgassing axis. The observed emissions must be optically thin. Furthermore the models assume that the outflow speed at large distance from the nucleus does not depend on direction. The value of the outflow speed is retrieved. The models are applied to CN images and HCN spectra of Comet 2P/Encke, obtained nearly simultaneously in November 2003 with the 2 m optical telescope on Mount Rozhen, Bulgaria, and with the 10 m Heinrich Hertz Submillimeter Telescope on Mount Graham, Arizona, USA. According to Sekanina (1988), Astron. J. 95, 911-924, at that time a single outgassing source was active. Input parameters to the models like the rotation period of the nucleus and a small correction to Sekanina’s rotation axis are determined from a simpler jet position angle model. The rotation is prograde with a sideric period of 11.056 ± 0.024 h, in agreement with literature values. The best fit model has an outflow speed of 0.95 ± 0.04 km s⁻¹. The same value has been derived from the corkscrew appearing in the CN images. The location of the outgassing axis is at colatitude δ_a = 7.4° ± 2.9° and longitude λ_a = 235° ± 17° (a definition of zero longitude is provided). Comet Encke’s outgassing corresponds approximately to the longitudinally averaged solar input on a spherical nucleus (i.e. very likely comes from deeper layers) but with some deficiency of outgassing at mid-latitudes and non-zero outgassing from the dark polar cap. The presence of gas flow from the dark polar cap is explained as evidence of gas flow across the terminator. The models rely mostly on the CN images. The HCN spectra are more noisy. They provide information how to determine the best fit outflow velocity and the sense of rotation. The model HCN spectra are distinctly non-Gaussian. Within error limits they are consistent with the observations. Models based solely on the HCN spectra are also presented but, because of the lower quality of the data and the unfavorable observing geometry, yield inferior results. As a by-product we determine the CN parent life time from our CN observations. The solar EUV and Lyα radiation field at the time of our observations is taken into account.     

Multi-wavelength simulations of atmospheric radiation from Io with a 3-D spherical-shell backward Monte Carlo radiative transfer model

Conflicting observations regarding the dominance of either sublimation or volcanism as the source of the atmosphere on Io and disparate reports on the extent of its spatial distribution and the absolute column abundance invite the development of detailed computational models capable of improving our understanding of Io’s unique atmospheric structure and origin. Improving upon previous models, Walker et al. (Walker, A.C., Gratiy, S.L., Levin, D.A., Goldstein, D.B., Varghese, P.L., Trafton, L.M., Moore, C.H., Stewart, B. [2009]. Icarus) developed a fully 3-D global rarefied gas dynamics model of Io’s atmosphere including both sublimation and volcanic sources of SO₂ gas. The fidelity of the model is tested by simulating remote observations at selected wavelength bands and comparing them to the corresponding astronomical observations of Io’s atmosphere. The simulations are performed with a new 3-D spherical-shell radiative transfer code utilizing a backward Monte Carlo method. We present: (1) simulations of the mid-infrared disk-integrated spectra of Io’s sunlit hemisphere at 19 μm, obtained with TEXES during 2001-2004; (2) simulations of disk-resolved images at Lyman-α obtained with the Hubble Space Telescope (HST), Space Telescope Imaging Spectrograph (STIS) during 1997-2001; and (3) disk-integrated simulations of emission line profiles in the millimeter wavelength range obtained with the IRAM-30 m telescope in October-November 1999. We found that the atmospheric model generally reproduces the longitudinal variation in band depth from the mid-infrared data; however, the best match is obtained when our simulation results are shifted ∼30° toward lower orbital longitudes. The simulations of Lyman-α images do not reproduce the mid-to-high latitude bright patches seen in the observations, suggesting that the model atmosphere sustains columns that are too high at those latitudes. The simulations of emission line profiles in the millimeter spectral region support the hypothesis that the atmospheric dynamics favorably explains the observed line widths, which are too wide to be formed by thermal Doppler broadening alone.     

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

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

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.     

Spatially-resolved high-resolution spectroscopy of Venus 1. Variations of CO2, CO, HF, and HCl at the cloud tops

Variations of the upper cloud boundary and the CO, HF, and HCl mixing ratios were observed using the CSHELL spectrograph at NASA IRTF. The observations were made in three sessions (October 2007, January 2009, and June 2009) at early morning and late afternoon on Venus in the latitude range of ±60°. CO₂ lines at 2.25 μm reveal variations of the cloud aerosol density (∼25%) and scale height near 65 km. The measured reflectivity of Venus at low latitudes is 0.7 at 2.25 μm and 0.028 at 3.66 μm, and the effective CO₂ column density is smaller at 3.66 μm than those at 2.25 μm by a factor of 4. This agrees with the almost conservative multiple scattering at 2.25 μm and single scattering in the almost black aerosol at 3.66 μm. The expected difference is just a factor of (1 − g)⁻¹ = 4, where g = 0.75 is the scattering asymmetry factor for Venus’ clouds. The observed CO mixing ratio is 52 ± 4 ppm near 08:00 and 40 ± 4 ppm near 16:30 at 68 km, and the higher ratio in the morning may be caused by extension of the CO morningside bulge to the cloud tops. The observed weak limb brightening in CO indicates an increase of the CO mixing ratio with altitude. HF is constant at 3.5 ± 0.2 ppb at 68 km in both morningside and afternoon observations and in the latitude range ±60°. Therefore the observations do not favor a bulge of HF, though HF is lighter than CO. Probably a source in the upper atmosphere facilitates the bulge formation. The recent measurements of HCl near 70 km are controversial (0.1 and 0.74 ppm) and require either a strong sink or a strong source of HCl in the clouds. The HCl lines of the (2-0) band are blended by the solar and telluric lines. Therefore we observed the P8 lines of the (1-0) band at 3.44 μm. These lines are spectrally clean and result in the HCl mixing ratio of 0.40 ± 0.03 ppm at 74 km. HCl does not vary with latitude within ±60°. Our observations support a uniformly mixed HCl throughout the Venus atmosphere.     

Orbital effects on Mercury’s escaping sodium exosphere

We present results from coronagraphic imaging of Mercury’s sodium tail over a 7° field of view. Several sets of observations made at the McDonald Observatory since May 2007 show a tail of neutral sodium atoms stretching more than 1000 Mercury radii (R_m) in length, or a full degree of sky. However, no tail was observed extending beyond 120 R_m during the January 2008 MESSENGER fly-by period, or during a similar orbital phase of Mercury in July 2008. Large changes in Mercury’s heliocentric radial velocity cause Doppler shifts about the Fraunhofer absorption features; the resultant change in solar flux and radiation pressure is the primary cause of the observed variation in tail brightness. Smaller fluctuations in brightness may exist due to changing source rates at the surface, but we have no explicit evidence for such changes in this data set. The effects of radiation pressure on Mercury’s escaping atmosphere are investigated using seven observations spanning different orbital phases. Total escape rates of atmospheric sodium are estimated to be between 5 and 13 × 10²³ atoms/s and show a correlation to radiation pressure. Candidate sources of Mercury’s sodium exosphere include desorption by UV sunlight, thermal desorption, solar wind channeled along Mercury’s magnetic field lines, and micro-meteor impacts. Wide-angle observations of the full extent of Mercury’s sodium tail offer opportunities to enhance our understanding of the time histories of these source rates.     

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.     

Meridional distribution of CH3C2H and C4H2 in Saturn’s stratosphere from CIRS/Cassini limb and nadir observations

Limb and nadir spectra acquired by Cassini/CIRS (Composite InfraRed Spectrometer) are analyzed in order to derive, for the first time, the meridional variations of diacetylene (C₄H₂) and methylacetylene (CH₃C₂H) mixing ratios in Saturn’s stratosphere, from 5 hPa up to 0.05 hPa and 80°S to 45°N. We find that the C₄H₂ and CH₃C₂H meridional distributions mimic that of acetylene (C₂H₂), exhibiting small-scale variations that are not present in photochemical model predictions. The most striking feature of the meridional distribution of both molecules is an asymmetry between mid-southern and mid-northern latitudes. The mid-southern latitudes are found depleted in hydrocarbons relative to their northern counterparts. In contrast, photochemical models predict similar abundances at north and south mid-latitudes. We favor a dynamical explanation for this asymmetry, with upwelling in the south and downwelling in the north, the latter coinciding with the region undergoing ring shadowing. The depletion in hydrocarbons at mid-southern latitudes could also result from chemical reactions with oxygen-bearing molecules.Poleward of 60°S, at 0.1 and 0.05 hPa, we find that the CH₃C₂H and C₄H₂ abundances increase dramatically. This behavior is in sharp contradiction with photochemical model predictions, which exhibit a strong decrease towards the south pole. Several processes could explain our observations, such as subsidence, a large vertical eddy diffusion coefficient at high altitudes, auroral chemistry that enhances CH₃C₂H and C₄H₂ production, or shielding from photolysis by aerosols or molecules produced from auroral chemistry. However, problems remain with all these hypotheses, including the lack of similar behavior at lower altitudes.Our derived mean mixing ratios at 0.5 hPa of (2.4 ± 0.3) × 10⁻¹⁰ for C₄H₂ and of (1.1 ± 0.3) × 10⁻⁹ for CH₃C₂H are compatible with the analysis of global-average ISO observations performed by Moses et al. (Moses, J.I., Bézard, B., Lellouch, E., Gladstone, G.R., Feuchtgruber, H., Allen, M. [2000a]. Icarus 143, 244-298). Finally, we provide values for the ratios [CH₃C₂H]/[C₂H₂] and [C₄H₂]/[C₂H₂] that can constrain the coupled chemistry of these hydrocarbons.     

The aftermath of the July 2009 impact on Jupiter: Ammonia, temperatures and particulates from Gemini thermal infrared spectroscopy

We obtained longitudinally resolved thermal infrared spectra (8-13 μm and 17-25 μm) of Jupiter’s impact debris at the Gemini South Telescope on July 24, 2009; five days after the July 19th collision. These were used to study the mechanisms responsible for the redistribution of thermal energy and material (ammonia and stratospheric particulates) following the impact. Upwelling of (8.5 ± 4.1) × 10¹⁴ g of tropospheric air was sufficient to deposit (6.7 ± 4.1) × 10¹² g of NH₃ over a 6° longitude range above the impact core. The NH₃ was distributed over the 20-80 mbar region with a peak abundance of 1.0 ± 0.6 ppm at 45 mbar. Only a 10th of this abundance was observed over the western ejecta, and it is unlikely that these observations were sensitive to NH₃ entrained in the ballistic plume itself. The pattern of excess thermal energy was markedly different from that of Shoemaker-Levy 9 (SL9), with a localized tropospheric perturbation of 2.0 ± 1.0 K at 200-300 mbar and a broader stratospheric warming of up to 3.5 ± 2.0 K at 10-30 mbar. We find no evidence of residual warmth at p < 1 mbar five days after the impact. The excess thermal energy places lower limits on the total energy of the impact (1.8-15.7 × 10²⁶ ergs), which limits the impactor diameter to 70-510 m (depending on the bulk density chosen for the material).The models of the Gemini spectra required three distinct aerosol features, indicative of the mineralogy of the dark particulate debris, centred at 9.1, 10.0 and 18.5 μm. The retrieved opacities for each of these features were distributed over a larger area (9-10° longitude) and at higher altitudes (above the 10-mbar level) than the stratospheric NH₃, and they are more spatially inhomogeneous. This implies the particulates were either entrained with the rising hot plume or created upon plume re-entry and are subsequently redistributed by stratospheric winds. The three particulate features were consistent with a mixture of amorphous iron and magnesium-rich silicates and silicas in the debris field. A broad 10-μm signature was coincident with peaks expected from material rich in amorphous olivines (but poor in pyroxenes), and similar to silicate features observed during SL9. A narrow 9.1-μm signature was interpreted as a combination of amorphous and crystalline silica. Finally, a broad 18.5-μm emitter was not adequately reproduced by a mixture of simple olivines and pyroxenes and remains to be identified.     

Titan’s atomic hydrogen corona

Based on measurements performed by the Hydrogen Deuterium Absorption Cell (HDAC) aboard the Cassini orbiter, Titan’s atomic hydrogen exosphere is investigated. Data obtained during the T9 encounter are used to infer the distribution of atomic hydrogen throughout Titan’s exosphere, as well as the exospheric temperature.The measurements performed during the flyby are modeled by performing Monte Carlo radiative transfer calculations of solar Lyman-α radiation, which is resonantly scattered on atomic hydrogen in Titan’s exosphere. Two different atomic hydrogen distribution models are applied to determine the best fitting density profile. One model is a static model that uses the Chamberlain formalism to calculate the distribution of atomic hydrogen throughout the exosphere, whereas the second model is a Particle model, which can also be applied to non-Maxwellian velocity distributions.The density distributions provided by both models are able to fit the measurements although both models differ at the exobase: best fitting exobase atomic hydrogen densities of n_H = (1.5 ± 0.5) × 10⁴ cm⁻³ and n_H = (7 ± 1) × 10⁴ cm⁻³ were found using the density distribution provided by both models, respectively. This is based on the fact that during the encounter, HDAC was sensitive to altitudes above about 3000 km, hence well above the exobase at about 1500 km. Above 3000 km, both models produce densities which are comparable, when taking into account the measurement uncertainty.The inferred exobase density using the Chamberlain profile is a factor of about 2.6 lower than the density obtained from Voyager 1 measurements and much lower than the values inferred from current photochemical models. However, when taking into account the higher solar activity during the Voyager flyby, this is consistent with the Voyager measurements. When using the density profile provided by the particle model, the best fitting exobase density is in perfect agreement with the densities inferred by current photochemical models.Furthermore, a best fitting exospheric temperature of atomic hydrogen in the range of T_H = (150-175) ± 25 K was obtained when assuming an isothermal exosphere for the calculations. The required exospheric temperature depends on the density distribution chosen. This result is within the temperature range determined by different instruments aboard Cassini. The inferred temperature is close to the critical temperature for atomic hydrogen, above which it can escape hydrodynamically after it diffused through the heavier background gas.     

Dynamics of Jupiter’s equatorial region at cloud top level from Cassini and HST images

We present a study of the equatorial region of Jupiter, between latitudes ∼15°S and ∼15°N, based on Cassini ISS images obtained during the Jupiter flyby at the end of 2000, and HST images acquired in May and July 2008. We examine the structure of the zonal wind profile and report the detection of significant longitudinal variations in the intensity of the 6°N eastward jet, up to 60 m s⁻¹ in Cassini and HST observations. These longitudinal variations are, in the HST case, associated with different cloud morphology. Photometric and radiative transfer analysis of the cloud features used as tracers in HST images show that at most there is only a small height difference, no larger than ∼0.5-1 scale heights, between the slow (∼100 m s⁻¹) and fast (∼150 m s⁻¹) moving features. This suggests that speed variability at 6°N is not dominated by vertical wind shears but instead we propose that Rossby wave activity is the responsible for the zonal variability. Removing this variability, we find that Jupiter’s equatorial jet is actually symmetric relative to equator with two peaks of ∼140-150 m s⁻¹ located at latitudes 6°N and 6°S and at a similar pressure level. We also study the local dynamics of particular equatorial features such as several dark projections associated with 5 μm hot spots and a large, long-lived feature called the White Spot (WS) located at 6°S. Convergent flow at the dark projections appears to be a characteristic which depends on the particular morphology and has only been detected in some cases. The internal flow field in the White Spot indicates that it is a weakly rotating quasi-equatorial anticyclone relative to the ambient meridionally sheared flow.     

Comparing Phoebe’s 2005 opposition surge in four visible light filters

We observed Phoebe for 13 nights over a period of 55 days before, during, and after the 2005 Saturn opposition with the New Mexico State University (NMSU) 1-m telescope at Apache Point Observatory (APO) in Sunspot, NM and characterized the width and magnitude of Phoebe’s opposition surge in BVRI filters. Our observations cover a phase angle range of 4.87° to 0.0509°. We use a Hapke reflectance model incorporating shadow hiding and coherent backscatter to investigate the wavelength dependence of Phoebe’s opposition surge. We find a significant opposition surge magnitude of 55-58% between phase angles of 5° and 0°. We find the strongest opposition surge for phase angles less than 2° in the I-band. The coherent backscatter angular width is on the order of 0.50°. We find Phoebe’s albedo to be spectrally flat within our error limits, with a B-band albedo of 0.0855 ± 0.0031, a V-band albedo of 0.0856 ± 0.0023, an R-band albedo of 0.0843 ± 0.0020, and an I-band albedo of 0.0839 ± 0.0023. We compare Phoebe’s albedo, color, and opposition surge magnitudes and slopes with those of other outer solar system bodies and find similarities to Centaurs, Nereid, Puck, and Comets 19P/Borrelly, 9P/Tempel 1, and 81P/Wild 2. We find that this comparison supports the idea that Phoebe originated in the Kuiper Belt. We also discuss the caveats of using results from a Hapke reflectance model to derive specific surface particle properties.