We use a climate model to simulate two intervals of Permian climate: the Sa
kmarian (ca. 280 Ma), at the end of the major Permo-Carboniferous glaciatio
n, and the Wordian (ca. 265 Ma). We explore the climate sensitivity to vari
ous levels of atmospheric CO2 concentration and to changes in geography and
topography between the two periods. The model simulates large seasonality
and high aridity in the continental interiors of both hemispheres for both
periods. The northern summer monsoon weakens and the southern monsoon stren
gthens between the Sakmarian and the Wordian, owing to changes in geography
and topography. The northern middle and high latitudes cool in winter, bet
ween the Sakmarian and Wordian, associated with northward shift of the cont
inents. This high-latitude cooling strengthens the winter westerlies and sh
ifts the maximum storm-track precipitation south. In the Southern Hemispher
e, the winter westerlies weaken from the Sakmarian to the Wordian. Starting
the simulations with no permanent ice fields (i.e., by assuming that the l
ate Sakmarian postdates deglaciation) and imposing increased levels of atmo
spheric CO2 four times the present level, we find no tendency for reinitiat
ion of major glaciation. Some permanent snow fields do develop in high sout
hern latitudes, but these are primarily at high elevation. However, the com
bination of low CO2 levels (such as present-day levels) and a cold summer o
rbital configuration produces expanded areas of permanent snow. The results
are based on statistics derived from the final 5 yr of 20-yr simulations.
Paleoenvironmental indicators such as coal, evaporite, phosphate, and eolia
n sand deposits agree qualitatively with the simulated climate. The extreme
cold simulated in high latitudes is inconsistent with estimates of high-la
titude conditions. Either the interpretation of observations is incorrect,
the model is incorrect, or both; a possible model deficiency that leads to
cold conditions in high latitudes is the relatively weak ocean-heat transpo
rt simulated by the heat diffusion parameterization of the upper ocean mode
l.