Glacial advance,occupation and retreat sediments associated with multi‐stage ice‐dammed lakes: north‐central Alberta,Canada

Ice sheets that advance upvalley, against the regional gradient, commonly block drainage and result in ice‐dammed proglacial lakes along their margins during advance and retreat phases. Ice‐dammed glacial lakes described in regional depositional models, in which ice blocks a major lake outlet, are often confined to basins in which the glacial lake palaeogeographical position generally remains semi‐stable (e.g. Great Lakes basins). However, in places where ice retreats downvalley, blocking regional drainage, the palaeogeographical position and lake level of glacial lakes evolve temporally in response to the position of the ice margin (referred to here as ‘multi‐stage’ lakes). In order to understand the sedimentary record of multi‐stage lakes, sediments were examined in 14 cored boreholes in the Peace and Wabasca valleys in north‐central Alberta, Canada. Three facies associations (FAI–III) were identified from core, and record Middle Wisconsinan ice‐distal to ice‐proximal glaciolacustrine (FAI) sediments deposited during ice advance, Late Wisconsinan subglacial and ice‐marginal sediments (FAII) deposited during ice‐occupation, and glaciolacustrine sediments (FAIII) that record ice retreat from the study area. Modelling of the lateral extent of FAs using water wells and gamma‐ray logs, combined with interpreted outlets and mapped moraines based on LiDAR imagery, facilitated palaeogeographical reconstruction of lakes and the identification of four major retreat‐phase lake stages. These lake reconstructions, together with the vertical succession of FAs, are used to develop a depositional model for ice‐dammed lakes during a cycle of glacial advance and retreat. This depositional model may be applied in other areas where meltwater was impounded by glacial ice advancing up the regional gradient, in order to understand the complex interaction between depositional processes, ice‐marginal position, and supply of meltwater and sediment in the lake basin. In particular, this model could be applied to decipher the genetic origin of diamicts previously interpreted to record strictly subglacial deposition or multiple re‐advances.     

The formation of natural cryogenic brines

The source of salts in the Ca-chloridic, hypersaline brines (up to 190 g Cl L⁻¹) occurring in crystalline basement rocks in the Canadian, Fennoscandian and Bohemian Shields and their evolution have been investigated and reported. The Cl-Br-Na relationship indicates that these waters have been concentrated from seawater, by freezing during glacial times. The Na/Cl ratio (0.25 to 0.35) in the more saline fluids is compatible with cooling down to −30°C, where the most saline waters have been concentrated by a factor of 25 to 30 relative to the parent seawater.The brines formed from seawater within cryogenic troughs, along the subarctic continental margins, around ice sheets. The depressions within which the brines formed are the cryogenic analogues of the classic, evaporitic lagoon. One million years suffice to saturate with brine a 2000km-radius by 1km-depth rock volume at an H₂O removal rate of only 2.8 mm/yr. Density-induced brine migration on a continental scale takes place via fissures below the ice.Our calculations, that were performed on a hypothetical ice sheet with dimensions compatible with the Laurentide ice sheet, demonstrate that during 1m.y., a 60m thick cryogenic sediment section could have formed. However, the precipitated minerals (mirabilite and hydrohalite) are repeatedly dispersed by the advance and retreat of the ice sheet, dissolved by melt water-seawater mixtures, and eroded during postglacial uplift, leaving almost no trace in the geological record.The cryogenic brines formed intermittently during and between glacial periods. The repeating advance and retreat of the ice sheets exerted a major control on the direction and intensity of brine flow. The cryogenic concentration of seawater and the migration of brine towards the center of the glaciostatic depression occurred mainly during the build up of the ice sheet, while reversal of the water flow from the center of the cryogenic basin outwards happened upon deglaciation. The flow of the waters in the subsurface was, inevitably, accompanied by significant dilution with melt water from the ice sheets.Using a “granitic” U concentration of 4 ppm and a (Ca-Mg mass balance based) rock/water ratio anywhere between 3.4 and 6.8 kg L⁻¹, a few hundred thousand years of brine-rock interaction are sufficient for the growth of ¹²⁹I in the most saline Canadian Shield brine to its present concentration (3.4×10⁸ atoms ¹²⁹I L⁻¹). Hence, both the formation of the saline fluids and their emplacement in their present sites occurred most likely within the Pleistocene.The young age calculated for cryogenic brines in crystalline shields and the dynamic water flow therein should raise concern about the planning and construction of high-grade nuclear waste repositories in such rocks, which are already under way.     

