Search

Subglacial processes exert a major control on ice streaming. Constraining subglacial conditions thus allows for more accurate predictions of ice mass loss. Due to the difficulty in observing large‐scale conditions of the modern subglacial environment, we turn to geological records of ice streaming in deglaciated environments. Morphometric values of streamlined subglacial bedforms provide valuable information about the relative speed, direction, and maturity of past ice streams and the relationship between ice streaming and subglacial erosion and deposition. However, manually identifying streamlined subglacial bedforms across deglaciated landscapes, sometimes in clusters of several thousand, is an arduous task with difficult‐to‐control sources of variability and human‐biased errors. This paper presents a new tool that utilizes a machine learning approach to automatically identify glacially derived streamlined features. Slope variations across a landscape, identified by topographic position index, undergo analysis from a series of supervised machine learning models trained from over 600 000 data points identified across the deglaciated Northern Hemisphere. A filtered data set produced through the combination of scientifically driven preprocessing and statistical downsampling improved the robustness of our approach. After cross‐validation, we found that Random Forest detected the most true positives, up to 94.5% on a withheld test set, and an ensemble average of machine learning models provided the highest stability when applied within the range of applicable data sets, performing at up to 79% identification of true positives on an out of distribution area of interest. We build these models into an open‐source Python package, bedfinder, and apply it to new data in the Green Bay Lobe region, USA, finding the general ice‐flow direction and average streamlined subglacial bedform elongation with minimal effort. This type of open, reproducible machine learning analysis is at the leading edge of glacial geomorphology research and will continue to improve with integration of newly acquired and previously collected data. [ABSTRACT FROM AUTHOR]

Reliable projections of future sea level rise from the polar ice sheets depend on the ability of ice sheet models to accurately reproduce flow dynamics in an evolving ice sheet system. Ice sheet models are sensitive to the choice of the basal sliding law, which remains a significant source of uncertainty. In this study we apply a range of sliding laws to a hindcast model of Jakobshavn Isbræ, western Greenland, from 2009 to 2018. We demonstrate that a linear viscous sliding law requires the assimilation of regular velocity observations into the model in order to reproduce the observed large seasonal and inter-annual variations in flow speed. This requirement introduces a major limitation for producing accurate future projections. A regularised Coulomb friction law, in which basal traction has an upper limit, is able to more accurately reproduce the range of speeds from 2012 to 2015, the period of peak flow and maximal retreat, without the requirement for assimilating regular observations. Additionally, we find evidence that the speed at which sliding transitions between power-law and Coulomb regimes may vary spatially and temporally. These results point towards the possible form of an ideal sliding parameterisation for accurately modelling fast-flowing glaciers and ice streams, although determining this is beyond the scope of this study. [ABSTRACT FROM AUTHOR]

Despite their recognized significance on global climate and extensive research efforts, the mechanism(s) driving Heinrich events remain(s) a subject of debate. Here, we use the 3D thermomechanically coupled glacial systems model (GSM) to examine the Hudson Strait ice stream surge cycling and the role of three factors previously hypothesized to play a critical role in Heinrich events: ice shelves, glacial isostatic adjustment, and sub-surface ocean temperature forcings. In contrast to all previous modeling studies examining HEs, the GSM uses a transient last glacial cycle climate forcing, global viscoelastic glacial isostatic adjustment model, and sub-glacial hydrology model. The results presented here are based on a high-variance sub-ensemble retrieved from North American history matching for the last glacial cycle. Over our comparatively wide sampling of the potential parameter space (52 ensemble parameters for climate forcing and process uncertainties), we find two modes of Hudson Strait ice streaming: classic binge–purge versus near-continuous ice streaming with occasional shutdowns and subsequent surge onset overshoot. Our model results indicate that large ice shelves covering the Labrador Sea during the last glacial cycle only occur when extreme calving restrictions are applied. The otherwise minor ice shelves provide insignificant buttressing for the Hudson Strait ice stream. While sub-surface ocean temperature forcing leads to minor differences regarding surge characteristics, glacial isostatic adjustment does have a significant impact. Given input uncertainties, the strongest controls on ice stream surge cycling are the poorly constrained deep geothermal heat flux under Hudson Bay and Hudson Strait and the basal drag law. Decreasing the geothermal heat flux within available constraints and/or using a Coulomb sliding law instead of a Weertman-type power law leads to a shift from the near-continuous streaming mode to the binge–purge mode. [ABSTRACT FROM AUTHOR]

We analyse ambient-noise seismic data from 23 three-component seismic nodes to study firn velocity structure and seismic anisotropy near the EastGRIP camp along the Northeast Greenland Ice Stream (NEGIS). Using nine-component correlation tensors, we derive dispersion curves of Rayleigh and Love wave group velocities from 3 to 40 Hz. These velocity distributions exhibit anisotropy along and across the flow. To assess these variations, we invert dispersion curves for shear wave velocities (Vsh and Vsv) in the top 150 m of the NEGIS using a Markov chain Monte Carlo approach. The reconstructed 1-D shear velocity model reveals radial anisotropy in the firn, with Vsh 12 %–15 % greater than Vsv , peaking at the critical density (550 kg m−3). We combine density data from firn cores drilled in 2016 and 2018 to create a new density parameterisation for the NEGIS, serving as a reference for our results. We link seismic anisotropy in the NEGIS to effective and intrinsic causes. Seasonal densification, wind crusts, and melt layers induce effective anisotropy, leading to faster Vsh waves. Changes in firn recrystallisation cause intrinsic anisotropy, altering the Vsv / Vsh ratio. We observe a shallower firn–ice transition across the flow (≈ 50 m) compared with along the flow (≈ 60 m), suggesting increased firn compaction due to the predominant wind direction and increased deformation towards the shear margin. We demonstrate that short-duration (9 d minimum), passive, seismic deployments and noise-based analysis can determine seismic anisotropy in firn, and we reveal 2-D firn structure and variability. [ABSTRACT FROM AUTHOR]

The grounding line position of glaciers and ice shelves is an essential observation for the study of the Earth's ice sheets. However, in some locations, such as the Antarctic Peninsula, where many grounding lines have not been mapped since the 1990s, remote sensing of grounding line position remains challenging. Here we present a tidal motion offset correlation (TMOC) method for measuring the grounding line position of tidewater glaciers and ice shelves, based on the correlation between tide amplitude and synthetic aperture radar offset tracking measurements. We apply this method to the Antarctic Peninsula Ice Sheet to automatically delineate a new grounding line position for 2019–2020, with near complete coverage along 9300 km of coastline, updating the 20-year-old record. A comparison of the TMOC grounding line to contemporaneous interferometrically measured grounding line position shows the method has a mean seaward offset compared to interferometry of 185 m and a standard deviation of 295 m. Our results show that over the last 24 years there has been grounding line retreat at a number of fast-flowing ice streams on the Antarctic Peninsula, with the most retreat concentrated in the north-eastern sector, where grounding lines have retreated following the collapse of ice shelves. We observe a maximum grounding line retreat since 1996 of 16.3 ± 0.5 km on Hektoria Glacier, with other notable glaciers retreating by 9.3 ± 0.5, 9.1 ± 0.5 and 3.6 ± 0.5 km. Our results document dynamic change on Antarctic Peninsula glaciers and show the importance of using an updated grounding line location to delineate the boundary between floating and grounded ice. [ABSTRACT FROM AUTHOR]

Regional subsurface mapping of glacial depositional systems preserved in buried bedrock paleovalleys, and quantitative analysis of new LiDAR imagery of surface glacial landforms using machine learning techniques, when combined, are powerful tools for assessing the dynamics of the Laurentide Ice Sheet (LIS) during the last (Wisconsinan) glaciation in southern Ontario. While age dating of deposits preserved below Last Glacial Maximum tills (LGM: marine isotope stage (MIS) 2 < c.24 000 years B.P. (ybp)) is still sparse, newly available sedimentological data derived by cored drilling, combined with legacy outcrop data, identify thick (100 m+) successions of glaciolacustrine sediments and a lack of till(s), indicating that the ice sheet margin did not extend beyond the Niagara Escarpment at the western end of Lake Ontario, during the earliest phases of the glaciation (MIS 4) or the ensuing mid-Wisconsinan (MIS 3). Ice was able to extend into New York State blocking the Rome outlet to the Hudson Valley ponding deep proglacial lakes in the glacio-isostatically depressed Huron–Ontario–Erie basins recorded by thick glaciolacustrine sediments in paleovalleys. These were cannibalized by an expanding Late Wisconsinan ice sheet after ∼24 000 ybp recorded by extensive till sheets resting on a marked erosional unconformity, with drumlinized surfaces. Analysis and visualization of LiDAR data identifies discrete statistically validated flow sets of highly elongated streamlined bedforms (mega-scale glacial lineations (MSGLs)). These provide key evidence of a major reorganization of the ice sheet margin during deglaciation into lobate paleo ice streams shortly after 17 400 ybp. MSGLs are cut across earlier LGM drumlinized tills creating widespread "palimpsest" surfaces. At least two principal phases of fast ice flow can be identified, marked by large fluxes of sediment and the rapid building of large gravel and sand-dominated moraine complexes within interlobate depocentres, the largest glacial landforms in southern Ontario. Analysis of LiDAR data further reveals the common presence of DeGeer moraines where ice margins retreated in water, and iceberg scours. Future work using LiDAR mapping has the objective of fully documenting the number, extent, and timing of ice streams to enhance glaciological modelling when the ice sheet rapidly lost mass. [ABSTRACT FROM AUTHOR]

