Ice-thickness measurements in Antarctic ice-shelf grounding zones are necessary for calculating the mass balance of individual catchments, but remain poorly constrained for most of the continent. We describe a new inverse modelling optimization approach to estimate ice thickness in the grounding zone of Antarctic outlet glaciers and ice shelves using spatial patterns of tide-induced flexure derived from differential interferometric synthetic aperture radar (InSAR). We demonstrate that the illposedness of the inverse formulation of the elastic-plate equations for bending can be controlled by regularization. In one dimension, the model recreates smooth, synthesized profiles of ice thickness from flexure information to within 1–2%. We test the method in two dimensions and validate it in the grounding zone of Beardmore Glacier, a major outlet glacier in the Transantarctic Mountains, using interferograms created from TerraSAR-X satellite imagery acquired in 2012. We compare our results with historic and modern ice-thickness data (radio-echo sounding from 1967 and ground-penetrating radar from 2010). We match both longitudinal and transverse thickness transects to within 50 m root-mean-square error using an effective Young’s modulus of 1.4 GPa. The highest accuracy is achieved close to the grounded ice boundary, where current estimates of thickness based on surface elevation measurements contain a systematic bias towards thicker ice.

This study presents the first complete glacier inventory of the Torngat Mountains, northern Labrador, Canada. In total, 195 glaciers and ice masses are identified, covering a total area of 24.5 ± 1.8 km2. Mapped ice masses range in size from 0.01 to 1.26 km2, with a median size of 0.08 km2. Ice masses have a median elevation of 776 m a.s.l. and span an altitudinal range of 290–1500 m a.s.l. Indications of ice flow suggest at least 105 active glaciers in the Torngat Mountains. Analysis of morphometric and topographic parameters suggests that the regional distribution of ice masses is linked to physiographic setting while the preservation of coastal ice masses at low elevation is related to local meteorological conditions. In the most coastal environments, ice masses are shown to exist below the regional glaciation level due to topographic shadowing, coastal proximity and widespread debris cover. This study provides a baseline for future change assessment.

We assess the runoff and surface mass balance (SMB) of the Greenland ice sheet in the Nuuk region (southwest) using output of two regional climate models (RCMs) evaluated by observations. The region encompasses six glaciers that drain into Godthåbsfjord. RCM data (1960–2012) are resampled to a high spatial resolution to include the narrow (relative to the native grid spacing) glacier trunks in the ice mask. Comparing RCM gridded results with automatic weather station (AWS) point measurements reveals that locally models can underestimate ablation and overestimate accumulation by up to tens of per cent. However, comparison with lake discharge indicates that modelled regional runoff totals are more accurate. Model results show that melt and runoff in the Nuuk region have doubled over the past two decades. Regional SMB attained negative values in recent high-melt years. Taking into account frontal ablation of the marine-terminating glaciers, the region lost 10–20 km3 w.e. a–1 in 2010–12. If 2010 melting prevails during the remainder of this century, a low-end estimate of sea-level rise of 5 mm is expected by 2100 from this relatively small section (2.6%) of the ice sheet alone.

Field observations, aerial photographs and satellite images are used to reconstruct the past surges of Lowell Glacier, Yukon, Canada, since 1948 based on the timing of terminus advances. A total of five surges occurred over this time, each with a duration of ~1–2 years. The time between successive surges ranged from 12 to 20 years, and appears to have been shortening over time. The relatively short advance and quiescent phases of Lowell Glacier, together with rapid increases in velocity during surges, suggest that the surging is controlled by a hydrological switch. The 2009–10 surge saw ablation area velocities increase by up to two orders of magnitude from quiescent velocities, and the terminus increase in area by 5.1 km2 and in length by up to 2.85 km. This change in area was the smallest since 1948, and follows the trend of decreasing surge extents over time. This decrease is likely driven by a strongly negative surface mass balance of Lowell Glacier since at least the 1970s, and means that the current town site of Haines Junction is very unlikely to be flooded by damming caused by any future advances of the glacier under the current climate regime.

