While numerous maps of Greenland ice flow velocity exist, most have gaps in coverage and/or accuracy is limited. We processed a large volume of synthetic aperture radar and Landsat 8 imagery collected between 1995 and 2015 to produce a nearly complete map of ice flow velocity for Greenland at a far greater accuracy than most prior products. We evaluated the accuracy of this map by comparing it with a variety of measured and estimated velocities. For the slow-moving interior of the ice sheet, where estimates are determined from interferometric phase, the errors are ~2 m a−1 or better. For coastal areas, where estimates are determined entirely from speckle- or feature-tracking methods, errors are 2–3 m a−1, which is in good agreement with the estimated formal errors. Especially for the slow-moving majority of the ice sheet, this map provides an important source of data for numerous types of glaciological studies.

Ice cores from mid-latitude mountain glaciers provide detailed information on past climate conditions and regional environmental changes, which is essential for placing current climate change into a longer term perspective. In this context, it is important to define guidelines and create dedicated maps to identify suitable areas for future ice-core drillings. In this study, the suitability for ice-core drilling (SICD) of a mountain glacier is defined as the possibility of extracting an ice core with preserved stratigraphy suitable for reconstructing past climate. Morphometric and climatic variables related to SICD are selected through literature review and characterization of previously drilled sites. A quantitative Weight of Evidence method is proposed to combine selected variables (i.e. slope, local relief, temperature and direct solar radiation) to map the potential drilling sites in mid-latitude mountain glaciers. The method was first developed in the European Alps and then applied to the Asian High Mountains. Model performances and limitations are discussed and first indications of new potential drilling sites in the Asian High Mountains are provided. Results presented here can facilitate the selection of future drilling sites especially on unexplored Asian mountain glaciers towards the understanding of climate and environmental changes.

The mechanics of snow friction are central to competitive skiing, safe winter driving and efficient polar sleds. For nearly 80 years, prevailing theory has postulated that self-lubrication accounts for low kinetic friction on snow: dry-contact sliding warms snow grains to the melting point, and further sliding produces meltwater layers that lubricate the interface. We sought to verify that self-lubrication occurs at the grain scale and to quantify the evolution of real contact area to aid modeling. We used high-resolution (15 µm) infrared thermography to observe the warming of stationary snow under a rotating polyethylene slider. Surprisingly, we did not observe melting at contacting snow grains despite low friction values. In some cases, slider shear failed inter-granular bonds and produced widespread snow movement with no persistent contacts to melt (μ < 0.03). When the snow grains did not move and persistent contacts evolved, the slider abraded rather than melted the grains at low resistance (μ < 0.05). Optical microscopy revealed that the abraded particles deposited in air pockets between grains and thereby carried heat away from the interface, a process not included in current models. Overall, our results challenge whether self-lubrication is indeed the dominant mechanism underlying low snow kinetic friction.

The morphology of englacial drainage networks and their temporal evolution are poorly characterised, particularly within cold ice masses. At present, direct observations of englacial channels are restricted in both spatial and temporal resolution. Through novel use of a terrestrial laser scanning (TLS) system, the interior geometry of an englacial channel in Austre Brøggerbreen, Svalbard, was reconstructed and mapped. Twenty-eight laser scan surveys were conducted in March 2016, capturing the glacier surface around a moulin entrance and the uppermost 122 m reach of the adjoining conduit. The resulting point clouds provide detailed 3-D visualisation of the channel with point accuracy of 6.54 mm, despite low (<60%) overall laser returns as a result of the physical and optical properties of the clean ice, snow, hoar frost and sediment surfaces forming the conduit interior. These point clouds are used to map the conduit morphology, enabling extraction of millimetre-to-centimetre scale geometric measurements. The conduit meanders at a depth of 48 m, with a sinuosity of 2.7, exhibiting teardrop shaped cross-section morphology. This improvement upon traditional surveying techniques demonstrates the potential of TLS as an investigative tool to elucidate the nature of glacier hydrological networks, through reconstruction of channel geometry and wall composition.

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.

We present a comprehensive review of the status and changes in glacier length (since the 1850s), area and mass (since the 1960s) along the Himalayan-Karakoram (HK) region and their climate-change context. A quantitative reliability classification of the field-based mass-balance series is developed. Glaciological mass balances agree better with remotely sensed balances when we make an objective, systematic exclusion of likely flawed mass-balance series. The Himalayan mean glaciological mass budget was similar to the global average until 2000, and likely less negative after 2000. Mass wastage in the Himalaya resulted in increasing debris cover, the growth of glacial lakes and possibly decreasing ice velocities. Geodetic measurements indicate nearly balanced mass budgets for Karakoram glaciers since the 1970s, consistent with the unchanged extent of supraglacial debris-cover. Himalayan glaciers seem to be sensitive to precipitation partly through the albedo feedback on the short-wave radiation balance. Melt contributions from HK glaciers should increase until 2050 and then decrease, though a wide range of present-day area and volume estimates propagates large uncertainties in the future runoff. This review reflects an increasing understanding of HK glaciers and highlights the remaining challenges.

