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.     

Compressible convection in the deep atmospheres of giant planets

Fast rotating giant planets such as Jupiter and Saturn possess alternate prograde and retrograde zonal winds which are stable over long periods of time. We consider a compressible model of convection in a spherical shell with rapid rotation, using the anelastic approximation, to explore the parameter range for which such zonal flows can be produced.We consider models with a large variation in density across the layer. Our models are based only on the molecular H/He region above the metallic hydrogen transition at about 2 Mbar, and we do not include the hydromagnetic effects which may be important if the electrical conductivity is significant. We find that the convective velocities are significantly higher in the low density regions of the shell, but the zonal flow is almost independent of the z-coordinate parallel to the rotation axis. We analyse how this behaviour is consistent with the Proudman-Taylor theorem.We find that deep prograde zonal flow near the equator is a very robust feature of our models. Prograde and retrograde jets alternating in latitude can occur inside the tangent cylinder in compressible as well as Boussinesq models, particularly at lower Prandtl numbers. However, the zonal jets inside the tangent cylinder are suppressed if a no-slip condition is imposed at the inner boundary. This suggests that deep high latitude jets may be suppressed if there is significant magnetic dissipation.Our compressible calculations include the viscous dissipation in the entropy equation, and we find this is comparable to, and in some cases exceeds, the total heat flux emerging from the surface. For numerical reasons, these simulations cannot reach the extremely low Ekman number found in giant planets, and they necessarily also have a much larger heat flux than planets. We therefore discuss how our results might scale down to give solutions with lower dissipation and lower heat flux.     

Scaling laws for convection and jet speeds in the giant planets

Three dimensional studies of convection in deep spherical shells have been used to test the hypothesis that the strong jet streams on Jupiter, Saturn, Uranus, and Neptune result from convection throughout the molecular envelopes. Due to computational limitations, these simulations must be performed at parameter settings far from jovian values and generally adopt heat fluxes 5-10 orders of magnitude larger than the planetary values. Several numerical investigations have identified trends for how the mean jet speed varies with heat flux and viscosity in these models, but no previous theories have been advanced to explain these trends. Here, we show using simple arguments that if convective release of potential energy pumps the jets and viscosity damps them, the mean jet speeds split into two regimes. When the convection is weakly nonlinear, the equilibrated jet speeds should scale approximately with F/ν, where F is the convective heat flux and ν is the viscosity. When the convection is strongly nonlinear, the jet speeds are faster and should scale approximately as (F/ν)^1/2. We demonstrate how this regime shift can naturally result from a shift in the behavior of the jet-pumping efficiency with heat flux and viscosity. Moreover, both Boussinesq and anelastic simulations hint at the existence of a third regime where, at sufficiently high heat fluxes or sufficiently small viscosities, the jet speed becomes independent of the viscosity. We show based on mixing-length estimates that if such a regime exists, mean jet speeds should scale as heat flux to the 1/4 power. Our scalings provide a good match to the mean jet speeds obtained in previous Boussinesq and anelastic, three-dimensional simulations of convection within giant planets over a broad range of parameters. When extrapolated to the real heat fluxes, these scalings suggest that the mass-weighted jet speeds in the molecular envelopes of the giant planets are much weaker—by an order of magnitude or more—than the speeds measured at cloud level.     

Generation of equatorial jets by large-scale latent heating on the giant planets

Three-dimensional numerical simulations show that large-scale latent heating resulting from condensation of water vapor can produce multiple zonal jets similar to those on the gas giants (Jupiter and Saturn) and ice giants (Uranus and Neptune). For plausible water abundances (3-5 times solar on Jupiter/Saturn and 30 times solar on Uranus/Neptune), our simulations produce ∼20 zonal jets for Jupiter and Saturn and 3 zonal jets on Uranus and Neptune, similar to the number of jets observed on these planets. Moreover, these Jupiter/Saturn cases produce equatorial superrotation whereas the Uranus/Neptune cases produce equatorial subrotation, consistent with the observed equatorial-jet direction on these planets. Sensitivity tests show that water abundance, planetary rotation rate, and planetary radius are all controlling factors, with water playing the most important role; modest water abundances, large planetary radii, and fast rotation rates favor equatorial superrotation, whereas large water abundances favor equatorial subrotation regardless of the planetary radius and rotation rate. Given the larger radii, faster rotation rates, and probable lower water abundances of Jupiter and Saturn relative to Uranus and Neptune, our simulations therefore provide a possible mechanism for the existence of equatorial superrotation on Jupiter and Saturn and the lack of superrotation on Uranus and Neptune. Nevertheless, Saturn poses a possible difficulty, as our simulations were unable to explain the unusually high speed (∼) of that planet’s superrotating jet. The zonal jets in our simulations exhibit modest violations of the barotropic and Charney-Stern stability criteria. Overall, our simulations, while idealized, support the idea that latent heating plays an important role in generating the jets on the giant planets.     

