A seasonal model of the Saturnian upper troposphere: Comparison with Voyager infrared measurements

An extension of the seasonal climate model of R. D. Cess and J. Caldwell (1979, Icarus, 38, 349–357) to Saturn's upper troposphere is presented. The ring-modulated latitudinal dependence of the insolation, the ring thermal emission, the oblateness of the planet, the orbit eccentricity, and the latitudinal variation of the internal heat flux are taken into account. Calculations agree closely with the temperature—latitude profiles retrieved from Voyager IRIS measurements at atmospheric levels located above the 0.2-bar pressure level; they reproduce the observed large-scale hemispheric asymmetry which is then shown to result from the seasonally variable insolation. Aerosol absorption is found to be the dominant source of atmospheric solar heating in the troposphere and the model suggests an aerosol mean unit optical depth around the 0.25-bar level in the equatorial region and around the 0.35-bar level at other latitudes. The model fails to predict the retrieved temperature—latitude profiles below the 0.3-bar level. This discrepancy is attributed to the existence of clouds at these levels which are responsible for an additional far-infrared opacity not taken into account in the temperature retrieval. The cloud-top altitude would be about 0.3 bar except in the 20 to 40°N region where these clouds would be confined below the 0.6-bar level. The poor correlation between infrared measurements and visible images is discussed and a possible model of Saturn's cloud structure is proposed.     

Zonal winds at high latitudes on Venus: An improved application of cyclostrophic balance to Venus Express observations

Recent retrievals of zonal thermal winds obtained in a cyclostrophic regime on Venus are generally consistent with cloud tracking measurements at mid-latitudes, but become unphysical in polar regions where the values obtained above the clouds are often less than or close to zero. Using a global atmospheric model, we show that the main source of errors that appear in the polar regions when retrieving the zonal thermal winds is most likely due to uncertainties in the zonal wind intensity in the choice of the lower boundary condition.Here we suggest a new and robust method to better estimate the lower boundary condition for high latitudes, thereby improving the retrieved zonal thermal winds throughout the high latitudes middle atmosphere. This new method is applied to temperature fields derived from Visible and Infrared Thermal Imaging Spectrometer (VIRTIS) data on board the Venus Express spacecraft. We obtain a zonal thermal wind field that is in better agreement with other, more direct methods based on either retrieving the zonal winds from cloud tracking or from direct measurements of the meridional slope of pressure surfaces.     

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.     

Wind variations in Jupiter's equatorial atmosphere: A QQO counterpart?

Jupiter's equatorial atmosphere, much like the Earth's, is known to show quasi-periodic variations in temperature, particularly in the stratosphere, but variations in other jovian atmospheric tracers have not been studied for any correlations to these oscillations. Data taken at NASA's Infrared Telescope Facility (IRTF) from 1979 to 2000 were used to obtain temperatures at two levels in the atmosphere, corresponding to the upper troposphere (250 mbar) and to the stratosphere (20 mbar). We find that the data show periodic signals at latitudes corresponding to the troposphere zonal wind jets, with periods ranging from 4.4 (stratosphere, 95% confidence at 4° S planetographic latitude) to 7.7 years (troposphere, 97% confidence at 6° N). We also discuss evidence that at some latitudes the troposphere temperature variations are out of phase from the stratosphere variations, even where no periodicity is evident. Hubble Space Telescope images were used, in conjunction with Voyager and Cassini data, to track small changes in the troposphere zonal winds from 20° N to 20° S latitude over the 1994-2000 time period. Oscillations with a period of 4.5 years are found near 7°-8° S, with 80-85% significance. Further, the strongest evidence for a QQO-induced tropospheric wind change tied to stratospheric temperature change occurs near these latitudes, though tropospheric temperatures show little periodicity here. Comparison of thermal winds and measured zonal winds for three dates indicate that cloud features at other latitudes are likely tracked at pressures that can vary by up to a few hundred millibar, but the cloud altitude change required is too large to explain the wind changes measured at 7° S.     

Dynamical properties of the Venus mesosphere from the radio-occultation experiment VeRa onboard Venus Express

