Ice Margin Research Articles

Total isostatic response to the complete unloading of the Greenland and Antarctic Ice Sheets

The land surface beneath the Greenland and Antarctic Ice Sheets is isostatically suppressed by the mass of the overlying ice. Accurate computation of the land elevation in the absence of ice is important when considering, for example, regional geodynamics, geomorphology, and ice sheet behaviour. Here, we use contemporary compilations of ice thickness and lithospheric effective elastic thickness to calculate the fully re-equilibrated isostatic response of the solid Earth to the complete removal of the Greenland and Antarctic Ice Sheets. We use an elastic plate flexure model to compute the isostatic response to the unloading of the modern ice sheet loads, and a self-gravitating viscoelastic Earth model to make an adjustment for the remaining isostatic disequilibrium driven by ice mass loss since the Last Glacial Maximum. Feedbacks arising from water loading in areas situated below sea level after ice sheet removal are also taken into account. In addition, we quantify the uncertainties in the total isostatic response associated with a range of elastic and viscoelastic Earth properties. We find that the maximum change in bed elevation following full re-equilibration occurs over the centre of the landmasses and is +783 m in Greenland and +936 m in Antarctica. By contrast, areas around the ice margins experience up to 123 m of lowering due to a combination of sea level rise, peripheral bulge collapse, and water loading. The computed isostatic response fields are openly accessible and have a number of applications for studying regional geodynamics, landscape evolution, cryosphere dynamics, and relative sea level change.

Chlorine-36 Surface Exposure Dating of Late Holocene Moraines and Glacial Mass Balance Modeling, Monte Sierra Nevada, South-Central Chilean Andes (38°S)

The development of robust chronologies of Neoglaciation from individual glaciers throughout the high-altitude Andes can provide fundamental knowledge of influences such as regional temperature and precipitation variability, and aid in predicting future changes in the Andean climate system. However, records of Late Holocene glaciation from the Central Chilean Andes are sparse, and often poorly constrained. Here, we present 36Cl surface exposure ages, dendrochronologic constraints, and glacial mass balance modeling simulations of Late Holocene glacier fluctuations in the Central-South Chilean Andes. A series of concentric moraine ridges were identified on Monte Sierra Nevada (38°S), where exposure dating of basaltic boulders was used to establish a chronology of ice recession. We infer that moraine abandonment of the most distal ridge in the valley commenced by ∼4.2 ka, and was followed by glacier margin retreat to an up-valley position. Exposure ages of the oldest Late Holocene boulders (∼2.5–0.8 ka) along the marginal extents of the moraine complex indicate fluctuations of the glacier terminus prior to ∼0.65 ka. A final expansion of the ice margin reoccupied the position of the 4.2 ka moraine, with abatement from the outermost composite moraine occurring by ∼0.70 ka, as constrained by tree-ring data from live Araucaria araucana trees. Finally, a series of nested moraines dating to ∼0.45–0.30 ka, formed from a pulsed ice recession during the latest Holocene when the lower reaches of the glacial snout was most likely debris mantled. A distributed temperature index model combined with a glacier flow model was used to quantify an envelope of possible climatic conditions of Late Holocene glaciation. The glacial modeling results suggest conditions were ∼1.5°C colder and 20% wetter during peak Neoglaciation relative to modern conditions. These records also suggest a near-coeval record of Late Holocene climate variability between the middle and high southern latitudes. Furthermore, this study presents some of the youngest 36Cl exposure ages reported for moraines in the Andes, further supporting this method as a valuable geochronologic tool for assessing Late Holocene landscape development.

Drag forces at the ice-sheet bed and resistance of hard-rock obstacles: the physics of glacial ripping

AbstractGlacial ripping involves glaciotectonic disintegration of rock hills and extensive removal of rock at the ice-sheet bed, triggered by hydraulic jacking caused by fluctuating water pressures. Evidence from eastern Sweden shows that glacial ripping caused significant subglacial erosion during the final deglaciation of the Fennoscandian ice sheet, distinct from abrasion and plucking (quarrying). Here we analyse the ice drag forces exerted onto rock obstacles at the base of an ice sheet, and the resisting forces of such rock obstacles: glaciotectonic disintegration requires that ice drag forces exceed the resisting forces of the rock obstacle. We consider rock obstacles of different sizes, shapes and fracture patterns, informed by natural examples from eastern Sweden. Our analysis shows that limited overpressure events, unfavourable fracture patterns, low-transmissivity fractures, slow ice and streamlined rock hamper rock hill disintegration. Conversely, under fast ice flow and fluctuating water pressures, disintegration is possible if the rock hill contains subhorizontal, transmissive fractures. Rock steps on previously smooth, abraded surfaces, caused by hydraulic jacking, also enhance drag forces and can cause disintegration of a rock hill. Glacial ripping is a physically plausible erosion mechanism, under realistic glaciological conditions prevalent near ice margins.

