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.     

Venus Doppler winds at cloud tops observed with ESPaDOnS at CFHT

We present new wind measurements in Venus’ lower mesosphere from visible spectroscopy during the 2007 worldwide coordinated ground campaign in support of ESA's Venus Express mission. These observations consisted of high-resolution spectra of Fraunhofer lines in the entire visible range (0.37-1.05 μm) to measure the winds near 68 km using the Doppler shift of solar radiation scattered by clouds toward the observer's direction. The observations included various points of the dayside hemisphere at a phase angle of ∼109°. We took advantage of two symmetrical elongations in July and September 2007 at Canada-France-Hawaii's 3.6-m telescope. Kinematical fits to the Doppler winds provide a mean equatorial velocity of (104±10) m s⁻¹ for the zonal retrograde flow. This velocity agrees quite well with the mean value obtained by tracking the UV markings from several spacecraft.     

Dynamics of Jupiter’s equatorial region at cloud top level from Cassini and HST images

We present a study of the equatorial region of Jupiter, between latitudes ∼15°S and ∼15°N, based on Cassini ISS images obtained during the Jupiter flyby at the end of 2000, and HST images acquired in May and July 2008. We examine the structure of the zonal wind profile and report the detection of significant longitudinal variations in the intensity of the 6°N eastward jet, up to 60 m s⁻¹ in Cassini and HST observations. These longitudinal variations are, in the HST case, associated with different cloud morphology. Photometric and radiative transfer analysis of the cloud features used as tracers in HST images show that at most there is only a small height difference, no larger than ∼0.5-1 scale heights, between the slow (∼100 m s⁻¹) and fast (∼150 m s⁻¹) moving features. This suggests that speed variability at 6°N is not dominated by vertical wind shears but instead we propose that Rossby wave activity is the responsible for the zonal variability. Removing this variability, we find that Jupiter’s equatorial jet is actually symmetric relative to equator with two peaks of ∼140-150 m s⁻¹ located at latitudes 6°N and 6°S and at a similar pressure level. We also study the local dynamics of particular equatorial features such as several dark projections associated with 5 μm hot spots and a large, long-lived feature called the White Spot (WS) located at 6°S. Convergent flow at the dark projections appears to be a characteristic which depends on the particular morphology and has only been detected in some cases. The internal flow field in the White Spot indicates that it is a weakly rotating quasi-equatorial anticyclone relative to the ambient meridionally sheared flow.     

Stratospheric superrotation in the TitanWRF model

TitanWRF general circulation model simulations performed without sub-grid-scale horizontal diffusion of momentum produce roughly the observed amount of superrotation in Titan’s stratosphere. We compare these results to Cassini-Huygens measurements of Titan’s winds and temperatures, and predict temperature and winds at future seasons. We use angular momentum and transformed Eulerian mean diagnostics to show that equatorial superrotation is generated during episodic angular momentum ‘transfer events’ during model spin-up, and maintained by similar (yet shorter) events once the model has reached steady state. We then use wave and barotropic instability analysis to suggest that these transfer events are produced by barotropic waves, generated at low latitudes then propagating poleward through a critical layer, thus accelerating low latitudes while decelerating the mid-to-high latitude jet in the late fall through early spring hemisphere. Finally, we identify the dominant waves responsible for the transfers of angular momentum close to northern winter solstice during spin-up and at steady state. Problems with our simulations include peak latitudinal temperature gradients and zonal winds occurring ∼60 km lower than observed by Cassini CIRS, and no reduction in zonal wind speed around 80 km, as was observed by Huygens. While the latter may have been due to transient effects (e.g. gravity waves), the former suggests that our low (∼420 km) model top is adversely affecting the circulation near the jet peak, and/or that we require active haze transport in order to correctly model heating rates and thus the circulation. Future work will include running the model with a higher top, and including advection of a haze particle size distribution.     

Hemispherical distribution of CO above the Venus’ clouds by ground-based 2.3 μm spectroscopy

