To understand the mechanism of simultaneous drainage event related to supraglacial lakes on a debris-covered glacier, we investigated water-level variations of supraglacial lakes on the southern Inylchek Glacier in Kyrgyzstan. To examine these variations, we used daily aerial images for 2017–2019 from an uncrewed aerial vehicle that were converted to 15 cm-digital surface models and ortho-images. Our main results are as follows: (1) When one lake drained, the water levels of other lakes simultaneously increased, indicating that drainage water is shared with several lakes through a main englacial conduit. In one drainage event, a branched off englacial conduit clearly connected to a main englacial conduit. (2) Sometimes several lakes discharged simultaneously, indicating that several lakes had connected to a main englacial conduit that had opened. Such cases can cause larger-scale drainage than that from the opening of a branched off englacial conduit. (3) Simultaneous drainage occurred twice in the same year, each time through a different conduit, indicating that the main englacial conduit can be abandoned and reused. (4) In some lakes, the water level on the hydraulic gradient line increased gradually with nearly the same increase rate just before drainage. Such an increase may be an indicator of a possible simultaneous drainage event.

This study investigates stagnation conditions of the Pensilungpa glacier, western Himalaya. Multiple glacier parameters (length, area, debris extent and thickness, snowline altitude (SLA), velocity, downwasting and ice cliffs) were studied using field measurements (2016–18), high-resolution imagery from GoogleEarth (2013–17) and spaceborne Landsat, ASTER and SRTM data (1993–2017) to comprehend the glacier's current state. Results show a moderate decrease in length (6.62 ± 2.11 m a−1) and area (0.11 ± 0.03% a−1), a marked increase in SLA (~6 m a−1) and debris cover (2.86 ± 0.29% a−1) and a slowdown of ~50% during 1993–2016. Notable thinning of −0.88 ± 0.04 m a−1 was observed between 2000 and 2017 showing a similar trend as field measurements during 2016–17 (−0.88 m) and 2017–18 (−1.54 m). Further, results reveal a stagnation of the lower ablation zone (LAZ). Less mass supply and heterogeneous debris growth (6.67 ± 0.41% a−1) over the previous decade resulted in slowdown, margin insulation and slope-inversion, leading to stagnation. Stagnation of LAZ caused bulging in the dynamic upper ablation zone and favored the development of supraglacial ponds and ice cliffs. Ice cliffs have grown significantly (48% in number; 41% in area during 2013–17) and their back-wasting now dominates the ablation process.

We investigate the errors caused by neglecting the crystal-orientation fabric when inferring the basal friction coefficient field, and whether such errors can be alleviated by inferring an isotropic enhancement factor field to compensate for missing fabric information. We calculate the steady states that arise from ice flowing over a sticky spot and a bedrock bump using a vertical-slab numerical ice-flow model, consisting of a Weertman sliding law and the anisotropic Johnson flow law, coupled to a spectral fabric model of lattice rotation and dynamic recrystallisation. Given the steady or transient states as input for a canonical adjoint-based inversion, we find that Glen's isotropic flow law cannot necessarily be used to infer the true basal drag or friction coefficient field, which are obscured by the orientation fabric, thus potentially affecting vertically integrated mass fluxes. By inverting for an equivalent isotropic enhancement factor, a more accurate mass flux can be recovered, suggesting that joint inversions for basal friction and the isotropic flow-rate factor may be able to compensate for mechanical anisotropies caused by the fabric. Thus, in addition to other sources of rheological uncertainty, fabric might complicate attempts to relate subglacial conditions to basal properties inferred from an inversion relying on Glen's law.

The overall trend of rapid retreat of Alpine glaciers contains considerable variability of responses at the scale of individual glaciers. As a step towards a regional assessment of glacier state that allows a detailed differentiation of single glaciers, we explore the potential of a self-organizing maps (SOM) algorithm to identify and cluster recurring patterns of thickness change at glaciers in western Austria. Using digital elevation models and glacier inventories for three time periods, we compute the frequency distribution of surface elevation change over the area of each glacier in the data set, for each period. The results of the SOM clustering show a distinct pattern shift over time: From 1969 to 1997, surface elevation change occurred at relatively uniform rates across a given glacier. Since 1997, the distribution of surface elevation change at individual glaciers has been far less uniform, indicating accelerated processes of disintegration. Tracking the evolution of individual glaciers throughout the time periods via the clusters highlights both the broader regional trend as well as glaciers that deviate from this trend, e.g. some very small, high elevation glaciers that have reverted to reduced and more uniform volume loss patterns.

