We investigate variations in the surface area of glaciers in Sagarmatha national park, Nepal, during the second half of the 20th century through comparison of a map applicable to the late 1950s with the official map of Nepal in the early 1990s. The comparison reveals a slight overall decrease in glacier area (by 4.9%, from 403.9 to 384.6 km2), a result which, though potentially subject to errors arising from cartographic interpretation, is in line with the area reductions found by other studies of Asian glaciers. We find that the areas of some individual glaciers, the largest situated at higher altitudes, increased during the study period. This was most apparent for the glaciers oriented to the south, with the increase occurring mainly in the glacier accumulation zones while the fronts tended to recede. Meanwhile, the smaller glaciers, situated lower and on steep basins, experienced a reduction. For the smaller glaciers, the sections most affected by change were the accumulation zones, and these glaciers showed a tendency for the front to advance. In this region there is a lack of climate data for high altitudes. Nevertheless, observations from stations situated around the park suggest that, alongside temperature variations which are often considered the primary factor eliciting glacier response, changes in precipitation play a significant role.

Increased mass losses from the Greenland ice sheet and inferred contributions to sea-level rise have heightened the need for hydrologic observations of meltwater exiting the ice sheet. We explore whether temporal variations in ice-sheet surface hydrology can be linked to the development of a downstream sediment plume in Kangerlussuaq Fjord by comparing: (1) plume area and suspended sediment concentration from Moderate Resolution Imaging Spectroradiometer (MODIS) imagery and field data; (2) ice-sheet melt extent from Special Sensor Microwave/Imager (SSM/I) passive microwave data; and (3) supraglacial lake drainage events from MODIS. Results confirm that the origin of the sediment plume is meltwater release from the ice sheet. Interannual variations in plume area reflect interannual variations in surface melting. Plumes appear almost immediately with seasonal surface-melt onset, provided the estuary is free of landfast sea ice. A seasonal hysteresis between melt extent and plume area suggests late-season exhaustion in sediment supply. Analysis of plume sensitivity to supraglacial events is less conclusive, with 69% of melt pulses and 38% of lake drainage events triggering an increase in plume area. We conclude that remote sensing of sediment plume behavior offers a novel tool for detecting the presence, timing and interannual variability of meltwater release from the ice sheet.

Rutford Ice Stream and Carlson Inlet are neighbouring glaciers in West Antarctica. Rutford Ice Stream flows at speeds greater than 350 m a−1, whereas Carlson Inlet, which has some similar dimensions and supports a similar driving stress, flows 10–50 times slower. We discuss a range of observations concerning Carlson Inlet, and conclude that there is good indirect evidence that it is a relict ice stream, which ceased streaming more than 240 years BP, but sufficiently recently that its surface morphology, basal water content and basal morphology still retain characteristics produced by streaming. An analysis of expected subglacial drainage pathways indicates that Carlson Inlet is not streaming because it is currently starved of subglacial water, which is currently directed beneath Rutford Ice Stream. This current state of water piracy by Rutford Ice Stream is, however, sensitive to minor thickness changes on the ice streams; a ∼120 m (<4%) thickening of Rutford Ice Stream would divert almost all the subglacial water in the system towards Carlson Inlet and could reactivate its flow. The result highlights the importance of subglacial drainage in controlling ice-stream evolution and the requirement for ice-sheet models to couple ice flow with subglacial drainage.

