Two small glaciers on James Ross Island, the north-eastern Antarctic Peninsula, experienced surface mass gain between 2009 and 2015 as revealed by field measurements. A positive cumulative surface mass balance of 0.57 ± 0.67 and 0.11 ± 0.37 m w.e. was observed during the 2009–2015 period on Whisky Glacier and Davies Dome, respectively. The results indicate a change from surface mass loss that prevailed in the region during the first decade of the 21st century to predominantly positive surface mass balance after 2009/10. The spatial pattern of annual surface mass-balance distribution implies snow redistribution by wind on both glaciers. The mean equilibrium line altitudes for Whisky Glacier (311 ± 16 m a.s.l.) and Davies Dome (393 ± 18 m a.s.l.) are in accordance with the regional data indicating 200–300 m higher equilibrium line on James Ross and Vega Islands compared with the South Shetland Islands. The mean accumulation-area ratio of 0.68 ± 0.09 and 0.44 ± 0.09 determined for Whisky Glacier and Davies Dome, respectively, is similar to the value reported for Vega Island and within the range of typical values for high-latitude glaciers.

Calving plays a key role in the recent rapid retreat of glaciers around the world. However, many processes related to calving are poorly understood since direct observations are scarce and challenging to obtain. When calving occurs at a glacier front, surface-water waves arise over the ocean or a lake in front of glaciers. To study calving processes from these surface waves, we performed field observations at Glaciar Perito Moreno, Patagonia. We synchronized time-lapse photography and surface waves record to confirm that glacier calving produces distinct waves compared with local noise. A total of 1074 calving events were observed over the course of 39 d. During austral summer, calving occurred twice more frequently than in spring. The cumulative distribution of calving-interevent time interval followed exponential model, implying random occurrence of events in time. We further investigated wave properties and found that source-to-sensor distance can be estimated from wave dispersion within ~20% error. We also found that waves produced by different calving types showed similar spectra in the same frequency range between 0.05–0.2 Hz, and that the amplitude of surface waves increased with the size of calving. This study demonstrates the potential of surface-wave monitoring for understanding calving processes.

Light-absorbing impurities (LAI) can darken snow and ice surfaces, reduce snow/ice albedo and accelerate melt. Efforts to allocate the relative contribution of different LAI to snow/ice albedo reductions have been limited by uncertainties in the optical properties of LAI. We developed a new method to measure LAI spectral reflectance at the submicron scale by modifying a Hyperspectral Imaging Microscope Spectrometer (HIMS). We present the instrument's internal calibration, and the overall small influence of a particle's orientation on its measured reflectance spectrum. We validated this new method through the comparison with a field spectroradiometer by measuring different standard materials. Measurements with HIMS at the submicron scale and the bulk measurements of the same standard materials with the field spectroradiometer are in good agreement with an average deviation between the spectra of 3.2% for the 400–1000 nm wavelength range. The new method was used (1) to identify BC (black carbon), mineral dust including hematite and the humic substances present in an environmental sample from Plaine Morte glacier and (2) to collect the individual reflectance spectra of each of these types of impurity. The results indicate that this method is applicable to heterogeneous samples such as the LAI found in snow and ice.

Increasing ship traffic and human activity in the Arctic has led to a growing demand for accurate Arctic weather forecast. High-quality forecasts obtained by models are dependent on accurate initial states achieved by assimilation of observations. In this study, a multi-variate nudging (MVN) method for assimilation of sea-ice variables is introduced. The MVN assimilation method includes procedures for multivariate update of sea-ice volume and concentration, and for extrapolation of observational information spatially. The MVN assimilation scheme is compared with the Ensemble Kalman Filter (EnKF) using the Los Alamos Sea Ice Model. Two multi-variate experiments are conducted: in the first experiment, sea-ice thickness from the European Space Agency's Soil Moisture and Ocean Salinity mission is assimilated, and in the second experiment, sea-ice concentration from the ocean and Sea Ice Satellite Application Facility is assimilated. The multivariate effects are cross-validated by comparing the model with non-assimilated observations. It is found that the simple and computationally cheap MVN method shows comparable skills to the more complicated and expensive EnKF method for multivariate update. In addition, we show that when few observations are available, the MVN method is a significant model improvement compared to the version based on one-dimensional sea-ice concentration assimilation.

Twin glaciers collapsed in 2016 near Aru Co, western Tibet and caused extreme loss to human beings. In this study, we attempted to track the dynamics of glaciers in the region, for example the glacier area and mass changes in Aru Co for the period 1971–2016, which were determined using topographic maps and Landsat images and ASTER-derived DEMs (2011–16), the Shuttle Radar Terrain Mission DEM (2000) and topographic maps (1971). Our results showed that the glacier area of Aru Co decreased by −0.4 ± 4.1% during 1971–2016. The geodetic mass-balance results showed that the glaciers in Aru Co lost mass at a rate of −0.15 ± 0.30 m w.e. a−1 during 1971–99, while they gained mass at a rate of 0.33 ± 0.61 m w.e. a−1 for the period 1999–2016. The twin glaciers experienced a larger negative mass budget than the others in the region before 1999. This process produced large amounts of meltwater, followed by a sustained increase in the meltwater on the pressure melting point, possibly in response to a period of positive mass balance (1999–2016) and then, transferred to the glacier bed until the glaciers collapsed.

