Search

The mixing of water in the ocean is a key influence on our climate and ecosystems, but its importance is often under-recognized. Mixing in the seas around Antarctica—a key process that redistributes heat, carbon, nutrients, plankton, and all other things in the sea, with profound consequences—also affects the stability of the continent's glaciers and ice sheets, with consequences for sea level rise globally. A recent discovery showed that when the fronts of glaciers disintegrate, they "calve" huge chunks of ice that can cause underwater tsunamis in the ocean, which can spread for many miles and cause strong bursts of mixing when they break. This surprising finding is changing the way we think about mixing close to Antarctica, what causes it, and how it matters. [ABSTRACT FROM AUTHOR]

The unknown and not observable basal drag distribution underneath the glaciers strongly impacts ice flow speeds and, with that, the ongoing mass loss of the Antarctic Ice Sheet. Therefore, basal drag is required for precise ice sheet modeling and accurate projections of future sea‐level rise. This can be achieved by applying an inverse method based on observed ice surface velocity data. The forward model equations, including boundary conditions, represent the ice dynamics in an approximated way. The ice‐base boundary condition is the main focus here, as it describes a non‐linear friction law. This law depends on the unknown basal drag parameter determined by the inversion and utilizes an effective pressure from a subglacial hydrology model. The inverse method minimizes a cost function consisting of the sum of a term quantifying the misfit between simulated and observed surface velocities and a regularization term to penalize unrealistic oscillations in the basal drag parameter. An L‐curve analysis determines the optimal weighting of both cost function terms. Here, we perform inversions for three different domains of Antarctica, comprising about 9.4 Mio km2${\rm km}^2$, to compare the variability of the resulting basal drag and the L‐curves. The results present a low basal drag, as well as a low variability, predominating over large parts of the interior of EAIS and WAIS. In contrast, some fast‐flowing glaciers reveal a patchy pattern of alternating high and low basal drag. In addition, parts of the grounding line exhibit a high basal drag, which potentially affects the future retreat behavior of the ice sheets. [ABSTRACT FROM AUTHOR]

Beginning in March 2022, the Antarctic Peninsula's Hektoria Glacier experienced an unprecedented retreat of ∼23 km over 1.5 years, one of the fastest observed glacier retreats on record. Improving constraints on the drivers of such extreme events is key to understanding glacier change around the continent and future sea‐level rise. We use satellite remote sensing and reanalysis data to characterize changes in Hektoria, a former Larsen B Ice Shelf tributary, over the last ∼20 years and document a period of retreat from 2002 to 2011, and readvancement from 2011 to 2022. We find that the long‐term ice front and velocity response (2002–2022) correlated more strongly with changes in modeled ocean temperatures compared to surface air temperatures. However, the acute loss of buttressing support following fast ice collapse paired with a near‐contemporaneous extreme atmospheric river in the region likely catalyzed the unprecedented 2022–2023 retreat. Plain Language Summary: The Antarctic Ice Sheet is one of the largest sources for future sea level rise, yet how much and how fast ice is lost to the ocean here remains relatively unknown. Ice shelves can buttress glaciers from flowing quickly into the ocean, stabilizing their movement and limiting mass discharge. As ice shelves retreat or break up, glaciers accelerate, adding mass to the ocean. In this study, we use imagery and elevation data collected from airborne studies and satellites to characterize how Hektoria Glacier—a marine‐terminating glacier located on the Eastern Antarctic Peninsula that was previously a Larsen B Ice Shelf tributary—has changed over the past 20 years. We compare these changes with available ocean and air temperatures in the region to determine how they influenced the observed fluctuations over time. We find that Hektoria retreated from 2002 to 2011 and readvanced from 2011 to 2022, followed by an unprecedented retreat of ∼23 km between March 2022 and August 2023. We find that abrupt changes in stress following buttressing loss drives glacier change, while modeled ocean temperatures wield influence on Hektoria's long‐term fluctuations and atmospheric temperatures drive shorter term changes in glacier response. Key Points: Hektoria Glacier retreated ∼23 km between March 2022 and August 2023—one of the fastest observed marine‐terminating glacier retreatsChanges in buttressing support and mid‐depth ocean temperatures served as primary drivers for change at Hektoria between 2002 and 2022Understanding long‐ and short‐term glacier response to ocean and atmospheric variability is key to improved sea level rise predictions [ABSTRACT FROM AUTHOR]

The surface elevation of the Greenland Ice Sheet is constantly changing due to the interplay between surface mass balance processes and ice dynamics, each exhibiting distinct spatiotemporal patterns. Here, we employ satellite and airborne altimetry data with fine spatial (1 km) and temporal (monthly) resolutions to document this spatiotemporal evolution from January 2003 to August 2023. To estimate elevation changes of the Greenland Ice Sheet (GIS), we utilize radar altimetry data from CryoSat-2 and EnviSat, laser altimetry data from the ICESat and ICESat-2, and laser altimetry data from NASA's Operation IceBridge Airborne Topographic Mapper. We produce continuous monthly ice surface elevation changes from January 2003 to August 2023 on a 1 km grid covering the entire GIS. We estimate cumulative ice loss of 4,352 Gt ± 315 Gt (12.1 ± 0.9 mm sea level equivalent) during this period, excluding peripheral glaciers. Between 2003 and 2023, the ice sheet land-terminating margin underwent a significant cumulative thinning of several meters. Ocean-terminating glaciers exhibited thinning between 20–40 m, with Jakobshavn Isbræ experiencing an exceptional thinning of nearly 70 m. This dataset of fine-resolution altimetry data in both space and time will support studies of ice mass loss and useful for GIS ice sheet modelling. To validate our monthly mass changes of the Greenland ice sheet, we use mass change from satellite gravimetry and mass change from the Input-Output method. On multiannual timescales, there is a strong correlation between the time series, with R values ranging from 0.88 to 0.92. [ABSTRACT FROM AUTHOR]

The Greenland Ice Sheet is highly sensitive to climate change, leading to significant retreat along its edges. This rapid ice loss contributes to rising sea levels and impacts the Earth's climate stability. Understanding the extent of recent glacier retreat is crucial in order to determine if it is unprecedented or within ranges of natural variability. Palaeoenvironmental studies aim to identify past glacial phases and landscape changes using advanced dating methods such as cosmic ray exposure (CRE) dating. In NE Greenland, CRE dating has helped establish the timing of glacial oscillations, yet a comprehensive understanding of glacial fluctuations during specific periods still needs to be developed. This study aims to chronologically constrain the postglacial landscape evolution of two NE Greenland valleys from the Young Sund–Tyrolerfjord area (74°N, 20–25°E) from the onset of deglaciation and throughout the Holocene to better understand glacial and postglacial changes. The chronological framework relies on 27 10Be cosmic‐ray exposure ages that constrain our interpretation of the geomorphological features in both valleys. Inconsistencies were observed in the ages dataset, highlighting potential bias associated with nuclide inheritance and post‐glacial dynamics. Despite limitations, the CRE results confirm the general pattern observed in NE Greenland: (i) major deglaciation and disconnection of glaciers from the main glacial systems during the Lateglacial and Early Holocene with a rapid but not homogeneous deglaciation within the range from ~14.3 to 11.9 ka; (ii) no evidence of glacial activity during the Middle Holocene, probably associated with the withdrawn position of the ice masses’ fronts; and (iii) glacier expansion during the Late Holocene, with a Little Ice Age advance as the last significant period of glacial regrowth. [ABSTRACT FROM AUTHOR]

Two-thirds of all glaciers worldwide are projected to disappear by 2100 CE. Large uncertainties however remain in maritime settings, where some glaciers have recently gained mass in response to increased snowfall. One of these regions is southern Patagonia, where increased precipitation since the 1980s seems to have attenuated glacier retreat. Whether this exceptional behavior will continue in a warmer future is unclear. Here, we use a numerical ice-flow model constrained by paleoglaciological data to simulate how climate variability influenced the evolution of three maritime outlet glaciers of the Southern Patagonian Icefield during the last 6000 years. Our experiments suggest that precipitation drove 67% of the centennial-scale fluctuations in the volume of the modeled glaciers. When applied to the temperature projected by 2100 CE, our simulations show that precipitation needs to increase by 10–50% to maintain present-day glacier volumes, depending on the climate scenario (SSP1-2.6 to SSP5-8.5). This implies that if greenhouse-gas emission cuts fail, these glaciers will enter a warmer regime unseen over the last 6000 years, where precipitation cannot offset glacier mass loss. Conversely, if emissions are curtailed, increased precipitation could halt mass loss of some of Patagonia's largest glaciers, and potentially of other maritime glaciers worldwide. [ABSTRACT FROM AUTHOR]