Three debris-free glaciers with strongly differing annual glaciological glacier-wide mass balances (MBs) are monitored in the Everest region (central Himalaya, Nepal). The mass budget of Mera Glacier (5.1 km2 in 2012), located in the southern part of this region, was balanced during 2007–15, whereas Pokalde (0.1 km2 in 2011) and West Changri Nup glaciers (0.9 km2 in 2013), ~30 km further north, have been losing mass rapidly with annual glacier-wide MBs of −0.69 ± 0.28 m w.e. a−1 (2009–15) and −1.24 ± 0.27 m w.e. a−1 (2010–15), respectively. An analysis of high-elevation meteorological variables reveals that these glaciers are sensitive to precipitation, and to occasional severe cyclonic storms originating from the Bay of Bengal. We observe a negative horizontal gradient of annual precipitation in south-to-north direction across the range (≤−21 mm km−1, i.e. −2% km−1). This contrasted mass-balance pattern over rather short distances is related (i) to the low maximum elevation of Pokalde and West Changri Nup glaciers, resulting in years where their accumulation area ratio is reduced to zero and (ii) to a steeper vertical gradient of MB for glaciers located in the inner arid part of the range.

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.

In glacial environments particle-size analysis of moraines provides insights into clast origin, transport history, depositional mechanism and processes of reworking. Traditional methods for grain-size classification are labour-intensive, physically intrusive and are limited to patch-scale (1 m2) observation. We develop emerging, high-resolution ground- and unmanned aerial vehicle-based ‘Structure-from-Motion’ (UAV-SfM) photogrammetry to recover grain-size information across a moraine surface in the Heritage Range, Antarctica. SfM data products were benchmarked against equivalent datasets acquired using terrestrial laser scanning, and were found to be accurate to within 1.7 and 50 mm for patch- and site-scale modelling, respectively. Grain-size distributions were obtained through digital grain classification, or ‘photo-sieving’, of patch-scale SfM orthoimagery. Photo-sieved distributions were accurate to <2 mm compared to control distributions derived from dry-sieving. A relationship between patch-scale median grain size and the standard deviation of local surface elevations was applied to a site-scale UAV-SfM model to facilitate upscaling and the production of a spatially continuous map of the median grain size across a 0.3 km2 area of moraine. This highly automated workflow for site-scale sedimentological characterization eliminates much of the subjectivity associated with traditional methods and forms a sound basis for subsequent glaciological process interpretation and analysis.

Permeability and meltwater flow have been studied in sea ice in the Siberian and central Arctic during the summers of 1995 and 1996. A bail-test technique has been adapted to allow for measurements of in situ permeability, found to range between 10−11 and 10−8 m2. Permeability varied by about a factor of 2 between 1995 (above-normal melt rates) and 1996 (below-normal melt rates). Release of fluorescent tracers (fluorescein, rhodamine) furthermore allowed the derivation of flow velocities and assessment of the relevant driving forces. Hydraulic gradients in rough ice and wind stress in ponded ice were found to be particularly important, driving meltwater over distances of several meters per day. The mid- to late-summer ice was found to be permeable enough to completely divert meltwater from the surface into the ice interior. It is shown, however, that lower permeabilities of the upper ice layers as well as refreezing of meltwater, particularly during the early melt season, foster the development of surface melt ponds.

Despite the important role of supraglacial debris in ablation, knowledge of debris thickness on Himalayan glaciers is sparse. A recently developed method based on reanalysis data and thermal band satellite imagery has proved to be potentially suitable for debris thickness estimation without the need for detailed field data. In this study, we further develop the method and discuss possibilities and limitations arising from its application to a glacier in the Himalaya with scarce in situ data. Surface temperature patterns are consistent for 13 scenes of Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER) and Landsat 7 imagery and correlate well with incoming shortwave radiation and air temperature. We use an energy-balance approach to subtract these radiation or air temperature effects, in order to estimate debris thickness patterns as a function of surface temperature. Both incoming shortwave and longwave radiation are estimated with reasonable accuracy when applying parameterizations and reanalysis data. However, the model likely underestimates debris thickness, probably due to incorrect representation of vertical debris temperature profiles, the rate of heat storage and turbulent sensible heat flux. Moreover, the uncertainty of the result was found to increase significantly with thicker debris, a promising result since ablation is enhanced by thin debris of 1–2 cm.

