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]

Numerical models predict that discharge from the polar ice sheets will become the largest contributor to sea-level rise over the coming centuries. However, the predicted amount of ice discharge and associated thinning depends on how well ice sheet models reproduce glaciological processes, such as ice flow in regions of large topographic relief, where ice flows around bedrock summits (i.e. nunataks) and through outlet glaciers. The ability of ice sheet models to capture long-term ice loss is best tested by comparing model simulations against geological data. A benchmark for such models is ice surface elevation change, which has been constrained empirically at nunataks and along margins of outlet glaciers using cosmogenic exposure dating. However, the usefulness of this approach in quantifying ice sheet thinning relies on how well such records represent changes in regional ice surface elevation. Here we examine how ice surface elevations respond to the presence of strong topographic relief that acts as an obstacle by modelling ice flow around and between idealised nunataks during periods of imposed ice sheet thinning. We find that, for realistic Antarctic conditions, a single nunatak can exert an impact on ice thickness over 20 km away from its summit, with its most prominent effect being a local increase (decrease) of the ice surface elevation of hundreds of metres upstream (downstream) of the obstacle. A direct consequence of this differential surface response for cosmogenic exposure dating is a delay in the time of bedrock exposure upstream relative to downstream of a nunatak summit. A nunatak elongated transversely to ice flow is able to increase ice retention and therefore impose steeper ice surface gradients, while efficient ice drainage through outlet glaciers produces gentler gradients. Such differences, however, are not typically captured by continent-wide ice sheet models due to their coarse grid resolutions. Their inability to capture site-specific surface elevation changes appears to be a key reason for the observed mismatches between the timing of ice-free conditions from cosmogenic exposure dating and model simulations. We conclude that a model grid refinement over complex topography and information about sample position relative to ice flow near the nunatak are necessary to improve data–model comparisons of ice surface elevation and therefore the ability of models to simulate ice discharge in regions of large topographic relief. [ABSTRACT FROM AUTHOR]

Studying the response of the Antarctic ice sheets during periods when climate conditions were similar to the present can provide important insights into current observed changes and help identify natural drivers of ice sheet retreat. In this context, the marine isotope substage 11c (MIS11c) interglacial offers a suitable scenario, given that during its later portion orbital parameters were close to our current interglacial. Ice core data indicate that warmer-than-present temperatures lasted for longer than during other interglacials. However, the response of the Antarctic ice sheets and their contribution to sea level rise remain unclear. We explore the dynamics of the Antarctic ice sheets during this period using a numerical ice sheet model forced by MIS11c climate conditions derived from climate model outputs scaled by three glaciological and one sedimentary proxy records of ice volume. Our results indicate that the East and West Antarctic ice sheets contributed 4.0–8.2 m to the MIS11c sea level rise. In the case of a West Antarctic Ice Sheet collapse, which is the most probable scenario according to far-field sea level reconstructions, the range is reduced to 6.7–8.2 m independently of the choices of external sea level forcing and millennial-scale climate variability. Within this latter range, the main source of uncertainty arises from the sensitivity of the East Antarctic Ice Sheet to a choice of initial ice sheet configuration. We found that the warmer regional climate signal captured by Antarctic ice cores during peak MIS11c is crucial to reproduce the contribution expected from Antarctica during the recorded global sea level highstand. This climate signal translates to a modest threshold of 0.4 ∘ C oceanic warming at intermediate depths, which leads to a collapse of the West Antarctic Ice Sheet if sustained for at least 4000 years. [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]

After more than a century of geological research, the Cordilleran ice sheet of North America remains among the least understood in terms of its former extent, volume, and dynamics. Because of the mountainous topography on which the ice sheet formed, geological studies have often had only local or regional relevance and shown such a complexity that ice-sheet-wide spatial reconstructions of advance and retreat patterns are lacking. Here we use a numerical ice sheet model calibrated against field-based evidence to attempt a quantitative reconstruction of the Cordilleran ice sheet history through the last glacial cycle. A series of simulations is driven by time-dependent temperature offsets from six proxy records located around the globe. Although this approach reveals large variations in model response to evolving climate forcing, all simulations produce two major glaciations during marine oxygen isotope stages 4 (62.2-56.9 ka) and 2 (23.2-16.9 ka). The timing of glaciation is better reproduced using temperature reconstructions from Greenland and Antarctic ice cores than from regional oceanic sediment cores. During most of the last glacial cycle, the modelled ice cover is discontinuous and restricted to high mountain areas. However, widespread precipitation over the Skeena Mountains favours the persistence of a central ice dome throughout the glacial cycle. It acts as a nucleation centre before the Last Glacial Maximum and hosts the last remains of Cordilleran ice until the middle Holocene (6.7 ka). [ABSTRACT FROM AUTHOR]