On the reconstruction of pleistocene ice sheets: A review

Pleistocene ice sheets can be reconstructed through three separate approaches: (1) Evidence based on glacial geological studies, such as erratic trains, till composition, crossing striations and exposures of multiple tills/nonglacial sediments. (2) Reconstructions based on glaciological theory and observations. These can be either two- or three-dimensional models; they can be constrained by ‘known’ ice margins at specific times; or they can be ‘open-ended’ with the history of growth and retreat controlled by parameters resting entirely within the model. (3) Glacial isostatic rebound after deglaciation provides a measure of the distribution of mass (ice) across a region. A ‘best fit’ ice sheet model can be developed that closely approximates a series of relative sea level curves within an area of a former ice sheet; in addition, the model should also provide a reasonable sea level fit to relative sea level curves at sites well removed from glaciation.This paper reviews some of the results of a variety of ice sheet reconstructions and concentrates on the various attempts to reconstruct the ice sheets of the last (Wisconsin, Weischelian, Würm, Devensian) glaciation. Evidence from glacial geology suggests flow patterns at variance with simple, single-domed ice sheets over North America and Europe. In addition, reconstruction of ice sheets from glacial isostatic sea level data suggests that the ice sheets were significantly thinner than estimates based on 18 ka equilibrium ice sheets (cf. Denton and Hughes, 1981). The review indicates it is important to differentiate between ice divides, which control the directions of glacial flow, and areas of maximum ice thickness, which control the glacial isostatic rebound of the crust upon deglaciation. Recent studies from the Laurentide Ice Sheet region indicate that the center of mass was not over Hudson Bay; that a major ice divide lay east of Hudson Bay so that flow across the Hudson Bay and James Bay lowlands was from the northeast; that Hudson Bay was probably open to marine invasions two or three times during the Wisconsin Glaciation; and that the Laurentide Ice Sheet was thinner than an equilibrium reconstruction would suggest.     

Late Wisconsinan glacial landsystems on Atlantic Canadian shelves: New evidence from multibeam and single-beam sonar data

Multibeam sonar surveys in the past decade, augmented by single-beam data from the OLEX charting system, reveal landsystems on Atlantic Canadian shelves that are diagnostic of Late Wisconsinan ice-sheet dynamics. Four landsystems are described. (1) The Bay of Fundy landsystem comprises two contrasting sets of bedforms, and is interpreted as evidence of topographically controlled fast-flowing ice adjacent to slower-moving ice. (2) The German Bank landsystem off southwest Nova Scotia is comprised of glacially fluted terrain overprinted by De Geer moraines and arcuate recessional moraines. We infer that a flow of grounded glacial ice out of the Bay of Fundy was followed by steady retreat, punctuated by at least one major re-advance. (3) The Placentia Bay landsystem consists of a convergent field of streamlined landforms with superimposed De Geer moraines, overprinted in one area by flutings. We infer that this landsystem was formed in the onset zone of fast-flowing ice, and that overprinting was due to a re-advance of ice from offshore. (4) The south coast of Newfoundland landsystem, which includes arcuate, fjord-mouth moraines and a coast-parallel, fluted moraine, indicates strong topographic control on a retreating marine ice margin as it reached a fjord coastline. These submarine glacial landsystems are not inconsistent with a conceptual model showing Late Wisconsinan ice advance to shelf edges, rapid calving retreat along deepwater channels and slower retreat of ice margins grounded in shallow water. The re-advances documented two of the study areas have parallels in the Last British Ice Sheet, confirming that the reorganization of marine-based ice sheets, caused by calving in embayments, led to internally forced re-advances.     

Sedimentation and Iithostratigraphy of the North Sea Drift and Lowestoft Till Formations in the Coastal cliffs of northeast Norfolk,England