Antarctica's Lambert Glacier drains about one-sixth of the ice from the East Antarctic Ice Sheet and is considered stable due to the strong buttressing provided by the Amery Ice Shelf. While previous projections of the sea-level contribution from this sector of the ice sheet have predicted significant mass loss only with near-complete removal of the ice shelf, the ocean warming necessary for this was deemed unlikely. Recent climate projections through 2300 indicate that sufficient ocean warming is a distinct possibility after 2100. This work explores the impact of parametric uncertainty on projections of the response of the Lambert–Amery system (hereafter "the Amery sector") to abrupt ocean warming through Bayesian calibration of a perturbed-parameter ice-sheet model ensemble. We address the computational cost of uncertainty quantification for ice-sheet model projections via statistical emulation, which employs surrogate models for fast and inexpensive parameter space exploration while retaining critical features of the high-fidelity simulations. To this end, we build Gaussian process (GP) emulators from simulations of the Amery sector at a medium resolution (4–20 km mesh) using the Model for Prediction Across Scales (MPAS)-Albany Land Ice (MALI) model. We consider six input parameters that control basal friction, ice stiffness, calving, and ice-shelf basal melting. From these, we generate 200 perturbed input parameter initializations using space filling Sobol sampling. For our end-to-end probabilistic modeling workflow, we first train emulators on the simulation ensemble and then calibrate the input parameters using observations of the mass balance, grounding line movement, and calving front movement with priors assigned via expert knowledge. Next, we use MALI to project a subset of simulations to 2300 using ocean and atmosphere forcings from a climate model for both low- and high-greenhouse-gas-emission scenarios. From these simulation outputs, we build multivariate emulators by combining GP regression with principal component dimension reduction to emulate multivariate sea-level contribution time series data from the MALI simulations. We then use these emulators to propagate uncertainty from model input parameters to predictions of glacier mass loss through 2300, demonstrating that the calibrated posterior distributions have both greater mass loss and reduced variance compared to the uncalibrated prior distributions. Parametric uncertainty is large enough through about 2130 that the two projections under different emission scenarios are indistinguishable from one another. However, after rapid ocean warming in the first half of the 22nd century, the projections become statistically distinct within decades. Overall, this study demonstrates an efficient Bayesian calibration and uncertainty propagation workflow for ice-sheet model projections and identifies the potential for large sea-level rise contributions from the Amery sector of the Antarctic Ice Sheet after 2100 under high-greenhouse-gas-emission scenarios. [ABSTRACT FROM AUTHOR]

Current climate warming is accelerating mass loss from glaciers and ice sheets. In Greenland, the rates of mass changes are now dominated by changes in surface mass balance (SMB) due to increased surface melting. To improve the future sea-level rise projections, it is therefore critical to have an accurate estimate of the SMB, which depends on the representation of the processes occurring within the snowpack. The Explicit Snow (ES) scheme implemented in the land surface model Organising Carbon and Hydrology In Dynamic Ecosystems (ORCHIDEE) has not yet been adapted to ice-covered areas. Here, we present the preliminary developments we made to apply the ES model to glaciers and ice sheets. Our analysis mainly concerns the model's ability to represent ablation-related processes. At the regional scale, our results are compared to the MAR regional atmospheric model outputs and to MODIS albedo retrievals. Using different albedo parameterizations, we performed offline ES simulations forced by the MAR model over the 2000–2019 period. Our results reveal a strong sensitivity of the modelled SMB components to the albedo parameterization. Results inferred with albedo parameters obtained using a manual tuning approach present very good agreement with the MAR outputs. Conversely, with the albedo parameterization used in the standard ORCHIDEE version, runoff and sublimation were underestimated. We also tested parameters found in a previous data assimilation experiment, calibrating the ablation processes using MODIS snow albedo. While these parameters greatly improve the modelled albedo over the entire ice sheet, they degrade the other model outputs compared to those obtained with the manually tuned approach. This is likely due to the model overfitting to the calibration albedo dataset without any constraint applied to the other processes controlling the state of the snowpack. This underlines the need to perform a "multi-objective" optimization using auxiliary observations related to internal snowpack processes. Although there is still room for further improvements, the developments reported in the present study constitute an important advance in assessing the Greenland SMB with possible extension to mountain glaciers or the Antarctic ice sheet. [ABSTRACT FROM AUTHOR]

In recent decades, glaciers in the Amundsen Sea Embayment in West Antarctica have made the largest contribution to mass loss from the entire Antarctic Ice Sheet. Glacier retreat and acceleration have led to concerns about the stability of the region and the effects of future climate change. Coastal thinning and near-synchronous increases in ice flux across neighbouring glaciers suggest that ocean-driven melting is one of the main drivers of mass imbalance. However, the response of individual glaciers to changes in ocean conditions varies according to their local geometry. One of the largest and fastest-flowing of these glaciers, Pine Island Glacier (PIG), underwent a retreat from a subglacial ridge in the 1940s following a period of unusually warm conditions. Despite subsequent cooler periods, the glacier failed to recover back to the ridge and continued retreating to its present-day position. Here, we use the ice-flow model Úa to investigate the sensitivity of this retreat to changes in basal melting. We show that a short period of increased basal melt was sufficient to force the glacier from its stable position on the ridge and undergo an irreversible retreat to the next topographic high. Once high melting begins upstream of the ridge, only near-zero melt rates can stop the retreat, indicating a possible hysteresis in the system. Our results suggest that unstable and irreversible responses to warm anomalies are possible and can lead to substantial changes in ice flux over relatively short periods of only a few decades. [ABSTRACT FROM AUTHOR]

Glaciers and ice caps surrounding the Greenland Ice Sheet are found to be sensitive to warming climate thus the knowledge of their thickness and internal structure is substantial to determine their future impact on sea level and local environment. Still, in situ glaciological measurements of such glaciers are very scarce. Here, we present the results of ground penetrating radar (GPR) and uncrewed aerial vehicle surveys conducted on the two southern outlet glaciers of Qaanaaq Ice Cap in NW Greenland. GPR measurements reveal up to 170 m thick ice and the lack of englacial hyperbolae indicating no developed en/subglacial drainage system. The glaciers consist mainly of radar transparent facies characteristic for cold ice, while limited scattering facies appear closer to the glacier's terminus beneath the thinnest ice and are attributed to debris inside the ice. Results show that the glaciers flow into narrow V‐shaped valleys suggesting spatio‐temporally limited subglacial erosion and restricted possible distribution of temperate ice in the past. The comparison of the ice thickness measurement data with global ice thickness model estimates shows considerable discrepancies emphasising the need of modelling improvements in the case of narrow valley and outlet glaciers. [ABSTRACT FROM AUTHOR]

Sediment erosion, transport, and deposition by glaciers and ice sheets play crucial roles in shaping landscapes, provide important nutrients to downstream ecosystems, and preserve key indicators of past climate conditions in the geologic record. While previous work has quantified sediment fluxes from subglacial meltwater, we also observe sediment entrained within basal ice, transported by the flow of the glacier itself. However, the formation and evolution of these debris‐rich ice layers remains poorly understood and rarely represented in landscape evolution models. Here, we identify a characteristic sequence of basal ice layers at Mendenhall Glacier, Alaska. We develop a numerical model of frozen fringe and regelation processes that describes the co‐evolution of this sequence and explore the sensitivity of the model to key properties of the subglacial sedimentary system, using the Instructed Glacier Model to parameterize ice dynamics. Then, we run numerical simulations over the spatial extent of Mendenhall Glacier, showing that the sediment transport model can predict the observed basal ice stratigraphy at the glacier's terminus. From the model results, we estimate basal ice layers transport between 23,300 m3 ${\mathrm{m}}^{3}$a−1 ${\mathrm{a}}^{-1}$ and 39,800 m3 ${\mathrm{m}}^{3}$a−1 ${\mathrm{a}}^{-1}$ of sediment, mostly entrained in the lowermost ice layers nearest to the bed, maximized by high effective pressures and slow, convergent flow fields. Overall, our results highlight the role of basal sediment entrainment in delivering eroded material to the glacier terminus and indicate that this process should not be ignored in broader models of landscape evolution. Plain Language Summary: Glaciers shape landscapes by eroding mountain ranges and moving the eroded sediment downstream to the lakes, fjords, and other ecosystems near the terminus. To quantify this supply of sediment, and thus understand the evolution of high alpine landscapes, we often use numerical models to simulate the processes of erosion and transport. Previous work has focused on the role of meltwater channels underneath glaciers to transport sediment, analogous to how rivers move sediment across much of Earth's surface, but we also know from field observations that a portion of the underlying sediment is picked up and incorporated into the ice itself. Here, we develop a new model to simulate where and how sediment is incorporated into ice beneath the glacier, and how those debris‐rich ice layers are modified as the glacier flows downstream. We use this model to predict the total amount of sediment reaching a glacier's snout, as well as the spatial distribution and thickness of the debris‐rich ice. Our results show that the amount of sediment transported in glacier ice likely represents a significant component of the total amount of sediment eroded or transported by the glacier, and thus should not be ignored in broader models of landscape evolution. Key Points: We develop a model of sediment entrainment beneath an alpine glacier to predict sediment fluxes and basal ice stratigraphy at the terminusThe model explains the co‐occurrence of solid, debris‐rich ice layers and overlying dispersed basal ice facies observed at glacier marginsAt a case study, we find sediment can be entrained up to 5 m in basal ice, maximized at high effective pressure and slow, convergent flow [ABSTRACT FROM AUTHOR]

