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.     

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.     

The effects of vigorous mixing in a convective model of zonal flow on the ice giants

Previous studies have used models of three-dimensional (3D) Boussinesq convection in a rotating spherical shell to explain the zonal flows on the gas giants, Jupiter and Saturn. In this paper we demonstrate that this approach can also generate flow patterns similar to those observed on the ice giants, Uranus and Neptune. The equatorial jets of Uranus and Neptune are often assumed to result from baroclinic cloud layer processes and have been simulated with shallow layer models. Here we show that vigorous, 3D convection in a spherical shell can produce the retrograde (westward) equatorial flows that occur on the ice giants as well as the prograde (eastward) equatorial flows of the gas giants. In our models, the direction of the equatorial jet depends on the ratio of buoyancy to Coriolis forces in the system. In cases where Coriolis forces dominate buoyancy, cylindrical Reynolds stresses drive prograde equatorial jets. However, as buoyancy forces approach and exceed Coriolis forces, the cylindrical nature of the flow is lost and 3D mixing homogenizes the fluid's angular momentum; the equatorial jet reverses direction, while strong prograde jets form in the polar regions. Although the results suggest that conditions involving strong atmospheric mixing are responsible for generating the zonal flows on the ice giants, our present models require roughly 100 and 10 times the internal heat fluxes observed on Uranus and Neptune, respectively.     

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.     

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.     

Constraints on deep-seated zonal winds inside Jupiter and Saturn

The atmospheres of Jupiter and Saturn exhibit strong and stable zonal winds. How deep the winds penetrate unabated into each planet is unknown. Our investigation favors shallow winds. It consists of two parts. The first part makes use of an Ohmic constraint; Ohmic dissipation associated with the planet's magnetic field cannot exceed the planet's net luminosity. Application to Jupiter (J) and Saturn (S) shows that the observed zonal winds cannot penetrate below a depth at which the electrical conductivity is about six orders of magnitude smaller than its value at the molecular-metallic transition. Measured values of the electrical conductivity of molecular hydrogen yield radii of maximum penetration of 0.96RJ and 0.86RS, with uncertainties of a few percent of R. At these radii, the magnetic Reynolds number based on the zonal wind velocity and the scale height of the magnetic diffusivity is of order unity. These limits are insensitive to difficulties in modeling turbulent convection. They permit complete penetration along cylinders of the equatorial jets observed in the atmospheres of Jupiter and Saturn. The second part investigates how deep the observed zonal winds actually do penetrate. As it applies heuristic models of turbulent convection, its conclusions must be regarded as tentative. Truncation of the winds in the planet's convective envelope would involve breaking the Taylor-Proudman constraint on cylindrical flow. This would require a suitable nonpotential acceleration which none of the obvious candidates appears able to provide. Accelerations arising from entropy gradients, magnetic stresses, and Reynolds stresses appear to be much too weak. These considerations suggest that strong zonal winds are confined to shallow, stably stratified layers, with equatorial jets being the possible exception.     

Thermal structure of Saturn from Voyager infrared measurements: Implications for atmospheric dynamics

Saturn atmospheric temperatures at the 150-mbar level retrieved from Voyager IRIS measurements indicate the presence of small-scale meridional gradients which are approximately symmetric with respect to the equator, but are superposed on a large-scale hemispheric thermal asymmetry. Under the assumption that the retrieved values at this atmospheric level represent kinetic temperatures on a constant pressure surface, it is suggested that the small-scale structure is produced by a meriodional circulation associated with the dissipative decay of the zonal winds with height, while the hemispheric asymmetry represents a thermal response to the seasonally varying insolation. The small-scale gradients are correlated with zonal winds derived from Voyager images at mid and high latitudes through the thermal wind relation; the calculated thermal wind shears suggest a decay with height of the jet system toward a state of uniform eatward flow. The existence of the approximately symmetric zonal winds and associated temperature gradients in the presence of a large-scale seasonal thermal response suggests that the jet system is driven at depths substantially below the levels where seasonally modulated insolation is important ($\text{p?0.5}\text{bar}$).     

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

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

Instability of Zonal Flows in Rotating Spherical Shells: An Application to Jupiter

