We have determined the in situ electric field attenuation length Lα (defined as the length over which the signal amplitude diminishes by a factor 1/e) for radio-frequency signals broadcast vertically through South Polar ice and reflected off the underlying bed. Conservatively assuming a bedrock field reflectivity for f = 380 MHz, and T = –50°C; the errors incorporate uncertainties in R. This value is consistent with previous estimates that the radiofrequency attenuation length exceeds the attenuation length at optical frequencies by an order of magnitude.

The dynamics of glacier motion are governed to a large extent by the properties of the basal interface. In this paper we address the interaction of a glacier with a layer of till at its bed in an attempt to test whether our physical understanding of till is sufficient to explain general features of the observed flow field and changes in geometry of Black Rapids Glacier, Alaska, U.S.A. We also investigate whether or not a till layer has a clear surface-observable signature in the dynamics of the glacier. Towards this end we use a finite-element ice-flow model with a Coulomb failure criterion within the basal till layer. We find that simple “till physics” can be used to describe decadal, seasonal and short-term (hours to days) velocity variations, and possibly uplift events. Mechanisms for each of these variations involve an increase in the extent of till at failure, a transfer of shear stress across the bed, and a consequent increase in ice deformation. “Effective shape factors” are calculated that permit a simple incorporation of this boundary condition into glacier response models. Our analyses, however, have not resulted in the identification of a clear and unique signature of a till layer in the surface dynamics of a glacier.

The ApRES (autonomous phase-sensitive radio-echo sounder) instrument is a robust, lightweight and relatively inexpensive radar that has been designed to allow long-term, unattended monitoring of ice-shelf and ice-sheet thinning. We describe the instrument and demonstrate its capabilities and limitations by presenting results from three trial campaigns conducted in different Antarctic settings. Two campaigns were ice sheet-based – Pine Island Glacier and Dome C – and one was conducted on the Ross Ice Shelf. The ice-shelf site demonstrates the ability of the instrument to collect a time series of basal melt rates; the two grounded ice applications show the potential to recover profiles of vertical strain rate and also demonstrate some of the limitations of the present system.

Snowmelt regions on Greenland ice are mapped daily with the SeaWinds wideswath Ku-band (13.4 GHz) scatterometer on the QuikSCAT satellite. The approach exploits the high temporal resolution of SeaWinds/QuikSCAT data for the melt mapping using diurnal backscatter change independent of the absolute calibration. The results reveal several pronounced melting and refreezing events, and effects of topography are evident in the melt patterns. The spatial resolution is sufficient to identify melt features on the Sukker-toppen Iskappe west of the main ice sheet. An anomalous warming event, caused by down-ward mixing of warm air, is detected in late September 1999 over the west flank of the southern Greenland ice sheet. Time-series images of melt regions are presented over the period from summer to the fall freeze-up. The satellite observations are verified with in situ measurements from the Greenland Climate Network stations.

Karakoram glaciers experienced balanced or slightly positive mass budgets since at least the 1970s. Here, we provide an update on the state of balance of Central and Eastern Karakoram glaciers (12 000 km2) between 2008 and 2016 by differencing DEMs derived from satellite optical images. The mass budget of Central Karakoram glaciers was slightly positive (0.12 ± 0.14 m w.e. a−1) while eastern Karakoram glaciers lost mass (−0.24 ± 0.12 m w.e. a−1). The glacier-wide mass balances of surge-type and nonsurge-type glaciers were not statistically different. Our elevation change data also depict the effect of a > 100 Mm3 rock avalanche on Siachen Glacier ablation area in September 2010. It covered a 4 km2 area with a thick debris layer that unexpectedly, led to locally enhanced glacier mass loss during the following years. Enhanced melt opened a > 100 m deep 2 km2 depression and contributed to 6% of the mass loss of Siachen Glacier from 2010 to 2016 (−0.39 m w.e. a−1). We hypothesize that sub- or englacial melt may be responsible for this intriguing behaviour. This study contributes to a better knowledge of the regional pattern of the Karakoram anomaly and of the influence of rock avalanches on glacier mass changes.

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.