Glaciers spanning large altitudinal ranges often experience different climatic regimes with elevation, creating challenges in acquiring mass-balance and climate observations that represent the entire glacier. We use mixed methods to reconstruct the 1991–2014 mass balance of the Kahiltna Glacier in Alaska, a large (503 km2) glacier with one of the greatest elevation ranges globally (264–6108 m a.s.l.). We calibrate an enhanced temperature index model to glacier-wide mass balances from repeat laser altimetry and point observations, finding a mean net mass-balance rate of −0.74 mw.e. a−1( ± σ = 0.04, std dev. of the best-performing model simulations). Results are validated against mass changes from NASA's Gravity Recovery and Climate Experiment (GRACE) satellites, a novel approach at the individual glacier scale. Correlation is strong between the detrended model- and GRACE-derived mass change time series (R2 = 0.58 and p ≪ 0.001), and between summer (R2 = 0.69 and p = 0.003) and annual (R2 = 0.63 and p = 0.006) balances, lending greater confidence to our modeling results. We find poor correlation, however, between modeled glacier-wide balances and recent single-stake monitoring. Finally, we make recommendations for monitoring glaciers with extreme altitudinal ranges, including characterizing precipitation via snow radar profiling.

We analyzed a 2-year time series of meteorological data (January 2011–December 2012) from three automatic weather stations on Laohugou glacier No. 12, western Qilian Mountains, China. Air temperature, humidity and incoming radiation were significantly correlated between the three sites, while wind speed and direction were not. In this work, we focus on the effects of clouds on other meteorological parameters and on glacier melt. On an average, ~18% of top-of-atmosphere shortwave radiation was attenuated by the clear-sky atmosphere, and clouds attenuated a further 12%. Most of the time the monthly average increases in net longwave radiation caused by clouds were larger than decreases in net shortwave radiation but there was a tendency to lose energy during the daytime when melting was most intense. Air temperature and wind speed related to turbulent heat flux were found to suppress glacier melt during cloudy periods, while increased water vapor pressure during cloudy days could enhance glacier melt by reducing energy loss by latent heat. From these results, we have increased the physical understanding of the significance of cloud effects on continental glaciers.

The recession of mountain glaciers worldwide is increasing global sea level and, in many regions, human activities will have to adapt to changes in surface hydrology. Thus, it is important to provide up-to-date analyses of glacier change and the factors modulating their response to climate warming. Here we report changes in the extent of >120 glaciers on the Lyngen Peninsula, northern Norway, where glacier runoff is utilised for hydropower and where glacial lake outburst floods have occurred. Glaciers covered at least 114 km2 in 1953 and we compare this inventory with those from 1988, 2001 and a new one from 2014, and previously-dated Little Ice Age (LIA) limits. Results show a steady reduction in area (~0.3% a−1) between their LIA maximum (~1915) and 1988, consistent with increasing summer air temperatures, but recession paused between 1988 and 2001, coinciding with increased winter precipitation. Air temperatures increased 0.5°C per decade from the 1990s and the rate of recession accelerated to ~1% a−1 between 2001 and 2014 when glacier area totalled ~95.7 km2. Small glaciers (<0.05 km2) with low maximum elevations (<1400 m) experienced the largest percentage losses and, if warming continues, several glaciers may disappear within the next two decades.

We present the glacier-wide summer surface mass balances determined by a detailed hydrological balance (sSMBhydro) and the quantification of the uncertainties of the calculations on the Argentière and Mer de Glace-Leschaux drainage basins, located in the upper Arve watershed (French Alps), over the period 1996–2004. The spatial distribution of precipitation within the study area was adjusted using in situ winter mass-balance measurements. The sSMBhydro performance was assessed via a comparison with the summer surface mass balances based on in situ glaciological observations (sSMBglacio). Our results show that the sSMBhydro has an uncertainty of ± 0.67 m w.e. a−1 at Argentière and ± 0.66 m w.e. a−1 at Mer de Glace-Leschaux. Estimates of the Argentière sSMBhydro values are in good agreement with the sSMBglacio values. These time series show almost the same interannual variability. From the marked difference between the sSMBhydro and sSMBglacio values for the Mer de Glace-Leschaux glacier, we suspect a significant role of groundwater fluxes in the hydrological balance. This study underlines the importance of taking into account the groundwater transfers to represent and predict the hydro-glaciological behaviour of a catchment.

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.

In a series of papers, Bartelt and co-workers developed novel snow-avalanche models in which random kinetic energy (RKE) RK (a.k.a. granular temperature) is a key concept. The earliest models were for a single, constant density layer, using a Voellmy model but with RK-dependent friction parameters. This was then extended to variable density, and finally a suspension layer (powder-snow cloud) was added. The physical basis and mathematical formulation of these models are critically reviewed here, with the following main findings: (i) Key assumptions in the original RKE model differ substantially from established results on dense granular flows; in particular, the effective friction coefficient decreases to zero with velocity in the RKE model. (ii) In the variable-density model, non-canonical interpretation of the energy balance leads to a third-order evolution equation for the flow depth or density, whereas the stated assumptions imply a first-order equation. (iii) The model for the suspension layer neglects gravity and disregards well-established theoretical and experimental results on particulate gravity currents. Some options for improving these aspects are discussed.

The critique by Issler and other (2017) of our avalanche dynamics model maintains that we disregard well-established results of particulate gravity flows. Here we show that the arguments of Issler and others (1) violate Newton's laws of motion in the avalanche core and (2) ignore size-dependent drag forces on snow particles (Stokes law) that lead to the formation of avalanche suspension layers. We explain why we cannot amend our model equations to accommodate their suggestions. The goal of our approach is to describe highly non-stationary processes in the avalanche core that lead to a wide range of avalanche flow types and therefore different flow behavior. This is important for practical applications.