Measurements from the Galileo probe suggest that the zonal winds are deep rooted. Jupiter's high rotation rate makes it likely that the whole outer molecular H/He layer is involved in these long-lived jet flows. Assuming that the primary flows are geostrophic, and that the banded surface structure stretches right through the molecular H/He layer, we examine the conditions for such flows to be stable. As a first step, the linear stability of some prescribed banded zonal flows in a rotating spherical shell is explored. Incompressibility is assumed for simplicity, and the boundary condition is stress-free. We compare solutions for two aspect ratios, appropriate for the molecular H/He layers of Jupiter and Saturn, and two Ekman numbers (E=10⁻² and E=10⁻⁴). Convective and shear flow instabilities compete in our system. The convective instabilities are of the well-known columnar structure. Shear flow instabilities for the larger Ekman number are similar to the Taylor-Couette instability in rotating annuli. At the lower Ekman number, shear flow instabilities adopt a geostrophic character, assuming the form of rotating columns, similar to the convective instabilities. While the convective instability always sets in outside the tangent cylinder, shear instability can become unstable inside the tangent cylinder. If even a weak zonal flow is present inside the tangent cylinder, the flow is unstable to shear instability. This offers an explanation why the jovian zonal jet structure is much weaker at the higher latitudes that correspond to locations inside the tangent cylinder.     

Discrete cloud activity in Saturn’s equator during 1995, 1996 and 1997

A regular extensive CCD imaging of Saturn allowed us to analyze the discrete cloud activity in the Equatorial Zone from 1995 to 1997. The large-scale storm observed in 1994 at +10° (Sanchez-Lavega et al., 1994, Sanchez-Lavega et al., 1996) was rediscovered in 1995, reaching a lifetime >1 year. Its slow motion characterized by a zonal velocity difference of −150 ms⁻¹ relative to background flow is confirmed. Our red and near infrared observations showed a strong increase of white cloud activity in the southern Equatorial Zone (latitude −13.5°) during 1996, declining later on during 1997. Cloud tracking of two prominent plumes and other features allowed us to measure zonal wind velocities and to compare them to the Voyager zonal flow velocity profile. We note that in general the 1995–1997 features have velocities lower than those measured with the Voyagers. Altitude differences in the clouds and hence different zonal velocities, or real changes in the zonal jet as a consequence of Saturn’s insolation cycle and ring-shadowing, can be the reason for such differences.     

Zonal winds near Venus' cloud top level: A model study of the interaction between the zonal mean circulation and the semidiurnal tide

A primitive equation wave-mean flow interaction model, designed by J. R. Holton and used originally to study Earth's middle atmosphere, has been adapted to Venus in order to clarify our understanding of the interaction between the semidiurnal tide and the thermally driven mean meridional circulation near the cloud top level. With or without the tide the model produces midlatitude jets whose structure is insensitive to vertical shear of the background angular velocity above and below the cloud top level, but it is sensitive to background angular velocity at the cloud top level. When this background angular velocity is close to that of Venus, the latitudes and speeds of these jets are similar to the latitudes and speeds of jets at the Venus cloud top level as inferred from observed temperatures and the cyclostrophic balance condition. In agreement with the hypothesis of Fels and Lindzen, the model tide accelerates the equatorial zonal wind near the cloud top level and decelerates it at higher levels. The tidal vertical wavelength, maximum amplitude, dissipative decay length, and zonal flow accelerations are sensitive functions of background angular velocity, in agreement with elementary gravity wave theory. In the equatorial cloud top region, tidal acceleration is comparable in magnitude to the decelerative effects of vertical advection and the model's Rayleigh friction damping. For sufficiently rapid initial zonal flow near the cloud top level, the area-weighted global mean cloud top level zonal wind increases with time over a 50-day model run as a result of tidal acceleration. Agreement between the model tide and the observed tide, or the tide determined in the more detailed calculations of Pechmann and Ingersoll, is best when the background angular velocity at the jet level is about 30% larger than that observed.     

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.     

Mapping zonal winds at Venus’s cloud tops from ground-based Doppler velocimetry

