The southern Patagonian glaciers are known for having extremely high ablation rates. Foehn winds are one of the suspected causes, however, their influence on the annual ablation, their interannual variations, and their relationship with climate change is not well understood. We analysed the in-situ meteorological data from 2003–2020 recorded at Glaciar Perito Moreno. Daily temperature lapse rates varied substantially, from −7.8°C km−1 to 10.4°C km−1, due to foehn, fog, and katabatic winds. We find that, on average, foehn events occurred 1073 hours per year, and accounted for 20% of the annual surface ablation. This increase in surface ablation rates during foehn events occurs as a result of the enhanced sensible heat flux and net shortwave radiation. The downglacier-directed foehn winds warm the air mass over the glacier, but because of the high humidity of the foehn here, they often release latent heat by condensation. Variations in the Amundsen Sea Low influence foehn occurrence by modulating the westerly winds, which is related to the hemispherical ocean and atmospheric variability. Our results show that the local climate play an important role in the surface melting of Patagonian glaciers.

We study the meteorological drivers of melt at two glaciers in Taylor Valley, Antarctica, using 22 years of weather station observations and surface energy fluxes. The glaciers are located only 30 km apart, but have different local climates; Taylor Glacier is generally drier and windier than Commonwealth Glacier, which receives more snowfall due to its proximity to the coast. Commonwealth Glacier shows more inter-annual melt variability, explained by variable albedo due to summer snowfall events. A significant increase in surface melt at Commonwealth Glacier is associated with a decrease in summer minimum albedo. Inter-annual variability in melt at both glaciers is linked to degree-days above freezing during föhn events, occurring more frequently at Taylor Glacier. At Taylor Glacier melt occurs most often with positive air temperatures, but föhn conditions also favour sublimation, which cools the surface and prevents melt for the majority of the positive air temperatures. At Commonwealth Glacier, most of the melt instead occurs with sub-zero air temperatures, driven by strong solar radiative heating. Future melt at Taylor Glacier will likely be more sensitive to changes in föhn events, while Commonwealth Glacier will be impacted more by changes in near coastal weather, where moisture inputs can drive cloud cover, snowfall and change albedo.

Over half of the recent mass loss from the Greenland ice sheet, and its associated contribution to global sea level rise, can be attributed to increased surface meltwater runoff, with the remainder a result of dynamical processes such as calving and ice discharge. It is therefore important to quantify the distribution of melting on the ice sheet if we are to adequately understand past ice sheet change and make predictions for the future. In this article, we present a novel semi-empirical approach for characterising ice sheet surface conditions using high-resolution synthetic aperture radar (SAR) backscatter data from the Sentinel-1 satellite. We apply a state-space model to nine sites within North-East Greenland to identify changes in SAR backscatter, and we attribute these to different surface types with reference to optical satellite imagery and meteorological data. A set of decision-making rules for labelling ice sheet melting states are determined based on this analysis and subsequently applied to previously unseen sites. We show that our method performs well in (1) recognising some of the ice sheet surface types such as snow and dark ice and (2) determining whether the surface is melting or not melting. Sentinel-1 SAR data are of high spatial resolution; thus, in developing a method to identify the state of the surface from these data, we improve our capability to understand the variation of ice sheet melting across time and space.

Pigmented microalgae inhabiting snow and ice environments lower the albedo of glacier and ice-sheet surfaces, significantly enhancing surface melt. Our ability to accurately predict their role in glacier and ice-sheet surface mass balance is limited by the current lack of empirical data to constrain their representation in predictive models. Here we present new empirical optical properties for snow and ice algae and incorporate them in a radiative transfer model to investigate their impact on snow and ice surface albedo. We found ice algal cells to be more efficient absorbers than snow algal cells, but their blooms had comparable impact on surface albedo due to the different photic conditions of their habitats. We then used the model to reconstruct the effect of ice algae on bare ice albedo spectra collected at our field site in southern Greenland, where blooms dropped the albedo locally by between 3 and 43%, equivalent to 1–10 L m$^{-2}$ d$^{-1}$ of melted ice. Using the newly parametrized model, future studies could investigate biological albedo reduction and algal quantification from remote hyperspectral and multispectral imagery.