Ground-penetrating radar (GPR) is becoming a commonly applied technique in geomorphology. However, its use in the study of subglacial bedforms has yet to be fully explored and exploited. This paper presents the results of a GPR feasibility study conducted on a drumlinized terrain in Cumbria, UK, where five drumlins were investigated using multiple radar antenna frequencies. The site was selected for the presence of nearby bedrock outcrops, suggesting a shallow drumlinized diamict–bedrock contact and a permeable lithology. Despite the clayey sediment and unfavourable weather conditions, a considerable penetration depth of ~12 m was achieved when using a 50 MHz antenna, with a separation of 1 m, trace spacing of 1 m and 128-fold vertical stack. Results indicate that the drumlinized diamict is in direct erosional contact with the bedrock. While the internal drumlin geometry is generally chaotic on the stoss side, evidence of layering dipping downflow at an angle greater than the drumlin surface profile was found on the lee side. The inter-drumlin areas comprise ~4 m of infill sediment that masks part of the original drumlin profile. Overall, this study indicates that GPR can be deployed successfully in the study of glacial bedform sedimentary architecture.

Borehole instrument records from a cirque glacier with an overdeepened bed are examined to assess the validity of widely held glacial hydrological assumptions. At this glacier, hydraulic-potential calculations suggest water below overburden pressure will flow into the overdeepening, where the steepness of the riegel causes water to pool in the basin and increase in pressure. Our subglacial water pressure data also show high consistent pressures in the overdeepening and the presence of an active, variable-pressure drainage system towards the margin of the cirque. Therefore, we find that although uniform hydraulic-potential calculations are not directly applicable, they can still be useful for interpretation of the subglacial hydrological system. We also examine supercooling assumptions under different pressure and temperature regimes for water flowing over a riegel, driven using our borehole records of subglacial water temperatures that are consistently above the pressure-melting point during the late melt season. Our results show that even a slight increase in basal temperatures relative to the local pressure-melting point is sufficient to prevent a reduction in basal hydraulic conductivity as a result of supercooling freeze-on.

Knowledge of the surface topography, velocity field and strain field at an ice-core site is critical to the accurate interpretation of ice-core records. At Dome Argus (Dome A), where a Chinese deep ice-core drilling project is being carried out, we have produced an accurate surface digital elevation model from GPS measurements in January 2013 at 47 sites. We identify two peaks at Dome A, with the northern peak ~7 cm higher than the southern peak. Repeat GPS measurements at 12 sites in 2008 and 2013 provide a surface velocity field around the dome. The surface velocity ranges from 3.1±2.6 to 29.4±1.2 cm a–1, with a mean of 11.1 ~2.4 cm a–1. The surface flow directions are near perpendicular to the surface elevation contours. Velocities from GPS are lower than derived from satellite radar interferometry (InSAR). From GPS velocities, the accuracy of velocity from the existing InSAR velocity field is determined, resulting in a standard deviation of 0.570 m a–1 in speed and 117.5º in direction. This result is consistent with the reported accuracy of InSAR, showing the value of in situ GPS measurements for assessing and correcting remote-sensing results. A surface strain field for the drilling site over Dome A is calculated from 24 strain triangles, showing north–south extension, east– west compression and vertical layer thinning.

A long-standing problem in avalanche science is to understand how forests stop avalanches. In this paper we quantify the effect of forests on small and medium avalanches, crucial for road and skirun safety. We performed field studies on seven avalanches where trees affected the runout. We gathered information concerning the release zone location and dimension, deposition patterns and heights, runout distance and forest structure. In these studies the trees were not destroyed, but acted as rigid obstacles. Wedge-like depositions formed behind (1) individual tree stems, (2) dense tree groups and (3) young trees with low-lying branches. Using the observations as a guide, we developed a one-parameter function to extract momentum corresponding to the stopped mass from the avalanche. The function was implemented in a depth-averaged avalanche dynamics model and used to predict the observed runout distances and mean deposition heights for the seven case studies. The approach differs from existing forest interaction models, which modify avalanche friction to account for tree breakage and debris entrainment. Our results underscore the importance of forests in mitigating the danger from small-to-medium avalanches.