Upernavik Isstrøm, the largest contributor to sea-level rise in northwest Greenland, has experienced complex and contrasting ice-flow-speed changes across its five outlets over the last two decades. In this study, we present a detailed remote-sensing analysis of the ice dynamics at Upernavik's outlets from 2000 to 2021 to evaluate the details of these changes. Previous research suggested that the presence or absence of floating ice tongues strongly influences Upernavik's ice dynamics. We use several lines of evidence to document the presence of floating ice tongues, and find that, while several outlets experienced ice-tongue formation and/or loss during the study period, these changes do not explain observed fluctuations in ice-flow velocity. Further exploration of ice-dynamic forcings using a flowline model suggests that changes in basal slipperiness near the terminus have a strong impact on upstream ice dynamics and can explain the velocity variations. Our results suggest that speed fluctuations at Upernavik's outlets may be seasonally and interannually controlled by bed conditions near the terminus, and highlight the need for further research on the influence of basal conditions on complex tidewater glacier dynamics. [ABSTRACT FROM AUTHOR]

Antarctica has an active subglacial hydrological system, with interconnected subglacial lakes fed by subglacial meltwater. Subglacial hydrology can influence basal sliding, inject freshwater into the sub‐ice‐shelf cavity, and impact sediment transport and deposition which can affect the stability of grounding lines (GLs). We used satellite altimetry data from the ICESat, ICESat‐2, and CryoSat‐2 missions to document the second recorded drainage of Engelhardt Subglacial Lake (SLE), which began in July 2021 and discharged more than 2.3 km3 of subglacial water into the Ross Ice Shelf cavity. We used differential synthetic aperture radar interferometry from RADARSAT‐2 and TerraSAR‐X alongside ICESat‐2 repeat‐track laser altimetry (RTLA) and REMA digital elevation model strips to detect 2–13 km of GL retreat since the previous drainage event in 2003–06. Combining these satellite observations, we evaluated the mechanism triggering SLE drainage, the cause of the observed GL retreat, and the interplay between subglacial hydrology and GL dynamics. We find that: (a) SLE drainage was initiated by influx from a newly identified upstream lake; (b) the observed GL retreat is mainly driven by the continued retreat of Engelhardt Ice Ridge and long‐term dynamic thinning that caused a grounded ice plain to reach flotation; and (c) SLE drainage and GL retreat were largely independent. We also discuss the possible origins and influence of a 27 km grounded promontory found to protrude seaward from the GL. Our observations demonstrate the importance of high‐resolution satellite data for improving the process‐based understanding of dynamic and complex regions around the Antarctic Ice Sheet margins. Plain Language Summary: Large volumes of water flow beneath the Antarctic Ice Sheet through an interconnected network of rivers and lakes. This water system impacts slipperiness at the base of the ice, affecting how fast it moves. It also delivers freshwater into the ocean, directly contributing to sea‐level rise and increasing melt beneath the floating ice shelves. In this study, we use satellite data to track the 2021–24 drainage of Engelhardt Subglacial Lake in West Antarctica. This lake is located close to the Ross Ice Shelf grounding line, the point at the edge of the ice sheet where the ice first lifts off the bedrock and starts to float on the ocean. This region of the grounding line has retreated by up to 13 km since the last lake drainage in 2003–06. Here, we investigate what caused the drainage, the reasons for the grounding line retreat, and whether the two processes are connected. We also report the growth of a grounded promontory extending 27 km out to sea, which may be evidence of a former ice stream moraine or melt channel. These findings help to improve our limited understanding of the relationship between subglacial hydrology and grounding line dynamics in Antarctica. Key Points: Satellite altimetry detects second recorded drainage of Engelhardt Subglacial Lake in July 2021, triggered by an influx from an upstream lakeThe grounding line retreated by 2–13 km since the last drainage in 2003, linked to the retreat of Engelhardt Ice Ridge and ice plain ungroundingSatellite observations suggest that the 2021–24 lake drainage and grounding line retreat were largely independent processes [ABSTRACT FROM AUTHOR]

Integrating ice‐sheet models with empirical data pertaining to palaeo‐ice sheets promotes advances in the models used in sea‐level predictions and can improve our understanding of past ice‐sheet behaviour. The large number of empirical constraints on the last British–Irish Ice Sheet make it ideal for model–data comparison experiments. Here, we present an ensemble of 600 model simulations, which are compared with data on former ice‐flow extent, flow geometry and deglaciation timing. Simulations which poorly recreate data were ruled out, allowing us to examine the remaining physically realistic simulations which capture the ice sheets' behavioural tendencies. Our results led to a novel reconstruction of behaviour in the data‐poor region of the North Sea, insights into the ice stream, potential ice‐shelf and readvance dynamics, and the potential locations of peripheral ice caps. We also propose that the asynchronous behaviour of the British–Irish Ice Sheet is a consequence of the geography of the British Isles and the merging and splitting of different bodies of ice through saddle merger and collapse. Furthermore, persistent model–data mismatches highlight the need for model development, especially regarding the physics of ice–ocean interactions. Thus, this work highlights the power of integrating models and data, a long‐held aim of palaeoglaciology. [ABSTRACT FROM AUTHOR]

Knowledge of the behaviour of marine‐based ice sheets during times of climatic warming, such as the last deglaciation, provides important information to understand how ice sheets respond to external forcing. We analysed swath bathymetric and acoustic sub‐bottom profiler data from Wrigley Gulf on the western Amundsen Sea shelf, West Antarctica, to identify glacial features and reconstruct past changes in the extent of the West Antarctic Ice Sheet (WAIS) and ice flow directions. Glacial bedforms mapped within a bathymetric cross‐shelf trough include features showing cross‐cutting and overprinting relationship and indicate changes in ice‐flow orientation. Here, we distinguish at least two phases of different ice‐flow patterns on the Wrigley Gulf shelf. During the earlier phase, seaward ice stream flow on the inner shelf was deflected towards the east due to the existence of an ice dome on the middle‐outer continental shelf. Retreat of grounded ice towards the centre of this dome is indicated by the asymmetric cross profile of recessional moraines mapped on the middle shelf. The later glaciation phase was characterized by fast, NNW‐directed ice flow across the shelf along a broad front and subsequent stepwise landward retreat, which is evident from the common occurrence and orientation of mega‐scale glaciation lineations and grounding zone wedges on the middle‐inner shelf. It is uncertain whether the two phases of glaciation recorded on the seafloor occurred during the last and penultimate glacial periods or at different times of the last glaciation. Reliable chronological constraints from sediment cores and additional geomorphological information are needed to understand the cause of the changes in WAIS dynamics reflected by the two ice‐flow phases. [ABSTRACT FROM AUTHOR]

Ice-core records from the interior of the Greenland ice sheet suggest widespread thinning during the Holocene. However, the recurring underestimation of this thinning in numerical models raises concerns about both the veracity of such reconstructions and the reliability of glaciological models. Recent work suggests the 8000-year-old Northeast Greenland Ice Stream (NEGIS), including a now-extinct northern tributary, may have been an early influence on Greenland ice-sheet dynamics. Yet, the inaccurate reproduction of NEGIS-like dynamics in most models hampers investigation of whether this feature played a role in Holocene ice-sheet thinning. Here we show that grounding-line retreat in northeast Greenland triggers elevation changes at the northern summit via ice-dynamic effects modulated by the paleo NEGIS system. In our simulations, fast ice-stream flow caused by transiently imposed reduced basal shear stress following the northeast retreat explains 55% (± 18%) of the estimated ice thinning, showing that ice-stream dynamics is one of the main drivers of the NGRIP Holocene surface elevation drop. Our findings show that the ice-flow in northeast Greenland plays a large role in ice-surface elevation changes in central Greenland. Using a three-dimensional ice sheet model, this work demonstrates that more than half of the Holocene surface elevation drop recorded in ice cores from north-central Greenland can be attributed to ice dynamic changes in northeast Greenland following ice-sheet retreat during the last deglaciation. [ABSTRACT FROM AUTHOR]

