Obliquity Driven Climate Change in Mars' Recent Past

We use the NASA/Ames grid Mars general circulation model with an updated radiation code and cloud nlicrophysics scheme. To speed up the simulations, we run the model at fairly coarse resolution (7.5" latitude x 22.5" longitude). Radiation Code Fluxes and heating rates are calculated from a radiation code based on the two-stream solution to radiative transfer that fully accounts for multiple scattering in the presence of gaseous absorp- tion. The model has 12 spectral intervals. Dust and water ice scattering properties are included. For dust, we use the Ockert-Bell [9] values in the visible, and Forget [IO] values in the infrared. For ice, we can either compute them online as the cloud evolves, or we can specify them. Gaseous opacities for water vapor and CO2 are calculated from correlated k-distributions taken from full line-by-line models. Cloud Microphysics Our cloud scheme is based on momentlorder scheme in which the mass mixing ratio and number density of the cloud ensemble are the advected species. From these we obtain a mean particle size and an estimate or the particle size distribution assuming a variance) which we then divide into 8 bins. Cloud microphysics is performed in each of these bins and includes nucleation, condensation, and gravitational settling. Dust is treated as a tracer and serves as condensation nuclei. The altered size distribution is then converted back into a mean size, a mixing ratio, and a particle number density.