To investigate the mechanisms behind Holocene regional climate trends from north of 58Â°N in the North Atlantic-Fennoscandian region Principal Component Analysis (PCA) was performed and a temperature anomaly stack produced from 81 proxy derived summer temperature time series from 74 sites. The PC results show distinctly different trends for near-surface versus surface temperatures, demonstrating the importance of handling these separately. The first PC of weighted sea surface summer temperature time series and continental time series explains 45Â Â±Â 8% of the variance, where the uncertainty is the standard deviation of the distribution of variance explained across the 1000 age-uncertain ensemble members. PC1 has a relatively uniform expression over the whole region, closely following the summer insolation at 65Â°N. The second PC explains 22Â Â±Â 4% of the variance and shows a non-uniform expression, with loadings in opposite directions in the northern and southeastern parts of the region. Comparing the PC time series with model runs and with the timing of the demise of the Laurentide Ice Sheet (LIS), suggest that this pattern reflects both topographic and albedo effects of the LIS as well as release of meltwater into the North Atlantic and Arctic Oceans. Comparing the stack of gridded records with published global stacks reveals an unusual Holocene temperature development in the North Atlantic-Fennoscandian region most likely resulting from the location relative to the decaying LIS.