Determining a tilt in Titan’s north-south albedo asymmetry from Cassini images

Analysis of Titan’s hemispheric brightness asymmetry from mapped Cassini images reveals an axis of symmetry that is tilted with respect to the rotational axis of the solid body. Twenty images taken from 2004 through 2007 show a mean axial offset of 3.8 ± 0.9° relative to the solid body’s pole, directed 79 ± 24° to the west of the sub-solar longitude. These values are consistent with recent measurements of an implied atmospheric spin axis determined from isothermal mapping by [Achterberg, R.K., Conrath, B.J., Gierasch, P.J., Flasar, F.M., Nixon, C.A., 2008. Icarus 197, 549-555].     

Compositional evidence for Titan's stratospheric tilt

Five years of Cassini CIRS infrared spectra have been used to determine the tilt of Titan's stratospheric symmetry axis with respect to the solid body rotation axis. Measurements of HCN abundance centred around 5 mbar (125 km altitude) at equatorial latitudes show the symmetry axis is tilted by 4.0±1.5^° in a direction 70±40^°W of the sub-solar point. This value is consistent with tilts determined from temperature and haze measurements by Achterberg et al. (2008a) and Roman et al. (2009). The consistency of results from three independent methods suggests that Titan's entire stratosphere is tilted and provides a powerful constraint on the underlying atmospheric dynamics.     

Latitudinal transport by barotropic waves in Titan's stratosphere.: I. General properties from a horizontal shallow-water model

We present a numerical study of barotropic waves in Titan's stratosphere based on a shallow-water model. The forcing of the zonal flow by the mean meridional circulation is represented by a relaxation towards a barotropically unstable wind profile. The relaxation profile is consistent with observations and with previous results from a 3D general circulation model. The time constant of the forcing that best matches the northward eddy-transport of zonal momentum from the 3D model is τ∼5 Titan days. The eddy wind field is a zonal wavenumber-2 wave with a peak amplitude about 10% of the mean wind speed. The latitudinal transport of angular momentum by the wave tends to keep the flow close to marginal stability by carrying momentum upgradient, from the core of the jets into the low latitudes. Although the strongest eddy motions occur at the latitudes of the wind maxima, the strongest mixing takes place at the barotropically unstable regions, close to ±30° and spanning about 30° in latitude. An eddy-mixing time constant of the order of 1 Titan day is inferred within these regions, and of a few tens of days within regions of stable flow. Horizontal gradients in transient tracer fields are less than 10% of the latitudinal gradient of the meridional tracer profile. Cassini's detection of such waves could provide a direct observation of wind speeds at stratospheric levels.     

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

Four climate regimes on a land planet with wet surface: Effects of obliquity change and implications for ancient Mars

Series of numerical experiments are performed using a general circulation model to gain insights on the hydrologic cycle on ancient Mars. Since the state of the ancient Mars atmosphere is not well constrained, we did not try to simulate an ancient Mars climate under warm and wet condition. In stead, we used an idealized model and tried to extract general features of the hydrologic cycle by modeling an ideal land planet that has no ocean on its surface. Four different climate regimes, “warm-upright,” “warm-oblique,” “frozen-upright,” and “frozen-oblique” regimes, are recognized depending on the inclination of the spin axis (obliquity) and average surface temperature. The period of active hydrologic cycle suggested from the geomorphology on Mars seems to be consistent with that at the “warm-oblique” regime, which appears at warm (above-freezing) environment with high-obliquity (higher than about 30°) condition.     

On the Origin of the Spin of Planets and Stars and its Connection with Gravitomagnetism

The origin of the spin of planets and stars is, to a certain extent, still unexplained. In general, we attribute their rotation to the swirl of their constituent primitive gases. In this paper, we try to show that the rotation of celestial bodies depends only on their mass, apparent radius and tilt of their spin axes. We reach this conclusion within the framework of gravitomagnetism, implied by the Einstein’s general relativity theory (GR). Our results show that it might possible, in principle, to calculate the mass of spinning objects by measuring their apparent radius, the speed of rotation and the tilt of the axis of rotation.     

A comprehensive numerical simulation of Io’s sublimation-driven atmosphere

