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.     

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

Modification of Jupiter’s stratosphere three weeks after the 2009 impact

Infrared spectroscopy sensitive to thermal emission from Jupiter’s stratosphere reveals effects persisting 23 days after the impact of a body in late July 2009. Measurements obtained on 2009 August 11 UT at the impact latitude of 56°S (planetocentric), using the Goddard Heterodyne Instrument for Planetary Wind and Composition mounted on the NASA Infrared Telescope Facility, reveal increased ethane abundance and the effects of aerosol opacity. An interval of reduced thermal continuum emission at 11.744 μm is measured ∼60-80° towards planetary east of the impact site, estimated to be at 305° longitude (System III). Retrieved stratospheric ethane mole fraction in the near vicinity of the impact site is enhanced by up to ∼60% relative to quiescent regions at this latitude. Thermal continuum emission at the impact site, and somewhat west of it, is significantly enhanced in the same spectra that retrieve enhanced ethane mole fraction. Assuming that the enhanced continuum brightness near the impact site results from thermalized aerosol debris blocking contribution from the continuum formed in the upper troposphere and indicating the local temperature, then continuum emission by a haze layer can be approximated by an opaque surface inserted at the 45-60 mbar pressure level in the stratosphere in an unperturbed thermal profile, setting an upper limit on the pressure and therefore a lower limit on the altitude of the top of the impact debris at this time. The reduced continuum brightness east of the impact site can be modeled by an opaque surface near the cold tropopause, which is consistent with a lower altitude of ejecta/impactor-formed opacity. The physical extent of the observed region of reduced continuum implies a minimum average velocity of 21 m/s transporting material prograde (planetary east) from the impact.     

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.     

Long-term evolution of the aerosol debris cloud produced by the 2009 impact on Jupiter

We present a study of the long-term evolution of the cloud of aerosols produced in the atmosphere of Jupiter by the impact of an object on 19 July 2009 (Sánchez-Lavega, A. et al. [2010]. Astrophys. J. 715, L155-L159). The work is based on images obtained during 5 months from the impact to 31 December 2009 taken in visible continuum wavelengths and from 20 July 2009 to 28 May 2010 taken in near-infrared deep hydrogen-methane absorption bands at 2.1-2.3 μm. The impact cloud expanded zonally from ∼5000 km (July 19) to 225,000 km (29 October, about 180° in longitude), remaining meridionally localized within a latitude band from 53.5°S to 61.5°S planetographic latitude. During the first two months after its formation the site showed heterogeneous structure with 500-1000 km sized embedded spots. Later the reflectivity of the debris field became more homogeneous due to clump mergers. The cloud was mainly dispersed in longitude by the dominant zonal winds and their meridional shear, during the initial stages, localized motions may have been induced by thermal perturbation caused by the impact’s energy deposition. The tracking of individual spots within the impact cloud shows that the westward jet at 56.5°S latitude increases its eastward velocity with altitude above the tropopause by 5-10 m s⁻¹. The corresponding vertical wind shear is low, about 1 m s⁻¹ per scale height in agreement with previous thermal wind estimations. We found evidence for discrete localized meridional motions with speeds of 1-2 m s⁻¹. Two numerical models are used to simulate the observed cloud dispersion. One is a pure advection of the aerosols by the winds and their shears. The other uses the EPIC code, a nonlinear calculation of the evolution of the potential vorticity field generated by a heat pulse that simulates the impact. Both models reproduce the observed global structure of the cloud and the dominant zonal dispersion of the aerosols, but not the details of the cloud morphology. The reflectivity of the impact cloud decreased exponentially with a characteristic timescale of 15 days; we can explain this behavior with a radiative transfer model of the cloud optical depth coupled to an advection model of the cloud dispersion by the wind shears. The expected sedimentation time in the stratosphere (altitude levels 5-100 mbar) for the small aerosol particles forming the cloud is 45-200 days, thus aerosols were removed vertically over the long term following their zonal dispersion. No evidence of the cloud was detected 10 months after the impact.     

Latitudinal variations of HCN, HC3N, and C2N2 in Titan's stratosphere derived from Cassini CIRS data