Constant-load creep tests were performed at −10°C at various compressive stresses from 0.05 to 0.75 MPa on specimens taken every 10 m along a firn core extracted at Summit, Greenland in June 2017. The microstructures before and after creep testing were examined using both X-ray microtomography (micro-CT) and optical images from thin sections. An Andrade-like equation was used to describe the primary creep behavior and yielded the time exponent k of 0.17–0.76. The onset of secondary creep occurred at strains of ~0.5–3% but was sometimes not observed at all in shallow firn specimens and at stresses ⩽0.43 MPa even for strain up to 32%. For the 50–80 m firn crept at stresses ⩾0.55 MPa, secondary creep occurred at strains of 2.6 ± 0.28%, and the stress exponent, n, in Glen's law, was found to range from 4.1 to 4.6, similar to those observed for fully dense ice. Micro-CT observations of crept specimens showed that in most cases, the specific surface area, the total porosity and the structure model index decreased, while the structure thickness increased with increasing density. These microstructural characteristics are consistent with the densification of the firn. Optical images from thin sections showed that recrystallization occurred in some specimens that had undergone secondary creep.

The ablation zones of debris-covered glaciers in Himalaya exhibit heterogeneous processes and melt patterns. Although sub-debris melt is measured at ablation stakes, the high variability of debris thickness necessitates distributed melt measurements at the glacier scale. Focusing on Annapurna III Glacier, we used uncrewed aerial system (UAS) photogrammetry to estimate total volume loss and slope-perpendicular glacier melt between May and November 2019 using flow-corrected point clouds. Results indicated the average elevation change was −1.10 ± 0.19 m, while the mean melt was −0.87 m w.e., equating to a mean melt rate of −0.47 cm w.e. d−1. However, the spatial pattern was highly variable due to complex local processes necessitating future study over short intervals. The evaluation of specific areas showed the interplay of debris thickness variability, subseasonal debris redistribution, supraglacial channel reconfiguration and the imprint of relict ice cliffs in leading to contemporary melt rates. Ice cliffs had higher melt distances (mean −3.9 ± 0.19 m) compared to non-cliff areas (mean −0.75 ± 0.19 m) and were the predominant control on the spatial patterns of seasonal melt rates. Crucially, the definition of ice cliff areas from thinning data has a profound impact on derived melt rates and melt enhancement. Our study demonstrates the possibility and utility of deriving fully-distributed slope-perpendicular melt measurements.

Glaciers in the Southern Patagonia Icefield (SPI) have been shrinking in recent decades, but due to a lack of field observations, understanding of the drivers of ablation is limited. We present a distributed surface energy balance model, forced with meteorological observations from a west–east transect located in the north of the SPI. Between October 2015 and June 2016, humid and warm on-glacier conditions prevailed on the western side compared to dry and cold conditions on the eastern side. Controls of ablation differ along the transect, although at glacier-wide scale sensible heat (mean of 72 W m−2 to the west and 51 W m−2 to the east) and net shortwave radiation (mean of 54 W m−2 to the west and 52 W m−2 to the east) provided the main energy sources. Net longwave radiation was an energy sink, while latent heat was the most spatially variable flux, being an energy sink in the east (−4 W m−2) and a source in the west (20 W m−2). Ablation was high, but at comparable elevations, it was greater to the west. These results provide new insights into the spatial variability of energy-balance fluxes and their control over the ablation of Patagonian glaciers.

Liquid water can exist at temperatures well below freezing beneath glaciers and ice sheets, where subglacial water systems, fresh and saline, have been shown to host unique microbial ecosystems. Geophysical techniques sensitive to fluid-content contrasts, e.g. electromagnetics, can characterize subglacial water and its salinity. Here, we assess the ground-based transient electromagnetic (TEM) method for deriving the resistivity and salinity of subglacial water. We adapt an existing open-source Bayesian inversion algorithm, which uses independent depth constraints, to output posterior distributions of resistivity and pore fluid salinity with depth. A variety of synthetic models, including a thin (5 m), conductive (0.16 Ωm), hypersaline (147 psu) subglacial lake, are used to evaluate the TEM method for imaging under 800 m-thick ice. The study demonstrates that TEM methods can resolve conductive, saline bodies accurately using external depth constraints, for example, from radar or seismic data. The depth resolution of TEM can be limited beneath deep (>800 m), thick (>50 m) conductive, water bodies and additional constraints from passive electromagnetic (EM) methods could be used to reduce ambiguities in the TEM results. Subsequently, non-invasive active and passive EM methods could provide profound insights into remote aqueous systems under glaciers and ice sheets.

The mechanics underlying ice–skate friction remain uncertain despite over a century of study. In the 1930s, the theory of self-lubrication from frictional heat supplanted an earlier hypothesis that pressure melting governed skate friction. More recently, researchers have suggested that a layer of abraded wear particles or the presence of quasi-liquid molecular layers on the surface of ice could account for its slipperiness. Here, we assess the dominant hypotheses proposed to govern ice–skate friction and describe experiments conducted in an indoor skating rink aimed to provide observations to test these hypotheses. Our results indicate that the brittle failure of ice under rapid compression plays a strong role. Our observations did not confirm the presence of full-contact water films and are more consistent with the presence of lubricating ice-rich slurries at discontinuous high-pressure zones (HPZs). The presence of ice-rich slurries supporting skates through HPZs merges pressure-melting, abrasion and lubricating films as a unified hypothesis for why skates are so slippery across broad ranges of speeds, temperatures and normal loads. We suggest tribometer experiments to overcome the difficulties of investigating these processes during actual skating trials.

