Bulk aerodynamic methods have been shown to perform poorly in computing turbulent heat fluxes at glacier surfaces during shallow katabatic winds. Katabatic surface layers have different wind shear and flux profiles to the surface layers for which the bulk methods were developed, potentially invalidating their use in these conditions. In addition, eddy covariance-derived turbulent heat fluxes are unlikely to be representative of surface conditions when eddy covariance data are collected close to the wind speed maximum (WSM). Here we utilize two months of eddy covariance and meteorological data measured at three different heights (1 m, 2 m, and 3 m) at Kaskawulsh Glacier in the Yukon, Canada, to re-examine the performance of bulk methods relative to eddy covariance-derived fluxes under different near-surface flow regimes. We propose a new set of processing methods for one-level eddy covariance data to ensure the validity of calculated fluxes during highly variable flows and low-level wind speed maxima, which leads to improved agreement between eddy covariance-derived and modelled fluxes across all flow regimes, with the best agreement (correlation >0.9) 1 m above the surface. Contrary to previous studies, these results show that adequately processed eddy covariance data collected at or above the WSM can provide valid estimates of surface heat fluxes.

Glaciers cover 132 900 km2 around the Antarctic Ice Sheet, but few are subject to annual mass-balance measurements. Lookalike Glacier and Davies Dome on James Ross Island have been monitored since 2009, providing the third longest mass-balance record for the northern Antarctic Peninsula. These glaciers had a balanced mass budget over the period 2009/10–2014/15 but started to lose their mass thereafter. Between 2014/15 and 2020/21, mass change rates were −0.15 ± 0.13 and −0.26 ± 0.11 m w.e. a−1 for Lookalike Glacier and Davies Dome, respectively. The mean equilibrium-line altitudes over this period at Lookalike Glacier (362 ± 18 m a.s.l.) and Davies Dome (>427 ± 22 m a.s.l.) are 51 and >34 m higher compared to the previous 6-year period. The mean accumulation area ratio values determined for the period 2014/15–2020/21 are lower than the balanced-budget ratio indicating that glaciers are out of balance with the current climate. The data confirm the transition from positive to negative mass-balance periods around 2014/15, which is attributed to the change in air temperature trends. The mean summer temperature increased by 0.9°C between the periods 2009/10–2014/15 and 2015/16–2020/21 and melt-season temperatures became predominantly positive.

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.

The retreat rates of Triangular Glacier since 1979 and its mass changes during the period 2014/15–2019/20 indicate the sensitive response of small ice masses on the eastern side of the Antarctic Peninsula to air temperature evolution. This cirque glacier in the northern part of James Ross Island receded rapidly during the period of regional warming in the late 20th century, losing 30.8% of its surface area between 1979 and 2006 (−1.7% a−1). The retreat rate then dropped to −0.3% a−1 following the regional cooling trend, but started to accelerate again (−0.8 to −2.3% a−1) with increasing air temperature since the summer 2014/15. Since the glaciological year 2015/16, Triangular Glacier has experienced enhanced snow melt, wind scour and permanent mass loss with annual mass balance ranging from −0.08 ± 0.35 to −0.56 ± 0.25 m w.e. The largest mass loss was observed in the glaciological year 2019/20, which included the warmest summer of the observation period. The cumulative mass balance of −1.66 ± 0.83 m w.e. over the years 2014/15–2019/20 is consistent with the termination of the positive mass-balance period that occurred in the north-eastern Antarctic Peninsula from 2009/10 to 2014/15.

