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.

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.

Accurately quantifying all the components of the surface energy balance (SEB) is a prerequisite for the reliable estimation of surface melt and the surface mass balance over ice and snow. This study quantifies the SEB closure by comparing the energy available for surface melt, determined from continuous measurements of radiative fluxes and turbulent heat fluxes, to the surface ablation measured on the Greenland ice sheet between 2003 and 2023. We find that the measured daily energy available for surface melt exceeds the observed surface melt by on average 18 ± 30 W m−2 for snow and 12 ± 54 W m−2 for ice conditions (mean ± SD), which corresponds to 46 and 10% of the average energy available for surface melt, respectively. When the surface is not melting, the daily SEB is on average closed within 5 W m−2. Based on the inter-comparison of different ablation sensors and radiometers installed on different stations, and on the evaluation of modelled turbulent heat fluxes, we conclude that measurement uncertainties prevent a better daily to sub-daily SEB closure. These results highlight the need and challenges in obtaining accurate long-term in situ SEB observations for the proper evaluation of climate models and for the validation of remote sensing products.

Global Navigation Satellite System (GNSS) observations and ground-based timelapse photography obtained over the record-high 2019/2020 melt season are combined to characterise the flexure and fracture behaviour of a previously formed doline on George VI Ice Shelf, Antarctica. The GNSS timeseries shows a downward vertical displacement of the doline centre with respect to the doline rim of ~60 cm in response to loading from a central meltwater lake. The GNSS data also show a tens-of-days episode of rapid-onset, exponentially decaying horizontal displacement, where the horizontal distance between the doline rim and its centre increases by ~70 cm. We interpret this event as the initiation and/or widening of a fracture, aided by stress perturbations associated with meltwater loading in the doline basin. Viscous flexure modelling indicates that the meltwater loading generates tensile surface stresses exceeding 75 kPa. This, together with our timelapse photos of circular fractures around the doline, suggests the first such documentation of meltwater-loading-induced ‘ring fracture’ formation on an ice shelf, equivalent to the fracture type proposed as part of the chain-reaction lake drainage process involved in the 2002 breakup of the Larsen B Ice Shelf.

New experiments have revealed that the addition of a small amount of wood ash doubles the angle of repose of granulated ice at −10°C, but has little effect at −30°C. It is hypothesized that this behavior may be explained in terms of the freezing of water that is formed at the higher temperature through the reaction of salt within the ash.

As warming intensifies across the Greenland ice sheet, an increasing number of shallow coring and radar studies are targeting the melt-impacted firn column to investigate meltwater processes. Highly inhomogeneous infiltration and refreezing, however, redistributes mass, distorting age–depth relationships and confounding comparisons between different cores. Here, we utilize a dynamic time warping algorithm for time series alignment to investigate and quantify the heterogeneous impact of melt processes on nine closely spaced (within 50 m) firn core-density profiles. The 10 m cores were collected relatively high in Greenland's percolation zone, where melt alteration is minimal compared to lower elevation. Our analysis demonstrates the effectiveness of dynamic time warping as a tool for assessing heterogeneity between ice core records. We find that the optimal alignment of density profiles in the nine cores requires vertical stretching and compression of individual profiles, ranging from, on average, <1 to ~16% of the core lengths. We identified four depth zones of mass redistribution that appear to coincide with observed ice layers. Further, ~75% of density measurements from each core do not align with an age model-derived density profile that assumes no mass redistribution of meltwater, indicating the pervasive impact of melt processes.

Glide-snow avalanches release due to a loss of friction at the snow–ground interface, which can result in large avalanches that endanger infrastructure in alpine regions. The underlying processes are still relatively poorly understood, in part due to the limited data available on glide processes. Here, we introduce a pixel-based algorithm to detect glide cracks in time-lapse photographs under changing illumination and shadow conditions. The algorithm was applied to 14 years of time-lapse images at Dorfberg (Davos, Switzerland). We analysed 947 glide-snow events at a high-spatial (0.5 m) and temporal (2–15 min) resolution. Avalanche activity and glide-crack opening dynamics were investigated across timescales ranging from seasonally to hourly. Events were separated into surface (meltwater percolation) and interface events (no meltwater percolation). The results show that glide activity is highly variable between and within seasons. Most avalanches released without crack formation or within 24 h after crack opening, and release was favoured in the afternoon hours. Glide rates often showed a stick–slip behaviour. The acceleration of glide rates and non-constant increases in glide crack aspect ratio were indicators for avalanche release. This comprehensive dataset provides the basis for further investigations into glide-snow avalanche drivers.

