Characteristics of planetary-scale waves simulated by a new venusian mesosphere and thermosphere general circulation model

We have developed a new general circulation model (GCM) for the venusian mesosphere and thermosphere (80-about 180 km). Our GCM simulations show that winds in the subsolar-to-antisolar direction (SS–AS) are predominant above about 90 km. A weak return flow of the SS–AS is seen below about 90 km. We performed GCM simulations imposing the planetary-scale waves (thermal tides, Rossby wave, and Kelvin wave) at the lower boundary. Although the diurnal and semidiurnal tides are damped below 95 km, the Rossby wave propagates up to around 130 km. However, the amplitude of the Rossby wave is too small (<1 m/s) to affect the general circulation. On the other hand, the Kelvin wave propagates up to about 130 km with a maximum zonal wind fluctuation of approximately 5.9 m/s on average. The amplitude of the Kelvin wave sometimes exceeds 10 m/s around the terminator. The Kelvin wave causes a temporal variation in the wind velocity at the altitude of the O₂-1.27 μm nightglow emission (about 95 km). Using a newly developed 1-D nightglow model and the composition distribution calculated from our GCM, we investigated the impact of the Kelvin wave on the nightglow distribution. Our results suggest that the Kelvin wave would cause temporal variations in the nightglow emission in the 23:50–00:20 LT region with an intensity of 1.1–1.3 MR and a period of approximately 4 days.     

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.     

Seasonal variation in the atmospheric activity of Jupiter

Periodicities of 22, 8 and 6 years have been found by a statistical analysis of the coefficient of atmospheric activity R of Jupiter for the time period 1910–1985.We have also found variations of the coefficient of asymmetry of the activity Z measured by the Observatory of Athens for the time period 1956–1985 with periodicities of 12, 8, 6 and 4 years by statistical analysis.The same periodicities have been recently found in an analysis of the total number of solar flares and in the number of high velocity solar wind streamers.     

Simulations of Titan's brightness by a two-dimensional haze model

We have used a 2-D microphysics model to study the effects of atmospheric motions on the albedo of Titan's thick haze layer. We compare our results to the observed variations of Titan's brightness with season and latitude. We use two wind fields; the first is a simple pole-to-pole Hadley cell that reverses twice a year. The second is based on the results of a preliminary Titan GCM. Seasonally varying wind fields, with horizontal velocities of about 1 cm sec-1 at optical depth unity, are capable of producing the observed change in geometric albedo of about 10% over the Titan year. Neither of the two wind fields can adequately reproduce the latitudinal distribution of reflectivity seen by Voyager. At visible wavelengths, where only haze opacity is important, upwelling produces darkening by increasing the particle size at optical depth unity. This is due to the suspension of larger particles as well as the lateral removal of smaller particles from the top of the atmosphere. At UV wavelengths and at 0.89 micrometers the albedo is determined by the competing effects of the gas the haze material. Gas is bright in the UV and dark at 0.89 micrometers. Haze transport at high altitudes controls the UV albedo and transport at low altitude controls the 0.89 micrometers albedo. Comparisons between the hemispheric contrast at UV, visible, and IR wavelengths can be diagnostic of the vertical structure of the wind field on Titan.     

Decadal variations in a Venus general circulation model

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

A sensitivity study of changes in Earth's rotation rate with an atmospheric general circulation model

A series of general circulation model simulations in which Earth's rotation rate has been increased is presented using the community climate model version 1 (CCM1) of the National Center for Atmospheric Research (NCAR). The rotation rate has been altered in order to simulate day-lengths of 24-, 22-, 20-, 18-, 16- and 14-h Earth days. This is a plausible range of Earth's day-length during the last 4 billion years. In an earlier study with a simple energy balance ocean, which does not store heat (sometimes referred to as a swamp), it was shown that reducing the day-length to a 14-h day caused a 20% reduction in the global mean cloud fraction. In this study however, using fixed sea surface temperatures (SSTs) with January solar forcing, a slight increase in clouds occurs with faster rotation, although changes in relative humidity are similar to the earlier study. Furthermore, as in an earlier study, there is more sinking in the mid-latitudes. This sinking is most prominent over the Pacific and Atlantic storm tracks, indicating that the baroclinic eddies have been weakened. With faster rotation rates, the storm tracks are defined by shorter waves as compared to the control simulation. A significant change in the large-scale zonally averaged circulation occurs when the day-length is reduced to a values less than 18-hours. The effects of faster rotation rates on stationary eddy heat transport may help to explain high latitude glaciation of the Ordovician some 440 Ma years ago.     

