Dynamical response of the southwestern Laurentide Ice Sheet to rapid Bølling–Allerød warming

The shift in climate that occurred between the Last Glacial Maximum (LGM) and the Early Holocene (ca. 18–12 kyr BP) displayed rates of temperature increase similar to present-day warming trends. The most rapid recorded changes in temperature occurred during the abrupt climate oscillations known as the Bølling–Allerød interstadial (14.7–12.9 kyr BP) and the Younger Dryas stadial (12.9–11.7 kyr BP). Reconstructing ice sheet dynamics during these climate oscillations provides the opportunity to assess long-term ice sheet evolution in reaction to a rapidly changing climate. Here, we use glacial geomorphological inversion methods (flowsets) to reconstruct the ice flow dynamics and the marginal retreat pattern of the southwestern sector of the Laurentide Ice Sheet (SWLIS). We combine our reconstruction with a recently compiled regional deglaciation chronology to depict ice flow dynamics that encompass the time period from pre-LGM to the Early Holocene. Our reconstruction portrays three macroscale reorganizations in the orientation and dynamics of ice streaming followed by regional deglaciation associated with rapid warming during the Bølling–Allerød interstadial. Initial westward flow is documented, likely associated with an early set of ice streams that formed during the advance to the LGM. During the LGM ice streaming displays a dominant north to south orientation. Ice sheet thinning at ∼15 ka is associated with a macroscale reorganization in ice stream flow, with a complex of ice streams recording south-eastward flow. A second macroscale reorganization in ice flow is then observed at ∼14 ka, in which southwestern ice flow is restricted to the Hay, Peace, Athabasca, and Churchill river lowlands. Rates of ice sheet retreat then slowed considerably during the Younger Dryas stadial; at this time, the ice margin was situated north of the Canadian Shield boundary and ice flow continued to be sourced from the northeast. Resulting from these changes in ice sheet dynamics, we recognize a three-part pattern of deglacial landform zonation within the SWLIS characterized by active ice margin recession, stagnation and downwasting punctuated by local surging (terrestrial ice sheet collapse): the outer deglacial zone contains large recessional moraines aligned with the direction of active ice margin retreat; the intermediate deglacial zone contains large regions of hummocky and stagnation terrain, in some areas crosscut by the signature of local surges, reflecting punctuated stagnation and downwasting; and the inner deglacial zone contains inset recessional moraines demarcating progressive regional ice margin retreat. We attribute these macroscale changes in ice flow geometry and associated deglacial behaviour to external climatic controls during the Bølling–Allerød and Younger Dryas but also recognize the role of internal (glaciological, lithological and topographic) controls in SWLIS dynamics

Reconstructing dynamics of the Baltic Ice Stream Complex during deglaciation of the Last Scandinavian Ice Sheet

Landforms left behind by the last Scandinavian Ice Sheet (SIS) offer an opportunity to investigate controls governing ice sheet dynamics. Terrestrial sectors of the ice sheet have received considerable attention from landform and stratigraphic investigations. In contrast, despite its geographical importance, the Baltic Sea remains poorly constrained due to limitations in bathymetric data. Both ice-sheet-scale investigations and regional studies at the southern periphery of the SIS have considered the Baltic depression to be a preferential route for ice flux towards the southern ice margin throughout the last glaciation. During the deglaciation the Baltic depression hosted the extensive Baltic Ice Lake, which likely exerted a considerable control on ice dynamics. Here we investigate the Baltic depression using newly available bathymetric data and peripheral topographic data. These data reveal an extensive landform suite stretching from Denmark in the west to Estonia in the east and from the southern European coast to the Åland Sea, comprising an area of 0.3 million km2. We use these landforms to reconstruct aspects of the ice dynamic history of the Baltic sector of the ice sheet. Landform evidence indicates a complex retreat pattern that changes from lobate ice margins with splaying lineations to parallel mega-scale glacial lineations (MSGLs) in the deeper depressions of the Baltic Basin. Ice margin still-stands on underlying geological structures indicate the likely importance of pinning points during deglaciation, resulting in a stepped retreat signal. Over the span of the study area we identify broad changes in the ice flow direction, ranging from SE–NW to N–S and then to NW–SE. MSGLs reveal distinct corridors of fast ice flow (ice streams) with widths of 30 km and up to 95 km in places, rather than the often-interpreted Baltic-wide (300 km) accelerated ice flow zone. These smaller ice streams are interpreted as having operated close behind the ice margin during late stages of deglaciation. Where previous ice-sheet-scale investigations inferred a single ice source, our mapping identifies flow and ice margin geometries from both Swedish and northern Bothnian sources. We anticipate that our landform mapping and interpretations may be used as a framework for more detailed empirical studies by identifying targets to acquire high-resolution bathymetry and sediment cores and also for comparison with numerical ice sheet modelling.</p