Chemical weathering signatures at Mt. Achernar, Central Transantarctic Mountains II: Surface exposed sediments

Mt Achernar Moraine is a high altitude, high latitude blue ice moraine where typical conditions preclude the presence of liquid water. Cosmogenic and salt accumulation dating indicate that the moraine’s surface is progressively older away from the active ice margin, with surface exposure ages up to 1 Ma. We analyze the chemical and mineralogical transformations in the <63 µm fraction along transects across the moraine. Data include bulk chemical composition, crystalline mineralogy by X-ray diffraction (XRD), and the composition of amorphous or low abundance products of chemical weathering by sequential extraction. These data are analyzed by multiple regression as a function of exposure age and as a function of composition of the moraine’s cobble and pebble-sized clasts. Change with exposure age is defined by the development of salts and carbonate minerals along with the input of detrital material, principally from sedimentary rocks. Clay minerals and amorphous cements breakdown as detrital material in proportions far above their abundance in the rock clasts, whereas framework silicates (i.e. feldspars and quartz) break down in relatively small proportions. Both the carbonate minerals and some of the salts form from atmospheric acids (i.e. H2CO3) that in turn react with other minerals. Mass balance shows that the input of these atmospheric acids balances with gains in authigenic smectites, zeolites, and amorphous material. Many of these minerals also form in the subglacial environment, but are poorly represented in the underlying rock, suggesting a similar chemical weathering regime in both the subglacial and surface environments of this hyper cold and arid setting. The rate of CO2 drawdown into carbonate minerals increases as the moraine progressively thickens, from 3 mg·m2·a−1 in freshly emerging sediments to ∼50 mg·m2·a−1 after 500 ka of exposure. Weathering from acidic aerosols is proportional to atmospheric flux documented in ice cores and does not vary with moraine thickness. The carbonate mineral formation rates are more than an order of magnitude below those of the subglacial environment and as much as two orders of magnitude below those found in warm desert soils. Nevertheless, the drawdown of atmospheric CO2 into carbonate minerals occurs in a terrestrial setting where water exists only in vapor form.

Unravelling the long-term, locally heterogenous response of Greenland glaciers observed in archival photography

Abstract. We present an approach for extracting quantifiable information from archival aerial photographs to extend the temporal record of change over a region of the central eastern Greenland Ice Sheet. The photographs we use were gathered in the 1930s as part of a surveying expedition, and so they were not acquired with photogrammetric analysis in mind. Nevertheless, we are able to make opportunistic use of this imagery, as well as additional, novel datasets, to explore changes at ice margins well before the advent of conventional satellite technology. The insights that a longer record of ice margin change bring is crucial for improving our understanding of how glaciers are responding to the changing climate. In addition, our work focuses on a series of relatively small and little studied outlet glaciers from the eastern margin of the ice sheet. We show that whilst air and sea surface temperatures are important controls on the rates at which these ice masses change, there is also significant heterogeneity in their responses, with non-climatic controls (such as the role of bathymetry in front of calving margins) being extremely important. In general, there is often a tendency to focus either on changes of the Greenland Ice Sheet as a whole, or on regional variations. Here, we suggest that even this approach masks important variability, and full understanding of the behaviour and response of the ice sheet requires us to consider changes that are taking place at the scale of individual glaciers.

Ice‐Marginal Proglacial Lakes Across Greenland: Present Status and a Possible Future

AbstractIce‐marginal lakes can affect glacier dynamics but are ignored in studies of the evolution of the Greenland ice sheet (GrIS) and of peripheral mountain glaciers and ice caps (PGICs). Here we show that lakes occupy 10% of the GrIS ice margin and occur on 5% of PGICs. Ice velocity at the GrIS margin is enhanced by ∼ 25% at lakes versus on land. Mean ice discharge into lakes is ∼4.9 Gt.yr, which is ∼1% of ice discharged through marine termini. We locate thousands of subglacial overdeepenings within which 7,404 km2 of future lakes could form, all of which will be ice‐marginal at some time. Future lakes in the west and east will be restricted to the margin of the GrIS and within alpine valleys, respectively. This status and possible future leads us to contend that lakes should be incorporated into projections of Greenland ice loss.