New upper limits for hydrogen halides on Saturn derived from Cassini-CIRS data

Far infrared spectra (10-600 cm⁻¹) from Cassini's Composite InfraRed Spectrometer (CIRS) were used to determine improved upper limits of hydrogen halides HF, HCl, HBr, and HI in Saturn's atmosphere. Three observations, comprising a total of 3088 spectra, gave 3σ upper limits on HF, HCl, HBr, and HI volume mole fractions of 8.0×¹⁰⁻¹², 6.7×¹⁰⁻¹¹, 1.3×¹⁰⁻¹⁰, and 1.4×¹⁰⁻⁹, respectively, at the 500 mbar pressure level. These upper limits confirm sub-solar abundances of halide species for HF, HCl, and HBr in Saturn's upper atmosphere—consistent with predictions from thermochemical models and influx of material from meteoroids. Our upper limit for HCl is 16 times lower than the tentative detection at 1.1×¹⁰⁻⁹ reported by Weisstein and Serabyn [Weisstein, E.W., Serabyn, E., 1996. Icarus 123, 23-36]. These observations are not sensitive to the deep halide abundance, which is expected to be enriched relative to the solar composition.     

Turbulent convection in rapidly rotating spherical shells: A model for equatorial and high latitude jets on Jupiter and Saturn

The origin of zonal jets on the jovian planets has long been a topic of scientific debate. In this paper we show that deep convection in a spherical shell can generate zonal flow comparable to that observed on Jupiter and Saturn, including a broad prograde equatorial jet and multiple alternating jets at higher latitudes. We present fully turbulent, 3D spherical numerical simulations of rapidly rotating convection with different spherical shell geometries. The resulting global flow fields tend to be segregated into three regions (north, equatorial, and south), bounded by the tangent cylinder that circumscribes the inner boundary equator. In all of our simulations a strong prograde equatorial jet forms outside the tangent cylinder, whereas multiple jets form in the northern and southern hemispheres, inside the tangent cylinder. The jet scaling of our numerical models and of Jupiter and Saturn is consistent with the theory of geostrophic turbulence, which we extend to include the effect of spherical shell geometry. Zonal flow in a spherical shell is distinguished from that in a full sphere or a shallow layer by the effect of the tangent cylinder, which marks a reversal in the sign of the planetary β-parameter and a jump in the Rhines length. This jump is manifest in the numerical simulations as a sharp equatorward increase in jet widths—a transition that is also observed on Jupiter and Saturn. The location of this transition gives an estimate of the depth of zonal flow, which seems to be consistent with current models of the jovian and saturnian interiors.     

Methane and its isotopologues on Saturn from Cassini/CIRS observations

High spectral resolution observations from the Cassini Composite Infrared Spectrometer [Flasar, F.M., and 44 colleagues, 2004. Space Sci. Rev. 115, 169-297] are analysed to derive new estimates for the mole fractions of CH₄, CH₃D and ¹³CH₄ of (4.7±0.2)×¹⁰⁻³, (3.0±0.2)×¹⁰⁻⁷ and (5.1±0.2)×¹⁰⁻⁵ respectively. The mole fractions show no hemispherical asymmetries or latitudinal variability. The analysis combines data from the far-IR methane rotational lines and the mid-IR features of methane and its isotopologues, using both the correlated-k retrieval algorithm of Irwin et al. [Irwin, P., and 9 colleagues, 2008. J. Quant. Spectrosc. Radiat. Trans. 109, 1136-1150] and a line-by-line approach to evaluate the reliability of the retrieved quantities. C/H was found to be enhanced by 10.9±0.5 times the solar composition of Grevesse et al. [Grevesse, N., Asplund, M., Sauval, A., 2007. Space Sci. Rev. 130 (1), 105-114], 2.25±0.55 times larger than the enrichment on Jupiter, and supporting the increasing fractional core mass with distance from the Sun predicted by the core accretion model of planetary formation. A comparison of the jovian and saturnian C/N, C/S and C/P ratios suggests different reservoirs of the trapped volatiles in a primordial solar nebula whose composition varies with distance from the Sun. This is supported by our derived D/H ratio in methane of (1.6±0.2)×¹⁰⁻⁵, which appears to be smaller than the jovian value of Lellouch et al. [Lellouch, E., Bézard, B., Fouchet, T., Feuchtgruber, H., Encrenaz, T., de Graauw, T., 2001. Astron. Astrophys. 370, 610-622]. Mid-IR emission features provided an estimate of , which is consistent with both the terrestrial ratio and jovian ratio, suggesting that carbon was accreted from a shared reservoir for all of the planets.     