In this study, we compare equilibrium-line altitudes (ELAs) calculated using the area–altitude balance ratio (AABR) and the accumulation–area ratio (AAR) methods, with measured ELAs derived from direct field observations. We utilise a GIS toolbox to calculate the ELA for 64 extant glaciers by applying the AABR and AAR methods to DEMs and polygons of their geometry. The calculated ELAs (c-ELAs) are then compared to measured zero-net balance ELAs (znb-ELAs) obtained from mass-balance time series held by the WGMS for the same glaciers. The correlation between znb-ELAs and AABR (1.56)/AAR (0.58) c-ELAs is very strong, with an r2 = 0.99. The smallest median difference between znb-ELAs and c-ELAs (i.e. 65.5 m) is obtained when a globally representative AABR of 1.56 is used. When applied to palaeoglacier-climate applications, this difference translates to ~0.42°C, well within the uncertainty of palaeotemperature proxies used to determine mean summer temperature at the ELA. The more widely used mean AABR of 1.75 is shown to be statistically invalid due to the skewness of the dataset. On this basis, when calculating glacier ELAs, we recommend the use of a global AABR value of 1.56.

Recent progress has been made in quantifying snowmelt in the Himalaya. Although the conditions are favorable for refreezing, little is known about the spatial variability of meltwater refreezing, hindering a complete understanding of seasonal snowmelt dynamics. This study aims to improve our understanding about how refreezing varies in space and time. We simulated refreezing with the seNorge (v2.0) snow model for the Langtang catchment, Nepalese Himalaya, covering a 5-year period. Meteorological forcing data were derived from a unique elaborate network of meteorological stations and high-resolution meteorological simulations. The results show that the annual catchment average refreezing amounts to 122 mm w.e. (21% of the melt), and varies strongly in space depending on elevation and aspect. In addition, there is a seasonal altitudinal variability related to air temperature and snow depth, with most refreezing during the early melt season. Substantial intra-annual variability resulted from fluctuations in snowfall. Daily refreezing simulations decreased by 84% (annual catchment average of 19 mm w.e.) compared to hourly simulations, emphasizing the importance of using sub-daily time steps to capture melt–refreeze cycles. Climate sensitivity experiments revealed that refreezing is highly sensitive to changes in air temperature as a 2°C increase leads to a refreezing decrease of 35%.

Changes in water pressure at the beds of glaciers greatly modify their sliding rate, affecting rates of ice mass loss and sea level change. However, there is still no agreement about the physics of subglacial sliding or how water affects it. Here, we present a new simplified physical model for the effect of transient subglacial hydrology on basal ice velocity. This model assumes that a fraction of the glacier bed is connected by an active hydrologic system that, when averaged over an appropriate scale, is governed by two parameters with limited spatial variability. The sliding model is reminiscent of Budd's empirical sliding law but with fundamental differences including a dependence on the fractional area of the active hydrologic system. With periodic surface meltwater forcing, the model displays classic diffusion-wave behavior, with a downstream time lag and decay of subglacial water pressure perturbations. Testing the model against Greenland observations suggests that, despite its simplicity, it captures key features of observed proglacial discharges and ice velocities with reasonable physical parameter values. Given these encouraging findings, including this sliding model in predictive ice-sheet models may improve their ability to predict time-evolving velocities and associated sea level change and reduce the related uncertainties.

Surface meltwater is becoming increasingly widespread on Antarctic ice shelves. It is stored within surface ponds and streams, or within firn pore spaces, which may saturate to form slush. Slush can reduce firn air content, increasing an ice-shelf's vulnerability to break-up. To date, no study has mapped the changing extent of slush across ice shelves. Here, we use Google Earth Engine and Landsat 8 images from six ice shelves to generate training classes using a k-means clustering algorithm, which are used to train a random forest classifier to identify both slush and ponded water. Validation using expert elicitation gives accuracies of 84% and 82% for the ponded water and slush classes, respectively. Errors result from subjectivity in identifying the ponded water/slush boundary, and from inclusion of cloud and shadows. We apply our classifier to the Roi Baudouin Ice Shelf for the entire 2013–20 Landsat 8 record. On average, 64% of all surface meltwater is classified as slush and 36% as ponded water. Total meltwater areal extent is greatest between late January and mid-February. This highlights the importance of mapping slush when studying surface meltwater on ice shelves. Future research will apply the classifier across all Antarctic ice shelves.