A new radioscopic imaging technique has been developed to measure firn density in unprecedented resolution and accuracy even when the porosity is low or the geometry of a core or piece of core is not perfect. The technique is based on an X-ray microfocus computer tomograph (ICE-CT) designed especially for ice-core applications. Applied on an archive piece of the Antarctic firn core B32 drilled in Dronning Maud Land in 1998, the obtained density profile shows a strong correlation with the calcium ion concentration as found previously in Greenland. Given the impurity–density relationship found previously in Greenland, our result suggests both improved accuracy of the new density measurements and an impurity–density relationship with a similar magnitude in Greenland to that on the Antarctic plateau. Our measurements provide first evidence that the impurity–density relationship is a universal feature of polar firn and that the calcium ion concentration can serve as a proxy to describe quantitatively the effect of the impurities on densification.

The generalized (Glen) flow relation for ice, involving the second invariants of the stress deviator and strain-rate tensors, is only expected to hold for isotropic polycrystalline ice. Previous single-stress experiments have shown that for the steady-state flow, which develops at large strains, the tertiary strain rate is greater than the minimum (secondary creep) value by an enhancement factor which is larger for shear than compression. Previous experiments combining shear with compression normal to the shear plane have shown that enhancement of the tertiary octahedral strain rate increases monotonically from compression alone to shear alone. Additional experiments and analyses presented here were conducted to further investigate how the separate tertiary shear and compression strain-rate components are related in combined stress situations. It is found that tertiary compression rates are more strongly influenced by the addition of shear than is given by a Glen-type flow relation, whereas shear is less influenced by additional compression. A scalar function formulation of the flow relation is proposed, which fits the tertiary creep data well and is readily adapted to a generalized form that can be extended to other stress configurations and applied in ice mass modelling.

The rate of mass loss from the Greenland ice sheet has increased over the past decade due, in large part, to changes in marine-terminating outlet glacier dynamics. These changes are attributed to increased submarine melt rates of floating ice tongues and submerged calving faces resulting from increased coastal ocean heat transport. We use remotely sensed data to calculate submarine melt rates for 13 marine-terminating outlet glaciers in Greenland on a semi-annual basis between 2000 and 2010. We assess temporal and spatial variability in the calculated submarine melt rates and compare those variabilities to concurrent glacier change and offshore ocean temperatures. Over the period of study, average melt rates ranged from 0.03 to 2.98 m d−1 and account for 5–85% of the total volume loss from the floating ice tongue, with no clear spatial pattern. Only four glaciers show substantial interannual variability in melt rate during the decade. Melt rates were uncorrelated with front retreat, speed and changes in ocean temperature. Although the small sample size limits our analysis of the relationship between oceanographic forcing and glacier response, these data suggest that the calving rate may vary with discharge but that submarine melt rates are independent of grounding line discharge.

Levoglucosan is a unique marker for biomass burning that can be transported in the atmosphere and preserved in archives such as ice cores. A new method to determine the concentrations of levoglucosan in Tibetan ice-core samples using high-performance liquid chromatography with electrospray ionization mass spectrometry (HPLC-ESI/MS) was developed. Levoglucosan was separated from coeluting water-soluble organic compounds using a C18 column with a gradient program from 50% to 90% methanol in ultrapure water. An external standard calibration curve (R2 = 0.9958) was established by plotting the ion m/z 163 [M+H]+ peak area versus the amount of analyte. The repeatability ranges between 11% and 2% at a concentration around 10 and 150 ng mL−1. The limit of detection was 10 ng mL−1 and the limit of quantification was 40 ng mL−1. Levoglucosan concentrations ranged from 10 to 718 ng mL−1 in the Muztagh Ata ice core and from 10 to 93 ng mL−1 in the Tanggula ice core. These concentrations, up to 1000 times higher than those measured in samples from Antarctic and Greenland, showed the higher vulnerability of the Tibetan Plateau glaciers to biomass burning events.