The most complete terrestrial sequence of Anglian (Elsterian) glacial sediments in western Europe was investigated in northeast Norfolk, England in order to reconstruct the evolution of the contemporary palaeoenvironments. Lithostratigraphically the glacial sediments in the northeast Norfolk coastal cliffs can be divided into the Northn Sea Drift and Lowestoft Till Formations. Three of the diamicton members of the North Sea Drift Formation (Happisburgh, Walcott and Cromer Diamictons) were deposited as lodgement and/or subglacial deformation till by grounded ice, but one, the Mundesley Diamicton, is waterlain and was deposited in an extensive glacial lake. Sands and fine sediments interbedded between the diamictons represent deltaic sands and glaciolacustrine sediments derived not solely from the melting ice in the north but also from extra-marginal rivers in the south. The Lowestoft Till Formation is not well preserved in the cliffs but includes lodgement till (Marly Drift till) and, most probably, associated meltwater deposits. Extensive glaciotectonism in the northern part of the area is shown to relate to oscillating ice that deposited the Cromer Diamicton and also partially to the ice sheet that deposited the Marly Drift till. It is suggested that during the Anglian Stage the present day northeast Norfolk coast was situated on the northwestern margin of an extensive glaciolacustrine basin. This basin was dammed by the Scandinavian ice sheet in the north and northeast. Because the grounding line of this ice sheet oscillated in space and time, part of the North Sea Drift diamictons were deposited directly by this ice. However, during ice retreat phases glaciolacustrine deposition comprised waterlain diamicton, sands and fines. When the Scandinavian ice sheet was situated in northernmost Norfolk, the British ice sheet (responsible for depositing the Marly Drift facies) entered the area from the west. This ice sheet partially deformed the North Sea Drift Formation sediments in the northern part of the area but not in the south, where the British ice sheet apparently terminated in water. The interplay of these two ice sheets on the northern and western margins of the glacial lake is thought to be the major determining factor for the accumulation of thick glacial deposits in this area during the Anglian glaciation.     

Radiocarbon Age Constraints on Rates of Advance and Retreat of the Puget Lobe of the Cordilleran Ice Sheet during the Last Glaciation

Calibrated radiocarbon dates of organic matter below and above till of the last (Fraser) glaciation provide limiting ages that constrain the chronology and duration of the last advance–retreat cycle of the Puget Lobe in the central and southeastern Puget Lowland. Seven dates for wood near the top of a thick proglacial delta have a weighted mean age of 17,420 ± 90 cal yr B.P., which is the closest limiting age for arrival of the glacier near the latitude of Seattle. A time–distance curve constructed along a flowline extending south from southwestern British Columbia to the central Puget Lowland implies an average glacier advance rate of ca. 135 m/yr. The glacier terminus reached its southernmost limit ca. 16,950 yr ago and likely remained there for ca. 100 yr. In the vicinity of Seattle, where the glacier reached a maximum thickness of 1000 m, ice covered the landscape for ca. 1020 yr. Postglacial dates constraining the timing of ice retreat in the central lowland are as old as 16,420 cal yr B.P. and show that the terminus had retreated to the northern limit of the lowland within three to four centuries after the glacial maximum. The average rate of retreat was about twice the rate of advance and was enhanced by rapid calving recession along flowline sectors where the glacier front crossed deep proglacial lakes.     

Deglaciation chronology inferred from marine sediments in a proglacial lake basin, western Spitsbergen, Svalbard

Sub-bottom sediment profiles and sediment cores show that the lacustrine sediments in lake Linnevatnet are underlain by marine sediments and a basal till that mantles the bedrock. The till was probably deposited by the glacier that during the Late Weichselian glacial maximum removed all pre-existing sediments from the basin. The cores were collected in closed basins, where continuous deposition is expected. The marine sediment in the studied cores is up to 8 m thick and consists of bioturbated clay and silt. Radiocarbon dates on shells from the base of the marine sequence suggest that glacial retreat from the lake basin occurred around 12,500BP. This is more than a thousand years older than basal shell dates from raised marine sediments on the slopes above the lake. Typical ice proximal litbofacies were not identified in the cores. stratigraphic record indicates both a rapid glacial retreat and that no younger glacial re-advances occurred. During the Younger Dryas local glaciers on western Svalbard were smaller than during the Little Ice Age. This is in sharp contrast to western Europe, where Younger Dryas glaciers were much larger than those the Little Ice Age.     

