We present an analysis of the active hydrologic system of MacAyeal Ice Stream (MacIS), West Antarctica, from a synthesis of multiple remote-sensing techniques: satellite laser altimetry; satellite image differencing; and hydrologic potential mapping (using a satellite-derived DEM and a bedrock DEM from airborne radio-echo sounding). Combining these techniques augments the information provided by each one individually, and allows us to develop a protocol for studying subglacial hydrologic systems in a holistic manner. Our study reveals five large active subglacial lakes under MacIS, the largest of which undergoes volume changes of at least 1.0 km3. We discuss the hydrologic properties of this system and present evidence for links between the lakes. At least three of the lakes are co-located with sticky spots, i.e. regions of high local basal shear stress. We also find evidence for surface elevation changes due to ice-dynamic effects (not just water movement) caused by changes in basal resistance. Lastly, we show that satellite radar altimetry is of limited use for monitoring lake activity on fast-flowing ice streams with surfaces that undulate on ∼10 km length scales.

Results from the Heinrich Event Intercomparison (HEINO) topic of the Ice-Sheet Model Intercomparison Project (ISMIP) are presented. ISMIP HEINO was designed to explore internal large-scale ice-sheet instabilities in different contemporary ice-sheet models. These instabilities are of interest because they are a possible cause of Heinrich events. A simplified geometry experiment reproduces the main characteristics of the Laurentide ice sheet, including the sedimented region over Hudson Bay and Hudson Strait. The model experiments include a standard run plus seven variations. Nine dynamic/thermodynamic ice-sheet models were investigated; one of these models contains a combination of the shallow-shelf (SSA) and shallow-ice approximation (SIA), while the remaining eight models are of SIA type only. Seven models, including the SIA–SSA model, exhibit oscillatory surges with a period of ∼1000 years for a broad range of parameters, while two models remain in a permanent state of streaming for most parameter settings. In a number of models, the oscillations disappear for high surface temperatures, strong snowfall and small sediment sliding parameters. In turn, low surface temperatures and low snowfall are favourable for the ice-surge cycles. We conclude that further improvement of ice-sheet models is crucial for adequate, robust simulations of cyclic large-scale instabilities.

This paper records the history of how and why the Journal of Glaciology was started and how it has grown into the Journal we know today.

We report the discovery in the Greenland ice sheet of a discrete layer of free nanodiamonds (NDs) in very high abundances, implying most likely either an unprecedented influx of extraterrestrial (ET) material or a cosmic impact event that occurred after the last glacial episode. From that layer, we extracted n-diamonds and hexagonal diamonds (lonsdaleite), an accepted ET impact indicator, at abundances of up to about 5×106 times background levels in adjacent younger and older ice. The NDs in the concentrated layer are rounded, suggesting they most likely formed during a cosmic impact through some process similar to carbon-vapor deposition or high-explosive detonation. This morphology has not been reported previously in cosmic material, but has been observed in terrestrial impact material. This is the first highly enriched, discrete layer of NDs observed in glacial ice anywhere, and its presence indicates that ice caps are important archives of ET events of varying magnitudes. Using a preliminary ice chronology based on oxygen isotopes and dust stratigraphy, the ND-rich layer appears to be coeval with ND abundance peaks reported at numerous North American sites in a sedimentary layer, the Younger Dryas boundary layer (YDB), dating to 12.9 ± 0.1 ka. However, more investigation is needed to confirm this association.

While performing statistical–dynamical simulations for avalanche predetermination, a propagation model must reach a compromise between precise description of the avalanche flow and computation times. Crucial problems are the choice of appropriate distributions describing the variability of the different inputs/outputs and model identifiability. In this study, a depth-averaged propagation model is used within a hierarchical Bayesian framework. First, the joint posterior distribution is estimated using a sequential Metropolis–Hastings algorithm. Details for tuning the estimation algorithm are provided, as well as tests to check convergence. Of particular interest is the calibration of the two coefficients of a Voellmy friction law, with model identifiability ensured by prior information. Second, the point estimates are used to predict the joint distribution of different variables of interest for hazard mapping. Recent developments are employed to compute pressure distributions taking into account the rheology of snow. The different steps of the method are illustrated with a real case study, for which all possible decennial scenarios are simulated. It appears that the marginal distribution of impact pressures is strongly skewed, with possible high values for avalanches characterized by low Froude numbers. Model assumptions and results are discussed.