Mid- and far-infrared spectra from the Composite InfraRed Spectrometer (CIRS) have been used to determine volume mixing ratios of nitriles in Titan's atmosphere. HCN, HC₃N, C₂H₂, and temperature were derived from 2.5 cm⁻¹ spectral resolution mid-IR mapping sequences taken during three flybys, which provide almost complete global coverage of Titan for latitudes south of 60° N. Three 0.5 cm⁻¹ spectral resolution far-IR observations were used to retrieve C₂N₂ and act as a check on the mid-IR results for HCN. Contribution functions peak at around 0.5-5 mbar for temperature and 0.1-10 mbar for the chemical species, well into the stratosphere. The retrieved mixing ratios of HCN, HC₃N, and C₂N₂ show a marked increase in abundance towards the north, whereas C₂H₂ remains relatively constant. Variations with longitude were much smaller and are consistent with high zonal wind speeds. For 90°-20° S the retrieved HCN abundance is fairly constant with a volume mixing ratio of around 1 × 10⁻⁷ at 3 mbar. More northerly latitudes indicate a steady increase, reaching around 4 × 10⁻⁷ at 60° N, where the data coverage stops. This variation is consistent with previous measurements and suggests subsidence over the northern (winter) pole at approximately 2 × 10⁻⁴ m s⁻¹. HC₃N displays a very sharp increase towards the north pole, where it has a mixing ratio of around 4 × 10⁻⁸ at 60° N at the 0.1-mbar level. The difference in gradient for the HCN and HC₃N latitude variations can be explained by HC₃N's much shorter photochemical lifetime, which prevents it from mixing with air at lower latitude. It is also consistent with a polar vortex which inhibits mixing of volatile rich air inside the vortex with that at lower latitudes. Only one observation was far enough north to detect significant amounts of C₂N₂, giving a value of around 9 × 10⁻¹⁰ at 50° N at the 3-mbar level.     

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.     

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.     

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.     

Characterising Saturn's vertical temperature structure from Cassini/CIRS

Thermal infrared spectra of Saturn from 10-1400 cm⁻¹ at 15 cm⁻¹ spectral resolution and a spatial resolution of 1°-2° latitude have been obtained by the Cassini Composite Infrared Spectrometer [Flasar, F.M., and 44 colleagues, 2004. Space Sci. Rev. 115, 169-297]. Many thousands of spectra, acquired over eighteen-months of observations, are analysed using an optimal estimation retrieval code [Irwin, P.G.J., Parrish, P., Fouchet, T., Calcutt, S.B., Taylor, F.W., Simon-Miller, A.A., Nixon, C.A., 2004. Icarus 172, 37-49] to retrieve the temperature structure and para-hydrogen distribution over Saturn's northern (winter) and southern (summer) hemispheres. The vertical temperature structure is analysed in detail to study seasonal asymmetries in the tropopause height (65-90 mbar), the location of the radiative-convective boundary (350-500 mbar), and the variation with latitude of a temperature knee (between 150 and 300 mbar) which was first observed in inversions of Voyager/IRIS spectra [Hanel, R., and 15 colleagues, 1981. Science 212, 192-200; Hanel, R., Conrath, B., Flasar, F.M., Kunde, V., Maguire, W., Pearl, J.C., Pirraglia, J., Samuelson, R., Cruikshank, D.P., Gautier, D., Gierasch, P.J., Horn, L., Ponnamperuma, C., 1982. Science 215, 544-548]. Uncertainties due to both the modelling of spectral absorptions (collision-induced absorption coefficients, tropospheric hazes, helium abundance) and the nature of our retrieval algorithm are quantified.Temperatures in the stratosphere near 1 mbar show a 25-30 K temperature difference between the north pole and south pole. This asymmetry becomes less pronounced with depth as the radiative time constant for the atmospheric response increases at deeper pressure levels. Hemispherically-symmetric small-scale temperature structures associated with zonal winds are superimposed onto the temperature asymmetry for pressures greater than 100 mbar. The para-hydrogen fraction in the 100-400 mbar range is greater than equilibrium predictions for the southern hemisphere and parts of the northern hemisphere, and less than equilibrium predictions polewards of 40° N.The temperature knee between 150-300 mbar is larger in the summer hemisphere than in the winter, smaller and higher at the equator, deeper and larger in the equatorial belts and small at the poles. Solar heating on tropospheric haze is proposed as a possible mechanism for this effect; the increased efficiency of ortho- to para-hydrogen conversion in the southern hemisphere is consistent with the presence of larger aerosols in the summer hemisphere, which we demonstrate to be qualitatively consistent with previous studies of Saturn's tropospheric aerosol distribution.     