LakeCC: a tool for efficiently identifying lake basins with application to palaeogeographic reconstructions of North America

Along the margins of continental ice sheets, lakes formed in isostatically depressed basins during glacial retreat. Their shorelines and extent are sensitive to the ice margin and the glacial history of the region. Proglacial lakes, in turn, also impact the glacial isostatic adjustment due to loading, and ice dynamics by posing a marine-like boundary condition at the ice margin. In this study we present a tool that efficiently identifies lake basins and the corresponding maximum water level for a given ice sheet and topography reconstruction. This algorithm, called the LakeCC model, iteratively checks the whole map for a set of increasing water levels and fills isolated basins until they overflow into the ocean. We apply it to the present-day Great Lakes and the results show good agreement (∼1−4%) with measured lake volume and depth. We then apply it to two topography reconstructions of North America between the Last Glacial Maximum and the present. The model successfully reconstructs glacial lakes such as Lake Agassiz, Lake McConnell and the predecessors of the Great Lakes. LakeCC can be used to judge the quality of ice sheet reconstructions. © 2019 The Authors Journal of Quaternary Science Published by John Wiley & Sons Ltd.     

Multiple re‐advances of a Lake Vättern outlet glacier during Fennoscandian Ice Sheet retreat,south‐central Sweden

Lake Vättern represents a critical region geographically and dynamically in the deglaciation of the Fennoscandian Ice Sheet. The outlet glacier that occupied the basin and its behaviour during ice‐sheet retreat were key to the development and drainage of the Baltic Ice Lake, dammed just west of the basin, yet its geometry, extent, thickness, margin dynamics, timing and sensitivity to regional retreat forcing are rather poorly known. The submerged sediment archives of Lake Vättern represent a missing component of the regional Swedish deglaciation history. Newly collected geophysical data, including high‐resolution multibeam bathymetry of the lake floor and seismic reflection profiles of southern Lake Vättern, are used here together with a unique 74‐m sediment record recently acquired by drill coring, and with onshore LiDAR‐based geomorphological analysis, to investigate the deglacial environments and dynamics in the basin and its terrestrial environs. Five stratigraphical units comprise a thick subglacial package attributed to the last glacial period (and probably earlier), and an overlying >120‐m deglacial sequence. Three distinct retreat–re‐advance episodes occurred in southern Lake Vättern between the initial deglaciation and the Younger Dryas. In the most recent of these, ice overrode proglacial lake sediments and re‐advanced from north of Visingsö to the southern reaches of the lake, where ice up to 400 m thick encroached on land in a lobate fashion, moulding crag‐and‐tail lineations and depositing till above earlier glacifluvial sediments. This event precedes the Younger Dryas, which our data reveal was probably restricted to north‐central sectors of the basin. These dynamics, and their position within the regional retreat chronology, indicate a highly active ice margin during deglaciation, with retreat rates on average 175 m a^?1. The pronounced topography of the Vättern basin and its deep proglacial‐dammed lake are likely to have encouraged the dynamic behaviour of this major Fennoscandian outlet glacier.     

Detrital zircon age and provenance constraints on late Paleozoic ice-sheet growth and dynamics in Western and Central Australia

U–Pb dating and Hf-isotope provenance analysis of detrital zircons from the glaciogenic lower Permian Grant Group of the Canning Basin indicate sources principally from basement terranes in central Australia, with subordinate components from terranes to the south and north. Integrating these data with field outcrop and subsurface evidence for ice sheets, including glacial valleys and striated pavements along the southern and northern margins of the basin, suggests that continental ice sheets extended over several Precambrian upland areas of western and central Australia during the late Paleozoic ice age (LPIA). The youngest zircons constrain the maximum age for contemporaneous ice sheet development to the late Carboniferous (Kasimovian), whereas palynology provides a minimum age of early Permian (Asselian–Sakmarian). Considering the palynological age of the Grant Group within the context of regional and global climate proxies, the main phase of continental ice sheet growth was possibly in the Ghzelian–Asselian. The presence of ice sheets older than Kasimovian in western and central Australia remains difficult to prove given a regional gap in deposition possibly covering the mid-Bashkirian to early Ghzelian within the main depocentres and even larger along basin margins, and the poor evidence for older Carboniferous glacial facies. There is also no evidence for extensive glacial facies younger than mid-Sakmarian in this region as opposed to eastern Australia where the youngest regional glacial phase was Guadalupian.     

