Accelerated discharge through marine-terminating outlet glaciers has been a key component of the rapid mass loss from Arctic glaciers since the 1990s. However, glacier retreat and its climatic controls have not been assessed at the pan-Arctic scale. Consequently, the spatial and temporal variability in the magnitude of retreat, and the possible drivers are uncertain. Here we use remotely sensed data acquired over 273 outlet glaciers, located across the entire Atlantic Arctic (i.e. areas potentially influenced by North Atlantic climate and/or ocean conditions, specifically: Greenland, Novaya Zemlya, Franz Josef Land and Svalbard), to demonstrate high-magnitude, accelerating and near-ubiquitous retreat between 1992 and 2010. Overall, mean retreat rates increased by a factor of 3.5 between 1992 and 2000 (−30.5 m a−1) and 2000–10 (−105.8 m a−1), with 97% of the study glaciers retreating during the latter period. The Retreat was greatest in northern, western and south-eastern Greenland and also increased substantially on the Barents Sea coast of Novaya Zemlya. Glacier retreat showed no significant or consistent relationship with summer air temperatures at decadal timescales. The rate of frontal position change showed a significant, but weak, correlation with changes in sea-ice concentrations. We highlight large variations in retreat rates within regions and suggest that fjord topography plays an important role. We conclude that marine-terminating Arctic outlet glaciers show a common response of rapid and accelerating retreat at decadal timescales.

Quantifying the biogeochemical cycling of carbon in glacial ecosystems is of great significance for regional, and potentially global, carbon flow estimations. The concentration and quality of organic carbon (OC) is an important indicator of biogeochemical and physical processes that prevail in an ice-sheet ecosystem. Here we determine the content and quality of OC in debris from the surface of the Greenland ice sheet (GrIS) using microscopic, chromatographic, spectrophotometric and high-temperature combustion techniques. The total OC content in the debris increased with distance from the edge of the ice sheet, from virtually zero to >6% dry weight at 50 km inland, and there was a peak in the carbohydrate proportion and the microbial abundance at ∼6km inland. The highest (galactose + mannose)/(arabinose + xylose) ratios, indicating maximum autochthonous microbial production, were found at >10km inland. We propose that three key processes influence the carbon cycling on the GrIS: aeolian input of microbial inoculum and nutrients, in situ biological C transformation and the wash-away of supraglacial debris by meltwaters. We show that all these processes have significant spatial variability. While the total OC content of the debris on the ice sheet is probably controlled by the physical processes of wind transport and wash-away by meltwater, the microbial abundance and the quantity of the labile cell-contained OC within the debris is likely to be driven by the balance between the wash-away and the microbial productivity.

Eddy covariance data collected over a horizontal surface on the largest ice body on Kilimanjaro, Tanzania, over 26–29 July 2005 were used to assess the uncertainty of calculating sublimation with a surface energy balance (SEB) model. Data required for input to the SEB model were obtained from an existing automatic weather station. Surface temperatures that were solved iteratively by the SEB model were used to compute emitted longwave radiation, turbulent heat fluxes using the aerodynamic bulk method and the subsurface heat flux. Roughness lengths for momentum and temperature, which were found to be the most important input parameters controlling the magnitude of modelled (bulk method) turbulent heat fluxes, were obtained using eddy covariance data. The roughness length for momentum was estimated to be 1.7×10–3 m, while the length for temperature was one order of magnitude smaller. Modelled sensible and latent heat fluxes (bulk method) compared well to eddy covariance data, with root-mean-square differences between 3.1 and 4.8 Wm–2 for both turbulent heat fluxes. Modelled sublimation accounted for about 90% of observed ablation, confirming that mass loss by melting is much less important than sublimation on the horizontal surfaces of the remaining plateau glaciers on Kilimanjaro.