The Cordilleran Ice Sheet (CIS) covered most of British Columbia and southern Yukon Territory at the local Last Glacial Maximum (lLGM) during Marine Oxygen Isotope Stage 2. However, its subsequent demise is not well understood, particularly at high elevations east of its ocean-terminating margin. We present 10Be exposure ages from two high-elevation sites in southern and central British Columbia that help constrain the time of initial deglaciation at these sites. We sampled granodiorite erratics at elevations of 2126-2230 m a.s.l. in the Marble Range and 1608-1785 m a.s.l. in the Telkwa Range at the western margin of the Interior Plateau. The erratics at both sites are near ice-marginal meltwater channels that delineate the local ice surface slope and thus the configuration of the ice sheet during deglaciation. The locations of the erratics and their relations to meltwater channels ensure that the resulting 10Be ages date CIS deglaciation and not the retreat of local montane glaciers. Our sample sites emerged above the surface of the CIS as its divide migrated westward from the Interior Plateau to the axis of the Coast Mountains. Two of the four samples from the summit area of the Marble Range yielded apparent exposure ages of 14.0 ± 0.7 and 15.2 ± 0.8 ka. These ages are 1.8-3.0 ka younger than the well-established lLGM age of ca 17 ka for the Puget lobe of the CIS in Washington State; they are 1.7 ka younger than the lLGM age for the Puget lobe if a snow-shielding correction to their uncertainty-weighted mean age is applied. The other two samples yielded much older apparent exposure ages (20.6 ± 1.4 and 33.0 ± 1.5 ka), indicating the presence of inherited isotopes. Four samples collected from the summit area of the Telkwa Range in the Hazelton Mountains yielded well clustered apparent exposure ages of 10.1 ± 0.6, 10.2 ± 0.7, 10.4 ± 0.5, and 11.5 ± 1.1 ka. Significant present-day snow cover introduces a large uncertainty in the apparent exposure ages from this site. A snow-shielding correction based on present-day snow cover data increases the uncertainty-weighted mean exposure age of the Telkwa Range erratics to 12.4 ± 0.7 ka, consistent with deglacial 14C ages from areas near sea level to the west. Our exposure ages show a thinning of the southern portion of the CIS shortly after the lLGM and persistence of a remnant mountain ice cap in the central Coast Mountains into the Younger Dryas Chronozone. Our data also show that the summit area of the Marble Range was ice-covered during the lLGM. The presence of an ice body of considerable dimension in north-central British Columbia until, or possibly even after, the Younger Dryas highlights the need for geomorphological and geochronological studies of the ice dispersal centre over the Skeena Mountains in northwest British Columbia and the need for better understanding of the response of the CIS to Lateglacial climate fluctuations. [ABSTRACT FROM AUTHOR]

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). [ABSTRACT FROM AUTHOR]

The Cordilleran Ice Sheet ( CIS) covered much of the mountainous northwestern part of North America at least several times during the Pleistocene. The pattern and timing of its growth and decay are, however, poorly understood. Here, we present a reconstruction of the pattern of ice-sheet retreat in central British Columbia at the end of the last glaciation based on a palaeoglaciological interpretation of ice-marginal meltwater channels, eskers and deltas mapped from satellite imagery and digital elevation models. A consistent spatial pattern of high-elevation (1600-2400 m a.s.l.), ice-marginal meltwater channels is evident across central British Columbia. These landforms indicate the presence of ice domes over the Skeena Mountains and the central Coast Mountains early during deglaciation. Ice sourced in the Coast Mountains remained dominant over the southern and east-central parts of the Interior Plateau during deglaciation. Our reconstruction shows a successive westward retreat of the ice margin from the western foot of the Rocky Mountains, accompanied by the formation and rapid evolution of a glacial lake in the upper Fraser River basin. The final stage of deglaciation is characterized by the frontal retreat of ice lobes through the valleys of the Skeena and Omineca Mountains and by the formation of large esker systems in the most prominent topographic lows of the Interior Plateau. We conclude that the CIS underwent a large-scale reconfiguration early during deglaciation and was subsequently diminished by thinning and complex frontal retreat towards the Coast Mountains. [ABSTRACT FROM AUTHOR]