The marginal areas of the Greenland ice sheet develop streams and lakes each summer, documenting that surface runoff of meltwater is a major component of ice-sheet mass balance. Here we map the slush limit, a proxy for the extent of surface runoff, using daily MODIS data for the years 2000–2021. We develop an automated algorithm capable of detecting daily slush limits, provided sufficient image quality. The algorithm is applied to the ice sheet's western flank (61.7 $^{\circ }$N to 76.5 $^{\circ }$N). We find significant increasing trends in maximum slush limits until the year 2012, but not thereafter. We show that the slush limit typically rises quickly early in the ablation season but stabilizes before melting ceases. The data provide evidence that upward migration of surface runoff in summer 2012 stopped early at the upper margin of the ice slabs. These thick and continuous ice layers are located close to the surface, in the firn, and impede percolation of melt into deeper pore space. Had the ice slabs extended higher, the summer 2012 provided sufficient energy to raise the slush limit by another $\sim$300 m in elevation.

Increased surface melt in the percolation zone of the Greenland ice sheet causes significant changes in the firn structure, directly affecting the amount and timing of meltwater runoff. Here we force an energy-balance model with automatic weather stations data at two sites in the percolation zone of southwest Greenland ($2040$ and 2360 m a.s.l.) between spring $2017$ and fall $2019$. Extensive model validation and sensitivity analysis reveal that the skin layer formulation used to compute the surface temperature by closing the energy balance leads to a consistent overestimation of melt by more than a factor of two or three depending on the site. In contrast, model results match the observations well when the model is forced by observed surface temperatures; however, unexplained residuals in the energy balance occur. The sensible and ground heat flux differ markedly in the two simulations accounting largely for the difference in modeled melt amounts. This indicates that the energy available for melt is highly sensitive to small changes in surface temperature. Thus, regional climate models that also use the skin layer formulation may have a bias in surface temperature and melt energy in the percolation zone of the ice sheet.

Marine-terminating glaciers lose mass through melting and iceberg calving, and we find that meltwater drainage systems influence calving timing at Helheim Glacier, a tidewater glacier in East Greenland. Meltwater feeds a buoyant subglacial discharge plume at the terminus of Helheim Glacier, which rises along the glacial front and surfaces through the mélange. Here, we use high-resolution satellite and time-lapse imagery to observe the surface expression of this meltwater plume and how plume timing and location compare with that of calving and supraglacial meltwater pooling from 2011 to 2019. The plume consistently appeared at the central terminus even as the glacier advanced and retreated, fed by a well-established channelized drainage system with connections to supraglacial water. All full-thickness calving episodes, both tabular and non-tabular, were separated from the surfacing plume by either time or by space. We hypothesize that variability in subglacial hydrology and basal coupling drive this inverse relationship between subglacial discharge plumes and full-thickness calving. Surfacing plumes likely indicate a low-pressure subglacial drainage system and grounded terminus, while full-thickness calving occurrence reflects a terminus at or close to flotation. Our records of plume appearance and full-thickness calving therefore represent proxies for the grounding state of Helheim Glacier through time.

Surface melting on Amery Ice Shelf (AIS), East Antarctica, produces an extensive supraglacial drainage system consisting of hundreds of lakes connected by surface channels. This drainage system forms most summers on the southern portion of AIS, transporting meltwater large distances northward, toward the ice front and terminating in lakes. Here we use satellite imagery, Landsat (1, 4 and 8), MODIS multispectral and Sentinel-1 synthetic aperture radar to examine the seasonal and interannual evolution of the drainage system over nearly five decades (1972–2019). We estimate seasonal meltwater input to one lake by integrating output from the regional climate model [Regional Atmospheric Climate Model (RACMO 2.3p2)] over its catchment defined using the Reference Elevation Model of Antarctica. We find only weak positive relationships between modeled seasonal meltwater input and lake area and between meltwater input and lake volume. Consecutive years of extensive melting lead to year-on-year expansion of the drainage system, potentially through a link between melt production, refreezing in firn and the maximum extent of the lakes at the downstream termini of drainage. These mechanisms are important when evaluating the potential of drainage systems to grow in response to increased melting, delivering meltwater to areas of ice shelves vulnerable to hydrofracture.

