The increasingly negative mass balance of the Greenland ice sheet (GrIS) over the last ~25 years has been associated with enhanced surface melt and increased ice loss from marine-terminating outlet glaciers. Accelerated retreat during 2000–2010 was concentrated in the southeast and northwest sectors of the ice sheet; however, there was considerable spatial and temporal variability in the timing and magnitude of retreat both within and between these regions. This behaviour has yet to be quantified and compared for all glaciers in both regions. Furthermore, it is unclear whether retreat has continued after 2010 in the northwest, and whether the documented slowdown in the southeast post-2005 has been sustained. Here, we compare spatial and temporal patterns of frontal change in the northwest and southeast GrIS, for the period 2000–2015. Our results show near-ubiquitous retreat of outlet glaciers across both regions for the study period; however, the timing and magnitude of inter-annual frontal position change is largely asynchronous. We also find that since 2010, there is continued terminus retreat in the northwest, which contrasts with considerable inter-annual variability in the southeast. Analysis of the role of glacier-specific factors demonstrates that fjord and bed geometry are important controls on the timing and magnitude of glacier retreat.

The stable oxygen isotope composition, major ions and isotopic compositions of strontium (Sr), neodymium (Nd) and lead (Pb) in insoluble dust from recent surface snow samples along the transect from the Zhongshan and Progress stations (located on the Amery Ice Shelf, East Antarctica) to Dome A (Summit, Antarctica) were analysed. No previous isotopic fingerprinting studies have been conducted for this transect. These data were used to document the dust provenances in Antarctica along the transect up to the highest site, Dome A, for the first time. The insoluble dust in snow samples along the coast displays an overall crust-line isotopic signature that is characterised by highly radiogenic 87Sr/86Sr values and less radiogenic 143Nd/144Nd values. These signatures are comparable with those of samples collected near the ice-free areas of the Zhongshan and Progress stations. Spatial differences are statistically significant along the transect, and the Sr, Nd and Pb isotope components in insoluble dust from two continuous snow samples at Dome A exhibit marked differences, indicating that additional dust reaches the East Antarctic Plateau. The isotopic characteristics of insoluble dust from this transect indicate that the long-distance natural dust and anthropogenic pollutants in these samples primarily originate from Australia.

Haxilegen Glacier No.51 (43.731°N, 84.391°E; CN5Y741C0051) is located in the Kuytun river basin, Erenharbirga range, eastern Tien Shan. This study presents the annual mass balance of Haxilegen Glacier No.51 for 7 hydrological years and uses a temperature-index and an accumulation model to reconstruct the annual mass balance from 1999 to 2015. The model is calibrated against annual altitudinal mass-balance measurements and then applied to the period with no measurements. We find an accumulated mass balance of −6.06 ± 0.88 m w.e.a−1 over the period of 16 hydrological years, with an average annual value of −0.32 ± 0.22 m w.e.a−1. The mean glacier-wide annual, summer and winter balances for 1999 to 2015 are −0.37, −0.54 and 0.16 ± 0.22 m w.e.a−1, respectively, with a high correlation coefficient (r = 0.95, p < 0.001) between annual balance and summer balance. The calculated mass-balance sensitivity of the glacier to temperature is −0.51 m w.e.a−1 °C−1 and to precipitation is 0.08 m w.e.a−1 for a 10% increase. The sensitivity of seasonal mass balance indicates that temperature during the melt season (May–August) and annual precipitation are the major contributors to mass-balance fluctuation.

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.

Accurate quantification of rates of glacier mass loss is critical for managing water resources and for assessing hazards at ice-clad volcanoes, especially in arid regions like southern Peru. In these regions, glacier and snow melt are crucial dry season water resources. In order to verify previously reported rates of ice area decline at Nevado Coropuna in Peru, which are anomalously rapid for tropical glaciers, we measured changes in ice cap area using 259 Landsat images acquired from 1980 to 2014. We find that Coropuna Ice Cap is presently the most extensive ice mass in the tropics, with an area of 44.1 km2, and has been shrinking at an average area loss rate of 0.409 km2 a−1 (~0.71% a−1) since 1980. Our estimated rate of change is considerably lower than previous studies (1.4 km2 a−1 or ~2.43% a−1), but is consistent with other tropical regions, such as the Cordillera Blanca located ~850 km to the NW (~0.68% a−1). Thus, if glacier recession continues at its present rate, our results suggest that Coropuna Ice Cap will likely continue to contribute to water supply for agricultural and domestic uses until ~2120, which is nearly 100 years longer than previously predicted.

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.