Supraglacial lakes on grounded regions of the Greenland and Antarctic ice sheets sometimes produce ‘lake ogives’ or banded structures that sweep downstream from the lakes. Using a variety of remote-sensing data, we demonstrate that lake ogives originate from supraglacial lakes that form each year in the same bedrock-fixed location near the equilibrium-line altitude. As the ice flows underneath one of these lakes, an ‘image’ of the lake is imprinted on the ice surface both by summer-season ablation and by superimposed ice (lake ice) formation. Ogives associated with a lake are sequenced in time, with the downstream ogives being the oldest, and with spatial separation equal to the local annual ice displacement. In addition, lake ogives can have decimeter- to meter-scale topographic relief, much like wave ogives that form below icefalls on alpine glaciers. Our observations highlight the fact that lake ogives, and other related surface features, are a consequence of hydrological processes in a bedrock-fixed reference frame. These features should arise naturally from physically based thermodynamic models of supraglacial water transport, and thus they may serve as fiducial features that help to test the performance of such models.

Glacier-dammed lakes can yield subglacial outburst floods (jökulhlaups) repeatedly. Predicting flood timing is crucial for hazard mitigation, but incomplete understanding of flood-initiation physics makes this challenging. Here we examine the predictability of the timing of jökulhlaups from Merzbacher Lake, Kyrgyzstan, using five flood-date prediction models of varying complexity. The simplest model, which offers a benchmark against which the other models are compared, assumes that floods occur on the same date each year. The other four models predict flood dates using a flood-initiation threshold approach and incorporate weather forcing (approximated by the output of two climate reanalyses) behind the meltwater input to the lake; the most complex of these models accounts for a moving subglacial water divide beneath the glacier that dams the lake. Each model is optimized against recorded flood dates to maximize its prediction ability. In terms of their flood prediction ability, our two best models are those that assume a variable outburst threshold governed by the rate of meltwater input to the lake and the rate of lake-level rise. They excel over the simplest and most complex models and correctly predict flood dates to within ±20 days 57.4% of the time. We also quantify the impact of weather uncertainty on prediction success. Our findings can inform practical flood-forecasting schemes and future investigations of flood-initiation physics.

Sea ice is generally covered with snow, which can vary in thickness from a few centimeters to >1 m. Snow cover acts as a thermal insulator modulating the heat exchange between the ocean and the atmosphere, and it impacts sea-ice growth rates and overall thickness, a key indicator of climate change in polar regions. Snow depth is required to estimate sea-ice thickness using freeboard measurements made with satellite altimeters. The snow cover also acts as a mechanical load that depresses ice freeboard (snow and ice above sea level). Freeboard depression can result in flooding of the snow/ice interface and the formation of a thick slush layer, particularly in the Antarctic sea-ice cover. The Center for Remote Sensing of Ice Sheets (CReSIS) has developed an ultra-wideband, microwave radar capable of operation on long-endurance aircraft to characterize the thickness of snow over sea ice. The low-power, 100 mW signal is swept from 2 to 8 GHz allowing the air/snow and snow/ ice interfaces to be mapped with 5 cm range resolution in snow; this is an improvement over the original system that worked from 2 to 6.5 GHz. From 2009 to 2012, CReSIS successfully operated the radar on the NASA P-3B and DC-8 aircraft to collect data on snow-covered sea ice in the Arctic and Antarctic for NASA Operation IceBridge. The radar was found capable of snow depth retrievals ranging from 10 cm to >1 m. We also demonstrated that this radar can be used to map near-surface internal layers in polar firn with fine range resolution. Here we describe the instrument design, characteristics and performance of the radar.

