Between 2003 and 2007, North Gasherbrum Glacier on the northeastern slope of the Karakoram mountains in Asia underwent a dramatic acceleration, during which a velocity wave propagated down the glacier. There was a significant transfer of ice from up-glacier downstream, which resulted in a strong surface elevation increase over the lower tongue, but only a moderate advance of the glacier snout. We interpret this behaviour as that of a glacier surge, and we explain the observations by means of a simple version of the Kamb drainage-switching theory.

A thin-blade snow hardness gauge was developed that measures penetration resistance over a length scale (on the order of 10–100 grain contacts) relevant to the fracture of slab avalanches. A thin blade was chosen to measure the ruptures of bonds and grain structures and minimize the effects of snow compaction during penetration. The apparatus consists of a 10 cm wide, 0.6 mm thick stainless-steel blade attached to a digital push–pull gauge. Blade penetration measurements are easy to conduct in the field and laboratory and required no post-processing or subjective interpretation. Measurements were conducted in snow pits to test the effects of penetration rate, blade orientation and blade width. The blade hardness index, defined as the maximum force of penetration, is a highly repeatable measure across observers compared to the hand hardness test. The blade hardness index was a better variable than the density for correlating with tensile strength measurements in a cold laboratory and with a cohesive strength measure in the field. As strength is one of the most important parameters in the fracture mechanics of slab avalanches, the strong correlation between thin-blade penetration and strength should benefit future slope stability evaluations using this gauge.

A new ring-shear device allows basal slip and related processes to be studied in laboratory experiments for the cases of hard or soft beds. The device rotates a confined ring of ice (0.9 m outside diameter) across a horizontal bed at a constant velocity or drag, while a vertical stress is applied and basal water pressure is controlled. A bath with circulating fluid regulated to ∼0.01°C surrounds the ice chamber and keeps the ice at its pressure-melting temperature. In a first experiment with a stepped rigid bed and zero basal water pressure, steady lengths of step cavities depended upon slip velocity raised to a power of 0.59, in general agreement with the square-root dependence of some models of sliding and linked-cavity hydraulics. Transient cavity growth after slip velocity increases was not monotonic, with damped volume oscillations that converged to a steady value. Once ice separated from lee surfaces, drag on the bed was constant and independent of slip velocity and cavity size, consistent with a shear-stress upper bound like that indicated by sliding models. Shear strains near the bed exceeded 30 and ice developed multiple-cluster c-axis fabrics similar to those of sheared ice in temperate glaciers.

Little quantitative information is available concerning the size of ice particles in the turbulent clouds of powder-snow avalanches. To quantify particle size distributions, we have developed an experimental device that collects particles in real-scale powder avalanches. The device was placed on the concrete bunker of the Swiss Vallée de la Sionne avalanche dynamics test site. On 31 January 2003, a large powder-snow avalanche struck the bunker and we were able to collect particle samples. The collected particles have been photographed and the pictures digitized. An image analysis tool allows us to determine an equivalent particle radius. The captured particles have a geometric mean of 0.16 mm; the largest particles were 0.8 mm in size and the smallest particles 0.03 mm.

We present a microstructural characterization of fine-grained layers from the top 90 m of firn from Summit, Greenland, performed using a combination of scanning electron microscopy techniques including secondary electron imaging, energy-dispersive spectroscopy and electron backscattered patterns, and X-ray microcomputed tomography. The impurities in the firn, both soluble impurities and dust particles, were found largely in the grain interiors. Both c- and a-axis pole figures do not show strong evidence of a preferred orientation of the grains even at the bottom of the firn column. The firn structure became increasingly anisotropic with vertical alignment in the top 3 m, probably due to vapor transport associated with dry-snow metamorphism. The anisotropy decreases below this level until at ∼50 m the average firn structure is close to isotropic. In the near surface, the level of anisotropy is weaker than at Hercules Dome, Antarctica, confirming that differences in accumulation rates and temperatures leave enduring evidence in the structure of the firn. The fraction of closed-off pores is relatively low until ∼65 m; below that it rises through the end of our sampling at 90 m. Our microstructure measurements on the microscale are consistent with in situ firn-air sampling measurements on a decimeter scale, both indicating the existence of the lock-in zone starting near 69 m depth, and pore close-off at 81 m at this site.