Melting glacier surfaces are unique ecosystems for specialized microbes, frequently harbouring blooms of microalgae with pigments contributing to the darkening of ice surfaces, reducing albedo and enhancing melt rates. The main cause of this phenomenon is algae of the genus Ancylonema. Prior investigation depended on field-collected material because these algae resisted cultivation. To enhance research on how these algae dominate melting ice, we established a strain of Ancylonema alaskanum from an alpine glacier and exposed to temperatures around the freezing point at irradiations of ~10% of full sunlight. The morphology of the culture changed, with the cells becoming longer and turning green by losing their brownish pigmentation, indicating that these dark phenols are crucial for survival in the cryosphere. Photophysiological comparisons of strain and glacial material showed adaptation of the photosynthetic apparatus to prevailing conditions. This laboratorial strain opens possibilities for a wide range of comparative ‘omics’ research.

Supraglacial debris cover regulates the melt rates of many glaciers in mountainous regions around the world, thereby modifying the availability and quality of downstream water resources. However, the influence of supraglacial debris is often poorly represented within glaciological models, due to the absence of a technique to provide high-precision, spatially continuous measurements of debris thickness. Here, we use high-resolution UAV-derived thermal imagery, in conjunction with local meteorological data, visible UAV imagery and vertically profiled debris temperature time series, to model the spatially distributed debris thickness across a portion of Llaca Glacier in the Cordillera Blanca of Peru. Based on our results, we simulate daily sub-debris melt rates over a 3-month period during 2019. We demonstrate that, by effectively calibrating the radiometric thermal imagery and accounting for temporal and spatial variations in meteorological variables during UAV surveys, thermal UAV data can be used to more precisely represent the highly heterogeneous patterns of debris thickness and sub-debris melt on debris-covered glaciers. Additionally, our results indicate a mean sub-debris melt rate nearly three times greater than the mean melt rate simulated from satellite-derived debris thicknesses, emphasising the importance of acquiring further high-precision debris thickness data for the purposes of investigating glacier-scale melt processes, calibrating regional melt models and improving the accuracy of runoff predictions.

A conceptual model, based on field observations and assumed physics of a perennial firn aquifer near Helheim Glacier (southeast Greenland), is evaluated via steady-state 2-D simulation of liquid water flow and energy transport with phase change. The simulation approach allows natural representation of flow and energy advection and conduction that occur in vertical meltwater recharge through the unsaturated zone and in lateral flow within the saturated aquifer. Agreement between measured and simulated aquifer geometry, temperature, and recharge and discharge rates confirms that the conceptual field-data-based description of the aquifer is consistent with the primary physical processes of groundwater flow, energy transport and phase change. Factors that are found to control simulated aquifer configuration include surface temperature, meltwater recharge rate, residual total-water saturation and capillary fringe thickness. Simulation analyses indicate that the size of perennial firn aquifers depends primarily on recharge rates from surface snowmelt. Results also imply that the recent aquifer expansion, likely due to a warming climate, may eventually produce lakes on the ice-sheet surface that would affect the surface energy balance.

A supraglacial debris layer controls energy transfer to the ice surface and moderates ice ablation on debris-covered glaciers. Measurements of vertical temperature profiles within the debris enables the estimation of thermal diffusivities and sub-debris ablation rates. We have measured the debris-layer temperature profiles at 16 locations on Satopanth Glacier (central Himalaya) during the ablation seasons of 2016 and 2017. Debris temperature profile data are typically analysed using a finite-difference method, assuming that the debris layer is a homogeneous one-dimensional thermal conductor. We introduce three more methods for analysing such data that approximate the debris layer as either a single or a two-layered conductor. We analyse the performance of all four methods using synthetic experiments and by comparing the estimated ablation rates with in situ glaciological observations. Our analysis shows that the temperature measurements obtained at equispaced sensors and analysed with a two-layered model improve the accuracy of the estimated thermal diffusivity and sub-debris ablation rate. The accuracy of the ablation rate estimates is comparable to that of the in situ observations. We argue that measuring the temperature profile is a convenient and reliable method to estimate seasonal to sub-seasonal variations of ablation rates in the thickly debris-covered parts of glaciers.