Moraines in the Austrian Alps record repeated phases of glacier stabilization through the Late Glacial and the Early Holocene

Climate is currently warming due to anthropogenic impact on the Earth’s atmosphere. To better understand the processes and feedbacks within the climate system that underlie this accelerating warming trend, it is useful to examine past periods of abrupt climate change that were driven by natural forcings. Glaciers provide an excellent natural laboratory for reconstructing the climate of the past as they respond sensitively to climate oscillations. Therefore, we study glacier systems and their behavior during the transition from colder to warmer climate phases, focusing on the period between 15 and 10 ka. Using a combination of geomorphological mapping and beryllium-10 surface exposure dating, we reconstruct ice extents in two glaciated valleys of the Silvretta Massif in the Austrian Alps. The mountain glacier record shows that general deglaciation after the Last Glacial Maximum (LGM) was repeatedly interrupted by glacier stabilization or readvance, perhaps during the Oldest Dryas to Bølling transition (landform age: 14.4 ± 1.0 ka) and certainly during the Younger Dryas (YD; 12.9–11.7 ka) and the Early Holocene (EH; 12–10 ka). The oldest landform age indicates a lateral ice margin that postdates the ‘Gschnitz’ stadial (ca. 17–16 ka) and predates the YD. It shows that local inner-alpine glaciers were more extensive until the onset of the Bølling warm phase (ca. 14.6 ka), or possibly even into the Bølling than during the subsequent YD. The second age group, ca. 80 m below the (pre-)Bølling ice margin, indicates glacier extents during the YD cold phase and captures the spatial and temporal fine structure of glacier retreat during this period. The ice surface lowered approximately 50–60 m through the YD, which is indicative of milder climate conditions at the end of the YD compared to its beginning. Finally, the third age group falls into a period of more substantial warming, the YD–EH transition, and shows discontinuous glacier retreat during the glacial to interglacial transition. The new geochronologies synthesized with pre-existing moraine records from the Silvretta Massif evidence multiple cold phases that punctuated the general post-LGM warming trend and illustrate the sensitive response of Silvretta glaciers to abrupt climate oscillations in the past.

Sedimentary Signatures of Persistent Subglacial Meltwater Drainage From Thwaites Glacier, Antarctica

Subglacial meltwater drainage can enhance localized melting along grounding zones and beneath the ice shelves of marine-terminating glaciers. Efforts to constrain the evolution of subglacial hydrology and the resulting influence on ice stability in space and on decadal to millennial timescales are lacking. Here, we apply sedimentological, geochemical, and statistical methods to analyze sediment cores recovered offshore Thwaites Glacier, West Antarctica to reconstruct meltwater drainage activity through the pre-satellite era. We find evidence for a long-lived subglacial hydrologic system beneath Thwaites Glacier and indications that meltwater plumes are the primary mechanism of sedimentation seaward of the glacier today. Detailed core stratigraphy revealed through computed tomography scanning captures variability in drainage styles and suggests greater magnitudes of sediment-laden meltwater have been delivered to the ocean in recent centuries compared to the past several thousand years. Fundamental similarities between meltwater plume deposits offshore Thwaites Glacier and those described in association with other Antarctic glacial systems imply widespread and similar subglacial hydrologic processes that occur independently of subglacial geology. In the context of Holocene changes to the Thwaites Glacier margin, it is likely that subglacial drainage enhanced submarine melt along the grounding zone and amplified ice-shelf melt driven by oceanic processes, consistent with observations of other West Antarctic glaciers today. This study highlights the necessity of accounting for the influence of subglacial hydrology on grounding-zone and ice-shelf melt in projections of future behavior of the Thwaites Glacier ice margin and marine-based glaciers around the Antarctic continent.

Timing of the last deglaciation phases in the southern Baltic area inferred from Bayesian age modeling

