To provide a new reconstruction of the deglaciation of the Fennoscandian Ice Sheet, in the form of calendar-year time-slices, which are particularly useful for ice sheet modelling, we have compiled and synthesized published geomorphological data for eskers, ice-marginal formations, lineations, marginal meltwater channels, striae, ice-dammed lakes, and geochronological data from radiocarbon, varve, optically-stimulated luminescence, and cosmogenic nuclide dating. This is summarized as a deglaciation map of the Fennoscandian Ice Sheet with isochrons marking every 1000 years between 22 and 13 cal kyr BP and every hundred years between 11.6 and final ice decay after 9.7 cal kyr BP. Deglaciation patterns vary across the Fennoscandian Ice Sheet domain, reflecting differences in climatic and geomorphic settings as well as ice sheet basal thermal conditions and terrestrial versus marine margins. For example, the ice sheet margin in the high-precipitation coastal setting of the western sector responded sensitively to climatic variations leaving a detailed record of prominent moraines and other ice-marginal deposits in many fjords and coastal valleys. Retreat rates across the southern sector differed between slow retreat of the terrestrial margin in western and southern Sweden and rapid retreat of the calving ice margin in the Baltic Basin. Our reconstruction is consistent with much of the published research. However, the synthesis of a large amount of existing and new data support refined reconstructions in some areas. For example, the LGM extent of the ice sheet in northwestern Russia was located far east and it occurred at a later time than the rest of the ice sheet, at around 17-15 cal kyr BP. We also propose a slightly different chronology of moraine formation over southern Sweden based on improved correlations of moraine segments using new LiDAR data and tying the timing of moraine formation to Greenland ice core cold stages. Retreat rates vary by as much as an order of magnitude in different sectors of the ice sheet, with the lowest rates on the high-elevation and maritime Norwegian margin. Retreat rates compared to the climatic information provided by the Greenland ice core record show a general correspondence between retreat rate and climatic forcing, although a close match between retreat rate and climate is unlikely because of other controls, such as topography and marine versus terrestrial margins. Overall, the time slice reconstructions of Fennoscandian Ice Sheet deglaciation from 22 to 9.7 cal kyr BP provide an important dataset for understanding the contexts that underpin spatial and temporal patterns in retreat of the Fennoscandian Ice Sheet, and are an important resource for testing and refining ice sheet models.

In this article, we present a map of the glacial geomorphology of the Altai andWestern Sayan Mountains, covering an area of almost 600,000 km2. Although numerous studies provide evidence for restricted Pleistocene glaciations in this area, others have hypothesized the past existence of an extensive ice sheet. To provide a framework for accurate glacial reconstructions of the Altai and Western Sayan Mountains, we present a map at a scale of 1:1,000,000 based on a mapping from 30 m resolution ASTER DEM and 15 m/30 mresolution Landsat ETM+ satellite imagery. Four landform classes have been mapped: marginal moraines, glacial lineations, hummocky terrain, and glacial valleys. Our mapping reveals an abundance of glacial erosional and depositional landforms. The distribution of these glacial landforms indicates that the Altai and Western Sayan Mountains have experienced predominantly alpine-style glaciations, with some small ice caps centred on the higher mountain peaks. Large marginal moraine complexes mark glacial advances in intermontane basins. By tracing the outer limits of present-day glaciers, glacial valleys, and moraines, we estimate that the past glacier coverage have totalled to 65,000 km2 (10.9% of the mapped area), whereas present-day glacier coverage totals only 1300 km2 (0.2% of the mapped area). This demonstrates the usefulness of remote sensing techniques for mapping the glacial geomorphology in remote mountain areas and for quantifying the past glacier dimensions. The glacial geomorphological map presented here will be used for further detailed reconstructions of the paleoglaciology and paleoclimate of the region.

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.

Paleoclimatic constraints from regions at the confluence of major climate systems are particularly important in understanding past climate change. Using geomorphic mapping based on remote sensing and field investigations, combined with in situ cosmogenic Be-10 and Al-26 dating of boulders associated with glacial landforms, we investigate the chronology of past glaciation in the Inylchek and Sary-Dzaz valleys in the eastern Kyrgyz Tian Shan, a tectonically active area with some of the highest peaks in the world outside of the Himalayas. Cosmogenic Be-10 and (26) Al exposure ages of boulders on moraines record up to five glacial advances including: Lateglacial age lateral moraine remnants and meltwater channels in the upper Inylchek Valley; Last Glacial Maximum (LGM, Marine Oxygen Isotope Stage [MIS] 2) moraines in the Sary-Dzaz Valley and in a terminal moraine complex at the west end of the Inylchek Valley, overriding older moraines; an MIS 4 or 5 moraine remnant above the Inylchek terminal moraine complex; and an older high moraine remnant down-valley from the confluence of the Inylchek and Sary-Dzaz valleys. The evidence for glacial extent in this study is consistent with a limited ice expansion hypothesis for Tian Shan glaciation. Published results from the western and central Kyrgyz Tian Shan do not show evidence for significant LGM glacier expansion, which in combination with the results presented here, indicate a spatial variation in glacier records along the Tian Shan. This may reflect either paleoclimatic gradients or the impact of local physiographic conditions on responses to regional climate change, or both.