Modeling North American Freshwater Runoff through the Last Glacial Cycle

The Northern Hemisphere ice sheets decayed rapidly during deglacial phases of the ice-age cycle, producing meltwater fluxes that may have been of sufficient magnitude to perturb oceanic circulation. The continental record of ice-sheet history is more obscured during the growth and advance of the last great ice sheets, ca. 120,000–20,000 yr B.P., but ice cores tell of high-amplitude, millennial-scale climate fluctuations that prevailed throughout this period. These climatic excursions would have provoked significant fluctuation of ice-sheet margins and runoff variability whenever ice sheets extended to mid-latitudes, giving a complex pattern of freshwater delivery to the oceans. A model of continental surface hydrology is coupled with an ice-dynamics model simulating the last glacial cycle in North America. Meltwater discharged from ice sheets is either channeled down continental drainage pathways or stored temporarily in large systems of proglacial lakes that border the retreating ice-sheet margin. The coupled treatment provides quantitative estimates of the spatial and temporal patterns of freshwater flux to the continental margins. Results imply an intensified surface hydrological environment when ice sheets are present, despite a net decrease in precipitation during glacial periods. Diminished continental evaporation and high levels of meltwater production combine to give mid-latitude runoff values that are highly variable through the glacial cycle, but are two to three times in excess of modern river fluxes; drainage to the North Atlantic via the St. Lawrence, Hudson, and Mississippi River catchments averages 0.356 Sv for the period 60,000–10,000 yr B.P., compared to 0.122 Sv for the past 10,000 yr. High-amplitude meltwater pulses to the Gulf of Mexico, North Atlantic, and North Pacific occur throughout the glacial period, with ice-sheet geometry controlling intricate patterns of freshwater routing variability. Runoff from North America is staged in the final deglaciation, with a stepped sequence of pulses through the Mississippi, St. Lawrence, Arctic, and Hudson Strait drainages.     

The evolution of the terrestrial-terminating Irish Sea glacier during the last glaciation

Here we reconstruct the last advance to maximum limits and retreat of the Irish Sea Glacier (ISG), the only land-terminating ice lobe of the western British Irish Ice Sheet. A series of reverse bedrock slopes rendered proglacial lakes endemic, forming time-transgressive moraine- and bedrock-dammed basins that evolved with ice marginal retreat. Combining, for the first time on glacial sediments, optically stimulated luminescence (OSL) bleaching profiles for cobbles with single grain and small aliquot OSL measurements on sands, has produced a coherent chronology from these heterogeneously bleached samples. This chronology constrains what is globally an early build-up of ice during late Marine Isotope Stage 3 and Greenland Stadial (GS) 5, with ice margins reaching south Lancashire by 30 ± 1.2 ka, followed by a 120-km advance at 28.3 ± 1.4 ka reaching its 26.5 ± 1.1 ka maximum extent during GS-3. Early retreat during GS-3 reflects piracy of ice sources shared with the Irish-Sea Ice Stream (ISIS), starving the ISG. With ISG retreat, an opportunistic readvance of Welsh ice during GS-2 rode over the ISG moraines occupying the space vacated, with ice margins oscillating within a substantial glacial over-deepening. Our geomorphological chronosequence shows a glacial system forced by climate but mediated by piracy of ice sources shared with the ISIS, changing flow regimes and fronting environments.     

Pattern and chronology of glacial Lake Peace shorelines and implications for isostacy and ice‐sheet configuration in northeastern British Columbia,Canada