A new chronology of the last Scandinavian Ice Sheet retreat in the southern Baltic basin is proposed. Based on Bayesian age modeling, we show that the most likely ages of particular deglaciation phases are 16.5 ± 0.5 ka for the Gardno Phase, 15.6 ± 0.6 ka for the Słupsk Bank Phase, and 13.9 ± 0.5 ka for the Southern Middle Bank Phase. The Gardno moraines are correlated with the Halland Coastal moraines in southern Sweden and the Middle Lithuanian moraines in Lithuania and Latvia. Ice margin stillstands of the Słupsk Bank Phase and Souhern Middle Bank Phase are correlated with the Göteborg and Vimmerby moraines, and with the North Lithuanian (Haanja) and Otepää moraines. The average retreat rates of the ice margin is about 55 m/yr between the Gardno Phase and the Słupsk Bank Phase, and about 40 m/yr between the Słupsk Bank Phase and the Southern Middle Bank Phase. This suggests that the last deglaciation did not accelerate after the Gardno Phase when an extensive ice-dammed lake was formed in front of the retreating ice sheet. The ice margin was probably grounded rather than floating, which prevented its more rapid retreat. The timing of the two main ice margin stillstands at the Słupsk Bank and at the Southern Middle Bank corresponds to the cool periods around 15 and 14 ka interpreted from paleotemperatures of Greenland based on ice core GISP2. This suggests that the main phases of the last deglaciation in the southern Baltic region were at least partly triggered by climatic fluctuations in the Northern Hemisphere.

Catastrophic ice shelf collapse along the NW Laurentide Ice Sheet highlights the vulnerability of marine-based ice margins

We document the unparalleled depositional and chronological record of an 80,000 km2 ice shelf from the former Laurentide Ice Sheet whose advance and retreat spanned only ∼400 years (11.6–11.2 cal ka BP). Its catastrophic breakup (∼150 years; 11.3–11.2 cal ka BP) coincided with a rapidly warming atmosphere and ocean, paralleling reported conditions impacting the ongoing evolution of marine-based margins of Antarctica and Greenland. Our record highlights the instability of the NW Laurentide Ice Sheet whose marine-based ice streams and ice shelves retreated in channels of intermediate depth (400–600 m) that lack reverse slopes, suggesting that current ice dynamics models may underestimate the sensitivity of similar margins in Antarctica to ongoing global warming. Our reconstruction is based on > 700 km of coastal field surveys complemented by radiocarbon dates stratigraphically encompassing the full lifespan of the former Viscount Melville Sound Ice Shelf. Notably, this detailed record demonstrates that large ice shelves can be removed rapidly, cautioning against dismissing predictions of possible 21st century mean-sea-level rise >1 m derived from the ongoing retreat of marine-based ice sheets. We emphasize that during the last deglaciation of the western Canadian Arctic Archipelago, prominent ice streams transitioned into ice shelves whose depositional landforms are pervasive, demonstrating their pivotal role in the evisceration of the NW Laurentide Ice Sheet. The VMSIS provides an exceptionally detailed geological record whose tightly constrained chronology of growth, maintenance and breakup serves as an instructive example for those modern marine-based margins whose instability is being investigated to address estimates of possible future sea level rise.

Review article: Existing and potential evidence for Holocene grounding line retreat and readvance in Antarctica

Abstract. Widespread existing geological records from above the modern ice sheet surface and outboard of the current ice margin show that the Antarctic Ice Sheet (AIS) was much more extensive at the Last Glacial Maximum (∼ 20 ka) than at present. However, whether it was ever smaller than present during the last few millennia, and (if so) by how much, is known only for a few locations because direct evidence lies within or beneath the ice sheet, which is challenging to access. Here, we describe how retreat and readvance (henceforth “readvance”) of AIS grounding lines during the Holocene could be detected and quantified using subglacial bedrock, subglacial sediments, marine sediment cores, relative sea-level (RSL) records, geodetic observations, radar data, and ice cores. Of these, only subglacial bedrock and subglacial sediments can provide direct evidence for readvance. Marine archives are of limited utility because readvance commonly covers evidence of earlier retreat. Nevertheless, stratigraphic transitions documenting change in environment may provide support for direct evidence from subglacial records, as can the presence of transgressions in RSL records, and isostatic subsidence. With independent age control, ice structure revealed by radar can be used to infer past changes in ice flow and geometry, and therefore potential readvance. Since ice cores capture changes in surface mass balance, elevation, and atmospheric and oceanic circulation that are known to drive grounding line migration, they also have potential for identifying readvance. A multidisciplinary approach is likely to provide the strongest evidence for or against a smaller-than-present AIS in the Holocene.