The jovian anticyclone BA: I. Motions and interaction with the GRS from observations and non-linear simulations

A study of the dynamics of the second largest anticyclone in Jupiter, Oval BA, and its red colour change that occurred in late 2005 is presented in a three part study. The first part, this paper, deals with its long-term kinematical and dynamical behaviour monitored since its formation in 2000 to September 2008 using ground-based observations archived at the public International Outer Planet Watch (IOPW) database. The vortex changed its zonal drift velocity from 1.8 m s⁻¹ in the period 2000-2002 to 0.8 m s⁻¹ in 2002-2003, and to 2.5 m s⁻¹ since late 2003. It also migrated southwards by 1.0 ± 0.5° in latitude between 2000 and 2004, remaining afterwards at an almost fixed latitude position. During the period 2000-2007, the oval also changed its triangular-like shape to a more symmetrical one. No latitudinal change was found in the months before the development of a red annulus in its interior. The colour change took place in less than 5 months in 2005-2006 and no red colour feature was observed to have been present or entrained by BA months before the annulus development. After detailed examination of the four encounters between BA and GRS that took place during this 9 year period, we did not detect any noticeable change in its drift rate or in apparent structure associated with the encounters at cloud level. Also, the area of BA did not significantly change in this period. Additionally, we found that BA displays a long-term oscillation of ∼160 days in its longitude position with peak to peak amplitude of 1.2°. Numerical experiments using the global circulation model EPIC reproduce accurately the shape, connecting it to its latitude migration, and morphology of the oval and confirm that no strong interaction between BA and the GRS is possible at least in the current situation.     

Seasonal change on Saturn from Cassini/CIRS observations, 2004-2009

Five years of thermal infrared spectra from the Cassini Composite Infrared Spectrometer (CIRS) are analyzed to determine the response of Saturn’s atmosphere to seasonal changes in insolation. Hemispheric mapping sequences at 15.0 cm⁻¹ spectral resolution are used to retrieve the variation in the zonal mean temperatures in the stratosphere (0.5-5.0 mbar) and upper troposphere (75-800 mbar) between October 2004 (shortly after the summer solstice in the southern hemisphere) and July 2009 (shortly before the autumnal equinox).Saturn’s northern mid-latitudes show signs of dramatic warming in the stratosphere (by 6-10 K) as they emerge from ring-shadow into springtime conditions, whereas southern mid-latitudes show evidence for cooling (4-6 K). The 40-K asymmetry in stratospheric temperatures between northern and southern hemispheres (at 1 mbar) slowly decreased during the timespan of the observations. Tropospheric temperatures also show temporal variations but with a smaller range, consistent with the increasing radiative time constant of the atmospheric response with increasing pressure. The tropospheric response to the insolation changes shows the largest magnitude at the locations of the broad retrograde jets. Saturn’s warm south-polar stratospheric hood has cooled over the course of the mission, but remains present.Stratospheric temperatures are compared to a radiative climate model which accounts for the spatial distribution of the stratospheric coolants. The model successfully predicts the magnitude and morphology of the observed changes at most latitudes. However, the model fails at locations where strong dynamical perturbations dominate the temporal changes in the thermal field, such as the hot polar vortices and the equatorial semi-annual oscillation (Orton, G., and 27 colleagues [2008]. Nature 453, 196-198). Furthermore, observed temperatures in Saturn’s ring-shadowed regions are larger than predicted by all radiative-climate models to date due to the incomplete characterization of the dynamical response to the shadow. Finally, far-infrared CIRS spectra are used to demonstrate variability of the para-hydrogen distribution over the 5-year span of the dataset, which may be related to observed changes in Saturn’s tropospheric haze in the spring hemisphere.     

Phosphine on Jupiter and Saturn from Cassini/CIRS

