A decadal variance decomposition method is applied to the Northern Hemisphere (NH) 500-hPa geopotential height (GPH) and the sea level pressure (SLP) taken from the last millennium (850–1850 AD) experiment with the coupled climate model CCSM4, to estimate the contribution of the intra-decadal variability to the inter-decadal variability. By removing the intra-decadal variability from the total inter-decadal variability, the residual variability is more likely to be associated with slowly varying external forcings and slow-decadal climate processes, and therefore is referred to as slow-decadal variability. The results show that the (multi-)decadal changes of the NH 500-hPa GPH are primarily dominated by slow-decadal variability, whereas the NH SLP field is primarily dominated by the intra-decadal variability. At both pressure levels, the leading intra-decadal modes each have features related to the El Nino–southern oscillation, the intra-decadal variability of the Pacific decadal oscillation (PDO) and the Arctic oscillation (AO); while the leading slow-decadal modes are associated with external radiative forcing (mostly with volcanic aerosol loadings), the Atlantic multi-decadal oscillation and the slow-decadal variability of AO and PDO. Moreover, the radiative forcing has much weaker effect to the SLP than that to the 500-hPa GPH.