Polar ice is known to be one of the most anisotropic natural materials. For a given fabric the polycrystal viscous response is strongly dependent on the actual state of stress and strain rate. Within an ice sheet, grounded-ice parts and ice shelves have completely different stress regimes, so one should expect completely different impacts of ice anisotropy on the flow. The aim of this work is to quantify, through the concept of enhancement factors, the influence of ice anisotropy on the flow of grounded ice and ice shelves. For this purpose, a full-Stokes anisotropic marine ice-sheet flowline model is used to compare isotropic and anisotropic diagnostic velocity fields on a fixed geometry. From these full-Stokes results, we propose a definition of enhancement factors for grounded ice and ice shelves, coherent with the asymptotic models used for these regions. We then estimate realistic values for the enhancement factors induced by ice anisotropy for grounded ice and ice shelves.

In the framework of the European Project for Ice Coring in Antarctica (EPICA), a comprehensive glaciological pre-site survey has been carried out on Amundsenisen, Dronning Maud Land, East Antarctica, in the past decade. Within this survey, four intermediate-depth ice cores and 13 snow pits were analyzed for their ionic composition and interpreted with respect to the spatial and temporal variability of volcanic sulphate deposition. The comparison of the non-sea-salt (nss)-sulphate peaks that are related to the well-known eruptions of Pinatubo and Cerro Hudson in AD 1991 revealed sulphate depositions of comparable size (15.8±3.4 kg km–2) in 11 snow pits. There is a tendency to higher annual concentrations for smaller snow-accumulation rates. The combination of seasonal sodium and annually resolved nss-sulphate records allowed the establishment of a time-scale derived by annual-layer counting over the last 2000 years and thus a detailed chronology of annual volcanic sulphate deposition. Using a robust outlier detection algorithm, 49 volcanic eruptions were identified between AD 165 and 1997. The dating uncertainty is ±3 years between AD 1997 and 1601, around ±5 years between AD 1601 and 1257, and increasing to ±24 years at AD 165, improving the accuracy of the volcanic chronology during the penultimate millennium considerably.

The flow of valley glaciers is examined in the light of recent laboratory experiments on the behaviour of ice under load. Simple expressions are given for the velocity distributions in some cases of laminar flow, and the modification of a pure laminar flow theory necessary to explain the formation of transverse crevasses and thrust planes is considered. The paper ends with some remarks about the formation of crevasse patterns on the surfaces of glaciers. The statical equilibrium of a circular ice cap is discussed in an appendix.

We present results from cold-laboratory observations of changes in isotopic (δ18O and δD) content by sublimation in snow and ice samples under nearly isothermal conditions. The results show large increases in observed δ18O and δD in snow samples within several centimeters of the surface. They contradict the assumption of a non-changing isotopic content due to layer-by-layer transport mechanisms driven by sublimation/desublimation processes. The data also do not support the idea that isotopic changes of snow and firn are limited by the possibility that the ice matrix incorporates the atmospheric water vapor and that forced water-vapor diffusion in the pore space (wind pumping) is a requirement for isotopic content change. The observations show that sublimation from ice samples results in much lower increases in heavy-isotope content in the first several millimetres near the sublimating surface over the same time period, despite sublimation intensities similar to those of the snow samples. The results suggest that continuous phase transitions inside snow (recrystallization) are the process responsible for the isotopic content change because they are the primary mass-exchange mechanism between the snow mass and the surrounding environment. Modeling the isotopic content of the ice matrix therefore requires inclusion of a two-stage process: fractionation at the ice-matrix surface due to repetitive phase transitions, and fractionation due to preferable diffusion of light water isotopes in the pore space. For interpretation of the observed natural isotopic profiles in snow, the first process can be linked to the time a snow layer undergoes recrystallization, while the second process is related to the total ice/snow mass gain/loss determined by the external environmental conditions.

X-ray diffraction photographs were taken of ice dendrites grown from water and from the vapour. In both cases [1120] was found to be the growth direction.

A method of growing cylindrical single crystals of any desired orientation has been devised. Some of these crystals were extended to several times their original length. Slip took place exclusively on the basal plane. Experiments designed to locate the glide direction gave inconclusive results and suggested that, in general, crystals near the melting point can glide in any direction.

Before deformation, ice crystals gave fairly sharp X-ray diffraction spots on Laue photographs. After deformation the spots are drawn out into streaks which later break up into rows of small spots. This indicates that bending of the crystals is followed by the formation of small crystallites of slightly differing orientations. Analogous effects are observed in metals.