The global distribution of phosphine (PH₃) on Jupiter and Saturn is derived using 2.5 cm⁻¹ spectral resolution Cassini/CIRS observations. We extend the preliminary PH₃ analyses on the gas giants [Irwin, P.G.J., and 6 colleagues, 2004. Icarus 172, 37-49; Fletcher, L.N., and 9 colleagues, 2007a. Icarus 188, 72-88] by (a) incorporating a wider range of Cassini/CIRS datasets and by considering a broader spectral range; (b) direct incorporation of thermal infrared opacities due to tropospheric aerosols and (c) using a common retrieval algorithm and spectroscopic line database to allow direct comparison between these two gas giants.The results suggest striking similarities between the tropospheric dynamics in the 100-1000 mbar regions of the giant planets: both demonstrate enhanced PH₃ at the equator, depletion over neighbouring equatorial belts and mid-latitude belt/zone structures. Saturn's polar PH₃ shows depletion within the hot cyclonic polar vortices. Jovian aerosol distributions are consistent with previous independent studies, and on Saturn we demonstrate that CIRS spectra are most consistent with a haze in the 100-400 mbar range with a mean optical depth of 0.1 at 10 μm. Unlike Jupiter, Saturn's tropospheric haze shows a hemispherical asymmetry, being more opaque in the southern summer hemisphere than in the north. Thermal-IR haze opacity is not enhanced at Saturn's equator as it is on Jupiter.Small-scale perturbations to the mean PH₃ abundance are discussed both in terms of a model of meridional overturning and parameterisation as eddy mixing. The large-scale structure of the PH₃ distributions is likely to be related to changes in the photochemical lifetimes and the shielding due to aerosol opacities. On Saturn, the enhanced summer opacity results in shielding and extended photochemical lifetimes for PH₃, permitting elevated PH₃ levels over Saturn's summer hemisphere.     

The evolving flow of Jupiter’s White Ovals and adjacent cyclones

We present results regarding the dynamical meteorology of Jupiter’s White Ovals at different points in their evolution. Starting from the era with three White Ovals FA, BC, and DE (Galileo), continuing to the post-merger epoch with only one Oval BA (Cassini), and finally to Oval BA’s current reddened state (New Horizons), we demonstrate that the dynamics of their flow have similarly evolved along with their appearance. In the Galileo epoch, Oval DE had an elliptical shape with peak zonal wind speeds of ∼90 m s⁻¹ in both its northern and southern peripheries. During the post-merger epoch, Oval BA’s shape was more triangular and less elliptical than Oval DE; in addition to widening in the north-south direction, its northern periphery was 20 m s⁻¹ slower, and its southern periphery was 20 m s⁻¹ faster than Oval DE’s flow during the Galileo era. Finally, in the New Horizons era, the reddened Oval BA had evolved back to a classical elliptical form. The northern periphery of Oval BA increased in speed by 20 m s⁻¹ from Cassini to New Horizons, ending up at a speed nearly identical to that of the northern periphery of Oval DE during Galileo. However, the peak speeds along the southern rim of the newly formed Oval BA were consistently faster than the corresponding speeds in Oval DE, and they increased still further between Cassini and New Horizons, ending up at ∼140-150 m s⁻¹. Relative vorticity maps of Oval BA reveal a cyclonic ring surrounding its outer periphery, similar to the ring present around the Great Red Spot. The cyclonic ring around Oval BA in 2007 appears to be moderately stronger than observed in 1997 and 2001, suggesting that this may be associated with the coloration of the vortex. The modest strengthening of the winds in Oval BA, the appearance of red aerosols, and the appearance of a turbulent, cyclonic feature to Oval BA’s northwest create a strong resemblance with the Great Red Spot from both a dynamical and morphological perspective.In addition to the White Ovals, we also measure the winds within two compact cyclonic regions, one in the Galileo data set and one in the Cassini data set. In the images, these cyclonic features appear turbulent and filamentary, but our wind field reveals that the flow manifests as a coherent high-speed collar surrounding relatively quiescent interiors. Our relative vorticity maps show that the vorticity likewise concentrates in a collar near the outermost periphery, unlike the White Ovals which have peak relative vorticity magnitudes near the center of the vortex. The cyclones contain several localized bright regions consistent with the characteristics of thunderstorms identified in other studies. Although less studied than their anticyclonic cousins, these cyclones may offer crucial insights into the planet’s cloud-level energetics and dynamical meteorology.     

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.     

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.     

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.     

