The subglacial observatory beneath the Argentière glacier provides a rare opportunity to study the interactions between glacier sliding velocity and subglacial runoff. The sliding velocity has been monitored in this cavity almost continuously since 1997 and the resulting data indicate a decrease in annual sliding velocities over the last two decades. We found close relationships between annual surface velocity, sliding velocity and ice thickness. These relationships indicate that the ice-flow velocity changes do not depend on subglacial water runoff changes at the annual timescale. The seasonal magnitudes of sliding also show a decrease over the last two decades. At the seasonal timescale, sliding velocity increases before or simultaneously with the large runoff increase in May, indicating a distributed drainage system. Conversely, at the end of the melt season, sliding velocity continues to decrease after the runoff returns to low winter values. The simultaneous increases of runoff and sliding velocity occur mainly before the spring transition. Later, sliding velocity generally appears not to be related to water inputs coming from the surface, except for some large accelerations after midAugust that are always associated with periods of rapidly increasing water inputs to the subglacial drainage system.

An extensive set of shore-parallel fractures are observed at the grounding line of Kamb Ice Stream (KIS) in West Antarctica. Seismicity measured in the grounding zone is, as elsewhere around Antarctica, tidally forced and moreover strand cracks propagate nearly exclusively on the falling tide. Measured surface deflection and a model of fracture propagation are used to conclude that bending on the falling tide favors propagation while bending on the rising tide suppresses propagation. Without the perturbation due to tidal bending, strand cracks would be rare and appear farther downstream than observed. We speculate that the very large number of cracks observed at KIS is due to the stagnant-to-floating transition at that grounding line, which allows cyclic bending of the same ice and relatively large stretching rates.

A 1-D ice cover model was developed to predict and constrain drivers of long-term ice thickness trends in chemically stratified lakes of Taylor Valley, Antarctica. The model is driven by surface radiative heat fluxes and heat fluxes from the underlying water column. The model successfully reproduced 16 a (between 1996 and 2012) of ice thickness changes for the west lobe of Lake Bonney (average ice thickness = 3.53 m) and Lake Fryxell (average ice thickness = 4.22 m). Long-term ice thickness trends require coupling with the thermal structure of the water column. The heat stored within the temperature maximum of lakes exceeding a liquid water column depth of 20 m can either impede or facilitate ice thickness change depending on the predominant climatic trend (cooling or warming). As such, shallow (<20 m deep water columns) perennially ice-covered lakes without deep temperature maxima are more sensitive indicators of climate change. The long-term ice thickness trends are a result of surface energy flux and heat flux from the deep temperature maximum in the water column, the latter of which results from absorbed solar radiation.

A short-term series of quantitative observations of calving activity of Fuerza Aérea Glacier (Greenwich Island, the South Shetland Islands, Antarctica) was conducted in order to test new methods of monitoring calving. The volume of single calving events was quantified by combining terrestrial laser scanning (TLS) surveys with continuous video recording of the ice front. An empirical formula for area/volume scaling of the calved ice block was proposed based on the TLS measured calved ice volume and the calved ice front area obtained by manual delineation on the images acquired with the video camera. This combination of methods proves to be a valuable tool for glacier monitoring, providing both high-temporal resolution and precise quantitative measurements of the calving volume. The size distribution of calving events is best approximated by a power law and within the short period of observations (14 d) calving was found to be an intrinsic process not dependent on environmental forcings. Over the period of 21 January–04 February 2013 the ice flow velocity at the terminus of Fuerza Aérea Glacier was 0.26 ± 0.07 m d−1 and the calving rate was 0.41 ± 0.07 m d−1.

For the Andes Cordillera, where observed mass-balance records are sparse, long-term glacier velocity measurements potentially represent a useful tool for assessing glacier health. Utilising manual and automatic feature-tracking techniques applied to Corona, Landsat and ASTER satellite imagery, this paper presents surface velocity fluctuations for Glaciar Universidad between 1967 and 1969, and 1985 and 2015, the longest such time series available for the Andes Cordillera, outside Patagonia. This time series reveals an increase in the surface velocities of the main glacier trunk between 1967 and 1987 (~90%) followed by a deceleration between 1987 and 2015 (~80%), with ice velocities observed between 2014 and 2015 possibly representing a 48 a low. In response to the surface velocity fluctuations, the glacier front advanced between 1985 and 1992 (cumulative change of 137 ± 14 m), and again to a lesser magnitude during the 1996–98 and 2004–08 periods. Although having exhibited possible surge behaviour during the 1940s, the synchrony of the glacier changes presented for Glaciar Universidad with those reported for nearby glaciers, suggests that this glacier is responding to climatic trends. If the above scenario is true, the results indicate a general pattern of increasingly negative glacier mass-balance conditions since the late 1980s.