Evolution of cold dry snow and firn plays important roles in glaciology; however, the physical formulation of a densification law is still an active research topic. We forced eight firn-densification models and one seasonal-snow model in six different experiments by imposing step changes in temperature and accumulation-rate boundary conditions; all of the boundary conditions were chosen to simulate firn densification in cold, dry environments. While the intended application of the participating models varies, they are describing the same physical system and should in principle yield the same solutions. The firn models all produce plausible depth-density profiles, but the model outputs in both steady state and transient modes differ for quantities that are of interest in ice core and altimetry research. These differences demonstrate that firn-densification models are incorrectly or incompletely representing physical processes. We quantitatively characterize the differences among the results from the various models. For example, we find depth-integrated porosity is unlikely to be inferred with confidence from a firn model to better than 2 m in steady state at a specific site with known accumulation rate and temperature. Firn Model Intercomparison Experiment can provide a benchmark of results for future models, provide a basis to quantify model uncertainties and guide future directions of firn-densification modeling.

We present newly acquired airborne radar data showing ice thickness and surface elevation for Pine Island Glacier, Antarctica. These data, when combined with earlier measurements, suggest the presence of a lightly grounded area immediately above the grounding line of Pine Island Glacier. We identify this region as an “ice plain”. It lies close to the centre line of the glacier, has an elevation above buoyancy of <50 m and extends inland for >28 km. The upstream edge of the ice plain is defined by a “coupling line”. The configuration of the ice plain implies that nearby thinning of the ice stream would result in substantial grounding-line retreat. We suggest that the grounding-line retreat of Pine Island Glacier, observed between 1992 and 1996, probably commenced sometime after 1981.

We investigate the elevation and mass-balance response of tributary glaciers following the loss of the Larsen A and B ice shelves, Antarctic Peninsula (in 1995 and 2002 respectively). Our study uses MODIS imagery to track ice extent, and ASTER and SPOT5 digital elevation models (DEMs) plus ATM and ICESat laser altimetry to track elevation changes, spanning the period 2001–09. The measured Larsen B tributary glaciers (Hektoria, Green, Evans, Punchbowl, Jorum and Crane) lost up to 160 m in elevation during 2001–06, and thinning continued into 2009. Elevation changes were small for the more southerly Flask and Leppard Glaciers, which are still constrained by a Larsen B ice shelf remnant. In the northern embayment, continued thinning of >3 m a−1 on Drygalski Glacier, 14 years after the Larsen A ice shelf disintegrated, suggests that mass losses for the exposed Larsen B tributaries will continue for years into the future. Grounded ice volume losses exceed 13 km3 for Crane Glacier and 30 km3 for the Hektoria–Green–Evans glaciers. The combined mean loss rate for 2001–06 is at least 11.2 Gt a−1. Our values differ significantly from published mass-budget-based estimates for these embayments, but are a reasonable fraction of GRACE-derived rates for the region (∼40 Gt a−1).

We propose a general framework for iceberg-calving models that can be applied to any calving margin. The framework is based on mass continuity, the assumption that calving rate and terminus velocity are not independent and the simple idea that terminus thickness following a calving event is larger than terminus thickness at the event onset. The theoretical, near steady-state analysis used to support and analyze the framework indicates that calving rate is governed, to first order, by ice thickness, thickness gradient, strain rate, mass-balance rate and backwards melting of the terminus; the analysis furthermore provides a physical explanation for a previously derived empirical relationship for ice-shelf calving (Alley and others, 2008). In the calving framework the pre- and post-calving terminus thicknesses are given by two unknown but related functions. The functions can vary independently of changes in glacier flow and geometry, and can therefore account for variations in calving behavior due to external forcings and/or self-sustaining calving processes (positive feedbacks). Although the calving framework does not constitute a complete calving model, any thickness-based calving criterion can easily be incorporated into the framework. The framework should be viewed as a guide for future attempts to parameterize calving.