Abstract: Deglaciation chronologies for some sectors of former ice sheets are relatively poorly constrained because of the paucity of features or materials traditionally used to constrain the timing of deglaciation. In areas without good deglaciation varve chronologies and/or without widespread occurrence of material that indicates the start of earliest organic radiocarbon accumulations suitable for radiocarbon dating, typically only general patterns and chronologies of deglaciation have been deduced. However, mid-latitude ice sheets that had warm-based conditions close to their margins often produced distinctive deglaciation landform assemblages, including eskers, deltas, meltwater channels and aligned lineation systems. Because these features were formed or significantly altered during the last glaciation, boulder or bedrock samples from them have the potential to yield reliable deglaciation ages using terrestrial cosmogenic nuclides (TCN) for exposure age dating. Here we present the results of a methodological study designed to examine the consistency of TCN-based deglaciation ages from a range of deglaciation landforms at a site in northern Norway. The strong coherence between exposure ages across several landforms indicates great potential for using TCN techniques on features such as eskers, deltas and meltwater channels to enhance the temporal resolution of ice-sheet deglaciation chronologies over a range of spatial scales. [Copyright &y& Elsevier]

Please clickhereto download the map associated with this article. The glacial geomorphology of the central Tibetan Plateau was mapped over 285,000 km2.Here we present a map covering 135,000 km at a scale of 1:660,000. The glacial geomorphology was mapped using 15 and 30 m resolution Landsat 7 ETM+ satellite imagery, a 90 m resolution SRTM digital elevation model, and satellite and aerial images displayed in Google Earth. Four landform categories were discernible and mapped; glacial valleys, marginal moraines, glacial lineations, and hummocky terrain. The mapped landforms indicate multiple glacial advances of valley and piedmont glaciers. The mapped landform record lends no support to individual ice centres having coalesced to form a plateau-wide ice sheet. [ABSTRACT FROM PUBLISHER]

Abstract: Exposure durations of glacial landforms in widely separated areas of central Yukon Territory affected by the northern sector of the Cordilleran Ice Sheet (CIS) and alpine glaciers have been determined using cosmogenic 10Be in quartz. The aim of our research is to test previous reconstructions of glacial history and to begin to address the paucity of chronological control for the lateral and vertical extent of the northern CIS. Chronological evidence for CIS expansion predating the Last Glacial Maximum comes from minimum surface exposure durations of c 100 ka for two bedrock samples within the Reid glacial limit, indicating a possible marine Oxygen Isotope Stage (OIS) 6 age for this event, and from minimum exposure durations of about 40 ka for boulders on moraines constructed by alpine glaciers on a nunatak within the McConnell glacial limit (OIS 2), indicating a possible OIS 4 age. High elevation minimum surface exposure durations within the McConnell limit indicate that some areas formerly mapped as nunataks were covered by cold-based ice prior to 30 ka. Montane glaciation in the Mackenzie Mountains, outside the McConnell glacial limit, was contemporaneous with nearby CIS advance at 17 ka, with CIS retreat by 15 ka. Deglaciation of the Tintina Trench, a major ice discharge route, was completed by 12 ka. At this time ice in an adjacent discharge route to the south was still entering higher-elevation valleys in the Pelly Mountains. A Lateglacial readvance may have peaked at ca 10 ka in the Ogilvie Mountains. Considerable variation in ages from individual landforms, and possible complex histories, require additional cosmogenic nuclide measurements to confirm interpretations. [ABSTRACT FROM AUTHOR]