Basal slip accounts for a large part of the flow of ice streams draining ice from Antarctica and Greenland into the ocean. Therefore, an appropriate representation of basal slip in ice flow models is a prerequisite for accurate sea level rise projections. Various friction laws have been proposed to describe basal slip in models. Here, we compare the influence on grounding line (GL) dynamics of four friction laws: the traditional Weertman law and three effective pressure-dependent laws, namely the Schoof, Tsai and Budd laws. It turns out that, even when they are tuned to a common initial reference state, the Weertman, Budd and Schoof laws lead to thoroughly different steady-state positions, although the Schoof and Tsai laws lead to much the same result. In particular, under certain circumstances, it is possible to obtain a steady GL located on a reverse slope area using the Weertman law. Furthermore, the predicted transient evolution of the GL as well as the projected contributions to sea level rise over a 100-year time horizon vary significantly depending on the friction law. We conclude on the importance of choosing an appropriate law for reliable sea level rise projections and emphasise the need for a coupling between ice flow models and physically based subglacial hydrological models.

The Himalaya has some of the largest glacier concentrations outside the polar regions. Despite this, long-term measurements detecting the impact of global warming and changing precipitation patterns on glaciers are rare. The Nanga Parbat massif in northern Pakistan is an exception. The cartographer and glaciologist R. Finsterwalder investigated glacier dynamics of this mountain massif in the 1930s, and several other studies document changes since then. The aim of this study is to detect and analyse the changes of Raikot Glacier over the past seven decades. We use a multitemporal and multiscale approach, based on repeat terrestrial images, additional historical data and remotely sensed imagery (Corona, ASTER, Landsat, QuickBird). The multitemporal approach covers the period 1934–2007. While the analyses show a total glacier retreat of ~200 m in 73 years, this general trend was interrupted by a significant glacier advance between the 1950s and 1980s. Although down-wasting processes can be inferred from an increase in debris-covered area, a general trend of reduced glacier thickness does not appear significant over the whole observation period.

Despite their relatively small total ice volume, mid-latitude valley glaciers are expected to make a significant contribution to global sea-level rise over the next century due to the sensitivity of their mass-balance systems to small changes in climate. Here we use a degree-day model to reconstruct the past century of mass-balance variation at ‘Ka Roimata o Hine Hukatere’ Franz Josef Glacier, New Zealand, and to predict how mass balance may change over the next century. Analysis of the relationship between temperature, precipitation and mass balance indicates that temperature is a stronger control than precipitation on the mass balance of Franz Josef Glacier. The glacier’s mass balance, relative to its 1986 geometry, has decreased at a mean annual rate of 0.02m a–1 w.e. between 1894 and 2005. We compare this reduction to observations of terminus advance and retreat, of which Franz Josef Glacier has the best record in the Southern Hemisphere. For the years 2000–05 the relative mass balance ranged from –0.75 to +1.50m a–1 w.e., with 2000/01 the only year showing a negative mass balance. In a regionally downscaled Intergovernmental Panel on Climate Change mean warming scenario, the annual relative mass balance will continue to decrease at 0.02m a–1 w.e. through the next century.

This study assesses the impact of a large 2010 calving event on the current and future stability of Petermann Glacier, Greenland, and ascertains the glacier’s interaction with different components of the climate and ocean system. We use a numerical ice-flow model that captures the major aspects of the glacier’s mass budget, the resistive forces controlling glacier flow, and includes dynamic calving. Satellite observations and model results show that the recent break-off of 25% of the floating tongue did not result in a significant glacier speed-up due to the low lateral resistance of this relatively wide and thin ice tongue. We demonstrate that seasonal speed-up at Petermann Glacier is mainly driven by meltwater lubrication rather than freeze-up conditions in the fjord. Results also show that sub-shelf ocean melt may have a profound effect on the future stability of Petermann Glacier, emphasizing the urgent need for more observations, and a better understanding of fjord temperature variability and circulation.