Io’s sublimation-driven atmosphere is modeled using the direct simulation Monte Carlo (DSMC) method. These rarefied gas dynamics simulations improve upon earlier models by using a three-dimensional domain encompassing the entire planet computed in parallel. The effects of plasma heating, planetary rotation, inhomogeneous surface frost, molecular residence time of SO₂ on the exposed (non-volatile) rocky surface, and surface temperature distribution are investigated. Circumplanetary flow is predicted to develop from the warm dayside toward the cooler nightside. Io’s rotation leads to a highly asymmetric frost surface temperature distribution (due to the frost’s high thermal inertia) which results in circumplanetary flow that is not axi-symmetric about the subsolar point. The non-equilibrium thermal structure of the atmosphere, specifically vibrational and rotational temperatures, is also examined. Plasma heating is found to significantly inflate the atmosphere on both the dayside and nightside. The plasma energy flux causes high temperatures at high altitudes but plasma energy depletion through the dense gas column above the warmest frost permits gas temperatures cooler than the surface at low altitudes. A frost map (Douté, S., Schmitt, B., Lopes-Gautier, R., Carlson, R., Soderblom, L., Shirley, J., and the Galileo NIMS Team [2001]. Icarus 149, 107-132) is used to control the sublimated flux of SO₂ which can result in inhomogeneous column densities that vary by nearly a factor of four for the same surface temperature. A short residence time for SO₂ molecules on the “rock” component is found to smooth lateral atmospheric inhomogeneities caused by variations in the surface frost distribution, creating an atmosphere that looks nearly identical to one with uniform frost coverage. A longer residence time is found to agree better with mid-infrared observations (Spencer, J.R., Lellouch, E., Richter, M.J., López-Valverde, M.A., Jessup, K.L, Greathouse, T.K., Flaud, J. [2005]. Icarus 176, 283-304) and reproduce the observed anti-jovian/sub-jovian column density asymmetry. The computed peak dayside column density for Io assuming a surface frost temperature of 115 K agrees with those suggested by Lyman-α observations (Feaga, L.M., McGrath, M., Feldman, P.D. [2009]. Icarus 201, 570-584). On the other hand, the peak dayside column density at 120 K is a factor of five larger and is higher than the upper range of observations (Jessup, K.L., Spencer, J.R., Ballester, G.E., Howell, R.R., Roesler, F., Vigel, M., Yelle, R. [2004]. Icarus 169, 197-215; Spencer et al., 2005).     

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.     

Power spectral analysis of Jupiter’s clouds and kinetic energy from Cassini

We present suggestive evidence for an inverse energy cascade within Jupiter’s atmosphere through a calculation of the power spectrum of its kinetic energy and its cloud patterns. Using Cassini observations, we composed full-longitudinal mosaics of Jupiter’s atmosphere at several wavelengths. We also utilized image pairs derived from these observations to generate full-longitudinal maps of wind vectors and atmospheric kinetic energy within Jupiter’s troposphere. We computed power spectra of the image mosaics and kinetic energy maps using spherical harmonic analysis. Power spectra of Jupiter’s cloud patterns imaged at certain wavelengths resemble theoretical spectra of two-dimensional turbulence, with power-law slopes near −5/3 and −3 at low and high wavenumbers, respectively. The slopes of the kinetic energy power spectrum are also near −5/3 at low wavenumbers. At high wavenumbers, however, the spectral slopes are relatively flatter than the theoretical prediction of −3. In addition, the image mosaic and kinetic energy power spectra differ with respect to the location of the transition in slopes. The transition in slope is near planetary wavenumber 70 for the kinetic energy spectra, but is typically above 200 for the image mosaic spectra. Our results also show the importance of calculating spectral slopes from full 2D velocity maps rather than 1D zonal mean velocity profiles, since at large wavenumbers the spectra differ significantly, though at low wavenumbers, the 1D zonal and full 2D kinetic energy spectra are practically indistinguishable. Furthermore, the difference between the image and kinetic energy spectra suggests some caution in the interpretation of power spectrum results solely from image mosaics and its significance for the underlying dynamics. Finally, we also report prominent variations in kinetic energy within the equatorial jet stream that appear to be associated with the 5 μm hotspots. Other eddies are present within the flow collar of the Great Red Spot, suggesting caution when interpreting snapshots of the flow inside these features as representative of a time-averaged state.     

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.     

Cosmological considerations regarding Uranus

The cosmogony of Uranus is discussed within the context of a picture in which solid condensed materials accumulate to form a large body, which then acquires significant amounts of gas from the primitive solar nebula. Of prime cosmogonical importance is the tilt of the equatorial plane of the planet and of the plane of tilt of the planet can easily occur as a result of a major collision during the formation process; it seems most likely that the tilt of the satellite orbits requires that they were formed from a gaseous disc rotating about the planet after the tilt of the planetary rotational axis had occurred. Possible methods for tilting this gaseous disc are discussed. A strong early magnetic field may have helped in this and may have played an essential role in showing down the spin of the planet to the present observed value. These processes may have produced significant compositional differences between the satellites of Uranus and those of Jupiter and Saturn.     

Quasi-periodic dust events in the summertime south polar region of Mars