We here present a detailed glacial geomorphological map covering 136,500 km² of the Bayan Har sector of the northeastern Tibetan Plateau - an area previously suggested to have nourished the most extensive Quaternary glaciers of the Tibetan Plateau. The map, presented at a scale of 1:650,000, is based on remote sensing of a 90 m SRTM digital elevation model and 15/30 m Landsat ETM+ satellite imagery. Seven landform types have been mapped; glacial valleys, glacial troughs, glacial lineations, marginal moraines, marginal moraine remnants, meltwater channels and hummocky terrain. A large number of glacial landforms exist, concentrated around mountain blocks protruding above the surrounding plateau area, testifying to former glacial activity. In contrast, large plateau areas of lower altitude lack glacial landforms. The mapped glacial geomorphology indicates multiple former glacial advances primarily by valley and piedmont glaciers, but lends no support to the hypothesis of ice sheet scale glaciation in the area. The presented glacial geomorphological map demonstrates the usefulness of remote sensing techniques for mapping the glacial geomorphology of the Tibetan Plateau, and it will be used for reconstructing the paleoglaciology of the Bayan Har sector of the northeastern Tibetan Plateau. [ABSTRACT FROM AUTHOR]

Abstract: This special issue presents a collection of papers that address a wide range of important challenges and exciting advances in the field of glacial landscape evolution. Primarily, these papers reflect persistent uncertainty that surrounds the mechanisms and timescales of glacial landscape evolution. For example, estimates of the duration of glacial occupancy required for the evolution of characteristic glacial valley forms from previously fluvial landscapes range from 100 kyrs for landscapes beneath large ice sheets (Jamieson et al.) to ~400–600 kyrs for glaciated alpine terrains (Brook et al.). Further, the mechanisms of glacial erosion are debated through analyses of the importance of ice thickness (Brocklehurst et al.; van der Beek and Bourbon), ice surface steepness (Vieira) and, in the case of large ice sheets, the co-evolution of ice sheet thermal regime, dynamics, and subglacial topography (Kleman et al.; Swift et al.). Debate concerning the potential climatic impacts of landscape evolution in alpine terrains is represented by van der Beek and Bourbon, who infer a significant increase in relief as a direct result of glacial erosion, and by Brocklehurst et al. and Heimsath and McGlynn, who demonstrate respectively that glacial relief production can be surprisingly modest and that rates of glacial erosion may be lower than those for fluvial incision. Further confirmation that valleys beneath large ice sheets evolve through selective linear erosion comes from studies that have combined geomorphological evidence with cosmogenic nuclide (Briner et al.) and apatite (U–Th)/He thermochronometry (Swift et al.), and the resulting style of landscape evolution is demonstrated by the antiquity of fjords in East Greenland (Swift et al.) and of deep erosion zones and thick drift covered zones in Fennoscandia (Kleman et al.), although the location of areal scouring zones may be subject to major alteration during single glacial events (Kleman et al.). Another set of papers shows that analyses of glacial lineation systems continue to provide important data on the dynamics of glacial landscape evolution, whether the lineations are formed underneath ice streams (Bradwell et al.; Andreassen et al.) or not (Jansson and Glasser), and whether they indicate intricate patterns of landscape modification (Andreassen et al.) or preservation (Jansson and Glasser). The final three papers address rarely-reported issues relating to landscapes of glacial deposition, including moraine degradation (Putkonen et al.), proglacial hydrogeology (Robinson et al.), and the evolution of hummocky-till topography (Clayton et al.). [Copyright &y& Elsevier]

One of the key elements in reconstructing the palaeoglaciology of the northeastern sector of the Fennoscandian Ice Sheet is the Keiva ice marginal zone (KIZ) along the southern and eastern coast of Kola Peninsula, including the Keiva I and II moraines. From detailed geomorphological mapping of the KIZ, primarily using aerial photographs and satellite images, combined with fieldwork, we observed the following. (1) The moraines display ice contact features on both the Kola side and the White Sea side along its entire length. (2) The Keiva II moraine is sloping along its length from c. 100 m a.s.l. in the west (Varzuga River) to c. 250 m a.s.l. in the east (Ponoy River). (3) The KIZ was partly overrun and fragmented by erosive White Sea-based ice after formation. From these observations we conclude that the KIZ is not a synchronous feature formed along the lateral side of a White Sea-based ice lobe. If it was, the moraines should have a reversed slope. Rather, we interpret it to be time transgressive, formed at a northeastward-migrating junction between a warm-based Fennoscandian Ice Sheet expanding from the west and southwest into the White Sea depression, and a sluggish cold-based ice mass centred over eastern Kola Peninsula. In contrast to earlier reconstructions, we find it unlikely that an ice expansion of this magnitude was a mere re-advance during the deglaciation. Instead, we propose that the KIZ was formed during a major expansion of a Fennoscandian Ice Sheet at a time pre-dating the Last Glacial Maximum. [ABSTRACT FROM AUTHOR]