The weakening and/or removal of floating ice shelves in Antarctica can induce inland ice flow acceleration. Numerical modelling suggests these processes will play an important role in Antarctica's future sea-level contribution, but our understanding of the mechanisms that lead to ice tongue/shelf collapse is incomplete and largely based on observations from the Antarctic Peninsula and West Antarctica. Here, we use remote sensing of structural glaciology and ice velocity from 2001 to 2020 and analyse potential ocean-climate forcings to identify mechanisms that triggered the rapid disintegration of ~2445 km2 of ice mélange and part of the Voyeykov Ice Shelf in Wilkes Land, East Antarctica between 27 March and 28 May 2007. Results show disaggregation was pre-conditioned by weakening of the ice tongue's structural integrity and was triggered by mélange removal driven by a regional atmospheric circulation anomaly and a less extensive latent-heat polynya. Disaggregation did not induce inland ice flow acceleration, but our observations highlight an important mechanism through which floating termini can be removed, whereby the break-out of mélange and multiyear landfast sea ice triggers disaggregation of a structurally-weak ice shelf. These observations highlight the need for numerical ice-sheet models to account for interactions between sea-ice, mélange and ice shelves.

The Kaskawulsh Glacier is an iconic outlet draining the icefields of the St. Elias Mountains in Yukon, Canada. We determine and attempt to interpret its catchment-wide mass budget since 2007. Using SPOT5/6/7 data we estimate a 2007–18 geodetic balance of −0.46 ± 0.17 m w.e. a−1. We then compute balance fluxes and observed ice fluxes at nine flux gates to examine the discrepancy between the climatic mass balance and internal mass redistribution by glacier flow. Balance fluxes are computed using a fully distributed mass-balance model driven by downscaled and bias-corrected climate-reanalysis data. Observed fluxes are calculated using NASA ITS_LIVE surface velocities and glacier cross-sectional areas derived from ice-penetrating radar data. We find the glacier is still in the early stages of dynamic adjustment to its mass imbalance. We estimate a committed terminus retreat of ~23 km under the 2007–18 climate and a lower bound of 46 km3 of committed ice loss, equivalent to ~15% of the total glacier volume.

Close-range (CR) L-band radiometry and quasi-simultaneous in situ snow characterizations were conducted in May 2019 at the Swiss Camp research site in the ablation zone of the western Greenland ice sheet. Snow liquid-water and its melt/refreeze are retrieved from microwave antenna temperatures measured with the ground-based L-band radiometer ELBARA-III. The emission model (EM) used in the retrieval is a two-layer configuration of the ‘L-Band Specific Microwave Emission Model of Layered Snowpack’. Consistent snow wetness retrievals were achieved from both single- and multi-angle CR observations of L-band antenna temperatures. This suggests that multi-angle observation is not a pre-requisite for snow wetness retrieval. Therefore, in addition to soil moisture and ocean salinity (SMOS) multi-angle measurements, snow wetness can be estimated from spaceborne L-band brightness temperatures measured at a single observation angle, such as from NASA's SMAP satellite. Our results provide partial validation of a recently presented snow wetness retrieval approach based on the same EM and applied over Greenland using multi-angle SMOS brightness temperatures. Agreement between measured CR antenna temperatures and SMOS brightness temperatures is found to be within the 95% confidence intervals of ELBARA-III and SMOS measurement uncertainties. Our measurements confirm the modeled response of antenna temperatures to diurnal variations of snow wetness.

Supraglacial debris is significant in many regions and complicates modeling of glacier melt, which is required for predicting glacier change and its influences on hydrology and sea-level rise. Temperature-index models are a popular alternative to energy-balance models when forcing data are limited, but their transferability among glaciers and inherent uncertainty have not been documented in application to debris-covered glaciers. Here, melt factors were compiled directly from published studies or computed from reported melt and MERRA-2 air temperature for 27 debris-covered glaciers around the world. Linear mixed-effects models were fit to predict melt factors from debris thickness and variables including debris lithology and MERRA-2 radiative exchange. The models were tested by leave-one-site-out cross-validation based on predicted melt rates. The best model included debris thickness (fixed effect) and glacier and year (random effects). Predictions were more accurate using MERRA-2 than on-site air temperature data, and pooling MERRA-2-derived and reported melt factors improved cross-validation accuracy more than including additional predictors such as shortwave or longwave radiation. At one glacier where monthly ablation was measured over 4 years, seasonal variation of melt factors suggested that heat storage significantly affected the relation between melt and energy exchange at the debris surface.