Thermal and hydrological dynamics and their impacts on the stability of a moraine dam were analyzed and simulated for the Longbasaba Lake in the Himalaya, based on soil temperature, moisture and heat flux data observed at different depths in the dam from 2012 to 2016. Annual average heat income is greater than heat expenditure on the dam surface. The mean annual temperature at observed the depths of 0–150 cm is >0°C, although the average annual air temperature was −3.6°C over the dam, indicating a relatively larger temperature difference between moraine dam and air. The volumetric soil moisture content is relatively low with an annual average of 5%, peaking after the snow cover melting and active layer thawing. Simulation results indicate that the average yearly maximum thawing depth has been ~0.3 m deeper than the average yearly maximum freezing depth during the observation period. In the past 55 years, the yearly maximum thawing depth has increased, while yearly maximum freezing depth has decreased, implying that the permafrost in the dam has been deteriorating. The annual surplus heat and increasing permafrost thawing depth will result in further deterioration of permafrost and melting of buried ice in the dam, thereby decreasing its stability.

Determining the position and stability of the grounding line of a marine ice sheet is a major challenge for ice-sheet models. Here, we investigate the role of lateral shear and ice-shelf buttressing in grounding line dynamics by extending an existing boundary layer theory to laterally confined marine ice sheets. We derive an analytic expression for the ice flux at the grounding line of confined marine ice sheets that depends on both local bed properties and non-local ice-shelf properties. Application of these results to a laterally confined version of the MISMIP 1a experiment shows that the boundary condition at the ice-shelf front (i.e. the calving law) is a major control on the location and stability of the grounding line in the presence of buttressing, allowing for both stable and unstable grounding line positions on downwards sloping beds. These results corroborate the findings of existing numerical studies that the stability of confined marine ice sheets is influenced by ice-shelf properties, in contrast to unconfined configurations where grounding line stability is solely determined by the local slope of the bed. Consequently, the marine ice-sheet instability hypothesis may not apply to buttressed marine ice sheets.

Microfossil records from ice archives allow vegetation, fire and land-use activity reconstructions on broad spatial scales. Samples typically contain low microfossil concentrations. Therefore, large ice volumes are often needed for palynology. Hence, it is crucial to extract maximum microfossil numbers through appropriate physical-chemical treatments. We compare six methods covering the main water reduction procedures: evaporation, filtration and centrifugation with snow samples. Adding a known number of Lycopodium marker spores prior to sample treatment and a second marker (Eucalyptus) after laboratory processing allows a quantitative microfossil loss assessment during pollen extraction. We applied the best-performing method (average loss of 22%) to high-alpine firn cores from Colle Gnifetti glacier for validation with a natural archive containing extremely low microfossil concentrations. We conclude that samples processed with different microfossil extraction protocols may give different results for pollen concentrations, percentages and ratios between different pollen types, especially if vesiculate conifer pollen is an important pollen assemblage component. We recommend a new evaporation-based method which delivers the smallest and least variable losses among the tested approaches. Since microfossil losses are inevitable during laboratory procedure, adding markers prior to sample processing is mandatory to achieve reliable microfossil concentration and influx estimates.

The permeability of sea ice can strongly affect the dissipation of wave energy into the ice pack. Sea-ice permeability is known to be impacted by the brine volume fraction and the blockage of flow pathways by the freezing of infiltrating lower salinity water. Here we investigate another process impacting sea-ice permeability, namely, inelastic deformation. We report the results of a first-of-its-kind field-scale deformation experiment to investigate the impact of compressive loading on sea-ice permeability. We observed that deformation decreased permeability by four orders of magnitude or more in some locations, while elsewhere permeability was unaffected or possibly increased. We show that the observed changes in permeability are consistent with expected changes in stress state and, as a result, in the mechanisms of deformation.

Repeat airborne laser altimetry measurements show widespread thinning (surface lowering) of glaciers in Canada's Queen Elizabeth Islands since 1995. Thinning rates averaged for 50 m elevation bins, were more than three times higher during the period 2005/06 to 2012/14 pentad than during the previous two pentads. Strongly negative thickness change (dh/dt) anomalies from 2005/06 to 2012/14, relative to the 1995–2012/14 mean, suggest that most of the measured thinning occurred during the most recent 5–6 year period when mean summer land surface temperatures (LSTs) were anomalously high and the mean summer black-sky shortwave broadband albedos (BSA) were anomalously low, relative to the 2000/01–15/16 period, and upper-air (700 hPa) and near surface (2 m) air temperatures were between 0.8°C and 1.5°C higher than 1995–2012 mean. Comparisons of dh/dt with mean summer LST and BSA measurements from the Moderate Resolution Imaging Spectroradiometer and with surface longitudinal strain rates computed from surface velocity fields derived from RADARSAT 1/2 and Landat-7 ETM + data suggest that surface elevation changes were driven mainly by changes in climate. An exception to this occurs along many fast-flowing outlet glaciers where ice dynamics appear also to have played an important role in surface elevation changes.