Annual balances of eight alpine glaciers were slightly negative for 1961–90 and highly negative for 1991–2018. We explain this by changes in positive degree-day sums and summer temperatures extrapolated to the median altitudes of the glaciers. We test a new way of calculating degree-day sums that performs better than the traditional method which used daily mean temperatures. Annual degree-day sums are highly correlated with May–September temperatures as suggested in 1866 by Karl von Sonklar. We find moderate correlations between annual balances and degree-day sums, and with May–September temperatures. Calculated degree-day factors for the eight glaciers cover the reported range for snow and ice ablation, while the temperature sensitivity of annual balance is from −0.4 to −1.0 m w.e. for a +1°C temperature change. We accurately predict mean balances for 1991–2018 using May–September temperatures in regression models calibrated for 1961–90. May–September temperatures in the Alps have already increased ~+3°C since 1880 and, if temperatures continue to rise, these glaciers will shrink rapidly. As annual balances are already negative for present-day temperatures, these glaciers will not be ‘safe’ under the further temperature increase permitted by the Paris Agreement.

We present the first direct comparison of turbulence conditions measured simultaneously over exposed ice and a 0.08 m thick supraglacial debris cover on Suldenferner, a small glacier in the Italian Alps. Surface roughness, sensible heat fluxes (~20–50 W m−2), latent heat fluxes (~2–10 W m−2), topology and scale of turbulence are similar over both glacier surface types during katabatic and synoptically disturbed conditions. Exceptions are sunny days when buoyant convection becomes significant over debris-covered ice (sensible heat flux ~ −100 W m−2; latent heat flux ~ −30 W m−2) and prevailing katabatic conditions are rapidly broken down even over this thin debris cover. The similarity in turbulent properties implies that both surface types can be treated the same in terms of boundary layer similarity theory. The differences in turbulence between the two surface types on this glacier are dominated by the radiative and thermal contrasts, thus during sunny days debris cover alters both the local surface turbulent energy fluxes and the glacier component of valley circulation. These variations under different flow conditions should be accounted for when distributing temperature fields for modeling applications over partially debris-covered glaciers.

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.

In this study we analyse a 15-year long time series of surface mass-balance (SMB) measurements performed between 2001 and 2016 in the ablation zone of the Morteratsch glacier complex (Engadine, Switzerland). For a better understanding of the SMB variability and its causes, multiple linear regressions analyses are performed with temperature and precipitation series from nearby meteorological stations. Up to 85% of the observed SMB variance can be explained by the mean May–June–July temperature and the total precipitation from October to March. A new method is presented where the contribution of each month's individual temperature and precipitation to the SMB can be examined in a total sample of 224 (16.8 million) combinations. More than 90% of the observed SMB can be explained with particular combinations, in which the May–June–July temperature is the most recurrent, followed by October temperature. The role of precipitation is less pronounced, but autumn, winter and spring precipitation are always more important than summer precipitation. Our results indicate that the length of the ice ablation season is of larger importance than its intensity to explain year-to-year variations. The widely used June–July–August temperature index may not always be the best option to describe SMB variability through statistical correlation.

The MODerate resolution Imaging Spectroradiometer (MODIS) land surface temperature (LST) data have been widely used for air temperature estimation in mountainous regions where station observations are sparse. However, the performance of MODIS LST in high-elevation glacierized areas remains unclear. This study investigates air temperature estimation in glacierized areas based on ground observations at four glaciers across the Tibetan Plateau. Before being used to estimate the air temperature, MODIS LST data are evaluated at two of the glaciers, which indicates that MODIS night-time LST is more reliable than MODIS daytime LST data. Then, linear models based on each of the individual MODIS LST products from two platforms (Terra and Aqua) and two overpasses (night-time and daytime) are built to estimate daily mean, minimum and maximum air temperatures in glacierized areas. Regional glacier surface (RGS) models (mean /-mean-square differences: 3.3, 3.0 and 4.8°C for daily mean, minimum and maximum air temperatures, respectively) show higher accuracy than local non-glacier surface models (mean root-mean-square differences: 4.2, 4.7 and 5.7°C). In addition, the RGS models based on MODIS night-time LST perform better to estimate daily mean, minimum and maximum air temperatures than using temperature lapse rate derived from local stations.

