The Martian Hadley circulation: Comparison of “viscous” model predictions to observations

The predictions of an analytical steady thermally forced viscous model for the zonally averaged global circulation of the Martian atmosphere are compared to observations of atmospheric temperatures and eolian feature directions. The temperature of the winter polar atmosphere is a sensitive indicator of the value of the mean eddy diffusivity v of the atmosphere. For global dust storm season, the observed large warming of the atmosphere over the winter (north) pole during the 1977b storm as well as the distribution and azimuths of Type I(b) dust streaks can be well reproduced if v = 10⁷−10⁸cm²sec⁻¹. Observed northeasterly surface winds at high northern latitudes could be produced by plausible fluxes of CO₂ from the northern polar cap. The thermal structure at southern latitudes during this storm is satisfactorily produced, but a major discrepancy exists between observations and predictions of the temperature structure at middle northern latitudes. By comparing these results with those from previous inviscid models of the zonal mean circulation on Mars, we suggest v > 10⁷cm²sec⁻¹ at heights above 30 km and v ⪡ 10⁷cm²sec⁻¹ at lower altitudes may yield a meridional circulation which produces the dynamical warming necessary to account for the observed thermal structure at 0.6 mbar at all latitudes. For late southern summer when the dust content of the atmosphere has decreased to levels characteristic of much of the Martian year, the lack of dynamical warming of the winter pole implies v ≤ 10⁶cm²sec⁻¹. In addition, the predicted surface wind azimuths at this season agree well with the observed azimuths of Type I(d) dust streaks, which are observed to form at this season. The observed latitudinal concentration of the streaks is not explained by the model and may require additional processes which have not been included.     

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.     

Mars surface and atmospheric temperature during the 2001 global dust storm

The dramatic growth and evolution of the 2001 martian global dust storm were captured using the Submillimeter Wave Astronomy Satellite (SWAS). While the lower and middle atmosphere (pressures greater than 50 μbar, up to ∼45 km altitude) showed rapid heating of up to 40 K, the average surface brightness temperature plummeted by ∼20 K at the peak of the storm. The storm appears to have had little impact on the global temperature structure at altitudes above ∼ 10 μbar (∼ 60 km).     

A numerical model of the Martian polar cap winds

An investigation of the Martian polar cap winds and their response to a variety of factors is carried out by a series of numerical experiments based on a zonally symmetric primitive equation model. These factors are the seasonal thermal forcing, mass exchange between polar caps and atmosphere, large-scale topography, and polar cap size. The thermal forcing sets up a circulation whose surface winds adjust to achieve angular momentum balance, with low-latitude easterlies and high-latitude westerlies. The maximum westerlies occur roughly where the horizontal temperature gradients are largest. This pattern changes when cap and atmosphere exchange mass. Corriolis forces acting on the net outflow or inflow produce easterlies at the surface during spring (outflow) and westerlies during winter (inflow). Topography appears to have a small effect, but cap size does play a role, the circulation intensity increasing with cap size. Peak surface winds occur when outflow or inflow is a maximum and are 20 m sec^?1 during spring and 30 m sec^?1 during winter for the northern hemisphere. The model results show that surface winds near the edge of a retreating polar cap are substantially enhanced, a result which is consistent with the Viking observations of local dust storm activity near the edge of the south polar cap during spring. The results also indicate that the surficial wind indicators near the south pole are formed during spring and those near the north pole during winter. The implication is that the high-latitude dune fields in the northern hemisphere are formed at a time when the terrain is being covered with frost. It is therefore suggested that the saltating particles are “snowflakes” which have formed by the mechanism proposed by Pollack etal. The model results for the winter simulation, which have formed by the mechanism transport by large-scale eddies, compare favorably with general circulation model (GCM) calculations. This suggests that the eddy transports may be less important than those associated with the net mass flow, and that 2-D climate modeling may be more succesful for Mars than Earth.     