We report on analysis of meteorological data for the period 27 May–20 August 2004, from two automatic weather stations on McCall Glacier, Alaska, USA, aimed at studying the relationship between climate and ablation. One station is located on a mountain ridge and the other in the ablation area where we also analyzed the energy balance. The weather station on the glacier measured an average temperature of 5.3°C (at 2 m height above surface) and wind speed of 3.1 m s−1 (at 3 m height). A sonic height ranger and ablation stakes indicate a specific mass balance of –1.94 ± 0.09 m w.e between 15 June and 20 August. The specific mass balance calculated from the surface energy balance, –2.06 ± 0.18 m w.e., is in close correspondence to this. The latter is the sum of 0.12 m w.e. of snowfall, 0.003 m w.e. of deposition and –2.18 m w.e. of melt. Net radiation contributes 74% of the melt energy. Compared to ablation measurements in the early 1970s, summer ablation was large. This increase is explained by a combination of a relatively higher net radiation, a lower albedo and larger turbulent heat fluxes that led to more energy being available for melting. No single meteorological variable can be isolated as being the principal reason for the high ablation, however. The lower ice albedo (0.19) is possibly due to ash deposits from forest fires.

The length of time TM over which a glacier responds to a prior change in climate is investigated with reference to the linearized theory of kinematic waves and to results from numerical models. We show the following: TM may in general be estimated by a volume time-scale describing the time required for a step change in mass balance to supply the volume difference between the initial and final steady states. The factor f in the classical estimate of τM = ƒl/u, where I is glacier length and u is terminus velocity, has a simple geometrical interpretation. Ft is the ratio of thickness change averaged over the full length I to the change at the terminus. Although both u and f relate to dynamic processes local to the terminus zone, the ratio f/u and, therefore, Tm are insensitive to details of the terminus dynamics, in contrast to conclusions derived from some simplified kinematic wave models. A more robust estimate of Tm independent of terminus dynamics is given by TM= h/(–b) where h is a thickness scale for the glacier and –b is the mass-balance rate (negative) at the terminus. We suggest that Tm for mountain glaciers can be substantially less than the 1O2–103 years commonly considered to be theoretically expected.

The genesis of some types of patterned ground, including hummocks, frost boils and sorted stone circles, has been attributed to differential frost heave (DFH). However, a theoretical model that adequately describes DFH has yet to be developed and validated. In this paper, we present a mathematical model for the initiation of DFH, and discuss how variations in physical (i.e. soil/vegetation properties) and environmental (i.e. ground/air temperatures) properties affect its occurrence and length scale. Using the Fowler and Krantz multidimensional frost-heave equations, a linear stability analysis anda quasi-steady-state real-time analysis are performed. Results indicate that the following conditions positively affect the spontaneous initiation of DFH: silty soil, small Young’s modulus, small non-uniform surface heat transfer or cold uniform surface temperatures, and small freezing depths. The initiating mechanism for DFH is multidimensional heat transfer within the freezing soil. Numerical integration of the linear growth rates indicates that expression of surface patterns can become evident on the 10–100 year time-scale.

We use optical (ASTER and Landsat) and radar (ERS-1 and ERS-2) satellite imagery to document changes in the Prince Gustav Ice Shelf, Antarctic Peninsula, and its tributary glaciers before and after its January 1995 collapse. The satellite image record captures the transition from an ice-shelf glacier system to a tidewater glacial system and the subsequent rapid retreat and inferred ‘fatal’ negative mass balances that occur as lower glacier elevations lead to higher ablation and tidewater-style calving collapse. Pre-1995 images show that the central ice shelf was fed primarily by Sjögren Glacier flowing from the Antarctic Peninsula and by Röhss Glacier flowing from James Ross Island. Numerous structural discontinuities (rifts and crevasses) and melt ponds were present on the ice shelf before the collapse. After the ice shelf collapsed, Röhss Glacier retreated rapidly, becoming a tidewater glacier in 2002 and receding a total of ∼15 km between January 2001 and March 2009, losing >70% of its area. Topographic profiles of Röhss Glacier from ASTER-derived digital elevation models show a thinning of up to ∼150 m, and surface speeds increased up to ninefold (0.1–0.9 m d−1) over the same period. The rates of speed increase and elevation loss, however, are not monotonic; both rates slowed between late 2002 and 2005, accelerated in 2006 and slowed again in 2008–09. We conclude that tributary glaciers react to ice-shelf removal by rapid (if discontinuous) recession, and that the response of tidewater glaciers on the Antarctic Peninsula to ice-shelf removal occurs over timescales ranging from sub-annual to decadal.