Kyagar glacier is located in the Chinese Karakoram mountains. The glacier tongue entirely blocks the riverbed in the upper Shaksgam valley and impounds a glacial lake, which was the source of several violent and disastrous glacial lake outburst floods (GLOFs). A GLOF early warning system was implemented between 2011 and 2013. We present an integrative analysis of the hazard potential of Kyagar lake, taking into account the ice flow dynamics of Kyagar glacier as well as the recent surface mass-balance response to climate change. Comparison of two high-resolution digital elevation models (DEMs) for the ice dam shows surface lowering rates of >5ma– 1 between 2002 and 2011, leading to a significant reduction in the maximum potential lake volume. However, two DEMs covering the entire glacier for the period 2000–10 indicate mass gains in its central part, and flow speed measurements show an acceleration in this region. This pattern of local ice-thickness changes combined with varying ice flow velocities is typical for surge-type glaciers. The velocity of the glacier surface and of the ice dam between 2011 and 2012 are analyzed at high temporal and spatial resolution, based on feature tracking of synthetic aperture radar (SAR) images.

The presence of leads with open water or thin ice is an important feature of the Arctic sea ice cover. Leads regulate the heat, gas and moisture fluxes between the ocean and atmosphere and are areas of high ice growth rates during periods of freezing conditions. Here, an algorithm providing an automatic lead detection based on synthetic aperture radar images is described that can be applied to a wide range of Sentinel-1 scenes. By using both the HH and the HV channels instead of single co-polarised observations the algorithm is able to classify more leads correctly. The lead classification algorithm is based on polarimetric features and textural features derived from the grey-level co-occurrence matrix. The Random Forest classifier is used to investigate the importance of the individual features for lead detection. The precision–recall curve representing the quality of the classification is used to define threshold for a binary lead/sea ice classification. The algorithm is able to produce a lead classification with more that 90% precision with 60% of all leads classified. The precision can be increased by the cost of the amount of leads detected. Results are evaluated based on comparisons with Sentinel-2 optical satellite data.

The RADARSAT-1 Antarctic Mapping Project (RAMP) is a collaboration between NASA and the Canadian Space Agency to map Antarctica using synthetic aperture radar (SAR). The first Antarctic Mapping Mission (AMM-1) was successfully completed in October 1997. Data from the acquisition phase of the 1997 campaign have been used to achieve the primary goal of producing the first high-resolution SAR image map of Antarctica. The Modified Antarctic Mapping Mission (MAMM) occurred during the fall of 2000. The acquisition strategy concentrated on collecting highest- resolution RADARSAT-1 data of Antarctica’s fast glaciers for change detection and feature-retracking estimates of surface velocity. Additionally, extensive data were acquired for interferometric analysis over the entire viewable region, which extends north of 80.1˚ S latitude. This paper summarizes the goals and strategy behind MAMM. It goes on to discuss ice-sheet margin changes observed on several ice shelves around the continent. Margin changes are documented by comparing AMM-1 and MAMM data with earlier datasets including European Remote-sensing Satellite-1 SAR imagery, Landsat imagery, the Antarctic Digital Database (version 1) coastline and Declassified Intelligence Satellite Photography. Analysis reveals a complex pattern of ice-margin advance and retreat without indicating any systematic behavior in ice-sheet extent about the ice-sheet perimeter.

The morphology, sediments, and processes associated with the construction of a moraine along the western margin of the ice shelf in George VI Sound, Antarctica, are discussed. The moraine occurs as a double ridge where the ice sheet grounds against promontories on Alexander Island and is approximately horizontal over a distance of 120 km. It consists of exotic rock debris carried into the ice shelf by Antarctic Peninsula glaciers and local rock debris derived from the grounding line on Alexander Island. As the coast steepens, so the proportion of exotic rocks increases. The transport of basal material from the peninsula implies that there can be little bottom melting beneath this part of the ice shelf. The moraine is modified by streams and marginal lakes which periodically drain into and through the ice shelf. Tidal lakes are impounded against the ice shelf in shallower embayments and consist of fresh water overlying sea-water. A conceptual model of the moraine is developed and may help to explain some features of puzzling horizontal moraines found in formerly glaciated areas.

Blue-ice areas (BIAs) and their geographical distribution in Antarctica were mapped using Landsat-7 ETM+ images with 15 m spatial resolution obtained during the 1999–2003 austral summers and covering the area north of 82.5° S, and a snow grain-size image of the MODIS-based Mosaic of Antarctica (MOA) dataset with 125 m grid spacing acquired during the 2003/04 austral summer from 82.5°S to the South Pole. A map of BIAs was created with algorithms of thresholds based on band ratio and reflectance for ETM+ data and thresholds based on snow grain size for the MOA dataset. The underlying principle is that blue ice can be separated from snow or rock by their spectral discrepancies and by different grain sizes of snow and ice. We estimate the total area of BIAs in Antarctica during the data acquisition period is 234 549 km2, or 1.67% of the area of the continent. Blue ice is scattered widely over the continent but is generally located in coastal or mountainous regions. The BIA dataset presented in this study is the first map covering the entire Antarctic continent sourced solely from ETM+ and MODIS data. This dataset can potentially benefit other studies in glaciology, meteorology, climatology and paleoclimate, meteorite collection and airstrip site selection.