On the Reconstruction of Palaeo-Ice Sheets: Recent Advances and Future Challenges

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

Exploring the ingredients required to successfully model the placement, generation, and evolution of ice streams in the British-Irish Ice Sheet

Ice stream evolution is a major uncertainty in projections of the future of the Greenland and Antarctic Ice sheets. Accurate simulation of ice stream evolution requires an understanding of a number of “ingredients” that control the location and behaviour of ice stream flow. Here, we test the influence of geothermal heat flux, grid resolution, and bed hydrology on simulated ice streaming. The palaeo-record provides snapshots of ice stream evolution, with a particularly well constrained ice sheet being the British-Irish Ice Sheet (BIIS). We implement a new basal sliding scheme coupled with thermo-mechanics into the BISICLES ice sheet model, to simulate the evolution of the BIIS ice streams. We find that the simulated location and spacing of ice streams matches well with the empirical reconstructions of ice stream flow in terms of position and direction when simple bed hydrology is included. We show that the new basal sliding scheme allows the accurate simulation for the majority of BIIS ice streams. The extensive empirical record of the BIIS has allowed the testing of model inputs, and has helped demonstrate the skill of the ice sheet model in simulating the evolution of the location, spacing, and migration of ice streams through millennia. Simulated ice streams also prompt new empirical mapping of features indicative of streaming in the North Channel region. Ice sheet model development has allowed accurate simulation of the palaeo record, and allows for improved modelling of future ice stream behaviour

An updated radiocarbon-based ice margin chronology for the last deglaciation of the North American Ice Sheet Complex

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 km2, 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

Constraining bedrock erosion during extreme flood events

The importance of high-magnitude, short-lived flood events in controlling the evolution of bedrock landscapes is not well understood. During such events, erosion processes can shift from one regime to another upon the passing of thresholds, resulting in abrupt landscape changes that can have a long lasting legacy on landscape morphology. Geomorphological mapping and topographic analysis document the evidence for, and impact of, extreme flood events within the Jökulsárgljúfur canyon (North-East Iceland). Surface exposure dating using cosmogenic 3He of fluvially sculpted bedrock surfaces determines the timing of the floods that eroded the canyon and helps constrain the mechanisms of bedrock erosion during these events. Once a threshold flow depth has been exceeded, the dominant erosion mechanism becomes the toppling and transportation of basalt lava columns and erosion occurs through the upstream migration of knickpoints. Surface exposure ages allow identification of three periods of rapid canyon cutting during erosive flood events about 9, 5 and 2 ka ago, when multiple active knickpoints retreated large distances (> 2 km), each leading to catastrophic landscape change within the canyon. A single flood event ~9 ka ago formed, and then abandoned, Ásbyrgi canyon, eroding 0.14 km3 of rock. Flood events ~5 and ~2 ka ago eroded the upper 5 km of the Jökulsárgljúfur canyon through the upstream migration of vertical knickpoints such as Selfoss, Dettifoss and Hafragilsfoss. Despite sustained high discharge of sediment-rich glacial meltwater (ranging from 100 to 500 m3 s-1); there is no evidence for a transition to an abrasion-dominated erosion regime since the last erosive flood: the vertical knickpoints have not diffused over time and there is no evidence of incision into the canyon floor. The erosive signature of the extreme events is maintained in this landscape due to the nature of the bedrock, the discharge of the river, large knickpoints and associated plunge pools. The influence of these controls on the dynamics of knickpoint migration and morphology are explored using an experimental study. The retreat rate of knickpoints is independent of both mean discharge, and temporal variability in the hydrograph. The dominant control on knickpoint retreat is the knickpoint form which is set by the ratio of channel flow depth to knickpoint height. Where the knickpoint height is five times greater than the flow depth, the knickpoints developed undercutting plunge pools, accelerating the removal of material from the knickpoint base and the overall retreat rate. Smaller knickpoints relative to the flow depth were more likely to diffuse from a vertical step into a steepened reach or completely as the knickpoint retreated up the channel. These experiments challenge the established assumption in models of landscape evolution that a simple relationship exists between knickpoint retreat and discharge/drainage area. In order to fully understand how bedrock channels, and thus landscapes, respond and recover to transient forcing, further detailed study of the mechanics of erosion processes at knickpoints is required

The glacial geomorphology of the Mackenzie Mountains region, Canada

