Search

In situ 14C in quartz provides a recently developed tool to date exposure of bedrock surfaces up to ~25 000 years. From outcrops located in east-central Sweden, we test the accuracy of in situ 14C dating against (i) a relative sea level (RSL) curve constructed from radiocarbon dating of organic material in isolation basins, and (ii) the timing of local deglaciation constructed from a clay varve chronology complemented with radiocarbon dating. Five samples of granitoid bedrock were taken along an elevation transect extending southwestwards from the Baltic Sea coast near Forsmark. Because these samples derive from bedrock outcrops positioned below the highest postglacial shoreline, they target the timing of progressive landscape emergence above sea level. In contrast, in situ 14C concentrations in an additional five samples taken from granitoid outcrops above the highest postglacial shoreline, located 100 km west of Forsmark, should reflect local deglaciation ages. The ten in situ 14C measurements provide robust age constraints that, within uncertainties, compare favorably with the RSL curve and with the local deglaciation chronology. These data demonstrate the utility of in situ 14C to accurately date ice sheet deglaciation, and durations of postglacial exposure, in regions where cosmogenic 10Be and 26Al routinely return complex exposure results. [ABSTRACT FROM AUTHOR]

In low relief Precambrian gneiss terrain in eastern Sweden, abraded bedrock surfaces were ripped apart by the Fennoscandian Ice Sheet. The resultant boulder spreads are covers of large, angular boulders, many with glacial transport distances of 1–100 m. Boulder spreads occur alongside partly disintegrated roches moutonnées and associated fracture caves, and are associated with disrupted bedrock, which shows extensive fracture dilation in the near surface. These features are distributed in ice-flow parallel belts up to 10 km wide and extend over distances of >500 km. Our hypothesis is that the assemblage results from (1) hydraulic jacking and bedrock disruption, (2) subglacial ripping and (3) displacement, transport and final deposition of boulders. Soft sediment fills indicate jacking and dilation of pre-existing bedrock fractures by groundwater overpressure below the ice sheet. Overpressure reduces frictional resistance along fractures. Where ice traction overcomes this resistance, the rock mass strength is exceeded, resulting in disintegration of rock surfaces and ripping apart into separate blocks. Further movement and deposition create boulder spreads and moraines. Short boulder transport distances and high angularity indicate that glacial ripping operated late in the last deglaciation. The depths of rock mobilized in boulder spreads are estimated as 1–4 m. This compares with 0.6–1.6 m depths of erosion during the last glaciation derived from cosmogenic nuclide inventories of samples from bedrock surfaces without evidence of disruption. Glacially disrupted and ripped bedrock is also made ready for removal by future ice sheets. Hence glacial ripping is a highly effective process of glacial erosion. [ABSTRACT FROM AUTHOR]

The most recent deglaciation of the North American Ice Sheet Complex (NAISC: comprising the Innuitian, Cordilleran, and Laurentide ice sheets) offers a broad perspective from which to analyze the timing and rate of ice retreat, deglacial sea-level rise, and abrupt climate change events. Previous efforts to portray the retreat of the NAISC have been focused largely on minimum-limiting radiocarbon ages and ice margin location(s) tied to deglacial landforms that were not, for the most part, chronologically constrained. Here, we present the first version of North American Deglaciation Isochrones (NADI-1) spanning 25 to 1 ka in calendar years before present. Key new features of this work are (i) the incorporation of cosmogenic nuclide data, which offer a direct constraint on the timing of ice recession; (ii) presentation of all data and time-steps in calendar years; (iii) optimal, minimum, and maximum ice extents for each time-step that are designed to capture uncertainties in the ice margin position, and; (iv) extensive documentation and justification for the placement of each ice margin. Our data compilation includes 2229 measurements of 10Be, 459 measurements of 26Al and 35 measurements of 36Cl from a variety of settings, including boulders, bedrock surfaces, cobbles, pebbles, and sediments. We also updated a previous radiocarbon dataset (n = 4947), assembled luminescence ages (n = 397) and gathered uranium-series data (n = 2). After scrutiny of the geochronological dataset, we consider >90% of data to be reliable or likely reliable. Key findings include (i) a highly asynchronous maximum glacial extent in North America, occurring as early as 27 ka to as late as 17 ka, within and between ice sheets. In most marine realms, extension of the ice margin to the continental shelf break at 25 ka is somewhat speculative because it is based on undated and spatially scattered ice stream and geomorphic evidence; (ii) detachment of the Laurentide and Cordilleran ice sheets took place gradually via southerly and northerly 'unzipping' of the ice masses, starting at 17.5 ka and ending around 14 ka; (iii) the final deglaciation of Hudson Bay began at 8.5 ka, with the collapse completed by 8 ka. The maximum extent of ice during the last glaciation occurred at 22 ka and covered 15,470,000 km2. All North American ice sheets merged at 22 ka for the first time in the Quaternary. The highly asynchronous Last Glacial Maximum in North America means that our isochrones (starting at 25 ka) capture ice advance across some areas, which is based on limited evidence and is therefore somewhat speculative. In the Supplementary Data, the complete NADI-1 chronology is available in PDF, GIF and shapefile format, together with additional visualizations and spreadsheets of geochronological data. The NADI-1 shapefiles are also available at https://doi.org/10.5281/zenodo.8161764. • Uncertainties captured in optimal, minimum, and maximum ice extents. • Detailed documentation for the placement of each ice margin. • NADI-1 ice margins available in.shp,.pdf and.gif format for 49 timesteps. • Supplementary Data: additional visualizations, geochronological data. [ABSTRACT FROM AUTHOR]

