In this paper, we explore geomorphological evidence allowing a first- order reconstruction of the extent and pattern of frozen-bed conditions under the last Fennoscandian ice sheet. We mapped relict landscapes, i.e. glacial landforms and subaerially developed ground surfaces predating the last ice sheet and marking sustained frozen- bed conditions, at four different spatial scales. At the ice-sheet scale, relict landscapes are most abundant between the Last Glacial Maximum ice divide and the elevation axis of the Scandinavian mountain range. The location of frozen-bed zones was mainly a function of dispersal centre location (low surface temperatures and small strain heating) and small ice thickness over the eastern flank of the mountain range. At the mesoscale (260 × 360 km map area), the pattern of relict surfaces is governed by inward-cutting ice-stream erosion.Topographical control was weak, but relation to flow pattern was strong, with the major frozen-bed zone located where ice flow was strongly divergent. At the regional scale (40 × 65 km map area) in hilly terrain, topographical control was strong with relict surfaces only appearing above a plane dipping in the up-ice direction. At the local scale (12 × 14 km map area), control by topography was likewise strong, but the detailed boundary pattern was irregular, with specific landforms occurring both up- and down-ice of frozen patches.

Thirteen annually resolved accumulation-rate records covering the last ~200 years from the Pine Island–Thwaites and Ross drainage systems and the South Pole are used to examine climate variability over West Antarctica. Accumulation is controlled spatially by the topography of the ice sheet, and temporally by changes in moisture transport and cyclonic activity. A comparison of mean accumulation since 1970 at each site to the long-term mean indicates an increase in accumulation for sites located in the western sector of the Pine Island–Thwaites drainage system. Accumulation is negatively associated with the Southern Oscillation Index (SOI) for sites near the ice divide, and periods of sustained negative SOI (1940–42, 1991–95) correspond to above-mean accumulation at most sites. Correlations of the accumulation-rate records with sea-level pressure (SLP) and the SOI suggest that accumulation near the ice divide and in the Ross drainage system may be associated with the mid-latitudes. The post-1970 increase in accumulation coupled with strong SLP–accumulation-rate correlations near the coast suggests recent intensification of cyclonic activity in the Pine Island– Thwaites drainage system.

Sea-ice freeboard heights for 17 ICESat campaign periods from 2003 to 2009 are derived from ICESat data. Freeboard is combined with snow depth from Advanced Microwave Scanning Radiometer for Earth Observing System (AMSR-E) data and nominal densities of snow, water and sea ice, to estimate sea-ice thickness. Sea-ice freeboard and thickness distributions show clear seasonal variations that reflect the yearly cycle of growth and decay of the Weddell Sea (Antarctica) pack ice. During October–November, sea ice grows to its seasonal maximum both in area and thickness; the mean freeboards are 0.33–0.41m and the mean thicknesses are 2.10–2.59 m. During February–March, thinner sea ice melts away and the sea-ice pack is mainly distributed in the west Weddell Sea; the mean freeboards are 0.35–0.46m and the mean thicknesses are 1.48–1.94 m. During May–June, the mean freeboards and thicknesses are 0.26–0.29m and 1.32–1.37 m, respectively. the 6 year trends in sea-ice extent and volume are (0.023±0.051)×106 km2 a–1 (0.45% a–1) and (0.007±0.092)×103 km3 a–1 (0.08% a–1); however, the large standard deviations indicate that these positive trends are not statistically significant.