The flow of glaciers and ice sheets is largely influenced by friction at the ice‐bed interface that can trigger rapid changes in glacier motion ranging from seasonal velocity variations to cyclic surge instabilities or even devastating glacier collapse. This wide range of transient glacier dynamics is currently not captured by models, and its implications for long‐term glacier evolution are uncertain. This highlights the need of developing improved descriptions for processes that occur at the glacier bed. Here, we present a model that describes the evolution of basal friction inspired by a "rate and state" approach, coupled to models of subglacial drainage and glacier flow, and investigate how these couplings affect the dynamics of glaciers. We show that a wide range of sliding behavior results from a feedback loop between subglacial drainage efficiency and friction which depends on the evolution of a frictional state that can be interpreted as the degree of cavitation or till porosity for hard and soft beds, respectively. In our simulations, we find that glaciers are susceptible to surging if they exhibit a transition to velocity weakening friction associated with a poor sensitivity of the drainage capacity to the frictional state. This potential materializes if the local topography and mass balance create the conditions for high water pressure to build up in an area sufficiently large to exceed a critical length. We advocate accounting for feedback loops between friction and drainage as a promising avenue for better understanding dynamical instabilities of glaciers and ice sheets. Plain Language Summary: Glaciers move through a combination of basal sliding and viscous deformation of the ice. Observed velocities range widely from glaciers creeping at a few meters per year, to continuously fast‐flowing ice‐streams at velocities of kilometers per year, all the way to devastating ice avalanches. Some glaciers can alter their velocity in quasi‐periodic intervals, known as surge‐type glaciers. These glaciers can move steadily for decades, then speed up for a short period of time by up to two orders of magnitude or more, before returning to quiescence. In this study, we show that a wide range of glacier behavior can be understood from the interaction between friction and the subglacial drainage system. Glacier surges are possible if the efficiency of the drainage system is not high enough to drain the water that enters and is produced at the glacier base. This allows water pressure to increase over time, which can potentially trigger a frictional instability. Depending on the configuration of the subglacial drainage and how it evolves in response to sliding, glaciers can be stable or surge‐type, in addition to exhibiting several different classes of velocity variations. Key Points: We present a model combining rate and state friction with subglacial hydrology processes that spontaneously generates glacier surgesA wide range of glacier behavior emerges from the model, including periodic and quasi‐periodic glacier surgesSurges arise if the drainage system is unable to sufficiently drain water from the glacier base, triggering a frictional instability [ABSTRACT FROM AUTHOR]

Glaciers and ice sheets are an integral part of Earth's system, advancing and retreating in response to changes in climate. Clues about the past, present and future behaviour of these ice masses are found throughout current and former glaciated landscapes. In this commentary, we outline recent scientific advances from a collection of articles in which geomorphological evidence is used to inform us about the behaviour of glaciers and ice sheets across a range of spatial (landform to continent) and temporal (seasons to millennia) scales. Through a diversity of approaches including field measurements, remote sensing and numerical modelling, these studies build on an extensive background literature to deepen our understanding of how ice flows, how glaciers and ice sheets respond to climate change, and of the processes of ice advance and retreat and the stability of the system. Further integration of knowledge across the fields of geomorphology and glaciology will have tangible benefits for managing the societal and environmental impacts of glacier change and for improved projections of sea‐level rise over the coming decades to centuries. [ABSTRACT FROM AUTHOR]

The negative surface mass balance of glaciers and ice caps under a warming climate impacts local ecosystems, influencing the volume and timing of water flow in local catchments while also contributing to global sea level rise. Peripheral glaciers distinct to the Greenland ice sheet respond faster to climate change than the main ice sheet. Accurate assessment of surface mass balance depends on in situ observations of near-surface climate and ice ablation, but very few in situ observations of near-surface climate and ice ablation are freely available for Greenland's peripheral glaciers. The transect of three automated weather and ablation stations on the peripheral A. P. Olsen Ice Cap in northeast Greenland is an example of these much needed data. The transect has been monitored since 2008, and in 2022, the old weather and ablation stations were replaced by a new standardized setup. In order to ensure comparable data quality of the old and new monitoring station setups, it is necessary to re-evaluate the data collected between 2008 and 2022. This paper presents the fully reprocessed near-surface climate and ablation data from the A. P. Olsen Ice Cap transect from 2008 to 2022, with a focus on data quality and the usability in ice ablation process studies. The usability of the data is exemplified by the data in an energy balance melt model for two different years. We show that the inherent uncertainties in the data result in an accurate reproduction of ice ablation for just one of the two years. A transect of three automatic ablation and weather stations of this length is unique to Greenland's peripheral glaciers, and it has a broad scale of usage from input to climate reanalysis and detailed surface ablation studies. The dataset can be downloaded at 10.22008/FK2/X9X9GN. [ABSTRACT FROM AUTHOR]

Greenland's marine‐terminating glaciers connect the ice sheet to the ocean and provide a critical boundary where heat, freshwater, and nutrient exchanges take place. Buoyant freshwater runoff from inland ice sheet melt is discharged at the base of marine‐terminating glaciers, forming vigorous upwelling plumes. It is understood that subglacial plumes modify waters near glacier fronts and increase submarine glacier melt by entraining warm ambient waters at depth. However, ocean observations along Greenland's coastal margins remain biased toward summer months which limits accurate estimation of ocean forcing on glacier retreat and acceleration. Here, we fill a key observational gap in northwest Greenland by describing seasonal hydrographic variation at glacier fronts in Melville Bay using in situ observations from moorings deployed year‐round, CTDs, and profiling floats. We evaluated local and remote forcing using remote sensing and reanalysis data products alongside a high‐resolution ocean model. Analysis of the year‐round hydrographic data revealed consistent above‐sill seasonality in temperature and salinity. The warmest, saltiest waters occurred in spring (April–May) and primed glaciers for enhanced submarine melt in summer when meltwater plumes entrain deep waters. Waters were coldest and freshest in early winter (November–December) after summer melt from sea ice, glacier ice, and icebergs provided cold freshwater along the shelf. Ocean variability was greatest in the summer and fall, coincident with increased freshwater runoff and large wind events before winter sea ice formation. Results increase our mechanistic understanding of Greenland ice‐ocean interactions and enable improvements in ocean model parameterization. Plain Language Summary: Many of Greenland's glaciers terminate in the ocean and form an important boundary between coastal waters and the ice sheet. During summer months, meltwater from the ice sheet flows out below glaciers that terminate in the ocean. The less‐dense meltwater rises to the surface as a plume, increasing near‐glacier ocean mixing and submarine glacier melt. Understanding processes at the glacier‐ocean interface, including how plumes modify nearshore waters, requires sustained observations near glacier termini. However, ocean measurements are largely taken during the summer and do not necessarily represent the total variability observed in the system. Here, we fill an important knowledge gap by presenting year‐round measurements of temperature and salinity near glacier fronts in Northwest Greenland. We found the warmest, saltiest waters in the spring and coldest, freshest waters in early winter at all sites. Warm waters in the spring prime deep glaciers for submarine melt during summer. Increased winds in the fall initiated vertical mixing in nearshore waters before winter sea ice created a surface barrier, stabilizing the water column. Results advance our understanding of the ocean's role in the acceleration and retreat of Greenland's glaciers, which is important for estimating large‐scale ocean circulation and Greenland's ice sheet mass loss. Key Points: Moored hydrographic observations at glacier fronts in NW Greenland revealed consistent above‐sill seasonality in temperature and salinityOcean temperatures reached a maximum in spring, priming deep glaciers for submarine melt during summerRemote forcing supplied the renewal of warm, salty water in the spring and local forcing controlled increased variability in the summer/fall [ABSTRACT FROM AUTHOR]

Marine‐terminating glacier fjords play a central role in the transport of oceanic heat toward ice sheets, regulating their melt. Mixing processes near glacial termini are key to this circulation but remain poorly understood. We present new summer measurements of circulation and mixing near a marine‐terminating glacier with active sub‐glacial discharge. 65% of the fjord's vertical overturning circulation is driven by the buoyant plume, however we newly report intense vertical and horizontal mixing in the plume's horizontal spreading phase, accounting for the remaining 35%. Buoyant plume theory supports 2%–5% of total glacial melt. Thus, most of the heat associated with vertical overturing short‐circuits the glacial front. We find however that turbulence in the horizontal spreading phase redistributes the short‐circuited heat back into the surface waters of the near‐glacial zone. Our findings highlight the need for further research on the complex mixing processes that occur near the glacier terminus. Plain Language Summary: Melting of glacial ice is the single largest contributor to global sea‐level rise. Many glaciers flow into the ocean where the near‐vertical ice wall is bathed in relatively warm sea water. Freshwater from ice surface melting seeps down through cracks and crevasses in the glacier to be discharged at the base of the ice wall, many tens to hundreds of meters below the sea surface. This fresh water rises from the depths as a highly turbulent plume, drawing in and pushing upwards the surrounding seawater, eventually spreading horizontally as a mixture of fresh discharge and mixed‐in seawater. Few measurements exist in the dangerous zone where iceberg often calve. Using data from a robotic platform we show that the vertical rise and the horizontal spreading of fresh water both play important roles in the total ocean‐induced melting of the glacial face in this type of system. Key Points: Entrainment into the buoyant plume drives ∼65% of fjord vertical overturning circulationIntense turbulent mixing in the horizontal spreading phase of the plume drives the remaining ∼35%95%–98% of the heat in the rising plume short‐circuits the glacier, to be redistributed into glacial proximal waters by horizontal mixing [ABSTRACT FROM AUTHOR]