Detailed maps of bed elevation and ice thickness are essential for understanding and projecting the evolution of the ice sheets. Such maps are traditionally obtained using airborne radar-sounding profiler data interpolated onto regular grids using geostatistical tools such as kriging. Here we compare three mapping techniques applied to a dense radar survey of Russell Glacier, West Greenland, by NASA Operation IceBridge: (1) radar tomography (RT) processing of the radar data to map the bed elevation, (2) interpolation of radar-derived thickness by ordinary kriging (KR) and (3) reconstruction of ice thickness based on the principles of mass conservation (MC) combining radar-sounding profiler and ice motion data. RT eliminates ambiguities caused by off-nadir reflections, but is spatially limited. KR yields a standard error in bed elevation of 35 m, but large errors (>300 m a−1) in flux divergence when combined with ice motion data. MC yields a comparable performance in bed elevation mapping, and errors smaller than 1 m a−1 in flux divergence. When the number of radar-sounding tracks is reduced, the performance of KR decreases more rapidly than for MC. Our study site shows that MC is capable of maintaining precision levels of 60 m at 400 m posting with flight tracks separated by 5 km.

We report sea-ice temperature and bulk salinity measurements as well as textural analysis from 33 first-year drift- and fast-ice stations sampled between November 2007 and June 2008 in the southern Beaufort Sea–Amundsen Gulf, Canadian Arctic, during the International Polar Year Circumpolar Flaw Lead (IPY-CFL) system study. We use this significant dataset to investigate the halo-thermodynamic evolution of sea ice from growth to melt. A strong desalination phase is observed over a small time window in the spring. Using calculated proxies of sea-ice permeability (brine volume fraction) and of the intensity of brine convection (Rayleigh number) we demonstrate that this phase corresponds to full-depth gravity drainage initiated by a restored connectivity of the brine network with warming in the spring. Most stations had a textural sequence typical of Arctic first-year ice, with granular ice overlying columnar ice. Unusual textural features were observed sporadically: sandwiched granular ice, platelet ice and draped platelet ice. We suggest that turbulence in leads and double diffusion in strong brine plumes following the refreeze of cracks are plausible mechanisms for the formation of these textures.

The contribution to sea level to 2200 from the grounded, mainland Antarctic Peninsula ice sheet (APIS) was calculated using an ice-sheet model initialized with a new technique computing ice fluxes based on observed surface velocities, altimetry and surface mass balance, and computing volume response using a linearized method. Volume change estimates of the APIS resulting from surface mass-balance anomalies calculated by the regional model RACMO2, forced by A1B and E1 scenarios of the global models ECHAM5 and HadCM3, predicted net negative sea-level contributions between −0.5 and −12 mm sea-level equivalent (SLE) by 2200. Increased glacier flow due to ice thickening returned ∼15% of the increased accumulation to the sea by 2100 and ∼30% by 2200. The likely change in volume of the APIS by 2200 in response to imposed 10 and 20 km retreats of the grounding line at individual large outlet glaciers in Palmer Land, southern Antarctic Peninsula, ranged between 0.5 and 3.5 mm SLE per drainage basin. Ensemble calculations of APIS volume change resulting from imposed grounding-line retreat due to ice-shelf break-up scenarios applied to all 20 of the largest drainage basins in Palmer Land (covering ∼40% of the total area of APIS) resulted in net sea-level contributions of 7–16 mm SLE by 2100, and 10–25 mm SLE by 2200. Inclusion of basins in the northern peninsula and realistic simulation of grounding-line movement for AP outlet glaciers will improve future projections.

Ice-core data are the only source of ice response on timescales of hundreds to thousands of years, which is necessary to describe polar ice-sheet flows and their interaction with changing climatic conditions. Ice properties, such as a grain growth relation, require knowledge of each ice particle’s history since first deposited at the surface, and not simply the present core state. This study is an analysis of 11 ‘present’ ice cores, using their age/depth data and present accumulation, to infer their past accumulation and mean strain rate and, in turn, the strain at each level of the core. These are the ingredients necessary to determine the particle paths and associated evolution of properties. The conventional assumptions of a purely vertical velocity distribution with a uniform compressive strain rate in the core are made, but unlike an earlier analysis which allowed independent surface accumulation and strain-rate histories, it is now assumed that the strain rate is proportional to the accumulation. Determination of the past accumulation and proportionality factor is accomplished by correlation with the present age/depth data, and while good correlations are obtained – a consistent solution – other close correlations can yield quite different accumulation histories and corresponding proportional strain-rate histories.