Abstract: The existence of sliding and frozen bed areas under ice sheets is significant in understanding basal thermal regimes, patterns of erosion and landform development, and in constraining boundary conditions for the reconstructions of ice sheets. Recognition of subglacial boundaries between sliding and frozen-bed areas for former ice sheets is typically based on distinct morphological contrasts between areas with glacial landform assemblages and relict areas showing little alteration of pre-existing features. Some of these boundaries, especially on continental shield areas, however, are clearly visible from air photos but have minimal topographic expression. Understanding the chronology and erosional development of such boundaries is important to provide insight into the pattern and persistence of basal conditions under ice sheets. Geomorphic evidence and cosmogenic radionuclide concentrations of bedrock outcrops on either side of two sliding boundaries on Ultevis and Arvestuottar, low-relief upland plateaus in northern Sweden, are consistent with negligible erosion in relict landscape (frozen bed) areas due to the last glaciation, but also indicate insignificant erosion in the sliding areas. Such a pattern and magnitude of landscape modification indicates that sliding was short lived in these areas, likely as a transient phase during deglaciation. These sites demonstrate that short periods of sliding are in some cases sufficient to produce landscapes that are recognized as ‘glacial’ from air photos. Thus, regions of sliding identified on shield areas must be viewed as the cumulative total area that has experienced sliding at any time during a glaciation. The actual extent of sliding areas during any single ice sheet phase is presumably considerably less than this cumulative total, which has important implications for establishing appropriate basal boundary conditions for ice sheet reconstructions. [Copyright &y& Elsevier]

Some areas within ice sheet boundaries retain pre-existing landforms and thus either remained as ice free islands (nunataks) during glaciation, or were preserved under ice. Differentiating between these alternatives has significant implications for paleoenvironment, ice sheet surface elevation, and ice volume reconstructions. In the northern Swedish mountains, in situ cosmogenic 10Be and 26Al concentrations from glacial erratics on relict surfaces as well as glacially eroded bedrock adjacent to these surfaces, provide consistent last deglaciation exposure ages (∼8–13 kyr), confirming ice sheet overriding as opposed to ice free conditions. However, these ages contrast with exposure ages of 34–61 kyr on bedrock surfaces in these same relict areas, demonstrating that relict areas were preserved with little erosion through multiple glacial cycles. Based on the difference in radioactive decay between 26Al and 10Be, the measured nuclide concentration in one of these bedrock surfaces suggests that it remained largely unmodified for a minimum period of 845−418+461 kyr. These results indicate that relict areas need to be accounted for as frozen bed patches in basal boundary conditions for ice sheet models, and in landscape development models. Subglacial preservation also implies that source areas for glacial sediments in ocean cores are considerably smaller than the total area covered by ice sheets. These relict areas also have significance as potential long-term subglacial biologic refugia. [Copyright &y& Elsevier]