Direct measurement of optical properties of glacier ice using a photon-counting diffuse LiDAR

AbstractThe production of meltwater from glacier ice, which is exposed at the margins of land ice during the summer, is responsible for a large proportion of glacier mass loss. The rate of meltwater production from glacier ice is especially sensitive to its physical structure and chemical composition which combine to determine the albedo of glacier ice. However, the optical properties of near-surface glacier ice are not well known since most prior work has focused on laboratory-grown ice or deep cores. Here, we demonstrate a measurement technique based on diffuse propagation of nanosecond-duration laser pulses in near-surface glacier ice that enables the independent measurement of the scattering and absorption coefficients, allowing for a complete description of the processes governing radiative transfer. We employ a photon-counting detector to overcome the high losses associated with diffuse optics. The instrument is highly portable and rugged, making it optimally suited for deployment in remote regions. A set of measurements taken on Crook and Collier Glaciers, Oregon, serves as a demonstration of the technique. These measurements provide insight into both physical structure and composition of near-surface glacier ice and open new avenues for the analysis of light-absorbing impurities and remote sensing of the cryosphere.

Hydro‐Mechanical Interactions of a Rock Slope With a Retreating Temperate Valley Glacier

AbstractRock slope failures often result from progressive rock mass damage which accumulates over long timescales. In deglaciating environments, rock slopes are affected by stress perturbations driven by mechanical unloading due to ice downwasting and concurrent changes in thermal and hydraulic boundary conditions. Since in‐situ data are rare, the different processes and their relative contribution to slope damage remain poorly understood. Here, we analyze borehole monitoring data from a rock slope adjacent to the retreating Great Aletsch Glacier (Switzerland) and compare it to englacial water levels, climate data, and decreasing ice levels. Rock slope pore pressures show a seasonal signal controlled by infiltration events as well as effects from the connectivity to the englacial hydrological system. We find that reversible and irreversible strains are driven by: (a) hydromechanical effects caused by englacial pressure fluctuations and infiltration events, (b) stress transfer related to changing mechanical glacial loads from short‐term englacial water level fluctuations and longer term ice downwasting, and (c) thermomechanical effects from annual temperature cycles penetrating the shallow subsurface, which primarily result in reversible deformation. We relate most observed irreversible strain (damage) to mechanical unloading from ice downwasting. Damage is strongest directly at the ice margin and moves through the slope at the pace of glacial retreat and advance. Locations with many retreat/advance cycles are very sensitive to landslide formation. The current climate warming impacts very sensitive valley sectors, which is confirmed by landslide distributions and activity in the study area.

A 10Be-dated record of glacial retreat in Connemara, Ireland, following the Last Glacial Maximum and implications for regional climate

Late Pleistocene stadials are traditionally viewed as periods of severe, year-round cooling centred on the North Atlantic, where the abrupt onset and termination of stadial conditions is widely, though not unequivocally, thought to reflect the impact on Atlantic Meridional Circulation of meltwater fluxes from fringing ice masses. Widespread subaerial melting of land-based ice, however, implies a state of glacier-climate disequilibrium that, according to some conceptual models, had to persist for the duration of each stadial event rather than immediately prior to it. Yet, to date, the spatial extent of any such melting remains unclear, as do the magnitude and duration of any cognate atmospheric warming. To help address this ambiguity, we reconstructed the timing and nature of deglaciation of the Connemara ice centre, Ireland, which, being located immediately downwind of the North Atlantic, affords an ideal location for assessing terrestrial expressions of stadial and non-stadial conditions. We report 15 new cosmogenic beryllium-10 ages of glacial erratic boulders that reveal the timing and magnitude of fluctuations in this sector of the former Irish ice sheet during Heinrich Stadial 1 (HS1). Coupled with geomorphic mapping, these results indicate rapid, widespread deglaciation of the Connemara ice centre at approximately 17 ka, under full stadial conditions, and that the ice margin retreated back to regional source areas rather than eastward towards the centre of the island. Recognising the potential divergence from the traditional view of cold North Atlantic stadials, we propose that the rate and nature of deglaciation in Connemara during HS1 likely reflects enhanced subaerial melting during the summer ablation season. • Fifteen new 10 Be dates constrain the deglaciation of an ice centre in Connemara, Ireland to Heinrich Stadial 1. • Our chronology contrasts with the current paradigm and suggests summer melting occurred despite anomalous winter cooling. • The Connemara deglacial record supports a broader pattern of enhanced seasonality in the North Atlantic during stadials.

