Mars General Circulation Model Research Articles

Global seasonal variations of the near-surface relative humidity levels on present-day Mars

We investigate the global seasonal variations of near-surface relative humidity and relevant attributes, like temperature and water vapor volume mixing ratio on Mars using calculations from modelled and measurement data. We focus on 2 AM local time snapshots to eliminate daily effects related to differences in insolation, and to be able to compare calculations based on modelling data from the Laboratoire de Météorologie Dynamique Mars General Circulation Model with the observations of Mars Global Surveyor Thermal Emission Spectrometer. We study the seasonal effects by examining four specific dates in the Martian year, the northern spring equinox, summer solstice, autumn equinox, and winter solstice. We identify three specific zones, where the near-surface relative humidity levels are systematically higher than in their vicinity regardless of season. We find that these areas coincide with low thermal inertia features, which control surface temperatures on the planet, and are most likely covered with unconsolidated fine dust with grain sizes smaller than ∼ 40 μm. By comparing the data of relative humidity, temperature and water vapor volume mixing ratio at three different heights (near-surface, ∼ 4 m and ∼ 23 m above the surface), we demonstrate that the thermal inertia could play an important role in determining near-surface humidity levels. We also notice that during the night the water vapor levels drop at ∼ 4 m above the surface. This, together with the temperature and thermal inertia values, shows that water vapor likely condenses in the near-surface atmosphere and on the ground during the night at the three aforementioned regions. This condensation may be in the form of brines, wettening of the fine grains by adsorption or deliquescence. This study specifies areas of interest on the surface of present day Mars for the proposed condensation, which may be examined by in-situ measurements in the future.

Flow Associated With the Condensation and Sublimation of Polar Ice Caps on Mars

AbstractThe flow associated with the sublimation and condensation of carbon dioxide ice in the polar regions of Mars is an important part of the general circulation. However, so far the pattern of the flow is generally not well known. In this study, a Mars general circulation model (MarsWRF) has been applied to derive the pattern of this flow based on some specially designed experiments. The results suggest that the equatorward sublimation and the poleward condensation flows in both hemispheres are rather shallow except near the poles, and the strength of the sublimation flow is significantly stronger than the condensation flow in both hemispheres. The results also show that the strength of the sublimation and condensation flow in the southern hemisphere is much stronger than that in the northern hemisphere. During the summer time, the surface wind in the high‐latitude region near the pole could be dominated by the sublimation flow.

Replication of the historic record of martian global dust storm occurrence in an atmospheric general circulation model

The MarsWRF Mars general circulation model (GCM) with radiatively active dust and orbit-spin coupling reproduces the observational record of Mars years with and without global dust storms (GDS) with a success rate of 77%. Atmospheric conditions diagnostic to the occurrence or non-occurrence of GDS were successfully simulated in 17 of the 22 Mars years of the available historic record. Statistical significance at the 99% level is obtained in a comparison of success rates between an occurrence model with stochastic forcing and those of the GCM forced by orbit-spin coupling. These results provide proof of concept for the orbit-spin coupling hypothesis as a factor contributing to the observed interannual variability of the weather and climate of Mars.

Numerical modeling of orbit-spin coupling accelerations in a Mars general circulation model: Implications for global dust storm activity

We employ the MarsWRF general circulation model (GCM) to test the predictions of a new physical hypothesis: a weak coupling of the orbital and rotational angular momenta of extended bodies is predicted to give rise to cycles of intensification and relaxation of circulatory flows within atmospheres. The dynamical core of MarsWRF has been modified to include the orbit-spin coupling accelerations due to solar system dynamics for the years 1920–2030. The modified GCM is subjected to extensive testing and verification. We compare forced and unforced model outcomes for large-scale zonal and meridional flows, and for near-surface wind velocities and surface wind stresses. The predicted cycles of circulatory intensification and relaxation within the modified GCM are observed. Most remarkably, the modified GCM reproduces conditions favorable for the occurrence of perihelion-season global-scale dust storms (GDSs) on Mars in years in which such storms were observed. A strengthening of the meridional overturning circulation during the dust storm season occurs in the GCM in all recorded years with perihelion-season global-scale dust storms. The increased upwelling produced in the southern hemisphere in southern summer may facilitate the transport of dust to high altitudes in the Mars atmosphere during the dust storm season, where radiative heating may further strengthen the circulation. Significantly increased surface winds and surface wind stresses are also obtained. These may locally facilitate dust lifting from the surface. Based on comparison to the historical record, there is a strong likelihood of a perihelion-season GDS in Mars year 33 and/or Mars year 34.