We have determined the ice mass evolution of the Antarctic and Greenland ice sheets (AIS and GIS) and Gulf of Alaska (GOA) glaciers from a new GRACE global solution of equal-area surface mass concentration parcels (mascons) in equivalent height of water. The mascons were estimated directly from the reduction of the inter-satellite K-band range-rate (KBRR) observations, taking into account the full noise covariance, and formally iterating the solution. The new solution increases signal recovery while reducing the GRACE KBRR observation residuals. The mascons were estimated with 10 day and 1 arcdeg equal-area sampling, applying anisotropic constraints. An ensemble empirical mode decomposition adaptive filter was applied to the mascon time series to compute annual mass balances. The details and causes of the spatial and temporal variability of the land-ice regions studied are discussed. The estimated mass trend over the total GIS, AIS and GOA glaciers for the time period 1 December 2003 to 1 December 2010 is −380 ± 31 Gt a−1, equivalent to −1.05 ± 0.09 mm a−1 sea-level rise. Over the same time period we estimate the mass acceleration to be −41 ± 27 Gt a−2 , equivalent to a −0.11 ± 0.08 mm a−2 sea-level acceleration. The trends and accelerations are dependent on significant seasonal and annual balance anomalies.

As part of the Ross Ice Shelf Geophysical and Glaciological Survey (RIGGS), ice velocities were measured on the Ross Ice Shelf (RIS) during 1973–78. Comparisons of these with velocity estimates at the same locations derived from RADARSAT synthetic aperture radar (SAR) measurements in 1997 and 2009 show velocity reduction in the southeast quadrant of the ice shelf by almost 200 m a−1, with deceleration rates increasing with time. Large areas of ice shelf in this region are lightly grounded, forming an ‘ice plain’ that increases local buttressing of the ice streams. ICESat measurements show this ice plain to be thickening. The observed decrease in ice-shelf velocities implies a total reduction in the mass of ice flowing into the RIS from the West Antarctic ice sheet (WAIS) by ~23 Gt a−1, shifting the mass balance of the WAIS drainage basin from strongly negative in the 1970s to strongly positive in 2009. The resulting decrease in ice advection should lead to ice-shelf thinning further seaward of the ice plain. This thinning would reduce the lateral drag and back-stress of the shelf ice, further contributing to thinning through an increase in spreading rate. ICESat measurements show recent thinning of most of the freely floating ice shelf.

The Glacier Bay region of southeast Alaska, USA, and British Columbia, Canada, has undergone major glacier retreat since the Little Ice Age (LIA). We used airborne laser altimetry elevation data acquired between 1995 and 2011 to estimate the mass loss of the Glacier Bay region over four time periods (1995–2000, 2000–05, 2005–09, 2009–11). For each glacier, we extrapolated from center-line profiles to the entire glacier to estimate glacier-wide mass balance, and then averaged these results over the entire region using three difference methods (normalized elevation, area-weighted method and simple average). We found that there was large interannual variability of the mass loss since 1995 compared with the long-term (post-LIA) average. For the full period (1995–2011) the average mass loss was 3.93 ± 0.89 Gt a−1 (0.6 ± 0.1 m w.e. a−1), compared with 17.8 Gt a−1 for the post-LIA (1770–1948) rate. Our mass loss rate is consistent with GRACE gravity signal changes for the 2003–10 period. Our results also show that there is a lower bias due to center-line profiling than was previously found by a digital elevation model difference method.

Consumer-grade digital cameras have become ubiquitous accessories of science. Particularly in glaciology, the recognized importance of short-term variability has motivated their deployment for increasingly time-critical observations. However, such devices were never intended for precise timekeeping, and their use as such needs to be accompanied by appropriate management of systematic, rounding and random errors in reported image times. This study describes clock drift, subsecond reporting resolution and timestamp precision as the major obstacles to precise camera timekeeping, and documents the subsecond capability of camera models from 17 leading manufacturers. We present a complete and accessible methodology to calibrate cameras for absolute timing and provide a suite of supporting scripts. Two glaciological case studies serve to illustrate how the methods relate to contemporary investigations: (1) georeferencing aerial photogrammetric surveys with camera positions time-interpolated from GPS tracklogs; and (2) coupling videos of glacier-calving events to synchronous seismic waveforms.