Knowledge of debris-free and debris-covered glaciers is important for understanding the varying response of glaciers to climate change. Measurements at the debris-free Parlung No. 4 Glacier and the debris-covered 24 K Glacier in the southeastern Tibetan Plateau were carried out to compare the meteorology and surface energy fluxes and to understand the factors controlling the melting process. The meteorological comparisons displayed temporally synchronous fluctuations in air temperature, relative humidity and incoming longwave radiation (Lin), but notable differences in precipitation, incoming shortwave radiation (Sin) and wind speed. Under the prevailing regional precipitation and debris conditions, more Lin (42 W m−2) was supplied from warmer and more humid air and more Sin (58 W m−2) was absorbed at the 24 K Glacier. The relatively high energy supply led mainly to an increased energy output via turbulent heat fluxes and outgoing longwave radiation, rather than glacier melting beneath the thick debris. The sensitivity experiment showed that melt rates were sensitive to energy supply at debris thicknesses <~10 cm. In contrast, energy supply to the Parlung No. 4 Glacier mainly resulted in snow/ice melting, the magnitude of which was significantly influenced by energy supplied by Sin and the sensible heat flux.

The spatio-temporal distribution of air temperature over mountain glaciers can demonstrate complex patterns, yet it is often represented simplistically using linear vertical temperature gradients (VTGs) extrapolated from off-glacier locations. We analyse a network of centreline and lateral air temperature observations at Tsanteleina Glacier, Italy, during summer 2015. On average, VTGs are steep (<−0.0065 °C m−1), but they are shallow under warm ambient conditions when the correlation between air temperature and elevation becomes weaker. Published along-flowline temperature distribution methods explain centreline observations well, including warming on the lower glacier tongue, but cannot estimate lateral temperature variability. Application of temperature distribution methods improves simulation of melt rates (RMSE) in an energy-balance model by up to 36% compared to the environmental lapse rate extrapolated from an off-glacier station. However, results suggest that model parameters are not easily transferable to glaciers with a small fetch without recalibration. Such methods have potential to improve estimates of temperature across a glacier, but their parameter transferability should be further linked to the glacier and atmospheric characteristics. Furthermore, ‘cold spots’, which can be >2°C cooler than expected for their elevation, whose occurrence is not predicted by the temperature distribution models, are identified at one-quarter of the measurement sites.

Previous estimates of melt and surface sublimation on glaciers of the subtropical semiarid Andes (29–34°S) have been obtained at few specific locations, but it is not clear how ablation components vary across the entire extent of a glacier in this dry environment. Here, we simulate the distributed energy and mass balance of Juncal Norte Glacier (33°S) during a 2-month summer period. Forcing fields of near-surface air temperature and wind speed are generated using two methods accounting for the main physical processes that shape their spatial variations. Simulated meteorological variables and ablation agree well with observations on the glacier tongue and reveal complex patterns of energy and mass fluxes. Ablation decreases from 70 mm w.e. d−1 at the low-albedo glacier terminus (~3000 m), where almost 100% of total ablation corresponds to melt, to <5 mm w.e. d−1 at wind-exposed, strong-radiated sites above 5500 m, where surface sublimation represents >75% of total ablation. Our simulations provide the first glacier-scale estimates of ablation components on a glacier in the study region and better reproduce the observed and expected spatial variations of melt and surface sublimation, in comparison with more simple assumptions, such as linear gradients and uniform wind speeds.

Glaciers in northern and southern Norway are subject to different daily and seasonal cycles of incoming solar radiation, which is presumably reflected in the importance of net solar radiation in their surface energy balance. We present a 3 year continuous record from an automatic weather station in the ablation zone of the ice cap Langfjordjøkelen, one of the most northerly glaciers of mainland Norway. Despite its location at 70º N, Langfjordjøkelen was found to have a maritime climate, with an annual mean air temperature of –1.08C, frequent cloud cover and end-of-winter snow depths over 3 m in all three years. The main melt season was May–October, but occasional melt events occurred on warm, cloudy winter days. Net solar and longwave radiation together accounted for 58% of the melt energy, with a positive contribution by net longwave radiation (7%). The sensible and latent heat fluxes supplied the remainder of the melt energy. Cloud optical thickness over Langfjordjøkelen was larger than on two glaciers in southern Norway, especially in the summer months. This resulted in a smaller contribution of net solar radiation to surface melt on Langfjordjøkelen; the effect of the higher latitude on net solar radiation was found to be small.