Information on bed topography and basal conditions is essential to developing the next-generation ice-sheet models needed to generate a more accurate estimate of ice-sheet contribution to sea-level rise. Synthetic aperture radar (SAR) images of the ice–bed can be analyzed to obtain information on bed topography and basal conditions. We developed a wideband SAR, which was used during July 2005 to perform measurements over a series of tracks between the GISP2 and GRIP cores near Summit Camp, Greenland. The wideband SAR included an eight-element receive-antenna array with multiple-phase centers. We applied the MUltiple SIgnal Classification (MUSIC) algorithm, which estimates direction of arrival signals, to single-pass multichannel data collected as part of this experiment to obtain fine-resolution bed topography. This information is used for producing fine-resolution estimates of bed topography over a large swath of 1600 m, with a 25 m posting and a relative accuracy of approximately 10 m. The algorithm-derived estimate of ice thickness is within 10 m of the GRIP ice-core length. Data collected on two parallel tracks separated by 500 m and a perpendicular track are compared and found to have difference standard deviations of 9.1 and 10.3 m for the parallel and perpendicular tracks, respectively.

Calving from tidewater glaciers and ice shelves accounts for around half the mass loss from both polar ice sheets, yet the process is not well represented in prognostic models of ice dynamics. Benn and others proposed a calving criterion appropriate for both grounded and floating glacier tongues or ice shelves, based on the penetration depth of transverse crevasses near the calving front, computed using Nye’s formula. The criterion is readily incorporated into glacier and ice-sheet models, but has not been fully validated with observations. We apply a three-dimensional extension of Benn and others’ criterion, incorporated into a full-Stokes model of glacier dynamics, to estimate the current position of the calving front of Johnsons Glacier, Antarctica. We find that two improvements to the original model are necessary to accurately reproduce the observed calving front: (1) computation of the tensile deviatoric stress opening the crevasse using the full-stress solution and (2) consideration of such a tensile stress as a function of depth. Our modelling results also suggest that Johnsons Glacier has a polythermal structure, rather than the temperate structure suggested by earlier studies.

Snow and weather conditions typical of exceptional cycles of fresh-snow avalanches in the northern Alps are investigated using the numerical avalanche-hazard forecasting procedure of Météo-France. Sensitivity tests are performed on the events of February 1999 in the Chamonix France region, and resulting snowpack instability modeled at the massif scale is compared using adapted new indices and maps. Our results complete conclusions of earlier observation-based studies by providing new insights into the snow and weather conditions of February 1999. The large avalanches mainly resulted from large and very unstable fresh-snow accumulations. Moreover, the snowpack instability was increased locally by wind transport of light and fresh snow in February. The mechanical weaknesses resulting from the weather conditions prior to February were a key factor in explaining the unusual volumes of these avalanches. This study suggests that the operational numerical SAFRAN/Crocus/ MÉPRA (SCM) chain provides reliable forecasts of extreme new-snow avalanche situations at the massif scale, but that local-scale simulations are still needed to improve the efficiency of risk mitigation and civil protection policies.

Women’s history in glaciology extends as far back in time as the discipline itself, although their contributions to the scientific discourse have for all of that history been constrained by the sociopolitical contexts of the times. The first Journal of Glaciology paper authored by a woman appeared in 1948, within a year of the founding of the Journal, but it was not until the 1980s that women produced more than a few percent of Journal and Annals of Glaciology papers. Here international perspectives on women’s participation in the sciences are presented in order to establish an economic and sociopolitical context for stories of women ‘pioneers’ in glaciology and a frame in which to discuss women’s persistent under-representation relative to men. We find that the experiences of individual glaciologists mirror women’s experiences in higher education and the sciences as a whole. The existence of both positive and negative trends in women’s participation in the sciences suggests caution in the interpretation of recent positive trends for women’s participation in glaciology.

