This work presents a reflection on the meaning and significance of knowledge coproduction in the field of glaciology. We start by invoking the paradigm of Structure–Form–Environment Interplay (SFEI) to formulate a generalised definition of glaciology, which highlights the relevance of knowledge coproduction. The adoption of a relational view of glaciological knowledge leads us to identify five core dimensions of knowledge coproduction: purpose, ethics, ambiguity, inclusion/exclusion, and relationships. Based on those dimensions, we delve into the decisive methodological aspects of the coproduction process, namely the definition of its purpose, the identification of participants, the organisation of the process, the recognition of ambiguity in Ways of Knowing (WoKs), and the consideration of ethical implications. In addition to the already known three stages of knowledge coproduction process (codesign, codevelopment, and codelivery), we propose the inclusion of an additional preparation stage, which entails the acknowledgment of the identity and involvement of all human and nonhuman participants, their positionality, and means to ensure their cultural and ontological safety. We reason that knowledge coproduction does not replace the scientific method, but rather complements it, eliciting the possibility to unveil deeper insights that might be difficult to attain through unilateral means.

We present the first simulations of Jakobshavn Isbrae (Sermeq Kujalleq), west Greenland, using a 3D Stokes calving model that permits unrestricted advance and retreat. Using the position-based crevasse-depth calving law, the model is applied to simulate the calving dynamics of 2016–2017 season when Jakobshavn Isbrae is assumed to be stable because of the presence of a strong proglacial ice mélange. The calving law needs to be adjusted to avoid an underestimation of calving, but once adjusted the calving model simulates seasonal calving dynamics that reflect observed calving-driven retreat very well. We find that a crevasse penetration threshold of 94.5% best matches observations from satellite imagery. Additional, 2-year transient simulations show that although ice mélange is essential to the glacier's winter readvance, when removed, the glacier only retreats a couple of kilometres before reaching a stable position. While the backstress provided by the ice mélange allows the glacier to advance beyond this point, the retreated terminus position is determined by a combination of bed geometry and glacier dynamics. Ultimately, while the ice mélange allows winter readvance, cessation of the well-documented rapid retreat of Jakobshavn Isbrae will be influenced by the bed geometry.

An array of information about the Antarctic ice sheet can be extracted from ice-sheet internal architecture imaged by airborne ice-penetrating radar surveys. We identify, trace and date three key internal reflection horizons (IRHs) across multiple radar surveys from South Pole to Dome A, East Antarctica. Ages of ~38 ± 2.2, ~90 ± 3.6 and ~162 ± 6.7 ka are assigned to the three IRHs, with verification of the upper IRH age from the South Pole ice core. The resultant englacial stratigraphy is used to identify the locations of the oldest ice, specifically in the upper Byrd Glacier catchment and the Gamburtsev Subglacial Mountains. The distinct glaciological conditions of the Gamburtsev Mountains, including slower ice flow, low geothermal heat flux and frozen base, make it the more likely to host the oldest ice. We also observe a distinct drawdown of IRH geometry around South Pole, indicative of melting from enhanced geothermal heat flux or the removal of deeper, older ice under a previous faster ice flow regime. Our traced IRHs underpin the wider objective to develop a continental-scale database of IRHs which will constrain and validate future ice-sheet modelling and the history of the Antarctic ice sheet.

Knowledge of ice-sheet catchments is critical for mass-balance assessments, especially glacier-scale input–output budgets. This study explores variations in the catchment of Sermeq Kujalleq, or Jakobshavn Isbrø, Greenland. Six observation-based catchment delineations are evaluated along with a 16-member catchment ensemble calculated from ice-sheet models within the Ice Sheet Model Intercomparison Project for CMIP6 (ISMIP6). The ‘present-day’ ISMIP6 ensemble mean area was found to be $\sim 6.3\%$ larger than the mean of the observed catchments. Ensemble spreads were comparable in size, $\pm 12.3\%$ and $\pm 15.4\%$, suggesting models are able to delineate the present-day catchment with the same degree of uncertainty as observational methods. The mean catchment area of a 13-member ISMIP6 ensemble shows temporal variation, increasing $\sim [ 2.7,\; \, 5.7,\; \, 9.1] \%$ under three ocean forcing scenarios and a RCP8.5 projection based on one GCM from 2015 to 2100, primarily as the southern catchment boundary migrates southward. This is interpreted as Sermeq Kujalleq exhibiting dynamic piracy, re-directing ice away from adjacent land terminating glaciers. For mass-balance assessments, present-day catchment delineation is more important than capturing the temporal evolution of individual catchments. However, the modeled temporal changes in catchment area are potentially underestimated, as the models exhibit insufficient acceleration of inland ice flow.