Circulation of the Venus upper mesosphere/lower thermosphere: Doppler wind measurements from 2001–2009 inferior conjunction,sub-millimeter CO absorption line observations

Sub-millimeter ¹²CO (346 GHz) and ¹³CO (330 GHz) line absorptions, formed within the mesospheric to lower thermospheric altitude (70–120 km) region of the Venus atmosphere, have been mapped across the nightside disk of Venus during 2001–2009 inferior conjunctions, employing the James Clerk Maxwell Telescope (JCMT). Radiative transfer analysis of these thermal line absorptions supports temperature and CO mixing profile retrievals, as described in a companion paper (Clancy et al., 2012). Here, we consider the analysis of the sharp line absorption cores of these CO spectra in terms of accurate Doppler wind profile measurements at 95–115 km altitudes versus local time (～8 pm–4 am) and latitude (～60N–60S). These Doppler wind measurements support determinations of the nightside zonal and subsolar-to-antisolar (SSAS) circulation components over a variety of timescales. The average behavior fitted from 21 retrieved maps of ¹²CO Doppler winds (obtained over hourly, daily, weekly, and interannual intervals) indicates stronger average zonal (85 m/s retrograde) versus SSAS (65 m/s) circulation at the 1 μbar pressure (108–110 km altitude) level. However, the absolute and relative magnitudes of these circulation components exhibit extreme variability over daily to weekly timescales. Furthermore, the individual Doppler wind measurements within each nightside mapping observation generally show significant deviations (20–50 m/s, averaged over 5000 km horizontal scales) from the simple zonal/SSAS solution, with distinct local time and latitudinal characters that are also time variable. These large scale residual circulations contribute 30–70% of the observed nightside Doppler winds at any given time, and may be most responsible for global variations in nightside lower thermospheric trace composition and temperatures, as coincidentally retrieved CO abundance and temperature distributions do not correlate with solution retrograde zonal and SSAS winds (see companion paper, Clancy et al., 2012). Limited comparisons of these nightside submillimeter results with dayside infrared Doppler wind measurements suggest distinct dayside versus nightside circulations, in terms of zonal winds in particular. Combined ¹²CO and ¹³CO Doppler wind mapping observations obtained since 2004 indicate that the average zonal and SSAS wind components increase by 50–100% between altitudes of 100 and 115 km. If gravity waves originating from the cloud levels are responsible for the extension of zonal winds into the thermosphere (Alexander, M.J. [1992]. Geophys. Res. Lett. 19, 2207–2210), such waves deposit substantial momentum (i.e., break) in the lower nightside thermosphere.     

Solar infrared occultation observations by SPICAM experiment on Mars-Express: Simultaneous measurements of the vertical distributions of H2O, CO2 and aerosol

The infrared AOTF spectrometer is a part of the SPICAM experiment onboard the Mars-Express ESA mission. The instrument has a capability of solar occultations and operates in the spectral range of 1-1.7 μm with a spectral resolution of ∼3.5 cm⁻¹. We report results from 24 orbits obtained during MY28 at Ls 130°-160°, and the latitude range of 40°-55° N. For these orbits the atmospheric density from 1.43 μm CO₂ band, water vapor mixing ratio based on 1.38 μm absorption, and aerosol opacities were retrieved simultaneously. The vertical resolution of measurements is better than 3.5 km. Aerosol vertical extinction profiles were obtained at 10 wavelengths in the altitude range from 10 to 60 km. The interpretation using Mie scattering theory with adopted refraction indices of dust and H₂O ice allows to retrieve particle size (reff∼0.5-1 μm) and number density (∼1 cm⁻³ at 15-30 km) profiles. The haze top is generally below 40 km, except the longitude range of 320°-50° E, where high-altitude clouds at 50-60 km were detected. Optical properties of these clouds are compatible with ice particles (effective radius reff=0.1-0.3 μm, number density N∼10 cm−3) distributed with variance νeff=0.1-0.2 μm. The vertical optical depth of the clouds is below 0.001 at 1 μm. The atmospheric density profiles are retrieved from CO₂ band in the altitude range of 10-90 km, and H₂O mixing ratio is determined at 15-50 km. Unless a supersaturation of the water vapor occurs in the martian atmosphere, the H₂O mixing ratio indicates ∼5 K warmer atmosphere at 25-45 km than predicted by models.     