A method has been devised and tested for measuring the c-axis orientation of crystal grains in thin sections of glacier ice. The crystal orientation and grain size of ice are of great interest to glaciologists since these parameters contain information on the prior thermal and flow history of the ice. The traditional method of determining c-axis orientation involves a transmission measurement through an ice sample, a process that is time-consuming and therefore impractical for obtaining a continuous record. A reflection- or backscatter-based method could potentially be used inside boreholes, with bubbles as reflectors to avoid such drawbacks. The concept demonstration of this paper is performed on ice slices, enabling a direct comparison of accuracy with traditional methods. Measurements of the crystal orientations (θ, ϕ) in 11 grains showed an average error of ±0.8° in ϕ, with no grain error >1.4°. Measurements of θ showed an average error of ±8.2° on ten grains, with unexplained disagreement on the remaining grain. Although the technique is applied specifically to glacier ice, it should be generally applicable to any transparent birefringent polycrystalline material.

Several different methodologies have previously been employed in the tracking of glacier terminus change, though a systematic comparison of these has not been undertaken. The frequent application of single methods to multiple glaciers over large geographical areas such as Greenland, raises the question of whether individual methodologies are robust. In this study we evaluate three existing methodologies that have been widely used to track terminus change (the centre-line, bow and box methods) against a full range of idealized glaciological scenarios and six examples of real glaciers. We also evaluate two new methodologies that aim to reduce measurement error compared with the existing methodologies. The first is a modification to the box method that can account for termini retreating through fjords that change orientation (termed the curvilinear box method), while the second determines the average terminus position relative to the glacier centre line using an inverse distance weighting extrapolation (termed the extrapolated centre-line method). No single method tested achieved complete accuracy for all scenarios, though the extrapolated centre-line method was able to successfully account for variable fjord orientation, width and terminus geometry with the least error.

Glacier roughness at sub-metre scales is an important control on the ice surface energy balance and has implications for scattering energy measured by remote-sensing instruments. Ice surface roughness is dynamic as a consequence of spatial and temporal variation in ablation. To date, studies relying on singular and/or spatially discrete two-dimensional profiles to describe ice surface roughness have failed to resolve common patterns or causes of variation in glacier surface morphology. Here we demonstrate the potential of close-range digital photogrammetry as a rapid and cost-effective method to retrieve three-dimensional data detailing plot-scale supraglacial topography. The photogrammetric approach here employed a calibrated, consumer-grade 5 Mpix digital camera repeatedly imaging a plot-scale (≤25 m2) ice surface area on Midtre Lovénbreen, Svalbard. From stereo-pair images, digital surface models (DSMs) with sub-centimetre horizontal resolution and 3 mm vertical precision were achieved at plot scales ≤4 m2. Extraction of roughness metrics including estimates of aerodynamic roughness length (z0) was readily achievable, and temporal variations in the glacier surface topography were captured. Close-range photogrammetry, with appropriate camera calibration and image acquisition geometry, is shown to be a robust method to record sub-centimetre variations in ablating ice topography. While the DSM plot area may be limited through use of stereo-pair images and issues of obliquity, emerging photogrammetric packages are likely to overcome such limitations.

The Randolph Glacier Inventory (RGI) is a globally complete collection of digital outlines of glaciers, excluding the ice sheets, developed to meet the needs of the Fifth Assessment of the Intergovernmental Panel on Climate Change for estimates of past and future mass balance. The RGI was created with limited resources in a short period. Priority was given to completeness of coverage, but a limited, uniform set of attributes is attached to each of the ~198 000 glaciers in its latest version, 3.2. Satellite imagery from 1999–2010 provided most of the outlines. Their total extent is estimated as 726 800 ± 34 000 km2. The uncertainty, about ±5%, is derived from careful single-glacier and basin-scale uncertainty estimates and comparisons with inventories that were not sources for the RGI. The main contributors to uncertainty are probably misinterpretation of seasonal snow cover and debris cover. These errors appear not to be normally distributed, and quantifying them reliably is an unsolved problem. Combined with digital elevation models, the RGI glacier outlines yield hypsometries that can be combined with atmospheric data or model outputs for analysis of the impacts of climatic change on glaciers. The RGI has already proved its value in the generation of significantly improved aggregate estimates of glacier mass changes and total volume, and thus actual and potential contributions to sea-level rise.