Creep tests provide invaluable data to better understand the physical properties of ice under various conditions. We describe here a simple, updated pneumatic apparatus for experimental studies of ice rheological properties. The apparatus is designed to perform two simultaneous compression creep tests either in a cold room or in atmospheric conditions when coupled to an external cooling circulator. We present results from calibration tests of the apparatus and provide new data from creep tests performed on temperate glacier ice samples. These calibration and creep results show that the apparatus is able to provide fast and reliable mechanical ice characterization. The secondary creep rates measured in this study range between 1.59 × 10−8 s−1 (at 0.21 MPa) and 4.38 × 10−7 s−1 (at 0.71 MPa) at −10°C for quasi-isotropic ice, which is consistent with former standard published data. The corresponding mean parameter, A, is 5.20 × 10−16 s−1 kPa−3, which also compares well with the range of reported other studies.

An approximation to the first-order momentum balance with consistent boundary conditions is derived using variational methods. Longitudinal and lateral stresses are treated as depth-independent, but vertical velocity gradients are accounted for both in the nonlinear viscosity and in the treatment of basal stress, allowing for flow over a frozen bed. A numerical scheme is presented that is significantly less computationally expensive than that of a fully three-dimensional (3-D) solver. The numerical solver is subjected to the ISMIP-HOM experiments and experiments involving nonlinear sliding laws, and results are compared with those of 3-D models. The agreement with first-order surface velocities is favorable down to length scales of 10 km for flow over a flat bed with periodic basal traction, and ∼40 km for flow over periodic basal topography.

SnowDrift3D, a high-resolution, atmospheric snow-transport model, is presented for the first time. In contrast to most state-of-the-art snowdrift models, atmospheric particle transport, i.e. saltation and suspension, is accounted for by one passive transport equation. The model uses unsteady wind fields (spatial resolution of up to 2 m) computed with an atmospheric computational fluid dynamics model that is directly connected to the numerical weather prediction model ALADIN. Sensitivity runs show that (1) the saltation mass flux is a function of cubic shear velocity, , (2) the model is marginally sensitive to the grid spacing at high resolutions (up to 2 m), (3) the model computes the redistribution of snow at high resolution in real time on dual core personal computers and (4) the changing topography of the snow cover should be included in cases of local erosion or deposition of a large amount of snow. Finally, we present a comparison of modeled and measured snow distributions obtained by terrestrial laser scanning showing area-wide linear correlation up to R = 0.33.

Debris cover is widely present in glacier ablation areas of the Tibetan Plateau, and its spatial distribution greatly affects glacier melt rates. High-resolution in situ measurements of debris thickness on Hailuogou glacier, Mount Gongga, southeastern Tibetan Plateau, show pronounced inhomogeneous debris distribution. An analysis of transverse and longitudinal profiles indicates that the ground-surveyed debris thicknesses and Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER)-derived thermal resistances of debris layers correlate strongly over the entire ablation area. Across- and along-glacier patterns of ASTER-derived thermal resistance correspond well with spatial patterns of debris thickness, which may reflect large-scale variations in the extent and thickness of the debris cover. The ice melt rate variability over the ablation area simulated by a surface energy-balance model that considered thermal resistance of the debris layer indicates clearly the crucial role of debris and its spatial continuity in modifying the spatial characteristics of melt rates. Because of the inhomogeneous distribution of debris thickness, about 67% of the ablation area on Hailuogou glacier has undergone accelerated melting, whereas about 19% of the ablation area has experienced inhibited melting, and the sub-debris melt rate equals the bare-ice melt rate in only 14% of the ablation area.

As ice streams flow into the Ross Ice Shelf, West Antarctica, their bed coupling transitions from weak to transient to zero as the ice goes afloat. Here we explore the nature of the bed across these crucial grounding zones using ice-penetrating radar. We collected several ground-based 2 MHz radar transects across the grounding zones of Whillans and Kamb Ice Streams and inferred bed-reflectivity changes from in situ measurements of depth-averaged dielectric attenuation, made possible by the observation of both primary and multiple bed echoes. We find no significant change in the bed reflectivity across either grounding zone. Combined with reflectivity modeling, this observation suggests that a persistent layer of subglacial water (>∼0.2 m) is widespread several kilometers upstream of the grounding zone. Our results are consistent with previous inferences of gradual grounding zones across this sector of the Ross Ice Shelf from airborne radar and satellite altimetry. Separately, the only clear bed-reflectivity change that we observed occurs ∼40 km downstream of the Kamb Ice Stream grounding zone, which we attribute to the onset of marine ice accretion onto the base of the ice shelf. This onset is much nearer to the grounding zone than previously predicted.