We present a newly developed 1-D numerical energy-balance and phase transition supraglacial lake model: GlacierLake. GlacierLake incorporates snowfall, in situ snow and ice melt, incoming water from the surrounding catchment, ice lid formation, basal freeze-up and thermal stratification. Snow cover and temperature are varied to test lake development through winter and the maximum lid thickness is recorded. Average wintertime temperatures of −2 to $-30^{\circ }{\rm C}$ and total snowfall of 0 to 3.45 m lead to a range of the maximum lid thickness from 1.2 to 2.8 m after ${\sim }250$ days, with snow cover exerting the dominant control. An initial ice temperature of $-15^{\circ }{\rm C}$ with simulated advection of cold ice from upstream results in 0.6 m of basal freeze-up. This suggests that lakes with water depths above 1.3 to 3.4 m (dependent on winter snowfall and temperature) upon lid formation will persist through winter. These buried lakes can provide a sizeable water store at the start of the melt season, expedite future lake formation and warm underlying ice even in winter.

We examine the spatial patterns of near-surface air temperature (Ta) over a melting glacier using a multi-annual dataset from McCall Glacier, Alaska. The dataset consists of a 10-year (2005–2014) meteorological record along the glacier centreline up to an upper glacier cirque, spanning an elevation difference of 900 m. We test the validity of on-glacier linear lapse rates, and a model that calculates Ta based on the influence of katabatic winds and other heat sources along the glacier flow line. During the coldest hours of each summer (10% of time), average lapse rates across the entire glacier range from −4.7 to −6.7°C km−1, with a strong relationship between Ta and elevation (R2 > 0.7). During warm conditions, Ta shows more complex, non-linear patterns that are better explained by the flow line-dependent model, reducing errors by up to 0.5°C compared with linear lapse rates, although more uncertainty might be associated with these observations due to occasionally poor sensor ventilation. We conclude that Ta spatial distribution can vary significantly from year to year, and from one glacier section to another. Importantly, extrapolations using linear lapse rates from the ablation zone might lead to large underestimations of Ta on the upper glacier areas.

Surface melt on the coastal Antarctic ice sheet (AIS) determines the viability of its ice shelves and the stability of the grounded ice sheet, but very few in situ melt rate estimates exist to date. Here we present a benchmark dataset of in situ surface melt rates and energy balance from nine sites in the eastern Antarctic Peninsula (AP) and coastal Dronning Maud Land (DML), East Antarctica, seven of which are located on AIS ice shelves. Meteorological time series from eight automatic and one staffed weather station (Neumayer), ranging in length from 15 months to almost 24 years, serve as input for an energy-balance model to obtain consistent surface melt rates and energy-balance results. We find that surface melt rates exhibit large temporal, spatial and process variability. Intermittent summer melt in coastal DML is primarily driven by absorption of shortwave radiation, while non-summer melt events in the eastern AP occur during föhn events that force a large downward directed turbulent flux of sensible heat. We use the in situ surface melt rate dataset to evaluate melt rates from the regional atmospheric climate model RACMO2 and validate a melt product from the QuikSCAT satellite.

Glaciological ablation is computed from point-scale data at a few ablation stakes that are usually regressed as a function of elevation and averaged over the area-elevation distribution of a glacier. This method is contingent on a tight control of elevation on local ablation. However, in debris-covered glaciers, systematic and random spatial variations of debris thickness modify the ablation rates. We propose and test a method to compute sub-debris ablation where stake data are interpolated as a function of debris-thickness alone and averaged over the debris-thickness distribution at different parts of the glacier. We apply this method on Satopanth Glacier located in Central Himalaya utilising ~1000 ablation measurements obtained from a network of up to 56 stakes during 2015–2017. The estimated mean sub-debris ablation ranges between 1.5±0.2 to 1.7±0.3 cm d−1. We show that the debris-thickness-dependent regression describes the spatial variability of the sub-debris ablation better than the elevation dependent regression. The uncertainties in ablation estimates due to the corresponding uncertainties in the measurement of ablation and debris-thickness distribution, and those due to interpolation procedures are estimated using Monte Carlo methods. Possible biases due to a finite number of stakes used are also investigated.

