Search

Because global sea level during the last interglacial (LIG; 130–115 ka) was higher than today, the LIG is a useful approximate analogue for improving predictions of future sea-level rise. Here, we synthesize sea-level proxies for the LIG in the glaciated Northern Hemisphere for inclusion in the World Atlas of Last Interglacial Shorelines (WALIS) database. We describe 82 sites from Russia, northern Europe, Greenland and North America from a variety of settings, including boreholes, riverbank exposures and along coastal cliffs. Marine sediments at these sites were constrained to the LIG using a variety of radiometric methods (radiocarbon, uranium–thorium, potassium–argon), non-radiometric methods (amino acid dating, luminescence methods, electron spin resonance, tephrochronology) as well as various stratigraphic and palaeo-environmental approaches. In general, the sites reported in this paper do not offer constraint on the global LIG highstand, but rather evidence of glacial isostatic adjustment (GIA)-influenced sea-level positions following the Marine Isotope Stage 6 glaciation (MIS 6; 191–130 ka). Most of the proxies suggest that sea level was much higher during the LIG than at the present time. Moreover, many of the sites show evidence of regression due to sea-level fall (owing to glacial isostatic uplift), and some also show fluctuations that may reflect regrowth of continental ice or increased influence of the global sea-level signal. In addition to documenting LIG sea-level sites in a large swath of the Northern Hemisphere, this compilation is highly relevant for reconstructing the size of MIS 6 ice sheets through GIA modelling. The database is available at https://doi.org/10.5281/zenodo.5602212 (Dalton et al., 2021).

The North American Ice Sheet Complex (comprising the Laurentide, Cordilleran and Innuitian ice sheets) was the largest ice mass in the Northern Hemisphere that grew towards and waned after the Last Glacial Maximum. The existing ice margin chronology available for the North American Ice Sheet Complex is based on radiocarbon data only and does not reflect other geochronometric information constraining the last deglaciation, such as cosmogenic exposure- or optically stimulated luminescence ages. Here we present a series of newly produced ice margin isochrones from 25 ka to present, in a time step of 500 years. For each isochron, we draw maximum, best estimate, and minimum ice margin position in an attempt to capture the existing uncertainty. The ice margin isochrones are based on (i) an up-to-date dataset of radiocarbon ages (~5000), (ii) 10Be and 26Al cosmogenic nuclide data that directly date ~80 ice-marginal features over North America, (iii) ~350 optically stimulated luminescence ages dating the deposition of an aeolian cover immediately post-deglaciation, (iv) the ice-sheet scale glacial geomorphology record. Our effort brings the information on the last North American Ice Sheet Complex deglaciation on par with that for the Eurasian Ice Sheets and should serve the broad community of Quaternary research from archaeology to numerical ice sheet modelling.

Because global sea level during the last interglacial (LIG; 130–115 ka) was higher than today, the LIG is a useful analogue for improving predictions of future sea level rise. Here, we synthesize sea level proxies for the LIG in the glaciated Northern Hemisphere for inclusion in the World Atlas of Last Interglacial Shorelines (WALIS) database. We describe 82 sites from Russia, northern Europe, Greenland and North America from a variety of settings, including boreholes, riverbank exposures and along coastal cliffs. Marine sediments at these sites were constrained to the LIG using a variety of radiometric methods (radiocarbon, U-Series dating, K-Ar dating), non-radiometric methods (amino acid dating, luminescence methods, and electron spin resonance, tephrochronology) as well as various stratigraphic and palaeoenvironmental approaches. As the areas in this database were covered by ice sheets from the penultimate glaciation and were affected by glacial isostatic adjustment (GIA), most of the proxies show that sea level was much higher than present during the LIG. Many of the sites show evidence of regression due sea level fall due to GIA uplift, and some also show fluctuations that may reflect regrowth of continental ice or increased influence of the global sea level signal. The database is available at https://doi.org/10.5281/zenodo.5602212 (Dalton et al., 2021).