Subglacial lakes have been mapped across Antarctica with two methods, radio‐echo sounding (RES) and ice‐surface deformation. At sites where both are coincident, these methods typically provide conflicting interpretations about the ice‐bed interface. With a single exception, active subglacial lakes identified by surface deformation do not display the expected flat, bright, and specular bed reflection in RES data, characteristic of non‐active lakes. This observational conundrum suggests that our understanding of Antarctic subglacial hydrology, especially beneath important fast‐moving ice streams, remains incomplete. Here, we use an airborne RES campaign that surveyed a well‐characterized group of active subglacial lakes on lower Mercer and Whillans ice streams, West Antarctica, to explore inconsistency between the two observational techniques. We test hypotheses of increased scattering and attenuation due to the presence of an active subglacial lake system that could suppress reflected bed‐echo power for RES observations in these locations, finding that entrained water is most plausible. Plain Language Summary: The bottom of an ice sheet is insulated from cold air temperatures, often warm enough to melt and pond liquid water into lakes. These lakes beneath the ice sheet have been identified by two independent measurements, first with radar methods and second with changes in height of the ice surface (altimetry). Interestingly, the two methods rarely identify the same lakes: radar generally detects lakes in the ice‐sheet interior, whereas altimetry detects active lakes near the ice‐sheet margins that fill and drain within the time series of repeated measurements (∼years). In this study, we investigate a group of active subglacial lakes at which both radar and altimetry data sets are available. We demonstrate that the radar returns from active lake reflections are much dimmer than expected based on non‐active lake signatures and investigate the physical processes controlling those dim reflections. We argue that water moves into the ice when the lake fills or drains and that is the most plausible explanation for the observational discrepancy. Key Points: Active subglacial lakes, identified by surface deformation, do not create the expected bright and specular radar reflectionEntrained water in basal ice suppresses radar power by scattering and attenuation, and it also likely alters the basal ice mechanicsUnderstanding the radar expression of subglacial water on Earth provides context for investigations of subsurface water on planetary bodies [ABSTRACT FROM AUTHOR]

The Last Glacial Maximum (LGM) was characterised by huge ice sheets covering the Northern Hemisphere, especially over North America, and by its cold climate. Previous authors have performed numerical simulations of the LGM to better understand coupled climate–ice sheet systems. However, the results of such simulations are sensitive to many model parameters. Here, we perform a 200-member ensemble of simulations of the North American and Greenland ice sheets and climate of the LGM with a coupled ice sheet–atmosphere–slab ocean model (FAMOUS-BISICLES) to explore sensitivities of the coupled climate–ice system to 16 uncertain parameters. In the ensemble of simulations, the global mean surface temperature is primarily controlled by the combination of parameters in the large-scale condensation scheme and the cumulus convection scheme. In simulations with plausible LGM global mean surface temperatures, we find that the albedo parameters have only a small impact on the Greenland ice volume due to the limited area of surface ablation associated with the cold climate. Instead, the basal sliding law controls the ice volume by affecting ice transport from the interior to the margin. On the other hand, like the Greenland ice sheet in future climate change, the LGM North American ice sheet volume is controlled by parameters in the snow and ice albedo scheme. Few of our simulations produce an extensive North American ice sheet when the global temperature is above 12 °C. Based on constraints on the LGM global mean surface temperature, the ice volume and the southern extent of the North American ice sheet, we select 16 acceptable simulations. These simulations lack the southern extent of ice compared to reconstructions, but they show reasonable performance on the ice sheet configuration and ice streams facing Baffin Bay and the Arctic Ocean. The strong sensitivities of the North American ice sheet to albedo at the LGM may imply a potential constraint on the future Greenland ice sheet by constraining the albedo schemes. [ABSTRACT FROM AUTHOR]

Glaciers and ice streams flowing over sediment beds commonly have a layer of ice‐rich debris adhered to their base, known as a "frozen fringe," but its impact on basal friction is unknown. We simulated basal slip over granular beds with a cryogenic ring shear device while ice infiltrated the bed to grow a fringe, and measured the frictional response under different effective stresses and slip speeds. Frictional resistance increased with increasing slip speed until it plateaued at the frictional strength of the till, closely resembling the regularized Coulomb slip law associated with clean ice over deformable beds. We hypothesize that this arises from deformation in a previously unidentified zone of weakly frozen sediments at the fringe's base, which is highly sensitive to temperature and stress gradients. We show how a rheologic model for ice‐rich debris coupled with the thermomechanics of fringe growth can account for the regularized Coulomb behavior. Plain Language Summary: Many glaciers move by sliding over sediment beds. As the glacier flows downslope, ice can infiltrate the underlying sediments, forming a layer of ice‐rich debris attached at the glacier's base. We investigated how this frozen fringe impacts glacier motion by simulating glacier slip in a cold‐room facility with a specialized ring shear device. We recreated glacier conditions, sliding ice over granular beds to form the fringe, and then assessed how frictional resistance at the slip interface varied under different stresses and ice speeds. We found that frozen fringe influences the relationship between ice speed and frictional resistance, known as the "slip law." As ice speed increased, basal friction increased to a threshold that matches the strength of the ice‐free sediment bed—mirroring the "regularized Coulomb slip law" inferred for clean ice over soft beds. We attribute this behavior to deformation in a weakly frozen zone at the base of the frozen fringe and show how this behavior can be incorporated into existing parameterizations of glacier slip. Key Points: Ring shear experiments show frozen fringe alters basal slip dynamics for soft‐bedded glaciersDeformation in a zone of weakly frozen sediments within the fringe leads to a regularized Coulomb slip response [ABSTRACT FROM AUTHOR]

The Amundsen Sea sector has some of the fastest-thinning ice shelves in Antarctica, caused by high, ocean-driven basal melt rates, which can lead to increased ice streamflow, causing increased sea level rise (SLR) contributions. In this study, we present the results of a new synchronously coupled ice-sheet–ocean model of the Amundsen Sea sector. We use the Wavelet-based, Adaptive-grid, Vertically Integrated ice sheet model (WAVI) to solve for ice velocities and the Massachusetts Institute of Technology general circulation model (MITgcm) to solve for ice thickness and three-dimensional ocean properties, allowing for full mass conservation in the coupled ice–ocean system. The coupled model is initialised in the present day and run forward under idealised warm and cold ocean conditions with a fixed ice front. We find that Thwaites Glacier dominates the future SLR from the Amundsen Sea sector, with a SLR that evolves approximately quadratically over time. The future evolution of Thwaites Glacier depends on the lifespan of small pinning points that form during the retreat. The rate of melting around these pinning points provides the link between future ocean conditions and the SLR from this sector and will be difficult to capture without a coupled ice–ocean model. Grounding-line retreat leads to a progressively larger Thwaites Ice Shelf cavity, leading to a positive trend in total melting, resulting from the increased ice basal surface area. Despite these important sensitivities, Thwaites Glacier retreats even in a scenario with zero ocean-driven melting. This demonstrates that a tipping point may have been passed in these simulations and some SLR from this sector is now committed. [ABSTRACT FROM AUTHOR]

The ice streams flowing into the Amundsen Sea, West Antarctica, are losing mass due to changes in oceanic basal melting of their floating ice shelves. Rapid ice‐shelf melting is sustained by the delivery of warm Circumpolar Deep Water to the ice‐shelf cavities, which is first supplied to the continental shelf by an undercurrent that flows eastward along the shelf break. Temporal variability of this undercurrent controls ice‐shelf basal melt variability. Recent work shows that on decadal timescales the undercurrent variability opposes surface wind variability. Using a regional model, we show that undercurrent variability is induced by sea‐ice freshwater fluxes, particularly those north of the shelf break, which affect the cross‐shelf break density gradient. This sea‐ice variability is linked to tropical Pacific variability impacting atmospheric conditions over the Amundsen Sea. Ice‐shelf melting also feeds back onto the undercurrent by affecting the on‐shelf density, thereby influencing shelf‐break density gradient anomalies. Plain Language Summary: The glaciers that flow toward the Amundsen Sea, West Antarctica, are losing ice faster than most others about the continent. Once these glaciers reach the coast, they extend out onto the ocean surface, forming ice shelves. The rapid loss of ice is caused by changes in melting by relatively warm ocean waters beneath the floating ice shelves. In the Amundsen Sea, a deep ocean current is responsible for delivering warm water from the deep ocean to the ice shelves. We present model results that show that this deep current varies on decadal timescales as a consequence of systematic sea‐ice melt and formation patterns. A faster current drives more rapid ice shelf melting which, via a feedback process, further accelerates the current. Climate variability originating in the tropical Pacific Ocean is responsible for the variability in the sea‐ice, and is therefore also responsible for the effects on melting of the ice shelves. Key Points: In the Amundsen Sea decadal variability of an undercurrent flowing along the shelf break drives decadal variability in ice‐shelf basal meltSea‐ice freshwater fluxes and positive feedbacks from ice‐shelf basal melt drive the undercurrent variabilityTropical Pacific teleconnections induce atmospheric anomalies over the Amundsen Sea which drive the sea‐ice freshwater flux variability [ABSTRACT FROM AUTHOR]

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