The most significant aspect of the general circulation of the atmosphere of Venus is its retrograde super-rotation. A complete characterization of this dynamical phenomenon is crucial for understanding its driving mechanisms. Here we report on ground-based Doppler velocimetry measurements of the zonal winds, based on high resolution spectra from the UV–Visual Echelle Spectrograph (UVES) instrument at ESO’s Very Large Telescope. Under the assumption of predominantly zonal flow, this method allows the simultaneous direct measurement of the zonal velocity across a range of latitudes and local times in the day side. The technique, based on long slit spectroscopy combined with the high spatial resolution provided by the VLT, has provided the first ground-based characterization of the latitudinal profile of zonal wind in the atmosphere of Venus, the first zonal wind field map in the visible, as well as new constraints on wind variations with local time. We measured mean zonal wind amplitudes between 106 ± 21 and 127 ± 14 m/s at latitudes between 18°N and 34°S, with the zonal wind being approximately uniform in 2.6°-wide latitude bands (0.3 arcsec at disk center). The zonal wind profile retrieved is consistent with previous spacecraft measurements based on cloud tracking, but with non-negligible variability in local time (longitude) and in latitude. Near 50° the presence of moderate jets is apparent in both hemispheres, with the southern jet being stronger by ～10 m/s. Small scale wind variations with local time are also present at low and mid-latitudes.     

Hadley circulations and Kelvin wave-driven equatorial jets in the atmospheres of Jupiter and Saturn

We propose a dynamical mechanism that can plausibly explain the origin of the broad prograde equatorial winds observed on Jupiter and Saturn, and examine the feasibility of this mechanism using two- (2D) and three-dimensional (3D) numerical simulation models. The idea is based on combining a narrow Gaussian jet peaking at the equator, which is induced by the momentum transfer from an upward propagating equatorial Kelvin-wave, and a pair of off-equatorial jets due to a meridional-vertical circulation similar to the tropical Hadley circulation on Earth. We employ for this feasibility study a 2D mechanistic mean-flow model which incorporates the influence of prescribed waves, and a 3D general circulation model, based on the generalised primitive equations of atmospheric motion. We then confirm that the dynamical models of both kinds can successfully reproduce theoretically expected flows of a reasonable magnitude, and that when two mechanisms are combined, a broad super-rotating jet is produced with off-equatorial maxima in zonal velocity for both Jupiter and Saturn, approximately in accordance with observations.     

Characteristics of EUV Coronal Jets Observed with STEREO/SECCHI

In this paper we present the first comprehensive statistical study of EUV coronal jets observed with the SECCHI (Sun Earth Connection Coronal and Heliospheric Investigation) imaging suites of the two STEREO spacecraft. A catalogue of 79 polar jets is presented, identified from simultaneous EUV and white-light coronagraph observations, taken during the time period March 2007 to April 2008, when solar activity was at a minimum. The twin spacecraft angular separation increased during this time interval from 2 to 48 degrees. The appearances of the coronal jets were always correlated with underlying small-scale chromospheric bright points. A basic characterization of the morphology and identification of the presence of helical structure were established with respect to recently proposed models for their origin and temporal evolution. Though each jet appeared morphologically similar in the coronagraph field of view, in the sense of a narrow collimated outward flow of matter, at the source region in the low corona the jet showed different characteristics, which may correspond to different magnetic structures. A classification of the events with respect to previous jet studies shows that amongst the 79 events there were 37 Eiffel tower-type jet events, commonly interpreted as a small-scale (～35 arc?sec) magnetic bipole reconnecting with the ambient unipolar open coronal magnetic fields at its loop tops, and 12 lambda-type jet events commonly interpreted as reconnection with the ambient field happening at the bipole footpoints. Five events were termed micro-CME-type jet events because they resembled the classical coronal mass ejections (CMEs) but on much smaller scales. The remaining 25 cases could not be uniquely classified. Thirty-one of the total number of events exhibited a helical magnetic field structure, indicative for a torsional motion of the jet around its axis of propagation. A few jets are also found in equatorial coronal holes. In this study we present sample events for each of the jet types using both, STEREO A and STEREO B, perspectives. The typical lifetimes in the SECCHI/EUVI (Extreme UltraViolet Imager) field of view between 1.0 to 1.7 R _⊙ and in SECCHI/COR1 field of view between 1.4 to 4 R _⊙ are obtained, and the derived speeds are roughly estimated. In summary, the observations support the assumption of continuous small-scale reconnection as an intrinsic feature of the solar corona, with its role for the heating of the corona, particle acceleration, structuring and acceleration of the solar wind remaining to be explored in more detail in further studies.     

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.     

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.     

