Three debris-free glaciers with strongly differing annual glaciological glacier-wide mass balances (MBs) are monitored in the Everest region (central Himalaya, Nepal). The mass budget of Mera Glacier (5.1 km2 in 2012), located in the southern part of this region, was balanced during 2007–15, whereas Pokalde (0.1 km2 in 2011) and West Changri Nup glaciers (0.9 km2 in 2013), ~30 km further north, have been losing mass rapidly with annual glacier-wide MBs of −0.69 ± 0.28 m w.e. a−1 (2009–15) and −1.24 ± 0.27 m w.e. a−1 (2010–15), respectively. An analysis of high-elevation meteorological variables reveals that these glaciers are sensitive to precipitation, and to occasional severe cyclonic storms originating from the Bay of Bengal. We observe a negative horizontal gradient of annual precipitation in south-to-north direction across the range (≤−21 mm km−1, i.e. −2% km−1). This contrasted mass-balance pattern over rather short distances is related (i) to the low maximum elevation of Pokalde and West Changri Nup glaciers, resulting in years where their accumulation area ratio is reduced to zero and (ii) to a steeper vertical gradient of MB for glaciers located in the inner arid part of the range.

We evaluate the variability in basal friction for Kronebreen, Svalbard, a fast-flowing tidewater glacier. We invert 3 years (2013–15) of surface velocities at high temporal resolution (generally 11 days), to estimate the changing basal properties of the glacier. Our results suggest that sliding behaviour of Kronebreen within a year is primarily influenced by changes in water input patterns during the meltwater season and basal friction is highly variable from a year to another. At present, models usually employ parameterisations to encompass the complex physics of glacier sliding by mathematically simulate their net effect. For such ice masses with strong seasonal variations of surface melt, the spatio-temporal patterns of basal friction imply that it is neither possible nor appropriate to use a parameterisation for bed friction that is fixed in space and/or time, at least in a timescale of a few years. Basal sliding may not only be governed by local processes such as basal topography or summer melt, but also be mediated by factors that vary over a larger distance and over a longer time period such as subglacial hydrology organisation, ice-thickness changes or calving front geometry.

Snow metamorphism and settlement change the microstructure of a snowpack simultaneously. Past experiments investigated snow deformation under isothermal conditions. In nature, temperature gradient metamorphism and settlement often occur together. We investigated snow settlement in the first days after the onset of temperature-gradient metamorphism in laboratory experiments by means of in-situ time-lapse micro-computed tomography. We imposed temperature gradients of up to 95 K m−1 on samples of rounded snow with a density of ~230 kg m−3 and induced settlement by applying 1.7 kPa stress with a passive load on the samples simultaneously. We found that snow settled about half as fast when a temperature gradient was present, compared with isothermal conditions. The change in specific surface area after 4 days caused by temperature-gradient metamorphism was only a few percent. The viscosity evolution correlated with the amount of the temperature gradient. Finite element simulations of the snow samples revealed that stress-bearing chains had developed in the snow structure, causing the large increase in viscosity. We could show that a small change in microstructure caused a large change in the mechanical properties. This explains the difficulty of predicting snow mechanical properties in applications such as firn compaction or snow avalanche formation.

Understanding the drivers of recent change at Greenlandic tidewater glaciers is of great importance if we are to predict how these glaciers will respond to climatic warming. A poorly constrained component of tidewater glacier processes is the near-terminus subglacial hydrology. Here we present a novel method for constraining near-terminus subglacial hydrology with application to marine-terminating Kangiata Nunata Sermia in South-west Greenland. By simulating proglacial plume dynamics using buoyant plume theory and a general circulation model, we assess the critical subglacial discharge, if delivered through a single compact channel, required to generate a plume that reaches the fjord surface. We then compare catchment runoff to a time series of plume visibility acquired from a time-lapse camera. We identify extended periods throughout the 2009 melt season where catchment runoff significantly exceeds the discharge required for a plume to reach the fjord surface, yet we observe no plume. We attribute these observations to spatial spreading of runoff across the grounding line. Persistent distributed drainage near the terminus would lead to more spatially homogeneous submarine melting and may promote more rapid basal sliding during warmer summers, potentially providing a mechanism independent of ocean forcing for increases in atmospheric temperature to drive tidewater glacier acceleration.