A Hovmöller diagram analysis of the dust optical depth measured by the Mars Global Surveyor Thermal Emission Spectrometer shows the occurrence of quasi-periodic westwardly-propagating disturbances with timescales of 10-20 sols during summer in the south polar region of Mars. Dust clouds emerge repeatedly around the region with a latitude of around 70-80°S and a longitude of 240-300°E, move westward at speeds of 3-6 m s⁻¹, reach the region with a longitude of 60-120°E, and finally disappear. This longitude range coincides with elevated terrains in the south polar region, and in this region an increase of dust optical depth encircling the south pole is also observed. This implies that the quasi-periodic dust events will contribute to the enhancement of the atmospheric dust loading in this region. These dust events might be related to baroclinic instability caused by the thermal contrast across the CO₂ cap edge, or the horizontal advection or vertical convection with radiative-dynamical feedback. The westward movement of the dust clouds suggests steady westward winds blowing in the near-surface layer, where the quasi-periodic dust lifting is expected to occur. Such a westward cap-edge flow will be created by the Coriolis force acting on the flow from the ice side to the regolith side.     

Forced waves in the martian atmosphere from MGS TES nadir data

We have analyzed the temperature retrievals from Mars Global Surveyor (MGS) Thermal Emission Spectrometer (TES) nadir spectra to yield latitude-height resolved maps of various atmospheric forced wave modes as a function of season for a full Mars year. Among the isolated wave modes is the zonal mean, time mean temperature, which we used to derive zonal mean zonal winds and stationary wave quasi-geostrophic indices of refraction, diagnostic of their propagation. The diurnal Kelvin wave was isolated in the data, with results roughly consistent with models (Wilson and Hamilton, 1996, J. Atmos. Sci. 33, 1290-1326). The s = 1 and s = 2 stationary waves were found to have significant amplitude in ducts extending up the winter polar jets, while the s = 3 stationary wave was found to be confined to near the surface. The s = 1 stationary wave was found to have little phase tilt with height during northern winter, but significant westward phase tilt with height in the southern winter. This indicates that the wave carries heat poleward, slightly more than that found in Barnes et al. (1996; J. Geophys. Res. 101, 12,753-12,776). The s = 1 stationary wave is likely the dominant mechanism for eddy meridional heat transport for the southern winter. We noted that the phase of the s = 2 stationary wave is nearly constant with time, but that the s = 1 stationary wave moved 90° of longitude from fall to winter and back in spring in the North. While interannual variability is not yet addressed, overall, these results provide the first comprehensive benchmark for forced waves in Mars’s atmosphere against which future atmospheric models of Mars can be compared.     

Secular spin dynamics of inner main-belt asteroids

Understanding the evolution of asteroid spin states is challenging work, in part because asteroids have a variety of orbits, shapes, spin states, and collisional histories but also because they are strongly influenced by gravitational and non-gravitational (YORP) torques. Using efficient numerical models designed to investigate asteroid orbit and spin dynamics, we study here how several individual asteroids have had their spin states modified over time in response to these torques (i.e., 951 Gaspra, 60 Echo, 32 Pomona, 230 Athamantis, 105 Artemis). These test cases which sample semimajor axis and inclination space in the inner main belt, were chosen as probes into the large parameter space described above. The ultimate goal is to use these data to statistically characterize how all asteroids in the main belt population have reached their present-day spin states. We found that the spin dynamics of prograde-rotating asteroids in the inner main belt is generally less regular than that of the retrograde-rotating ones because of numerous overlapping secular spin-orbit resonances. These resonances strongly affect the spin histories of all bodies, while those of small asteroids (?40 km) are additionally influenced by YORP torques. In most cases, gravitational and non-gravitational torques cause asteroid spin axis orientations to vary widely over short (?1 My) timescales. Our results show that (951) Gaspra has a highly chaotic rotation state induced by an overlap of the s and s₆ spin-orbit resonances. This hinders our ability to investigate its past evolution and infer whether thermal torques have acted on Gaspra's spin axis since its origin.     

The atmospheric signature of Io's Prometheus plume and anti-jovian hemisphere: evidence for a sublimation atmosphere