Decadal variations in a Venus general circulation model

The Community Atmosphere Model (CAM), a 3-dimensional Earth-based climate model, has been modified to simulate the dynamics of the Venus atmosphere. The most current finite volume version of CAM is used with Earth-related processes removed, parameters appropriate for Venus introduced, and some basic physics approximations adopted. A simplified Newtonian cooling approximation has been used for the radiation scheme. We use a high resolution (1° by 1° in latitude and longitude) to take account of small-scale dynamical processes that might be important on Venus. A Rayleigh friction approach is used at the lower boundary to represent surface drag, and a similar approach is implemented in the uppermost few model levels providing a ‘sponge layer’ to prevent wave reflection from the upper boundary. The simulations generate superrotation with wind velocities comparable to those measured in the Venus atmosphere by probes and around 50-60% of those measured by cloud tracking. At cloud heights and above the atmosphere is always superrotating with mid-latitude zonal jets that wax and wane on an approximate 10 year cycle. However, below the clouds, the zonal winds vary periodically on a decadal timescale between superrotation and subrotation. Both subrotating and superrotating mid-latitude jets are found in the approximate 40-60 km altitude range. The growth and decay of the sub-cloud level jets also occur on the decadal timescale. Though subrotating zonal winds are found below the clouds, the total angular momentum of the atmosphere is always in the sense of superrotation. The global relative angular momentum of the atmosphere oscillates with an amplitude of about 5% on the approximate 10 year timescale. Symmetric instability in the near surface equatorial atmosphere might be the source of the decadal oscillation in the atmospheric state. Analyses of angular momentum transport show that all the jets are built up by poleward transport by a meridional circulation while angular momentum is redistributed to lower latitudes primarily by transient eddies. Possible changes in the structure of Venus’ cloud level mid-latitude jets measured by Mariner 10, Pioneer Venus, and Venus Express suggest that a cyclic variation similar to that found in the model might occur in the real Venus atmosphere, although no subrotating winds below the cloud level have been observed to date. Venus’ atmosphere must be observed over multi-year timescales and below the clouds if we are to understand its dynamics.     

A tidal explanation for the Titan haze layers

Though Titan is in synchronous rotation around Saturn, it experiences gravitational tides as a consequence of its eccentric orbit. It is proposed that the vertical transport of aerosols by these tides produces the haze layers in Titan's upper atmosphere. Analysis shows that the zonal winds in Titan's superrotating atmosphere have a profound influence on which tidal components are effective in establishing the multiple detached-haze layers. If the Huygens Doppler winds are representative of the equatorial global superrotation, then the westward propagating s=2 mode is the responsible tidal component even though its forcing is significantly weaker than that of the s=0 and eastward s=2 components. The eastward s=2 tidal mode is eliminated by critical levels while the s=0 mode is viscously damped in the strong high altitude winds. At polar latitudes, however, the gravest s=0 mode is the one most likely to produce layering. It is also suggested that the atmospheric gravitational tides could be responsible for decelerating the superrotating atmosphere as seen in the Huygens Doppler wind velocity profile at about 80 km altitude.     

Saturn's cloud morphology and zonal winds before the Cassini encounter

Analyses of Hubble Space Telescope (HST) images of Saturn obtained from August 2003 to March 2004, with extensive support from ground-based telescopes, have been used to characterize the cloud morphology and motions in its atmosphere few months before the Cassini encounter. We present data on the major meteorological features as potential targets for Cassini observations. We extend our previous measurements of the zonal winds during the 1996-2002 period (A. Sánchez-Lavega et al., 2003, Nature, 423, 623-625), confirming the strong change in the equatorial jet, and the high hemispheric symmetry of the zonal wind pattern.