The abundance of carbon monoxide in the Venus’ dayside atmosphere above the clouds was measured by ground-based 2.3 μm spectroscopy for 4 days. The hemispherical distributions found show no significant latitudinal or longitudinal structure. The disc-averaged mixing ratio of 58 ± 17 ppm found at a representative height of 62-67 km is consistent with previous measurements. Such a flat distribution of CO abundance above the clouds seems to be controlled by an efficient horizontal eddy diffusion with a time scale of 30 days or shorter although the CO distribution below the clouds seems to be controlled by the meridional circulation. The pole-ward wind speed of the meridional circulation above the clouds is estimated to be 0.2 m s⁻¹ or less based on the difference between the CO mixing ratios above and below the clouds.     

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.     

Changes in Jupiter’s zonal velocity between 1979 and 2008

We show that the peak velocity of Jupiter’s visible-cloud-level zonal winds near 24°N (planetographic) increased from 2000 to 2008. This increase was the only change in the zonal velocity from 2000 to 2008 for latitudes between ±70° that was statistically significant and not obviously associated with visible weather. We present the first automated retrieval of fast (∼130 m s⁻¹) zonal velocities at 8°N planetographic latitude, and show that some previous retrievals incorrectly found slower zonal winds because the eastward drift of the dark projections (associated with 5-μm hot spots) “fooled” the retrieval algorithms.We determined the zonal velocity in 2000 from Cassini images from NASA’s Planetary Data System using a global method similar to previous longitude-shifting correlation methods used by others, and a new local method based on the longitudinal average of the two-dimensional velocity field. We obtained global velocities from images acquired in May 2008 with the Wide Field Planetary Camera 2 (WFPC2) on the Hubble Space Telescope (HST). Longer-term variability of the zonal winds is based on comparisons with published velocities based on 1979 Voyager 2 and 1995-1998 HST images. Fluctuations in the zonal wind speeds on the order of 10 m s⁻¹ on timescales ranging from weeks to months were found in the 1979 Voyager 2 and the 1995-1998 HST velocities. In data separated by 10 h, we find that the east-west velocity uncertainty due to longitudinal fluctuations are nearly 10 m s⁻¹, so velocity fluctuations of 10 m s⁻¹ may occur on timescales that are even smaller than 10 h. Fluctuations across such a wide range of timescales limit the accuracy of zonal wind measurements. The concept of an average zonal velocity may be ill-posed, and defining a “temporal mean” zonal velocity as the average of several zonal velocity fields spanning months or years may not be physically meaningful.At 8°N, we use our global method to find peak zonal velocities of ∼110 m s⁻¹ in 2000 and ∼130 m s⁻¹ in 2008. Zonal velocities from 2000 Cassini data produced by our local and global methods agree everywhere, except in the vicinity of 8°N. There, the local algorithm shows that the east-west velocity has large variations in longitude; vast regions exceed ∼140 m s⁻¹. Our global algorithm, and all of the velocity-extraction algorithms used in previously-published studies, found the east-west drift velocities of the visible dark projections, rather than the true zonal velocity at the visible-cloud level. Therefore, the apparent increase in zonal winds between 2000 and 2008 at 8°N is not a true change in zonal velocity.At 7.3°N, the Galileo probe found zonal velocities of 170 m s⁻¹ at the 3-bar level. If the true zonal velocity at the visible-cloud level at this latitude is ∼140 m s⁻¹ rather than ∼105 m s⁻¹, then the vertical zonal wind shear is much less than the currently accepted value.     

Effects of a large convective storm on Saturn's equatorial jet