Orbital change experiments with a Mars general circulation model

We use a Mars general circulation model to examine the effect of orbital changes on the planet’s general circulation and climate system. Experiments are performed for obliquities ranging from 0° to 60° for two different longitudes of perihelion. Each experiment simulates a full Mars year assuming a fixed atmospheric dust distribution and fixed amount of CO₂ in the atmosphere/cap system. We find that global mean surface temperatures and pressures decline with increasing obliquity due to the increasing extent of the winter polar caps. The seasonal CO₂ cycle and intensity of the solstice circulation amplify considerably with increasing obliquity such that global dust storms are likely at both solstices. The most significant feature of the high obliquity solstice circulations is the development of an intense low-level jet associated with the return branch of the Hadley circulation.Model surface stresses are used to map regions of preferred dust lifting, which are defined in terms of an annual deflation potential. For the present obliquity, the model-predicted regions of high deflation potential are in good agreement with Cantor et al.’s (2001, J. Geophys. Res.106, 23653-23688) observations, which gives us some confidence in the model’s ability to predict where lifting might occur when Mars’ orbit parameters are different than they are today. In general we find that the dust lifting potential increases sharply with obliquity and is greatest at times of high obliquity when perihelion coincides with northern summer solstice. Over an obliquity cycle, the model global annual deflation potential ranges from several tenths of a millimeter at 0° obliquity to almost 15 mm at 60° obliquity. Much higher values are possible when the atmosphere is very dusty.We find a strong correlation between the deflation potential and surface thermal inertia: regions of high deflation potential correspond to regions of high thermal inertia (high rock abundance), and regions of low deflation potential correspond to regions of low thermal inertia (high dust/sand abundance). Furthermore, while the regions of preferred lifting (high deflation potential) expand somewhat with increasing obliquity and dust loading, the central parts of Tharsis, Arabia, and Elysium show no tendency for significant lifting at any obliquity or longitude of perihelion. These regions may therefore be very old and represent net long-term sinks for atmospheric dust. It is the topography of the planet, through its influence on surface pressure and wind systems, which ultimately determines where dust accumulates.Finally, as was found by Fenton and Richardson (2001, J. Geophys. Res.106, 32885-32909), we find no tendency for the development of east-southeasterly winds at the Pathfinder site for any of our orbital change experiments. This suggests that the ancient wind regime discussed by Greeley et al. (2000, J. Geophys. Res.105, 1829-1840) was produced by other factors, such as polar wander.     

Three-dimensional on-line chemical modeling in a Mars general circulation model

In this paper we present the design and some results of an on-line chemical model of the atmosphere of Mars using a three-dimensional dynamical-chemical model. A chemistry module has been added to the 3D Mars Global Multiscale Model (GM3), a general circulation model based on the Canadian weather forecast model. As water is a major fuel for chemistry on Mars, we have taken care to simulate the atmospheric water vapor content, taking into account evaporation from and condensation to the polar caps and in the atmosphere we allow for freezing and formation of a bulk cloud and sedimentation. The chemical module includes the main neutral carbon-oxygen-hydrogen-nitrogen species and associated relevant reactions. The outputs of the model are analyzed in the light of our current understanding of the chemistry on Mars. The results are in general accordance with 1D models and with two other 3D models with more limited chemistry. The model does not require elevated diffusion coefficients to fit observational data in the thermosphere as was suggested in some 1D studies. The large scale transport and the downward thermospheric transport in the winter pole appear to have a significant role in the concentrations of some trace species.     

Seasonal variation of Titan's stratospheric ethylene (C2H4) observed

All previous observations of seasonal change on Titan have been of physical phenomena such as clouds and haze. We present here the first observational evidence of chemical change in Titan's atmosphere. Images taken during 1999-2002 (late southern spring on Titan) with the W.M. Keck I 10-meter telescope at 8-13 μm show a significant accumulation of ethylene (C₂H₄) in the south polar stratosphere as well as north-south stratospheric temperature variation (colder at poles). Our observations restrict this newly discovered south polar ethylene accumulation to latitudes south of 60° S. The only other observations of the spatial distribution of C₂H₄ were those of Voyager I, which found a significant north polar accumulation in early northern spring. We see no build-up in the north, although the highest northern latitudes are obstructed from view in the current season. Our observations constrain any unobserved north polar accumulation of C₂H₄ to north of 50° N latitude. Comparison of the Voyager I results with our new results show seasonal chemical change has occurred in Titan's atmosphere.     