Recognition of positions of glacial lakes along the margin of continental ice sheets is critical in reconstructing ice configuration during deglaciation. Advances in remote sensing technology (e.g. LiDAR) have enabled the generation of accurate digital‐elevation models (DEMs) that reveal unprecedented geomorphic detail. Combined with geographical information systems, these tools have considerably advanced the mapping and correlation of geomorphic features such as relict shorelines. Shorelines of glacial Lake Peace (GLP) developed between the Laurentide and Cordilleran ice sheets in northeastern British Columbia and northwestern Alberta. Shoreline mapping from high resolution DEMs produced more than 55 500 elevation data points from 3231 shorelines, enabling the identification of four major phases of GLP: Phase I (altitude 960–990 m a.s.l.); Phase II (890–915 m a.s.l.); Phase III (810–865 m a.s.l.); and Phase IV (724–733 m a.s.l.). The timing of Phase II of GLP is estimated by two optical ages of <16.0±2.5 and 14.2±0.5 ka BP. Extensive mapping of the shorelines allows for measuring of glacial isostatic adjustment as ice retreated. Shorelines currently dip to the northeast at around 0.4–0.5 m km^?1. This slope reflects the asynchronous retreat of the Cordilleran (CIS) and Laurentide (LIS) ice sheets. The relative uplift in the southwest of the study area within the Rocky Mountains and foothills suggests that the Late Wisconsinan (MIS 2) CIS persisted in the foothill after the LIS lost mass and retreated, or that the Late Wisconsinan CIS was very thick and caused deep crustal loading, which resulted in more uplift in the southwest before reaching equilibrium during, or shortly after deglaciation.     

Seafloor glacial features reveal the extent and decay of the last British Ice Sheet,east of Scotland

Three‐dimensional (3D) seismic datasets, 2D seismic reflection profiles and shallow cores provide insights into the geometry and composition of glacial features on the continental shelf, offshore eastern Scotland (58° N, 1–2° W). The relic features are related to the activity of the last British Ice Sheet (BIS) in the Outer Moray Firth. A landsystem assemblage consisting of four types of subglacial and ice marginal morphology is mapped at the seafloor. The assemblage comprises: (i) large seabed banks (interpreted as end moraines), coeval with the Bosies Bank moraine; (ii) morainic ridges (hummocky, push and end moraine) formed beneath, and at the margins of the ice sheet; (iii) an incised valley (a subglacial meltwater channel), recording meltwater drainage beneath former ice sheets; and (iv) elongate ridges and grooves (subglacial bedforms) overprinted by transverse ridges (grounding line moraines). The bedforms suggest that fast‐flowing grounded ice advanced eastward of the previously proposed terminus of the offshore Late Weichselian BIS, increasing the size and extent of the ice sheet beyond traditional limits. Complex moraine formation at the margins of less active ice characterised subsequent retreat, with periodic stillstands and readvances. Observations are consistent with interpretations of a dynamic and oscillating ice margin during BIS deglaciation, and with an extensive ice sheet in the North Sea basin at the Last Glacial Maximum. Final ice margin retreat was rapid, manifested in stagnant ice topography, which aided preservation of the landsystem record. Copyright © 2008 John Wiley & Sons, Ltd.     

Evidence from Ar/Ar Ages of Individual Hornblende Grains for Varying Laurentide Sources of Iceberg Discharges 22,000 to 10,500 yr B.P.

The abundance and lithic content of ice rafted detritus in glacial North Atlantic sediment cores vary abruptly on millennial time scales that have been correlated to Dansgaard-Oeschger cycles in the Greenland ice cores. There is growing evidence that various ice sheet outlets contributed increased iceberg fluxes at multiple discrete intervals, and the relative timing of iceberg discharges from different sources is important for understanding interactions between oceans and ice sheets. We present a provenance study based on ⁴⁰Ar/³⁹Ar dates of individual hornblende grains from 20 samples taken at 600 to 700 yr spacing between 10,500 and 22,000 yr B.P., from Orphan Knoll core EW9303-GGC31. Heinrich layers are characterized by a dominant Paleoproterozoic hornblende provenance consistent with published studies. A change in provenance between Heinrich events H2 and H1 indicates contributions of iceberg calving from the Newfoundland and southern Labrador margins. Between H1 and the Younger Dryas interval, Paleoproterozoic ice rafted grains remained dominant. The dominance of Baffin Island (or Greenland?) sources to the ice rafted detritus is ascribed to the retreat of the southern Laurentide ice sheet at about the time of H1—a retreat that isolated Newfoundland and southern Labrador ice from the shelf-slope boundary.     

