The VIF sea ice/snow and land snow models developed in a companion paper are implemented in a seasonal dynamic zonally averaged climate model. The simulation of the seasonal cycle of the present climate is compared to climatological data and to results of other models. The addition of ice/snow is found to have interesting consequences for the relative dominance of the opposing winter and summer effects of variations in the Earth's orbital elements, as already implied by Held (1982). In response to a change between extreme high and low summer insolation orbits, it is found that for sea ice (but not for land snow) the winter orbital change dominates over the opposite summer change, causing sea ice shifts which are contrary to the direction predicted by the Milankovitch theory. This is explained by the orbitally induced changes in the length of the seasons, causing insolation anomalies that peak at high latitudes in October and March, coupled with the lag effect of sea ice accretion/melting/heat storage, large thermal capacity of seawater, and sea ice area-thermal inertial feedback. On the other hand, for snow over land the shorter but more pronounced summer orbital change dominates, causing shifts in the region of positive annual snow budgets that agree in direction with the Milankovitch theory. ¿ American Geophysical Union 1988