Ice Phenology in Eurasian Lakes over Spatial Location and Altitude

Eurasian freezing lakes cover an almost 180° wide longitude sector between the latitudes 30° and 75° N, and their altitudes range from below the sea surface level up to 5 km elevation. Ice phenology varies widely in this region. However, these variations and their influence factors have been little studied. Analytic models are applied here to examine these variations supported by historical ice and weather data. These models are forced by a linear air–lake heat exchange formula based on local empirical fits. The weather brings latitude–longitude–altitude patterns to the large-scale lake ice phenology. Freezing and breakup dates are forced by the local air temperature and solar radiation, and their rates of change are also important. In addition, freezing depends on lake depth and breakup depends on accumulated ice thickness. Lake depth provides a lag and radiation balance provides a shift with respect to the air temperature in cooling of the lake, and breakup is dictated by spring warming conditions and ice thickness. Due to solar radiation forcing, the common degree-day approach is biased for modelling ice phenology, especially in low latitudes. Analytic models provide a first-order tool for climate sensitivity of ice seasons. The freezing date and breakup date both change by around five days per one-degree shift in air temperature away from the climatological ice margin; however, at this margin, the sensitivity is higher.

Meighen Ice Cap: changes in geometry, mass, and climatic response since 1959

Remote sensing (satellite and airborne) and in situ monitoring measurements indicate that the Meighen Ice Cap has thinned by 10 ± 3 m between 1960 and 2016, resulting in areal shrinkage by 32 km2 (38%) and total mass loss of 0.71 ± 0.2 gigatonnes. Retreat of the ice margin along the North basin by up to 2.5 km accounted for 50% of the total area loss over the 56-year period of this study. A strong inverse relationship between atmospheric summer temperature anomalies, and in situ mass balance measurements (r = −0.74) indicates that accelerated rates of mass loss from the Meighen Ice Cap are largely driven by regional scale summer warming. Increases in summer temperature anomalies by 1.9 °C from 1960–2004 to 2005–2016 coincided with a five-fold increase in surface melt across the Meighen Ice Cap. Since the 1990s, the Equilibrium Line Altitude has been positioned above the summit of the Meighen Ice Cap, resulting in loss of the accumulation zone that persisted for the first three decades of this study. Since the early 2000s, the ice cap summit has thinned by more than 4 m, which is unprecedented in the period of record. Response time calculations based on 2006–2016 average mass balance data indicate that the Meighen Ice Cap, which has been in existence for ∼3500 years, will completely disappear in ∼175 years.

PISM-LakeCC: Implementing an adaptive proglacial lake boundary in an ice sheet model

Abstract. During the Late Pleistocene and Holocene retreat of paleo-ice sheets in North America and Europe, vast proglacial lakes existed along the land terminating margins. These proglacial lakes impacted ice sheet dynamics by imposing boundary conditions analogous to a marine terminating margin. Such lacustrine boundary conditions cause changes in the ice sheet geometry, stress balance and frontal ablation and therefore affect the mass balance of the entire ice sheet. Despite this, dynamically evolving proglacial lakes have rarely been considered in detail in ice sheet modeling endeavors. In this study, we describe the implementation of an adaptive lake boundary in the Parallel Ice Sheet Model (PISM), which we call PISM-LakeCC. We test our model with a simplified glacial retreat setup of the Laurentide Ice Sheet (LIS). By comparing the experiments with lakes to control runs with no lakes, we show that the presence of proglacial lakes locally enhances the ice flow, which leads to a lowering of the ice sheet surface. In some cases, this also results in an advance of the ice margin and the emergence of ice lobes. In the warming climate, increased melting on the lowered ice surface drives the glacial retreat. For the LIS, the presence of lakes triggers a process similar to marine ice sheet instability, which caused the collapse of the ice saddle over Hudson Bay. In the control experiments without lakes, Hudson Bay is still glaciated when the climate reaches present-day (PD) conditions. The results of our study demonstrate that glacio-lacustrine interactions play a significant role in the retreat of land terminating ice sheet margins.

Rapid ice sheet response to deglacial and Holocene paleoenvironmental changes in eastern Prydz Bay, East Antarctica