ASIS v1.0: an adaptive solver for the simulation of atmospheric chemistry

Abstract. This article reports on the development and tests of the adaptive semi-implicit scheme (ASIS) solver for the simulation of atmospheric chemistry. To solve the ordinary differential equation systems associated with the time evolution of the species concentrations, ASIS adopts a one-step linearized implicit scheme with specific treatments of the Jacobian of the chemical fluxes. It conserves mass and has a time-stepping module to control the accuracy of the numerical solution. In idealized box-model simulations, ASIS gives results similar to the higher-order implicit schemes derived from the Rosenbrock's and Gear's methods and requires less computation and run time at the moderate precision required for atmospheric applications. When implemented in the MOCAGE chemical transport model and the Laboratoire de Météorologie Dynamique Mars general circulation model, the ASIS solver performs well and reveals weaknesses and limitations of the original semi-implicit solvers used by these two models. ASIS can be easily adapted to various chemical schemes and further developments are foreseen to increase its computational efficiency, and to include the computation of the concentrations of the species in aqueous-phase in addition to gas-phase chemistry.

The variability, structure and energy conversion of the northern hemisphere traveling waves simulated in a Mars general circulation model

Investigations of the variability, structure and energetics of the m=1−3 traveling waves in the northern hemisphere of Mars are conducted with the MarsWRF general circulation model. Using a simple, annually repeatable dust scenario, the model reproduces many general characteristics of the observed traveling waves. The simulated m=1 and m=3 traveling waves show large differences in terms of their structures and energetics. For each representative wave mode, the geopotential signature maximizes at a higher altitude than the temperature signature, and the wave energetics suggests a mixed baroclinic-barotropic nature. There is a large contrast in wave energetics between the near-surface and higher altitudes, as well as between the lower latitudes and higher latitudes at high altitudes. Both barotropic and baroclinic conversions can act as either sources or sinks of eddy kinetic energy. Band-pass filtered transient eddies exhibit strong zonal variations in eddy kinetic energy and various energy transfer terms. Transient eddies are mainly interacting with the time mean flow. However, there appear to be non-negligible wave-wave interactions associated with wave mode transitions. These interactions include those between traveling waves and thermal tides and those among traveling waves.

New observations of molecular nitrogen in the Martian upper atmosphere by IUVS on MAVEN

AbstractWe identify molecular nitrogen (N2) emissions in the Martian upper atmosphere using the Imaging Ultraviolet Spectrograph (IUVS) on NASA's Mars Atmosphere and Volatile EvolutioN (MAVEN) mission. We report the first observations of the N2 Lyman‐Birge‐Hopfield (LBH) bands at Mars and confirm the tentative identification of the N2 Vegard‐Kaplan (VK) bands. We retrieve N2 density profiles from the VK limb emissions and compare calculated limb radiances between 90 and 210 km against both observations and predictions from a Mars general circulation model (GCM). Contrary to earlier analyses using other satellite data, we find that N2 abundances exceed GCM results by about a factor of 2 at 130 km but are in agreement at 150 km. The analysis and interpretation are enabled by a linear regression method used to extract components of UV spectra from IUVS limb observations.

Modeling the effects of martian surface frost on ice table depth

Ground ice has been observed in small fresh craters in the vicinity of the Viking 2 lander site (48°N, 134°E). To explain these observations, current models for ground ice invoke levels of atmospheric water of 20 precipitable micrometers – higher than observations. However, surface frost has been observed at the Viking 2 site and surface water frost and snow have been shown to have a stabilizing effect on Antarctic subsurface ice. A snow or frost cover provides a source of humidity that should reduce the water vapor gradient and hence retard the sublimation loss from subsurface ice. We have modeled this effect for the Viking 2 landing site with combined ground ice and surface frost models. Our model is driven by atmospheric output fields from the NASA Ames Mars General Circulation Model (MGCM). Our modeling results show that the inclusion of a thin seasonal frost layer, present for a duration similar to that observed by the Viking Lander 2, produces ice table depths that are significantly shallower than a model that omits surface frost. When a maximum frost albedo of 0.35 was permitted, seasonal frost is present in our model from Ls=182° to Ls=16°, resulting in an ice table depth of 64cm – which is 24cm shallower than the frost-free scenario. The computed ice table depth is only slightly sensitive to the assumed maximum frost albedo or thickness in the model.