We report previously unpublished evidence for a Marine Isotope Stage 3 (MIS 3; 60–26 ka) glacial outburst flood in the Torngat Mountains (northern Quebec/Labrador, Canada). We present 10Be cosmogenic exposure ages from legacy fieldwork for a glacial lake shoreline with evidence for outburst flooding in the Torngat Mountains, with a minimum age of 36 ± 3 ka (we consider the most likely age, corrected for burial, to be ~56 ± 3 ka). This shoreline position and age can potentially constrain the Laurentide Ice Sheet margin in the Torngat Mountains. This region, considered a site of glacial inception, has no published dated geologic constraints for high-elevation MIS 3 ice margins. We estimate the freshwater flux associated with the inferred glacial outburst flood using high-resolution digital elevation maps corrected for glacial isostatic adjustment. Using assumptions about the ice-dammed locations we find that a freshwater flood volume of 1.14 × 1012 m3 could have entered the Hudson Strait. This glacial outburst flood volume could have contributed to surface ocean freshening to cause a measurable meltwater signal in δ18O records, but would not necessarily have been associated with substantial ice rafted debris. Future work is required to refine estimates of the size and timing of such a glacial outburst flood. Nevertheless, we outline testable hypotheses about the Laurentide Ice Sheet and glacial outburst floods, including possible implications for Heinrich events and glacial inception in North America, that can be assessed with additional fieldwork and cosmogenic measurements.

Accurately reconstructing the paleogeography of the Laurentide Ice Sheet (LIS) during Marine Isotope Stage 3 (MIS 3; ca. 57,000 to ca. 29,000 yr B.P.) is critical for understanding glacial growth toward the Last Glacial Maximum (LGM), refining sea-level histories, and studying the Earth system response to rapid climate change events. Here, we present a geochronological data set useful for testing hypotheses of global sea level and refining ice sheet configuration through this interval. Data (n = 735) span the entire MIS 3 interval and consist of 14C determinations (n = 651), cosmogenic exposure ages (n = 52), and optically stimulated luminescence dates (n = 32). On that basis, we hypothesize that the central region of the LIS underwent a dramatic reduction in ice from ca. 52 to 40 ka. Key to this hypothesis are geological records at sites in the Hudson Bay Lowlands (east central Canada) that suggest a marine incursion and development of terrestrial landscapes. We show that these landscapes are consistent with recently published glacial isostatic adjustment predictions that include widespread deglaciation of the eastern (Labrador) sector of the LIS with ice buildup over the western (Keewatin) sector at 42 ka. Ice growth from this minimum toward the LGM is likely to have been rapid. The agreement between this data set and modeling predictions prompts the reassessment of key Late Pleistocene records, including Heinrich events, loess deposition in the continental United States, and sedimentological records from the Gulf of Mexico.

Stratigraphic records from formerly glaciated regions are critical for detailed study of the timing, onset and dynamics of past ice sheets and the palaeoecology of previous ice-free intervals. We examined three stratigraphic sections from an 18-km stretch of the Albany River, Hudson Bay Lowlands, Canada, located at the geographic center for many Late Pleistocene ice sheets. Till characterization and correlation suggest that at least three glacial advances from shifting ice centers within the Labrador sector of the Laurentide Ice Sheet were preserved in these stratigraphic records. Non-glacial units (fluvial, organic-bearing sediments) were constrained via optically stimulated luminescence to two possible periods at ca. 73,000 to 68,000 yr BP and ca. 60,000 yr BP. Boreal and peatland taxa (Picea, Pinus, Poaceae, Betula, Cyperaceae, Sphagnum) dominated the pollen record at each site, whereas plant macrofossils analyzed at one site confirm the local presence of conifer trees (bark, needles, seed wings), bryophytes (largely Scotpidium spp), herbaceous plants (Caryophyllaceae, Carex, Poaceae), and an aquatic setting (e.g. Potamogeton, ephippia of Daphnia spp). Pollen-derived average summer temperature reconstructions suggested that local temperatures at the Albany sites were between 12 and 15 degrees C, which is similar to present-day estimates for the region (14.2 degrees C). Reconstructed annual precipitation estimates were 580-640 mm, which is similar to slightly higher than present-day estimates (564 mm). Non-glacial intervals at the Albany sites likely represent abandoned fluvial environments that supported water-logged peatland biota. Results from this research contribute toward ongoing efforts to constrain ice sheet dynamics over North America during the last glacial cycle (e.g. 71,000-14,000 yr BP) and provide insight into the complex Late Pleistocene palaeoclimate record at the innermost area of the glaciated region.