This paper presents an assessment of biological activity associated with ice surface debris (cryoconite) at the ice-sheet scale. Estimates of the mass distribution of cryoconite over the Greenland ice sheet (GIS) and the biological activity associated with it are presented and then coupled with a surface mass-balance model to estimate total carbon fluxes due to respiration and photosynthesis. We find an average loading of 66gm−2 at Kangerlussuaq, southwest Greenland, which compares well with recent estimates from Kronprins Christians Land (17–440 gm−2: Bøggild and others, 2010) in northeast Greenland. We also report a significant microbial biomass in cryoconite at both these places (103–104 cells mg−1) and carbon fluxes of the order of 1–3 μM C g−1d−1 for both respiration and photosynthesis. The modelling indicates that total respiration and photosynthesis fluxes are likely to be ∼101–102 GgCa−1 and thus far from trivial. However, estimation of the net ecosystem impact across the entire ice sheet on atmospheric CO2 concentrations is problematic because photosynthesis rates were almost certainly low during our field campaign. Therefore, like its water balance, the carbon balance of the GIS is now known to be important, but its accurate quantification will remain elusive until more data are forthcoming.

Radar sounding is a powerful geophysical approach for characterizing the subsurface conditions of terrestrial and planetary ice masses at local to global scales. As a result, a wide array of orbital, airborne, ground-based, and in situ instruments, platforms and data analysis approaches for radioglaciology have been developed, applied or proposed. Terrestrially, airborne radar sounding has been used in glaciology to observe ice thickness, basal topography and englacial layers for five decades. More recently, radar sounding data have also been exploited to estimate the extent and configuration of subglacial water, the geometry of subglacial bedforms and the subglacial and englacial thermal states of ice sheets. Planetary radar sounders have observed, or are planned to observe, the subsurfaces and near-surfaces of Mars, Earth's Moon, comets and the icy moons of Jupiter. In this review paper, and the thematic issue of the Annals of Glaciology on ‘Five decades of radioglaciology’ to which it belongs, we present recent advances in the fields of radar systems, missions, signal processing, data analysis, modeling and scientific interpretation. Our review presents progress in these fields since the last radio-glaciological Annals of Glaciology issue of 2014, the context of their history and future prospects.

Air temperatures at the Marambio (MAR), Esperanza (ESP) and Matienzo (MAT) stations have been analyzed to investigate recent climate change on the eastern part of the Antarctic Peninsula. They are compared with data from the Oreadas station on the South Orkney Islands, the longest record available in Antarctica, and from the Faraday (FAR) station on the western coast of the Peninsula. Though the interannual variability is comparatively high and the stations are located in different climatic regimes, a pronounced warming trend shows up in all records. At MAR a temperature increase of 1.5°C has been observed since the beginning of the record in 1971. This is of similar magnitude to the increase at FAR on the west coast, which was 2.5°C for the longer period since 1945. The steady retreat and collapse of the northern Larsen Ice Shelf (LIS) coincided with this warming trend. of particular importance for the ice-shelf mass balance in this region are the summer temperatures which show a statistically significant warming trend at MAR and ESP. The representativity of the summer temperatures of MAR for northern LIS is confirmed by intercomparison with the parallel measurements at MAT which is located on the ice shelf.

The Tibetan Plateau interior area (TPIA), often termed the Qangtang Plateau, is distinguished by many dome-like mountains higher than 6000 ma.s.l. These mountains provide favourable conditions for the development of ice caps and glaciers of extreme continental/subpolar type. According to historical topographic maps (1959–80) and recent Landsat images (2004–11), continuous retreat was observed and the glacierized part of this area decreased by 9.5% (0.27% a–1) with respect to the total glacier area of 8036.4 km2 in the 1970s. Glaciers in the Zhari Namco basin have experienced the highest area shrinkage, with a reduction rate of 0.72% a–1, while the smallest reduction occurred in Bangong Co (0.12% a–1) and Dogai Coying basins (0.11% a–1). A regional gradient of area loss was found, with a larger decrease in the south and a smaller decrease in the north of the plateau. Comparisons indicate glaciers have experienced smaller shrinkage in the TPIA than in surrounding regions. Glacier shrinkage in the TPIA is mainly attributed to an increase in air temperature, while precipitation, glacier size and positive difference of glaciation also played an important role.