A reanalysis of Venus winds at two cloud levels from Galileo SSI images

The global circulation of the Venus atmosphere is characterized at cloud level by a zonal super rotation studied over the years with data from a battery of spacecrafts: orbiters, balloons and probes. Among them, the Galileo spacecraft monitored the Venus atmosphere in a flyby in February 1990 in its route toward Jupiter. Since the flyby was almost equatorial, published analysis of zonal winds obtained from displacements of cloud elements on images obtained by the SSI camera [Belton, M.J.S., and 20 colleagues, 1991. Science 253, 1531-1536] stop at latitudes 50° north and south. In this paper we present new results on Venus winds based on a reanalysis of an extended set of images obtained at two wavelengths, 418 nm (violet) and 986 nm (near infrared), that sense different altitude levels in the upper cloud. Our main result is that we have been able to extend the zonal wind profile up to the polar latitudes: 70° N and 70° S at 418 nm and 70° N at 986 nm. Binned and smoothed profiles are given in tabular form. We show that the zonal winds drop in their velocity poleward of latitudes 45° N and 50° S where an intense meridional wind shear develops at the two cloud levels. Our data confirm the magnitude of this shear, retrieved previously from radio occultation data, but disagrees with it in the latitudinal location of the sheared region. The new wind data can be used to recalibrate the zonal winds retrieved from the previous measurements of the temperature field and the cyclostrophic balance assumption. The meridional profiles of the zonal winds at the two cloud levels are used to assess the vertical wind shear in the upper cloud layer as a function of latitude and locate the most unstable region.     

Winds and tidal oscillations in the upper middle atmosphere at Saskatoon (52°N, 107°W, L = 4.3) during the year June 1982–May 1983

The M.F. radar (2.2 MHz) operating at Saskatoon, Canada (52°N, 107°W) has been used to obtain continuous winds data for 1982/3: heights are ～60–110 km in daylight hours and ～80–110 km at night. The monthly data are compared with means of prevailing winds, tides (12-, 24-h) and planetary waves for the previous 4 years: 1978–1982. A regular climatology emerges that strongly suggests that the theoretical modelling of the upper middle atmosphere is an achievable goal: also that updated global data models are required.     

Validation of martian meteorological data assimilation for MGS/TES using radio occultation measurements

We describe an assimilation of thermal profiles below about 40 km altitude and total dust opacities into a general circulation model (GCM) of the martian atmosphere. The data were provided by the Thermal Emission Spectrometer (TES) on board the Mars Global Surveyor (MGS) spacecraft. The results of the assimilation are verified against an independent source of contemporaneous data represented by radio occultation measurements with an ultra-stable radio oscillator, also aboard MGS. This paper describes a comparison between temperature profiles retrieved by the radio occultation experiments and the corresponding profiles given by both an independent, carefully tuned GCM simulation and by an assimilation of TES observations performed over the period of time from middle, northern summer in martian year 24, corresponding to May 1999, until late, northern spring in martian year 27, corresponding to August 2004. This study shows that the assimilation of TES measurements improves the overall agreement between radio occultation observations and the GCM analysis, in particular below 20 km altitude, where the radio occultation measurements are known to be most accurate. Discrepancies still remain, mostly during the global dust storm of year 2001 and at latitudes around 60° N in northern winter-early spring. These are the periods of time and locations, however, for which discrepancies between TES and radio occultation profiles are also shown to be the largest. Finally, a further direct validation is performed, comparing stationary waves at selected latitudes and time of year. Apart from biases at high latitudes in winter time, data assimilation is able to represent the correct wave behaviour, which is one major objective for martian assimilation.     