Over recent decades, Greenland ice sheet surface melt has shown an increase both in intensity and spatial extent. Part of this water probably reaches the bed and can enhance glacier speed, advecting a larger volume of ice into the ablation area. In the context of a warming climate, this mechanism could contribute to the future rate of thinning and retreat of land-terminating glaciers of Greenland. These changes in ice flow conditions will in turn influence surface crevassing and thus the ability of water to reach the bed at higher elevations. Here, using a coupled basal hydrology and prognostic ice flow model, the evolution of a Greenland-type glacier subject to increasing surface melt is studied over a few decades. For different scenarios of surface melt increase over the next decades, the evolution of crevassed areas and the ability of water to reach the bed is inferred. Our results indicate that the currently observed crevasse distribution is likely to extend further upstream which will allow water to reach the bed at higher elevations. This will lead to an increase in ice flux into the ablation area which, in turn, accelerates the mass loss of land-terminating glaciers.

Stable isotopes of water (δ18O and δ2H) were measured in the debris-laden ice underlying an Antarctic blue ice moraine, and in adjoining Law Glacier in the central Transantarctic Mountains. Air bubble content and morphology were assessed in shallow ice core samples. Stable isotope measurements plot either on the meteoric waterline or are enriched from it. The data cluster in two groups: the ice underlying the moraine has a δ2H:δ18O slope of 5.35 ± 0.92; ice from adjoining portions of Law Glacier has a slope of 6.69 ± 1.39. This enrichment pattern suggests the moraine's underlying blue ice entrained sediment through refreezing processes acting in an open system. Glaciological conditions favorable to warm-based sediment entrainment occur 30–50 km upstream. Basal melting and refreezing are further evidenced by abundant vapor figures formed from internal melting of the ice crystals. Both the moraine ice and Law Glacier are sufficiently depleted of heavy isotopes that their ice cannot be sourced locally, but instead must be derived from far-field interior regions of the higher polar plateau. Modeled ice flow speeds suggest the ice must be at least 80 ka old, with Law Glacier's ice possibly dating to OIS 5 and moraine ice older still.

Autonomous temperature data loggers were used to measure the temperature profile within a growing ice cover and in the water below. The ice formed under natural conditions over the pond. We observed the presence of distinct layers of gas bubbles throughout the ice thickness. Temperature measurements allowed us to determine growth rates (μm s−1) and cooling rates (°C s−1) of the ice and demonstrated that these bubble layers formed during the peak ice growth rates from 0.58 to 0.92 µm s−1. The growth rates, leading to the formation of layers of bubbles, were more than an order of magnitude lower than for bubbles produced in controlled laboratory conditions (from 3 to 80 µm s−1). This observation introduces the possibility that solid impurities play a role in natural waters and that they must lower the limit of growth rates required for bubble occurrence. Data revealed a decrease in ice growth rates while cooling rates increased. We interpret this observation as an effect of the heat flux from the water to the ice (8.34–34.11 W m−2), and of gas concentration changes in the water below. Calculations of the ice thickness using traditional methods showed the necessity to include the heat flux from the water to the ice and the effect of gas bubbles within the ice and near the ice–water interface.

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.

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 non-linearly viscous ice flow in the vicinity of a glacier terminus, an observation region, depends crucially on the upstream flow as well as on the local surface and bed conditions. The former requires a likely complex solution of the balance laws and boundary conditions for the complete glacier. However, if the profile and downstream surface tangential velocity in the observation region are measured at an observation time t = 0, and a two-dimensional flow approximation is satisfactory, the complete stress and velocity fields satisfying local reduced model equations in the observation region at time t = 0 can be determined by asymptotic expansions in upstream distance from the (moving) terminus. Thus the full strain-rate and stress tensors are determined without prescribing the basal conditions. The terminus velocity is determined in terms of the net accumulation or melt flux and surface velocity at the terminus, with bounds for advance or retreat. The analysis and illustration are presented for a plane flow approximation.

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 thermal state at the bed of a large ice sheet is a critical boundary condition governing its future evolution. Radar surveys provide an opportunity for direct but remote observation of the ice-sheet bed, and therefore offer a means of constraining numerical ice-sheet models at the ice–bed interface. Here we have processed results of radar surveys of the Greenland Ice Sheet undertaken by the Program for Arctic Regional Climate Assessment (PARCA) between 1999 and 2003, to explore this opportunity. We consider the robustness of the measurements in the context of uncertain dielectric losses in the ice sheet, concluding that the observed radar signal characteristics reflect the character of the bed itself rather than that of uncertain englacial absorption. However, the identification of thaw is restricted to areas where subglacial water has sufficient depth to influence the radar reflection. We derive a map of inferred areas of subglacial thaw, and compare our results with other studies predicting regions with temperate bed. We show that in many areas the radar inferences of ponded water lie within areas predicted to be thawed by modelling and radiostratigraphy. There is clear disagreement in certain areas, suggesting the presence of high geothermal flux anomalies.