The dynamics of Venus’ mesosphere (60–100 km altitude) was investigated using data acquired by the radio-occultation experiment VeRa on board Venus Express. VeRa provides vertical profiles of density, temperature and pressure between 40 and 90 km of altitude with a vertical resolution of few hundred meters of both the Northern and Southern hemisphere. Pressure and temperature vertical profiles were used to derive zonal winds by applying an approximation of the Navier–Stokes equation, the cyclostrophic balance, which applies well on slowly rotating planets with fast zonal winds, like Venus and Titan. The main features of the retrieved winds are a midlatitude jet with a maximum speed up to 140 ± 15 m s^?1 which extends between 20°S and 50°S latitude at 70 km altitude and a decrease of wind speed with increasing height above the jet. Cyclostrophic winds show satisfactory agreement with the cloud-tracked winds derived from the Venus Monitoring Camera (VMC/VEx) UV images, although a disagreement is observed at the equator and near the pole due to the breakdown of the cyclostrophic approximation. Knowledge of both temperature and wind fields allowed us to study the stability of the atmosphere with respect to convection and turbulence. The Richardson number Ri was evaluated from zonal field of measured temperatures and thermal winds. The atmosphere is characterised by a low value of Richardson number from ～45 km up to ～60 km altitude at all latitudes that corresponds to the lower and middle cloud layer indicating an almost adiabatic atmosphere. A high value of Richardson number was found in the region of the midlatitude jet indicating a highly stable atmosphere. The necessary condition for barotropic instability was verified: it is satisfied on the poleward side of the midlatitude jet, indicating the possible presence of wave instability.     

Thermospheric winds in the cusp: Dependence of the latitude of the cusp

The dayside thermospheric wind pattern as observed from Spitsbergen generally shows moderately strong westward winds with a small poleward component. The flow is almost zonal, frequently with sufficient westward velocity that parcels of air cross the noon meridian travelling towards the morning before turning antisunward towards the nightside across the polar cap. There have been some exceptions which are characterized by much weaker winds having been increased in the poleward direction but with a very much reduced westward component. Making use of the meridian scanning photometer data obtained simultaneously on the same site, it is clearly shown that the normal behaviour occurs when the cusp, as indicated by the region of high $\text{630}\text{428}\text{nm}$ and $\text{630}\text{558}\text{nm}$ photometric intensity ratios, is to the North of the station. Just below the latitude of the cusp, the strong thermospheric flow generated by neutral coupling to the strong westward convection in the dusk sector continues across the dayside. It is maintained in the zonal direction because of the balance between the poleward Coriolis force and the equatorward pressure force caused by cusp heating. Poleward of the high pressure region at the cusp the flow is diverted northward and initially makes much slower progress across the Polar Cap. When the auroral oval has expanded such that the cusp is well to the South of our Spitsbergen station, the thermosphere in the sampled region has been found to be within this slow flow zone. On such occasions, the nightside speeds are well in excess of those on the dayside, in contrast to the normal behavior of comparable dayside and nightside wind speeds.     

4-Day waves in the Venus atmosphere

Ultraviolet albedo contrasts in the Venus atmosphere are probably large-scale atmospheric waves propagating slowly with respect to the rapid cloud-top zonal winds. Using a simple theoretical model and profiles of mean wind and thermal structure based on Pioneer Venus data, we find planetary-scale gravity waves with phase velocities matching the speeds of the uv markings. We propose an upward-propagating wave and waves trapped at cloud levels as candidates to explain the observed uv features.     

Venus mesosphere and thermosphere temperature structure: II. Day-night variations

A two-dimensional nonlinear hydrodynamic model has been developed for studying the global scale winds, temperature, and compositional structure of the mesosphere and thermosphere of Venus. The model is driven by absorption of solar radiation. Ultraviolet radiation produces both heating and photodissociation. Infrared solar heating and thermal cooling are also included with an accurate NLTE treatment. The most crucial uncertainty in determining the solar drive is the efficiency by which $\text{λ < 1080}\text{A}\text{?}$ solar radiation is converted to heat. This question was analyzed in Part I, where it was concluded that essentially all hot atom and O(¹D) energy may be transferred to vibrational-rotational energy of CO₂ molecules. If this is so, the minimum possible euv heating occurs and is determined by the quenching of the resulting excess rotational energy. The hydrodynamic model is integrated with this minimum heating and neglecting any small-scale vertical eddy mixing. The results are compared with predictions of another model with the same physics except that it assumes that 30% of $\text{λ < 1080}\text{A}\text{?}$ radiation goes into heat and that the heating from longer-wavelength radiation includes the O(¹D) energy. For the low-efficiency model, exospheric temperatures are ?300°K on the dayside and drop to < 180°K at the antisolar point. For the higher-efficiency model, the day-to-night temperature variation is from ?600°K to ?250°K. Both versions of the model predict a wind of several hundred meters per second blowing across the terminator and abruptly weakening to small values on the nightside with the mass flow consequently going into a strong tongue of downward motion on the nightside of the terminator. The presence of this circulation could be tested observationally by seeing if its signature can be found in temperature measurements. Both versions of the model indicate that a self-consistent large-scale circulation would maintain oxygen concentrations with ?5% mixing ratios near the dayside F-1 ionospheric peak but ?40% at the antisolar point at the same pressure level.     