Global climate change is causing widespread glacier retreat, with many small glaciers disappearing from the world's mountain ranges. We obtained the annual mass balance of a small glacier (the Ningchan No.1 Glacier) located on the northeastern Tibetan Plateau, from the years 2010 to 2015 using glaciological and geodetic methods. We also measured the glacier's thickness in 2014 using ground-penetrating radar. Employing topographical maps and ZY-3 images, we obtained Digital Elevation Models for 1972 and 2014. Our results showed that the mean annual mass balance from 2010 to 2015 was ~−0.9 ± 0.5 m w.e. The mean equilibrium line altitude was ~4680 m in the period 2010–15, which exceeds the maximum elevation of the glacier. The glacier has lost area and mass across its elevation range. The mean ice thickness was 24.0 ± 2.5 m in 2014. From 1972 to 2014, the glacier's area shrank from 0.77 ± 0.05 to 0.39 ± 0.04 km2, and the ice volume decreased by (14.96 ± 0.97) × 106 m3, equivalent to (12.72 ± 0.82) × 106 t w.e. over the same period.

Supraglacial ponds play a key role in absorbing atmospheric energy and directing it to the ice of debris-covered glaciers, but the spatial and temporal distribution of these features is not well documented. We analyse 172 Landsat TM/ETM+ scenes for the period 1999–2013 to identify thawed supraglacial ponds for the debris-covered tongues of five glaciers in the Langtang Valley of Nepal. We apply an advanced atmospheric correction routine (Landcor/6S) and use band ratio and image morphological techniques to identify ponds and validate our results with 2.5 m Cartosat-1 observations. We then characterize the spatial, seasonal and interannual patterns of ponds. We find high variability in pond incidence between glaciers (May–October means of 0.08–1.69% of debris area), with ponds most frequent in zones of low surface gradient and velocity. The ponds show pronounced seasonality, appearing in the pre-monsoon as snow melts, peaking at the monsoon onset at 2% of debris-covered area, then declining in the post-monsoon as ponds drain or freeze. Ponds are highly recurrent and persistent, with 40.5% of pond locations occurring for multiple years. Rather than a trend in pond cover over the study period, we find high interannual variability for each glacier after controlling for seasonality.

We examine repeat surface altimetry and radio echo observations of two supraglacial lakes in the percolation zone of the Greenland ice sheet to investigate the changes in firn conditions leading to lake formation and implications for meltwater storage within firn. Both lakes formed in 2011, when an anomalously high melt season was followed by low winter accumulation, resulting in reduced infiltration and storage in the near surface. The lakes expanded during the 2012 record melt season and retained liquid meltwater through the following winter. The lakes then contracted, with one lake slowly draining and refreezing and another rapidly draining to the subsurface. The lack of observable change in firn conditions surrounding the lakes indicates increased run-off in the near surface firn, likely along low-permeability ice layers formed during the previous melt seasons. This implies a reduced ability of the firn to absorb increased meltwater.

In this study, a modelling approach was used to investigate the cause of the diversion of erratic boulders from Mont Blanc and southern Valais by the Valais Glacier to the Solothurn lobe during the Last Glacial Maximum (LGM). Using the Parallel Ice Sheet Model, we simulated the ice flow field during the LGM, and analyzed the trajectories taken by erratic boulders from areas with characteristic lithologies. The main difficulty in this exercise laid with the large uncertainties affecting the paleo climate forcing required as input for the surface mass-balance model. In order to mimic the prevailing climate conditions during the LGM, we applied different temperature offsets and regional precipitation corrections to present-day climate data, and selected the parametrizations, which yielded the best match between the modelled ice extent and the geomorphologically-based ice-margin reconstruction. After running a range of simulations with varying parameters, our results showed that only one parametrization allowed boulders to be diverted to the Solothurn lobe during the LGM. This precipitation pattern supports the existing theory of preferential southwesterly advection of moisture to the alps during the LGM, but also indicates strongly enhanced precipitation over the Mont Blanc massif and enhanced cooling over the Jura Mountains.