The collapse of ice shelves could expose tall ice cliffs at ice sheet margins. The marine ice cliff instability (MICI) is a hypothesis that predicts that, if these cliffs are tall enough, ice may fail structurally leading to self-sustained retreat. To date, projections that include MICI have been performed with a single model based on a simple parameterization. Here, we implement a physically motivated parameterization in three ice sheet models and simulate the response of the Amundsen Sea Embayment after a hypothetical collapse of floating ice. All models show that Thwaites Glacier would not retreat further in the 21st century. In another set of simulations, we force the grounding line to retreat into Thwaites' deeper basin to expose a taller cliff. In these simulations, rapid thinning and velocity increase reduce the calving rate, stabilizing the cliff. These experiments show that Thwaites may be less vulnerable to MICI than previously thought, and model projections that include this process should be re-evaluated. [ABSTRACT FROM AUTHOR]

Numerical modeling of ice sheet dynamics is a critical tool for projecting future sea level rise. Among all the processes responsible for the loss of mass of the ice sheets, enhanced ice discharge triggered by the retreat of marine-terminating glaciers is one of the key drivers. Numerical models of ice sheet flow are therefore required to include ice front migration in order to reproduce today's mass loss and to be able to predict their future. However, the discontinuous nature of calving poses a significant numerical challenge for accurately capturing the motion of the ice front. In this study, we explore different stabilization techniques combined with varying reinitialization strategies to enhance the numerical stability and accuracy of solving the level-set function, which tracks the position of the ice front. Through rigorous testing on an idealized domain with a semicircular and a straight-line ice front, including scenarios with diverse front velocities, we assess the performance of these techniques. The findings contribute to advancing our ability to model ice sheet dynamics, specifically calving processes, and provide valuable insights into the most effective strategies for simulating and tracking the motion of the ice front. [ABSTRACT FROM AUTHOR]

Terminus change is a complex outcome of ice‐ocean boundary processes and poses challenges for ice sheet models due to inadequate calving laws, creating uncertainty in sea level change projections. To address this, we quantify glacier termini sinuosity and convexity, testing the hypothesis that terminus morphology reflects dominant calving processes. Using 10 glaciers with diverse calving styles in Greenland over the period from 1985 to 2021, we establish a supervised classification of calving style by comparing morphology and literature‐derived calving observations. Validation with four of these glaciers and flotation conditions and subglacial discharge routing observations confirms concave, smooth termini indicate buoyant flexure dominated‐calving, while convex, sinuous termini suggest serac failure dominated‐calving. We also identify a mixed style where both calving types may occur. We use these classes to label calving style from 1985 to 2021 for all 10 glaciers and explore how this changes over time as glaciers retreat. Plain Language Summary: The ice‐ocean boundary (terminus) of the Greenland ice sheet changes over space and time due to the melt and calving of ice. There is a range of calving behaviors due to many environmental and geometric factors, including changing ocean temperatures, glacier melting, and the geometry of the glacier. This makes representing how these boundary changes in ice sheet models are difficult, leading to uncertainties in predicting sea level rise from retreating glaciers. The study focuses on the shape (morphology) of 10 glaciers from 1985 to 2021, finding that smooth, concave termini indicate a type of calving related to the flotation of the glacier, while sinuous, convex termini indicate melt‐dominated calving. The research helps classify calving styles and understand how they change over time as glaciers retreat. Key Points: Two calving styles are linked to distinct terminus morphologies but both can exist for a single glacier (over time or space)Mixed styles of calving are more difficult to classifyClassifying terminus morphology over time can aid in identifying how terminus calving mechanisms change over time [ABSTRACT FROM AUTHOR]

Glaciers, distinct from the Greenland and Antarctic ice sheets, play a crucial role in Earth's climate system by affecting global sea levels, freshwater availability, nutrient and energy budgets and regional climate patterns. Accurate measurements of glacier mass changes are needed to understand and project glacier evolution and its related impacts on the climate system. Two distinct methods allow to measure glacier mass changes at high spatial resolution. Remotely sensed surface elevation data provides volume change estimates over large glacierized regions for multi-annual to decadal time periods. Field glaciological measurements provide annually to seasonally resolved information on glacier mass change for a small sample of the world's glaciers. By combining the two methods we provide annual time series of individual glacier mass changes and related uncertainties spanning the hydrological years from 1976 to 2023. The per-glacier time series can then be seamlessly integrated into annually resolved global regular grids of glacier mass changes at user-specified spatial resolution. Our results undergo a leave-one-out cross-validation confirming uncertainty estimates at the glacier level to be in the conservative side. Our dataset provides a new baseline for future glacier change modelling assessments and their impact on the world's energy, water, and sea-level budget. The present annual mass change time-series for the individual glaciers and the derived global gridded annual mass change product at a spatial resolution of 0.5° latitude and longitude will be made available from the WGMS webpage. During the review process, the dataset is temporarily available from URL: https://user.geo.uzh.ch/idussa/Dussaillant_etal_ESSD_data/. [ABSTRACT FROM AUTHOR]

This review of the dynamics of the flow of glaciers and ice sheets focusses on the mathematical models which have been developed to explain a number of observations in the behaviour of these large ice masses. It takes a personal view, starting with an account of the historical development of the subject, followed by a tour of some of the observational puzzles and theoretical solutions which I have found most intriguing. Most of the discourse involves basal processes: glacier sliding, subglacial hydrology, bedform generation. I describe in some detail wave ogives, surging glaciers, subglacial floods and subglacial bedforms. I variously explore issues in sliding theory: sub-temperate sliding, slip or till deformation; the rheology of till; the status of the instability theory of drumlin formation. And in an effort to broaden the field of view, I include discussion of Dansgaard-Oeschger and Heinrich events, grounding line stability, and the formation of eskers and, particularly, flutes. [ABSTRACT FROM AUTHOR]

Radio-echo sounding is a standard technique for imaging the englacial stratigraphy of glaciers and ice sheets. In most cases, internal reflection horizons (IRHs) represent former glacier surfaces, comprise information about past accumulation and ice deformation, and enable the linking of ice core chronologies. IRHs in the lower third of the ice column are often difficult to detect or coherently trace. In the polar ice sheets, progress in IRH detection has been made by using multistatic, phase-coherent radars, enabling focused synthetic aperture radar (SAR) processing. However, these radar systems are often not suitable for deployment on mountain glaciers. We present a proof-of-concept study for a lightweight, phase-coherent and ground-based radar system, based on the phase-sensitive radio-echo sounder (pRES). To improve the detectability of IRHs we additionally adapted a layer-optimized SAR processing scheme to this setup. We showcase the system capability at Colle Gnifetti, Switzerland and Italy, where specular reflections are now apparent down to the base of the glacier. Compared to previously deployed impulse radar systems, with the mobile pRES the age of the oldest continuously traceable IRH could be increased from 78±12 to 288±35 a. Corresponding reflection mechanisms for this glacier are linked to stratified acidic impurities which in the upper part were deposited at a higher rate due to increased industrial activity in the area. Possible improvements to the system are discussed. If successfully implemented, these may provide a new way to map the deep internal structure of Colle Gnifetti and other mountain glaciers more extensively in future deployments. [ABSTRACT FROM AUTHOR]

In this paper, we present a reconstruction of climate conditions during the Last Glacial Maximum on a karst plateau Snežnik, which lies in Dinaric Mountains (southern Slovenia) and bears evidence of glaciation. The reconstruction merges geomorphological ice limits, classified as either clear or unclear, and a computer modelling approach based on the Parallel Ice Sheet Model (PISM). Based on extensive numerical experiments where we studied the agreements between simulated and geomorphological ice extent, we propose using a combination of a high-resolution precipitation model that accounts for orographic precipitation combined with a simple elevation-based temperature model. The geomorphological ice extent can be simulated with climate to be around 6 °C colder than the modern day and with a lower-than-modern-day amount of precipitation, which matches other state-of-the art climate reconstructions for the era. The results indicate that an orographic precipitation model is essential for the accurate simulation of the study area, with moist southern winds from the nearby Adriatic Sea having a predominant effect on the precipitation patterns. Finally, this study shows that transforming climate conditions towards wetter and warmer or drier and colder does not significantly change the conditions for glacier formation. [ABSTRACT FROM AUTHOR]

