We employed the modern analog technique to quantitatively reconstruct temperature and precipitation over the past 2500 yr based on fossil pollen records from six high-elevation lakes in northern Colorado. Reconstructed annual temperatures for the study area did not deviate significantly from modern over the past 2500 yr despite hemispheric expressions of Medieval Climate Anomaly warmth and Little Ice Age cooling. Annual precipitation, however, shifted from lower than modern rates from 2500 to 1000 cal yr BP to higher than modern rates after 1000 cal yr BP, a greater than 100 mm increase in precipitation. Winter precipitation accounts for the majority of the change in annual precipitation, while summer precipitation rates did not change significantly over the past 2500 yr. The large change in winter precipitation rates from the first to second millennium of the Common Era is inferred from a shift in fossil pollen assemblages dominated by subalpine conifers, which have southern sites as modern analogs, to assemblages representing open subalpine vegetation with abundant Artemisia spp. (sagebrush), which have more northern modern analogs. The change helps to explain regional increases in lake levels and shifts in some isotopic and tree-ring data sets, highlighting the risk of large reductions in snowpack and water supplies in the Intermountain West.

A precisely 230Th-dated stalagmite δ13C profile from Hulu Cave, China, is presented to characterize the frequency and pattern of millennial-scale Asian monsoon (AM) variability from 160.6 to 132.5 ka. Evidence for an antiphased relationship of the δ13C and δ18O on the millennial scale suggests that the δ13C is indicative of the local hydrological cycle associated with changes in AM strength. Owing to the δ13C responding to AM changes more sensitively than the δ18O, we could identify 15 strong AM events that correlate to cold intervals recorded in Antarctic ice cores within 230Th dating uncertainty. This result supports a dynamic link of AM strength and southern hemispheric climates via the cross-equatorial airflows. Power spectrum analysis shows a predominant periodicity of 1.5–2.5 ka for the δ13C profile, similar to the Dansgaard-Oeschger frequency during the last glacial period. Moreover, the AM events are characterized by rapid transitions at the onset, suggesting that the observed millennial-scale AM variability is likely forced by northern high-latitude climates via north–south shifts of the Intertropical Convergence Zone associated with the bipolar seesaw mechanism. As evidence for a common mechanism for ice age terminations, a strong AM event (~134 ka) surrounding Termination II is analogous to the Bølling-Allerød warming interval.

A 500,000-yr-long deep-sea sediment-color record from the Northeast Atlantic was investigated to reconstruct the evolution of late Pleistocene climate variability on millennial time scales. Variations of the red–green color intensity are probably caused by climatically induced changes in the ice-rafted input of red-colored iron-bearing terrigenous material to the core site. The resolution of the age model impedes the detection of distinct spectral features at sub-Milankovitch periodicities. Hence, millennial-scale climate variability is quantified as time-dependent variance of the high-pass filtered color time series. The course of the estimated variance shows distinct patterns, which can be linked to continental ice mass. During the past 500,000 yr, large-amplitude millennial-scale climate variability occurs only if continental ice mass exceeds a threshold level, equivalent to sea level at approximately 40% of the lowering during the last glacial maximum.