The intermixed thermal and structural framework of cold firn, water-saturated firn and ice layers in Greenland's percolation zone can be challenging to penetrate with core drills. Here, we present our experiences using a hot water drill for research on the firn layer of the percolation zone. We built and deployed a lightweight and easily transportable system for drilling a transect of ~15 cm diameter boreholes through the full firn column thickness, to depths exceeding 100 m. An instrumented drill stem provides a scientific measurement of the firn properties while drilling. The system was successful at gaining rapid access to the firn column with mixed wet and cold conditions, was easily transported to the site and across the glacier surface, and required a small field crew to operate. The boreholes are well suited for in situ investigations of firn processes in Greenland percolation zone.

Hydrofracturing can enhance the depth to which crevasses propagate and, in some cases, allow full depth crevasse penetration and iceberg detachment. However, many existing crevasse models either do not fully account for the stress field driving the hydrofracture process and/or treat glacier ice as elastic, neglecting the non-linear viscous rheology. Here, we present a non-local continuum poro-damage mechanics (CPDM) model for hydrofracturing and implement it within a full Stokes finite element formulation. We use the CPDM model to simulate the propagation of water-filled crevasses in idealized grounded glaciers, and compare crevasse depths predicted by this model with those from linear elastic fracture mechanics (LEFM) and zero stress models. We find that the CPDM model is in good agreement with the LEFM model for isolated crevasses and with the zero stress model for closely-spaced crevasses, until the glacier approaches buoyancy. When the glacier approaches buoyancy, we find that the CPDM model does not allow the propagation of water-filled crevasses due to the much smaller size of the tensile stress region concentrated near the crevasse tip. Our study suggests that the combination of non-linear viscous and damage processes in ice near the tip of a water-filled crevasse can alter calving outcomes.

The retention of meltwater in the accumulation area of the Greenland ice sheet and other Arctic ice masses buffers their contribution to sea level change. However, sustained warming also results in impermeable ice layers or ‘ice slabs’ that seal the underlying pore space. Here, we use a 1-D, physically based, high-resolution model to simulate the surface mass balance (SMB), percolation, refreezing, ice layer formation and runoff from across the high-elevation area of Devon Ice Cap, Canada, from 2001 to 2016. We vary the thickness of the ‘impermeable’ ice layer at which underlying firn becomes inaccessible to meltwater. Thick near-surface ice layers are established by an initial deep percolation, the formation of decimetre ice layers and the infilling of interleaving pore space. The cumulative SMB increases by 48% by varying impermeable layer thickness between 0.01 and 5 m. Within this range we identify narrower range (0.25–1 m) that can simulate both the temporal variability in SMB and the observed near-surface density structure. Across this range, cumulative SMB variation is limited to 6% and 45–49% of mass retention takes place within the annually replenished snowpack. Our results indicate cooler summers after intense mid-2000s warming have led to a partial replenishment of pore space.

The positive degree-day (PDD) model provides a particularly simple approach to estimate surface melt from land ice based solely on air temperature. Here, we use a climate and snow pack simulation of the Greenland ice sheet (Modèle Atmosphérique Régional, MAR) as a reference, to analyze this scheme in three realizations that incorporate the sub-monthly temperature variability differently: (i) by local values, (ii) by local values that systematically overestimate the dampened variability associated with intense melting or (iii) by one constant value. Local calibrations reveal that incorporating local temperature variability, particularly resolving the dampened variability of melt areas, renders model parameters more temperature-dependent. This indicates that the negative feedback between surface melt and temperature variability introduces a non-linearity into the temperature – melt relation. To assess the skill of the individual realizations, we hindcast melt rates from MAR temperatures for each realization. For this purpose, we globally calibrate Greenland-wide, constant parameters. Realization (i) exhibits shortcomings in the spatial representation of surface melt unless temperature-dependent instead of constant parameters are calibrated. The other realizations perform comparatively well with constant parametrizations. The skill of the PDD model primarily depends, however, on the consistent calibration rather than on the specific representation of variability.