The century-long debate over the origins of inner gorges that were repeatedly covered by Quaternary glaciers hinges upon whether the gorges are fluvial forms eroded by subaerial rivers, or subglacial forms cut beneath ice. Here we apply cosmogenic nuclide exposure dating to seven inner gorges along similar to 500 km of the former Fennoscandian ice sheet margin in combination with a new deglaciation map. We show that the timing of exposure matches the advent of ice-free conditions, strongly suggesting that gorges were cut by channelized subglacial meltwater while simultaneously being shielded from cosmic rays by overlying ice. Given the exceptional hydraulic efficiency required for meltwater channels to erode bedrock and evacuate debris, we deduce that inner gorges are the product of ice sheets undergoing intense surface melting. The lack of postglacial river erosion in our seven gorges implicates subglacial meltwater as a key driver of valley deepening on the Baltic Shield over multiple glacial cycles.

The Tibetan Plateau holds an ample record of past glaciations, and there is an extensive set of glacial deposits dated by exposure dating. Here a compilation is presented of ¹⁰Be exposure ages from 485 glacial deposits with 1855 individual samples on the Tibetan Plateau, and ELA depression estimates for the glacial deposits based on a simple toe to headwall ratio approach. To recalculate the Tibetan Plateau exposure ages, ¹⁰Be production rates from 24 calibration sites across the world are compiled and recalibrated yielding an updated global reference ¹⁰Be production rate. The recalculated exposure ages from the Tibetan Plateau glacial deposits are then divided into three groups based on exposure age clustering, to discriminate good (well-clustered) from poor (scattered) deglaciation ages. A major part of the glacial deposits have exposure ages affected by prior or incomplete exposure, complicating exposure age interpretations. The well-clustered deglaciation ages are primarily from mountain ranges along the margins of the Tibetan Plateau with a main peak between 10 and 30 ka, indicating glacial advances during the global LGM. A large number of deglaciation ages older than 30 ka indicates maximum glaciation predating the LGM, but the exposure age scatter generally prohibits accurate definition of the glacial chronology. The ELA depression estimates scatter significantly, but the main part is remarkably low. Average ELA depressions of 337 ± 197 m for the LGM and 494 ± 280 m for the pre-LGM indicate restricted glacier expansion.

The Cordilleran Ice Sheet (CIS) covered most of British Columbia and southern Yukon Territory at the local Last Glacial Maximum (ILGM) 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 Be-10 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 Be-10 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 ILGM age of ca 17 ka for the Puget lobe of the CIS in Washington State; they are 1.7 ka younger than the ILGM 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 C-14 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 ILGM 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 ILGM. 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.

Glacial erosion has created distinctive types of landscapes reflecting the extent, duration, and processes of the parent glaciers. Alpine landscapes are representative of pervasive erosion at the local scale. Landscapes formerly covered by larger scale glaciations display a wide range of appearances, from intensively eroded to preserved. Landscapes of selective linear erosion were formed, where subglacial melting occurred along certain corridors that were flanking regions of subglacial freezing. Landscapes of areal scouring were formed where subglacial melting on low-relief surfaces allowed spatially extensive subglacial stripping to dominate. Landscapes of little or no erosion indicate a dominance of subglacial freezing conditions.

Glacial geomorphological mapping from satellite imagery and field investigations provide the basis for a reconstructionof the extent and style of glaciation of the Shaluli Shan, a mountainous area on the southeastern TibetanPlateau. Our studies provide evidence for multiple glaciations, including the formation of regional ice caps andvalley glaciers. The low-relief topographywithin the Shaluli Shan, the Haizishan Plateau, and Xinlong Plateau displayzonal distributions of glacial landforms that is similar to those imprinted by Northern Hemisphere ice sheetsduring the last glacial cycle, indicating the presence of regional, polythermal ice caps. Abundant alpine glaciallandforms occur on high mountain ranges. The pattern of glaciated valleys centered on high mountain rangesand ice-scoured low relief granite plateaus with distinctive patterns of glacial lineations indicate a strong topographiccontrol on erosional and depositional patterns by glaciers and ice caps. In contrast to the Shaluli Shan,areas farther north and west on the Tibetan Plateau have not yielded similar landform evidence for regionalice capswith complex thermal basal conditions. Such spatial differences across the Tibetan Plateau are the resultof variations in climate and topography that control the extent and style of glaciations and that reinforce the importanceof detailed geomorphological mapping for understanding paleoclimate variations and characteristics offormer glaciations.

The glacial geomorphology of the Tian Shan has been mapped, with the study area covering almost 638,000km(2). The map, designed to be printed at A0 size due to the elongated shape of the mountain range, is presented at a scale of 1:1,100,000. Five glacial landform categories are presented; glacial valleys, marginal moraines, glacial lineations, hummocky terrain and meltwater channels. These landform categories were mapped using the Shuttle Radar Topography Mission (SRTM) digital elevation model (90m resolution), Landsat 7 ETM+ satellite imagery (30m resolution), and images contained in Google Earth. The mapped landforms were created by glaciers that were restricted to mountain areas and their immediate surroundings.