Recent advances in remote sensing have produced global glacier surface elevation change data. Parsing these elevation change signals into contributions from the climate (i.e. climatic mass balance) and glacier dynamics (i.e. flux divergence) is critical to enhance our process-based understanding of glacier change. In this study, we evaluate three approaches for direct, continuous measurements of the climatic mass balance, flux divergence and elevation change at a site on Gulkana Glacier in Alaska using low-cost global navigation satellite system (GNSS) sensors, GNSS interferometric reflectometry (GNSS-IR), banded ablation stakes with time-lapse cameras and combinations thereof. Cumulative climatic mass balance over the season was 4.85 m and the three approaches were within 0.08 m through early July before the snowpack melted, and within 0.28 m through mid-August. The flux divergence increased from 0.52 ± 0.03 cm d−1 before June 3 to about 0.73 cm d−1 after June 27. We demonstrate a single GNSS system fixed atop an ablation stake can measure contemporaneous climatic mass balance, flux divergence and elevation change based on the antenna's position and GNSS-IR techniques. The ability of these systems to measure glacier mass balance and flux divergence offers unique opportunities for year-round observations on mountain glaciers in the future.

Müller Ice Cap sits on Umingmat Nunaat (Axel Heiberg Island), Nunavut, Canada, ~ 80°N. Its high latitude and elevation suggest it experiences relatively little melt and preserves an undisturbed paleoclimate record. Here, we present a suite of field measurements, complemented by remote sensing, that constrain the ice thickness, accumulation rate, temperature, ice-flow velocity, and surface-elevation change of Müller Ice Cap. These measurements show that some areas near the top of the ice cap are more than 600 m thick, have nearly stable surface elevation, and flow slowly, making them good candidates for an ice core. The current mean annual surface temperature is −19.6 °C, which combined with modeling of the temperature profile indicates that the ice is frozen to the bed. Modeling of the depth-age scale indicates that Pleistocene ice is likely to exist with measurable resolution (300–1000 yr m−1) 20–90 m from the bed, assuming that Müller Ice Cap survived the Holocene Climatic Optimum with substantial ice thickness (~400 m or more). These conditions suggest that an undisturbed Holocene climate record could likely be recovered from Müller Ice Cap. We suggest 91.795°W, 79.874°N as the most promising drill site.

The sensitivity of glacier mass balance to temperature and precipitation variations is crucial for informing models that simulate glaciers’ response to climate change. In this study, we simulate the glacier-wide mass balance of Mera Glacier with a surface energy-balance model, driven by in situ meteorological data, from 2016 to 2020. The analysis of the share of the energy fluxes of the glacier shows the radiative fluxes account for almost all the energy available during the melt season (May–October). However, turbulent fluxes are significant outside the monsoon (June–September). On an annual scale, melt is the dominant mass flux at all elevations, but 44% of the melt refreezes across the glacier. By reshuffling the available observations, we create 180 synthetic series of hourly meteorological forcings to force the model over a wide range of plausible climate conditions. A +1 (−1)°C change in temperature results in a −0.75 ± 0.17 (+0.93 ± 0.18) m w.e. change in glacier-wide mass balance and a +20 (−20)% change in precipitation results in a +0.52 ± 0.10 (−0.60 ± 0.11) m w.e. change. Our study highlights the need for physical-based approaches to produce consistent forcing datasets, and calls for more meteorological and glaciological measurements in High Mountain Asia.