Initial results from radio occultation measurements with the Mars Reconnaissance Orbiter: A nocturnal mixed layer in the tropics and comparisons with polar profiles from the Mars Climate Sounder

The Mars Reconnaissance Orbiter (MRO) performs radio occultation (RO) measurements on selected orbits, generally once per day. We have retrieved atmospheric profiles from two subsets of data, yielding a variety of new results that illustrate the scientific value of the observations. One set of measurements sounded the tropics in northern summer at a local time ∼1h before sunrise. Some of these profiles contain an unexpected layer of neutral stability with a depth of ∼4km and a pressure at its upper boundary of ∼160Pa. The mixed layer is bounded above by a temperature inversion and below by another strong inversion adjacent to the surface. This type of structure is observed near Gale Crater, in the Tharsis region, and at a few other locations, whereas profiles in Amazonis Planitia and Elysium Planitia show no sign of a detached mixed layer with an overlying inversion. We supplemented the measurements with numerical simulations by the NASA Ames Mars General Circulation Model, which demonstrate that water ice clouds can generate this distinctive type of temperature structure through their influence on radiative transfer at infrared wavelengths. In particular, the simulations predict the presence of a nocturnal cloud layer in the Tharsis region at a pressure of ∼150Pa (∼10km above the surface), and the nighttime radiative cooling at cloud level is sufficient to produce a temperature inversion above the cloud as well as convective instability below the cloud, consistent with the observations. The second set of measurements sounded mid-to-high northern latitudes in spring, when carefully coordinated observations by the MRO Mars Climate Sounder (MCS) are also available. The differences between the RO and MCS temperature profiles are generally consistent with the expected performance of the two instruments. Within this set of 21 comparisons the average temperature difference is less than 1K where the aerosol opacities are smaller than 10-3km-1, at pressures of 10–50Pa, whereas it increases to ∼2K where the aerosol opacities exceed this threshold, at pressures of 50–300Pa. The standard deviation of the temperature difference is ∼2K, independent of pressure. The second set of RO measurements also provides unique information about the stability of the annual CO2 cycle and the dynamics near the edge of the seasonal CO2 ice cap.

Local Dynamics of Baroclinic Waves in the Martian Atmosphere

Abstract The paper investigates the processes that drive the spatiotemporal evolution of baroclinic transient waves in the Martian atmosphere by a simulation experiment with the Geophysical Fluid Dynamics Laboratory (GFDL) Mars general circulation model (GCM). The main diagnostic tool of the study is the (local) eddy kinetic energy equation. Results are shown for a prewinter season of the Northern Hemisphere, in which a deep baroclinic wave of zonal wavenumber 2 circles the planet at an eastward phase speed of about 70° Sol−1 (Sol is a Martian day). The regular structure of the wave gives the impression that the classical models of baroclinic instability, which describe the underlying process by a temporally unstable global wave (e.g., Eady model and Charney model), may have a direct relevance for the description of the Martian baroclinic waves. The results of the diagnostic calculations show, however, that while the Martian waves remain zonally global features at all times, there are large spatiotemporal changes in their amplitude. The most intense episodes of baroclinic energy conversion, which take place in the two great plain regions (Acidalia Planitia and Utopia Planitia), are strongly localized in both space and time. In addition, similar to the situation for terrestrial baroclinic waves, geopotential flux convergence plays an important role in the dynamics of the downstream-propagating unstable waves.

Three‐dimensional Martian ionosphere model: I. The photochemical ionosphere below 180 km

We describe the Mars ionosphere with unprecedented detail in 3‐D, as simulated by a Mars general circulation model (the Laboratoire de Météorologie Dynamique Mars GCM), and compare it with recent measurements. The model includes a number of recent extensions and improvements. Different simulations for a full Martian year have been performed. The electron density at the main ionospheric peak is shown to vary with the Sun‐Mars distance and with the solar variability, both in the long‐term (11 year solar cycle) and on shorter temporal scales (solar rotation). The main electronic peak is shown to be located at the same pressure level during all the Martian year. As a consequence, its altitude varies with latitude, local time, and season according to the natural expansions and fluctuations of the neutral atmosphere, in agreement with previous models. The model predicts a nighttime ionosphere due only to photochemistry. The simulated ionosphere close to the evening terminator is in agreement with observations. No effort has been made to explain the patchy ionosphere observed in the deep nightside. We have compared the modeled ionosphere with Mars Global Surveyor and Mars Advanced Radar for Subsurface and Ionosphere Sounding data. The model reproduces the solar zenith angle variability of the electron density and the altitude of the peak, although it underestimates the electron density at the main peak by about 20%. The electron density at the secondary peak is strongly underestimated by the model, probably due to a very crude representation of the X‐ray solar flux. This is one of the aspects that needs a revision in future versions of the model.