Ice cliffs are melt hot spots that contribute disproportionately to melt on debris-covered glaciers. In this study, we investigate the impact of supraglacial stream hydrology on ice cliffs using in situ and remote sensing observations, streamflow measurements, and a conceptual geomorphic model of ice cliff backwasting applied to ice cliffs on Kennicott Glacier, Alaska. We found that 33 % of ice cliffs (accounting for 69 % of the ice cliff area) are actively influenced by streams, while half are nearer than 10 m from the nearest stream. Supraglacial streams contribute to ice cliff formation and maintenance by horizontal meandering, vertical incision, and debris transport. These processes produce an undercut lip at the ice cliff base and transport clasts up to tens of centimeters in diameter, preventing reburial of ice cliffs by debris. Stream meander morphology reminiscent of sedimentary river channel meanders and oxbow lakes produces sinuous and crescent ice cliff shapes. Stream avulsions result in rapid ice cliff collapse and local channel abandonment. Ice cliffs abandoned by streams are observed to be reburied by supraglacial debris, indicating a strong role played by streams in ice cliff persistence. We also report on a localized surge-like event at the glacier's western margin which drove the formation of ice cliffs from crevassing; these cliffs occur in sets with parallel linear morphologies contrasting with the crescent planform shape of stream-driven cliffs. The development of landscape evolution models may assist in quantifying the total net effect of these processes on steady-state ice cliff coverage and mass balance, contextualizing them with other drivers including supraglacial ponds, differential melt, ice dynamics, and collapse of englacial voids. [ABSTRACT FROM AUTHOR]

Enhanced sedimentation at glacial margins can produce submarine megaslides (>10,000 km3). We report a single megaslide in the Surveyor Fan, Gulf of Alaska. Minimum extant size is ∼16,124 km2 in area and ∼9,080 km3 in volume. Slope failure occurred ∼1.2 Ma at the onset of the mid‐Pleistocene transition (MPT). With accretion along the Aleutian‐Alaska Trench, the original volume is conservatively ∼16,280 km3, with only a 140 km run‐out due to its blocky, high shear strength nature. We suggest the megaslide was triggered by a major sediment influx at the onset of the MPT, when glacial‐interglacial cycles shifted from 41 to 100 Kyr. The absence of repeat megaslides may reflect a changing balance between seismic strengthening and sediment flux, where later sedimentation driven by cross‐shelf ice streams results in thin, fluidized, non‐cohesive slides. Continued accretion of the Surveyor Fan and megaslide also reduces critical wedge taper, further inhibiting major failure. Plain Language Summary: Sediment flux at glacial margins can produce submarine slides >10,000 km3 in size (megaslides). We report a single megaslide in the Surveyor Fan, Gulf of Alaska. Minimum extant size is ∼16,124 km2 in area and ∼9,080 km3 in volume. Failure occurred ∼1.2 Ma at the onset of the mid‐Pleistocene transition (MPT). Due to accretion along the Aleutian‐Alaska Trench, the original volume is conservatively ∼16,280 km3, but with only a 140 km run‐out due to its blocky, high shear strength nature. We suggest the megaslide was caused by the initial major sediment flux at the onset of the MPT, and that afterward only fluidized, thin, and non‐cohesive slides occurred. This is due to changing balance between seismic strengthening and sediment flux, and accretion of the Surveyor Fan and megaslide which reduces the critical wedge taper inhibiting major failure. Key Points: World's Fifth largest mapped megaslide documented beneath Surveyor Fan Gulf of AlaskaTiming of slope failure linked to onset of mid‐Pleistocene glacial intensificationAbsence of later failures due to changing balance of sediment flux/seismic strengthening and negative feedbacks from critical wedge processes [ABSTRACT FROM AUTHOR]

Ice shelves are assumed to flow steadily from their grounding lines to the ice front. We report the detection of ice‐propagating extensional Lamb (plate) waves accompanied by pulses of permanent ice shelf displacement observed by co‐located Global Navigation Satellite System receivers and seismographs on the Ross Ice Shelf. The extensional waves and associated ice shelf displacement are produced by tidally triggered basal slip events of the Whillans Ice Stream, which flows into the ice shelf. The propagation velocity of 2,800 m/s is intermediate between shear and compressional ice velocities, with velocity and particle motions consistent with predictions for extensional Lamb waves. During the passage of the Lamb waves the entire ice shelf is displaced about 60 mm with a velocity more than an order of magnitude above its long‐term flow rate. Observed displacements indicate a peak dynamic strain of 10−7, comparable to that of earthquake surface waves that trigger ice quakes. Plain Language Summary: Ice shelves normally flow steadily toward their boundaries with the open ocean at the ice front. However, seismographs and Global Navigation Satellite System receivers deployed on the Ross Ice Shelf record guided elastic plate waves traveling in the ice as well as permanent displacement of the ice shelf. The elastic waves and ice shelf displacement originate from basal slip events of the Whillans Ice Stream, which flows into the Ross Ice Shelf. The velocity of the elastic waves is about 2,800 m/s, as expected for guided plate waves propagating in an ice shelf. During the passage of the elastic waves, the entire ice shelf with an area of 500,000 square kilometers is displaced about 60 mm in a direction away from the Whillans Ice Stream. These observations show that the strain imparted to the ice shelf by the once or twice daily Whillans Ice Stream basal slip events is sufficient to trigger ice quakes and perhaps enhance the deformation of the ice shelf. Key Points: Extensional Lamb waves propagate across the Ross Ice Shelf, radiated from slip events at the base of the Whillans Ice StreamDuring the passage of the Lamb waves, the entire ice shelf is displaced about 60 mm, with a velocity an order of magnitude above its long‐term flow rateThe displacement pulses produce a peak dynamic strain of 10−7, suggesting that they could trigger icequakes in the ice shelf [ABSTRACT FROM AUTHOR]

Given the high number and diversity of events in a typical cryoseismic dataset, in particular those recorded on ice sheet margins, it is desirable to use a semi-automated method of grouping similar events for reconnaissance and ongoing analysis. We present a workflow for employing semi-unsupervised cluster analysis to inform investigations of the processes occurring in glaciers and ice sheets. In this demonstration study, we make use of a seismic event catalogue previously compiled for the Whillans Ice Stream, for the 2010–2011 austral summer outlined in Part 1,. We address the challenges of seismic event analysis for a complex wave field by clustering similar seismic events into groups using characteristic temporal, spectral, and polarization attributes of seismic time series with the k -means++ algorithm. This provides the basis for a reconnaissance analysis of a seismic wave field that contains local events (from the ice stream) set in an ambient wave field that itself contains a diversity of signals (mostly from the Ross Ice Shelf). As one result, we find that two clusters include stick-slip events that diverge in terms of length and initiation locality (i.e., central sticky spot and/or the grounding line). We also identify a swarm of high-frequency signals on 16–17 January 2011 that are potentially associated with a surface melt event from the Ross Ice Shelf. Used together with the event detection presented in Part 1, the semi-automated workflow could readily be generalized to other locations and, as a possible benchmark procedure, could enable the monitoring of remote glaciers over time and comparisons between locations. [ABSTRACT FROM AUTHOR]

