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.

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.

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.

Identification of the transient snowline (TSL) from high spatial resolution Landsat imagery on Lemon Creek Glacier (LCG), southeast Alaska, USA, and Mittivakkat Gletscher (MG), southeast Greenland, is used to determine snow ablation rates, the equilibrium-line altitude (ELA) and the accumulation-area ratio (AAR). The rate of rise of the TSL during the ablation season on a glacier where the balance gradient is known provides a measure of the snow ablation rate. On both LCG and MG, snow pits were completed in regions that the TSL subsequently transects. This further provides a direct measure of the snow ablation rates for a particular year. TSL observations from multiple dates during the ablation season from 1998 to 2011 at LCG and 1999 to 2012 at MG were used to explore the consistency of the TSL rise and snow ablation rate. On LCG and MG the satellite-derived mean TSL migration rates were 3.8 ± 0.6 and 9.4 ± 9.1 m d−1, respectively. The snow ablation rates were 0.028 ± 0.004 m w.e. d−1 for LCG and 0.051 ± 0.018 m w.e. d−1 for MG estimated by applying a TSL–mass-balance-gradient method, and 0.031 ± 0.004 and 0.047 ± 0.019 m w.e. d−1 by applying a snow-pit–satellite method, illustrating significant agreement between the two different approaches for both field sites. Also, satellite-derived ELA and AAR, and estimated net mass-balance (Ba) conditions were in agreement with observed ELA, AAR and Ba conditions for LCG and MG.

We document changes for Mittivakkat Gletscher, the peripheral glacier in Greenland with the longest field-based observed mass-balance and surface velocity time series. Between 1986 and 2011, this glacier changed by −15% in mean ice thickness and −30% in volume. We attribute these changes to summer warming and lower winter snow accumulation. Vertical strain compensated for ∼60% of the elevation change due to surface mass balance (SMB) in the lower part, and ∼25% in the upper part. The annual mean ice surface velocity changed by −30%, which can be fully explained by the dynamic effect of ice thinning, within uncertainty. Mittivakkat Gletscher summer surface velocities were on average 50–60% above winter background values, and up to 160% higher during peak velocity events. Peak velocity events were accompanied by uplift of a few centimeters.

We determined length, area, elevation and volume change of the Columbia Icefield using Interprovincial Boundary Commission Survey maps from 1919, eight sets of aerial photographs from 1948 to 1993, and satellite data from 1999 to 2009. Over the period 1919–2009, glaciers on average retreated 1150 ± 34 m and shrank by 2.4 ± 0.2 km2. Total area loss was 59.6 ± 1.2 km2 (23 ± 5%), and mean elevation change was −49 ± 25 m w.e., resulting in a total volume loss of 14.3 ± 2.0 km3 w.e. Large outlet glaciers experienced the greatest absolute ice loss, while small, detached glaciers lost the most relative length and area. Thinning rates of debris-covered ice were 30–60% lower than those for clean ice. All glacier changes were significantly correlated with each other (p < 0.01), with r values ranging from 0.54 to 0.82. Temperature is correlated with length and area change over periods lagged 1–5 years (p < 0.05), and with elevation and volume change over periods lagged 9–18 years (p < 0.05). Precipitation is correlated with glacier change over periods lagged 1–10 years (p < 0.05).

We present satellite-derived velocity patterns for the two contrasting melt seasons of 2009–10 across Russell Glacier catchment, a western, land-terminating sector of the Greenland ice sheet which encompasses the K(angerlussuaq)-transect. Results highlight great spatial heterogeneity in flow, indicating that structural controls such as bedrock geometry govern ice discharge into individual outlet troughs. Results also reveal strong seasonal flow variability extending 57 km up-glacier to 1200 m elevation, with the largest acceleration (100% over 11 days) occurring within 10 km of the margin coincident with spring melt. By late July 2010, 2 weeks before peak melt and runoff, 48 % of the 2400 km2 catchment had slowed to less than the winter mean. This observation supports the hypothesis that the subglacial hydrological system evolves from an inefficient distributed to an efficient drainage system, regulating flow dynamics. Despite this, the cumulative surface flux over the record melt year of 2010 was still greater compared with the perturbation over the average melt year of 2009. This study supports the proposition that local surface meltwater runoff couples to basal hydrology driving ice-sheet dynamics, and although the effect is nonlinear, our observations indicate that greater meltwater runoff yields increased net flux over this sector of the ice sheet.