We report impact pressures exerted by three wet-snow avalanches on a pylon equipped with piezoelectric load cells. These pressures were considerably higher than those predicted by conventional avalanche engineering guidelines. The time-averaged pressure linearly increased with the immersion depth of the load cells and it was about eight times larger than the hydrostatic snow pressure. At the same immersion depth, the pressures were very similar for all three avalanches and no dependency between avalanche velocity and pressure was apparent. The pressure time series were characterized by large fluctuations. For all immersion depths and for all avalanches, the standard deviations of the fluctuations were, on average, about 20% of the mean value. We compare our observations with results of slow-drag granular experiments, where similar behavior has been explained by formation and destruction of chain structures due to jamming of granular material around the pylon, and we propose the same mechanism as a possible microscale interpretation of our observations.

Firn microstructure is accurately characterized using images obtained from scanning electron microscopy (SEM). Visibly etched grain boundaries within images are used to create a skeleton outline of the microstructure. A pixel-counting utility is applied to the outline to determine grain area. Firn grain sizes calculated using the technique described here are compared to those calculated using the techniques of Gow (1969) and Gay and Weiss (1999) on samples of the same material, and are found to be substantially smaller. The differences in grain size between the techniques are attributed to sampling deficiencies (e.g. the inclusion of pore filler in the grain area) in earlier methods. The new technique offers the advantages of greater accuracy and the ability to determine individual components of the microstructure (grain and pore), which have important applications in ice-core analyses. The new method is validated by calculating activation energies of grain boundary diffusion using predicted values based on the ratio of grain-size measurements between the new and existing techniques. The resulting activation energy falls within the range of values previously reported for firn/ice.

Here we report a laboratory study of the effects of debris thickness, diurnally cyclic radiation and rainfall on melt rates beneath rock-avalanche debris and sand (representing typical highly permeable supraglacial debris). Under continuous, steady-state radiation, sand cover >50 mm thick delays the onset of ice-surface melting by >12 hours, but subsequent melting matches melt rates of a bare ice surface. Only when diurnal cycles of radiation are imposed does the debris reduce the longterm rate of ice melt beneath it. This is because debris >50 mm thick never reaches a steady-state heat flux, and heat acquired during the light part of the cycle is partially dissipated to the atmosphere during the nocturnal part of the cycle, thereby continuously reducing total heat flux to the ice surface underneath. The thicker the debris, the greater this effect. Rain advects heat from high-permeability supraglacial debris to the ice surface, thereby increasing ablation where thin, highly porous material covers the ice. In contrast, low-permeability rock-avalanche material slows water percolation, and heat transfer through the debris can cease when interstitial water freezes during the cold/night part of the cycle. This frozen interstitial water blocks heat advection to the ice–debris contact during the warm/day part of the cycle, thereby reducing overall ablation. The presence of metre-deep rock-avalanche debris over much of the ablation zone of a glacier can significantly affect the mass balance, and thus the motion, of a glacier. The length and thermal intensity of the diurnal cycle are important controls on ablation, and thus both geographical location and altitude significantly affect the impact of debris on glacial melting rates; the effect of debris cover is magnified at high altitude and in lower latitudes.