The deep wind structure of the giant planets: Results from an anelastic general circulation model

The giant gas planets have hot convective interiors, and therefore a common assumption is that these deep atmospheres are close to a barotropic state. Here we show using a new anelastic general circulation model that baroclinic vorticity contributions are not negligible, and drive the system away from an isentropic and therefore barotropic state. The motion is still aligned with the direction of the axis of rotation as in a barotropic rotating fluid, but the wind structure has a vertical shear with stronger winds in the atmosphere than in the interior. This shear is associated with baroclinic compressibility effects. Most previous convection models of giant planets have used the Boussinesq approximation, which assumes the density is constant in depth; however, Jupiter's actual density varies by four orders of magnitude through its deep molecular envelope. We therefore developed a new general circulation model (based on the MITgcm) that is anelastic and thereby incorporates this density variation. The model's geometry is a full 3D sphere down to a small inner core. It is nonhydrostatic, uses an equation of state suitable for hydrogen-helium mixtures (SCVH), and is driven by an internal heating profile. We demonstrate the effect of compressibility by comparing anelastic and Boussinesq cases. The simulations develop a mean state that is geostrophic and hydrostatic including the often neglected, but significant, vertical Coriolis contribution. This leads to modification of the standard thermal wind relation for a deep compressible atmosphere. The interior flow organizes in large cyclonically rotating columnar eddies parallel to the rotation axis, which drive upgradient angular momentum eddy fluxes, generating the observed equatorial superrotation. Heat fluxes align with the axis of rotation, and provide a mechanism for the transport of heat poleward, which can cause the observed flat meridional emission. We address the issue of over-forcing which is common in such convection models and analyze the dependence of our results on this; showing that the vertical wind structure is not very sensitive to the Rayleigh number. We also study the effect of rotation, showing how the transition from a rapidly to a slowly rotating system affects the dynamics.     

A physical model of Titan's aerosols

Microphysical simulations of Titan's stratospheric haze show that aerosol microphysics is linked to organized dynamical processes. The detached haze layer may be a manifestation of 1 cm sec-1 vertical velocities at altitudes above 300 km. The hemispherical asymmetry in the visible albedo may be caused by 0.05 cm sec-1 vertical velocities at altitudes of 150 to 200 km, we predict contrast reversal beyond 0.6 micrometer. Tomasko and Smith's (1982, Icarus 51, 65-95) model, in which a layer of large particles above 220 km altitude is responsible for the high forward scattering observed by Rages and Pollack (1983, Icarus 55, 50-62), is a natural outcome of the detached haze layer being produced by rising motions if aerosol mass production occurs primarily below the detached haze layer. The aerosol's electrical charge is critical for the particle size and optical depth of the haze. The geometric albedo, particularly in the ultraviolet and near infrared, requires that the particle size be near 0.15 micrometer down to altitudes below 100 km, which is consistent with polarization observations (Tomasko and Smith 1982, West and Smith 1991, Icarus 90, 330-333). Above about 400 km and below about 150 km Yung et al.'s (1984, Astrophys. J. Suppl. Ser. 55, 465-506) diffusion coefficients are too small. Dynamical processes control the haze particles below about 150 km. The relatively large eddy diffusion coefficients in the lower stratosphere result in a vertically extensive region with nonuniform mixing ratios of condensable gases, so that most hydrocarbons may condense very near the tropopause rather than tens of kilometers above it. The optical depths of hydrocarbon clouds are probably less than one, requiring that abundant gases such as ethane condense on a subset of the haze particles to create relatively large, rapidly removed particles. The wavelength dependence of the optical radius is calculated for use in analyzing observations of the geometric albedo. The lower atmosphere and surface should be visible outside of regions of methane absorption in the near infrared. Limb scans at 2.0 micrometers wavelength should be possible down to about 75 km altitude.     

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.     

Near-surface winds at the Huygens site on Titan: Interpretation by means of a general circulation model