Hubble Space Telescope observations revealed that Saturn's equatorial jet at the cloud level blows at ∼275 m s⁻¹ today, approximately half the ∼470 m s⁻¹ wind during the Voyager flybys in 1980-1981. Radiative transfer calculations estimate the clouds to be significantly higher today than in 1980. The higher clouds make it difficult to observationally isolate any true slowdown from the vertical wind shear because Voyager and Cassini observations show that the winds become slower with altitude. Here, we test the hypothesis that the large equatorial storm in 1990 called the Great White Spot (GWS) decelerated the equatorial jet. We first use order of magnitude estimates to show: (1) if the GWS triggers vertical momentum redistribution, a minor speed change in the troposphere can lead to a substantial stratospheric wind speed change; (2) storm-triggered turbulent mixing slows a prograde equatorial jet; and (3) a prograde equatorial jet inhibits turbulent mixing in latitude. To test whether a GWS-like large storm decelerates the equatorial jet, we perform numerical experiments using the Explicit Planetary Isentropic Coordinate (EPIC) atmosphere model. Our simulation results are consistent with our order of magnitude predictions. We show that the storm excites waves, and the waves transport westward momentum from the troposphere to the stratosphere and decelerate the equatorial jet by as much as ∼40 m s⁻¹ at the 10-mbar level. However, our results show that the storm's effect is too weak at the cloud levels to halve the jet's speed from ∼470 m s⁻¹. Our results suggest that a combination of higher clouds and a true slowdown is necessary to explain the apparent equatorial jet slowdown. We also analyze the effect of waves on the apparent cloud motions, and show that waves can influence cloud-tracking wind speed measurements.     

Experimental determination of the coefficient of restitution for meter-scale granite spheres

We present results of a series of large-scale experiments to measure the coefficient of restitution for 1-m-diameter rocky bodies in impacts with collision speeds up to ∼1.5 m s⁻¹. The experiments were conducted in an outdoor setting, with two 40-ton cranes used to suspend the ∼1300-kg granite spheres pendulum-style in mutual contact at the bottoms of their respective paths of motion. The spheres were displaced up to ∼1 m from their rest positions and allowed to impact each other in normal-incidence collisions at relative speeds up to ∼1.5 m s⁻¹. Video data from 66 normal-incidence impacts suggest a value for the coefficient of restitution of 0.83 ± 0.06 for collisions between ∼1-m-scale spheres at speeds of order 1 m s⁻¹. No clear trend of coefficient of restitution with impact speed is discernable in the data.     

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.     

In situ observation of penetration process in silica aerogel: Deceleration mechanism of hard spherical projectiles

A large number of cometary dust particles were captured with low-density silica aerogels by NASA’s Stardust Mission. Knowledge of the details of the capture mechanism of hypervelocity particles in silica aerogel is needed in order to correctly derive the original particle features from impact tracks. However, the mechanism has not been fully understood yet. We shot hard spherical projectiles of several different materials into silica aerogel of density 60 mg cm⁻³ and observed their penetration processes using an image converter or a high-speed video camera. In order to observe the deceleration of projectiles clearly, we carried out impact experiments at two velocity ranges; ∼4 km s⁻¹ and ∼200 m s⁻¹. From the movies we took, it was indicated that the projectiles were decelerated by hydrodynamic force which was proportional to v² (v: projectile velocity) during the faster penetration process (∼4 km s⁻¹) and they were merely overcoming the aerogel crushing strength during the slower penetration process (∼200 m s⁻¹). We applied these deceleration mechanisms for whole capture process to calculate the track length. Our model well explains the track length in the experimental data set by Burchell et al. (Burchell, M.J., Creighton, J.A., Cole, M.J., Mann, J., Kearsley, A.T. [2001]. Meteorit. Planet. Sci. 36, 209-221).     

In situ measurements of particle load and transport in dust devils

In situ (mobile) sampling of 33 natural dust devil vortices reveals very high total suspended particle (TSP) mean values of 296 mg m⁻³ and fine dust loadings (PM10) mean values ranging from 15.1 to 43.8 mg m⁻³ (milligrams per cubic meter). Concurrent three-dimensional wind profiles show mean tangential rotation of 12.3 m s⁻¹ and vertical uplift of 2.7 m s⁻¹ driving mean vertical TSP flux of 1689 mg m⁻³ s⁻¹ and fine particle flux of ∼1.0 to ∼50 mg m⁻³ s⁻¹. Peak PM10 dust loading and flux within the dust column are three times greater than mean values, suggesting previous estimates of dust devil flux might be too high. We find that deflation rates caused by dust devil erosion are ∼2.5-50 μm per year in dust devil active zones on Earth. Similar values are expected for Mars, and may be more significant there where competing erosional mechanisms are less likely.     