Tidewater glaciers in Greenland experienced widespread retreat during the last century. Information on their behaviour prior to this is often poorly constrained due to lack of observations, while determining the drivers prior to instrumental records is also problematic. Here we present a record of the dynamics of Kangiata Nunaata Sermia (KNS), southwest Greenland, from its Little Ice Age maximum (LIAmax) to 1859 – the period before continuous air temperature observations began at Nuuk in 1866. Using glacial geomorphology, historical accounts, photographs and GIS analyses, we provide evidence KNS was at its LIAmax by 1761, had retreated by ~5 km by 1808 and a further 7 km by 1859. This predates retreat at Jakobshavn Isbræ by 43–113 years, demonstrating the asynchroneity of tidewater glacier terminus response following the LIA. We use a one-dimensional flowband model to determine the relative sensitivity of KNS to atmospheric and oceanic climate forcing. Results demonstrate that terminus forcing rather than surface mass balance drove the retreat. Modelled glacier sensitivity to submarine melt rates is also insufficient to explain the retreat observed. However, moderate increases in crevasse water depth, driving an increase in calving, are capable of causing terminus retreat of the observed magnitude and timing.

The potential for future dynamical instability of Pine Island Glacier, West Antarctica, has been addressed in a number of studies, but information on its past remains limited. In this study we use airborne radio-echo sounding (RES) data acquired over Pine Island Glacier to investigate past variations in accumulation pattern. In the dataset a distinctive pattern of layers was identified in the central part of the glacier basin. We use these layers as chronological identifiers in order to construct elevation maps of the internal stratigraphy. The observed internal layer stratigraphy is then compared to calculated stratigraphy from a three-dimensional ice-flow model that has been forced with different accumulation scenarios. The model results indicate that the accumulation pattern is likely to have changed at least twice since the deposition of the deepest identified layer. Additional RES data linked to the Byrd ice core provide an approximate timescale. This timescale suggests that the layers were deposited at the beginning of or during the Holocene period. Thus the widespread changes occurring in the coastal extent of the West Antarctic ice sheet at the end of the last glacial period could have been accompanied by changes in accumulation pattern.

We provide a high-resolution map of elevation change rates at the Juneau Icefield (JIF), southeastern Alaska, in order to quantify its contribution to sea-level rise between 2000 and 2009/2013. We also produce the first high-resolution map of ice speeds at the JIF, which we use to constrain flux and look for acceleration. We calculate using stacked digital elevation models (DEMs) from the Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER) instrument and the Shuttle Radar Topography Mission (SRTM), taking into account SRTM C-band penetration via comparison with SRTM X-band elevations. Overall, the JIF is losing mass less rapidly (0.13 ± 0.12 m w.e. a–1) than other Alaskan icefields (0.79 m w.e. a–1). We determine glacier speeds using pixel-tracking on optical image pairs acquired from 2001 to 2010 by ASTER, from radar image pairs acquired between 2007 and 2011 and from radar interferometry in 1995. We detect seasonal speed variations but no interannual acceleration, ruling out dynamics as the cause of the observed thinning. Thinning must therefore be due to the documented warming in the region. Flux measurements confirm this for Mendenhall Glacier, showing that calving constitutes only 2.5–5% of mass loss there.

In polar ice sheets, the average grain size varies with depth. Ice grain size increases due to several factors including ice temperature and impurity content, which in turn varies with climate. The effect of impurities on grain growth is thought to be crucial but has never been observed experimentally. Using a methodology recently developed at Royal Holloway University of London, in situ chemical analysis of frozen ice at sub-ppm concentrations with unprecedented spatial resolution (~150 μm) is achievable using ultraviolet laser ablation inductively coupled plasma mass spectrometry (UV-LA-ICPMS) featuring a two-volume cryo-LA-cell. Following surface cleaning with a custom-built vice equipped with a ceramic blade, NGRIP ice slabs (~86 ka before AD 2000) have been analysed using a series of one-dimensional profiles and two-dimensional maps of laser spots at a resolution of 200–300 μm. Results demonstrate that cation impurities are not uniformly distributed in ice layers and show significant variations in concentration on a sub-millimetre scale. Furthermore, a different pattern of elemental distribution between clear ice and layers enriched in impurities (cloudy bands) has been identified: while concentration differences for cloudy bands are not resolvable between boundaries and inner grain domains, within clear ice, grain boundaries and junctions are significantly (up to 100 times) impurity-enriched relative to corresponding grain interiors.