This study aims at interpreting the zonal and meridional wind in Titan's troposphere measured by the Huygens probe by means of a general circulation model. The numerical simulation elucidates the relative importance of the seasonal variation in the Hadley circulation and Saturn's gravitational tide in affecting the actual wind profile. The observed reversal of the zonal wind at two altitudes in the lower troposphere can be reproduced with this model only if the near-surface temperature profile is asymmetric about the equator and substantial seasonal redistribution of angular momentum by the variable Hadley circulation takes place. The meridional wind near the surface is mainly caused by the meridional pressure gradient and is thus a manifestation of the Hadley circulation. Southward meridional wind in the PBL (planetary boundary layer) is consistent with the near-surface temperature at the equator being lower than at mid southern latitudes. Even small changes in the radiative heating profile in the troposphere can substantially affect the mean zonal and meridional wind including their direction. Saturn's gravitational tide is rather weak at the Huygens site due to the proximity to the equator, and does not clearly manifest itself in the instantaneous vertical profile of wind. Nevertheless, the simulated descent trajectory is more consistent with the observation if the tide is present. Because of a different force balance in Titan's atmosphere from terrestrial conditions, PBL-specific wind systems like on Earth are unlikely to exist on Titan.     

A coupled dynamics-microphysics model of Titan's atmosphere

We have developed a coupled general circulation model of Titan's atmosphere in which the aerosol haze is treated with a microphysical model and is advected by the winds. The radiative transfer accounts for the non uniform haze distribution and, in turn, drives the dynamics. We analyze the GCM results, especially focusing on the difference between a uniform haze layer and a haze layer coupled to the dynamics. In the coupled simulation the aerosols tend to accumulate at the poles, at latitudes higher than ±60°. During winter, aerosols strongly radiate at thermal infrared wavelengths enhancing the cooling rate near the pole. Since this tends to increase the latitudinal gradients of temperature the direct effect of this cooling excess, in contrast to the uncoupled haze case, is to increase the strength of the meridional cells as well as the strength of the zonal winds and profile. This is a positive feedback of the haze on dynamics. The coupled model reproduces observations about the state of the atmosphere better than the uniform haze model, and in addition, the northern polar hood and the detached haze are qualitatively reproduced.     

Seasonal evolution of Titan's dark polar hood: midsummer disappearance observed by the Hubble Space Telescope

Titan, Saturn's largest moon, has a dense organic-laden atmosphere that displays dramatic seasonal variations in composition and appearance. Here we document the evolution of the dark polar hood, first seen in 1980 by Voyager 1 around the north pole, and report quantitative measurements of the hood's disappearance from the south pole in 2002–2003 using previously unpublished observations with the Hubble Space Telescope Advanced Camera for Surveys ( HST /ACS). These data support a model of the hood as a transient structure associated with downwelling during polar winter.     

A model of Titan's haze of fractal aerosols constrained by multiple observations

We use Titan's geometric albedo to constrain the vertical distribution of the haze. Microphysical models incorporating fractal aggregates do not readily fit the methane features at 0.62 μm band and the dark 0.88 μm of the albedo spectrum simultaneously. We take advantage of this apparent discrepancy to constrain the haze vertical profile.We used the geometric albedo and several results and constraints from other works to better constrain the vertical haze extinction profile, especially in the low stratosphere. The objective of this model is to give a solution that simultaneously fits the main constraints known to apply to the haze.We find that the haze extinction increases with decreasing altitude with a scale height about equal to the atmospheric scale height down to 100 km. Below this altitude, extinction must decrease down to 30 km. This is necessary in order to have enough haze to sustain a relatively high albedo (0.076) in the dark 0.88 μm methane band and to show the 0.62 μm band in the haze continuum. We set the haze production rate around 7×10⁻¹⁴ kgm⁻² s⁻¹, and the aerosols production altitude around 400 km (or at pressure 1.5 Pa).The physical processes which generate such a profile are not clear. However, purely one-dimensional effects such as condensation, sedimentation, and rainout can be ruled out, and we believe that this relative clearing in Titan's troposphere and lower stratosphere is due to particle horizontal transport by the mean circulation.     

A photochemical model of Titan's atmosphere and ionosphere