Reconstruction of the trajectory of the Huygens probe using the Huygens Atmospheric Structure Instrument (HASI)

The Huygens probe returned scientific measurements from the atmosphere and surface of Titan on 14 January 2005. Knowledge of the trajectory of Huygens is necessary for scientific analysis of those measurements. We use measurements from the Huygens Atmospheric Structure Instrument (HASI) to reconstruct the trajectory of Huygens during its mission. The HASI Accelerometer subsystem measured the axial acceleration of the probe with errors of 3E−6 m s⁻², the most accurate measurements ever made by an atmospheric structure instrument on another planetary body. The atmosphere was detected at an altitude of 1498 km. Measurements of the normal acceleration of the probe, which are important for determining the probe's attitude during hypersonic entry, were significantly less accurate and limited by transverse sensitivity of the piezo sensors. Peak acceleration of 121.2 m s⁻² occurred at 234.9 km altitude. The parachute deployment sequence started at 157.1 km and a speed of 342.1 m s⁻¹. Direct measurements of pressure and temperature began shortly afterwards. The measured accelerations and equations of motion have been used to reconstruct the trajectory prior to parachute deployment. Measured pressures and temperatures, together with the equation of hydrostatic equilibrium and the equation of state, have been used to reconstruct the trajectory after parachute deployment. Uncertainties in the entry state of Huygens at the top of the atmosphere are significant, but can be reduced by requiring that the trajectory and atmospheric properties be continuous at parachute deployment.     

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.     

Spatially-resolved high-resolution spectroscopy of Venus 1. Variations of CO2, CO, HF, and HCl at the cloud tops

Variations of the upper cloud boundary and the CO, HF, and HCl mixing ratios were observed using the CSHELL spectrograph at NASA IRTF. The observations were made in three sessions (October 2007, January 2009, and June 2009) at early morning and late afternoon on Venus in the latitude range of ±60°. CO₂ lines at 2.25 μm reveal variations of the cloud aerosol density (∼25%) and scale height near 65 km. The measured reflectivity of Venus at low latitudes is 0.7 at 2.25 μm and 0.028 at 3.66 μm, and the effective CO₂ column density is smaller at 3.66 μm than those at 2.25 μm by a factor of 4. This agrees with the almost conservative multiple scattering at 2.25 μm and single scattering in the almost black aerosol at 3.66 μm. The expected difference is just a factor of (1 − g)⁻¹ = 4, where g = 0.75 is the scattering asymmetry factor for Venus’ clouds. The observed CO mixing ratio is 52 ± 4 ppm near 08:00 and 40 ± 4 ppm near 16:30 at 68 km, and the higher ratio in the morning may be caused by extension of the CO morningside bulge to the cloud tops. The observed weak limb brightening in CO indicates an increase of the CO mixing ratio with altitude. HF is constant at 3.5 ± 0.2 ppb at 68 km in both morningside and afternoon observations and in the latitude range ±60°. Therefore the observations do not favor a bulge of HF, though HF is lighter than CO. Probably a source in the upper atmosphere facilitates the bulge formation. The recent measurements of HCl near 70 km are controversial (0.1 and 0.74 ppm) and require either a strong sink or a strong source of HCl in the clouds. The HCl lines of the (2-0) band are blended by the solar and telluric lines. Therefore we observed the P8 lines of the (1-0) band at 3.44 μm. These lines are spectrally clean and result in the HCl mixing ratio of 0.40 ± 0.03 ppm at 74 km. HCl does not vary with latitude within ±60°. Our observations support a uniformly mixed HCl throughout the Venus atmosphere.     

Ion flow and momentum transfer in the Venus plasma environment