We assess the importance of basal boundary conditions for transient simulations of Basin 3, Austfonna ice cap between January 1995 and December 2011 and for the surge starting in 2012 by carrying out simulations with the full-Stokes model Elmer/Ice and the vertically-integrated model BISICLES. Time-varying surface mass-balance data from the regional climate model HIRHAM5 are downscaled according to elevation. Basal friction coefficient is varied through time by interpolating between two data-constrained inversions of surface velocity fields, from 1995 and 2011. Evolution of the basal boundary condition appears to be much more important for mass discharge and the dynamic response of the fast flowing unit in Basin 3 than either model choice or the downscaling method for the surface mass balance. In addition, temporally linear extrapolation of the evolution of basal friction coefficient beyond the 2011 distribution could not reproduce the expansion of the acceleration observed in southern Basin 3 between January 2012 and June 2013. This implies that changes in basal friction patterns, and in turn basal processes that are not currently represented in either model, are among the most important factors for the 2012 acceleration.

Flow-frictional resistance at the base of glaciers and ice sheets is strongly linked to subglacial water pressure. Understanding the physical mechanisms that govern meltwater fluxes in subglacial channels is hence critical for constraining variations in ice flow. Previous mathematical descriptions of soft-bed subglacial channels assume a viscous till rheology, which is inconsistent with laboratory data and the majority of field studies. Here, we use a grain-scale numerical formulation coupled to pore-water dynamics to analyze the structural stability of channels carved into soft beds. Contrary to the soft-bed channel models assuming viscous till rheology, we show that the flanks of till channels can support substantial ice loads without creep closure of the channel, because the sediment has finite frictional strength. Increased normal stress on the channel flanks causes plastic failure of the sediment, and the channel rapidly shrinks to increase the ice-bed contact area. We derive a new parameterization for subglacial channelized flow on soft beds and show that channel dynamics are dominated by fluvial erosion and deposition processes with thresholds linked to the plastic rheology of subglacial tills. We infer that the described limits to channel size may cause subglacial drainage to arrange in networks of multiple closely spaced channels.

Petermann Gletscher drains ~4% of the Greenland ice sheet (GrIS) area, with ~80% of its mass loss occurring by basal melting of its ice shelf. We use a high-resolution coupled ocean and sea-ice model with a thermodynamic glacial ice shelf to diagnose ocean-controlled seasonality in basal melting of the Petermann ice shelf. Basal melt rates increase by ~20% in summer due to a seasonal shift in ocean circulation within Nares Strait that is associated with the transition from landfast sea ice to mobile sea ice. Under landfast ice, cold near-surface waters are maintained on the eastern side of the strait and within Petermann Fjord, reducing basal melt and insulating the ice shelf. Under mobile sea ice, warm waters are upwelled on the eastern side of the strait and, mediated by local instabilities and eddies, enter Petermann Fjord, enhancing basal melt down to depths of 200 m. The transition between these states occurs rapidly, and seasonal changes within Nares Strait are conveyed into the fjord within the same season. These results suggest that long-term changes in the length of the landfast sea-ice season will substantially alter the structure of Petermann ice shelf and its contribution to GrIS mass loss.

New experiments reveal that the flexural strength of fresh water, columnar-grained ice loaded normal to the columns may be increased by about a factor of two upon reversed cyclic loading at a frequency of ~0.1 Hz at stress amplitudes of 1.3–2.6 MPa. The effect is explained in terms of a reduction in deformation-induced stress concentration through the activation of grain boundary sliding, which is evident through boundary decohesion.