Flow stripes along the surface of rapidly moving ice streams are shown to be an expected reaction of a viscous medium flowing over an irregular bed, whenever the velocity at the bed is large compared to shearing through the thickness. The principal features of the process are as follows. At high basal speeds the ice acts as a strongly selective band-pass filter transmitting basal undulations on spatial scales of a few ice thicknesses very effectively to the surface. The decay of the short-scale features by diffusion is strongly retarded by horizontal stress gradients. Consequently, localized disturbances at the bed produce topographic effects on the surface that are advected long distances downstream before decaying away. This mechanism may explain many of the flow stripes on active ice streams moving by rapid basal motion.

Two spectral snow-reflectance models that account for the effects of grain-size and liquid-water fraction are described and initial validation results presented. The models are based upon the spectral complex refractive index of liquid water and ice in the region from 400 to 2500 nm. Mie scattering calculations are used to specify the essential optical properties of snow in the models. Two approaches are explored to model the effect of liquid water in the snow. The first accounts for the liquid water as separate spheres interspersed with ice spheres in the snow layer. The second accounts for the liquid water as coatings on ice grains in the snow layer. A discrete-ordinate radiative transfer code is used to model the spectral reflectance of the snow for the Mie-calculated optical properties. Both the interspersed- and coated-sphere models show that the snow-absorption feature at 1030 nm shifts to shorter wavelength as the liquid-water content increased. The expression of these shifts is different for the two models. A comparison of the models with a spectral measurement of frozen and melting snow shows better agreement with the coated-sphere model. A spectral fitting algorithm was developed and tested with the coated-sphere model to derive the grain-size and liquid-water fraction from snow spectral reflectance measurements. Consistent values of grain-size and liquid water were retrieved from the measured snow spectra. This research demonstrates the use of spectral models and spectral measurements to derive surface snow grain-size and liquid-water fraction. The results of this research may be extended to regional and greater scales using data acquired by airborne and spaceborne imaging spectrometers for contributions to energy balance and hydrological modeling.

The majority of glaciers in central Chile have receded in recent decades, from >50m to only a few meters per year, mainly in response to an increase in the 0˚C isotherm altitude. The Aconcagua river basin (33˚ S) is one of the major glaciated basins in central Chile, with 121 km2 of ice in 2003. An earlier inventory using 1955 aerial photographs yielded a total surface area of 151 km2, implying a reduction in glacier area of 20% (0.63km2 a–1) over the 48 years. Photographic stereo models, high-resolution satellite images (Landsat, ASTER) and SRTM data have been used to delineate glacier basins. A focus on Glaciar Juncal Norte, one of the largest glaciers in the basin, allows a more detailed analysis of changes. The glacier has exhibited a smaller reduction (14%) between 1955 and 2006, and the resulting elevation changes over this smaller period are not significant. The above reduction rates are lower than in other glaciers of central Chile and Argentina. This trend emphasizes water runoff availability in a river where most of the water in the dry summers is generated by glaciers and snowpack, and where most of the superficial water rights are already allocated. Ongoing hydrological research including modelling of future water runoff will improve our understanding.

There has been a marked decline in the summer extent of Arctic sea ice over the past few decades. Data from autonomous ice mass-balance buoys can enhance our understanding of this decline. These buoys monitor changes in snow deposition and ablation, ice growth, and ice surface and bottom melt. Results from the summer of 2008 showed considerable large-scale spatial variability in the amount of surface and bottom melt. Small amounts of melting were observed north of Greenland, while melting in the southern Beaufort Sea was quite large. Comparison of net solar heat input to the ice and heat required for surface ablation showed only modest correlation. However, there was a strong correlation between solar heat input to the ocean and bottom melting. As the ice concentration in the Beaufort Sea region decreased, there was an increase in solar heat to the ocean and an increase in bottom melting.