Albedo variations are presented at two on-glacier sites in the semi-arid Andes, Chile, with >3 years of continuous measurements. Although <2 km apart and at similar elevations, the sites show significantly different albedo cycles. Whereas Toro 1 exhibits a clear seasonal cycle, Guanaco reveals a more complicated pattern, as exposed ice can occur in any month of the year. Daily albedo values are as low as 0.18 for debris-covered Toro 1, while minima are higher on Guanaco (0.38). A method is presented to discern cloud-free from cloudy conditions using measured incoming shortwave radiation only. A cloud climatology is provided showing very low cloudiness values. We see that effective cloud cover relates inversely to cloud occurrence (i.e. either more but thin or fewer but thick clouds). The cloud-free diurnal albedo cycle is pronounced, with lowest values around noon. Clouds increase albedo by 0.04 as a median hourly value, and 0.20 for the 95% quantile. There is a positive relationship between effective cloud cover and resulting albedo rise. Calculations as to whether the diurnal albedo cycle or the effect of clouds on albedo should be considered in energy-balance estimations show that the former is necessary whereas the latter can be neglected in the semi-arid Andes.

The continued use of the positive degree-day (PDD) method to predict ice-sheet melt is generally favoured over surface energy-balance methods partly due to the computational efficiency of the algorithm and the requirement of only one input variable (temperature). In this paper, we revisit some of the assumptions governing the application of the PDD method. Using hourly temperature data from the GC-Net network we test the assumption that monthly PDD total (PDDM) can be represented by a Gaussian distribution with fixed standard deviation of monthly temperature (σM). The results presented here show that the common assumption of fixed σM does not hold, and that σM may be represented more accurately as a quadratic function of average monthly temperature. For Greenland, the mean absolute error in predicting PDDM using our methodology is 3.9°C d, representing a significant improvement on current methods (7.8°C d, when σM = 4.5°C). Over a range of glaciated settings, our method reproduces PDDM, on average, to within 1.5–8.5°C d, compared to 4.4–15.7°C d when σM = 4.5°C. The improvement arises because we capture the systematic reduction in temperature variance that is observed over melting snow and ice, when surface temperatures cannot warm above 0°C.

Multi-annual records of glacier surface meteorology and energy balance are necessary to resolve glacier–climate interactions but remain sparse, especially in the Southern Hemisphere. To address this, we present a record from the ablation zone of Brewster Glacier, New Zealand, between October 2010 and September 2012. The mean air temperature was 1.2°C at 1760 m a.s.l., with only a moderate temperature difference between the warmest and coldest months (∼8°C). Long-term annual precipitation was estimated to exceed 6000 mm a−1, with the majority of precipitation falling within a few degrees of the freezing level. The main melt season was between November and March (83% of annual ablation), but melt events occurred during all months. Annually, net radiation was positive (a source of energy) and supplied 64% of the melt energy, driven primarily by net shortwave radiation. Net longwave radiation was often positive during cloudy conditions in summer, demonstrating the radiative importance of clouds during melt. Turbulent sensible and latent heat fluxes were directed towards the surface in the summer months, accounting for just over a third of the energy for melt (34%). The energy gain associated with rainfall was small except during heavy events in summer.