Cryoseismology is a powerful toolset for progressing the understanding of the structure and dynamics of glaciers and ice sheets. It can enable the detection of hidden processes such as brittle fracture, basal sliding, transient hydrological processes, and calving. Addressing the challenge of detecting signals from many different processes, we present a novel approach for the semi-automated detection of events and event-like noise, which is well-suited for use as Part 1 of a workflow where unsupervised machine learning will be used as Part 2 to facilitate the main reconnaissance of diverse detected event types. Implemented in the open-source and widely used ObsPy Python package, the multi-STA/LTA algorithm constructs a hybrid characteristic function from a set of short-term average (sta)–long-term average (lta) pairs (refer to Sect. in the main text for an explanation of how uppercase and lowercase STA/sta and LTA/lta abbreviations are differentiated). We apply the algorithm to data from a seismic array deployed on the Whillans Ice Stream (WIS) in West Antarctica (austral summer 2010–2011) to form a "catch-all" catalogue of events and event-like noise. The new algorithm compares favorably with standard approaches, yielding a diversity of seismic events, including all previously identified stick-slip events , teleseisms, and other noise-type signals. In terms of a catalogue overview, we investigate a partial association of seismicity with the tidal cycle and a slight association with ice temperature changes of the Antarctic summer. The new algorithm and workflow will assist in the comparison of different glacier environments using seismology, the identification of process change over time, and the targeting of possible subsequent high-resolution studies. [ABSTRACT FROM AUTHOR]

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

Ice surface speed increases dramatically from upstream to downstream in many ice streams and glaciers. This speed-up is thought to be associated with a transition from internal distributed deformation to highly localized deformation or sliding at the ice-bedrock interface. The physical processes governing this transition remain unclear. Here, we argue that highly localized deformation does not necessarily initiate at the ice-bedrock interface, but could also take the form of an internal shear band inside the ice flow that connects topographic highs. The power-law exponent n in the ice rheology amplifies the feedback between shear heating and shear localization, leading to the spontaneous formation of an internal shear band that can create flow separation within the ice. We model the thermomechanical ice flow over a sinusoidal basal topography by building on the high-resolution Stokes solver FastICE v1.0. We compile a regime diagram summarizing cases in which a sinusoidal topography with a given amplitude and wavelength leads to shear band formation for a given rheology. We compare our model results to borehole measurements from Greenland and find evidence to support the existence of an internal shear band. Our study highlights the importance of re-evaluating the degree to which internal deformation contributes to total deformation in the ice column and to the flow-to-sliding transition. Plain Language Summary On its way toward the ocean, ice speeds up dramatically from less than 1 m/year upstream to more than a kilometer per year downstream. In this paper, we investigate the physical processes controlling this speed-up. Specifically, we focus on the role of the bedrock topography and rheology in facilitating the transition from this slow to rapid motion. We use a two-dimensional numerical model to simulate the flow field within a slab of ice flowing down a ramp over a simplified topography. We find that including the bedrock topography can lead to a zone of highly localized deformation within the ice above topographic highs. We also find that a non-linear rheology amplifies this localization. We compare our model results to borehole measurements from Greenland and find evidence that supports the existence of this highly localized deformation zone. This study indicates that the localized deformation induced by bedrock topography and amplified by non-linear rheology could be one physical mechanism that governs the speed-up of the ice flow. [ABSTRACT FROM AUTHOR]

We observe the evacuation of 11-year-old landfast sea ice in the Larsen B embayment on the East Antarctic Peninsula in January 2022, which was in part triggered by warm atmospheric conditions and strong offshore winds. This evacuation of sea ice was closely followed by major changes in the calving behaviour and dynamics of a subset of the ocean-terminating glaciers in the region. We show using satellite measurements that, following a decade of gradual slow-down, Hektoria, Green, and Crane glaciers sped up by approximately 20 %–50 % between February and the end of 2022, each increasing in speed by more than 100 ma-1. Circumstantially, this is attributable to their transition into tidewater glaciers following the loss of their ice shelves after the landfast sea ice evacuation. However, a question remains as to whether the landfast sea ice could have influenced the dynamics of these glaciers, or the stability of their ice shelves, through a buttressing effect akin to that of confined ice shelves on grounded ice streams. We show, with a series of diagnostic modelling experiments, that direct landfast sea ice buttressing had a negligible impact on the dynamics of the grounded ice streams. Furthermore, we suggest that the loss of landfast sea ice buttressing could have impacted the dynamics of the rheologically weak ice shelves, in turn diminishing their stability over time; however, the accompanying shifts in the distributions of resistive stress within the ice shelves would have been minor. This indicates that this loss of buttressing by landfast sea ice is likely to have been a secondary process in the ice shelf disaggregation compared to, for example, increased ocean swell or the drivers of the initial landfast sea ice disintegration. [ABSTRACT FROM AUTHOR]

The retreat of the West Antarctic Ice Sheet (WAIS) in the Ross Sea after the Last Glacial Maximum (LGM) was more significant than for any other Antarctic sector. Here we combined the available chronology of retreat with new mapping of seismically resolvable grounding zone wedges (GZWs). Mapping GZWs is important because they record the locations of former stillstands in the extent of grounded ice for individual ice streams during the overall retreat. Our analysis shows that the longest stillstands occurred early in the deglacial period and had millennial durations. Stillstands ended abruptly with retreat distances measured in the tens to hundreds of kilometers creating deep embayments in the extent of grounded ice across the Ross Sea. The location of embayments shifted through time. The available chronological data show that cessation of WAIS and East Antarctic Ice Sheet (EAIS) stillstands was highly asynchronous across at least 5000 years. There was a general shift to shorter stillstands throughout the deglacial period. The asynchronous collapse of individual catchments during the deglacial period suggests that the Ross Sea sector would have contributed to multiple episodes of relatively small-amplitude sea-level rise as the WAIS and EAIS retreated from the region. The high sinuosity of the modern grounding zone in the Ross Sea suggests that this style of retreat persists. [ABSTRACT FROM AUTHOR]

A realistic initialization of ice flow models is critical for predicting future changes in ice sheet mass balance and their associated contribution to sea level rise. The initial thermal state of an ice sheet is particularly important, as it controls ice viscosity and basal conditions, thereby influencing the overall ice velocity. Englacial and subglacial conditions, however, remain poorly understood due to insufficient direct measurements, which complicate the initialization and validation of thermal models. Here, we investigate the impact of using different geothermal heat flux (GHF) datasets and vertical velocity profiles on the thermal state of the Antarctic ice sheet and compare our modeled temperatures to in situ measurements from 15 boreholes. We find that the temperature profile is more sensitive to vertical velocity than to GHF. The basal temperature of grounded ice and the amount of basal melting are influenced by both selection of GHF and vertical velocity. More importantly, we find that the standard approach, which consists of combining basal sliding speed and incompressibility to derive vertical velocities, provides reasonably good results in fast-flow regions (ice velocity >50 m yr -1) but performs poorly in slow-flow regions (ice velocity <50 m yr -1). Furthermore, the modeled temperature profiles in ice streams, where bed geometry is generated using a mass conservation approach, show better agreement with observed borehole temperatures compared to kriging-based bed geometry. [ABSTRACT FROM AUTHOR]

Only a few localised ice streams drain most of the ice from the Greenland Ice Sheet. Thus, understanding ice stream behaviour and its temporal variability is crucially important to predict future sea-level change. The interior trunk of the 700 km-long North-East Greenland Ice Stream (NEGIS) is remarkable due to the lack of any clear bedrock channel to explain its presence. Here, we present a 3-dimensional analysis of the folding and advection of its stratigraphic horizons, which shows that the localised flow and shear margins in the upper NEGIS were fully developed only ca 2000 years ago. Our results contradict the assumption that the ice stream has been stable throughout the Holocene in its current form and show that upper NEGIS-type development of ice streaming, with distinct shear margins and no bed topography relationship, can be established on time scales of hundreds of years, which is a major challenge for realistic mass-balance and sea-level rise projections. The flow of ice streams leaves traces in the stratigraphy of the ice sheets. Made visible by radar, they reveal the history of the upper North East Greenland Ice Stream. The ice stream is found to have existed in its current form for only about the last 2000 years. [ABSTRACT FROM AUTHOR]

The process of icing involves the freezing of overflow layers, on ground or within streams, and results in ice bodies called "aufeis" that are common in most northern landscapes. Knowledge about aufeis is still limited despite the cold region engineering challenge they represent. Understanding the causes of overflow events leading to aufeis development represents a key for the prediction, mitigation, and management of this geohazard and can also support the planning and design of infrastructure in the North. This paper introduces a practical classification for the diverse range of overflow processes that generate aufeis, including under-represented processes, such as the instability of winter streamflow. Importantly, it distinguishes flow conveyance from water supply overflow processes and describes the temporal aspect of icing intensity. Finally, research topics are proposed to improve our understanding of aufeis, including their predictability, the impact of climate change on their occurrence and extent, and stream morphology–aufeis interaction. [ABSTRACT FROM AUTHOR]