We have recently completed an analysis that examines in detail the spatial and temporal variations in global sea-ice coverage from 26 October 1978, through 20 August 1987. The sea-icemeasurements we analyzed are derived from data collected by a multispectral, dual-polarized, constant incidence-angle microwave imager, the Scanning Multichannel Microwave Radiometer (SMMR) on board the NASA Nimbus 7 satellite. The characteristics of the SMMR have permitted a more accurate calculation of total sea-ice concentrations (fraction of ocean area covered by sea ice) than earlier single-channel instruments and, for the first time, a determination of both multiyear sea-ice concentrations and physical temperatures of the sea-ice pack. An estimate of the SMMR wintertime total ice concentration accuracy of ± 7% in both hemispheres has been obtained. As this is an improvement over the estimated accuracies of previous microwave sensors, we are able to present improved calculations of the sea-ice extents (areas enclosed by the 15% ice concentration boundaries), sea-ice concentrations, and open-water areas within the ice margins. This analysis will be published in a book, Arctic and Antarctic sea ice, 1978–1987: satellite passive microwave observations and analysis, due for publication in1992. Some highlights from the analysis are presented in this paper.

Polynyas and leads are key elements of the wintertime Arctic sea-ice cover. They play a crucial role in surface heat loss, potential ice formation and consequently in the seasonal sea-ice budget. While polynyas are generally sufficiently large to be observed with passive microwave satellite sensors, the monitoring of narrow leads requires the use of data at a higher spatial resolution. We apply and evaluate different lead segmentation techniques based on sea-ice surface temperatures as measured by the Moderate Resolution Imaging Spectroradiometer (MODIS). Daily lead composite maps indicate the presence of cloud artifacts that arise from ambiguities in the segmentation process and shortcomings in the MODIS cloud mask. A fuzzy cloud artifact filter is hence implemented to mitigate these effects and the associated potential misclassification of leads. The filter is adjusted with reference data from thermal infrared image sequences, and applied to daily MODIS data from January to April 2008. The daily lead product can be used to deduct the structure and dynamics of wintertime sea-ice leads and to assess seasonal divergence patterns of the Arctic Ocean.

The distribution of basal traction on a transect along Pine Island Glacier, West Antarctica, is estimated by inverting observed surface velocities using a control method and a simple numerical stream-flow model. This model calculates the horizontal flow along a transect, based on the assumptions that the horizontal flow is independent of ice depth and that the driving stresses are balanced by resistive forces at the glacier bed and margin and by gradients in longitudinal stress. Basal traction is assumed to be linearly related to the basal velocity. For the lateral shear traction a parameterization based on an inversion of Glen’s flow law is used. The application of the control method allows us to calculate the set of model parameters (e.g. the basal friction coefficient) that gives the best fit between modelled andobserved surface velocities. The model is used to investigate the stress regime of Pine Island Glacier, in particular to estimate the importance of basal, lateral and longitudinal stresses relative to each other. In the flat region just behind the grounding line, basal drag, lateral drag and the longitudinal stress gradient are the same order of magnitude. In the steep region up-glacier from the grounding line, the driving stresses are highest and balanced predominantly by basal resistive stresses. Further upstream, in the trunk of the glacier, lateral and basal drag predominate.

Several hundred visible and thermal infrared satellite images of Antarctica’s southeast Amundsen Sea from 1986 to 2011, combined with aerial observations in 2009, show a strong inverse relation between prominent curvilinear surface depressions and the underlying basal morphology of the outer Pine Island Glacier ice shelf. Shipboard measurements near the calving front reveal positive temperature, salinity and current anomalies indicative of melt-laden, deep-water outflows near and above the larger channel termini. These buoyant plumes rise to the surface and are expressed as small polynyas in the sea ice and thermal signatures in the open water. The warm upwellings also trace the cyclonic surface circulation in Pine Island Bay. The satellite coverage suggests changing modes of ocean/ice interactions, dominated by leads along the ice shelf through 1999, fast ice and polynyas from 2000 to 2007, and larger areas of open water since 2008.