A simple quasi one-dimensional model of flowing avalanches is presented. It is a further development of that used in the Swiss Guidelines for practitioners. It is shown that shearing in avalanche movement is concentrated near the ground and that, due to the geometrical roughness of the ground, a flow resistance proportional to the square of velocity must be taken into account in addition to dry friction. For the change of flow on changing slope angles it is demonstrated that under certain conditions for internal friction a “normal” flow on a flat lower part can no longer be attained; the avalanche behaves like a rigid body. The runout distance is in fair agreement with the Guidelines if a larger internal friction is used. The main differences lie in much smaller deposition depths and smaller velocities during runout.

Different types of fast ice flow (both spatial and temporal) in valley glaciers (surging glaciers, tidewater glaciers and deforming-bed glaciers) and ice sheets (ice streams and deforming-bed ice-sheet flow) are discussed briefly. Although there are unlikely to be any specific individual landforms associated with fast ice flow, there may be landform assemblages.

At valley glacier scale, it is suggested that there are two landform assemblages: (1) an ice-thrust type, dominated by bulldozed push moraines and hummocky moraines (associated with glaciers with a high supraglacial sediment supply, a coarse-grained substrate and a coarse-grained proglacial sediment wedge); and (2) a bed-flow type dominated by “squeeze” push moraines, flutes and drumlins (associated with glaciers with a low supra- glacial sediment supply and fine-grained substrate). The ice-thrust type alone is only associated with discontinuous fast flow (on both rigid and deforming beds); whilst the bed-flow type is associated with both continuous and discontinuous fast flow.

It is suggested that these two landform assemblages may also be indicative of fast ice flow at ice-sheet scale, in particular the bed-flow style. If that is the case, then discontinuous fast ice flow may be indicated by the ice-thrust landform assemblage and the bed-flow style where drumlins are present.

It is also suggested that specific evidence for ice streams includes the distinctive land-form assemblages within valley or fan-like locations, and a predictable pattern of velocity reflected by drumlin elongation ratios.

An improved understanding of fresh water stored in the Himalaya is crucial for water resource management in South Asia and can be inferred from glacier mass-balance estimates. However, field investigations in the rugged Himalaya are limited to a few individual glaciers and short duration. Therefore, we have recently developed an approach that combines satellite-derived snowlines, a temperature-index melt model and the accumulation-area ratio method to estimate annual mass balance of glaciers at basin scale and for a long period. In this investigation, the mass balance of 146 glaciers in the Chandra basin, western Himalaya, is estimated from 1984 to 2012. We estimate the trend in equilibrium line altitude of the basin as +113 m decade−1 and the mean mass balance as −0.61 ± 0.46 m w.e. a−1. Our basin-wide mass-balance estimates are in agreement with the geodetic method during 1999–2012. Sensitivity analysis suggests that a 20% increase in precipitation can offset changes in mass balance for a 1 °C temperature rise. A water loss of 18% of the total basin volume is estimated, and 67% for small and low-altitude glaciers during 1984–2012, indicating a looming water scarcity crisis for villages in this valley.

Submarine melting at the ice–ocean interface is a significant term in the mass balance of marine-terminating outlet glaciers. However, obtaining direct measurements of the submarine melt rate, or the ocean heat transport towards the glacier that drives this melting, has been difficult due to the scarcity of observations, as well as the complexity of oceanic flows. Here we present a method that uses synoptic velocity and temperature profiles, but accounts for the dominant mode of velocity variability, to obtain representative heat transport estimates. We apply this method to the Sermilik Fjord–Helheim Glacier system in southeastern Greenland. Using lowered acoustic Doppler current profiler (LADCP) and hydrographic data collected in summer 2009, we find a mean heat transport towards the glacier of 29 × 109W, implying a submarine melt rate at the glacier face of 650 ma–1. The resulting adjusted velocity profile is indicative of a multilayer residual circulation, where the meltwater mixture flows out of the fjord at the surface and at the stratification maximum.