Basal slip accounts for a large part of the flow of ice streams draining ice from Antarctica and Greenland into the ocean. Therefore, an appropriate representation of basal slip in ice flow models is a prerequisite for accurate sea level rise projections. Various friction laws have been proposed to describe basal slip in models. Here, we compare the influence on grounding line (GL) dynamics of four friction laws: the traditional Weertman law and three effective pressure-dependent laws, namely the Schoof, Tsai and Budd laws. It turns out that, even when they are tuned to a common initial reference state, the Weertman, Budd and Schoof laws lead to thoroughly different steady-state positions, although the Schoof and Tsai laws lead to much the same result. In particular, under certain circumstances, it is possible to obtain a steady GL located on a reverse slope area using the Weertman law. Furthermore, the predicted transient evolution of the GL as well as the projected contributions to sea level rise over a 100-year time horizon vary significantly depending on the friction law. We conclude on the importance of choosing an appropriate law for reliable sea level rise projections and emphasise the need for a coupling between ice flow models and physically based subglacial hydrological models.

This study investigates the evolution of Greenland ice sheet flux focusing on five of the main fast flowing regions (Petermann glacier, North East Greenland Ice Stream, Kangerdlugssuaq glacier, Helheim glacier and Jakobshavn glacier) in response to 20th and 21st century climate change. A hybrid (shallow ice and shallow shelf) ice-sheet model (ISM) is forced with the combined outputs of a set of seven CMIP5 models and the regional climate model MAR. The ISM simulates the present-day ice velocity pattern, topography and surface mass balance (SMB) in good agreement with observations. Except for the Kangerdlugssuaq glacier, over the 21st century all the fast-flowing areas have exhibited a decrease in ice flux as a result of a negative SMB rather than dynamical changes. Only the fronts of Kangerdlugssuaq and Helheim glaciers have shown an interannual variability driven by dynamical rather than climate changes. Finally, the results predict a substantial inland ice margin retreat by the end of the 21st century, especially along the northern coasts.

Glacial earthquakes are slow earthquakes of magnitude M~5 associated with major calving events at near-grounded marine-terminating glaciers. These globally detectable earthquakes provide information on the grounding state of outlet glaciers and the timing of large calving events. Seismic source modeling of glacial earthquakes provides information on the size and orientation of forces associated with calving events. We compare force orientations estimated using a centroid-single-force technique with the calving-front orientations of the source glaciers at or near the time of earthquake occurrence. We consider earthquakes recorded at four glaciers in Greenland – Kangerdlugssuaq Glacier, Helheim Glacier, Kong Oscar Glacier, and Jakobshavn Isbræ – between 1999 and 2010. We find that the estimated earthquake force orientations accurately represent the orientation of the calving front at the time of the earthquake, and that seismogenic calving events are produced by a preferred section of the calving front, which may change with time. We also find that estimated earthquake locations vary in a manner consistent with changes in calving-front position, though with large scatter. We conclude that changes in glacial-earthquake source parameters reflect true changes in the geometry of the source glaciers, providing a means for identifying changes in glacier geometry and dynamics that complements traditional remote-sensing techniques.

We present the first dedicated study into the phenomenon of ice sails. These are clean ice structures that protrude from the surface of a small number of debris-covered glaciers and can grow to heights of over 25 m. We draw together what is known about them from the academic/exploration literature and then analyse imagery. We show here that ice sails can develop by one of two mechanisms, both of which require clean ice to become surrounded by debris-covered ice, where the debris layer is shallow enough for the ice beneath it to melt faster than the clean ice. Once formed, ice sails can persist for decades, in an apparently steady state, before debris layer thickening eventually causes a reversal in the relative melt rates and the ice sails decay to merge back with the surrounding glacier surface. We support our image-based analysis with a surface energy-balance model and show that it compares well with available observations from Baltoro Glacier in the Karakoram. A sensitivity analysis of the model is performed and confirms the results from our empirical study that ice sails require a relatively high evaporative heat flux and/or a relatively low sensible heat flux in order to exist.

