An understanding of the growth and demise of ice sheets over North America is essential to inform future climate models. One poorly studied subject is the glacial dynamics during interstadial Marine Isotope Stage (MIS) 3 (57–29 ka). To better constrain the southern margin of the Laurentide Ice Sheet during this time period, we re-examined a stratigraphic sequence in southeast Manitoba, Canada, and provide robust evidence for advance and retreat of ice. Around 46.6 ± 5.1 ka (1σ error), fluvial sands were deposited under similar precipitation and significantly cooler summer temperatures than present-day. Ice then advanced south over the area, before retreating once again and a return to boreal forest and grassland conditions. The area was then covered by proglacial Lake Vita, dammed by ice to the north. Geochronology constraints indicate Lake Vita existed from ca. 44.3 ± 3.6 to 30.4 ± 2.3 ka (1σ error), although gaps in the optical and finite radiocarbon ages suggest either a lack of data or plausible temporary ice-margin advances during this time period. Ice covering most of Manitoba during MIS 3 is in line with global δ18O records, and glacially influenced sediment deposition in the Mississippi River basin.

Global glaciations have varied in size and magnitude since the Early–Middle Pleistocene transition (~773 ka), despite the apparent regular and high-amplitude 100 ka pacing of glacial–interglacial cycles recorded in marine isotope records. The evidence on land indicates that patterns of glaciation varied dramatically between different glacial–interglacial cycles. For example, Marine Isotope Stages (MIS) 8, 10, and 14 are all noticeably absent from many terrestrial glacial records in North America and Europe. However, globally, the patterns are more complicated, with major glaciations recorded in MIS 8 in Asia and in parts of the Southern Hemisphere, such as Patagonia, for example. This spatial variability in glaciation between glacial–interglacial cycles is likely to be driven by ice volume changes in the West Antarctic Ice Sheet and associated interhemispheric connections through ocean–atmosphere circulatory changes. The weak global glacial imprint in some glacial–interglacial cycles is related to the pattern of global ice buildup. This is caused by feedback mechanisms within glacial systems themselves that partly result from long-term orbital changes driven by eccentricity.

The age of the Great Sand Dunes has been debated for nearly 150 yr. Seven ages ranging from Miocene to late Holocene have been proposed for them. This paper presents new information–chiefly subsurface stratigraphic data, OSL dates, and geomorphic evidence–that indicates that the Great Sand Dunes began to form in the latter part of the middle Pleistocene. The dunes overlie a thick wedge of piedmont-slope deposits, which in turn overlies sediment of Lake Alamosa, a paleolake that began to drain about 440 ka. The wedge of piedmont-slope deposits extends westward for at least 23 km and is as much as 60 m thick at a distance of 10 km from the Sangre de Cristo Range. Ostracodes from one well indicate that the eastern shoreline of Lake Alamosa extended to within 4.3 km of where the Great Sand Dunes eventually formed. The time represented by the wedge of piedmont-slope deposits is not known exactly, but the wedge post-dates 440 ka and was in place prior to 130 ka because by then the dunes overlying it were sufficiently close and tall enough to obstruct streams draining from the Sangre de Cristo Range.