New nitric oxide (NO) nightglow measurements with SPICAM/MEx as a tracer of Mars upper atmosphere circulation and comparison with LMD‐MGCM model prediction: Evidence for asymmetric hemispheres

We report observations of NO nightglow with the Spectroscopy for the Investigation of the Characteristics of the Atmosphere of Mars (SPICAM) experiment on board the Mars Express (MEx) spacecraft. NO molecules emit an ultraviolet photon when N and O atoms (produced at high altitude in the thermosphere) recombine. Therefore, this emission is a tracer of the atmospheric dynamics in the lower thermosphere where O and N atoms are produced, and below, in the altitude region 50–100 km where the emission is detected. A new retrieval method has been developed to analyze the measurements from this instrument in the stellar occultation mode without slit and retrieve the absolute brightness of the emission. We present the results from the processing of more than 2000 orbits, providing the first global latitude‐season distribution of the emission, established over three Martian years. The results are globally consistent with previously available measurements of dedicated limb nightglow obtained during the first Martian year of MEx (MY27). We compared the ensemble of both data sets with the predictions of the Laboratoire de Météorologie Dynamique Mars General Circulation Model (LMD‐MGCM), with the addition of the full chemistry of N atoms. We find an overall agreement between the observed and modeled airglow intensities, but discrepancies are also found. The frequency and magnitude of the NO airglow observations show important asymmetries between the Northern and the Southern Hemispheres. There is no detection of emission near the poles during equinox conditions, while the model predicts that it should be most intense because of a circulation with two descending branches at the poles.

The early Martian atmosphere: Investigating the role of the dust cycle in the possible maintenance of two stable climate states

AbstractC. Leovy (personal communication, 2007) speculated that two stable climate states on early Mars could have resulted from interactions between the dust and CO2 cycles. In one state, a highly active dust cycle would prevent atmospheric collapse, and in the second, the collapsed atmosphere would not maintain an active dust cycle. An initial assessment of this idea is presented based on a Mars general circulation model parameter study. A range of global dust loadings, CO2 ice albedos, and obliquities are investigated to explore conditions in which increasing the atmospheric dust content stabilizes an otherwise unstable atmosphere. We find that dust only stabilizes the atmosphere at high obliquity and when the CO2 ice albedo is high. Although results suggest that two stable states could have existed on early Mars under limited conditions, further work is needed to know if the conditions necessary are physically plausible.

The semidiurnal tide in the middle atmosphere of Mars

AbstractAtmospheric thermal tides are global oscillations in atmospheric fields that are subharmonics of a solar day. While atmospheric tides on Earth are mainly relevant in the upper atmosphere, on Mars, they dominate temperature variations and winds throughout the atmosphere. Observations and model simulations to date have suggested that the migrating diurnal tide is the predominant mode in the Martian atmosphere, and that the semidiurnal tide is only relevant in the tropical middle atmosphere during conditions of high dust loading. New comprehensive observations by the Mars Climate Sounder in a geometry that allows coverage of multiple local times show that the semidiurnal tide is a dominant response of the Martian atmosphere throughout the Martian year. The maximum semidiurnal amplitude of ~ 16 K is found at southern winter high latitudes, which makes it the largest tidal amplitude observed in the Martian middle atmosphere outside of dust storm conditions. The semidiurnal tide can be successfully modeled due to recent advances of Mars General Circulation Models (MGCMs) that include the radiatively active treatment of water ice clouds. Tidal forcing occurs through absorption of radiation by aerosols and points to the vertical structure of dust and clouds and their radiative effects as being essential for our understanding of the thermal structure and the general circulation of the Martian atmosphere. As with terrestrial GCMs trying to quantify mechanisms affecting climate, future Mars modeling efforts will require microphysical schemes to control aerosol distributions, and vertically and temporally resolved measurements of temperature and aerosols will be essential for their validation.