The presence of well-developed tors, boulder fields, and weathering mantles in the Parkajoki area of northeastern Sweden, near the centre of Fennoscandian glaciation, has been used to suggest that these landscapes were preserved during all glacial cycles since ice-sheet initiation in the late Cainozoic. This implies that all successive large-scale glaciations must have had frozen bed conditions across this area to allow for subglacial landscape preservation. Cosmogenic 10Be and 26Al data from three tors and a meltwater channel in the Parkajoki area were used to test this hypothesis of landscape preservation through multiple glacial cycles. Apparent exposure ages of tor summit bedrock surfaces ranging between 79 and 37 ka in an area deglaciated at ∼11 ka are consistent with the interpretation of these features as relict landforms that have survived glaciation with little or no erosion. Single nuclide minimum exposure age data require that the tors have survived at least two complete glacial cycles. This estimate is based on (i) the approximate duration of periods of ice sheet cover versus ice free conditions as deduced from the DSDP 607 marine benthic foraminifer oxygen isotope record, in conjunction with (ii) a record of Fennoscandian ice sheet flow traces (and hence, ice sheet extent), and (iii) noting that cosmogenic nuclides are accumulated only during ice free periods. In addition, mean cosmogenic 26Al/10Be concentration ratios from two of the sites indicate a minimum model total history of 605 ka and a maximum erosion rate of 1.6 m Ma−1. Thus, the numerical ages confirm the overall qualitative interpretation of landscape preservation through multiple glacial cycles. [Copyright &y& Elsevier]

The Parkajoki area in northeastern Sweden is situated near the central area of Fennoscandian glaciation. Despite its location, the area is dominated by landforms induced by subaerial weathering and erosion processes, such as well-developed tors and associated saprolites, boulder fields, and boulder depressions. The glacial geomorphology is dominated by lateral and proglacial meltwater channels. Subglacial imprints indicative of thawed-bed conditions and reshaping by glacier sliding (e.g., fluting, drumlins, striae) are lacking. Hence, most of the landscape still exhibits a preglacial appearance. Because of its location, near the central area of glaciation, we attribute preservation to frozen-bed conditions of overriding ice sheets. The widespread distribution of well-developed tors and boulder fields, and the degree of chemical weathering of the bedrock, indicate that the area has been protected from glacial erosion during all glacial cycles since ice sheet initiation in the late Cenozoic. Unlike most other areas with tors in glaciated regions, the Parkajoki area is uniquely situated in the lowlands at 150–400 m.a.s.l. Moreover, this relict landscape is surrounded by glacial landscapes (including drumlins, ribbed moraine, and eskers) of varying age and at similar elevation. Hence, topographical reasons for this area being persistently cold-based cannot be invoked. By inference, we conclude that strain heat release never managed to cancel the initial subglacial permafrost conditions. We attribute this to divergent ice flow towards the convex outline of the ice sheet margin during deglaciations and to the relative roughness of the area compared to its surroundings. The implication is that to explain preservation throughout the Quaternary, all large-scale glaciations must have had a similar evolution concerning ice sheet configuration and internal dynamics as the last glaciation for which good constraints on evolution and outline are available. [Copyright &y& Elsevier]

Abstract: The Tibetan Plateau is a topographic feature of extraordinary dimension and has an important impact on regional and global climate. However, the glacial history of the Tibetan Plateau is more poorly constrained than that of most other formerly glaciated regions such as in North America and Eurasia. On the basis of some field evidence it has been hypothesized that the Tibetan Plateau was covered by an ice sheet during the Last Glacial Maximum (LGM). Abundant field- and chronological evidence for a predominance of local valley glaciation during the past 300,000 calendar years (that is, 300 ka), coupled to an absence of glacial landforms and sediments in extensive areas of the plateau, now refute this concept. This, furthermore, calls into question previous ice sheet modeling attempts which generally arrive at ice volumes considerably larger than allowed for by field evidence. Surprisingly, the robustness of such numerical ice sheet model results has not been widely queried, despite potentially important climate ramifications. We simulated the growth and decay of ice on the Tibetan Plateau during the last 125 ka in response to a large ensemble of climate forcings (90 members) derived from Global Circulation Models (GCMs), using a similar 3D thermomechanical ice sheet model as employed in previous studies. The numerical results include as extreme end members as an ice-free Tibetan Plateau and a plateau-scale ice sheet comparable, in volume, to the contemporary Greenland ice sheet. We further demonstrate that numerical simulations that acceptably conform to published reconstructions of Quaternary ice extent on the Tibetan Plateau cannot be achieved with the employed stand-alone ice sheet model when merely forced by paleoclimates derived from currently available GCMs. Progress is, however, expected if future investigations employ ice sheet models with higher resolution, bidirectional ice sheet-atmosphere feedbacks, improved treatment of the surface mass balance, and regional climate data and climate reconstructions. [Copyright &y& Elsevier]