The Laurentide Ice Sheet was the largest global ice mass to grow and decay during the last glacial cycle (~115 ka to ~10 ka). Despite its importance for driving major changes in global mean sea level, long-term landscape evolution, and atmospheric circulation patterns, the history of the Laurentide (and neighbouring Innuitian) Ice Sheet is poorly constrained owing to sporadic preservation of stratigraphic records prior to the Last Glacial Maximum (LGM; ~25 ka) and a case-study approach to the dating of available evidence. Here, we synthesize available geochronological data from the glaciated region, together with published stratigraphic and geomorphological data, as well as numerical modelling output, to derive 19 hypothesised reconstructions of the Laurentide and Innuitian ice sheets from 115 ka to 25 ka at 5-kyr intervals, with uncertainties quantified to include best, minimum, and maximum ice extent estimates at each time-step. Our work suggests that, between 115 ka and 25 ka, some areas of North America experienced multiple cycles of rapid ice sheet growth and decay, while others remained largely ice-free, and others were continuously glaciated. Key findings include: (i) the growth and recession of the Laurentide Ice Sheet from 115 ka through 80 ka; (ii) significant build-up of ice to almost LGM extent at ~60 ka; (iii) a potentially dramatic reduction in North American ice at ~45 ka; (iv) a rapid expansion of the Labrador Dome at ~38 ka; and (v) gradual growth toward the LGM starting at ~35 ka. Some reconstructions are only loosely constrained and are therefore speculative (especially prior to 45 ka). Nevertheless, this work represents our most up-to-date understanding of the build-up of the Laurentide and Innuitian ice sheets during the last glacial cycle to the LGM based on the available evidence. We consider these ice configurations as a series of testable hypotheses for future work to address and refine. These results are important for use across a range of disciplines including ice sheet modelling, palaeoclimatology and archaeology and are available digitally.

The North American Ice Sheet Complex (NAISC; consisting of the Laurentide, Cordilleran and Innuitian ice sheets) was the largest ice mass to repeatedly grow and decay in the Northern Hemisphere during the Quaternary. Understanding its pattern of retreat following the Last Glacial Maximum is critical for studying many facets of the Late Quaternary, including ice sheet behaviour, the evolution of Holocene landscapes, sea level, atmospheric circulation, and the peopling of the Americas. Currently, the most up-to-date and authoritative margin chronology for the entire ice sheet complex is featured in two publications (Geological Survey of Canada Open File 1574 [Dyke et al., 2003]; 'Quaternary Glaciations - Extent and Chronology, Part II' [Dyke, 2004]). These often-cited datasets track ice margin recession in 36 time slices spanning 18 ka to 1 ka (all ages in uncalibrated radiocarbon years) using a combination of geomorphology, stratigraphy and radiocarbon dating. However, by virtue of being over 15 years old, the ice margin chronology requires updating to reflect new work and important revisions. This paper updates the aforementioned 36 ice margin maps to reflect new data from regional studies. We also update the original radiocarbon dataset from the 2003/2004 papers with 1541 new ages to reflect work up to and including 2018. A major revision is made to the 18 ka ice margin, where Banks and Eglinton islands (once considered to be glacial refugia) are now shown to be fully glaciated. Our updated 18 ka ice sheet increased in areal extent from 17.81 to 18.37 million km(2), which is an increase of 3.1% in spatial coverage of the NAISC at that time. Elsewhere, we also summarize, region-by-region, significant changes to the deglaciation sequence. This paper integrates new information provided by regional experts and radiocarbon data into the deglaciation sequence while maintaining consistency with the original ice margin positions of Dyke et al. (2003) and Dyke (2004) where new information is lacking; this is a pragmatic solution to satisfy the needs of a Quaternary research community that requires up-to-date knowledge of the pattern of ice margin recession of what was once the world's largest ice mass. The 36 updated isochrones are available in PDF and shapefile format, together with a spreadsheet of the expanded radiocarbon dataset (n = 5195 ages) and estimates of uncertainty for each interval. (C) 2020 Elsevier Ltd. All rights reserved.