Simulating the interannual variability of major dust storms on Mars using variable lifting thresholds

The redistribution of a finite amount of martian surface dust during global dust storms and in the intervening periods has been modelled in a dust lifting version of the UK Mars General Circulation Model. When using a constant, uniform threshold in the model’s wind stress lifting parameterisation and assuming an unlimited supply of surface dust, multiannual simulations displayed some variability in dust lifting activity from year to year, arising from internal variability manifested in surface wind stress, but dust storms were limited in size and formed within a relatively short seasonal window. Lifting thresholds were then allowed to vary at each model gridpoint, dependent on the rates of emission or deposition of dust. This enhanced interannual variability in dust storm magnitude and timing, such that model storms covered most of the observed ranges in size and initiation date within a single multiannual simulation. Peak storm magnitude in a given year was primarily determined by the availability of surface dust at a number of key sites in the southern hemisphere. The observed global dust storm (GDS) frequency of roughly one in every 3years was approximately reproduced, but the model failed to generate these GDSs spontaneously in the southern hemisphere, where they have typically been observed to initiate. After several years of simulation, the surface threshold field—a proxy for net change in surface dust density—showed good qualitative agreement with the observed pattern of martian surface dust cover. The model produced a net northward cross-equatorial dust mass flux, which necessitated the addition of an artificial threshold decrease rate in order to allow the continued generation of dust storms over the course of a multiannual simulation. At standard model resolution, for the southward mass flux due to cross-equatorial flushing storms to offset the northward flux due to GDSs on a timescale of ∼3years would require an increase in the former by a factor of 3–4. Results at higher model resolution and uncertainties in dust vertical profiles mean that quasi-periodic redistribution of dust on such a timescale nevertheless appears to be a plausible explanation for the observed GDS frequency.

Development of a fast, accurate radiative transfer model for the Martian atmosphere, past and present

We present details of an approach to creating ak‐distribution radiative transfer model (KDM) for use in the Martian atmosphere. Such models preserve the accuracy of more rigorous line‐by‐line models, but are orders of magnitude faster, and can be effectively implemented in 3‐D general circulation models. The approach taken here is sufficiently generalized that it can be employed for atmospheres of any arbitrary composition and mass, and demonstrations are provided for simulated atmospheres with a present‐day Martian surface pressure (∼6 mb) and a putative thick early Mars atmosphere (∼500 mb), both with and without atmospheric water vapor. KDM‐derived absorption coefficients are placed into a look‐up table at a set of gridded points in pressure, temperature and atmospheric composition, and a tri‐linear interpolation scheme is used to obtain the coefficients appropriate for the local atmospheric conditions. These coefficients may then be used within any of a variety of commonly used flux solvers to obtain atmospheric heating rates. A series of validation tests are performed with the KDM for both present‐day and early Mars atmospheric conditions, and the model is compared against several other widely used radiative transfer schemes, including several used in contemporary general circulation models. These validation results identify weaknesses in some other approaches and demonstrate the efficacy of the KDM, providing a rigorous test of these types of models for use in the Martian atmosphere. A demonstration of results obtained by implementing the KDM in a Mars general circulation model is provided.

Identifying Martian atmospheric instabilities and their physical origins using bred vectors

AbstractThe bred vector (BV) technique applied to the Geophysical Fluid Dynamics Laboratory (GFDL) Mars General Circulation Model (MGCM) identifies regions and seasons of instability of the Martian atmosphere, and a kinetic energy equation applied to the control and perturbed states elucidates their physical origins. Instabilities prominent in the late autumn to early spring seasons of each hemisphere along the polar temperature front result from baroclinic conversions from BV potential to BV kinetic energy near the surface, whereas both baroclinic and barotropic conversions play a role for the westerly jets aloft. The low‐level tropics and the northern hemisphere summer are relatively stable, with negative bred vector growth rates. Bred vector growth precedes initiation of travelling wave activity in the midlatitudes during the transition seasons, and their structure relates to the eddy field. Topography plays a role in determining favoured locations for near‐surface instabilities. Bred vectors are also linked to forecast ensemble spread in data assimilation and help explain the growth of forecast errors. We finally note that the ability to use breeding to identify instabilities as well as their physical origin depends on the fact that both the control and the perturbed solutions that give rise to bred vectors satisfy exactly the model's governing equations. As a result, this approach can be used with any dynamical system represented by a model. Copyright © 2012 Royal Meteorological Society