We construct a high‐resolution shear‐wave velocity (VS) model for the uppermost 100 m using ambient noise tomography near the West Antarctic Ice Sheet Divide camp. This is achieved via joint inversion of Rayleigh wave phase velocity and H/V ratio, whose signal‐to‐noise ratios are boosted by three‐station interferometry and phase‐matched filtering, respectively. The VS shows a steep increase (0.04–0.9 km/s) in the top 5 m, with sharp interfaces at ∼8–12 m, followed by a gradual increase (1.2–1.8 km/s) between 10 and 45 m depth, and to 2 km/s at ∼65 m. The compressional‐wave velocity and empirically‐obtained density profile compares well with the results from Herglotz–Wiechert inversion of diving waves in active‐source shot experiments and ice core analysis. Our approach offers a tool to characterize high‐resolution properties of the firn and shallow ice column, which helps to infer the physical properties of deeper ice sheets, thereby contributes to improved understanding of Earth's cryosphere. Plain Language Summary: Accurately determining the physical properties of Antarctic ice sheets enhances our understanding of climate change impacts on the Earth's cryosphere and improves sea level rise predictions. The uppermost ∼100 m of ice sheets are often covered by firn, an intermediate material between fresh granular snow and glacier ice. Firn acts as a buffer between the atmosphere and deep ice, playing a critical role in the mass balance, flow, and dynamics of the ice sheet. In this study, we use advanced techniques to measure seismic velocities of the top ∼100 m of the ice sheet and thus characterize the properties of the firn layer. Using a week‐long record of seismic ambient noise, we obtain velocity models that reveal the detailed structure of the firn layer, including rapid velocity increases over depth that indicate rapid changes in compaction and slight lateral velocity variations. These models were validated through favorable comparisons to conventional analyses and ice core records, and can be used to estimate other important ice properties (e.g., density) via empirical relationships. Our approach could easily be up‐scaled to long‐term monitoring of larger areas, to help inform predictive models of ice sheet dynamics. Key Points: We denoise fundamental‐mode Rayleigh waves extracted from ambient noise recorded in West Antarctic via three‐station interferometryWe improve the quality of Rayleigh wave H/V ratio measurements with phase‐matched filteringThe high‐resolution Vs model obtained by jointly inverting phase velocities and H/V ratios reveals high‐resolution firn structures [ABSTRACT FROM AUTHOR]

Understanding grounding line dynamics is critical for projecting glacier evolution and sea level rise. Observations from satellite radar interferometry reveal rapid grounding line migration forced by oceanic tides that are several kilometers larger than predicted by hydrostatic equilibrium, indicating the transition from grounded to floating ice is more complex than previously thought. Recent studies suggest seawater intrusion beneath grounded ice may play a role in driving rapid ice loss. Here, we investigate its impact on the evolution of Petermann Glacier, Greenland, using an ice sheet model. We compare model results with observed changes in grounding line position, velocity, and ice elevation between 2010 and 2022. We match the observed retreat, speed up, and thinning using 3‐km‐long seawater intrusion that drive peak ice melt rates of 50 m/yr; but we cannot obtain the same agreement without seawater intrusion. Including seawater intrusion in glacier modeling will increase the sensitivity to ocean warming. Plain Language Summary: Relatively warm seawater melts marine‐terminating glaciers from below. Recent observations suggest that seawater flows below grounded ice at high tide. The presence of seawater at this boundary, referred to as seawater intrusion, has the potential to increase glacier mass loss. We test this hypothesis on Petermann Glacier, Greenland, using an ice sheet flow model. We run the model to reconstruct the glacier's behavior from 2010 to 2022 with and without seawater intrusion. We compare the results with satellite observations of velocity, grounding line position, and ice thinning. When we use enhanced ice melt rates from kilometer‐scale seawater intrusion, we match the observed retreat, speed up, and thinning. When we do not, the model fails to replicate the observations. Seawater intrusion may play a critical role in glacier evolution. Adding this process to ice flow models will increase their sensitivity to ocean warming and projections of ice mass loss and sea level rise. Key Points: Ice melt caused by seawater intrusion in the grounding zone explains the observed grounding line retreat of Petermann GlacierWithout seawater intrusion in the grounding zone, we do not replicate the full extent of observed retreatIncluding seawater intrusion in the grounding zone increases glacier mass loss [ABSTRACT FROM AUTHOR]

Piedmont glaciers (lobes), typically found in high latitudes and large mountainous regions, extend from ice sheets and ice caps to lower altitudes. However, they can also occur, although less commonly, on mid‐latitude mountains. When these fan‐like glaciers retreat, they leave behind hummocky moraines scattered in a chaotic pattern. In this study, we have mapped one of these mid‐latitude sites and established a Terrestrial cosmogenic nuclide (TCN) glacial chronology on Mount Davraz, namely Davraz hummocky moraine field (37°46′00″N, 30°43′15″E). Our findings indicate that the glaciers in this area started receding from the early local Last Glacial Maximum (LGM) period (21.8 ± 2.4 ka) to the early Late‐glacial period (17.7 ± 2.2 ka), and eventually disappearing. The deglaciation of the Mt. Davraz palaeoglacier matches nearby mountains, supported by southerly winds as significant for regional glaciation. Our discoveries reveal a robust connection between southerly winds and nearby glaciation, contributing to our understanding of how climate influences glaciers. Likewise, the glacial timelines of the neighbouring mountains align with the glacial history of Mt. Davraz. [ABSTRACT FROM AUTHOR]

Glacier flow modulates sea level and is governed largely by the viscous deformation of ice. Multiple molecular-scale mechanisms facilitate viscous deformation, but it remains unclear how each contributes to glacier-scale deformation. Here, we present a model of ice deformation that bridges laboratory and glacier scales, unifies existing estimates of the viscous parameters, and provides a framework for estimating the parameters from observations and incorporating flow laws derived from laboratory observations into glacier-flow models. Our results yield a map of the dominant deformation mechanisms in the Antarctic Ice Sheet, showing that, contrary to long-standing assumptions, dislocation creep, characterized by a value of the stress exponent n = 4, likely dominates in all fast-flowing areas. This increase from the canonical value of n = 3 dramatically alters the climate conditions under which marine ice sheets may become unstable and drive rapid rates of sea-level rise. [ABSTRACT FROM AUTHOR]

The thickness of the Antarctic ice sheet is a crucial parameter for inferring glacier mass and its evolution process. In the literature, the gravity method has been proven to be one of the effective means for estimating ice sheet thickness. And it is a preferred approach when direct measurements are not available. However, few gravity inversion methods are valid in rugged terrain areas with undulating observation surfaces (UOSs). To solve this problem, this paper proposes an improved high-precision 3D density interface inversion method considering terrain and UOSs simultaneously. The proposed method utilizes airborne gravity data at their flight altitudes, instead of the continued data yield from the unstable downward continuation procedure. In addition, based on the undulating right rectangular prism model, the large reliefs of the terrain are included in the iterative inversion. The proposed method is verified on two synthetic examples and is successfully applied to real data in East Antarctica. [ABSTRACT FROM AUTHOR]

Over the past decade, a wealth of research has been devoted to the detection of crevasses in glaciers and ice sheets via remote sensing and machine learning techniques. It is often argued that remotely sensed damage maps can function as early warning signals for shifts in ice shelf conditions from intact to damaged states and can serve as an important tool for ice sheet modellers to improve future sea level rise predictions. Here, we provide evidence for the Filchner–Ronne and Pine Island ice shelves that remotely sensed damage maps are only weakly related to the ice rate factor field A derived by an ice flow model when inverting for surface velocities. This technique is a common procedure in ice flow models, as it guarantees that any inferred changes in A relate to changes in ice flow measured through observations. The weak relationship found is improved when investigating heavily damaged shear margins, as observed on the Pine Island Ice Shelf; however, even in this setting, this association remains modest. Our findings suggest that many features identified as damage through remote sensing methods are not of direct relevance to present-day ice shelf flow. While damage can clearly play an important role in ice shelf processes and thus be relevant for ice sheet behaviour and sea level rise projections, our results imply that mapping ice damage directly from satellite observations may not directly help improve the representation of these processes in ice flow models. [ABSTRACT FROM AUTHOR]

The Greenland Ice Sheet contributed 10.6 mm to global sea level rise between 1992 and 2018, and it is projected to be the largest glacial contributor to sea level rise by 2100. Here we assess the relative importance of two major sources of uncertainty in 21st century ice loss projections: (1) the choice of sliding law and (2) the surface mass balance (SMB) forecast. Specifically, we used the ice flow model Úa to conduct an ensemble of runs for 48 combinations of sliding law and SMB forecast for three major Greenland outlet glaciers (Kangerlussuaq (KG), Humboldt (HU) and Petermann (PG) glaciers) with differing characteristics and evaluated how the sensitivity to these factors varied between the study glaciers. Overall, our results show that SMB forecasts were responsible for 4.45 mm of the variability in sea level rise by 2100 compared with 0.33 mm sea level equivalent (SLE) due to sliding law. HU had the largest absolute contribution to sea level rise and the largest range (2.16–7.96 mm SLE), followed by PG (0.84–5.42 mm SLE), and these glaciers showed similar patterns of ice loss across the SMB forecasts and sliding laws. KG had the lowest range and absolute values (- 0.60 to 3.45 mm SLE) of sea level rise, and the magnitude of mass loss by SMB forecast differed markedly between HU and PG. Our results highlight SMB forecasts as a key focus for improving estimates of Greenland's contribution to 21st century sea level rise. [ABSTRACT FROM AUTHOR]

