The VAMPERS (Vrije Universiteit Amsterdam Permafrost Snow Model) has been coupled within iLOVECLIM, an earth system model. This advancement allows the thermal coupling between permafrost and climate to be examined from a millennial timescale using equilibrium experiments during the Last Glacial Maximum (21 ka) and transient experiments for the subsequent deglaciation period (21-11 ka). It appears that the role of permafrost during both stable and transitional (glacial-interglacial) climate periods is seasonal, resulting in cooler summers and warmer winters by approximately +/- 2 degrees C maximum. This conclusion reinforces the importance of including the active layer within climate models. In addition, the coupling of VAMPERS also yields a simulation of transient permafrost conditions, not only for estimating areal changes in extent but also total permafrost gain/loss.