The Eurasian ice sheet complex (EISC) was the third largest ice mass during the Last Glacial Maximum with a span of over 4500 km and responsible for around 20 m of eustatic sea-level lowering. Whilst recent terrestrial and marine empirical insights have improved understanding of the chronology, pattern and rates of retreat of this vast ice sheet, a concerted attempt to model the deglaciation of the EISC honouring these new constraints is conspicuously lacking. Here, we apply a first-order, thermomechanical ice sheet model, validated against a diverse suite of empirical data, to investigate the retreat of the EISC after 23 ka BP, directly extending the work of Patton et al. (2016) who modelled the build-up to its maximum extent. Retreat of the ice sheet complex was highly asynchronous, reflecting contrasting regional sensitivities to climate forcing, oceanic influence, and internal dynamics. Most rapid retreat was experienced across the Barents Sea sector after 17.8 ka BP when this marine-based ice sheet disintegrated at a rate of ∼670 gigatonnes per year (Gt a −1 ) through enhanced calving and interior dynamic thinning, driven by oceanic/atmospheric warming and exacerbated by eustatic sea-level rise. From 14.9 to 12.9 ka BP the EISC lost on average 750 Gt a −1 , peaking at rates >3000 Gt a −1 , roughly equally partitioned between surface melt and dynamic losses, and potentially contributing up to 2.5 m to global sea-level rise during Meltwater Pulse 1A. Independent glacio-isostatic modelling constrained by an extensive inventory of relative sea-level change corroborates our ice sheet loading history of the Barents Sea sector. Subglacial conditions were predominately temperate during deglaciation, with over 6000 subglacial lakes predicted along with an extensive subglacial drainage network. Moreover, the maximum EISC and its isostatic footprint had a profound impact on the proglacial hydrological network, forming the Fleuve Manche mega-catchment which had an area of ∼2.5 × 10 6 km 2 and drained the present day Vistula, Elbe, Rhine and Thames rivers through the Seine Estuary. During the Bølling/Allerød oscillation after c. 14.6 ka BP, two major proglacial lakes formed in the Baltic and White seas, buffering meltwater pulses from eastern Fennoscandia through to the Younger Dryas when these massive proglacial freshwater lakes flooded into the North Atlantic Ocean. Deglaciation temporarily abated during the Younger Dryas stadial at 12.9 ka BP, when remnant ice across Svalbard, Franz Josef Land, Novaya Zemlya, Fennoscandia and Scotland experienced a short-lived but dynamic re-advance. The final stage of deglaciation converged on present day ice cover around the Scandes mountains and the Barents Sea by 8.7 ka BP, although the phase-lagged isostatic recovery still continues today. [ABSTRACT FROM AUTHOR]

Reconstructing the growth and decay of palaeo-ice sheets is critical to understanding mechanisms of global climate change and associated sea-level fluctuations in the past, present and future. The significance of palaeo-ice sheets is further underlined by the broad range of disciplines concerned with reconstructing their behaviour, many of which have undergone a rapid expansion since the 1980s. In particular, there has been a major increase in the size and qualitative diversity of empirical data used to reconstruct and date ice sheets, and major improvements in our ability to simulate their dynamics in numerical ice sheet models. These developments have made it increasingly necessary to forge interdisciplinary links between sub-disciplines and to link numerical modelling with observations and dating of proxy records. The aim of this paper is to evaluate recent developments in the methods used to reconstruct ice sheets and outline some key challenges that remain, with an emphasis on how future work might integrate terrestrial and marine evidence together with numerical modelling. Our focus is on pan-ice sheet reconstructions of the last deglaciation, but regional case studies are used to illustrate methodological achievements, challenges and opportunities. Whilst various disciplines have made important progress in our understanding of ice-sheet dynamics, it is clear that data-model integration remains under-used, and that uncertainties remain poorly quantified in both empirically-based and numerical ice-sheet reconstructions. The representation of past climate will continue to be the largest source of uncertainty for numerical modelling. As such, palaeo-observations are critical to constrain and validate modelling. State-of-the-art numerical models will continue to improve both in model resolution and in the breadth of inclusion of relevant processes, thereby enabling more accurate and more direct comparison with the increasing range of palaeo-observations. Thus, the capability is developing to use all relevant palaeo-records to more strongly constrain deglacial (and to a lesser extent pre-LGM) ice sheet evolution. In working towards that goal, the accurate representation of uncertainties is required for both constraint data and model outputs. Close cooperation between modelling and data-gathering communities is essential to ensure this capability is realised and continues to progress. [ABSTRACT FROM AUTHOR]