A new 10 year surface mass balance (SMB) record of Hurd and Johnsons Glaciers, Livingston Island, Antarctica, is presented and compared with earlier estimates on the basis of local and regional meteorological conditions and trends. Since Johnsons is a tidewater glacier, we also include a calving flux calculation to estimate its total mass balance. The average annual SMB over the 10 year observation period 2002–11 is −0.15 ± 0.10 m w.e. for Hurd Glacier and 0.05 ± 0.10 m w.e. for Johnsons Glacier. Adding the calving losses to the latter results in a total mass balance of −0.09 ± 0.10 m w.e. There has been a deceleration of the mass losses of these glaciers from 1957–2000 to 2002–11, which have nearly halved for both glaciers. We attribute this decrease in the mass losses to a combination of increased accumulation in the region and decreased melt. The increased accumulation is attributed to larger precipitation associated with the recent deepening of the circumpolar pressure trough, while the melt decrease is associated with lower summer surface temperatures during the past decade.

Knowledge of iceberg calving is important for understanding instabilities of tidewater glaciers and ice sheets. Since 1995 the terminus of LeConte Glacier, Alaska, U.S.A., has retreated about 2 km and the glacier has thinned approximately 120 m at its 1999 terminus position. Our focus is short-term (hours to weeks) variability of the frequency and magnitude of calving events and calving flux. Both photogrammetric and visual observations are employed in a temporal analysis over a several-week period. We combined these data with measurements of ice speed, tide level, surface water input and water-storage estimates in an attempt to better understand the calving process. Contrary to results obtained over longer time-scales on other glaciers, our results show no correlation between ice speed and the frequency of calving. However, calving events do not appear to occur randomly; often they are a response to measurable changes in other parameters within the terminus region. Caclving can often be attributed to buoyancy perturbations and possibly flexure of the nearly floating terminus. Given the multiple possibilities for buoyancy perturbations, we have found no simple relationship between any specific forcing parameter and calving at short time-scales.

Structural glaciological, sedimentological and geophysical techniques are used to provide new insight concerning the formation of band ogives and associated structures at Bas Glacier d’Arolla, Switzerland. Sedimentary stratification, crevasse traces and transverse foliation are identified as planar structures in the lower icefall and glacier tongue. Stratification and crevasse traces are progressively deformed into, and enhance, the transverse foliation found in the glacier tongue. Three-dimensional geometry has been defined using ground-penetrating radar, which portrays four main characteristics: (i) deep reflectors interpreted as the ice/bed interface, (ii) alternating reflection-rich and reflection-poor zones interpreted as ogives, (iii) up-glacier-dipping reflectors, interpreted as planar structures, and (iv) down-glacier-dipping reflectors of uncertain origin. At the glacier surface, each band ogive consists of a light and dark band. The dark bands contain more intense foliation which, on differential weathering, traps fine debris. Clasts and clear ice of basal character within dark ogive bands suggest that basal ice has been raised to the glacier surface. The most applicable model for the formation of band ogives at Bas Glacier d’Arolla is a refinement of Posamentier’s (1978) “reverse faulting” hypothesis. In this context, multiple shear zones are formed, through which basal ice is uplifted to the glacier surface to produce the dark, foliated ogive bands. This model fits observations reported from other glaciers with band ogives.