Interaction of moist convection with zonal jets on Jupiter and Saturn

Observations suggest that moist convection plays an important role in the large-scale dynamics of Jupiter's and Saturn's atmospheres. Here we use a reduced-gravity quasigeostrophic model, with a parameterization of moist convection that is based on observations, to study the interaction between moist convection and zonal jets on Jupiter and Saturn. Stable jets with approximately the same width and strength as observations are generated in the model. The observed zonal jets violate the barotropic stability criterion but the modeled jets do so only if the flow in the deep underlying layer is westward. The model results suggest that a length scale and a velocity scale associated with moist convection control the width and strength of the jets. The length scale and velocity scale offer a possible explanation of why the jets of Saturn are stronger and wider than those of Jupiter.     

Jupiter and Saturn rotation periods

Anderson and Schubert [2007. Saturn's Gravitational field, internal rotation, and interior structure. Science 317, 1384-1387 (paper I)] proposed that Saturn's rotation period can be ascertained by minimizing the dynamic heights of the 100 mbar isosurface with respect to the geoid; they derived a rotation period of 10 h 32 m 35 s. We investigate the same approach for Jupiter to see if the Jovian rotation period is predicted by minimizing the dynamical heights of its isobaric (1 bar pressure level) surface using zonal wind data. A rotation period of 9 h 54 m 29.7 s is found. Further, we investigate the minimization method by fitting Pioneer and Voyager occultation radii for both Jupiter and Saturn. Rotation periods of 9 h 55 m 30 s and 10 h 32 m 35 s are found to minimize the dynamical heights for Jupiter and Saturn, respectively. Though there is no dynamical principle requiring the minimization of the dynamical heights of an isobaric surface, the successful application of the method to Jupiter lends support to its relevance for Saturn.We derive Jupiter and Saturn rotation periods using equilibrium theory to explain the difference between equatorial and polar radii. Rotation periods of 9 h 55 m 20 s and 10 h 31 m 49 s are found for Jupiter and Saturn, respectively. We show that both Jupiter's and Saturn's shapes can be derived using solid-body rotation, suggesting that zonal winds have a minor effect on the planetary shape for both planets.The agreement in the values of Saturn's rotation period predicted by the different approaches supports the conclusion that the planet's period of rotation is about 10 h 32 m.     

UV and IR auroral emission model for the outer planets: Jupiter and Saturn comparison

Planetary aurora display the dynamic behavior of the plasma gas surrounding a planet. The outer planetary aurora are most often observed in the ultraviolet (UV) and the infrared (IR) wavelengths. How the emissions in these different wavelengths are connected with the background physical conditions are not yet well understood. Here we investigate the sensitivity of UV and IR emissions to the incident precipitating auroral electrons and the background atmospheric temperature, and compare the results obtained for Jupiter and Saturn. We develop a model which estimates UV and IR emission rates accounting for UV absorption by hydrocarbons, ion chemistry, and non-LTE effects. Parameterization equations are applied to estimate the ionization and excitation profiles in the H₂ atmosphere caused by auroral electron precipitation. The dependences of UV and IR emissions on electron flux are found to be similar at Jupiter and Saturn. However, the dependences of the emissions on electron energy are different at the two planets, especially for low energy (<10 keV) electrons; the UV and IR emissions both decrease with decreasing electron energy, but this effect in the IR is less at Saturn than at Jupiter. The temperature sensitivity of the IR emission is also greater at Saturn than at Jupiter. These dependences are interpreted as results of non-LTE effects on the atmospheric temperature and density profiles. The different dependences of the UV and IR emissions on temperature and electron energy at Saturn may explain the different appearance of polar emissions observed at UV and IR wavelengths, and the differences from those observed at Jupiter. These results lead to the prediction that the differences between the IR and UV aurora at Saturn may be more significant than those at Jupiter. We consider in particular the occurrence of bright polar infrared emissions at Saturn and quantitatively estimate the conditions for such IR-only emissions to appear.     

The deep wind structure of the giant planets: Results from an anelastic general circulation model

The giant gas planets have hot convective interiors, and therefore a common assumption is that these deep atmospheres are close to a barotropic state. Here we show using a new anelastic general circulation model that baroclinic vorticity contributions are not negligible, and drive the system away from an isentropic and therefore barotropic state. The motion is still aligned with the direction of the axis of rotation as in a barotropic rotating fluid, but the wind structure has a vertical shear with stronger winds in the atmosphere than in the interior. This shear is associated with baroclinic compressibility effects. Most previous convection models of giant planets have used the Boussinesq approximation, which assumes the density is constant in depth; however, Jupiter's actual density varies by four orders of magnitude through its deep molecular envelope. We therefore developed a new general circulation model (based on the MITgcm) that is anelastic and thereby incorporates this density variation. The model's geometry is a full 3D sphere down to a small inner core. It is nonhydrostatic, uses an equation of state suitable for hydrogen-helium mixtures (SCVH), and is driven by an internal heating profile. We demonstrate the effect of compressibility by comparing anelastic and Boussinesq cases. The simulations develop a mean state that is geostrophic and hydrostatic including the often neglected, but significant, vertical Coriolis contribution. This leads to modification of the standard thermal wind relation for a deep compressible atmosphere. The interior flow organizes in large cyclonically rotating columnar eddies parallel to the rotation axis, which drive upgradient angular momentum eddy fluxes, generating the observed equatorial superrotation. Heat fluxes align with the axis of rotation, and provide a mechanism for the transport of heat poleward, which can cause the observed flat meridional emission. We address the issue of over-forcing which is common in such convection models and analyze the dependence of our results on this; showing that the vertical wind structure is not very sensitive to the Rayleigh number. We also study the effect of rotation, showing how the transition from a rapidly to a slowly rotating system affects the dynamics.     

The hydrogen ortho-to-para ratio in the stratospheres of the giant planets

Observations of the H₂ S(0) and S(1) quadrupole lines in the four giant planets by the short-wavelength spectrometer of the Infrared Space Observatory are analyzed. These lines probe pressure levels located between 10 and 1 mbar and allow us to determine the stratospheric hydrogen para fraction for the first time. In Jupiter and Saturn, the stratospheric para fraction is close to its tropopause value. In the stratosphere of these planets as well as in Neptune’s, the para fraction presents a significant departure from thermodynamic equilibrium. This situation results from a lagged conversion between the ortho and the para states as molecular hydrogen is transported upward under the influence of turbulent eddy diffusion. In contrast, the uranian stratosphere lies close to thermodynamic equilibrium. The magnitude of the departure from thermodynamic equilibrium appears to be anti-correlated with the amount of stratospheric aerosols. To validate this assumption, we estimate the hydrogen equilibration time with a one-dimensional diffusion model for different conversion processes in the gas phase or on aerosols. The comparison between our results and the tropospheric estimates from Conrath et al. (1998, Icarus,135, 501-517) shows that paramagnetic conversion on aerosols matches the estimated tropospheric and stratospheric relaxation times in the four giant planets. In contrast, paramagnetic conversion in the gas phase can only explain the relaxation times measured in Jupiter and Saturn atmospheres. This situation provides quantitative evidence for an equilibration mechanism dominated by conversion on aerosols.     

Stability of jets on Jupiter and Saturn

The linear stability of a zonal jet that decays with depth is investigated under the assumption that the thermal stratification is very small. A westerly cosine jet is found to be more stable than it is in a thin fluid shell with two-dimensional flow. This is in agreement with observations of Jupiter and Saturn, where jet curvature exceeds the barotropic stability criterion. This result constitutes an alternative hypothesis to that of Ingersoll and Pollard [Icarus 52 (1982) 62], who showed that deep jets extending through the interior are also more stable than thin shell jets. The flow regime assumed in the present work requires that a small stratification can exist and persist even in the presence of horizontal temperature gradients. Further work will be needed to test whether this is realistic.     

A three-dimensional model of moist convection for the giant planets II: Saturn's water and ammonia moist convective storms

Moist convective storms constitute a key aspect in the global energy budget of the atmospheres of the giant planets. Among them, Saturn is known to develop the largest scale convective storms in the Solar System, the Great White Spots (GWS) which occur rarely and have been detected once every 30 years approximately. On the average, Saturn seems to show much less convective storms than Jupiter with smaller size and reduced frequency and intensity. Here we present detailed simulations of the onset and development of storms at the Equator and mid-latitudes of Saturn. These are the regions where most of the recent convective activity of the planet has been observed. We use a 3D anelastic model with parameterized microphysics (Hueso and Sánchez-Lavega, 2001, Icarus 151, 257) studying the onset and evolution of water and ammonia moist convective storms up to sizes of a few hundred km. Water storms, while more difficult to initiate than in Jupiter, can be very energetic, arriving to the 150 mbar level and developing vertical velocities on the order of 150 m s⁻¹. Ammonia storms develop easier but with a much smaller intensity unless very large abundances of ammonia (10 times solar) are present in Saturn's atmosphere. The Coriolis forces play a major role in the morphology and properties of water based storms.     

Latitudinal variation of Saturn photochemistry deduced from spatially-resolved ultraviolet spectra

We obtained spatially-resolved ultraviolet spectra of Saturn in 1994 with the Faint Object Spectrometer and Goddard High Resolution Spectrograph of the Hubble Space Telescope. We observed four areas on the planet at 15° N, 33° S, 41° S, and 52° S, with a field-of-view of less than 2 × 2 arcsec², compared to the 16-arcsec planet diameter. The wavelength range, 1550-2300 Å, encompasses absorption from major hydrocarbons (C₂H₆, C₂H₄, C₂H₂, CH₃C₂H, C₄H₂) and water. We find global hydrocarbon abundances and a C₂H₂ vertical distribution compatible with infrared observations, in contrast with previous analyses of ultraviolet spectra. The stratospheric haze opacity decreases from polar region to the equator. Saturn mid-latitudes are photochemically distinct from the rest of the planet. At 33° S, the spectrum requires either (1) a distinctly different C₂H₂ vertical distribution or (2) a locally enhanced water abundance. At 41° S, the hydrocarbon abundance exhibits a local minimum, within a global trend of increasing abundance from equator to pole. This global trend may result from an increased abundance of short-lived hydrocarbons such as C₄H₂. Photochemical models predict a depletion of hydrocarbon molecules in the presence of stratospheric water [Moses et al., 2000. Icarus 143, 166-202]. These results are consistent with a localized influx of water, in the form of high charge to mass ratio particles, flowing into Saturn's atmosphere at latitudes magnetically linked to the rings.     

Deep jets on gas-giant planets

Three-dimensional numerical simulations of the atmospheric flow on giant planets using the primitive equations show that shallow thermal forcing confined to pressures near the cloud tops can produce deep zonal winds from the tropopause all the way down to the bottom of the atmosphere. These deep winds can attain speeds comparable to the zonal jet speeds within the shallow, forced layer; they are pumped by Coriolis acceleration acting on a deep meridional circulation driven by the shallow-layer eddies. In the forced layer, the flow reaches an approximate steady state where east-west eddy accelerations balance Coriolis accelerations acting on the meridional flow. Under Jupiter-like conditions, our simulations produce 25 to 30 zonal jets, similar to the number of jets observed on Jupiter and Saturn. The simulated jet widths correspond to the Rhines scale; this suggests that, despite the three-dimensional nature of the dynamics, the baroclinic eddies energize a quasi-two-dimensional inverse cascade modified by the β effect (where β is the gradient of the Coriolis parameter). In agreement with Jupiter, the jets can violate the barotropic and Charney-Stern stability criteria, achieving curvatures ^∂2u/∂y² of the zonal wind u with northward distance y up to 2β. The simulations exhibit a tendency toward neutral stability with respect to Arnol'd's second stability theorem in the upper troposphere, as has been suggested for Jupiter, although deviations from neutrality exist. When the temperature varies strongly with latitude near the equator, our simulations can also reproduce the stable equatorial superrotation with wind speeds greater than . Diagnostics show that barotropic eddies at low latitudes drive the equatorial superrotation. The simulations also broadly explain the distribution of jet-pumping eddies observed on Jupiter and Saturn. While idealized, these simulations therefore capture many aspects of the cloud-level flows on Jupiter and Saturn.     

Retrievals of atmospheric variables on the gas giants from ground-based mid-infrared imaging

Thermal-infrared imaging of Jupiter and Saturn using the NASA/IRTF and Subaru observatories are quantitatively analyzed to assess the capabilities for reproducing and extending the zonal mean atmospheric results of the Cassini/CIRS experiment. We describe the development of a robust, systematic and reproducible approach to the acquisition and reduction of planetary images in the mid-infrared (7-25 μm), and perform an adaptation and validation of the optimal estimation, correlated-k retrieval algorithm described by Irwin et al. [Irwin, P., Teanby, N., de Kok, R., Fletcher, L., Howett, C., Tsang, C., Wilson, C., Calcutt, S., Nixon, C., Parrish, P., 2008. J. Quant. Spectrosc. Radiat. Trans. 109 (6), 1136-1150] for channel-integrated radiances. Synthetic spectral analyses and a comparison to Cassini results are used to verify our abilities to retrieve temperatures, haze opacities and gaseous abundances from filtered imaging. We find that ground-based imaging with a sufficiently high spatial resolution is able to reproduce the three-dimensional temperature and para-H₂ fields measured by spacecraft visiting Jupiter and Saturn, allowing us to investigate vertical wind shear, pressure and, with measured cloud-top winds, Ertel potential vorticity on potential temperature surfaces. Furthermore, by scaling vertical profiles of NH₃, PH₃, haze opacity and hydrocarbons as free parameters during thermal retrievals, we can produce meridional results comparable with CIRS spectroscopic investigations. This paper demonstrates that mid-IR imaging instruments operating at ground-based observatories have access to several dynamical and chemical diagnostics of the atmospheric state of the gas giants, offering the prospect for quantitative studies over much longer baselines and often covering much wider areas than is possible from spaceborne platforms.     

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.     

A 2D study of the effects of the size of a solid core on the equatorial flow in giant planets

Zonal winds simulated in two-dimensional computer models of turbulent convection in the equatorial plane of giant planets have greater surface amplitudes for cases with smaller solid cores, and therefore larger buoyancy driving, all other properties being equal. This differential rotation in radius is maintained by the convergence of angular momentum flux, which occurs because of the convective flow that develops due to the effects of planetary rotation and density stratification. The superposition of the convective flow and the stronger zonal flow produces wave-like, instead of cellular convection.     

Emergence of polar-jet polygons from jet instabilities in a Saturn model

Voyager flybys of Saturn in 1980-1981 revealed a circumpolar wave at ≈78° north planetographic latitude. The feature had a dominant wavenumber 6 mode, and has been termed the Hexagon from its geometric appearance in polar-projected mosaics. It was also noted for being stationary with respect to Saturn’s Kilometric Radiation (SKR) rotation rate. The Hexagon has persisted for over 30 years since the Voyager observations until now. It has been observed from ground based telescopes, Hubble Space Telescope and multiple instruments onboard Cassini in orbit around Saturn. Measurements of cloud motions in the region reveal the presence of a jet stream whose path closely follows the Hexagon’s outline. Why the jet stream takes the characteristic six-sided shape and how it is stably maintained across multiple saturnian seasons are yet to be explained. We present numerical simulations of the 78.3°N jet using the Explicit Planetary Isentropic-Coordinate (EPIC) model and demonstrate that a stable hexagonal structure can emerge without forcing when dynamic instabilities in the zonal jet nonlinearly equilibrate. For a given amplitude of the jet, the dominant zonal wavenumber is most strongly dependent on the peak curvature of the jet, i.e., the second north-south spatial derivative of the zonal wind profile at the center of the jet. The stable polygonal shape of the jet in our simulations is formed by a vortex street with cyclonic and anticyclonic vortices lining up towards the polar and equatorial side of the jet, respectively. Our result is analogous to laboratory experiments of fluid motions in rotating tanks that develop polygonal flows out of vortex streets. However, our results also show that a vortex street model of the Hexagon cannot reproduce the observed propagation speed unless the zonal jet’s speed is modified beyond the uncertainties in the observed zonal wind speed, which suggests that a vortex street model of the Hexagon and the observed zonal wind profile may not be mutually compatible.