Synectic Aperture Radar (SAR) backscatter coefficient is sensitive to glacier surface physical characteristic changes, including the states of melting and refreezing, but it is also sensitive to incidence angle variation. This study explores the capability of monitoring Greenland Ice Sheet (GrIS) melting status with Sentinel-1 dual-polarized images by referring to Automatic Weather Station (AWS) records. Sentinel-1 SAR images at five coastal regions of the GrIS are obtained from 2017 to 2021. The backscatter coefficients are normalized to an incidence angle of 30° with an empirical model. Time series of five backscatter coefficients profiles covering AWS illustrates different patterns of the ice surface dielectric constant dynamics in different elevations. The wet snow radar zone shows clear backscatter coefficients decreasing during the melting seasons, but the bare ice radar zone behaves more complexly during the melting seasons. The numbers of melting days at different elevations are also derived for each profile based on −3 dB backscatter coefficient decrease of HH and/or HV polarization, showing the heterogeneous ablation processes over the GrIS. The daily maximum 2 m air temperature on two consecutive days (before and on the SAR acquisition day) exceeds 0°C, and the daily average 2 m air temperature exceeds −0.5°C on the SAR acquisition day that was recorded by the AWS finds good agreements with the −3 dB decrease of the backscatter coefficients, suggesting the GrIS surface melting can be well captured by dual-polarized Sentinel-1 C-band SAR images. The overall agreement and Kappa coefficients are mostly better than 0.85 and 0.70, respectively, for HH images and 0.80 and 0.60, respectively, for HV images, suggesting a better performance of the co-polarized image. High temporal resolution and wide-swath SAR sequence imagery provide suitable data sources for monitoring glacier surface melting-refreezing stats; further analysis is requested to quantitatively link the volume of melting with backscatter coefficient and other SAR data sources. [ABSTRACT FROM AUTHOR]

Warm water from the Southern Ocean has a dominant impact on the evolution of Antarctic glaciers and in turn on their contribution to sea level rise. Using a continuous time series of daily-repeat satellite synthetic-aperture radar interferometry data from the ICEYE constellation collected in March-June 2023, we document an ice grounding zone, or region of tidally controlled migration of the transition boundary between grounded ice and ice afloat in the ocean, at the main trunk of Thwaites Glacier, West Antarctica, a strong contributor to sea level rise with an ice volume equivalent to a 0.6-m global sea level rise. The ice grounding zone is 6 km wide in the central part of Thwaites with shallow bed slopes, and 2 km wide along its flanks with steep basal slopes. We additionally detect irregular seawater intrusions, 5 to 10 cm in thickness, extending another 6 km upstream, at high tide, in a bed depression located beyond a bedrock ridge that impedes the glacier retreat. Seawater intrusions align well with regions predicted by the GlaDS subglacial water model to host a high-pressure distributed subglacial hydrology system in between lower-pressure subglacial channels. Pressurized seawater intrusions will induce vigorous melt of grounded ice over kilometers, making the glacier more vulnerable to ocean warming, and increasing the projections of ice mass loss. Kilometer-wide, widespread seawater intrusion beneath grounded ice may be the missing link between the rapid, past, and present changes in ice sheet mass and the slower changes replicated by ice sheet models. [ABSTRACT FROM AUTHOR]

The Greenland Ice Sheet is losing mass at increasing rates. Substantial amounts of this mass loss occur by ice discharge which is influenced by ocean thermal forcing. The ice sheet is surrounded by thousands of peripheral, dynamically decoupled glaciers. The mass loss from these glaciers is disproportionately high considering their negligible share in Greenland' overall ice mass. We study the relevance of ocean thermal forcing for ice discharge evolution in the context of this contrasting behaviour. Our estimate of ice discharge from the peripheral tidewater glaciers yields a rather stable Greenland-wide mean of 5.40 ± 3.54 Gt a−1 over 2000–2021. The evolutions of ice discharge and ocean thermal forcing are heterogeneous around Greenland. We observe a significant sector-wide increase of ice discharge in the East and a significant sector-wide decrease in the Northeast. Ocean thermal forcing shows significant increases along the northern/eastern coast, while otherwise unchanged conditions or decreases prevail. For East Greenland, this implies a clear influence of ocean thermal forcing on ice discharge. Similarly, we find clear influences at peripheral tidewater glaciers with thick termini that are similar to ice sheet outlet glaciers. At the peripheral glaciers in Northeast Greenland ice discharge evolution opposes ocean thermal forcing for unknown reasons. [ABSTRACT FROM AUTHOR]

Herein, fast fracture initiation in glacier ice is modeled using a Material Point Method and a simplified constitutive law describing tensile strain softening. Relying on a simple configuration where ice flows over a vertical step, crevasse patterns emerge and are consistent with previous observations reported in the literature. The model's few parameters allows identification of a single dimensionless number controlling fracture spacing and depth. This scaling law delineates two regimes. In the first one, ice thickness does not play a role and only ice tensile strength controls the spacing, giving rise to numerous surface crevasses, as observed in crevasse fields. In this regime, scaling can recover classical values for ice tensile strength from macroscopic field observations. The second regime, governed by ice bending, produces large‐scale, deep fractures resembling serac falls or calving events. Plain Language Summary: In ice sheets and alpine glaciers, fast‐flowing sections are often characterized by crevasse fields that play a significant role in the cryo‐hydrologic system by facilitating meltwater flow, enhancing basal sliding, weakening the ice, and impacting glacier thermodynamics. Modeling these fractures at the glacier scale remains challenging and often necessitates integrating diverse models which hinders the straightforward consideration of physical issues associated with crevasse fields on a large scale. Here, a new numerical framework allows us to conduct field‐scale experiments and paves the way for a scaling law to elucidate the macroscopic factors influencing fracture fields and to easily incorporate crevasse depth and spacing into large‐scale models. A newly discovered scaling law highlights the transition between a mechanical behavior where the regular crevasse spacing is unaffected by geometry to a regime where geometry plays a significant role, particularly in large‐scale fracture processes like glacier calving. While the numerical experiments in this paper focus on glaciers, the model and conceptual framework is versatile and can address the mechanical behavior of fractures in broader geophysical contexts such as snow, rock or ice avalanches, tectonics and landslides. Key Points: Fractures in glacier flow are modeled using material point method with elastoplasticity and tensile strain softeningA dimensional analysis reveals a key dimensionless number characterizing two different regimes of fast fractureOne regime predicts acknowledged ice tensile strength from field observations and characterizes the regular crevasse spacing [ABSTRACT FROM AUTHOR]

Cryoconite is a granular structure present on the glaciers and ice sheets found in polar regions including the Himalayas. It is composed of organic and inorganic matter which absorb solar radiations and reduce ice surface albedo, therefore impacting the melting and retreat of glaciers. Though climate warming has a serious impact on Himalayan glaciers, the biodiversity of sub-glacier ecosystems is poorly understood. Moreover, cryoconite holes are unique habitats for psychrophile biodiversity hotspots in the NW Himalayas, but unfortunately, studies on the microbial diversity of such habitats remain elusive. Therefore, the current study was designed to explore the bacterial diversity of the Hamtah Glacier Himalaya using both culturable and non-culturable approaches. The culturable bacterial count ranged from 2.0 × 10³ to 8.8 × 105 colony-forming units (CFUs)/g at the different locations of the glacier. A total of 88 bacterial isolates were isolated using the culturable approach. Based on the 16S ribosomal RNA gene (16S rRNA), the identified species belong to seven genera, namely, Cryobacterium, Duganella, Janthinobacterium, Pseudomonas, Peribacillus, Psychrobacter, and Sphingomonas. In the non-culturable approach, high-throughput sequencing of 16S rRNA genes (using MiSeq) showed unique bacterial community profiles and represented 440 genera belonging to 20 phyla, namely, Proteobacteria, Actinobacteria, Firmicutes, Bacteroidetes, Chloroflexi, Acidobacteria, Planctomycetes, Cyanobacteria, Verrucomicrobia, Spirochaetes, Elusimicrobia, Armatimonadetes, Gemmatimonadetes, Deinococcus-Thermus, Nitrospirae, Chlamydiae, Chlorobi, Deferribacteres, Fusobacteria, Lentisphaerae, and others. High relative abundances of Proteobacteria, Actinobacteria, Firmicutes, and Bacteroidetes were observed in the samples. Phototrophic (Cyanobacteria and Chloroflexi) and nitrifier (Nitrospirae) in bacterial populations indicated sustenance of the micro-ecosystem in the oligotrophic glacier environment. The isolates varied in their phenotypic characteristics, enzyme activities, and antibiotic sensitivity. Furthermore, the fatty acid profiles of bacterial isolates indicate the predominance of branched fatty acids. Iso-, anteiso-, unsaturated and saturated fatty acids together constituted a major proportion of the total fatty acid composition. High cold-adapted enzyme activities such as lipase and cellulase expressed by Cryobacterium arcticum (KY783365) and protease and cellulase activities by Pseudomonas sp. strains (KY783373, KY783377-79, KY783382) provide evidence of the possible applications of these organisms. Additionally, antibiotic tests indicated that most isolates were sensitive to antibiotics. In conclusion, the present study contributed for the first time to bacterial diversity and biopotentials of cryoconites of Hamtah Glacier, Himalayas. Furthermore, the cold-adapted enzymes and polyunsaturated fatty acids (PUFAs)may provide an opportunity for biotechnology in the Himalayas. Inductively coupled plasma mass spectrometry (ICPMS) analyses showed the presence of several elements in cryoconites, providing a clue for the acceleratingmelting and retreating of theHamtah glacier. [ABSTRACT FROM AUTHOR]