We revisit the results of the ISMIP-HEINO benchmark by first analyzing the differences in various model outputs using a wavelet-based spectral technique. Second, the ISMIP-HEINO benchmark experiments are recomputed with a novel numerical ice-sheet model based on the SIA-I algorithm that enables both the shallow-ice and a higher-order approximation of the ice-flow equations to be performed. To assess the significance of the higher-order approximation in the ISMIP-HEINO experiment, a numerical sensitivity study for the shallow-ice approximation (SIA) simulations is also carried out. A high sensitivity of the SIA model response to surface temperature perturbations is found. We conclude that the variations in ISMIP-HEINO results are due to the differences in (1) simulated basal temperatures and (2) numerical treatment of the basal sliding condition.

Time series of observed microwave brightness temperatures at Dome C, East Antarctic plateau, were modeled over 27 months with a multilayer microwave emission model based on dense-medium radiative transfer theory. The modeled time series of brightness temperature at 18.7 and 36.5 GHz were compared with Advanced Microwave Scanning Radiometer–EOS observations. The model uses in situ high-resolution vertical profiles of temperature, snow density and grain size. The snow grain-size profile was derived from near-infrared (NIR) reflectance photography of a snow pit wall in the range 850–1100 nm. To establish the snow grain-size profile, from the NIR reflectance and the specific surface area of snow, two empirical relationships and a theoretical relationship were considered. In all cases, the modeled brightness temperatures were overestimated, and the grain-size profile had to be scaled to increase the scattering by snow grains. Using a scaling factor and a constant snow grain size below 3 m depth (i.e. below the image-derived snow pit grain-size profile), brightness temperatures were explained with a root-mean-square error close to 1 K. Most of this error is due to an overestimation of the predicted brightness temperature in summer at 36.5 GHz.

The Greenland ice sheet (GrIS) subglacial hydrological system may undergo a seasonal evolution, with significant geophysical and biogeochemical implications. We present results from a new isotope-mixing model to quantify the relative contributions of surface snow, glacial ice and delayed flow to the bulk meltwater discharge from a small (∼5 km2) land-terminating GrIS outlet glacier during melt onset (May) and at peak melt (July). We use radioactive (222Rn) and stable isotopes (18O, deuterium) to differentiate the water source contributions. Atmospherically derived 7Be further constrains meltwater transit time from the glacier surface to the ice margin. We show that (1) 222Rn is a promising tracer for glacial waters stored at the bed and (2) a quantitative chemical mixing model can be constructed by combining 222Rn and the stable water isotopes. Applying this model to the bulk subglacial outflow from our study area, we find a constant delayed-flow (stored) component from melt onset through peak melt. This component is diluted first by snowmelt and then by increasing glacial ice melt as the season progresses. Results from this pilot study are consistent with the hypothesis that subglacial drainage beneath land-terminating sections of the GrIS undergoes a seasonal evolution from a distributed to a channelized system.

A novel ‘selection curve’ method is developed to interpret annual layers in the West Antarctic ice sheet (WAIS) Divide ice core based on dielectric properties (DEP). Because dielectric measurements are non-contact and represent the integrated response of the ice volume, they are particularly useful for the brittle zone of the core. Seasonal differences in ice chemistry create an annual signal in DEP, though multiple peaks of varying strength within a year may complicate the interpretation of annual layers. The selection curve algorithm uses a spline curve whose shape selects successive annual peaks in plots of DEP. This spline curve was scaled to the average annual-layer thickness at a given depth, where the layer thickness was best estimated using the fast Fourier transform (FFT) power spectrum within a sliding 10 m window. To explore the accuracy and stability of the method, several spline curves were generated from varying lengths of calibration data taken from multiple depths in the WAIS core. Using 50 m of manually interpreted calibration data, the selection curve method matched a manual interpretation throughout the entire 1200 m dataset to within 2% root-mean-square error (RMSE). This method is equally applicable to glaciochemical and other time/depth series measurements.