The dependence of oxygen isotope fractionation on ice growth rate during the freezing of sea water is investigated based on laboratory experiments and field observations in McMurdo Sound, Antarctica. The laboratory experiments were performed in a tank filled with sea water, with sea ice grown under calm conditions at various room temperatures ranging from −5°C to −20°C. In McMurdo Sound, the ice growth rate was monitored using thermistor probes for first-year landfast ice that grew to ∼2 m in thickness. Combining these datasets allows, for the first time, examination of fractionation at a wide range of growth rates from 0.8 × 10−7 to 9.3 × 10−7 m s−1. In the analysis a stagnant boundary-layer model is parameterized using these two independent datasets. As a result, the optimum values of equilibrium pure-ice fractionation factor and boundary-layer thickness are estimated. It is suggested that a regime shift may occur at a growth rate of ∼2.0 × 10−7 m s−1. A case study on sea ice in the Sea of Okhotsk, where the growth rate is modeled by coupling the thermodynamic properties of the sea ice with meteorological data, demonstrates the utility of the fitted models.

The results of measuring relative brine content (Qlm) in the liquid phase of frozen sea water at temperatures between −2°C and −43°C with variable rate and direction of temperature change are presented. Various pulsed and steady-state nuclear magnetic resonance (NMR) methods are used and the results are compared with data in the literature. Temperature dependences of Qlm in sea water with dry NaCl added are given. Quasi-hysteresis corresponding to the crystallization range of certain seawater salts, mainly NaCl, which partially precipitates as the crystalline hydrate NaCl·2H2O at temperatures below −23°C, is detected in the temperature dependence of Qlm. The reasons why Qlm values, under various experimental conditions, differ from data in the literature are explained. Brine salinity (S) is calculated from the relative intensity of NMR signals of Ql, and agrees well with previously published data over a wide temperature range. It is shown that determining Qlm from base measurements of Ql in ice does not require information about brine salinity. Empirical equations for calculating brine content at different sea-water temperatures and salinity are given.

Results of determinations of proton magnetic relaxation time (T1) in the liquid phase of frozen sea water at temperatures between −2°C and −38°C, with variable rate and direction of temperature change, are presented. Temperature hysteresis corresponding to the crystallization range of certain sea-water salts, mainly NaCl which partially precipitates as the crystalline hydrate NaCl.°2H2O at temperatures below −23°C, is detected in the temperature dependence of T1, as in the analogous dependence of brine content. It is shown that formation of crystalline hydrates in brine is preceded by complete solvation of ions at −5 to −6°C. The main causes of hysteresis in temperature dependence of T1 and the extension of loop limits compared with the analogous loop for liquid brine content are indicated. Additional potentialities of using the nuclear magnetic resonance method for investigating sea-ice properties are presented.

The roughness of the bed beneath ice streams exerts an important control on their dynamics. Here the first in-depth analysis of roughness beneath the Greenland ice sheet is presented. Much of the interior is underlain by a smooth bed, while the margins are much rougher; this is particularly pronounced in the east where more mountainous topography dominates and fast-flow features are laterally constrained in deep narrow valleys. In contrast, fast flow in the west is much less laterally constrained and areas of high roughness are less extensive and pronounced. It is proposed that there is a major geological control on the distribution of bed variability since a significant thrust-fault coincides approximately with the boundary between rough terrain in the east and the smooth central region. Furthermore there is an abrupt change in roughness approximately coinciding with the crossing of this fault line. This suggests a limiting factor on the extent of fast flow in the east, which is lacking in the west. The size of many glaciers draining Greenland makes their local bed conditions difficult to determine with great confidence. However, the much larger Petermann Glacier lies in a deep trough characterized by a smooth bed that extends some distance into the ice sheet. The smooth bed may be due to deformable marine sediments which facilitate faster flow, although smoothing could also result from ice dynamics and subglacial erosion. That the smooth bed of Petermann Glacier extends some considerable distance raises concerns about the possible stability of this feature, and perhaps others in Greenland too.

Physically based projections of the Greenland ice sheet contribution to future sea-level change are subject to uncertainties of the atmospheric and oceanic climatic forcing and to the formulations within the ice flow model itself. Here a higher-order, three-dimensional thermomechanical ice flow model is used, initialized to the present-day geometry. The forcing comes from a high-resolution regional climate model and from a flowline model applied to four individual marine-terminated glaciers, and results are subsequently extended to the entire ice sheet. The experiments span the next 200 years and consider climate scenario SRES A1B. The surface mass-balance (SMB) scheme is taken either from a regional climate model or from a positive-degree-day (PDD) model using temperature and precipitation anomalies from the underlying climate models. Our model results show that outlet glacier dynamics only account for 6–18% of the sea-level contribution after 200 years, confirming earlier findings that stress the dominant effect of SMB changes. Furthermore, interaction between SMB and ice discharge limits the importance of outlet glacier dynamics with increasing atmospheric forcing. Forcing from the regional climate model produces a 14–31 % higher sea-level contribution compared to a PDD model run with the same parameters as for IPCC AR4.