The largest uncertainty in future sea-level rise is loss of ice from the Greenland and Antarctic Ice Sheets. Ice shelves, freely floating platforms of ice that fringe the ice sheets, play a crucial role in restraining discharge of grounded ice into the ocean through buttressing. However, since the 1990s, several ice shelves have thinned, retreated, and collapsed. If this pattern continues, it could expose thick cliffs that become structurally unstable and collapse in a process called marine ice cliff instability (MICI). However, the feedbacks between calving, retreat, and other forcings are not well understood. Here we review observed modes of calving from ice shelves and marine-terminating glaciers, and their relation to environmental forces. We show that the primary driver of calving is long-term internal glaciological stress, but as ice shelves thin they may become more vulnerable to environmental forcing. This vulnerability—and the potential for MICI—comes from a combination of the distribution of preexisting flaws within the ice and regions where the stress is large enough to initiate fracture. Although significant progress has been made modeling these processes, theories must now be tested against a wide range of environmental and glaciological conditions in both modern and paleo conditions. Ice shelves, floating platforms of ice fed by ice sheets, shed mass in a near-instantaneous fashion through iceberg calving. Most ice shelves exhibit a stable cycle of calving front advance and retreat that is insensitive to small changes in environmental conditions. Some ice shelves have retreated or collapsed completely, and in the future this could expose thick cliffs that could become structurally unstable called ice cliff instability. The potential for ice shelf and ice cliff instability is controlled by the presence and evolution of flaws or fractures within the ice. [ABSTRACT FROM AUTHOR]

Most climatic and glaciological processes exhibit internal variability, which is omitted from many ice sheet model simulations. Prior studies have found that climatic variability can change ice sheet sensitivity to the long-term mean and trend in climate forcing. In this study, we use an ensemble of simulations with a stochastic large-scale ice sheet model to demonstrate that variability in frontal ablation of marine-terminating glaciers changes the mean state of the Greenland Ice Sheet through noise-induced drift. Conversely, stochastic variability in surface mass balance does not appear to cause noise-induced drift in these ensembles. We describe three potential causes for noise-induced drift identified in prior statistical physics literature: noise-induced bifurcations, multiplicative noise, and nonlinearities in noisy processes. Idealized simulations and Reynolds decomposition theory show that for marine ice sheets in particular, noise-induced bifurcations and nonlinearities in variable ice sheet processes are likely the cause of the noise-induced drift. We argue that the omnipresence of variability in climate and ice sheet systems means that the state of real-world ice sheets includes this tendency to drift. Thus, the lack of representation of such noise-induced drift in spin-up and transient ice sheet simulations is a potentially ubiquitous source of bias in ice sheet models. [ABSTRACT FROM AUTHOR]

In recent decades, the Greenland ice sheet has been losing mass through glacier retreat and ice flow acceleration. This mass loss is linked with variations in submarine melt, yet existing ocean models are either coarse global simulations focused on decadal‐scale variability or fine‐scale simulations for process‐based investigations. Here, we unite these scales with a framework to downscale from a global state estimate (15 km) into a regional model (3 km) that resolves circulation on the continental shelf. We further downscale into a fjord‐scale model (500 m) that resolves circulation inside fjords and quantifies melt. We demonstrate this approach in Scoresby Sund, East Greenland, and find that interannual variations in submarine melt at Daugaard‐Jensen glacier induced by ocean temperature variability are consistent with the decadal changes in glacier ice dynamics. This study provides a framework by which coarse‐resolution models can be refined to quantify glacier submarine melt for future ice sheet projections. Plain Language Summary: Over the past several decades, the Greenland ice sheet has been losing ice and contributing to sea‐level rise. About half of this ice loss is induced by melt that occurs where glaciers meet the ocean. Using coarse‐scale ocean models that simulate circulation around the globe, previous studies have noted a strong link between ocean temperature and enhanced glacier ice loss. However, due to the small scale of Greenland's fjords, coarse models are unable to directly quantify circulation in these fjords and melt on submerged glaciers. In this study, we develop a new framework to "zoom in" on a fjord, using high‐resolution models driven by larger coarse‐resolution models. In this approach, we simulate melt on one of Greenland's biggest glaciers and find that periods of higher melt coincide with more ice loss as observed from satellites. Since this framework is adaptable to other regions, it could also be used to simulate melt on other glaciers and support estimates of future sea‐level rise. Key Points: Subsurface temperature variability is simulated in a narrow fjord network using regional models downscaled from a global state estimateModeled increases in ocean melt at Daugaard‐Jensen glacier coincide with the onset of acceleration in 2005 and retreat and thinning in 2011Model variations in shelf‐to‐fjord ocean properties match with observations, providing a basis to estimate ocean forcing in ice projections [ABSTRACT FROM AUTHOR]

Sunlight penetrates the ice surfaces of glaciers and ice sheets, forming a water‐bearing porous ice matrix known as the weathering crust. This crust is home to a significant microbial community. Despite the potential implications of microbial processes in the weathering crust for glacial melting, biogeochemical cycles, and downstream ecosystems, there have been few explorations of its microbial communities. In our study, we used 16S rRNA gene sequencing and shotgun metagenomics of a Svalbard glacier surface catchment to characterise the microbial communities within the weathering crust, their origins and destinies, and the functional potential of the weathering crust metagenome. Our findings reveal that the bacterial community in the weathering crust is distinct from those in upstream and downstream habitats. However, it comprises two separate micro‐habitats, each with different taxa and functional categories. The interstitial porewater is dominated by Polaromonas, influenced by the transfer of snowmelt, and exported via meltwater channels. In contrast, the ice matrix is dominated by Hymenobacter, and its metagenome exhibits a diverse range of functional adaptations. Given that the global weathering crust area and the subsequent release of microbes from it are strongly responsive to climate projections for the rest of the century, our results underscore the pressing need to integrate the microbiome of the weathering crust with other communities and processes in glacial ecosystems. [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]

Understanding drivers of marine-terminating ice sheet behavior is important for constraining ice contributions to global sea level rise. In part, the stability of marine-terminating ice is influenced by solid Earth conditions at the grounded-ice margin. While the Cordilleran Ice Sheet (CIS) contributed significantly to global mean sea level during its final post-Last-Glacial-Maximum (LGM) collapse, the drivers and patterns of retreat are not well constrained. Coastal outcrops in the deglaciated Puget Lowland of Washington State – largely below sea level during glacial maxima, then uplifted above sea level via glacial isostatic adjustment (GIA) – record the late Pleistocene history of the CIS. The preservation of LGM glacial and post-LGM deglacial sediments provides a unique opportunity to assess the variability in marine ice sheet behavior of the southernmost CIS. Based on paired stratigraphic and geochronological work, with a newly developed marine reservoir correction for this region, we identify that the late-stage CIS experienced stepwise retreat into a marine environment between 15 000 and 14 000 years before present, consistent with timing of marine incursion into the region reported in earlier works. Standstill of marine-terminating ice for at least 500 years, paired with rapid vertical landscape evolution, was followed by continued retreat of ice in a subaerial environment. These results suggest rapid rates of solid Earth uplift and topographic support (e.g., grounding zone wedges) stabilized the ice margin, supporting final subaerial ice retreat. This work leads to a better understanding of shallow-marine and coastal-ice-sheet retreat and is relevant to sectors of the contemporary Antarctic and Greenland ice sheets and marine-terminating outlet glaciers. [ABSTRACT FROM AUTHOR]

The heat transfer between the warm oceanic water and the floating portion of the Antarctic ice sheet (the ice shelves) occurs in a dynamic environment with year‐to‐year changes in the distribution of icebergs and fast‐ice (the “icescape”). Dramatic events such as the collapse of glacier tongues are apparent in satellite images but oceanographic observations are insufficient to capture the synoptic impact of such events on the supply of oceanic heat to ice shelves. This study uses a 3D numerical model and semi‐ idealized experiments to examine whether the current high melting rates of ice shelves in the Amundsen Sea could be mitigated by certain icescape configurations. Specifically, the experiments quantify the impacts on oceanic heat supply of presence/absence of the Thwaites Glacier Tongue, Bear Ridge Iceberg Chain, tabular iceberg B22, and fast‐ice cover seaward of Pine Island Ice Shelf (PIS). The experiments reveal that future changes in the coastal icescape are unlikely to reverse the high ice shelf melting rates of the Amundsen Sea, and that icescape changes between 2011 and 2022 actually enhanced them slightly. Ice shelves such as Crosson and Thwaites are found to have multiple viable sources of oceanic heat whose relative importance may shift following icescape reconfigurations but the overall heat supply remains high. Similarly, the formation of a fast‐ice cover seaward of PIS slows down the cavity circulation (by 7%) but does not reduce its heat supply. [ABSTRACT FROM AUTHOR]