Polar warming in the middle atmosphere of Mars

We report ground-based laser heterodyne spectroscopy of non-thermal emission in the cores of the 10.33-μmR(8) and 10.72-μmP(32) lines of ¹²C¹⁶O₂, obtained at 23 locations on the disk of Mars during the 1984 opposition, at L_s = 130°. The data were obtained at a sub-Doppler spectral resolution, and the temperature of the middle Martian atmosphere (50–85 km) is derived from the frequency width and intensity of the R(8) emission, and from the total intensity of the P(32) emission. We find that the temperature of the middle Martian atmosphere varies with latitude. Near the subsolar latitude, the average 50- to 85-km temperature is close to the radiative equilibrium value for a CO₂ atmosphere. However, at high latitudes in both the northern (summer) and southern (winter) hemispheres the 50- to 85-km temperature exceeds the CO₂ radiative equilibrium value; a meridional gradient in the range of 0.4 – 0.9°K per degree of latitude is indicated by our data. The highest temperatures are seen at high latitudes in the winter hemisphere, reminiscent of the seasonal effects seen at the Earth's mesopause. As in the terrestrial case, this winter polar warming in the Martian middle atmosphere necessitates departures from radiative equilibrium; dynamical heating of order 4 × 10² ergs g⁻¹ sec⁻¹ is required at the edge of the winter polar night. A comparison with 2-D circulation models shows that the presence of atmospheric dust may enhance this dynamical heating at high winter latitudes, and may also account for heating at high latitudes in the summer hemisphere.     

Thermospheric changes shortly after the onset of daytime joule heating

In the first few tens of minutes after the onset of widespread Joule heating, the motion of the ionospheric atmosphere can be approximated by the one-dimensional motion of a gas in a gravity field—a problem that is easily solved because the motion takes place at constant pressure. The solution provides an estimate of time for which the model is applicable to the physical situation. Seasonal variations of the early effects are examined by using ion profiles appropriate to each season. The results show that the atmosphere above 100 km is strongly modified within a few tens of minutes after the onset of widespread heating: the density can double, the temperature can increase several hundred degrees, and the molecular nitrogen concentration can quadruple. Vertical winds exceeding 100 m/sec at 400km altitude are possible for a brief period after the onset of electric fields of 100 mV/m—rare but observed events. In the first few tens of minutes after the onset of a given electric field, the greatest power is deposited in the thermosphere around summer solstice, while the greatest winds occur at 200 km altitude in the summer and at 400km in the winter. These differing seasonal effects show primarily that a given level of change occurs sooner for one season than another, not that long term seasonal differences exist. Once a magnetic storm is in progress, the quiet-day ion profiles change to the non-seasonal storm profile ; for this ion distribution, F-region effects are minimum regardless of season. Joule heating effects in the upper thermosphere are therefore concluded to be self-limiting.     

Features of the upper atmosphere revealed by analysis of the orbit of cosmos 1009 rocket, 1978-50B

COSMOS 1009 rocket was launched on 19 May 1978 into an orbit with initial perigee height 150 km and apogee 1100 km: its lifetime was only 17 days. The orbit has been determined daily during the final 14 days of its life, using the RAE orbit refinement program PROP6,with about 1100 observations supplied by NORAD. An average accuracy of about 60 m, radial and cross-track, was achieved.The orbits were analysed to reveal three features of the upper atmosphere at heights between 125 and 175 km. From the decrease in perigee height, five values of density scale height, accurate to ±4%, were obtained. The first three were within 10% of those from CIRA 1972; the fourth, after a magnetic storm, was higher than expected; the fifth gave evidence of the decrease in drag coefficient at heights below 130 km.Atmospheric oblateness produced a change of 4° in perigee position during the last four days of the life. Analysis showed that the ellipticity of the upper atmosphere was approximately equal to that of the Earth, f, for the first two of the four days, and about $\text{1}\text{2}\text{f}$ in the last two.The orbital inclination decreased during the 14 days by about 50 times its standard deviation, and the observed variation was analysed to determine zonal winds at heights of 150–160 km at latitudes near 47° north. The zonal wind was very weak (0±30 m/s) for 23–28 May at local times near 03h; and 90±30 m/s east-to-west for 29 May to 4 June at local times near 01 h.     

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.     