An aerodynamic roughness length map derived from extended Martian rock abundance data

Many boundary layer processes simulated within a Mars General Circulation Model (MGCM), including the description of the processes controlling dust rising from the Martian surface, are highly sensitive to the aerodynamic roughness length z0. On the basis of rock‐size frequency distributions inferred from different Martian landing sites and Earth analog sites, we have first established that lognormal‐modeled rock‐size frequency distributions are able to reproduce correctly the observed Martian rock populations. We have validated the hypothesis that the rock abundance ζ of a given area could be estimated at a first order from its thermophysical properties, namely its thermal inertia I and its albedo α. We have demonstrated the possibility of using rock abundance ζ to estimate the roughness density λ on Mars and to retrieve subsequently the aerodynamic roughness length by using semi‐empirical relationships based on terrestrial wind‐tunnel and field measurements. By combining our methodology with remote sensing measurements of the Thermal Emission Spectrometer aboard Mars Global Surveyor, we have derived a global map of the aeolian aerodynamic roughness length with a 1/8° × 1/8° resolution over the entire Martian surface. Contrary to what is often assumed, the Martian aeolian aerodynamic roughness length is spatially highly heterogeneous. At the fullest resolution, the Martian aerodynamic roughness length varies from 10−3 cm to 2.33 cm. About 84% of the Martian surface seems to be characterized by an aeolian aerodynamic roughness length value lower than 1 cm, the spatially uniform value that most of the MGCMs simulations have assumed recently. Since the aerodynamic roughness length z0 is a key parameter in deriving the erosion threshold wind velocities, we anticipate a significant impact of our findings on the efficiencies for lifting dust in future MGCMs.

The Ashima/MIT Mars GCM and argon in the martian atmosphere

We investigate the ability of modern general circulation models (GCMs) to simulate transport in the martian atmosphere using measurements of argon as a proxy for the transport processes. Argon provides the simplest measure of transport as it is a noble gas with no sinks or sources on seasonal timescales. Variations in argon result solely from ‘freeze distillation’, as the atmosphere condenses at the winter poles, and from atmospheric transport. Comparison of all previously published models when rescaled to a common definition of the argon enhancement factor (EF) suggest that models generally do a poor job in predicting the peak enhancement in southern winter over the winter pole – the time when the capability of the model transport approaches are most severely tested. Despite observed peak EF values of ∼6, previously published model predictions peaked at EF values of only 2–3. We introduce a new GCM that provides a better treatment of mass conservation within the dynamical core, includes more sophisticated tracer transport approaches, and utilizes a cube–sphere grid structure thus avoiding the grid-point convergence problem at the pole that exists for most current Mars GCMs. We describe this model – the Ashima Research/Massachusetts Institute of Technology Mars General Circulation Model (Ashima/MIT Mars GCM) and use it to demonstrate the significant sensitivity of peak EF to the choices of transport approach for both tracers and heat. We obtain a peak EF of 4.75 which, while over 50% higher than any prior model, remains well short of the observed value. We show that the polar EF value in winter is primarily determined by the competition between two processes: (1) mean meridional import of lower-latitude air not enriched in argon and (2) the leakage of enriched argon out of the polar column by eddies in the lowest atmospheric levels. We suggest possibilities for improving GCM representation of the CO2 cycle and the general circulation that may further improve the simulation of the argon cycle. We conclude that current GCMs may be insufficient for detailed simulation of transport-sensitive problems like the water cycle and potentially also the dust cycle.

Modeling the hydrological cycle on Mars

The study provides a detailed analysis of the hydrological cycle on Mars simulated with a newly developed microphysical model, incorporated in a spectral Mars General Circulation Model. The modeled hydrological cycle is compared well with simulations of other global climate models. The simulated seasonal migration of water vapor, circulation instability, and the high degree of temporal variability of localized water vapor outbursts are shown closely consistent with recent observations. The microphysical parameterization provides a significant improvement in the modeling of ice clouds evolved over the tropics and major ancient volcanoes on Mars. The most significant difference between the simulations presented here and other GCM results is the level at which the water ice clouds are found. The model findings also support interpretation of observed thermal anomalies in the Martian tropics during northern spring and summer seasons.