Ice cliffs and ponds on debris-covered glaciers have received increased attention due to their role in amplifying local melt. However, very few studies have looked at these features on the catchment scale to determine their patterns and changes in space and time. We have compiled a detailed inventory of cliffs and ponds in the Langtang catchment, central Himalaya, from six high-resolution satellite orthoimages and DEMs between 2006 and 2015, and a historic orthophoto from 1974. Cliffs cover between 1.4% (± 0.4%) in the dry and 3.4% (± 0.9%) in the wet seasons and ponds between 0.6% (± 0.1%) and 1.6% (± 0.3%) of the total debris-covered tongues. We find large variations between seasons, as cliffs and ponds tend to grow in the wetter monsoon period, but there is no obvious trend in total area over the study period. The inventory further shows that cliffs are predominately north-facing irrespective of the glacier flow direction. Both cliffs and ponds appear in higher densities several hundred metres from the terminus in areas where tributaries reach the main glacier tongue. On the largest glacier in the catchment ~10% of all cliffs and ponds persisted over nearly a decade.

To understand the geological origins of minerals in cryoconite and the nutrients sources for microbes on glaciers, we analyzed the Sr–Nd isotopic ratios of the four mineral fractions in cryoconites including saline, carbonate, phosphate, silicate and the organic fraction obtained from Gulkana Glacier in Alaska. The isotopes in the silicate mineral fraction exhibited spatial variation within the glacier (87Sr/86Sr: 0.704533–0.709563, εNd (0): −16.0 to 0.5), which can be explained by the different mixing ratios of the two distinct sources: one of the sources is lateral and terminal moraines or soil, and the other is the medial moraine of the glacier. The minerals in the cryoconite at the lower sites in the glacier are likely derived from the former source, whereas those at the upper sites are from latter sources. The mineralogical and elemental compositions also support mixing of the silicate minerals from the two local sources. The Sr isotopic ratios of the organic fraction also showed spatial variation on the glacier in the middle sites – a trend similar to those of the phosphate fraction. The results suggest that the organic matter is mostly the byproducts of microbes using the phosphate minerals as a nutrient source.

A geochemical method to characterize post-depositional melting and elution is demonstrated using inductively coupled plasma mass spectrometry to measure concentrations of dissolved and insoluble fractions of major crustal elements in snow samples collected from March 2006 to January 2010 at Urumqi Glacier No. 1, Tien Shan. Dust from these samples has compositional homogeneity, suggesting that dust has a stable dissolved fraction percentage (DFP, calculated as dissolved/(dissolved + insoluble)%). Calcium has the highest DFP (averaging 61.5 ± 19.4%), followed by Na (30.4 ± 19.6%), Mg (13.2 ± 9.8%), and K (7.9 ± 9.8%). Acid input can affect dissolution of Na and Ca. Taking DFP values of unmelted samples as the reference, the higher DFPs refer to strengthened dissolution from acid input, while the lower ones refer to elution. Based on the DFP difference between unmelted and eluted states, an elution sequence Ca > Na > Mg > K is obtained. Some details such as the beginning and the ending stages of elution can be found by DFP and acid input index, while using ion concentration is not capable of this. Our results reveal that acid input is an important mechanism for DFP changes, that the DFP index can provide an effective assessment of snow elution, and that this will aid in understanding low latitude ice cores.

The detection and monitoring of meltwater within firn presents a significant monitoring challenge. We explore the potential of small wireless sensors (ETracer+, ET+) to measure temperature, pressure, electrical conductivity and thus the presence or absence of meltwater within firn, through tests in the dry snow zone at the East Greenland Ice Core Project site. The tested sensor platforms are small, robust and low cost, and communicate data via a VHF radio link to surface receivers. The sensors were deployed in low-temperature firn at the centre and shear margins of an ice stream for 4 weeks, and a ‘bucket experiment’ was used to test the detection of water within otherwise dry firn. The tests showed the ET+ could log subsurface temperatures and transmit the recorded data through up to 150 m dry firn. Two VHF receivers were tested: an autonomous phase-sensitive radio-echo sounder (ApRES) and a WinRadio. The ApRES can combine high-resolution imaging of the firn layers (by radio-echo sounding) with in situ measurements from the sensors, to build up a high spatial and temporal resolution picture of the subsurface. These results indicate that wireless sensors have great potential for long-term monitoring of firn processes.