Glacial advance and retreat sequences in a Permo-Carboniferous section, central Transantarctic Mountains

Episodes of glacial advance and retreat can be recognized through analysis of vertical facies sequences in the Permo-Carboniferous Pagoda Formation of the Beardmore Glacier area, Antarctica. The formation includes a remarkably complete record of continental sedimentation near the terminus of a temperate glacier. Facies sequence is pre-eminent for inferring glacial advance and retreat. Other important criteria are abundance and geometry of sandstone interbedded with diamictite, diamictite character and nature of bed contacts. Using these characteristics advance and retreat sequences 5–60 m thick are recognized. A sharp contact, with a striated surface and erosional relief, overlain by structureless diamictite (lodgement till) is typical of grounded ice advance. Grounded ice retreat is characterized by structureless diamictite (lodgement till), overlain by crudely stratified diamictite (melt-out till) and then by diamictite interbedded with sandstone and conglomerate (flow till and glacio-fluvial or glacio-lacustrine deposits). Gradational contacts between shale overlain by diamictite and diamictite overlain by shale characterize advance and retreat, respectively, in subaqueous settings. Pauses in sediment accumulation, minor(?) fluctuations of the ice margin, and/or changes in subglacial dynamics are indicated by specific features within diamictite units such as probable frost-wedge casts, single layer boulder beds, sharp sedimentary contacts and changes in diamictite character. These minor(?) events are superimposed upon the main advance-retreat cycles. Study of both the overall facies sequence and of individual diagnostic structures, albeit in an incomplete stratigraphic record, permits a distinction between major and minor advance-retreat events. As many as six major advance-retreat cycles exist in some Pagoda sections, but the number of cycles present varies in different sections.     

Drainage beneath ice sheets: groundwater–channel coupling,and the origin of esker systems from former ice sheets

The nature of the drainage system beneath ice sheets is crucial to their dynamic behaviour but remains problematic. An experimentally based theory of coupling between groundwater and major channel systems is applied to the esker systems in the area occupied the last ice sheet in Europe, which we regard as a fossil imprint of major longitudinal drainage channels. We conclude that the large-scale distribution and spacing of major eskers is consistent with the theory of groundwater control, in which esker spacing is partly controlled by the transmissivity of the bed.It is concluded that esker patterns reflect the large-scale organisation of the subglacial drainage pattern in which channel development is coupled to groundwater flow and to the ice sheet's dynamic regime. The theory is then used to deduce: basal meltwater recharge rates and their spatial variability from esker spacing in an area in which the ice sheet was actively streaming during its final retreat; patterns of palaeo-groundwater flow and head distribution; and the seasonally varying magnitude of discharge from stream tunnels at the retreating ice sheet margin. Major channel/esker systems appear to have been stable at least over several hundred of years during the retreat of the ice sheet, although major dynamic events are demonstrably associated with major shifts in the hydraulic regime.Modelling suggests: that glaciation can stimulate deep groundwater circulation cells that are spatially linked to channel locations, with groundwater flow predominantly transverse to ice flow; that the circulation pattern has the potential to create large-scale anomalies in groundwater chemistry; and that the spacing of channels will change through the glacial cycle, influencing water pressures in stream tunnels, subglacial hydraulic gradients and effective pressure. If the latter is reduced sufficiently, it could trigger enhanced bed deformation, thus coupling drainage to ice sheet movement. It suggests the possibility of distinctive phases of sediment deformation and drumlin mobilisation during a glacial cycle.     

Modeling the deglaciation of the Green Bay Lobe of the southern Laurentide Ice Sheet