This study developed a novel, detailed sedimentological analysis for the complex interactions between rainout, iceberg rafting, tractional underflows, and settling of fines along a glacially influenced basin margin. The glaciomarine interval of the El Imperial Formation (Pennsylvanian, Serpukhovian–Bashkirian) in the San Rafael basin comprises massive to stratified diamictites, interpreted as rainout tills, thinly bedded diamictites, associated with cohesive debris flows, and mudstones containing ice-rafted debris (IRD), all capped by postglacial, transgressive, fine-grained sediments. The rhythmic intercalation of IRD-bearing (dropstone mudstones) and IRD-free (mudstones) intervals likely indicates variations in debris content within the ice margins, the on-and-off switching of ice streams, or dynamic oscillations of the ice terminus. The glaciomarine deposits exhibit soft sediment deformation on both large (metric to decametric) and small (centimetric) scales. This contribution refines previous interpretations of the soft sediment deformation, discerning between loading and slope triggered deformation. Large-scale deformation is characterized by coherent slump folds with low dispersion in the orientations of fold axial plane vergence and fold b-axes. Downslope-verging folds indicate a northward paleoslope, consistent with paleoflow indicators from flute casts found in sandstone turbidite beds. The diamictites affected by the large-scale soft sediment deformation are interpreted as rainout tills with a variable degree of gravity remobilization. Their association with thinly bedded diamictites and laminated mudstones with dropstones suggests that ice rafting played a significant role in the deposition of this succession. [ABSTRACT FROM AUTHOR]

Reconstructions of palaeo-ice stream activity provide an insight into the processes governing ice stream evolution over millennial timescales. The northwestern sector of the Laurentide Ice Sheet experienced a period of rapid retreat driven by warming during the Bølling–Allerød (14.7 – 12.9 ka) which may have contributed significantly to global mean sea level rise during this time. It therefore provides an opportunity to investigate ice sheet dynamics during a phase of rapid ice sheet retreat. Here, we classify coherent groups of ice flow parallel lineations into 326 flowsets and then categorise them as ice stream, deglacial, inferred deglacial or event type flowsets. Combined with ice marginal landforms and a new ice margin chronology (Dalton et al., 2023), we present the first reconstruction of ice flow dynamics of the northwestern Laurentide Ice Sheet at 500-year timesteps through the last deglaciation. At the local Last Glacial Maximum (17.5 ka), the ice stream network was dominated by large, marine-terminating ice streams (>1000 km long) that were fed by the Laurentide-Cordilleran ice saddle to the south and the Keewatin Ice Dome to the east. As the ice margin retreated onshore, the drainage network was characterised by shorter, land-terminating ice streams (<200 km long), with the exception of the Bear Lake and Great Slave Lake ice streams (~600 km long) that terminated in large glacial lakes. Rapid reorganisation of the ice drainage network, from predominantly northerly ice flow to westerly ice flow, occurred over ~2000 years, coinciding with a period of rapid ice sheet surface lowering in the ice saddle region. We note a peak in ice stream activity during the Bølling–Allerød that we suggest is a result of increased ablation and a steepening of the ice surface slope in ice stream onset zones and the increase in driving stresses which contributed to rapid ice drawdown. The subsequent cessation of ice stream activity by the end of the Bølling–Allerød was a result of ice drawdown lowering the ice surface profile, reducing driving stresses and leading to widespread ice stream shut-down. [ABSTRACT FROM AUTHOR]

The grounding line position of glaciers and ice shelves is an essential observation for the study of the Earth's ice sheets. However, in some locations, such as the Antarctic Peninsula, where many grounding lines have not been mapped since the 1990s, remote sensing of grounding line position remains challenging. Here we present a tidal motion offset correlation (TMOC) method for measuring the grounding line position of tidewater glaciers and ice shelves, based on the correlation between tide amplitude and synthetic aperture radar offset tracking measurements. We apply this method to the Antarctic Peninsula Ice Sheet to automatically delineate a new grounding line position for 2019–2020, with near complete coverage along 9,300 km of coastline, updating the 20-year-old record. A comparison of the TMOC grounding line to contemporaneous interferometrically-measured grounding line position shows the method has a mean seaward offset compared to interferometry of 185 m and a standard deviation of 295 m. Our results show that over the last 24 years there has been grounding line retreat at a number of fast flowing ice streams on the Antarctic Peninsula, with the most retreat concentrated in the north-eastern sector, where grounding lines have retreated following the collapse of ice shelves. We observe a maximum grounding line retreat since 1996 of 16.3 km on Hektoria Glacier, with other notable glaciers retreating by 9.3 km, 9.1 km, and 3.6 km respectively. Our results document dynamic change on Antarctic Peninsula glaciers and show the importance of using an updated grounding line location to delineate the boundary between floating and grounded ice. [ABSTRACT FROM AUTHOR]

The Baltic Ice Stream, a large fast‐flowing sector of the last Fennoscandian Ice Sheet that occupied the present‐day Baltic Sea basin, was first conceptualized in the earliest days of glacial geological research in Scandinavia. Landform and sedimentological evidence from the terrestrial margins support the concept and numerical ice‐sheet models demonstrate its existence and possible evolution. However, with evidence for the Baltic Ice Stream thus far limited to the terrestrial periphery, its true form, scale, function, and role in deglaciation have proven enigmatic. Here we present geomorphological evidence directly from the Baltic seabed that confirms the existence of and sheds light on the behaviour of the Baltic Ice Stream. Based on an extensive collection of high‐, moderate‐ and low‐resolution bathymetric terrain models covering a large proportion of the Baltic Sea floor, and complemented by LiDAR‐data for the Baltic islands, we have identified and mapped >20 000 individual subglacial bedforms, meltwater landforms and grounding line landforms. We reconstruct a six‐stage sequence of ice flow and retreat, finding that streaming was persistent in the Baltic but that pathways were variable in extent, timing and duration: different sectors of the Baltic exhibit asynchronous streaming and out‐of‐phase grounding line changes. During deglaciation, grounding line re‐advances occurred in both the southwestern and the northern Baltic Proper, and, while abundant iceberg ploughmarks attest to calving as a significant ice loss mechanism, lobate margins suggest supply to the Baltic catchment was consistently high. Our reconstruction is limited by a fragmentary geomorphic record. Here we put forward a first hypothesis for how the Baltic Ice Stream evolved, and hope it stimulates new geomorphic, stratigraphical and core data collection to extend the landform record, provide insights into subglacial and grounding line processes, and constrain the chronology for Baltic Ice Stream flow and retreat. [ABSTRACT FROM AUTHOR]

We estimate the seismic attenuation of P and S waves in the polar firn and underlying ice by spectral analysis of diving, refracted, and reflected waves from active-source three-component seismic signals obtained in 2010 on the Whillans Ice Stream (WIS), a fast-flowing ice stream in West Antarctica. The resulting quality factors are then successfully modeled using a rock-physics theory of wave propagation that combines White's mesoscopic attenuation theory of interlayer flow with that of Biot/squirt flow. The first theory describes an equivalent viscoelastic medium consisting of a stack of two alternating thin porous layers, both of which have thicknesses that are much greater than the pore size but smaller than the wavelength. On the other hand, in the so-called Biot/squirt-flow model, there are two loss mechanisms, namely the global Biot flow and the local flow from fluid-filled microcracks (or grain contacts) to the pore space and back, where the former is dominant over the latter. The fluid saturating the pores is assumed to be fluidized snow, defined as a mixture of snow particles and air, such as powder, with a rigidity modulus of zero. [ABSTRACT FROM AUTHOR]

To better constrain meltwater transport and ice viscosity in temperate glaciers, particularly in ice stream shear margins, we use a custom permeameter to study the untested model relationship between the permeability of temperate ice and its liquid water content. The permeability of lab-made ice of two mean grain diameters (1.8 and 4.2 mm) is measured, and water content is controlled with the ice salinity and measured calorimetrically. Fluorescein dye is added to through-flowing, chilled water to highlight flow pathways through the ice after experiments. As predicted by a simple model, permeability increases with approximately the square of the water content and by about three orders of magnitude across water contents of 0.1–4.4%. However, permeability values are less than those of the model by average factors of 2.6 and 4.1 for the finer and coarser ice, respectively. This discrepancy is likely due to tortuous, truncated or air-clogged veins. The order-of-magnitude agreement between measured and modeled values may indicate that reduced permeability from these factors is nearly compensated by preferential flow in oversized veins that are isolated or arborescent. Both kinds of preferred flow pathways are observed but the latter only in fine-grained ice at water contents > 2%.