Basal crevasses threaten the stability of ice shelves through the potential to form rifts and calve icebergs. Furthermore, it is important to determine the dependence of crevasse stability on temperature due to large vertical temperature variations on ice shelves. In this work, considering the vertical temperature profile through ice viscosity, we compare (1) the theoretical crack depths and (2) the threshold stress causing the transition from basal crevasses to full thickness fractures in several fracture theories. In the Zero Stress approximation, the depth-integrated force at the crevassed and non-crevassed location are unbalanced, violating the volume-integrated Stokes equation. By incorporating a Horizontal Force Balance (HFB) argument, recent work showed analytically that the threshold stress for rift initiation is only half of that predicted by the Zero Stress approximation. We generalize the HFB theory to show that while the temperature profile influences crack depths, the threshold rifting stress is insensitive to temperature. We compare with observations and find that HFB best matches observed rifts. Using HFB instead of Zero Stress for cracks in an ice-sheet model would substantially enlarge the predicted fracture depth, reduce the threshold rifting stress and potentially increase the projected rate of ice shelf mass loss.

Tidewater glaciers frequently advance and retreat in ways uncoupled from climate forcing. This complicates the task of forecasting the evolution of individual glaciers and the overall Greenland ice sheet, much of which is drained by tidewater glaciers. Past observational research has identified a set of processes collectively known as the tidewater glacier cycle (TGC) to describe tidewater glacier evolution in four stages: the advancing stage, the extended stage, the retreating stage and the retreated stage. Once glacier retreat is initiated, the TGC is thought to depend largely on the glacier's calving rate, which is controlled by fjord geometry. However, there has been little modeling or systematic observational work on the topic. Measuring calving rates directly is challenging and thus we developed an averaged von Mises stress state at the glacier terminus as a calving rate proxy that can be estimated from surface velocities, ice thickness, a terminus position and subglacial topography. We then analyzed 44 tidewater glaciers in Greenland and assessed the current state in the TGC for them. Of the 44 glaciers, we find that fjord geometry is causing instability in ten cases, vs stability in seven, with 11 in rapid retreat and 16 have been historically stable.

Knowledge of snow and firn-density change is needed to use elevation-change measurements to estimate glacier mass change. Additionally, firn-density evolution on glaciers is closely connected to meltwater percolation, refreezing and runoff, which are key processes for glacier mass balance and hydrology. Since 2016, the U.S. Geological Survey Benchmark Glacier Project has recovered firn cores from a site on Wolverine Glacier in Alaska's Kenai Mountains. We use annual horizons in repeat cores to track firn densification and meltwater retention over seasonal and interannual timescales, and we use density measurements to quantify how the firn air content (FAC) changes through time. The results suggest the firn is densifying due primarily to compaction rather than refreezing. Liquid-water retention in the firn is transient, likely due to gravity-fed drainage and irreducible-water-content decreases that accompany decreasing porosity. We show that the uncertainty (±60 kg m−3) in the commonly used volume-to-mass conversion factor of 850 kg m−3 is an underestimation when glacier-wide FAC variability exceeds 12% of the glacier-averaged height change. Our results demonstrate how direct measurements of firn properties on mountain glaciers can be used to better quantify the uncertainty in geodetic volume-to-mass conversions.

At high elevations on the Greenland ice sheet meltwater percolates and refreezes in place, and hence does not contribute to mass loss. However, meltwater generation and associated surface runoff is occurring from increasingly higher altitudes, causing changes in firn stratigraphy that have led to the presence of near-surface ice slabs. These ice slabs force meltwater to flow laterally instead of percolating downwards. Here we present a simple, physics-based quasi-2-D model to simulate lateral meltwater runoff and superimposed ice (SI) formation on top of ice slabs. Using an Eulerian Darcy flow scheme, the model calculates how far meltwater can travel within a melt season and when it appears at the snow surface. Results show that lateral flow is a highly efficient runoff mechanism, as lateral outflow exceeds locally generated meltwater in all model gridcells, with total meltwater discharge sometimes reaching more than 30 times the average amount of in situ generated melt. SI formation, an important process in the formation and thickening of the ice slabs, can retain up to 40% of the available meltwater, and generally delays the appearance of visible runoff. Validating the model against field- or remote-sensing data remains challenging, but the results presented here are a first step towards a more comprehensive understanding and description of the hydrological system in the accumulation zone of the southwestern Greenland ice sheet.