Past atmospheric pollution can be reconstructed from ice core trace element records retrieved from mountain glaciers. However, the current global temperature increase can result in post-depositional melt processes, significantly altering the originally stored information. Here, we present a comprehensive study on the behaviour of 35 trace elements (TEs) during meltwater percolation in a high-Alpine ice core segment from upper Grenzgletscher, Switzerland. Some TEs revealed significant concentration depletion, whereas others were well preserved depending on their water solubility and location at the grain scale. TEs present in insoluble minerals, typically enriched at grain boundaries, were found to be mostly preserved because their insolubility in water results in immobility with meltwater percolation. Water-soluble TEs revealed a variable meltwater-mobility. Whereas ultra-TEs tend to be preserved, likely due to incorporation into the ice lattice, abundant TEs are prone to relocation from grain-boundary regions. We propose that at Alpine sites, Ag, Al, Bi, Cu, Cs, Fe, Li, Mo, Pb, Rb, Sb, Th, Tl, U, V, W, Zr and the rare-earth elements may still be applicable as robust environmental proxies even if partial melting occurred, whereas Ba, Ca, Cd Co, Mg, Mn, Na, Ni, Sr and Zn are prone to significant depletion.

To obtain the effect of velocity and structural natural frequency (structural stiffness) on ice failure, an extended dynamic Van der Pol-based single degree-of-freedom ice/structure interaction model is developed. Three basic modes of response were reproduced: intermittent crushing, frequency lock-in and continuous crushing. Further analysis on the physical mechanism of ice/structure interaction is presented on the basis of feedback mechanism and energy mechanism, respectively. Internal effect and external effect from ice and structure were both explained in the feedback branch. Based on reproduced results, energy exchanges at different configurations are computed from the energy conservation using the first law of thermodynamics. A general conclusion on the predominant type of vibration when the ice velocity increases during the interaction process is forced, self-excited and forced in each of the three modes of responses. Ice force variations also show that there is more impulse energy during the lock-in range. Moreover, ice-induced vibration demonstrates an analogy of friction-induced self-excited vibration. Finally, the similarity between strain-stress curve and Stribeck curve shows that static and kinetic friction force variations are attributed to ice force characteristic, and can be used to explain the lower effective pressure magnitude during continuous crushing than the peak pressure during intermittent crushing.

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.

Snow slab avalanches are caused by cracks forming and propagating in a weak snow layer below a cohesive slab. The gradual damage process leading to the formation of the initial failure within the weak layer (WL) is still not entirely understood. To this end, we designed a novel test apparatus that allows performing loading experiments with large snow samples (0.25 m2) including a WL at different loading rates and simultaneously monitoring the acoustic emissions (AE) response. By analyzing the AE generated by micro-cracking, we studied the evolution of the damage process preceding snow failure. At fast loading rates, the exponent of the AE energy distribution (b-value) gradually changed, and both the energy rate and the inverse waiting time increased exponentially with increasing load. These changes in AE signature indicate a transition from small to large events and an acceleration of the damage processes leading to brittle failure. For the experiments at slow loading rate, these changes in the AE signature were not or only partially present, even if the sample failed, indicating a different evolution of the damage process. The observed characteristics in AE response provide new insights on how to model snow failure as a critical phenomenon.

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.

The controls on rapid surface lake drainage on the Greenland ice sheet (GrIS) remain uncertain, making it challenging to incorporate lake drainage into models of GrIS hydrology, and so to determine the ice-dynamic impact of meltwater reaching the ice-sheet bed. Here, we first use a lake area and volume tracking algorithm to identify rapidly draining lakes within West Greenland during summer 2014. Second, we derive hydrological, morphological, glaciological and surface-mass-balance data for various factors that may influence rapid lake drainage. Third, these factors are used within Exploratory Data Analysis to examine existing hypotheses for rapid lake drainage. This involves testing for statistical differences between the rapidly and non-rapidly draining lake types, as well as examining associations between lake size and the potential controlling factors. This study shows that the two lake types are statistically indistinguishable for almost all factors investigated, except lake area. Thus, we are unable to recommend an empirically supported, deterministic alternative to the fracture area threshold parameter for modelling rapid lake drainage within existing surface-hydrology models of the GrIS. However, if improved remotely sensed datasets (e.g. ice-velocity maps, climate model outputs) were included in future research, it may be possible to detect the causes of rapid drainage.