This paper reports a comparative analysis performed on a fraction-image time series of the Antarctic Peninsula from the period 1999–2009 generated by multiresolution remote-sensing images (SSM/I and SSMI/S with 25 km and QuikSCAT with 2.225 km spatial resolutions) for snow-melt detection. Our method is based on the (a) preprocessing of multitemporal remote-sensing data, (b) subpixel mixture analysis of SSMI and QuikSCAT image time series, and (c) evaluation of subpixel analysis, including an assessment of fraction images of wet snow using an independent ASAR dataset and sensitivity analysis on the melt metrics measured by these images. The temporal dynamics of the melt indices derived from the wet-snow fraction images presented a more realistic pattern than the traditional melt metrics measured by Boolean snow-melt detection approaches. Because the snow melt actually occurs at the pixel fractions, the multiscale analysis that was performed suggests an overestimation of the melt metrics calculated using Boolean approaches (which assume that the entire area of the detected pixel shows snow melt). The melt metrics measurements show an overestimation according to the decrease in spatial resolution related to the multiplicative effect of a larger pixel area.

Tropical glaciers constitute an important source of water for downstream populations. However, our understanding of glacial melt processes is still limited. One observed process that has not yet been quantified for tropical glaciers is the enhanced melt caused by the longwave emission transfer. Here, we use high-resolution surface temperatures obtained from the thermal infrared imagery of the Cuchillacocha Glacier, in the Cordillera Blanca, Peru in June 2014 to calculate a margin longwave flux. This longwave flux, reaching the glacier margin from the adjacent exposed rock, varies between 81 and 120 W m−2 daily. This flux is incorporated into a physically-based melt model to assess the net radiation budget at the modeled glacier margin. The simulation results show an increase in the energy available for melt by an average of 106 W m−2 during the day when compared with the simulation where the LWmargin flux is not accounted for. This value represents an increase in ablation of ~1.7 m at the glacier margin for the duration of the dry season. This study suggests that including the quantification of the glacier margin longwave flux in physically-based melt models results in an improved assessment of tropical glacier energy budget and meltwater generation.

Surface melt in the ablation zone is dominated by atmospheric temperature and surface albedo. We developed a surface mass-balance model with a dynamic component of glacier ice albedo which includes surface properties, clouds and the angle of the sun. The ice albedo reduction is mainly caused by impurity accumulation of non-biological origin such as dust and black carbon (BC), which is currently not included in other surface mass-balance models. Simulations show that dust from meltout is the main source of impurity mass at the melting glacier ice surface, and current rates of atmospheric deposition of dust play only a minor role. However, for BC the atmospheric deposition is the main source where ice melt rates are below 1 m, and atmospheric deposition is most likely from intercontinental transport due to the scarce population and lack of forests in Greenland.

The inland advance of supraglacial lakes (SGLs) towards the interior regions of the Greenland ice sheet (GrIS) may have implications for the water volumes reaching the subglacial drainage system, and could consequently affect long-term ice-sheet dynamics. Here, we investigate changes to the areas, volumes and elevation distributions of over 8000 manually delineated SGLs using 44 Landsat images of a 6200 km2 sector of north-west Greenland over three decades (1985–2016). Our results show that SGLs have advanced to higher maximum (+418 m) and mean (+299 m) elevations, and that there has been a near-doubling of total regional SGL areas and volumes over the study period, accelerating after 2000. These changes were primarily caused by an increased SGL area and volume at high (≥1200 m a.s.l.) elevations, where SGL coverage increased by over 2750% during the study period. Many of the observed changes, particularly the post-2000 accelerations, were driven by changes to regional surface-temperature anomalies. This study demonstrates the past and accelerating response of the GrIS's hydrological system due to climatic warming, indicating an urgent need to understand whether the increasingly inland SGLs will be capable of hydrofracture in the future, thus determining their potential implications for ice-sheet dynamics.

Understanding the climate response of the Antarctic Peninsula ice sheet is vital for accurate predictions of sea-level rise. However, since climate models are typically too coarse to capture spatial variability in local scale meteorological processes, our ability to study specific sectors has been limited by the local fidelity of such models and the (often sparse) availability of observations. We show that a high-resolution (5.5 km × 5.5 km) version of a regional climate model (RACMO2.3) can reproduce observed interannual variability in the Larsen B embayment sufficiently to enable its use in investigating long-term changes in this sector. Using the model, together with automatic weather station data, we confirm previous findings that the year of the Larsen B ice shelf collapse (2001/02) was a strong melt year, but discover that total annual melt production was in fact ~30% lower than 2 years prior. While the year before collapse exhibited the lowest melting and highest snowfall during 1980–2014, the ice shelf was likely pre-conditioned for collapse by a series of strong melt years in the 1990s. Melt energy has since returned to pre-1990s levels, which likely explains the lack of further significant collapse in the region (e.g. of SCAR Inlet).