Spatial variability in bed topography, characterized as bed roughness, impacts ice-sheet flow and organization and can be used to infer subglacial conditions and processes, yet is difficult to quantify due to sparse observations. Paleo-subglacial beds of formerly expanded glaciers found across the Antarctic continental shelf are well preserved, have relatively limited post-glacial sediment cover and contain glacial landforms that can be resolved at sub-meter vertical scales. We analyze high-resolution bathymetry offshore of Pine Island and Thwaites glaciers in the Amundsen Sea to explore spatial variability of bed roughness where streamlined subglacial landforms allow for the determination of ice-flow direction. We quantify bed roughness using std dev. and Fast Fourier Transform methods, each employed at local (100 km) and regional (101–2 km) scales and in along- and across-flow orientations to determine roughness expressions across spatial scales. We find that the magnitude of roughness is impacted by the parameters selected – which are often not sufficiently reported in studies – to quantify roughness. Important spatial patterns can be discerned from high-resolution bathymetry, highlighting both its usefulness in identifying patterns of streaming ice flow and underscores the need for a standardized way of characterizing topographic variability.

Ice streams are regions of rapid ice sheet flow characterised by a high degree of sliding over a deforming bed. The shallow shelf approximation (SSA) provides a convenient way to obtain closed-form approximations of the velocity and flux in a rapidly sliding ice stream when the basal drag is much less than the driving stress. However, the validity of the SSA approximation breaks down when the magnitude of the basal drag increases. Here we find a more accurate expression for the velocity and flux in this transitional regime before vertical deformation fully dominates, in agreement with numerical results. The closed-form expressions we derive can be incorporated into wider modelling efforts to yield a better characterisation of ice stream dynamics, and inform the use of the SSA in large-scale simulations.

Investigations of the time-dependent behavior of marine ice sheets and their sensitivity to basal conditions require numerical models because existing theoretical analyses focus only on steady-state configurations primarily with a power-law basal shear stress. Numerical results indicate that the choice of the sliding law strongly affects ice-sheet dynamic behavior. Although observed or simulated grounding-line retreat is typically interpreted as an indication of marine ice sheet instability introduced by Weertman (1974), this (in)stability is a characteristic of the ice sheet's steady states – not time-variant behavior. To bridge the gap between theoretical and numerical results, we develop a framework to investigate grounding line dynamics with generalized basal and lateral stresses (i.e. the functional dependencies are not specified). Motivated by observations of internal variability of the Southern Ocean conditions we explore the grounding-line response to stochastic variability. We find that adding stochastic variability to submarine melt rates that produced stable steady-state configurations leads to intermittently advancing and retreating grounding lines. They can also retreat in an unstoppable manner on time-scales significantly longer than the stochastic correlation time-scales. These results suggest that at any given time of their evolution, the transient behavior of marine ice sheets cannot be described in terms of ‘stable’ or ‘unstable’.

Increases in ocean temperatures in the Filchner Ronne region of Antarctica are likely to result in increased ice mass loss and sea level rise. We constrain projections of the 21st century sea level contribution of this region using process-based ice-sheet modeling, with model parameters controlling ice dynamics calibrated using observed surface speeds and Markov-chain Monte Carlo sampling. We use climate forcing from Representative Concentration Pathway (RCP) scenarios as well as a set of hypothetical scenarios of deep ocean warming to evaluate the sensitivity of this region to ocean temperatures. Projected changes in regional ice mass correspond to a decrease in global mean sea level of 24±7 mm over 2015–2100 under RCP 2.6 and 28±9 mm under RCP 8.5. Increased regional inland surface accumulation related to higher warming levels in RCP 8.5 leads to more ice above flotation, offsetting increased ice shelf basal melt. The tests involving step changes in ocean temperatures with constant surface forcing show that one degree of ocean warming from present results in an additional +11 mm contribution to sea level by 2100 and 1% of the ice-covered area in the domain becomes ungrounded (23 200 km2). The rate of mass loss with temperature increases at higher temperatures.

Upernavik Isstrøm, the largest contributor to sea-level rise in northwest Greenland, has experienced complex and contrasting ice-flow-speed changes across its five outlets over the last two decades. In this study, we present a detailed remote-sensing analysis of the ice dynamics at Upernavik's outlets from 2000 to 2021 to evaluate the details of these changes. Previous research suggested that the presence or absence of floating ice tongues strongly influences Upernavik's ice dynamics. We use several lines of evidence to document the presence of floating ice tongues, and find that, while several outlets experienced ice-tongue formation and/or loss during the study period, these changes do not explain observed fluctuations in ice-flow velocity. Further exploration of ice-dynamic forcings using a flowline model suggests that changes in basal slipperiness near the terminus have a strong impact on upstream ice dynamics and can explain the velocity variations. Our results suggest that speed fluctuations at Upernavik's outlets may be seasonally and interannually controlled by bed conditions near the terminus, and highlight the need for further research on the influence of basal conditions on complex tidewater glacier dynamics.

The dynamics of the Greenland Ice Sheet are affected by surface meltwater reaching the base of the ice, altering ice contact with the bedrock. Lack of understanding of this evolution hampers the ability to predict the effects of increasing temperatures on the Greenland Ice Sheet mass balance. Here we present a unique high-resolution study of ice velocity response to surface melting based on data from a COSMO-SkyMed satellite campaign over Upernavik Isstrøm (Northwest Greenland) for two months around the end of the 2014 melt season. We show that the velocity variations, due to both short-term (days) and seasonal variations in surface melt rates, are increasing in relative strength farther from the glacier terminus. Furthermore, we observe how ice dynamic response to frontal retreat, reaching several kilometres inland, can obscure the meltwater-induced velocity change close to the terminus. Future studies should consider the flow velocity dependence on the distance to the terminus, and local geometry, to distinguish subglacial hydrologic system changes from frontal processes and local basal conditions.

To better constrain meltwater transport and ice viscosity in temperate glaciers, particularly in ice stream shear margins, we use a custom permeameter to study the untested model relationship between the permeability of temperate ice and its liquid water content. The permeability of lab-made ice of two mean grain diameters (1.8 and 4.2 mm) is measured, and water content is controlled with the ice salinity and measured calorimetrically. Fluorescein dye is added to through-flowing, chilled water to highlight flow pathways through the ice after experiments. As predicted by a simple model, permeability increases with approximately the square of the water content and by about three orders of magnitude across water contents of 0.1–4.4%. However, permeability values are less than those of the model by average factors of 2.6 and 4.1 for the finer and coarser ice, respectively. This discrepancy is likely due to tortuous, truncated or air-clogged veins. The order-of-magnitude agreement between measured and modeled values may indicate that reduced permeability from these factors is nearly compensated by preferential flow in oversized veins that are isolated or arborescent. Both kinds of preferred flow pathways are observed but the latter only in fine-grained ice at water contents > 2%. [ABSTRACT FROM AUTHOR]

Spatial variability in bed topography, characterized as bed roughness, impacts ice-sheet flow and organization and can be used to infer subglacial conditions and processes, yet is difficult to quantify due to sparse observations. Paleo-subglacial beds of formerly expanded glaciers found across the Antarctic continental shelf are well preserved, have relatively limited post-glacial sediment cover and contain glacial landforms that can be resolved at sub-meter vertical scales. We analyze high-resolution bathymetry offshore of Pine Island and Thwaites glaciers in the Amundsen Sea to explore spatial variability of bed roughness where streamlined subglacial landforms allow for the determination of ice-flow direction. We quantify bed roughness using std dev. and Fast Fourier Transform methods, each employed at local (100 km) and regional (101–2 km) scales and in along- and across-flow orientations to determine roughness expressions across spatial scales. We find that the magnitude of roughness is impacted by the parameters selected – which are often not sufficiently reported in studies – to quantify roughness. Important spatial patterns can be discerned from high-resolution bathymetry, highlighting both its usefulness in identifying patterns of streaming ice flow and underscores the need for a standardized way of characterizing topographic variability. [ABSTRACT FROM AUTHOR]

Ice streams are regions of rapid ice sheet flow characterised by a high degree of sliding over a deforming bed. The shallow shelf approximation (SSA) provides a convenient way to obtain closed-form approximations of the velocity and flux in a rapidly sliding ice stream when the basal drag is much less than the driving stress. However, the validity of the SSA approximation breaks down when the magnitude of the basal drag increases. Here we find a more accurate expression for the velocity and flux in this transitional regime before vertical deformation fully dominates, in agreement with numerical results. The closed-form expressions we derive can be incorporated into wider modelling efforts to yield a better characterisation of ice stream dynamics, and inform the use of the SSA in large-scale simulations. [ABSTRACT FROM AUTHOR]