Between austral late winter 1993 and austral autumn 1998, during five cruises aboard the research vessel Nathaniel B. Palmer, almost 300 m of core was obtained from first-year ice floes in the Ross, Amundsen and Bellingshausen Seas. Analysis of the texture, stratigraphy and stable-isotopic composition of the ice was used to assess the magnitude of the role of flooding and snow-ice formation at the base of the snowpack in the thickening of the ice cover and the thinning of the snow cover. Snow ice occurred in all ice-thickness categories and made a significant contribution to the total ice mass (12−36%) in both autumn and winter. Although the amount of snow ice was often exceeded by the amount of frazil ice and congelation ice, the thickness of individual layers of each ice type indicated that snow ice often made a greater contribution to the thermodynamic thickening of the ice cover than the other ice types. The larger quantities of frazil ice and congelation ice were primarily the result of dynamic thickening. Flooding and snow-ice formation reduced the snow cover to 42−70% of the total snow accumulation depending on time and location. On the basis of this information, ship-based snow-depth estimates were adjusted to estimate the total snow accumulation on different ice-thickness categories.

Radar profiles of bed echo intensity can survey conditions at the ice–bed interface and test for the presence or absence of water. However, extracting information about basal conditions from bed echo intensities requires an estimate of the attenuation loss through the ice. We used the relationship between bed echo intensities from constant-offset radar data and ice thickness to estimate depth-averaged attenuation rates at several locations on and near Kamb Ice Stream (KIS), West Antarctica. We found values varying from 29 dBkm–1 at Siple Dome to 15 dBkm–1 in the main trunk region of KIS, in agreement with a previous measurement and models. Using these attenuation-rate values, we calculated the relative bed reflectivity throughout our KIS surveys and found that most of the bed in the trunk has high basal reflectivities, similar to those obtained in the location of boreholes that found water at the bed. Areas of lower bed reflectivity are limited to the sticky spot, where a borehole found a dry bed, and along the margins of KIS. These results support previous hypotheses that the basal conditions at locations like the sticky spot on KIS control its stagnation and possible reactivation.

We present a series of benchmark experiments designed for testing and comparing numerical ice-sheet models. Following the outcome of two EISMINT workshops organized to intercompare large-scale ice-sheet models currently in operation, model benchmark experiments ate described for ice sheets under fixed and moving margin conditions. These address both steady-state and time-dependent behaviour under schematic boundary conditions and with prescribed physics. A comparison was made of each model’s prediction of basic geophysical variables such as ice thickness, velocity and temperature. Consensus achieved in the model inter-comparison provides reference solutions against which the accuracy and consistency of ice-sheet modelling codes can be assessed.

Two- and three-dimensional avalanche dynamics models are being increasingly used in hazard-mitigation studies. These models can provide improved and more accurate results for hazard mapping than the simple one-dimensional models presently used in practice. However, two- and three-dimensional models generate an extensive amount of output data, making the interpretation of simulation results more difficult. To perform a simulation in three-dimensional terrain, numerical models require a digital elevation model, specification of avalanche release areas (spatial extent and volume), selection of solution methods, finding an adequate calculation resolution and, finally, the choice of friction parameters. In this paper, the importance and difficulty of correctly setting up and analysing the results of a numerical avalanche dynamics simulation is discussed. We apply the two-dimensional simulation program RAMMS to the 1968 extreme avalanche event In den Arelen. We show the effect of model input variations on simulation results and the dangers and complexities in their interpretation.

Using the finite-element code Elmer, we show that the full Stokes modeling of the ice-sheet/ice-shelf transition we propose can give consistent predictions of grounding-line migration. Like other marine ice-sheet models our approach is highly sensitive to the chosen mesh resolution. However, with a grid size down to <5 km in the vicinity of the grounding line, predictions start to be robust because: (1) whatever the grid size (<5 km) the steady-state grounding-line position is sensibly the same (6 km standard deviation), and (2) with a grid-size refinement in the vicinity of the grounding line (200 m), the steady-state solution is independent of the applied perturbation in fluidity, provided this perturbation remains monotonic.