Basal crevasses may play an important precursory role in determining both the location and propagation of rifts and iceberg dimensions. For example, icebergs calved recently from Thwaites Glacier, Antarctica, have the same width as surface undulations, strengthening the connection between basal crevasses, rifting and calving. We explore a novel method for estimating the heights of basal crevasses formed at the grounding lines of ice shelves and ice streams. We employ a thin-elastic beam (TEB) formulation and tensional yielding criterion to capture the physics of flexed ice at grounding lines. Observations of basal crevasse heights compare well with model predictions in the Siple Coast region of the Ross Ice Shelf. We find that the TEB method is most accurate in areas of low strain rate. We also test the method in other areas of Antarctica to produce order-of-magnitude maps of grounding-line basal crevasses and find general agreement with reported observations assuming basal crevasses develop in spatio-temporal sequence and are advected downstream. This method is computationally cheap and could be relatively easy to implement into damage-oriented large-scale ice models which aim at physically simulating calving and fracture processes.

Many of Greenland’s marine-terminating outlet glaciers have undergone rapid retreat in the past decade, accompanied by accelerated flow and dynamic thinning. Superimposed on this pattern of retreat, these glaciers undergo seasonal variations in terminus position, corresponding roughly to wintertime advance and summertime retreat. We compiled near-daily time series of terminus position for five of Greenland’s largest outlet glaciers (Daugaard Jensen, Kangerdlugssuaq and Helheim glaciers in East Greenland, and Jakobshavn Isbræ and Rink Isbræ in West Greenland) using Moderate Resolution Imaging Spectroradiometer (MODIS) satellite imagery. There are spatial differences in the timing of the onset of seasonal retreat among all the glaciers in our study, as well as variability in terminus behavior for individual glaciers from year to year. We examine whether this spatial and temporal variability is linked to above-freezing air temperatures or high sea surface temperatures, but find no simple relationship. Instead, we hypothesize that terminus geometry (ice thickness, subglacial topography, fjord bathymetry) exerts an important control on the response of marine-terminating glaciers to climate perturbations. Models for predicting outlet glacier response to climate change need to include this complex interaction between geometry and environmental forcing.

The values of the Voellmy friction parameters of 735 historical avalanches that have occurred along 26 paths in the Chamonix valley, France, since 1958 are back-analyzed with a depth-averaged hydraulic model, including sub-models for erosion, entrainment and deposition. For each path, the longitudinal and crosswise topographic profiles were derived from a high-resolution digital elevation model acquired by laser scanning. The initial snow depth and snow cohesion, as well as various physical properties of snow, were computed from numerical simulations of the detailed snowpack model Crocus fed by the SAFRAN meteorological analysis. For each event, the full ranges of the two friction parameters were scanned and the pairs of friction parameters for which the run-out altitude is found close enough to the observed one (with an uncertainty of ±5 m) were retained. Statistical class analysis was used to investigate the correlation between the obtained friction coefficients and the snow physical properties. No evident trend with the snow parameters was found for the inertial friction coefficient. For the static friction coefficient, an increasing trend with temperature and density was observed, as well as a decreasing trend with liquid water content and initial snow depth. Although modeling assumptions and limitations regarding data and the calibration procedure should be kept in mind, these trends are worth noting, allowing avalanche simulations to be refined to take into account prevailing weather and snow conditions.

El Niño and La Niña affect global climate and atmospheric circulation to determine winter temperature and precipitation patterns. Both winter temperatures and the associated precipitation patterns have effects on mountain snow deposition and snow avalanche occurrences. Approximately 25 000 slab avalanches from 30 winters were analyzed in relation to snowfall patterns contrasted for El Niño and La Niña winters for two avalanche areas with different snow climates in British Columbia (BC), Canada. La Niña winters were shown to produce more snow, more avalanches and a higher percentage of dry avalanches than wet avalanches. The data and analysis show that the avalanche patterns depend on the altitude and snow climate. Analysis of snowfall and accident data from the Andes of Chile suggests behavior opposite to BC. El Niño winters in central Chile produced the most snow and, by inference, the most avalanches. This paper is the first to show the links of El Niño and La Niña to snow avalanche activity.