The temporal variability of surface snow and glacier melt flux and runoff are investigated for the ablation area of Jakobshavn Isbræ, West Greenland. High-resolution meteorological observations both on and outside the Greenland ice sheet were used as model input. SnowModel, a physically based spatially distributed meteorological and snow evolution modeling system, is used to simulate the temporal variability of Jakobshavn Isbræ accumulation and ablation processes for 2000/01–2006/07. Winter snow depth observations and MODIS satellite-derived summer melt observations are used for model validation of accumulation and ablation. The modeled interannual runoff variability varied from 1.81 × 109 m3 (2001/02) to 5.21 × 109 m3 (2004/05), yielding a cumulative runoff at the Jakobshavn Glacier terminus of ∼2.25 to ∼4.5 m w.e. The average modeled Jakobshavn runoff of ∼3.4 km3 a−1 was merged with previous estimates of Jakobshavn ice discharge to quantify the freshwater flux to Illulissat Icefjord. For both runoff and ice discharge the average trends are similar, indicating increasing (insignificant) influx of fresh water to Ilulissat Icefjord for the period 2000/01–2006/07. This study suggests that surface runoff forms a minor part of the overall Jakobshavn freshwater flux to the fjord: about 7% (∼3.4 km3 a−1) of the average annual freshwater flux of ∼51.0 km3 a−1 originates from the surface runoff.

Uncertainties in estimates of glacier and ice-cap contribution to sea-level rise exist in part due to poor quantification of mass-balance errors, particularly those resulting from extrapolation of sparse measurements. Centre-line data are often assumed to be representative of the glacier as a whole, with little attention paid to extrapolation errors or their effect on mass-balance estimates. Here we present detailed digital elevation model (DEM) measurements of glacier-wide elevation changes over the last ~40 years at two glaciers on Svalbard, Norwegian Arctic. Austre Br0ggerbreen and Midtre Lovenbreen are shown to have lost 27.54 ± 0.98 and 9.65 ± 0.76 × 107m3 of ice, respectively, between 1966 and 2005, findings that we relate to trends in average summer air temperatures and winter accumulation. These volume losses correspond to geodetic balances of -0.58 ± 0.03 and -0.41 ± 0.03 mw.e. a-1, respectively. Our analysis revealed high spatial complexity in patterns of elevation change, varying between glaciers, between measurement intervals and within and between elevation bins. Balances from extrapolated centre-line geodetic data were the same (within errors) as those from full-coverage DEM differencing in the majority of comparisons, yet significantly underestimated balance in three instances. Additionally, field mass balance from centre-line ablation stake data underestimated balances from full-coverage geodetic measurements during three of six measurement periods. These findings may support the hypothesis that field measurements underestimate Svalbard glacier mass loss, at least partly as a result of the failure of centre-line measurements to account for glacier-wide variations in ablation. Our results demonstrate the importance of deriving accurate interpolation functions and constraining extrapolation errors from sparse measurements.

This paper reviews the seminal contributions (and tussles) of Weertman and Lliboutry to the theory of glacier sliding in the 1950s and 1960s, and summarizes later developments up to the present day.

Temperate glaciers are more sensitive to climate changes than polar or continental glaciers, and can drive remarkable runoff variation in local water catchments. Here we present recent glacier shrinkage and runoff change for Hailuogou glacier, a typical monsoon temperate glacier on the east slope of Mount Gongga (Minya Konga), China. The surface area of Hailuogou glacier has decreased by 3.5% (0.92 km2) between 1966 (aerial photographs) and 2007 (ASTER images). Flow measurements at a stream gauge about 500 m down-glacier commencing in 1994 display a remarkable increase in annual runoff (mostly during July–September) since 1999. Annual runoff over the same period in a nonglacierized but forested subcatchment (9.17 km2) did not experience significant change. By separating the daily rainfall component from the daily total discharge, monthly catchment water-balance series were calculated for the period 1994–2005, which shows an increasing trend of glacier storage loss. We concluded that air-temperature rise (with a trend of +0.2°C (10 a)−1 between 1988 and 2005, recorded at nearby weather stations) has had an increased effect on glacier mass loss and river runoff change during the past 20 years.