Glacier surface mass-balance measurements on Greenland started more than a century ago, but no compilation exists of the observations from the ablation area of the ice sheet and local glaciers. Such data could be used in the evaluation of modelled surface mass balance, or to document changes in glacier melt independently from model output. Here, we present a comprehensive database of Greenland glacier surface mass-balance observations from the ablation area of the ice sheet and local glaciers. The database spans the 123 a from 1892 to 2015, contains a total of ~3000 measurements from 46 sites, and is openly accessible through the PROMICE web portal (http://www.promice.dk). For each measurement we provide X, Y and Z coordinates, starting and ending dates as well as quality flags. We give sources for each entry and for all metadata. Two thirds of the data were collected from grey literature and unpublished archive documents. Roughly 60% of the measurements were performed by the Geological Survey of Denmark and Greenland (GEUS, previously GGU). The data cover all regions of Greenland except for the southernmost part of the east coast, but also emphasize the importance of long-term time series of which there are only two exceeding 20 a. We use the data to analyse uncertainties in point measurements of surface mass balance, as well as to estimate surface mass-balance profiles for most regions of Greenland.

We present results of a regional comparative study of surface mass changes from 2004 to 2008 based on Gravity Recovery and Climate Experiment (GRACE), The Ice, Cloud and Land Elevation Satellite (ICESat) and CHINARE observations over the Lambert Glacier/Amery Ice Shelf system (LAS). Estimation of the ICESat mass change rates benefitted from the density measurements along the CHINARE traverse and a spatial density adjustment method for reducing the effect of spatial density variations. In the high-elevation inland region, a positive trend was estimated from both ICESat and GRACE data, which is in line with the CHINARE accumulation measurements. In the coastal region, there were areas with high level accumulations in both ICESat and GRACE trend maps. In many high flow-speed glacier areas, negative mass change rates may be caused by dynamic ice flow discharges that have surpassed the snow accumulation. Overall, the mass change rate estimate in the LAS of 2004–2008 from the GRACE, ICESat and CHINARE data is 5.41 ± 4.59 Gt a−1, indicating a balanced to slightly positive mass trend. Along with other published results, this suggests that a longer-term positive mass trend in the LAS may have slowed in recent years.

The dominant mass-loss process on the Antarctic Peninsula has been ice-shelf collapse, including the Larsen A Ice Shelf in early 1995. Following this collapse, there was rapid speed up and thinning of its tributary glaciers. We model the impact of this ice-shelf collapse on upstream tributaries, and compare with observations using new datasets of surface velocity and ice thickness. Using a two-horizontal-dimension shallow shelf approximation model, we are able to replicate the observed large increase in surface velocity that occurred within Drygalski Glacier, Antarctic Peninsula. The model results show an instantaneous twofold increase in flux across the grounding line, caused solely from the reduction in backstress through ice shelf removal. This demonstrates the importance of ice-shelf buttressing for flow upstream of the grounding line and highlights the need to explicitly include lateral stresses when modelling real-world settings. We hypothesise that further increases in velocity and flux observed since the ice-shelf collapse result from transient mass redistribution effects. Reproducing these effects poses the next, more stringent test of glacier and ice-sheet modelling studies.

We explore a new method to retrieve seasonal glacier mass balances (MBs) from low-resolution optical remote sensing. We derive annual winter and summer snow maps of the Alps during 1998–2014 using SPOT/VEGETATION 1 km resolution imagery. We combine these seasonal snow maps with a DEM to calculate a ‘mean regional’ altitude of snow (Z) in a region surrounding a glacier. Then, we compare the interannual variation of Z with the observed winter/summer glacier MB for 55 Alpine glaciers over 1998–2008, our calibration period. We find strong linear relationships in winter (mean R2 = 0.84) and small errors for the reconstructed winter MB (mean RMSE = 158 mm (w.e.) a−1). This is lower than errors generally assumed for the glaciological MB measurements (200–400 mm w.e. a−1). Results for summer MB are also satisfying (mean R2 and RMSE, respectively, 0.74 and 314 mm w.e. a−1). Comparison with observed seasonal MB available over 2009–2014 (our evaluation period) for 19 glaciers in winter and 13 in summer shows good agreement in winter (RMSE = 405 mm w.e. a−1) and slightly larger errors in summer (RMSE = 561 mm w.e. a−1). These results indicate that our approach might be valuable for remotely determining the seasonal MB of glaciers over large regions.

Empirical data on the volume/area scaling of glaciers and ice caps are compared with theoretical models of longitudinal glacier profiles. Ice caps are described well by theoretical solutions of the differential equations modeling longitudinal profiles, but glaciers for which the surface slope is approximated by the ice thickness divided by the length are not. A power-law profile recently derived analytically, reproduces exactly the thickness-length scaling, first predicted by scaling analysis for ice caps.

Supraglacial deposits are known for their influence on glacier ablation. The magnitude of this influence depends on the thickness and the type of the deposited material. The effects of thin layers of atmospheric black carbon and of thick moraine debris have been intensively studied. Studies related to regional-scale deposits of volcanic tephra with thicknesses varying between millimetres and metres and thus over several orders of magnitude are scarce. We present results of a field experiment in which we investigated the influence of supraglacial deposits of tephra from Grímsvötn volcano on bare-ice ablation at Svínafelsjökull, Iceland. We observed that the effective thickness at which ablation is maximized ranges from 1.0 to 2.0 mm. At ~10 mm a critical thickness is reached where sub-tephra ablation equals bare-ice ablation. We calibrated two empirical ablation models and a semi-physics-based ablation model that all account for varying tephra-layer thicknesses. A comparison of the three models indicates that for tephra deposits in the lower-millimetre scale the temperature/radiation-index model performs best, but that a semi-physics-based approach could be expected to yield superior results for tephra deposits of the order of decimetres.

We investigate an internal surge of Karayaylak Glacier, which was reported by the media in May 2015. To differentiate the May 2015 glacier surge from other glacier advances, we surveyed changes in velocity, crevasses and glacier area using Landsat 8 OLI L1T, ZY-1-02C and Gaofen-1 images from October 2014 to July 2015. The velocity, measured by automatic feature extraction and tracking during the active phase, was 10–100 times the velocity during the quiescent phase, with a maximum of (20.2 ± 0.9) m d−1 (mean ± standard error) from 8 to 15 May 2015 in the west branch of the glacier. The surge initiation and termination took place from 13 April to 16 June 2015. Ice in the west branch (length, 7 km; area, 6.8 km2) of Karayaylak Glacier accelerated down to the east branch, leading to the development of crevasses and ice covering an additional 0.1 km2 of summer pasture on the northwestern side. However, we detected no advance of the glacier's terminus.

Recent studies have suggested a physical link between Ca++ ions and an increase in the ductility or ‘softening’ of polycrystalline ice. In order to investigate the potential effects of Ca++ on deformation, we created sets of both undoped and CaSO4-doped specimens of polycrystalline ice for testing in uniaxial tension or compression. Deformation tests in tension were carried out under a constant load at an initial stress of 0.75 MPa and a temperature of −6°C. Compression tests were carried out at −10 and −20°C at constant strain rates of 1×10−4 s−1, 1 × 10−5 s−1 and 1 × 10−6 s−1 and taken to 5% strain. Our results show that CaSO4 increases the strength of polycrystalline ice at higher strain rates and lower temperatures, an effect that decreases with decreasing strain rate and higher temperatures. A microstructural analysis of the post-test compression specimens reveals mean grain diameters much larger in the CaSO4-doped specimens tested at the lowest applied strain rate of 1 × 10−6 s−1. Precipitates were found to have formed along grain boundaries in some doped specimens and evidence of intergranular fracture was observed in all specimens tested at 1 × 10−4 and 1 × 10−5 s−1. In tension-tested specimens, there was no difference in the mean grain diameter between doped and undoped specimens at 25% strain.

Canada's Queen Elizabeth Islands (QEI) contain ~14% of the world's glacier and ice-cap area. Sparse in-situ measurements indicate that interannual variability in glacier surface mass balance in this region is driven primarily by variations in summer melt, and that the annual surface mass balance of four index glaciers has become increasingly negative since 2007. Here, we use a 16-a record of satellite-derived mean summer (June–August) land surface temperatures (LST) from NASA's Moderate Resolution Imaging Spectroradiometer to investigate large-scale spatial and temporal variability in the duration and intensity of summer melt across glaciated surfaces in the QEI from 2000 to 2015. During this period, QEI mean summer glacier surface temperatures increased at an average rate of 0.06 ± 0.04°C a−1, for a total of nearly 1°C. Most of this increase occurred between 2005 and 2012, when mean summer near-surface (2 m) and upper-air (700 hPa) temperatures were 1.0–1.2°C higher than the 1948–2015 mean. There is a strong correlation between the glacier LST and 700 hPa air temperature records (r> 0.8). The period 2005–12, when mean summer LSTs were anomalously high, was likely the warmest period in the region since at least 1948.

The McMurdo Dry Valleys are a cold, dry polar desert and the alpine glaciers therein exhibit small annual and seasonal mass balances, often <±0.06 m w.e. Typically, winter is the accumulation season, but significant snow storms can occur any time of year occasionally making summer the accumulation season. The yearly equilibrium line altitude is poorly correlated with mass balance because the elevation gradient of mass balance on each glacier can change dramatically from year to year. Most likely, winds redistribute the light snowfall disrupting the normal gradient of increasing mass balance with elevation. Reconstructed cumulative mass balance shows that the glaciers have lost <2 m w.e. over the past half century and area changes show minimal retreat. In most cases these changes are less than the uncertainty and the glaciers are considered in equilibrium. Since 2000, however, the glaciers have lost mass despite relatively stable summer air temperatures suggesting a different mechanism in play. Whether this trend is a harbinger of future changes or a temporary excursion is unclear.