In situ observations of precipitation days (days when snow or rain was reported in routine synoptic observations) from Faraday/Vernadsky station on the western side of the Antarctic Peninsula, and fields from the 40 year European Centre for Medium-Range Weather Forecasts re-analysis (ERA-40) project are used to investigate precipitation and atmospheric circulation changes around the Antarctic Peninsula. It is shown that the number of precipitation days is a good proxy for mean sea-level pressure (MSLP) over the Amundsen–Bellingshausen Sea. The annual total of precipitation days at the station has been increasing at a statistically significant rate of +12.4 days decade–1 since the early 1950s, with the greatest increase taking place during the summer and autumn. This is the time of year when the Southern Annular Mode (SAM) has experienced its greatest shift to a positive phase, with MSLP values decreasing in the Antarctic coastal zone. The lower pressures in the circumpolar trough have resulted in greater ascent and increased precipitation at Faraday/Vernadsky.

Hydrologic and hydrochemical studies have been conducted in Green Lakes Valley, Colorado Front Range, USA, above 3550 m since 1982. They show a classic seasonal hydrograph dominated by snowmelt and an earlier date for the start of spring flow and for peak flow over the period of record. This is consistent with patterns found at lower elevations in Colorado and throughout western North America. They also show an increasing trend in flows in September and October of 2.6 ±0.7 mm a−1 which is not found elsewhere and cannot be accounted for by increased autumn precipitation and the melting of surface ice. Because this late-season increase is not found at the highest elevations or in basins in which there is no evidence of permafrost, it seems best explained by the thawing of alpine permafrost at intermediate elevations. This is corroborated by an increase in the concentration of base cations and silica, and particularly in Ca2+ and SO4–, in the stream discharge starting in 2000. As with the physical hydrology, the geochemical signals have not been detected at the higher elevations in the basin, though they have previously been associated with streamflow from a small rock glacier in the valley. The combined evidence suggests the degradation of ice-rich permafrost on the north-facing slopes of the valley below 3700 m, where it has been detected at 3 m depth by geophysical surveys.

Glaciers in the Qilian mountains, located in the northeastern part of the Tibetan Plateau, constitute an important freshwater resource for downstream populations and natural systems. To enhance our understanding of the variability of the glaciers, temporally and spatially comprehensive information on them is needed. In this study, the glacier outlines of ~1990, ~2000 and ~2010 for the whole area were delineated in a semi-automated manner using band TM3/TM5 ratio images of Landsat ETM+ or TM scenes with the help of a merged ASTER GDEM/SRTM v4.1 digital elevation model. Combining our own results with those of previously published studies that span the period back to 1956, we found that the glacier area shrank by 30 ±8% from 1956 to 2010 and the shrinkage accelerated remarkably in the past two decades. The linear trends of annual air temperature and precipitation measured at weather stations within the glacierized areas were 0.03-0.05°Ca-1 (significant only after 2000) and 0.37-1.58 mm a–1 (not significant) respectively from 1961 to 2010. Glaciers shrank mainly due to the increasing temperature. Glaciers in the Qilian mountains are very unlikely to have experienced positive mass balance over the past decade. Moreover, given the trend toward higher temperatures, the glaciers in this region will continue to shrink.

Aerial photography has been used as a mapping tool in the Antarctic Peninsula region sinre the late 1920s. Following pioneering work by Wilkins in 1928, Ellsworth in 1934 and the British Graham Land Expedition in 1934-37, the Falkland Islands and Dependencies Aerial Survey Expedition carried out extensive aerial photography during the period 1955-57. Since then, many other aerial surveys have been carried out, and the result is an archive of aerial photography that, for some localities, spans 40 years. The production of maps both from different generations of photographs and satellite images has revealed many changes in the extent of ice cover with time. For example, changes in ice shelves such as the Wordie Ice Shelf, Larsen Ice Shelf and Müller Ice Shelf, are well recorded, and the termini of some glaciers have retreated. However, the most pervasive change is the consistent decline in the extent of small bodies of snow and ice. This paper shows how perennial snow or ice cover has decreased in the northern Marguerite Bay area, at 68°S. The correlation of the change with elevation and with climate records from Adelaide and Rothera research stations in the Antarctic Peninsula region is examined.