We present a novel and computationally efficient subglacial hydrological model that represents in a unified way both hard and soft bed rheologies as well as a dynamic switch between efficient and inefficient subglacial discharge. The subglacial model is dynamically linked to a regularized Coulomb friction law, allowing for a coupled evolution of the ice sheet on decadal to centennial time scales. The hydrological model is tested on an idealized marine ice sheet and subsequently applied to the drainage basin of Thwaites Glacier, West Antarctica, that is composed of a heterogeneous (hard/soft) bed. We find that subglacial hydrology embedded in the sliding law accelerates the grounding line retreat of Thwaites Glacier under present-day climatic conditions. Highest retreat rates are obtained for hard bed configurations and/or inefficient drainage systems.We show that the sensitivity is particularly driven by large gradients in effective pressure, more so than the value of effective pressure itself in the vicinity of the grounding line. Clearly, a better understanding of the subglacial system is needed with respect to both the spatial and temporal variability in effective pressure and bed rheological conditions. [ABSTRACT FROM AUTHOR]

The regular and robust mapping of grounding lines is essential for various applications related to the mass balance of marine ice sheets and glaciers, especially in Antarctica and Greenland. With Differential Interferometric Synthetic Aperture Radar (DInSAR) interferograms, it is possible to accurately capture the tide-induced bending of the ice shelf at a continent-wide scale and a temporal resolution of a few days. While current processing chains typically automatically generate differential interferograms, grounding lines are still primarily identified and delineated on the interferograms by a human operator. This method is time-consuming and inefficient, considering the volume of data from current and future SAR missions. We developed a pipeline that utilizes the Holistically-Nested Edge Detection (HED) neural network to delineate DInSAR interferograms automatically. We trained HED in a supervised manner using 421 manually annotated grounding lines for outlet glaciers and ice shelves on the Antarctic Ice Sheet. We also assessed the contribution of non-interferometric features like elevation, ice velocity and differential tide levels towards the delineation task. Our best-performing network generated grounding lines with a median distance of 186 m from the manual delineations. Additionally, we applied the network to generate grounding lines for undelineated interferograms, demonstrating the network's generalization capabilities and potential to generate high-resolution temporal and spatial mappings. [ABSTRACT FROM AUTHOR]

Glacier algal blooms dominate the surfaces of glaciers and ice sheets during summer melt seasons, with larger blooms anticipated in years that experience the greatest melt. Here, we characterize the glacier algal bloom proliferating on Morteratsch glacier, Switzerland, during the record 2022 melt season, when the Swiss Alps lost three times more ice than the decadal average. Glacier algal cellular abundance (cells ml-1), biovolume (µm3 cell-1), photophysiology (Fv/Fm, rETRmax), and stoichiometry (C:N ratios) were constrained across three elevations on Morteratsch glacier during late August 2022 and compared with measurements of aqueous geochemistry and outputs of nutrient spiking experiments. While a substantial glacier algal bloom was apparent during summer 2022, abundances ranged from 1.78 × 104 to 8.95 × 105 cells ml-1 of meltwater and did not scale linearly with themagnitude of the 2022melt season. Instead, spatiotemporal heterogeneity in algal distribution across Morteratsch glacier leads us to proposemelt-water-redistribution of (larger) glacier algal cells down-glacier and presumptive export of cells from the system as an important mechanism to set overall bloom carrying capacity on steep valley glaciers during high melt years. Despite the paradox of abundant glacier algae within seemingly oligotrophic surface ice, we found no evidence for inorganic nutrient limitation as an important bottom-up control within our study site, supporting our hypothesis above. Fundamental physical constraints may thus cap bloom carrying-capacities on valley glaciers as 21st century melting continues. [ABSTRACT FROM AUTHOR]

High subglacial heat flow and volcanic activity in West Antarctica contribute to instability and accelerated flow into the ocean of the West Antarctic ice sheet. In this case, a catastrophic rise in sea level by tens of centimeters – the first meters can occur in a very short geological time (years-decades) due to the rapid sliding of large masses of ice in West Antarctica into the ocean. If the Pine Island (50 cm sea level rise) or Thwaites (65 cm sea level rise) glaciers slide into the ocean, the West Antarctic Ice Sheet will lose support from these glaciers and may begin to collapse. In this case, the sea level will rise by a few meters. Based on Glen's rheological law for a two-dimensional model of the movement of ice as a nonlinear viscous fluid, the flow velocities of a 3000 m thick glacier were calculated under conditions of adhesion to the bed (~20 m/year) and under conditions of sliding along the bedrock when the lower edge of the glacier melts due to increased heat flow from below (~3000 m/year). These velocities are in good agreement with the velocities of the Pine Island, Thwaites, Amery, Denman and Totten glaciers. The rapid movement of some outlet glaciers in East Antarctica is also likely caused by melting of their bases, suggesting increased subglacial heat flow in these areas of East Antarctica. [ABSTRACT FROM AUTHOR]

Greenlandic fjords, located between the ice sheet and the ocean, are dynamic systems that can sustain highly variable levels of primary productivity and are sensitive to climate change. In our current climate trajectory, meltwater discharge is expected to significantly increase but its long‐term effects on fjord productivity are still poorly constrained. Paleo‐archives can offer valuable insights into long‐term effects. Here, we present two marine sediment core records from Nuup Kangerlua, Southwest Greenland. Our goal is to better understand to what extent, and on what time‐scales, climate fluctuations and associated glacier dynamic changes have impacted fjord productivity over the past ca. 3300 years. Our multiproxy records include diatom fluxes and assemblage composition, sediment biogeochemistry, and grain‐size distribution. Our study reveals that fjord productivity is tightly linked to regional climate variability; relatively higher productivity levels coincided with mild climate periods whereas the climate cooling of the last millennium led to a decrease in productivity. The diatom records suggest that lower productivity is associated with shorter or less intense summer blooms, increased sea‐ice cover and/or a stratified water column. Diatom assemblages demonstrate cold sea‐surface conditions around 1600 CE that might be linked to local advance of glaciers. Cold conditions and decreasing productivity culminated at 1850 CE, when glaciers in the fjord retreated and high glacial meltwater discharge would have altered the fjord hydrography, likely leading to limited nutrient availability. Our long‐term records support the idea that changing climate and cryosphere conditions have a non‐linear impact on the productivity of Greenlandic fjords. Key Points: Diatom assemblages record multiannual to multidecadal changes in the timing, magnitude, and composition of phytoplankton bloomsFjord productivity is linked to climate variability and to associated changes in marine‐terminating glaciersHigher fjord productivity coincided with warmer climate periods while the lowest productivity was recorded at the end of the Little Ice Age [ABSTRACT FROM AUTHOR]

The Antarctic Ice Sheet is the largest potential contributor to global sea level rise due to global climate warming, making it important to monitor changes in the marginal outlet glacier dynamics and surface topography. The unmanned aerial vehicle (UAV), a new platform, has been proven to be a useful tool in glaciological applications because it can collect data in a variety of ways and has a high spatial and temporal resolution compared to traditional ground-based measurements and remote sensing from space. In this study, we used the combined UAV platform composed of the DJI Mavic 3 Enterprise (M3E) and M300 (L1) versions to perform high-accuracy on-site investigations of Dalk Glacier, a typical outlet glacier near the Zhongshan Station in East Antarctica. During the 39th Chinese Antarctic Expedition (CHINARE) in 2022–2023, UAV surveys with multiple sensors were performed. The M3E optical camera makes it possible to reconstruct high-resolution ortho-mosaic and digital elevation models (DEMs) of Dalk Glacier with high accuracy based on photogrammetric principles, helping us to extract the features of glacier surface topography, while dense point clouds derived from the L1 LiDAR camera are used to validate and improve the three-dimensional accuracy. To further monitor the dynamics and stability of Dalk Glacier, we generate ice flow velocity maps using ortho-mosaics for the years 2019, 2020, and 2023. The spatiotemporal changes in ice flow velocities were further compared with the products derived from satellite remote sensing data. [ABSTRACT FROM AUTHOR]

The new photon-counting laser altimetry satellite ICESat-2 was successfully launched on September 15, 2018 with an unprecedented ice surface elevation measurement accuracy of 2–4 cm. It is meaningful for accurately estimating volumetric changes in the Antarctic Ice Sheet. Cross-validation of different types of data, especially comparison with in-situ data, is important for the ice sheet mass balance results of a new satellite. This paper proposes an elevation-difference method of grid elevation change rate model, based on ICESat-2 ATL11 elevation time series data, to distinguish the linear change trend of ice sheet surface elevation from the elevation change resulting from periodic precipitation. In order to compare the estimated elevation-change rate with the flux computed from in-situ snow stake velocity measurements and GRACE-FO gravity survey data, we made corrections for firn air content, elastic, and glacier isostatic adjustment. Based on the ATL11 data from 2019 to 2023, our results show that the ice sheet change in Basin 11 along the CHINARE traverse is from ~0.019 m yr−1 to ~0.121 m yr−1, and the mass balance in the upstream of the traverse in Basin 11 is ~1.9±0.2 Gt yr−1. It is comparable to ~2.4±1.2 Gt yr−1 from GRACE-FO during the same time period. Furthermore, the flux across traverse-11.9±1.1 Gt yr−1 is comparable to that of ~-9.7±0.9 Gt yr−1 across the same flux gate during 1997–2009 which is calculated based on GNSS-derived ice velocity observations, considering the time period difference and uncertainties. [ABSTRACT FROM AUTHOR]