Winter polar warmings and the meridional transport on Mars simulated with a general circulation model

Winter polar warmings in the middle atmosphere of Mars occur due to the adiabatic heating associated with the downward branch of the cross-equatorial meridional circulation. Thus, they are the manifestation of the global meridional transport rather than of local radiative effects. We report on a series of numerical experiments with a recently developed general circulation model of the martian atmosphere to examine the relative roles of the mechanical and thermal forcing in the meridional transport. The experiments were focused on answering the question of whether the martian circulation is consistent with the thermally driven nearly inviscid Hadley cell, as was pointed out by some previous studies, or it is forced mainly by zonally asymmetric eddies. It is demonstrated that, under realistic conditions in the middle atmosphere, the meridional transport is maintained primarily by dissipating large-scale planetary waves and solar tides. This mechanism is similar to the “extratropical pump” in the middle atmosphere on Earth. Only in the run with artificially weak zonal disturbances, was the circulation reminiscent of thermally induced Hadley cells. In the experiment with an imposed dust storm, the modified atmospheric refraction changes the vertical propagation of the eddies. As the result, the Eliassen-Palm fluxes convergence increases in high winter latitudes of the middle atmosphere, the meridional transport gets stronger, and the polar temperature rises. Additional numerical experiments demonstrated that insufficient model resolution, increased numerical dissipation, and, especially, neglect of non-LTE effects for the 15 μm CO₂ band could weaken the meridional transport and the magnitude of polar warmings in GCMs.     

Comparative studies of E-region ionospheric drifts and meteor winds

Ionospheric drifts using total reflections from the E-region have been compared with neutral winds measured by meteor radar. Close agreement was found when both measurements were made in a common volume of atmosphere. Even with a separation of 700 km between the measuring regions the results were very similar. It is concluded that the drift technique does measure the movement of the neutral atmosphere in the altitude range 95–120 km. The agreement between measurements from widely separated regions indicates the horizontal scale of the wind structure is at least 700 km.     

Stellar occultation measurements of molecular oxygen in the lower thermosphere

Stellar ultraviolet light near 1500 Å is attenuated in the Earth's upper atmosphere due to strong absorption in the Schumann-Runge continuum of molecular oxygen. The intensity of stars in the Schumann-Runge continuum region has been monitored by the University of Wisconsin stellar photometers aboard the OAO-2 satellite during occultation of the star by the Earth's atmosphere. These data have been used to determine the molecular oxygen number density profile at the occultation tangent point. The results of 14 stellar occultations obtained in low and middle latitudes are presented giving the night-time vertical number density profile of molecular oxygen in the 140–200 km region. In general, the measured molecular oxygen number density is about a factor of 2 lower than the number densities predicted by the CIRA 1965 model. Also, the number density at a given height appears to decrease with decreasing solar activity. Measurements taken at low latitudes during the August 1970 geomagnetic storm showed a decrease in the molecular oxygen number density at a given height several days after the peak of the storm followed by a slow recovery to pre-storm densities.     

Interannual Variability of Martian Global Dust Storms: Simulations with a Low-Order Model of the General Circulation