We modelled the flow of the Larsen C and northernmost Larsen D ice shelves, Antarctic Peninsula, using a model of continuum mechanics of ice flow, and applied a fracture criterion to the simulated velocities to investigate the ice shelf’s present-day stability. Constraints come from satellite data and geophysical measurements from the 2008/09 austral summer. Ice-shelf thickness was derived from BEDMAP and ICESat data, and the density–depth relationship was inferred from our in situ seismic reflection data. We obtained excellent agreements between modelled and measured ice-flow velocities, and inferred and observed distributions of rifts and crevasses. Residual discrepancies between regions of predicted fracture and observed crevasses are concentrated in zones where we assume a significant amount of marine ice and therefore altered mechanical properties in the ice column. This emphasizes the importance of these zones and shows that more data are needed to understand their influence on ice-shelf stability. Modelled flow velocities and the corresponding stress distribution indicate that the Larsen C ice shelf is stable at the moment. However, weakening of the elongated marine ice zones could lead to acceleration of the ice shelf due to decoupling from the slower parts in the northern inlets and south of Kenyon Peninsula, leading to a velocity distribution similar to that in the Larsen B ice shelf prior to its disintegration.

Ice originating near the inland ice divide of the ice sheet can reappear as marginal ice at the surface near the ice terminal in the ablation area. We have analyzed δ18O values and ion concentrations of the Skallen, Skarvsnes and Hamna terminal ice sections, located along the estuary line in the Sôya drainage basin, East Antarctica. The data suggest that the upper part of the Skallen terminal ice section could have originated from inland precipitation on the Shirase drainage basin during marine isotope stage (MIS) 5e, while the upper part of Skarvsnes and Hamna terminal ice sections could have originated from inland precipitation on the Sôya drainage basin. We calculate past elevation maps for the Antarctic ice sheet using the three-dimensional model, SICOPOLIS. This model suggests that the upstream portion of the Sôya drainage basin during the glacial period (MIS 2, 3 or 4) was located to the northeast of its present location. A flow history is proposed wherein ice from the inland Shirase drainage area flowed over the present ice-divide line from the Shirase to the Sôya drainage basin during the glacial period. The ice in the Sôya drainage basin then flowed to the marginal part of the sheet after the ice divide had assumed its present position.

The volume- and area-change evolution of glaciers can be obtained by employing the volume–area scaling approach during mass-balance modelling. This method usually requires information on the initial surface area and ice volume to adjust the volume–area relation to the specific ice body. However, absolute volumetric data on glaciers are very rare, so the applicability of volume–area scaling is limited. In order to use volume–area scaling on glaciers for which only limited information is available, a new method is presented to calibrate the volume–area relation without prior knowledge of this relation by using glacier extent information from different times. To validate the method and illustrate its practicability, we model the range of probable future changes in ice volume and surface area of ‘Glaciar Noroeste’, an outlet glacier of Gran Campo Nevado ice cap, southern Chilean Patagonia, during the 21st century, based on IPCC SRES scenarios B1 and A2.

For the first time, to our knowledge, a scientific study of the formation and evolution of waterfall ice, the ephemeral ice structures that form from the freezing of liquid water seeping on steep rock, was performed. We surveyed and analysed three waterfall ice structures near Glacier d’Argentière, Mont Blanc massif, France, between winter 2007 and spring 2009. We reconstruct the global evolution of two vertical ice structures using automatic digital cameras, while the internal ice microstructure was analysed using ice coring and sampling. Macro- and microstructural observations are considered, along with temperature conditions recorded at a nearby meteorological station and directly within the ice structure. They reveal that vertical structures initially grow rapidly from the aggregation of stalactites with microstructures indicative of temperature conditions during their crystallization. After this initial stage, the volume of the ice structure reaches an asymptotic value, as water continues to flow inside the structure, isolated from the outside cold ice; the outer surface remains dry. At the end of the season, the collapse of the free-standing structure does not occur by progressive melting, but is initiated by a horizontal crack propagation at the top. The initiation of this crack seems to be triggered by a drastic temperature decrease.