Dispersion in Neptune’s zonal wind velocities from NIR Keck AO observations in July 2009

We report observations of Neptune made in H-(1.4–1.8 μm) and K’-(2.0–2.4 μm) bands on 14 and 16 July 2009 from the 10-m W.M. Keck II Telescope using the near-infrared camera NIRC2 coupled to the Adaptive Optics (AO) system. We track the positions of 54 bright atmospheric features over a few hours to derive their zonal and latitudinal velocities, and perform radiative transfer modeling to measure the cloud-top pressures of 50 features seen simultaneously in both bands. We observe one South Polar Feature (SPF) on 14 July and three SPFs on 16 July at ～65?°S. The SPFs observed on both nights are different features, consistent with the high variability of Neptune’s storms. There is significant dispersion in Neptune’s zonal wind velocities about the smooth Voyager wind profile fit of Sromovsky et al. (Icarus, 105:140, 1993), much greater than the upper limit we expect from vertical wind shear, with the largest dispersion seen at equatorial and southern mid-latitudes. Comparison of feature pressures vs. residuals in zonal velocity from the smooth Voyager wind profile also directly reveals the dominance of mechanisms over vertical wind shear in causing dispersion in the zonal winds. Vertical wind shear is not the primary cause of the difference in dispersion and deviation in zonal velocities between features tracked in H-band on 14 July and those tracked in K’-band on 16 July. Dispersion in the zonal velocities of features tracked over these short time periods is dominated by one or more mechanisms, other than vertical wind shear, that can cause changes in the dispersion and deviation in the zonal velocities on timescales of hours to days.     

Wind measurements in Mars' middle atmosphere: IRAM Plateau de Bure interferometric CO observations

The IRAM Plateau de Bure Interferometer has been used to map the CO(1-0) rotational line in Mars' middle atmosphere. Absolute winds and thermal profiles were retrieved during the 1999, 2001, 2003 and 2005 planet's oppositions. The observations sampled various seasons (Ls=143, 196, 262, 317 and 322), and different dust situations (clear, global storm, regional storm). The absolute winds were derived by measuring directly the Doppler lineshifts. The main zonal circulation near 50 km is dominated by strong retrograde winds, with typical velocities of 70-170 m/s, strongly varying seasonally, latitudinally, and longitudinally (in particular between morning and evening). Comparison of the retrieved temperature with a general circulation model indicates that the model often underestimates the temperatures in the middle (20-50 km) atmosphere, and overestimates them above 50 km.     

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.     

Uranus and Neptune: Shape and rotation

Both Uranus and Neptune are thought to have strong zonal winds with velocities of several 100 m s⁻¹. These wind velocities, however, assume solid-body rotation periods based on Voyager 2 measurements of periodic variations in the planets’ radio signals and of fits to the planets’ magnetic fields; 17.24 h and 16.11 h for Uranus and Neptune, respectively. The realization that the radio period of Saturn does not represent the planet’s deep interior rotation and the complexity of the magnetic fields of Uranus and Neptune raise the possibility that the Voyager 2 radio and magnetic periods might not represent the deep interior rotation periods of the ice giants. Moreover, if there is deep differential rotation within Uranus and Neptune no single solid-body rotation period could characterize the bulk rotation of the planets. We use wind and shape data to investigate the rotation of Uranus and Neptune. The shapes (flattening) of the ice giants are not measured, but only inferred from atmospheric wind speeds and radio occultation measurements at a single latitude. The inferred oblateness values of Uranus and Neptune do not correspond to bodies rotating with the Voyager rotation periods. Minimization of wind velocities or dynamic heights of the 1 bar isosurfaces, constrained by the single occultation radii and gravitational coefficients of the planets, leads to solid-body rotation periods of ∼16.58 h for Uranus and ∼17.46 h for Neptune. Uranus might be rotating faster and Neptune slower than Voyager rotation speeds. We derive shapes for the planets based on these rotation rates. Wind velocities with respect to these rotation periods are essentially identical on Uranus and Neptune and wind speeds are slower than previously thought. Alternatively, if we interpret wind measurements in terms of differential rotation on cylinders there are essentially no residual atmospheric winds.     

Titan: Far-infrared and microwave remote sensing of methane clouds and organic haze