We present simulations of the interannual variability of martian global dust storms (GDSs) with a simplified low-order model (LOM) of the general circulation. The simplified model allows one to conduct computationally fast long-term simulations of the martian climate system. The LOM is constructed by Galerkin projection of a 2D (zonally averaged) general circulation model (GCM) onto a truncated set of basis functions. The resulting LOM consists of 12 coupled nonlinear ordinary differential equations describing atmospheric dynamics and dust transport within the Hadley cell. The forcing of the model is described by simplified physics based on Newtonian cooling and Rayleigh friction. The atmosphere and surface are coupled: atmospheric heating depends on the dustiness of the atmosphere, and the surface dust source depends on the strength of the atmospheric winds. Parameters of the model are tuned to fit the output of the NASA AMES GCM and the fit is generally very good.Interannual variability of GDSs is possible in the LOM, but only when stochastic forcing is added to the model. The stochastic forcing could be provided by transient weather systems or some surface process such as redistribution of the sand particles in storm generating zones on the surface. The results are sensitive to the value of the saltation threshold, which hints at a possible feedback between saltation threshold and dust storm activity. According to this hypothesis, erodable material builds up as a result of a local process, whose effect is to lower the saltation threshold until a GDS occurs. The saltation threshold adjusts its value so that dust storms are barely able to occur.     

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.     

Diagnostics of Titan's stratospheric dynamics using Cassini/CIRS data and the 2-dimensional IPSL circulation model

The dynamics of Titan's stratosphere is discussed in this study, based on a comparison between observations by the CIRS instrument on board the Cassini spacecraft, and results of the 2-dimensional circulation model developed at the Institute Pierre-Simon Laplace, available at http://www.lmd.jussieu.fr/titanDbase [Rannou, P., Lebonnois, S., Hourdin, F., Luz, D., 2005. Adv. Space Res. 36, 2194-2198]. The comparison aims at both evaluating the model's capabilities and interpreting the observations concerning: (1) dynamical and thermal structure using temperature retrievals from Cassini/CIRS and the vertical profile of zonal wind at the Huygens landing site obtained by Huygens/DWE; and (2) vertical and latitudinal profiles of stratospheric gases deduced from Cassini/CIRS data. The modeled thermal structure is similar to that inferred from observations (Cassini/CIRS and Earth-based observations). However, the upper stratosphere (above 0.05 mbar) is systematically too hot in the 2D-CM, and therefore the stratopause region is not well represented. This bias may be related to the haze structure and to misrepresented radiative effects in this region, such as the cooling effect of hydrogen cyanide (HCN). The 2D-CM produces a strong atmospheric superrotation, with zonal winds reaching 200 m s⁻¹ at high winter latitudes between 200 and 300 km altitude (0.1-1 mbar). The modeled zonal winds are in good agreement with retrieved wind fields from occultation observations, Cassini/CIRS and Huygens/DWE. Changes to the thermal structure are coupled to changes in the meridional circulation and polar vortex extension, and therefore affect chemical distributions, especially in winter polar regions. When a higher altitude haze production source is used, the resulting modeled meridional circulation is weaker and the vertical and horizontal mixing due to the polar vortex is less extended in latitude. There is an overall good agreement between modeled chemical distributions and observations in equatorial regions. The difference in observed vertical gradients of C₂H₂ and HCN may be an indicator of the relative strength of circulation and chemical loss of HCN. The negative vertical gradient of ethylene in the low stratosphere at 15° S, cannot be modeled with simple 1-dimensional models, where a strong photochemical sink in the middle stratosphere would be necessary. It is explained here by dynamical advection from the winter pole towards the equator in the low stratosphere and by the fact that ethylene does not condense. Near the winter pole (80° N), some compounds (C₄H₂, C₃H₄) exhibit an (interior) minimum in the observed abundance vertical profiles, whereas 2D-CM profiles are well mixed all along the atmospheric column. This minimum can be a diagnostic of the strength of the meridional circulation, and of the spatial extension of the winter polar vortex where strong descending motions are present. In the summer hemisphere, observed stratospheric abundances are uniform in latitude, whereas the model maintains a residual enrichment over the summer pole from the spring cell due to a secondary meridional overturning between 1 and 50 mbar, at latitudes south of 40-50° S. The strength, as well as spatial and temporal extensions of this structure are a difficulty, that may be linked to possible misrepresentation of horizontally mixing processes, due to the restricted 2-dimensional nature of the model. This restriction should also be kept in mind as a possible source of other discrepancies.