In this analysis, we studied the response of the GCM to a wider range of orbital configurations and water ice reservoirs, and show that during periods of moderate obliquity ( = 25–35°) and high dust opacity (τdust = 1.5–2.5), broad-scale glaciation in the northern mid-latitudes occurs if water ice deposited on the flanks of the Tharsis volcanoes at higher obliquity is available for sublimation. We find that high dust contents of the atmosphere increase its water vapor holding capacity, thereby moving the saturation region to the northern mid-latitudes. Precipitation events are then controlled by topographic forcing of stationary planetary waves and transient weather systems, producing surface ice distribution and amounts that are consistent with the geological record. Ice accumulation rates of 10 mm yr−1 lead to the formation of a 500–1000 m thick regional ice sheet that will produce glacial flow patterns consistent with the geological observations.