Flow at ice divides, their shape, size and internal structure depend not only on local conditions, but also on the flow regimes and past histories of the surrounding ice masses. Here we use field data from Kealey Ice Rise, Ellsworth Land, West Antarctica, in combination with flow modelling to investigate any possible signs of transients in the flow of the surrounding ice masses. Kealey Ice Rise shows linear surface features running parallel to its ridge in satellite imagery and a conspicuous layering in the ground-penetrating radar data known as double-peaked Raymond bumps. Through numerical modelling, by using an anisotropic full-Stokes thermomechanical flow solver, we analyse the evolution of Kealey Ice Rise and the timescales involved. We conclude that the features observed in the stratigraphy of Kealey Ice Rise require at least 3 ka of near-stationary flow conditions. However, we cannot exclude the possibility of a recent flow reorganization in the last century. We stress that the signs of stationary flow in radar stratigraphy observed in Kealey Ice Rise have been observed in other ice divides in the East Ellsworth Land area, suggesting stationary flow conditions over a millennial timescale in the region.

We investigate the performance of five glacier melt models over a multi-decadal period in order to assess their ability to model future glacier response. The models range from a simple degree-day model, based solely on air temperature, to more-sophisticated models, including the full shortwave radiation balance. In addition to the empirical models, the performance of a physically based energy-balance (EB) model is examined. The melt models are coupled to an accumulation and a surface evolution model and applied in a distributed manner to Rhonegletscher, Switzerland, over the period 1929–2012 at hourly resolution. For calibration, seasonal mass-balance measurements (2006–12) are used. Decadal ice volume changes for six periods in the years 1929–2012 serve for model validation. Over the period 2006–12, there are almost no differences in performance between the models, except for EB, which is less consistent with observations, likely due to lack of meteorological in situ data. However, simulations over the long term (1929–2012) reveal that models which include a separate term for shortwave radiation agree best with the observed ice volume changes, indicating that their melt relationships are robust in time and thus suitable for long-term modelling, in contrast to more empirical approaches that are oversensitive to temperature fluctuations.

This study evaluates the impact of a recent upgrade in the physics package of the regional atmospheric climate model RACMO2 on the simulated surface mass balance (SMB) of the Antarctic ice sheet. The modelled SMB increases, in particular over the grounded ice sheet of East Antarctica (+44 Gt a–1), with a small change in West Antarctica. This mainly results from an increase in precipitation, which is explained by changes in the cloud microphysics, including a new parameterization for ice cloud supersaturation, and changes in large-scale circulation patterns, which alter topographically forced precipitation. The spatial changes in SMB are evaluated using 3234 in situ SMB observations and ice-balance velocities, and the temporal variability using GRACE satellite retrievals. The in situ observations and balance velocities show a clear improvement of the spatial representation of the SMB in the interior of East Antarctica, which has become considerably wetter. No improvements are seen for West Antarctica and the coastal regions. A comparison of model SMB temporal variability with GRACE satellite retrievals shows no significant change in performance.

During summer 2011, seismic activity in the ablation zone of the western Greenland ice sheet (GrIS) was monitored using a network of three-component seismometers. The seismic record includes a large variety of icequakes and seismic tremors that demonstrate a clear correlation with subglacial water flow. We verified the existence of well-known shallow icequakes (related to surface crevasse formation), deep icequakes (located at 100–160 m depth) and narrow-banded short-term seismic tremors (tens of seconds in duration). In addition, we present previously unreported long-term tremors lasting several hours. Using attenuation of the measured tremor amplitude, we locate the epicentre of this long-term tremor to a large moulin within our study area. Between 3 and 11 Hz, our continuous seismic record is dominated by this ‘moulin tremor’ and shows strong correlation with the water level of the generating moulin. We argue that monitoring of icequake and glacial tremor sources bears high potential for investigating glacier hydraulics and dynamics, and is thus an ideal supplement to traditional glaciological measurements.