Ice shelves regulate the flow of the Antarctic ice sheet toward the ocean and its contribution to sea-level rise. Accurately monitoring the basal and surface melting of ice shelves is therefore essential for predicting the ice sheet's response to climatic warming. In this study, we utilize Sentinel-1A synthetic aperture radar satellite imagery combined with shipboard measurements of water temperature and salinity to investigate the presence of surficial meltwater plumes along the Antarctic coastline. Our approach reveals a strong correlation between areas of pronounced low radar backscatter extending from ice shelves and significant decreases in water temperature and salinity, suggesting meltwater-enriched ocean waters. We propose that the low radar backscatter signature of meltwater outflows is caused by stable stratification of the upper water column, driven by density contrasts from buoyant, low-salinity meltwater and surface current shear that reduce Bragg scattering waves. The resulting smooth water surfaces were observed adjacent to the surface expression of deep basal channels, documented in a helicopter survey along part of the Bellingshausen Sea ice edge. We present high-temporal resolution satellite radar as a tool for identifying meltwater release from beneath ice shelves, capable of all-weather, day-and-night imaging.

Outlet glaciers in Greenland are undergoing retreat and diffusive thinning in response to external forcings, but the rates and magnitudes of these responses differ from glacier to glacier for unclear reasons. We test how changes in ice overburden pressure and basal lubrication affect diffusive thinning rates and their spatial patterns by conducting numerical experiments over various idealized Greenland-like glacier domains. We find that ~10 km frontal retreat over a decade can produce sustained thinning rates as large as 16 m a−1 due to ice overburden pressure changes, at outlet glaciers with high basal drag (>60 kPa) and lateral resistive stress (>70 kPa). Localized basal lubrication perturbations induce upstream thinning and downstream thickening up to 12 m a−1; the duration of the lubrication forcing generally has a greater effect than its intensity on induced thickness changes. Lastly, episodic grounding line retreats over a rough bed produce a stepped time series of thinning broadly consistent with observations of dynamic elevation change on multiple Greenland glaciers. Our findings highlight the critical role of the total grounding zone – not ice front position – through the resistive stress change in relation to total glacier thinning.

We conducted a series of experiments that revealed the formation of mm-scale penitente structures in ice illuminated by broadband light under moderate vacuum conditions between 50 and 2000 Pa. The experimental apparatus consists of a 0.3 m diameter cylindrical vacuum chamber with a cooling jacket surrounding the outer radius and bottom surface. Light shines in through an optical window at the top to illuminate most of the ice surface. We observe penitente-like structures at temperatures between −15$^\circ$C and $-2^\circ$C and pressures close to the equilibrium vapor pressure at the ice surface temperature. The formation of these structures is very sensitive to slight changes in background pressure, and the structures tend to vanish with significant deviations away from the equilibrium curve, resulting in a smooth sublimated crater formation instead of penitentes. Application of the physical model by Claudin and others (2015, doi: 10.1103/PhysRevE.92.033015) at experimental conditions generally agrees with observations for penitente spacing.

Ruth Glacier is situated in the Central Alaska Range, with the Don Sheldon Amphitheater comprising much of its broad accumulation area, directly adjacent to North America's tallest mountain, Denali. From there it funnels through the ‘Great Gorge,’ flanked by steep valley walls reaching over 1500 m. We combine airborne and ground-based radar measurements of ice thickness with satellite-derived surface velocities to constrain ice flux above and below the gorge, and employ a mass conservation approach to estimate the glacier's thickness within the gorge. We measure ice thickness in the amphitheater to reach 950 m, and estimate centerline thickness in the gorge to range from 610 to 960 m. Our estimates are up to two times greater than those suggested by global models, and allow us to confirm that the Great Gorge rivals Hells Canyon as the deepest gorge in North America. We found that the geometry of the gorge prevents radar measurements of ice thickness there since returns from the subglacial valley walls would precede and potentially occlude nadir bed returns. The same may be true of other unmapped mountain glaciers; however, thickness may be determined using appropriately located flux gates where radar sounding is feasible, combined with mass conservation methods.