Using the Hubble Space Telescope's Space Telescope Imaging Spectrograph we have obtained for the first time spatially resolved 2000-3000 Å spectra of Io's Prometheus plume and adjoining regions on Io's anti-jovian hemisphere in the latitude range 60° N-60° S, using a 0.1″ slit centered on Prometheus and tilted roughly 45° to the spin axis. The SO₂ column density peaked at 1.25×10¹⁷ cm⁻² near the equator, with an additional 5×10¹⁶ cm⁻² enhancement over Prometheus corresponding to a model volcanic SO₂ output of 10⁵ kg s⁻¹. Apart from the Prometheus peak, the SO₂ column density dropped fairly smoothly away from the subsolar point, even over regions that included potential volcanic sources. At latitudes less than ±30°, the dropoff rate was consistent with control by vapor pressure equilibrium with surface frost with subsolar temperature 117.3±0.6 K, though SO₂ abundance was higher than predicted by vapor pressure control at mid-latitudes, especially in the northern hemisphere. We conclude that, at least at low latitudes on the anti-jovian hemisphere where there are extensive deposits of optically-thick SO₂ frost, the atmosphere is probably primarily supported by sublimation of surface frost. Although the 45° tilt of our slit prevents us from separating the dependence of atmospheric density on solar zenith angle from its dependence on latitude, the pattern is consistent with a sublimation atmosphere regardless of which parameter is the dominant control. The observed drop in gas abundance towards higher latitudes is consistent with the interpretation of previous Lyman alpha images of Io as indicating an atmosphere concentrated at low latitudes. Comparison with previous disk-resolved UV spectroscopy, Lyman-alpha images, and mid-infrared spectroscopy suggests that Io's atmosphere is denser and more widespread on the anti-jovian hemisphere than at other longitudes. SO₂ gas temperatures were in the range of 150-250 K over the majority of the anti-jovian hemisphere, consistent with previous observations. SO was not definitively detected in our spectra, with upper limits to the SO/SO₂ ratio in the range 1-10%, roughly consistent with previous observations. S₂ gas was not seen anywhere, with an upper limit of 7.5×10¹⁴ cm⁻² for the Prometheus plume, confirming that this plume is significantly poorer in S₂ than the Pele plume (S₂ /SO₂<0.005, compared to 0.08-0.3 at Pele). In addition to the gas absorption signatures, we have observed continuum emission in the near ultraviolet (near 2800 Å) for the first time. The brightness of the observed emission was directly correlated with the SO₂ abundance, strongly peaking in the equatorial region over Prometheus. Emission brightness was modestly anti-correlated with the jovian magnetic latitude, decreasing when Io intersected the torus centrifugal equator.     

Cassini CIRS retrievals of ammonia in Jupiter's upper troposphere

The zonal mean ammonia abundance on Jupiter between the 400- and 500-mbar pressure levels is inferred as a function of latitude from Cassini Composite Infrared Spectrometer data. Near the Great Red Spot, the ammonia abundance is mapped as a function of latitude and longitude. The Equatorial Zone is rich in ammonia, with a relative humidity near unity. The North and South Equatorial Belts are depleted relative to the Equatorial Zone by an order of magnitude. The Great Red Spot shows a local maximum in the ammonia abundance. Ammonia abundance is highly correlated with temperature perturbations at the same altitude. Under the assumption that anomalies in ammonia and temperature are both perturbed from equilibrium by vertical motion, we find that the adjustment time constant for ammonia equilibration is about one third of the radiative time constant.     

A fully 3-dimensional thermal model of a comet nucleus

A 3-D numerical model of comet nuclei is presented. An implicit numerical scheme was developed for the thermal evolution of a spherical nucleus composed of a mixture of ice and dust. The model was tested against analytical solutions, simplified numerical solutions, and 1-D thermal evolution codes. The 3-D code was applied to comet 67P/Churyumov-Gerasimenko; surface temperature maps and the internal thermal structure was obtained as function of depth, longitude and hour angle. The effect of the spin axis tilt on the surface temperature distribution was studied in detail. It was found that for small tilt angles, relatively low temperatures may prevail on near-pole areas, despite lateral heat conduction. A high-resolution run for a comet model of 67P/Churyumov-Gerasimenko with low tilt angle, allowing for crystallization of amorphous ice, showed that the amorphous/crystalline ice boundary varies significantly with depth as a function of cometary latitude.     

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.     

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.     

Latitudinal transport by barotropic waves in Titan's stratosphere.: II. Results from a coupled dynamics-microphysics-photochemistry GCM

We present a 2D general circulation model of Titan's atmosphere, coupling axisymmetric dynamics with haze microphysics, a simplified photochemistry and eddy mixing. We develop a parameterization of latitudinal eddy mixing by barotropic waves based on a shallow-water, longitude-latitude model. The parameterization acts locally and in real time both on passive tracers and momentum. The mixing coefficient varies exponentially with a measure of the barotropic instability of the mean zonal flow. The coupled GCM approximately reproduces the Voyager temperature measurements and the latitudinal contrasts in the distributions of HCN and C₂H₂, as well as the main features of the zonal wind retrieved from the 1989 stellar occultation. Wind velocities are consistent with the observed reversal time of the North-South albedo asymmetry of 5 terrestrial years. Model results support the hypothesis of a non-uniform distribution of infrared opacity as the cause of the Voyager temperature asymmetry. Transport by the mean meridional circulation, combined with polar vortex isolation may be at the origin of the latitudinal contrasts of trace species, with eddy mixing remaining restricted to low latitudes most of the Titan year. We interpret the contrasts as a signature of non-axisymmetric motions.