The locations, shapes and masses of evolving ice sheets during interglacial to glacial transitions remain elusive, but need to be clarified for identification of the climatic feedbacks that amplified low-amplitude insolation changes to become the major global climatic shift of the last glacial maximum. The best way to explore probable spatial patterns of ice sheet inception and build-up is through ice sheet modelling, and validating model output against geomorphological and geological data. We have modelled the 120–55 kyr BP evolution of the Laurentide Ice Sheet, focusing on the locations, outline, flow patterns, and basal thermal regimes of its precursor ice centres. For this, we employed a time-dependent mass-balance-driven finite element model, forced with GRIP ice core oxygen isotope values calibrated to temperature. The model achieves a good fit to relict ice flow systems in Keewatin, the Interior Plains, and the Hudson Bay lowland. Thus, our model indicates that the first coherent ice centre in North America was situated over Baffin Island and the Melville and Boothia Peninsulas, and expanded onto mainland Canada at 90–75 kyr BP. A previously unknown 200-km-long moraine zone near the Thlewiaza River, Keewatin, may demarcate the marginal position of this early Central Arctic Ice Sheet. Cordilleran and Keewatin ice merged, and the Hudson Bay lowland was inundated from the east (from Quebec), by 65 kyr BP. [Copyright &y& Elsevier]

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]

The Eurasian ice-sheet complex (EISC) was the third largest ice mass during the Last Glacial Maximum (LGM), after the Antarctic and North American ice sheets. Despite its global significance, a comprehensive account of its evolution from independent nucleation centres to its maximum extent is conspicuously lacking. Here, a first-order, thermomechanical model, robustly constrained by empirical evidence, is used to investigate the dynamics of the EISC throughout its build-up to its maximum configuration. The ice flow model is coupled to a reference climate and applied at 10 km spatial resolution across a domain that includes the three main spreading centres of the Celtic, Fennoscandian and Barents Sea ice sheets. The model is forced with the NGRIP palaeo-isotope curve from 37 ka BP onwards and model skill is assessed against collated flowsets, marginal moraines, exposure ages and relative sea-level history. The evolution of the EISC to its LGM configuration was complex and asynchronous; the western, maritime margins of the Fennoscandian and Celtic ice sheets responded rapidly and advanced across their continental shelves by 29 ka BP, yet the maximum aerial extent (5.48 × 10 6 km 2 ) and volume (7.18 × 10 6 km 3 ) of the ice complex was attained some 6 ka later at c. 22.7 ka BP. This maximum stand was short-lived as the North Sea and Atlantic margins were already in retreat whilst eastern margins were still advancing up until c. 20 ka BP. High rates of basal erosion are modelled beneath ice streams and outlet glaciers draining the Celtic and Fennoscandian ice sheets with extensive preservation elsewhere due to frozen subglacial conditions, including much of the Barents and Kara seas. Here, and elsewhere across the Norwegian shelf and North Sea, high pressure subglacial conditions would have promoted localised gas hydrate formation. [ABSTRACT FROM AUTHOR]

The advance of the Cordilleran Ice Sheet (CIS) towards its Last Glacial Maximum (LGM) configuration and its subsequent retreat remain poorly understood. We use the glacial landform record to determine ice dynamics for the central sector of the CIS in northern British Columbia, Canada, beneath the LGM ice divide. We classify seventy ice-flow indicator flowsets based on morphology, elevation, orientation and cross-cutting relationships into one of three stages, whereby stage 1 is oldest and stage 3 youngest. Combined with ice-contact geomorphology, our reconstruction highlights complex changes in ice flow over time as a result of ice divide migrations through the LGM and deglacial phases. The orientation and distribution of landforms indicates active post-LGM ice retreat westward through the Cassiar and Omineca mountains. We map the regional distribution of independent mountain glaciers, ice caps, and ice fields that regrew during a cooling event in the Late Glacial and show that some of these readvance glaciers were subsequently overrun by advancing outlet glaciers of the CIS. We use the cross-cutting relationship between readvance glaciers and CIS outlet glaciers and available chronological data to reconstruct the eastern CIS margin during the Late Glacial for the first time. • Glacial flowsets document Cordilleran Ice Sheet advance and retreat. • Ice contact landforms reveal ice retreat towards the Skeena and Coast mountains. • Cross-cutting relationships suggest the CIS advanced during the Late Glacial. • We define the eastern margin of the Late Glacial CIS for the first time. [ABSTRACT FROM AUTHOR]