Rock glaciers support alpine biodiversity and may respond more slowly to warming than snow or glaciers. While responses of snow and glaciers to climate change are relatively well understood, a robust assessment of rock glacier environmental niche, future distributions of rock glaciers and potential for development of rock glaciers from current glaciers is lacking. Using process-relevant, high-resolution environmental descriptors, we develop a species distribution model of the topographic, geologic and hydroclimatic niche of rock glaciers that provides novel estimates of potential rock glacier distributions for different climates. We identify mean annual air temperature and headwall area as the dominant controls on rock glacier spatial distributions, with rock glaciers more likely to be found in areas with mean annual temperatures close to −5°C, little rain, northern aspects and broad headwalls. While rock glacier climate equilibration may take hundreds of years, we find that equilibration to present climate will result in a 50% reduction in rock glacier habitat and equilibration to late 21st-century climate under a high-end warming scenario will result in a 99% reduction in rock glacier habitat across the western USA. Under future conditions, we find limited potential for glacier to rock glacier transformation (3% of glacierized area), concentrated in cold, high elevation, moderate precipitation areas.

The time-evolution of glacier basal motion remains poorly constrained, despite its importance in understanding the response of glaciers to climate warming. Athabasca Glacier provides an ideal site for observing changes in basal motion over long timescales. Studies from the 1960s provide an in situ baseline dataset constraining ice deformation and basal motion. We use two complementary numerical flow models to investigate changes along a well-studied transverse profile and throughout a larger study area. A cross-sectional flow model allows us to calculate transverse englacial velocity fields to simulate modern and historical conditions. We subsequently use a 3-D numerical ice flow model, Icepack, to estimate changes in basal friction by inverting known surface velocities. Our results reproduce observed velocities well using standard values for flow parameters. They show that basal motion declined significantly (30–40%) and this constitutes the majority (50–80%) of the observed decrease in surface velocities. At the same time, basal resistive stress has remained nearly constant and now balances a much larger fraction of the driving stress. The decline in basal motion over multiple decades of climate warming could serve as a stabilizing feedback mechanism, slowing ice transport to lower elevations, and therefore moderating future mass loss rates.

Supraglacial lakes on Antarctic ice shelves can have far-reaching implications for ice-sheet stability, highlighting the need to understand their dynamics, controls and role in the ice-sheet mass budget. We combine a detailed satellite-based record of seasonal lake evolution in Dronning Maud Land with a high-resolution simulation from the regional climate model Modèle Atmosphérique Régional to identify drivers of lake variability between 2014 and 2021. Correlations between summer lake extents and climate parameters reveal complex relationships that vary both in space and time. Shortwave radiation contributes positively to the energy budget during summer melt seasons, but summers with enhanced longwave radiation are more prone to surface melting and ponding, which is further enhanced by advected heat from summer precipitation. In contrast, previous winter precipitation has a negative effect on summer lake extents, presumably by increasing albedo and pore space, delaying the accumulation of meltwater. Downslope katabatic or föhn winds promote ponding around the grounding zones of some ice shelves. At a larger scale, we find that summers during periods of negative southern annular mode are associated with increased ponding in Dronning Maud Land. The high variability in seasonal lake extents indicates that these ice shelves are highly sensitive to future warming or intensified extreme events.

Glaciers in High-Mountain Asia are experiencing varying rates and patterns of mass loss due to a complex interplay between glacier surface processes, local conditions and climate forcing. Spatially distributed surface mass balance (SMB) estimates can provide valuable insight into these drivers, but observations are currently limited in both space and time. We used very-high-resolution optical stereo images acquired by commercial satellites to prepare time series of digital elevation models (DEMs), and derived contemporaneous surface velocity and elevation change products for six debris-covered glaciers in Nepal. We developed new methods to produce flow-corrected Lagrangian SMB maps to isolate local surface ablation signals with enough detail to study individual ice cliffs. Our results show reduced ablation under thick debris cover and enhanced ablation over ice cliffs. Ablating ice cliffs were responsible for $10\!-\!38\%$ of the total ablation over debris-covered areas, even though they covered $\leq \!11\%$ of the total area. Seasonal SMB products reveal the timing and patterns of summer accumulation and ablation, underscoring the importance of snow avalanches for low-elevation debris-covered glaciers in the region. Our approach can be applied to other glaciers with repeat high-resolution DEM coverage and extended for regional analyses of SMB on seasonal to interannual timescales.