We assess the runoff and surface mass balance (SMB) of the Greenland ice sheet in the Nuuk region (southwest) using output of two regional climate models (RCMs) evaluated by observations. The region encompasses six glaciers that drain into Godthåbsfjord. RCM data (1960–2012) are resampled to a high spatial resolution to include the narrow (relative to the native grid spacing) glacier trunks in the ice mask. Comparing RCM gridded results with automatic weather station (AWS) point measurements reveals that locally models can underestimate ablation and overestimate accumulation by up to tens of per cent. However, comparison with lake discharge indicates that modelled regional runoff totals are more accurate. Model results show that melt and runoff in the Nuuk region have doubled over the past two decades. Regional SMB attained negative values in recent high-melt years. Taking into account frontal ablation of the marine-terminating glaciers, the region lost 10–20 km3 w.e. a–1 in 2010–12. If 2010 melting prevails during the remainder of this century, a low-end estimate of sea-level rise of 5 mm is expected by 2100 from this relatively small section (2.6%) of the ice sheet alone.

Windpumping has been identified as a process that could potentially enhance sublimation of surface snow at high forcing frequency and spawn air movement deeper in firn at lower frequencies. We performed an experiment to examine the relationship between high-frequency wind and pressure measurements within the top meter of an alpine snowpack and compared experimental results with two theoretical predictions. We find that both theoretical predictions underestimate high-frequency perturbation pressure attenuation with depth in the near-surface snowpack and the discrepancy between theory and measurement increases with perturbation pressure frequency. The impact of this result for near-surface snow is that potential enhanced sublimation will occur over a shallower snow depth than these two theories predict. Correspondingly, interstitial air mixing at depth in firn will be driven by lower frequencies than these two theories predict. While direct measurement of these energy-rich lower frequencies is beyond the scope of this paper, stationary pressure measurements validate the presence of a pressure field that could drive near-surface circulation.

The skill of degree-day glacier melt models is highly dependent on the choice of degree-day factor (DDF), which is often assumed to remain constant in time. Here we explore the validity of this assumption in a changing climate for two locations on Vestari Hagafellsjökull (1979–2012) using a surface energy-balance (SEB) approach that isolates the effect of changes in the prevailing weather on the DDF. At lower elevation, we observe stable DDF during the period of study; however, at higher elevation, DDF is noted to be more variable and a statistically-significant downward trend is observed. This is found to result from an inappropriate threshold air temperature (Tcrit) from which to initiate the positive-degree-day sum, and is removed by setting Tcrit to −1.83°C, rather than the usual value of 0°C used in degree-day melt models. The stationarity of DDF once Tcrit is adjusted contradicts previous research and lends support to the use of constant DDF for projecting future glacier melt. Optimizing Tcrit also improves the skill of melt simulations at our study sites. This research thus highlights the importance of Tcrit for both melt model performance and the evaluation of DDF stationarity in a changing climate.

Air temperature is a key control of processes affecting snow and glaciers in high-elevation catchments, including melt, snowfall and sublimation. It is therefore a key input variable to models of land–surface–atmosphere interaction. Despite this importance, its spatial variability is poorly understood and simple assumptions are made to extrapolate it from point observations to the catchment scale. We use a dataset of 2.75 years of air temperature measurements (from May 2012 to November 2014) at a network of up to 27 locations in the Langtang River, Nepal, catchment to investigate air temperature seasonality and consistency between years. We use observations from high elevations and from the easternmost section of the basin to corroborate previous findings of shallow lapse rates. Seasonal variability is strong, with shallowest lapse rates during the monsoon season. Diurnal variability is also strong and should be taken into account since processes such as melt have a pronounced diurnal variability. Use of seasonal lapse rates seems crucial for glacio-hydrological modelling, but seasonal lapse rates seem stable over the 2–3 years investigated. Lateral variability at transects across valley is high and dominated by aspect, with south-facing sites being warmer than north-facing sites and deviations from the fitted lapse rates of up to several degrees. Local factors (e.g. topographic shading) can reduce or enhance this effect. The interplay of radiation, aspect and elevation should be further investigated with high-elevation transects.