Prydz Bay lies at the terminus of one of East Antarctica's largest glacial systems and is a key region for understanding the response of the ice sheet to past and future environmental changes. In this study, we explore the dynamics and paleo-geometry of the ice sheet in eastern Prydz Bay, using a combination of bathymetric features on the seafloor to delineate past flow patterns, and cosmogenic nuclide dating of glacial deposits on land to constrain ice sheet terminus chronologies. Large streamlined bedforms on the sea floor record the existence of a primary ice stream in Svenner Channel, collecting ice from multiple tributary ice streams and discharging into the main trunk stream of the Lambert Glacier-Amery Ice Shelf system in Prydz Channel. The location and orientation of the north-eastern tributary to this ice stream also provide evidence for a substantial independent ice dome on the outer shelf at Four Ladies Bank. Exposure of ice-free areas at outer Rauer Group and Vestfold Hills indicates that grounding line retreat across eastern Prydz Bay was largely complete by ∼14 ka BP, and the ice margin had retreated to within ∼1 km of its present position by ∼10 ka BP. Onshore moraines record periods of ice margin retreat during the middle (∼6 ka BP) and very late (∼0.5 ka BP) Holocene coinciding with local warm periods. Subsidence recorded in modern Global Navigation Satellite System (GNSS) observations suggests the regional ice sheet was smaller than at present in the period between the middle and very late Holocene advances. The unusually dynamic ice sheet behaviour in this area is attributed to the smooth, reverse-slope bed characteristics of eastern Prydz Bay, which enabled rapid retreat of the ice sheet margin several hundred kilometres during deglacial events.

Blue-ice moraines formation in the Heritage Range, West Antarctica: Implications for ice sheet history and climate reconstruction

Blue ice is found in areas of Antarctica where katabatic winds, focussed by steep surface slopes or by topography around nunataks, cause enhanced surface ablation. This process draws up deeper, older ice to the ice sheet surface, often bringing with it englacial sediment. Prevailing theories for dynamically stable moraine surfaces in East Antarctica suggest that: (i) it is this material, once concentrated, that forms blue-ice moraines (BIM), (ii) that the moraine formation can be dated using cosmogenic isotope approaches, and that, (iii) since we expect an increase in exposure age moving away from the ice margin towards bedrock, dating across the moraine can be used to constrain ice-sheet history. To test this lateral accretion model for BIM formation we visited Patriot, Marble and Independence Hills in the southern Heritage Range, West Antarctica. Detailed field surveys of surface form, sediment and moraine dynamics were combined with geophysical surveys of the englacial structure of the moraines and cosmogenic nuclide analysis of surface clasts. Results suggest sediment is supplied mainly by basal entrainment, supplemented by debris-covered valley glaciers transferring material onto the ice sheet surface, direct deposition from rock-fall and slope processes from nunataks. We find that once sediment coalesces in BIM, significant reworking occurs through differential ablation, slope and periglacial processes. We bring these processes together in a conceptual model, concluding that many BIM in West Antarctica are dynamic and, whilst they persist through glacial cycles, they do not always neatly record ice sheet retreat patterns since linear distance from the ice margin does not always relate to increased clast exposure age. Understanding the dynamic processes involved in moraine formation is critical to the effective interpretation of the typically large scatter of cosmogenic nuclide exposure ages, opening a deep window into the million-year history of the West Antarctic Ice Sheet.

Millennial-scale migration of the frozen/melted basal boundary, western Greenland ice sheet

AbstractThe geometry and thermal structure of western Greenland ice sheet are known to have undergone relatively substantial change over the Holocene. Evolution of the frozen and melted fractions of the bed associated with the ice-sheet retreat over this time frame remains unclear. We address this question using a thermo-mechanically coupled flowline model to simulate a 11 ka period of ice-sheet retreat in west central Greenland. Results indicate an episode of ~100 km of terminus retreat corresponded to ~16 km of upstream frozen/melted basal boundary migration. The majority of migration of the frozen area is associated with the enhancement of the frictional and strain heating fields, which are accentuated toward the retreating ice margin. The thermally active bedrock layer acts as a heat sink, tending to slow contraction of frozen-bed conditions. Since the bedrock heat flux in our region is relatively low compared to other regions of the ice sheet, the frozen region is relatively greater and therefore more susceptible to marginward changes in the frictional and strain heating fields. Migration of melted regions thus depends on both geometric changes and the antecedent thermal state of the bedrock and ice, both of which vary considerably around the ice sheet.