Ice streams are major contributors to ice sheet mass loss and sea level rise. Effects of their dynamic behaviour are imprinted into seismic properties, such as wave speeds and anisotropy. Here, we present results from a distributed acoustic sensing (DAS) experiment in a deep ice-core borehole in the onset region of the Northeast Greenland Ice Stream, with focus on phenomenological and methodological aspects. A series of active seismic surface sources produced clear recordings of the P and S wavefield, including internal reflections, along a 1500 m long fibre-optic cable that was placed into the borehole. The combination of nonlinear traveltime tomography with a firn model constrained by multimode surface wave data, allows us to invert for P and S wave speeds with depth-dependent uncertainties on the order of only 10 m s−1, and vertical resolution of 20–70 m. The wave speed model in conjunction with the regularly spaced DAS data enable a straightforward separation of internal upward reflections followed by a reverse-time migration that provides a detailed reflectivity image of the ice. While the differences between P and S wave speeds hint at anisotropy related to crystal orientation fabric, the reflectivity image seems to carry a pronounced climatic imprint caused by rapid variations in grain size. Further improvements in resolution do not seem to be limited by the DAS channel spacing. Instead, the maximum frequency of body waves below ∼200 Hz, low signal-to-noise ratio caused by poor coupling, and systematic errors produced by the ray approximation, appear to be the leading-order issues. Among these, only the latter has a simple existing solution in the form of full-waveform inversion. Improving signal bandwidth and quality, however, will likely require a significantly larger effort in terms of both sensing equipment and logistics. [ABSTRACT FROM AUTHOR]

While subglacial hydrology is known to play a role in glacial dynamics on sub-annual to decadal scales, it remains unclear whether subglacial hydrology plays a critical role in ice sheet evolution on centennial or longer timescales. Furthermore, several drainage systems have been inferred, but it is unclear which is most applicable at the continental/glacial scale. More fundamentally, it is even unclear if the structural choice of subglacial hydrology truly matters for this context. Here we compare the contribution to the surge behaviour of an idealized Hudson Strait-like ice stream from three subglacial hydrology systems. We use the newly updated BAsal Hydrology Model – BrAHMs2.0 – and provide model verification tests. BrAHMs2.0 incorporates two process-based representations of inefficient drainage dominant in the literature (linked cavity and poro-elastic) and a non-mass-conserving zero-dimensional form (herein termed leaky bucket) coupled to an ice sheet systems model (the Glacial Systems Model, GSM). The linked-cavity and poro-elastic configurations include an efficient drainage scheme while the leaky bucket does not. All three systems have a positive feedback on ice velocity, whereby faster basal velocities increase melt supply. The poro-elastic and leaky-bucket systems have diagnostic effective pressure relationships – only the linked-cavity system has an additional negative feedback, whereby faster basal ice velocities increase the dynamical effective pressure due to higher cavity opening rates. We examine the contribution of mass transport, efficient drainage, and the linked-cavity negative feedback to surging. We also assess the likely bounds on poorly constrained subglacial hydrology parameters and adopt an ensemble approach to study their impact and interactions within those bounds. We find that subglacial hydrology is an important system inductance for realistic ice stream surging but that the three formulations all exhibit similar surge behaviour within parametric uncertainties. Even a detail as fundamental as mass-conserving transport of subglacial water is not necessary for simulating a full range of surge frequency and amplitude. However, one difference is apparent: the combined positive and negative feedbacks of the linked-cavity system yield longer duration surges and a broader range of effective pressures than its poro-elastic and leaky-bucket counterparts. [ABSTRACT FROM AUTHOR]

The Antarctic Ice Sheet represents the largest source of uncertainty in future sea level rise projections, with a contribution to sea level by 2100 ranging from -5 to 43 cm of sea level equivalent under high carbon emission scenarios estimated by the recent Ice Sheet Model Intercomparison for CMIP6 (ISMIP6). ISMIP6 highlighted the different behaviors of the East and West Antarctic ice sheets, as well as the possible role of increased surface mass balance in offsetting the dynamic ice loss in response to changing oceanic conditions in ice shelf cavities. However, the detailed contribution of individual glaciers, as well as the partitioning of uncertainty associated with this ensemble, have not yet been investigated. Here, we analyze the ISMIP6 results for high carbon emission scenarios, focusing on key glaciers around the Antarctic Ice Sheet, and we quantify their projected dynamic mass loss, defined here as mass loss through increased ice discharge into the ocean in response to changing oceanic conditions. We highlight glaciers contributing the most to sea level rise, as well as their vulnerability to changes in oceanic conditions. We then investigate the different sources of uncertainty and their relative role in projections, for the entire continent and for key individual glaciers. We show that, in addition to Thwaites and Pine Island glaciers in West Antarctica, Totten and Moscow University glaciers in East Antarctica present comparable future dynamic mass loss and high sensitivity to ice shelf basal melt. The overall uncertainty in additional dynamic mass loss in response to changing oceanic conditions, compared to a scenario with constant oceanic conditions, is dominated by the choice of ice sheet model, accounting for 52 % of the total uncertainty of the Antarctic dynamic mass loss in 2100. Its relative role for the most dynamic glaciers varies between 14 % for MacAyeal and Whillans ice streams and 56 % for Pine Island Glacier at the end of the century. The uncertainty associated with the choice of climate model increases over time and reaches 13 % of the uncertainty by 2100 for the Antarctic Ice Sheet but varies between 4 % for Thwaites Glacier and 53 % for Whillans Ice Stream. The uncertainty associated with the ice–climate interaction, which captures different treatments of oceanic forcings such as the choice of melt parameterization, its calibration, and simulated ice shelf geometries, accounts for 22 % of the uncertainty at the ice sheet scale but reaches 36 % and 39 % for Institute Ice Stream and Thwaites Glacier, respectively, by 2100. Overall, this study helps inform future research by highlighting the sectors of the ice sheet most vulnerable to oceanic warming over the 21st century and by quantifying the main sources of uncertainty. [ABSTRACT FROM AUTHOR]

The development of the Fennoscandian Ice Sheet (FIS) in its SW Baltic Sea sector during MIS 2 led to the formation of the Jasmund Glacitectonic Complex (JGC). The evolutional stages of the latter reflect small-dimensioned ice sheet oscillations. A glacitectonically formed piggyback basin is in focus to give direct information on the evolution of the JGC. We present a set of five OSL ages from syn-kinematic glacilacustrine and alluvial fan sediments deposited between ice margin-parallel thrust-bounded ridges from the southern margin of the JGC. The ice front proximal depositional environment and short transport distances of the sediments implicate a possible insufficient bleaching of the luminescence signal. Nonetheless, we were able to gain reliable ages of ~22 ka for the growth of the southern JGC. The samples pre-date the Oder ice stream reaching its maximum extent, the Pomeranian ice marginal belt at ~18-20 ka. A comparison of our data with the existing age data related to the evolution of the JGC will contribute to a better understanding of the FIS dynamics in the research area and the timing of the Pomeranian phase in NE Germany. [ABSTRACT FROM AUTHOR]

Reconstructions of the southwestern Laurentide Ice Sheet (LIS) from geomorphology have revealed complex cross‐cutting ice streams, indicative of surging and flow switching. This flow pattern and behaviour is distinct from the rest of the ice sheet where ice streams flowed radially to the ice sheet margin. Many ice streams in the southwestern sector originate around the Alberta Oil Sands (AOS). Previous reports have detected oil sands material in surficial and glacial sediments south of the AOS, demonstrating potential glacial mobilisation of the oil sands. In this study, we use geochemical fingerprinting to systematically investigate surficial sediments from the former Central Alberta Ice Stream (CAIS) flow track. We compare the geochemical signatures of 82 sediments from within and outside the CAIS limits with those of the AOS (from mines and natural exposures), using gas chromatography–mass spectrometry oil–oil correlation techniques. Our results provide geochemical evidence of glacial erosion and long‐distance mobilisation of AOS producing contamination signatures in sediments throughout the CAIS flow track. The strength of the AOS signature is particularly strong in sediments along the terminating margins, to the east of Calgary and in the Cooking Lake area to the southeast of Edmonton. These results inform theoretical models of enhanced slipperiness, inspired by slippery liquid infused porous surfaces (SLIPS), whereby oil lubrication of the basal sediment influences the degree of sliding and basal deformation in the ice stream. We hypothesise that naturally occurring hydrocarbons at the basal interface exerted control on the location of the onset zone of the CAIS and surrounding ice streams in Alberta. The enhanced slipperiness caused by oil contamination may also explain ice streaming and surging in other ice sheets, such as the Barents Sea Ice Sheet, where hydrocarbons are known to have been driven to the sedimentary interface during glaciation. [ABSTRACT FROM AUTHOR]