It is shown that Titan's surface and plausible atmospheric thermal opacity sources—gaseous N₂, CH₄, and H₂, CH₄ cloud, and organic haze—are sufficient to match available Earth-based and Voyager observations of Titan's thermal emission spectrum. Dominant sources of thermal emission are the surface for wavelenghts λ ? 1 cm, atmospheric N₂ for 1 cm ? λ ? 200 μm,, condensed and gaseous CH₄ for 200 μm ? λ ? 20 μm, and molecular bands and organic haze for λ ? 20 μm. Matching computed spectra to the observed Voyager IRIS spectra at 7.3 and 52.7° emission angles yields the following abundances and locations of opacity sources: CH₄ clouds: 0.1 g cm^? at a planetocentric radius of 2610–2625 km, 0.3 g cm^?2 at 2590–2610 km, total 0.4 ± 0.1 g cm^–2 above 2590 km; organic haze: $\text{4 ± 2 × 10}{}^{\mathrm{?6}}\text{,}\text{g cm}\text{,}{}^{\mathrm{?2}}$ above 2750 km; tropospheric H₂: 0.3 ± 0.1 mol%. This is the first quantitative estimate of the column density of condensed methane (or CH₄/C₂H₆) on Titan. Maximum transparency in the middle to far IR occurs at 19 μm where the atmospheric vertical absorption optical depth is $\text{?0.6}$ A particle radius $\text{r}\text{? 2 μ}\text{m}$ in the upper portion of the CH₄ cloud is indicated by the apparent absence of scattering effects.     

Martian winds and dust clouds

Previous calculations of the surface wind stress required to raise dust on Mars are reconsidered and the threshold friction velocity is found to be about 2.0 m sec^?1 with particles of 200–300 μm being the most easily lifted. With this friction velocity, the planetary resistance law yields a corresponding wind at the top of the Ekman layer of 60 m sec^?1, and the logarithmic wind law yields a corresponding wind at the top of the Prandtl layer of 38 m sec^?1. These speeds are somewhat lower than those used by previous investigators.Various mechanisms for producing such strong winds are examined and it is concluded that the general circulation, thermal effects of topography, mechanical effects of topography and dust devils are all capable of doing so.Dust storms associated with small-scale disturbances are found to be incapable of growth. A scaling analysis of the equations of horizontal motion and of hydrostatic balance shows that a dust cloud at least 10 km thick and several tens of km in radius can, by absorption of sunlight, generate temperature gradients that, in turn, produce winds capable of raising more dust. Thus, a feedback mechanism is suggested in which an initial dust cloud exceeding certain critical dimensions can grown to planetary size. The preference of large dust storms to occur at southern hemisphere summer solstice is attributed to the maximum of insolation at that time. It is suggested that the frequent origin in the Noachis-Hellas region may be due to orographic features of the right scale and to low height in that area.     

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.     

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.     

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.     

On the maintenance of thermal wind balance and equatorial superrotation in Titan's stratosphere

Titan's atmospheric winds, like those on Venus, exhibit superrotation at high altitudes. Titan general circulation models have yielded conflicting results on whether prograde winds in excess of 100 m s⁻¹ at the 1 mbar level are possible based on known physical processes that drive wind systems. A comprehensive two-dimensional (2D) model for Titan's stratosphere was constructed to systematically explore the physical mechanisms that produce and maintain stratospheric wind systems. To ensure conservation of angular momentum in the limit of no net exchange of atmospheric angular momentum with the solid satellite and no external sources and sinks, the zonal momentum equation was solved in flux form for total angular momentum. The relationships among thermal wind balance, meridional circulation, and zonal wind were examined with numerical experiments over a range of values for fundamental input parameters, including planetary rotation rate, radius, internal friction due to wave stresses, and net radiative drive. The magnitude of mid-latitude jets is most sensitive to a single parameter, the planetary rotation rate and results from the conversion of planetary angular momentum to relative angular momentum by the meridional circulation, whereas the strength of meridional circulation is mainly determined by the magnitude of the radiative drive. For Titan's slowly rotating atmosphere, the meridional temperature gradient is vanishingly small, even when the radiative drive is enhanced beyond reasonable magnitudes, and can be inferred from zonal winds in gradient/thermal wind balance. In our 2D model large equatorial superrotation in Titan's stratosphere can be only produced through internal drag forcing by eddy momentum fluxes, which redistribute angular momentum within the atmosphere, while still conserving the total angular momentum of the atmosphere with time. We cannot identify any waves, such as gravitational or thermal tides, that are sufficiently capable of generating the required eddy forcing of >50 m s⁻¹ Titan-day⁻¹ to maintain peak prograde winds in excess of 100 m s⁻¹ at the 1 mbar level.