The specific surface area (SSA), defined as the surface area of ice per unit mass, is an important variable characterizing the complex microstructure of snow. Its application range covers the physical evolution of snow (metamorphism), photochemistry and optical and microwave remote sensing. This paper presents a new device, POSSSUM (Profiler Of Snow Specific Surface area Using SWIR reflectance Measurement), designed to allow the rapid acquisition of SSA profiles down to ∼20 m depth and with an effective vertical resolution of 10–20 mm. POSSSUM is based on the infrared (IR) reflectance technique: A laser diode operating at 1310 nm illuminates the snow at nadir incidence angle along the face of a drilled hole. The reflected radiance is measured at three zenith angles (20°, 40° and 60°) each for two azimuth angles (0° and 180°). A second laser operating at a shorter wavelength (635 nm), which is almost insensitive to SSA, allows the distance to the snow face to be estimated. The reflected IR radiance and the distance are combined to estimate bidirectional reflectances. These reflectances are converted into hemispherical reflectances and in turn into SSA using a theoretical formulation based on an asymptotic solution of the radiative transfer equation. The evaluation and validation of POSSSUM’s SSA measurements took place in spring 2009 in the French Alps. The new method was compared with the methane adsorption technique and DUFISSS, another well-validated instrument based on the IR technique. The overall measurement error is in the range 10–15%.

Avalanche deposits consist of rounded granules composed of aggregates of snow and ice particles. The size of the granules is related to vertical shear gradients within the flow; studying the granule-size distribution may be useful in understanding the flow and stopping of avalanches. We applied a sediment-size sampling method to measure snow granule-size distributions at different depositional environments on two dry and two wet avalanche deposits at three field sites. The granule-size distributions are approximately log-normal, similar to many natural sediment deposits. The median granule size in the wet and dry avalanches varies between 65 and 162 mm. Wet avalanches tend to produce more large granules than dry avalanches, indicating both smaller flow velocities and near-surface shear gradients. Granule size is similar in frontal lobes and levee deposits, suggesting that levee formation occurs independently of the size segregation at the avalanche front.

The altitude of the 0°C isotherm obtained from radiosonde data of the aerological Chilean stations Antofagasta, Quintero/Santo Domingo, Puerto Montt and Punta Arenas are analyzed, along with surface temperature and precipitation records from nearby stations. The strong effect of the 1976/77 climate shift due to a change in the Pacific Decadal Oscillation is evident in the temperature and precipitation data. The data are used as input for an empirical model which reconstructs annually the equilibrium-line altitude (ELA) for the last 49 years on the western side of the southern Andes. The model takes air temperature, precipitation and altitude as the main parameters, and was first developed by Fox (1993) and applied by Condom and others (2007). From the radiosonde data, a significant positive trend of the 0°C isotherm has occurred in the northern, central and southern regions, indicating an ELA rise due to regional warming. General glacier retreat, ice thinning and negative mass balance observed during the past few decades in virtually all the Chilean Andes concur with the observed ELA reconstruction. In the Punta Arenas radiosonde record there is slight evidence for precipitation increase but no evidence for significant warming in the past few decades. This results in a slight lowering of the ELA according to the model reconstruction, which does not agree with the strong and increased glacier retreat observed in recent decades in Patagonia.

Nearly 100 days of hourly glacier motion, hydrology and hydrochemistry measurements on Kennicott Glacier, Alaska, USA, demonstrate the complicated relationship between water and motion at the glacier bed. Our observations capture the transient glacier response to seasonal and daily melt cycles, and to a jökulhlaup that prompts a sixfold increase in glacier speed. Sliding is promoted whenever the water inputs to the glacier exceed the capacity of the subglacial hydrologic system to transmit the water. Sensitivity of sliding to daily meltwater inputs varies strongly through the season, implying that the state of the hydrologic system governs the sensitivity of basal sliding. A numerical model constructed to explore these relationships reveals: the roles of the effective pressure; the exponent to which this is taken in the ‘sliding law’ (0.1 <γ< 0.6); glacier macroporosity (φ < 2%); and the ‘cavity-generating capacity’ of the glacier bed, which encapsulates the sizes and spacing of bed roughness elements. Temporal changes in the effective pressure associated with evolution of both water inputs and subglacial water transmission capacity can explain the varying strength of diurnal velocity fluctuations of Kennicott Glacier. Spatial patterns of glacier macroporosity and of basal roughness govern variation in sensitivity of sliding to water inputs.