Observations show that the Greenland ice sheet has been losing mass at an increasing rate over the past few decades, which makes it a major contributor to sea-level rise. Here we use a three-dimensional higher-order ice-flow model, adaptive mesh refinement and inverse methods to accurately reproduce the present-day ice flow of the Greenland ice sheet. We investigate the effect of the ice thermal regime on (1) basal sliding inversion and (2) projections over the next 100 years. We show that steady-state temperatures based on present-day conditions allow a reasonable representation of the thermal regime and that both basal conditions and century-scale projections are weakly sensitive to small changes in the initial temperature field, compared with changes in atmospheric conditions or basal sliding. We conclude that although more englacial temperature measurements should be acquired to validate the models, and a better estimation of geothermal heat flux is needed, it is reasonable to use steady-state temperature profiles for short-term projections, as external forcings remain the main drivers of the changes occurring in Greenland.

Inclinometry with embedded probes is analyzed with a Stokes model of a solid body floating in a fluid with much smaller viscosity for a two-dimensional flow field. The assumption that such a probe behaves like a Lagrangian unit vector is only justified for probes embedded in a Newtonian fluid with lengths at least four times their width. A fluid with Glen-type rheology results in a slightly smaller rotation rate of the probe compared to Newtonian fluids.

We explore the small-scale spatial and temporal transferability of model parameters between two points in the ablation zone of tropical Glaciar Shallap, Cordillera Blanca, Peru (9°S, −77° W; ~4800 m a.s.l.) in order to provide a robust assessment of the performance of a process-based glacier mass-balance model. Relative surface height change is calculated at hourly time-steps, and cumulative values are compared to surface height measurements made at irregular intervals (14–64 days) over the course of two continuous hydrological years (August 2006–August 2008). Best-performing parameter combinations were determined for each point from the outcome of 1000 model simulations for which parameters were varied randomly within a defined range. With these parameter combinations measurements for a specific location and time-span are well reproduced. Transferring the parameter combination as optimized for one location to the other location in the ablation zone increases the errors of modeled cumulative mass balance by 5–1326 mm ice eq.a−1. Transferring the parameter combinations as optimized for one year to the other year increases the modeled errors in cumulative mass balance by 18–3179 mm ice eq.a−1. Model errors generally increase during periods with frequent snowfall and snow cover. This could reflect either the inherent difficulty of modeling complex snow processes, or the inability of the model to correctly capture the pattern of albedo evolution at this site. The magnitude of errors associated with parameter transfer in space and time highlights the need for improving model performance for robust climatological and/or hydrological analyses on tropical glaciers.

Marine-terminating glaciers may experience seasonal and short-term flow variations, which can impact rates of ice flux through the glacier terminus. We explore the relationship between variability in the flow of a large tidewater glacier (Belcher Glacier, Nunavut, Canada), the seasonal cycle of surface meltwater production and the rapid drainage of supraglacial lakes. We demonstrate a novel method for analyzing time-lapse photography to quantify lake area change rates (a proxy for net filling and drainage rates) and develop a typology of lake drainage styles. GPS records of ice motion reveal four flow acceleration events which can be linked to lake drainage events discovered in the time-lapse photography. These events are superimposed on a longer pattern of velocity variation that is linked to seasonal variation in surface melting. At the terminus of the glacier, the ice displacement associated with the lake drainage events constitutes ∼10% of the seasonally accelerated displacement or 0.4% of the total annual ice displacement (336 m a−1). While the immediate ice response to these individual perturbations may be small, these drainage events may enhance overall seasonal acceleration by opening and/or sustaining meltwater conduits to the glacier bed.