The temporal and spatial coarseness of megafaunal fossil records complicates attempts to to disentangle the relative impacts of climate change, ecosystem restructuring, and human activities associated with the Late Quaternary extinctions. Advances in the extraction and identification of ancient DNA that was shed into the environment and preserved for millennia in sediment now provides a way to augment discontinuous palaeontological assemblages. Here, we present a 30,000-year sedimentary ancient DNA (sedaDNA) record derived from loessal permafrost silts in the Klondike region of Yukon, Canada. We observe a substantial turnover in ecosystem composition between 13,500 and 10,000 calendar years ago with the rise of woody shrubs and the disappearance of the mammoth-steppe (steppe-tundra) ecosystem. We also identify a lingering signal of Equus sp. (North American horse) and Mammuthus primigenius (woolly mammoth) at multiple sites persisting thousands of years after their supposed extinction from the fossil record.

The collapse of the steppe-tundra biome (mammoth steppe) at the end of the Pleistocene is used as an important example of top-down ecosystem cascades, where human hunting of keystone species led to profound changes in vegetation across high latitudes in the Northern Hemisphere. Alternatively, it is argued that this biome transformation occurred through a bottom-up process, where climate-driven expansion of shrub tundra (Betula, Salix spp.) replaced the steppe-tundra vegetation that grazing megafauna taxa relied on. In eastern Beringia, these differing hypotheses remain largely untested, in part because the precise timing and spatial pattern of Late Pleistocene shrub expansion remains poorly resolved. This uncertainty is caused by chronological ambiguity in many lake sediment records, which typically rely on radiocarbon (14C) dates from bulk sediment or aquatic macrofossils—materials that are known to overestimate the age of sediment layers. Here, we reexamine Late Pleistocene pollen records for which 14C dating of terrestrial macrofossils is available and augment these data with 14C dates from arctic ground-squirrel middens and plant macrofossils. Comparing these paleovegetation data with a database of published 14C dates from megafauna remains, we find the postglacial expansion of shrub tundra preceded the regional extinctions of horse (Equus spp.) and mammoth (Mammuthus primigenius) and began during a period when the frequency of 14C dates indicates large grazers were abundant. These results are not consistent with a model of top-down ecosystem cascades and support the hypothesis that climate-driven habitat loss preceded and contributed to turnover in mammal communities.

The timing of Laurentide Ice Sheet deglaciation along its southwestern margin controlled the evolution of large glacial lakes and has implications for human migration into the Americas. Accurate reconstruction of the ice sheet’s retreat also constrains glacial isostatic adjustment models and is important for understanding ice-sheet sensitivity to climate forcing. Despite its significance, retreat of the southwestern Laurentide Ice Sheet (SWLIS) is poorly constrained by minimum-limiting 14C data. We present 26 new cosmogenic 10Be exposure ages spanning the western Interior Plains, Canada. Using a Bayesian framework, we combine these data with geomorphic mapping, 10Be, and high-quality minimum-limiting 14C ages to provide an updated chronology. This dataset presents an internally consistent retreat record and indicates that the initial detachment of the SWLIS from its convergence with the Cordilleran Ice Sheet began by ca. 15.0 ka, concurrent with or slightly prior to the onset of the Bølling-Allerød interval (14.7–12.9 ka) and retreated >1200 km to its Younger Dryas (YD) position in ~2500 yr. Ice-sheet stabilization at the Cree Lake Moraine facilitated a meltwater drainage route to the Arctic from glacial Lake Agassiz within the YD, but not necessarily at the beginning. Our record of deglaciation and new YD constraints demonstrate deglaciation of the Interior Plains was ~60% faster than suggested by minimum 14C constraints alone. Numerical modeling of this rapid retreat estimates a loss of ~3.7 m of sea-level equivalent from the SWLIS during the Bølling-Allerød interval.