Understanding the sequence of events occurring during the last major glacial to interglacial transition (21 ka bp to 9 ka bp) is a challenging task that has the potential to unveil the mechanisms behind large scale climate changes. Though many studies have focused at a complex understanding of the sequence of rapid climatic change that accompanied or interrupted the deglaciation, few have analysed it in a more theoretical framework with simple forcings. In the following, we address when and where the first significant temperature anomalies appear when using slow varying forcing of the last deglaciation. We use here coupled transient simulations of the last deglaciation, including ocean, atmosphere and vegetation components to analyse the spatial timing of the deglaciation. To keep the analysis in a simple framework, we do not include rapid freshwater forcings that have led to rapid climate shifts during that time period. We aim to disentangle the direct and subsequent response of the climate system to slow forcing and moreover the location where those changes are more clearly expressed. In a data-modelling comparison perspective, this could help understanding the physically plausible phasing between known forcings and recorded climatic changes. Our analysis of climate variability could also help to distinguish deglacial warming signals from internal climate variability. We thus are able to better pinpoint the onset of local deglaciation, as defined by the first significant local warming, and further show that there is a large regional variability associated with it, even with the set of slow forcings used here.