Significance During the Holocene (11,600 y ago to present), northern peatlands accumulated significant C stocks over millennia. However, virtually nothing is known about peatlands that are no longer in the landscape, including ones formed prior to the Holocene: Where were they, when did they form, and why did they disappear? We used records of peatlands buried by mineral sediments for a reconstruction of peat-forming wetlands for the past 130,000 y. Northern peatlands expanded across high latitudes during warm periods and were buried during periods of glacial advance in northern latitudes. Thus, peat accumulation and burial represent a key long-term C storage mechanism in the Earth system., Glacial−interglacial variations in CO2 and methane in polar ice cores have been attributed, in part, to changes in global wetland extent, but the wetland distribution before the Last Glacial Maximum (LGM, 21 ka to 18 ka) remains virtually unknown. We present a study of global peatland extent and carbon (C) stocks through the last glacial cycle (130 ka to present) using a newly compiled database of 1,063 detailed stratigraphic records of peat deposits buried by mineral sediments, as well as a global peatland model. Quantitative agreement between modeling and observations shows extensive peat accumulation before the LGM in northern latitudes (>40°N), particularly during warmer periods including the last interglacial (130 ka to 116 ka, MIS 5e) and the interstadial (57 ka to 29 ka, MIS 3). During cooling periods of glacial advance and permafrost formation, the burial of northern peatlands by glaciers and mineral sediments decreased active peatland extent, thickness, and modeled C stocks by 70 to 90% from warmer times. Tropical peatland extent and C stocks show little temporal variation throughout the study period. While the increased burial of northern peats was correlated with cooling periods, the burial of tropical peat was predominately driven by changes in sea level and regional hydrology. Peat burial by mineral sediments represents a mechanism for long-term terrestrial C storage in the Earth system. These results show that northern peatlands accumulate significant C stocks during warmer times, indicating their potential for C sequestration during the warming Anthropocene.

Well-dated paleorecords from periods prior to the Last Glacial Maximum (LGM) are important for validating models of ice sheet build-up and growth. However, owing to glacial erosion, most Late Pleistocene records lie outside of the previously glaciated region, which limits their ability to inform about the dynamics of paleo-ice sheets. Here, we evaluate new and previously published chronology data from the Missinaibi Formation, a Pleistocene-aged deposit in the Hudson Bay Lowlands (HBL), Canada, located near the geographic center of the Laurentide Ice Sheet (LIS). Available radiocarbon (AMS = 44, conventional = 36), amino acid (n = 13), uranium-thorium (U-Th, n = 14), thermoluminescence (TL, n = 15) and optically stimulated luminescence (OSL, n = 5) data suggest that an ice-free HBL may have been possible during parts of Marine Isotope Stage 7 (MIS 7; ca. 243,000 to ca. 190,000 yr BP), MIS 5 (ca. 130,000 to ca. 71,000 yr BP) and MIS 3 (ca. 29,000 to ca. 57,000). While MIS 7 and MIS 5 are well-documented interglacial periods, the development of peat, forest bed and fluvial deposits dating to MIS 3 (n = 20 radiocarbon dates; 4 TL dates, 3 OSL dates), suggests that the LIS retreated and remained beyond, or somewhere within, the boundaries of the HBL during this interstadial. Ice sheet models approximate the margin of the LIS to Southern Ontario during this time, which is 700 km south of the HBL. Therefore, if correct, our data help constrain a significantly different configuration and dynamicity for the LIS than previously modelled. We can find no chronological basis to discount the MIS 3 age assignments. However, since most data originate from radiocarbon dates lying close to the reliable limit of this geochronometer, future work on dating the Missinaibi Formation using other geochronological methods (e.g. U-Th, OSL) is necessary in order to confirm the age estimates and strengthen the boundaries of the LIS during this period.

We examined late Holocene (ca. 3300 yr BP to present-day) climate variability in the central Northwest Territories (Canadian Subarctic) using a diatom and sedimentological record from Danny’s Lake (63.48oN, 112.54oW), located 40 km southwest of the modern-day tree- line. High-resolution sampling paired with a robust age model (25 radiocarbon dates) allowed for the examination of both lake hydroecological conditions (30-year intervals; dia- toms) and sedimentological changes in the watershed (12-year intervals; grain size records) over the late Holocene. Time series analysis of key lake ecological indicators (diatom spe- cies Aulacoseira alpigena, Pseudostaurosira brevistriata and Achnanthidium minutissimum) and sedimentological parameters, reflective of catchment processes (coarse silt fraction), suggests significant intermittent variations in turbidity, pH and light penetration within the lake basin. In the diatom record, we observed discontinuous periodicities in the range of ca. 69, 88–100, 115–132, 141–188, 562, 750 and 900 years (>90% and >95% confidence inter- vals), whereas the coarse silt fraction was characterized by periodicities in the >901 and 95% confidence interval). Periodicities in the proxy data from the Danny’s Lake sediment core align with changes in total solar irradiance over the past ca. 3300 yr BP and we hypothesize a link to the Suess Cycle, Gleissberg Cycle and Pacific Decadal Oscilla- tion via occasional inland propagation of shifting air masses over the Pacific Ocean. This research represents an important baseline study of the underlying causes of climate vari- ability in the Canadian Subarctic and provides details on the long-term climate variability that has persisted in this region through the past three thousand years.