An analysis of ion data from 390 Venus Express, VEX, orbits demonstrates that the flow of solar wind- and ionospheric ions near Venus is characterized by a marked asymmetry. The flow asymmetry of solar wind H⁺ and ionospheric O⁺ points steadily in the opposite direction to the planet’s orbital motion, and is most pronounced near the Pole and in the tail/nightside region. The flow asymmetry is consistent with aberration forcing, here defined as lateral forcing induced by the planet’s orbital motion. In addition to solar wind forcing by the radial solar wind expansion, Venus is also subject a lateral/aberration forcing induced by the planet’s orbital motion transverse to the solar wind flow.The ionospheric response to lateral solar wind forcing is analyzed from altitude profiles of the ion density, ion velocity and ion mass-flux. The close connection between decreasing solar wind H⁺ mass-flux and increasing ionospheric O⁺ mass-flux, is suggestive of a direct/local solar wind energy and momentum transfer to ionospheric plasma. The bulk O⁺ ion flow is accelerated to velocities less than 10 km/s inside the dayside/flank Ionopause, and up to 6000 km in the tail. Consequently, the bulk O⁺ outflow does not escape, but remains near Venus as a fast (km/s) O⁺ zonal wind in the Venus polar and nightside upper ionosphere. Furthermore, the total O⁺ mass-flux in the Venus induced magnetosphere, increases steadily downward to a maximum of 2 × 10⁻¹⁴ kg/(m² s) at ≈400 km altitude, suggesting a downward transport of energy and momentum. The O⁺, and total mass-flux, decay rapidly below 400 km. With no other plasma mass-flux as replacement, we argue that the reduction of ion mass-flux is caused by ion-neutral drag, a transfer of ion energy and momentum to neutrals, implying that the O⁺ plasma wind is converted to a neutral (thermosphere) wind at Venus. Incidentally, such a neutral wind would go in the same direction as the Venus atmosphere superrotation.     

Scale-dependent infrared radiative damping rates on Mars and their role in the deposition of gravity-wave momentum flux

Using a Curtis-matrix model of 15 μm CO₂ radiative cooling rates for the martian atmosphere, we have computed vertical scale-dependent IR radiative damping rates from 0 to 200 km altitude over a broad band of vertical wavenumbers ∣m∣ = 2π(1-500 km)⁻¹ for representative meteorological conditions at 40°N and average levels of solar activity and dust loading. In the middle atmosphere, infrared (IR) radiative damping rates increase with decreasing vertical scale and peak in excess of 30 days⁻¹ at ∼50-80 km altitude, before gradually transitioning to scale-independent rates above ∼100 km due to breakdown of local thermodynamic equilibrium. We incorporate these computed IR radiative damping rates into a linear anelastic gravity-wave model to assess the impact of IR radiative damping, relative to wave breaking and molecular viscosity, in the dissipation of gravity-wave momentum flux. The model results indicate that IR radiative damping is the dominant process in dissipating gravity-wave momentum fluxes at ∼0-50 km altitude, and is the dominant process at all altitudes for gravity waves with vertical wavelengths ?10-15 km. Wave breaking becomes dominant at higher altitudes only for “fast” waves of short horizontal and long vertical wavelengths. Molecular viscosity plays a negligible role in overall momentum flux deposition. Our results provide compelling evidence that IR radiative damping is a major, and often dominant physical process controlling the dissipation of gravity-wave momentum fluxes on Mars, and therefore should be incorporated into future parameterizations of gravity-wave drag within Mars GCMs. Lookup tables for doing so, based on the current computations, are provided.     

The rate of granule ripple movement on Earth and Mars

The rate of movement for 3- and 10-cm-high granule ripples was documented in September of 2006 at Great Sand Dunes National Park and Preserve during a particularly strong wind event. Impact creep induced by saltating sand caused ∼24 granules min⁻¹ to cross each cm of crest length during wind that averaged ∼9 m s⁻¹ (at a height well above 1 m), which is substantially larger than the threshold for saltation of sand. Extension of this documented granule movement rate to Mars suggests that a 25-cm-high granule ripple should require from hundreds to thousands of Earth-years to move 1 cm under present atmospheric conditions.     

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

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

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.