ABSTRACTDuring the Last Glacial Maximum, the Mackenzie Mountains region was glaciated by three distinct ice sources; the Laurentide Ice Sheet, the Cordilleran Ice Sheet, and independent montane glaciers. Rapid ice sheet thinning of the Laurentide-Cordilleran ice saddle in the south of this region contributed to rapid sea level rise events and influenced the style of deglaciation to the north. The current understanding of the glacial history of the broader region has been established through mapping from aerial imagery and early surveys between the early 1970s to the 2010s. The central portions of the Mackenzie Mountains have not yet been mapped. We present a new glacial geomorphological map for the Mackenzie Mountains region covering over 220,000 km2. This updated geomorphological map will form the basis of future work to reconstruct the former maximum ice extents, flow dynamics, and retreat pattern

Ice streams in the Laurentide Ice Sheet: a new mapping inventory

Rapidly flowing ice streams dominate the drainage of continental ice sheets and are a key component of their mass balance. Due to their potential impact on sea level, their activity in the Antarctic and Greenland Ice Sheets has undergone detailed scrutiny in recent decades. However, these observations only cover a fraction of their ‘life-span’ and the subglacial processes that facilitate their rapid flow are very difficult to observe. To circumvent these problems, numerous workers have highlighted the potential of investigating palaeo-ice streams tracks, preserved in the landform and sedimentary record of former ice sheets. As such, it is becoming increasingly important to know where and when palaeo-ice streams operated. In this paper, we present a new map of ice streams in the North American Laurentide Ice Sheet (LIS; including the Innuitian Ice Sheet), which was the largest of the ephemeral Pleistocene ice sheets and where numerous ice streams have been identified. We compile previously published evidence of ice stream activity and complement it with new mapping to generate the most complete and consistent mapping inventory to date. The map depicts close to three times as many ice streams (117 in total) compared to previous inventories, and categorises them according to the evidence they left behind, with some locations more speculative than others. The map considerably refines our understanding of LIS dynamics, but there is a clear requirement for improved dating of ice stream activity

Investigating absolute chronologies of glacial advances in the NW sector of the Cordilleran Ice Sheet with terrestrial in situ cosmogenic nuclides

Geologic mapping in Yukon Territory, Canada, over the past 100 years has revealed a consistent pattern of diminishing Cordilleran Ice Sheet (CIS) extent during successively younger glaciations. Although this pattern is generally accepted, there is still much uncertainty about the number of glaciations, their ages, and the dynamics of the different lobes that constituted the digitate ice sheet margin, their subglacial thermal regimes, and ice thicknesses. We address uncertainties in the timing of glaciation using cosmogenic nuclide exposure dating at key localities that straddle several major lobes of the CIS in west- central Yukon Territory. Differences in exposure duration within what are thought to be the same map units are perhaps due to inheritance (older than expected), but more likely result from postglacial shielding (younger than expected) or surface erosion. Despite a significant spread in exposure durations on moraines and within map units, and tending to rely on longest exposure durations on moraines due to postglacial degradation and shielding, our results indicate that the McConnell glacial advance occurred during Marine Oxygen Isotope Stage (MIS) 2, judging from oldest minimum apparent exposure ages of 15.7 ± 1.5 and 17.7 ± 1.6 ka, a Gladstone glacial advance occurred before 51.8 ` 4.7 ka (MIS 4) and Reid glacial advances before 79.8 ± 7.3 and 82.8 ± 7.5 ka (consistent with MIS 6). Traces of even older glacial advances predate 100 ka (107.5 ± 9.9 - 154.3 ± 14.2 ka)

Glacial geomorphology of the northwest Laurentide Ice Sheet on the northern Interior Plains and western Canadian Shield, Canada

The majority of the Northwest Territories of mainland Canada was covered by the Laurentide Ice Sheet during the Last Glacial Maximum. The increasing coverage of high resolution remotely sensed data provides new opportunities to map the glacial geomorphology and study the glacial history of this remote location. Here we present a comprehensive map of glacial landforms within the northern Interior Plains and adjacent areas of the Canadian Shield, comprising around 6% of the Laurentide Ice Sheet bed. Twelve landform types were mapped from the high resolution ArcticDEM: ice flow parallel lineations, subglacial ribs, crevasse-squeeze ridges, major and minor moraine crests, hummocky terrain complexes and ridges, shear margin moraines, major, minor and lateral and submarginal meltwater channels, esker ridges and complexes, glaciofluvial complexes, perched deltas, raised shorelines and aeolian dunes. Together, these landforms provide a record of the highly dynamic behaviour of the northwest sector of the Laurentide Ice Sheet