Water content, water distribution and variations with time of these quantities in the natural snow cover are important parameters for microwave remote-sensing studies, snow stability investigations and snow hydrology studies. The most promising method for snow wetness determination is the measurement of the dielectric constant at radio frequencies. Recent developments of electronic devices for long-term recording of snow wetness are reported. The measurement system consists of two parts: the tuning and display unit with the battery pack and the flat capacitive sensor. The system is controlled by a built-in low-power microcomputer, the operation frequency for the sensor is 20.00 MHz. Measurement data are displayed on an LCD and are also stored in a RAM; data transfer for a PC for further data processing is possible by a built-in V24 interface. Up to four individual sensors can be operated simultaneously. Field measurements of long-term variations in snow wetness during diurnal cycles and registrations of penetrating meltwater or rain-water waves with the formation of water shock fronts are presented.

The geomorphology of palaeo-ice-stream beds and the internal structure of underlying tills can provide important information about the subglacial conditions during periods of fast flow and quiescence. This paper presents observations from three-dimensional seismic data, revealing the geomorphology of buried beds of the Bjørnøyrenna (Bear Island Trough) ice stream, the main drainage outlet of the former Barents Sea ice sheet. Repeated changes in ice dynamics are inferred from the observed successions of geomorphic features. Megablocks, aligned in long chains parallel to inferred ice-stream flowlines, and forming dipping plates that are thrust one on top of another, are taken as evidence for conditions of compressive ice flow. Mega-scale glacial lineations (MSGL) and pull-apart of underlying sediment blocks suggest extensional flow. The observed pattern of megablocks and rafts overprinted by MSGL indicates a change in ice dynamics from a compressional to an extensional flow regime. Till stiffening, due to subglacial freezing, is the favoured mechanism for creating switches in sub-ice-stream conditions. The observed pattern of geomorphic features indicates that periods of slowdown or quiescence were commonly followed by reactivation and fast flow during several glaciations, suggesting that this may be a common behaviour of marine ice streams.

The model used by Lingle and Clark (1985) to approximate the deformation of the Earth under a single ice stream is adapted to the purposes of continent-scale ice-sheet simulation. The model combines a layered elastic spherical Earth (Farrell, 1972) with a viscous half-space overlain by an elastic plate lithosphere (Cathles, 1975). For the half-space model we identify a new mathematical formulation, essentially a time-dependent partial differential equation, which generalizes and improves upon the standard elastic plate lithosphere with relaxing asthenosphere model widely used in ice-sheet simulation. The new formulation allows a significantly faster numerical strategy, a spectral collocation method based directly on the fast Fourier transform. We verify this method by comparing to an integral formula for a disk load. We also demonstrate that the magnitudes of numerical errors made in approximating coupled ice-flow/Earth-deformation systems are significantly smaller than pairwise differences between several Earth models. Our implementation of the Lingle and Clark (1985) model offers important features of spherical, layered, self-gravitating, viscoelastic Earth models without the computational expense.

We show that the ice plain in the mouth of Whillans Ice Stream (formerly Ice Stream B), Antarctica, moves by stick–slip motion. During a spring-tide period, rapid motions regularly occur near high tide and during falling tide. This correlation is weaker during a neap-tide period when the tidal magnitudes are less. Precise timing of these motion events suggests that they propagate through the region with a mean velocity of 88 m s−1.We hypothesize that this speed is associated with the propagation of shear waves through a wet subglacial till. Motion events are also seen on more smoothly flowing floating ice. Event delays are very short between grounded and floating stations, suggesting the events propagate through the ice shelf as an elastic wave. We further hypothesize the events are caused by the interaction of a sticky bed, the accumulation of stored elastic strain through the compression of ice by upstream inflow, and tidal forcing. Motion events seem to be triggered either by reduction of vertical normal stresses at high tide or by the increase of shear stresses from sub-shelf ocean currents during falling tide. Event magnitudes are not related to the length of the preceding quiescent period, suggesting significant viscous dissipation within the till.