Changes in Jupiter’s Great Red Spot (1979-2006) and Oval BA (2000-2006)

We analyze velocity fields of the Great Red Spot (GRS) and Oval BA that were previously extracted from Cassini, Galileo, and Hubble Space Telescope images (Asay-Davis, X.S., Marcus, P.S., Wong, M., de Pater, I. [2009]. Icarus 203, 164-188). Our analyses use reduced-parameter models in which the GRS, Oval BA, and surrounding zonal (east-west) flows are assumed to have piece-wise-constant potential vorticity (PV), but with finite-sized transition regions between the pieces of constant PV rather than sharp steps. The shapes of the regions of constant PV are computed such that the flow is a steady, equilibrium solution of the 2D quasigeostrophic equations when viewed in a frame translating uniformly in the east-west direction. All parameter values of the models, including the magnitudes of the PV, areas of the regions with constant PV, locations of the transition regions, widths of the transition regions, and the value of the Rossby deformation radius, are found with a genetic algorithm such that the velocity produced by the equilibrium solution is a “best-fit” to the observed velocity fields. A Monte Carlo method is used to estimate the uncertainties in the best-fit parameter values.The best-fit results show that there were significant changes (greater than the uncertainties) in the PV of the GRS between Galileo in 1996 and Hubble in 2006. In particular, the shape of the PV anomaly of the GRS became rounder, and the area of the PV anomaly of the GRS decreased by 18%, although the magnitudes of PV in the anomaly remained constant. In contrast, neither the area nor the magnitude of the PV anomaly of the Oval BA changed from 2000, when its cloud cover was white, to 2006, when its cloud cover was red. The best-fit results also show that the areas of the PV anomalies of the GRS and of the Oval BA are smaller than the areas of their corresponding cloud covers at all times. Using the best-fit values of the Rossby deformation radius, we show that the Brunt-Väisälä frequency is 15% larger at 33°S than at 23°S. As expected (Marcus, 1993), the best-fit results show that the PV of the zonal flow has “jumps” at the latitudes of the maxima of the eastward-going jet streams. However, a surprising result is that a large “jump” in the PV of the zonal flow occurs at the location of a maximum of the westward going jet stream neighboring the GRS. Another surprise is that the jumps in the PV of the zonal flow do not all have the same sign, which implies that there is not a monotonic “staircase” of zonal PV from north to south as was anticipated ( [Marcus, 1993] and [McIntyre, 2008]).     

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.     

Jupiter’s shrinking Great Red Spot and steady Oval BA: Velocity measurements with the ‘Advection Corrected Correlation Image Velocimetry’ automated cloud-tracking method

We show that between 1996 and 2006, the area circumscribed by the high-speed collar of the Great Red Spot (GRS) shrunk by 15%, while the peak velocities within its collar remained constant. This shrinkage indicates a dynamical change in the GRS because the region circumscribed by the collar is nearly coincident with the location of the potential vorticity anomaly of the GRS. It was previously observed that the area of the clouds associated with the GRS has been shrinking. However, the cloud cover of the GRS is not coincident with the location of its potential vorticity anomaly or any other of its known dynamical features. We show that the peak velocities of the Oval BA were nearly the same in 2000, when the Oval was white, and in 2006, when it was red, as were all of the other features of the two velocities fields. To measure temporal changes in the GRS and Oval, we extracted velocities from images taken with Galileo, Cassini, and the Hubble Space Telescope using a new iterative method called Advection Corrected Correlation Image Velocimetry (ACCIV). ACCIV finds correlations over image pairs with 10-h time separations when other automated velocity-extraction methods are limited to time separations of 2 h or less. Typically, ACCIV velocities produced from images separated by 10 h had errors that are 3-6 times smaller than similar velocities extracted from images separated by 2 h or less. ACCIV produces velocity fields containing hundreds of thousands of independent correlation vectors (tie-points). Dense velocity fields are needed to locate the loci of peak velocities and other features.