We use a time-dependent two-dimensional ice-flow model to explore the development of the Green Bay Lobe, an outlet glacier of the southern Laurentide Ice Sheet, leading up to the time of maximum ice extent and during subsequent deglaciation (c. 30 to 8 cal. ka BP). We focus on conditions at the ice-bed interface in order to evaluate their possible impact on glacial landscape evolution. Air temperatures for model input have been reconstructed using the GRIP δ ¹⁸ O record calibrated to speleothem records from Missouri that cover the time periods of c. 65 to 30 cal. ka BP and 13.25 to 12.4 cal. ka BP. Using that input, the known ice extents during maximum glaciation and early deglaciation can be reproduced reasonably well. The model fails, however, to reproduce short-term ice margin retreat and readvance events during later stages of deglaciation. Model results indicate that the area exposed after the retreat of the Green Bay Lobe was characterized by permafrost until at least 14 cal. ka BP. The extensive drumlin zones that formed behind the ice margins of the outermost Johnstown phase and the later Green Lake phase are associated with modeled ice margins that were stable for at least 1000 years, high basal shear stresses (c. 100 kPa) and permafrost depths of 80-200 m. During deglaciation, basal meltwater and sliding became more important.     

Milankovitch theory of ice ages: Hypothesis of ice-sheet linkage between regional insolation and global climate

J. D. Hays, J. Imbrie, and N. J. Shackleton (1976, Science194, 1121–1132) showed that the astronomical theory explained many features of late Quaternary ice-age climates, but they did not specify the physical mechanisms involved. Here it is proposed that interlocked variations of ice-sheet heat sinks in both polar hemispheres amplified and transmitted Milankovitch summer half-year insolation changes (a version of the astronomical theory) between 45° and 75°N into the globally synchronous climate changes recorded in geologic records. It is suggested that late Quaternary ice sheets had terrestrial components (grounded above sea level, melting margins, fluctuations controlled by climate) and marine components (grounded below sea level, drained largely by ice streams, limited melting margins, fluctuations controlled primarily by sea level and secondarily by climate, interior surface elevations coupled to downdraw through ice streams). Northern Hemisphere ice sheets were largely marine (with minor melting margins) in the Arctic and terrestrial (with major melting margins) in the midlatitudes. West Antarctic and peripheral East Antarctic ice was marine-based and lacked melting margins. Because of their geographic array, these terrestrial and marine components formed an ice-sheet system whose variations were coupled on a global scale. Milankovitch summer isolation changes near midlatitude Northern Hemisphere melting margins controlled most variations of this system, because advance or retreat of melting margins initiated concurrent eustatic sea-level change. Such sea-level change afforded the critical interlocking mechanism between terrestrial and marine components because it forced simultaneous expansion or contraction of marine margins in both polar hemispheres. This initiated an amplifying feedback loop among all marine components and influenced interior downdraw through ice streams. Arctic summer insolation change was less important because northern melting margins were relatively minor. Its greatest influence was on surface ablation of ice streams that controlled interior downdraw. This affected eustatic sea level and activated global linkage of marine sectors. By analogy with present-day Antarctica, late Quaternary ice sheets were enormous planetary heat sinks due to their reflective and radiative surface characteristics. It is suggested that the effectiveness of these ice-sheet heat sinks varied with their areal extent and interior surface elevation. Thus, it is postulated that concurrent growth or decay of these interlocked ice-sheet heat sinks in both polar hemispheres served as the global amplifier of regional Milankovitch summer insolation.     

Dating bedrock gorge incision in the French Western Alps (Ecrins-Pelvoux massif) using cosmogenic 10Be

We report in-situ produced ¹⁰Be data from the Gorge du Diable (French Western Alps) to date and quantify bedrock gorge incision into a glacial hanging valley. We sampled gorge sidewalls and the active channel bed to derive both long-term and present-day incision rates. ¹⁰Be ages of sidewall profiles reveal rapid incision through the late Holocene (ca 5 ka) at rates ranging from 6.5 to 13 mm yr⁻¹. Present-day incision rates are significantly lower and vary from 0.5 to 3 mm yr⁻¹ within the gorge. Our data imply either delayed initiation of gorge incision after final ice retreat from internal Alpine valleys at ca 12 ka, or post-glacial surface reburial of the gorge. Our results suggest that fluvial incision rates >1 cm yr⁻¹ into crystalline bedrock may be encountered in transient landscape features induced by glacial-interglacial transitions.