Abstract: The last glacial maximum (LGM) outline and subsequent retreat pattern (21–7kyr) of North American ice sheets are reasonably well established. However, the evolution of the ice sheets during their build-up phase towards the LGM between 115 and 21kyr has remained elusive, making it difficult to verify numerical ice sheet models for this important time interval. In this paper we outline the pre-LGM ice sheet evolution of the Laurentide and Cordilleran ice sheets by using glacial geological and geomorphological records to make a first-order reconstruction of ice sheet extent and flow pattern. We mapped the entire area covered by the Laurentide and Cordilleran ice sheets in Landsat MSS images and approximately 40% of this area in higher resolution Landsat ETM+ images. Mapping in aerial photographs added further detail primarily in Quebec-Labrador, the Cordilleran region, and on Baffin Island. Our analysis includes the recognition of approximately 500 relative-age relationships from crosscutting lineations. Together with previously published striae and till fabric data, these are used as the basis for relative-age assignments of regional flow patterns. For the reconstruction of the most probable ice sheet evolution sequence we employ a stepwise inversion scheme with a clearly defined strategy for delineating coherent landforms swarms (reflecting flow direction and configuration), and linking these to previously published constraints on relative and absolute chronology. Our results reveal that ice-dispersal centres in Keewatin and Quebec were dynamically independent for most of pre-LGM time and that a massive Quebec dispersal centre, rivalling the LGM in extent, existed at times when the SW sector of the ice sheet had not yet developed. The oldest flow system in eastern Quebec-Labrador (Atlantic swarm had an ice divide closer to the Labrador coast than later configurations). A northern Keewatin-Central Arctic Ice Sheet existed prior to the LGM, but is poorly chronologically constrained. There is also evidence for older and more easterly Cordilleran Ice Sheet divide locations than those that prevailed during the Late Wisconsinan. In terms of ice sheet build-up dynamics, it appears that “residual” ice caps after warming phases may have played an important role. In particular, the location and size of remnant ice masses at the end of major interstadials, i.e. OIS 5c and 5a, must have been critical for subsequent build-up patterns, because such remnant “uplands” may have fostered much more rapid ice sheet growth than what would have occurred on a fully deglaciated terrain. The ice-sheet configuration during stadials would also be governed largely by the additional topography that such “residual” ice constitutes because of inherent mass balance–topography feedbacks. [Copyright &y& Elsevier]

Studying the behaviour of the Antarctic Ice Sheet during periods with a climate that wassimilar or warmer than today can yield analogues for understanding its present-day changesand future trajectory. Among these, Marine Isotope Stage 11 (MIS 11) is one of thebest candidates given its limited insolation variability, which is very close to ourcurrent interglacial. However, Antarctic ice core data indicates that CO2 levels weremerely close to pre industrial, but that warmer-than-present temperatures (about 2 K)lasted for much longer than other Quaternary interglacials. Sea level was 6 to 13 mhigher than present, implying a potential substantial contribution from Antarctica. While substantial work has been conducted regarding the behaviour of the Greenland IceSheet during MIS 11, the Antarctic Ice Sheet history and dynamics during this period remainpoorly constrained, both in terms of its configuration and its potential contribution to sealevel. Thus, in this study we present an ensemble of continental scale transient simulations ofthe Antarctic Ice Sheet between 420 ka and 394 ka using the ice sheet model SICOPOLISforced by climate simulations from CCSM3. We evaluate the possible range of icesheet geometries, the Antarctic contribution to sea level rise and the associatedrange of uncertainties during the longest interglacial of the Quaternary period. Ourresults aim to provide insights into the style of the Antarctic deglaciation towards thepeak warming of MIS 11 and help set the bounds on the minimum volumes andextents of the western and eastern ice sheet counterparts as possible future analogues. [ABSTRACT FROM AUTHOR]