A global-mean model of coupled neutral and ion chemistry on Titan has been developed. Unlike the previous coupled models, the model involves ambipolar diffusion and escape of ions, hydrodynamic escape of light species, and calculates the H₂ and CO densities near the surface that were assigned in some previous models. We tried to reduce the numbers of species and reactions in the model and remove all species and reactions that weakly affect the observed species. Hydrocarbon chemistry is extended to C₁₂H₁₀ for neutrals and C₁₀H⁺₁₁ for ions but does not include PAHs. The model involves 415 reactions of 83 neutrals and 33 ions, effects of magnetospheric electrons, protons, and cosmic rays. UV absorption by Titan's haze was calculated using the Huygens observations and a code for the aggregate particles. Hydrocarbon, nitrile, and ion chemistries are strongly coupled on Titan, and attempt to calculate them separately (e.g., in models of ionospheric composition) may result in significant error. The model densities of various species are typically in good agreement with the observations except vertical profiles in the stratosphere that are steeper than the CIRS limb data. (A model with eddy diffusion that facilitates fitting to the CIRS limb data is considered as well.) The CO densities are supported by the O⁺ flux from Saturn's magnetosphere. The ionosphere includes a peak at 80 km formed by the cosmic rays, steplike layers at 500-700 and 700-900 km and a peak at 1060 km (SZA = 60°). Nighttime densities of major ions agree with the INMS data. Ion chemistry dominates in the production of bicyclic aromatic hydrocarbons above 600 km. The model estimates of heavy positive and negative ions are in reasonable agreement with the Cassini results. The major haze production is in the reactions C₆H + C₄H₂, C₃N + C₄H₂, and condensation of hydrocarbons below 100 km. Overall, precipitation rate of the photochemical products is equal to 4-7 kg cm⁻² Byr⁻¹ (50-90 m Byr⁻¹ while the global-mean depth of the organic sediments is ∼3 m). Escape rates of methane and hydrogen are 2.9 and 1.4 kg cm⁻² Byr⁻¹, respectively. The model does not support the low C/N ratio observed by the Huygens ACP in Titan's haze.     

The water cycle in the general circulation model of the martian atmosphere

Within the numerical general-circulation model of the Martian atmosphere MAOAM (Martian Atmosphere: Observation and Modeling), we have developed the water cycle block, which is an essential component of modern general circulation models of the Martian atmosphere. The MAOAM model has a spectral dynamic core and successfully predicts the temperature regime on Mars through the use of physical parameterizations typical of both terrestrial and Martian models. We have achieved stable computation for three Martian years, while maintaining a conservative advection scheme taking into account the water–ice phase transitions, water exchange between the atmosphere and surface, and corrections for the vertical velocities of ice particles due to sedimentation. The studies show a strong dependence of the amount of water that is actively involved in the water cycle on the initial data, model temperatures, and the mechanism of water exchange between the atmosphere and the surface. The general pattern and seasonal asymmetry of the water cycle depends on the size of ice particles, the albedo, and the thermal inertia of the planet’s surface. One of the modeling tasks, which results from a comparison of the model data with those of the TES experiment on board Mars Global Surveyor, is the increase in the total mass of water vapor in the model in the aphelion season and decrease in the mass of water ice clouds at the poles. The surface evaporation scheme, which takes into account the turbulent rise of water vapor, on the one hand, leads to the most complete evaporation of ice from the surface in the summer season in the northern hemisphere and, on the other hand, supersaturates the atmosphere with ice due to the vigorous evaporation, which leads to worse consistency between the amount of the precipitated atmospheric ice and the experimental data. The full evaporation of ice from the surface increases the model sensitivity to the size of the polar cap; therefore, the increase in the latter leads to better results. The use of a more accurate dust scenario changes the model temperatures, which also strongly affects the water cycle.     

Seasonal buffering of atmospheric pressure on Mars

An isothermal reservoir of carbon dioxide in gaseous contact with the Martian atmosphere would reduce the amplitude and advance the phase of global atmospheric pressure fluctuations caused by seasonal growth and decline of polar CO₂ frost caps. Adsorbed carbon dioxide in the upper ～10 m of Martian regolith is sufficient to buffer the present atmosphere on a seasonal basis. Available observations and related polar cap models do not confirm or refute the operation of such a mechanism. Implications for the amplitude and phase of seasonal pressure fluctuations are subject to direct test by the upcoming Viking mission to Mars.