Glacier outlines are mapped for the upper Bhagirathi and Saraswati/Alaknanda basins of the Garhwal Himalaya using Corona and Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER) satellite images acquired in 1968 and 2006, respectively. A subset of glaciers was also mapped using Landsat TM images acquired in 1990. Glacier area decreased from 599.9 ± 15.6 km2 (1968) to 572.5 ± 18.0 km2 (2006), a loss of 4.6 ± 2.8%. Glaciers in the Saraswati/Alaknanda basin and upper Bhagirathi basin lost 18.4 ± 9.0 km2 (5.7 ± 2.7%) and 9.0 ± 7.7 km2 (3.3 ± 2.8%), respectively, from 1968 to 2006. Garhwal Himalayan glacier retreat rates are lower than previously reported. More recently (1990–2006), recession rates have increased. The number of glaciers in the study region increased from 82 in 1968 to 88 in 2006 due to fragmentation of glaciers. Smaller glaciers (<1 km2) lost 19.4 ± 2.5% (0.51 ± 0.07% a−1) of their ice, significantly more than for larger glaciers (>50 km2) which lost 2.8 ± 2.7% (0.074 ± 0.071 % a−1). From 1968 to 2006, the debris-covered glacier area increased by 17.8 ± 3.1% (0.46 ± 0.08% a−1) in the Saraswati/Alaknanda basin and 11.8 ± 3.0% (0.31 ± 0.08% a−1) in the upper Bhagirathi basin. Climate records from Mukhim (∼1900 m a.s.l.) and Bhojbasa (∼3780 m a.s.l.) meteorological stations were used to analyze climate conditions and trends, but the data are too limited to make firm conclusions regarding glacier–climate interactions.

The existence of solutions of the non-sliding shallow-ice-sheet equation on a flat horizontal bed with a mass balance linearly depending on altitude is proven for fixed margins. Free-margin solutions for the same mass balance do not exist. Fixed-margin solutions show unbounded shear stress and nonzero mass flux at the margin. Steady-state solutions with realistic margins, vanishing ice flux and vanishing shear stress are found numerically for ice sheets with Weertman-type sliding.

Snow slab avalanches primarily release by propagation of shear fractures within thin weak layers under much thicker slabs. In some cases, the weak layer is on the order of 1 mm thick and such fractures may be considered to be a mode II shear fracture at initiation. In the cases analysed in this paper, the weak layer has finite thickness, and slope-normal effects may be present. Field data from >500 snow shear fracture tests are analysed and applied to the problem of weak-layer fracture. The paper contains a detailed analysis using a simple analytical model to estimate the critical length prior to an unstable shear fracture. The model contains the assumption of a finite fracture process zone which may be a significant fraction of the slab depth, D, or the critical length, L, for weak-layer shear fracture. The field results show that the L/D ratio varies from ∼0.1 to just over 2, and the model results are close to the same range. The analysis also shows that both the field and model results for L/D follow a Gumbel probability density function. Since the experimental field data contain rate-dependent (viscoelastic) and slope-normal effects, it is imperative to account for these in the model and for snow slab instability evaluation when using test data. Detailed evaluations considering both these effects are given. The applicability of the test data to avalanche release is discussed.

Static (or ‘normal’) grain growth, i.e. grain boundary migration driven solely by grain boundary energy, is considered to be an important process in polar ice. Many ice-core studies report a continual increase in average grain size with depth in the upper hundreds of metres of ice sheets, while at deeper levels grain size appears to reach a steady state as a consequence of a balance between grain growth and grain-size reduction by dynamic recrystallization. The growth factor k in the normal grain growth law is important for any process where grain growth plays a role, and it is normally assumed to be a temperature-dependent material property. Here we show, using numerical simulations with the program Elle, that the factor k also incorporates the effect of the microstructure on grain growth. For example, a change in grain-size distribution from normal to log-normal in a thin section is found to correspond to an increase in k by a factor of 3.5.

Avalanche hazard evaluation relies in part on representative snowpack stability observations. Thus, understanding the spatial patterns of snowpack instabilities and their environmental determinants is crucial. This case study integrates intensive field observations with spatial modeling to identify associations between incoming radiation, surface hoar development and its subsequent shear strength across an inclined forest opening. We examined a buried surface hoar layer in southwest Montana, USA, over five sampling days, collecting 824 SnowMicroPen resistance profiles and performing 352 shear frame tests. Spatial models of incoming long- and shortwave radiation were generated for the surface hoar formation period using modeled hemispheric sky visibility, physically based parameters and the Bird Clear Sky Radiation Model in a Geographic Information System. Before burial, the surface hoar persisted despite moderate winds and relatively high air temperatures. The buried surface hoar layer thickness varied between 3 and 21 mm within a distance of 30 m. Modeled incoming radiation explained spatial variations in layer thickness and shear strength. In areas exposed to large amounts of radiation, the surface hoar layer was strong and thin, while areas with limited incoming radiation (due to high sky visibility and shading) possessed a thicker surface hoar layer that sheared more easily. This demonstrates the usefulness of microclimate modeling for slope-scale avalanche hazard evaluation. We also identify that over the 3 week sample period, strengthening occurred without thinning of the surface hoar layer.