An anomalously slight glacier mass gain during 2000 to the 2010s has recently been reported in the Karakoram region. However, to date, no investigations of the region-wide glacier mass balance in the Karakoram prior to 2000 have been reported, leaving a knowledge gap for assessing glacier responses to climate change. We calculated elevation and mass change using DEMs generated from KH-9 images acquired during 1973–1980 and the 1 arc-second SRTM DEM. We find a slight mass loss of −0.09 ± 0.03 m w.e. a−1 (12 366 km2) for 1973–2000, which is less negative than the global average rate for 1971–2009 (−0.31 ± 0.19 m w.e. a−1). Within the Karakoram, the glacier change patterns are spatially and temporally heterogeneous. In particular, a nearly stable state in the central Karakoram (−0.04 ± 0.05 m w.e. a−1 during the period 1974–2000) implies that the Karakoram anomaly dates back to the 1970s. Combined with the previous studies, we conclude that the Karakoram glaciers as a whole were in a nearly balanced state during the 1970s to the 2010s.

Temporal variations of the radionuclide 10Be are broadly synchronous across the globe and thus provide a powerful tool to synchronize ice core chronologies from different locations. We compared the 10Be record of the Akademii Nauk (AN) ice core (Russian Arctic) for the time period CE 1590–1950 to the 10Be records of two well-dated Greenland ice cores (Dye3 and NGRIP). A high correlation (r = 0.59) was found between the AN and Dye3 records whereas the correlation with NGRIP was distinctly lower (r = 0.45). Sources of deviations may include local fluctuations in the deposition of 10Be due to changes in the precipitation patterns, and artefacts due to the core-sampling strategy. In general, the existing age model was validated, confirming the AN ice core to be a unique and well-dated source of palaeoclimate parameters for the Russian Arctic. We further used numerical simulations to test the influence of the core-sampling strategy on the results and derived an optimized sampling strategy for the deeper parts of the ice core.

Light absorbing impurities (LAI) initiate powerful snow albedo feedbacks, yet due to a scarcity of observations and measurements, LAI radiative forcing is often neglected or poorly constrained in climate and hydrological models. To support physically-based modeling of LAI processes, daily measurements of dust and black carbon (BC) stratigraphy, optical grain size, snow density and spectral albedo were collected over the 2013 ablation season in the Rocky Mountains, CO. Surface impurity concentrations exhibited a wide range of values (0.02–6.0 mg g−1 pptw) with 98% of mass being deposited by three episodic dust events in April. Even minor dust loading initiated albedo decline, and the negative relationship between dust concentrations and albedo was log-linear. As melt progressed, individual dust layers coalesced and emerged at the snow surface, with minimal mass loss to meltwater scavenging. The observations show that the convergence of dust layers at the surface reduced albedo to 0.3 and snow depth declined ~50% faster than other years with similar depth but less dust. The rapid melt led to an unexpected reduction in both grain size and density in uppermost surface layers. BC concentrations co-varied with dust concentrations but were several orders of magnitude lower (<1–20 ppb).

Discharge and water temperature measurements in the Skaftá river and measurements of the lowering of the ice over the subglacial lake at the western Skaftá cauldron, Vatnajökull, Iceland, were made during a rapidly rising glacial outburst flood (jökulhlaup) in September 2006. Outflow from the lake, flood discharge at the glacier terminus and the transient subglacial volume of floodwater during the jökulhlaup are derived from these data. The 40 km long initial subglacial path of the jökulhlaup was mainly formed by lifting and deformation of the overlying ice, induced by water pressure in excess of the ice overburden pressure. Melting of ice due to the heat of the floodwater from the subglacial lake and frictional heat generated by the dissipation of potential energy in the flow played a smaller role. Therefore this event, like other rapidly rising jökulhlaups, cannot be explained by the jökulhlaup theory of Nye (1976). Instead, our observations indicate that they can be explained by a coupled subglacial-sheet–conduit mechanism where essentially all of the initial flood path is formed as a sheet by the propagation of a subglacial pressure wave.