Unlike in most other regions, Karakoram glaciers are either stable or advancing, a phenomenon known as the Karakoram anomaly. Despite studies of glacier surges and the derivation of surface velocity maps, the spatiotemporal variability of glacier dynamics still remains poorly understood, particularly in the Eastern Karakoram Range. We use Advanced Land Observing Satellite/the Phased Array type L-band Synthetic Aperture Radar (ALOS/PALSAR)-1/2 data from 2007 to 2011 and 2014 to 2015 to examine detailed surface velocity patterns of the Siachen, Baltoro, Kundos, Singkhu and Gasherbrum Glaciers. The first three glaciers show considerable velocity variability (20–350 m a−1), with clear seasonal patterns. Although all glaciers, except for Baltoro, flow slowest in 2015, the velocity structures are individual and vary in space and time. In Gasherbrum Glacier, peak surge-phase velocities are seasonally modulated, with maxima in summers 2006 and 2007, suggesting surface melt plays an important role in maintaining the active phase. Given the relatively close proximity of these glaciers, we assume that surface melt timing and rates are comparable. We therefore argue that the observed spatiotemporal and interannual velocity patterns are determined by local and internal mechanisms, including englacial and subglacial hydrology, thermal processes and tributary configuration of each individual glacier.

The Holocene climatic optimum was a period 8–5 kyr ago when annual mean surface temperatures in Greenland were 2–3°C warmer than present-day values. However, this warming left little imprint on commonly used temperature proxies often used to derive the climate forcing for simulations of the past evolution of the Greenland ice sheet. In this study, we investigate the evolution of the Greenland ice sheet through the Holocene when forced by different proxy-derived temperature histories from ice core records, focusing on the effect of sustained higher surface temperatures during the early Holocene. We find that the ice sheet retreats to a minimum volume of ~0.15–1.2 m sea-level equivalent smaller than present in the early or mid-Holocene when forcing an ice-sheet model with temperature reconstructions that contain a climatic optimum, and that the ice sheet has continued to recover from this minimum up to present day. Reconstructions without a warm climatic optimum in the early Holocene result in smaller ice losses continuing throughout the last 10 kyr. For all the simulated ice-sheet histories, the ice sheet is approaching a steady state at the end of the 20th century.

The climatic condition that prevails at a glacier equilibrium line altitude (ELA) is often parameterized in terms of summer air temperature (T) and annual precipitation (P). This simple parameterization was initially proposed by Hans W:son Ahlmann. The physical background of the relationship between P and T on the equilibrium line, however, has been left unexplained since Ahlmann first questioned the mathematical form of the relationship. This relationship can be explained when the thermal and hydrological processes of the ELA formation are investigated. The present authors studied the energy exchange processes that prevail on the ELA during the melt season. The inclusion of solar radiation brings Ahlmann's hypothesis closer to energy balance, and improves his P/T diagram. By comparing the observed fluxes from the polar through the mid-latitude to the equatorial glaciers, it was found that these glaciers in different climatic regions share important similarities at the ELA. Further, it was found that the classic P/T curve originally proposed by Ahlmann in the early 20th century is a concise expression of the conservation principle of energy and mass at the ELA of glaciers, and takes the form of a polynomial of the fourth order.

The Tibetan Plateau (TP) lakes are sensitive to climate change due to ice-albedo feedback, but almost no study has paid attention to the ice albedo of TP lakes and its potential impacts. Here we present a recent field experiment for observing the lake ice albedo in the TP, and evaluate the applicability of the Moderate Resolution Imaging Spectroradiometer (MODIS) products as well as ice-albedo parameterizations. Most of the observed lake ice albedos on TP are <0.12, and the clear blue ice albedo is only 0.075, much lower than reported in the previous studies. Even that of ice covered with snow patches is only 0.212. MOD10A1 albedo product has the best agreement with observations, followed by those of MYD10A1. MCD43A3 product is consistently higher than the observations. Due to an error of snow flag and inconsistent time windows in MCD43A2 and MCD43A3, at certain times, the albedo of the ice without snow is even higher than that covered with snow. When the solar zenith angle is not considered, there is no significant correlation between the albedo and the ice surface temperature. None of the existing ice-albedo parameterizations can reproduce well the observed relationship of the albedo and surface temperature.