We investigate the transferability of an enhanced temperature-index melt model that was developed and tested on Haut Glacier d’Arolla, Switzerland, in the 2001 season. The model’s empirical parameters (temperature factor, TF, and shortwave radiation factor, SRF) are recalibrated for: (1) other locations on Haut Glacier d’Arolla; (2) subperiods of distinct meteorological conditions; (3) different years on Haut Glacier d’Arolla; and (4) other glaciers in different years. The model parameters are optimized against simulations of an energy-balance model validated against ablation observations. Results are compared with those obtained with the original parameters. The model works very well when applied to other sites, seasons and glaciers, with the exception of overcast conditions. Differences are due to underestimation of high melt rates. The parameter values are associated with the prevailing energy-balance conditions, showing that high SRF are obtained on clear-sky days, whereas higher TF are typical of locations where glacier winds prevail and turbulent fluxes are high. We also provide a range of parameters clearly associated with the site’s location and its meteorological characteristics that could help to assign parameter values to sites where few data are available.

This study independently quantifies geodetic elevation change and models surface mass balance to solve the continuity equation. The approach is tested on two dynamically different glaciers, Kongsvegen and Kronebreen in northwest Svalbard, through two time epochs, 1966-1990/95 (I) and 1990/95-2007 (II). On Kongsvegen, a dynamically inactive glacier, the residual term represents an error associated with determining elevation changes and surface mass balance. It is apparent that centerline mass-balance estimates are not representative of the entire glacier, which we relate to center-line accumulation being larger than the elevation bin average. On Kronebreen, a fast-flowing and actively calving glacier, a significant part of the residual is identified with the long-term calving flux. For both glaciers, the cumulative surface mass balance remained close to zero during the first epoch but became increasingly negative in the second epoch. The long-term calving flux of Kronebreen is estimated to be -0.14 ± 0.03 km3 w.e.a-1 during epoch I and-0.20 ± 0.05 km3 w.e.a-1 in epoch II.

We describe an instrument developed to make in situ measurements of salinity and solid- fraction profiles in growing sea ice. The vertical resolution of the measurements is up to a few millimeters, with a temporal resolution of up to fractions of a second. The technique is based on impedance measurements between platinum wires around which sea ice grows. Data obtained using this instrument in laboratory experiments are in good agreement with theoretical predictions. In a field test in the Arctic, the bulk salinity of growing sea ice has been measured in situ throughout the whole depth of the ice layer. The data are compared with bulk salinities obtained from ice cores, and confirm the general understanding that the bulk salinity in ice-core studies is significantly underestimated in the lower parts of the cores. The approach can also be used in other glaciological applications and for general studies of two-phase, two-component porous media.

Jakobshavn Isbræ, West Greenland, underwent a large, rapid and not well understood change in flow dynamics in 1998, leading to a doubling of its ice discharge rates. We calculate the width-averaged forces controlling flow of Jakobshavn Isbræ in 1995, 2000 and 2005 to elucidate processes responsible for this change in flow speed. In contrast to earlier suggestions, we conclude that the observed acceleration was not caused by the loss of back-stress due to weakening and subsequent break-up of the floating ice tongue alone. Gradients in longitudinal stress are small at all times considered (∼3% of the driving stress) and basal and lateral drag provide resistance to flow. Over the 10 year period considered, the average driving stress increased by 20 kPa, which was balanced by a comparable increase in lateral drag. We surmise that the velocity changes resulted from weakening of the ice in the lateral shear margins and perhaps a change in properties at the bed. Possible mechanisms for weakening of ice in the lateral shear margins include cryo-hydrologic warming of subsurface ice in the ablation zone and hydraulic weakening due to higher water content of ice in the shear margins.

In mountainous regions, snow transport due to wind significantly influences snow distribution and, as a result, avalanche danger. A physically based numerical two-layer model is developed to simulate blowing and drifting snow in Alpine terrain. One layer describes the driving-wind field and the transport in suspension. The description is based on the atmospheric boundary-layer equations, using an e − ∊ model for the turbulent closure. The second layer describes the transport due to saltation, including erosion and deposition of snow. Here, conservation equations for mass and momentum are formulated for the mixture of snow and air. Particle trajectory calculations are used to parameterize quantities characterizing the saltation layer. Both layers are mutually coupled by boundary conditions. A two-way coupling between particles and airflow is taken into account. Comparisons between simulation results and field measurements around an Alpine crest show encouraging results.

Analysis of the cross-flow transmission of force from the central parts of a well-lubricated ice stream to its margins shows that there is a corresponding shift in the lateral location of motion-induced heat generation. The rate of basal heat generation in the center can be substantially smaller than the local rate of potential energy loss given by driving stress times the speed of downslope motion. The basal heating is a maximum for an intermediate level of lubrication for which speed is about 40% of the speed over a friction-less bed and base stress is about 25% of the driving stress. Stable and unstable balances between meltwater production and drainage on the bed are identified. A stable steady state with a speed less (more) than that giving maximum heat generation is termed drainage-(production-) limited, since an increase in speed would lead to increased (decreased) basal melting and must (need not) be balanced by increased drainage. It is shown that gradual evolution of the basal water drainage system and the factors affecting basal melting can cause discontinuous jumps between fast- and slow-moving states. A simplified analysis applied to six cross-sections of West Antarctic Ice Streams B, D, E and Rutford Ice Stream shows them to be diverse in the level of support from the sides and corresponding shift of mechanical heating sideward from their central parts. The cross-sections of Ice Stream B near “Upstream B” may be production-limited, because of especially high lubrication and related support from the sides. Cross-sections in the upper part of Ice Stream D, Ice Stream E and Rutford Ice Stream are in a drainage-limited condition. Substantial reduction of basal heat generation by side drag (in most cases) and expected high heat flow into the basal ice associated with low thickness (in some cases) tends to favor basal freezing. Nevertheless, all of the examined cross-sections except one are expected to experience basal melting with a modest geothermal heat-flux density of 60 m W m−1 or less in some cases. The lower part of Ice Stream B is an exception, where the analysis indicates that geothermal flux density must exceed 80–100 m W−1 m to maintain melting. If this high geothermal flux is not present, then the base of the lower part of Ice Stream B may be freezing, which would suggest continued deceleration of this part of Ice Stream B.

We use an ice-flow model to demonstrate how flow variations initiated in the marginal zone of an ice sheet affect flow farther inland through longitudinal (along-flow) coupling. Our findings allow for an alternate interpretation of seasonal accelerations observed near the equilibrium line of the Greenland ice sheet (Zwally and others, 2002). We demonstrate that these observations can be explained by accelerations initiated up to 12 km closer to the margin where the ice is ∼40% thinner, is heavily crevassed, experiences a seasonal doubling of velocity, and where the ablation rate, surface meltwater flux and ice temperature are likely higher. Our modeling and observations suggest that conditions and processes normally found near ice-sheet margins are adequate for explaining the observations of Zwally and others (2002). This and considerations of the likely subglacial hydrology in the marginal zone lead us to suggest that seasonal accelerations may have limited impact on ice-sheet mass balance even in the face of climate warming.

We provide an assessment of the supraglacial water budget of a moulin basin on the western margin of the Greenland ice sheet for 15 days in August 2009. Meltwater production, the dominant input term to the 1.14 ± 0.06 km2 basin, was determined from in situ ablation measurements. The dominant water-output terms from the basin, accounting for 52% and 48% of output, respectively, were moulin discharge and drainage into crevasses. Moulin discharge exhibits large diurnal variability (0.017–0.54 m3 s−1) with a distinct late-afternoon peak at 16:45 local time. This lags peak meltwater production by ∼2.8 ± 4.2 hours. An Extreme Ice Survey time-lapse photography sequence complements the observations of moulin discharge. We infer, from in situ observations of moulin geometry, previously published borehole water heights and estimates of the temporal lag between meltwater production and observed local ice surface uplift (‘jacking’), that the transfer of surface meltwater to the englacial water table via moulins is nearly instantaneous (<30 min). We employ a simple crevasse mass-balance model to demonstrate that crevasse drainage could significantly dampen the surface meltwater fluctuations reaching the englacial system in comparison to moulin discharge. Thus, unlike crevasses, moulins propagate meltwater pulses to the englacial system that are capable of overwhelming subglacial transmission capacity, resulting in enhanced basal sliding.

The surface velocity field of Devon Ice Cap, Nunavut, Canada, was mapped using interferometric synthetic aperture radar (InSAR). Ascending European Remote-sensing Satellite 1 and 2 (ERS-1/-2) tandem mode data were used for the western and southeast sectors, and 3 day repeat pass ERS-1 imagery for the northeast sector. Speckle-tracking procedures were used with RADARSAT 1 imagery to obtain surface velocities over the terminus of Belcher Glacier (a major calving front) where decorrelation between ERS data occurred. The InSAR data highlight a significant contrast in ice-flow dynamics between the east and west sides of the ice cap. Ice movement west of the main north–south divide is dominated by relatively uniform ‘sheet’ flow, but three fast-flowing outlet glaciers that extend 14–23km beyond the ice-cap margin also drain this region. Several outlet glaciers that extend up to 60 km inland from the eastern margin drain the eastern side of the ice cap. The dominant ice-flow regimes were classified based on the relationship between the driving stress (averaged over a length scale of ten ice thicknesses) and the ratio of surface velocity to ice thickness. The mapped distribution of flow regimes appears to depict the spatial extent of basal sliding across the ice cap. This is supported by a close relationship between the occurrence of flow stripes on the ice surface and flow regimes where basal sliding was found to be an important component of the glacier motion. Iceberg calving rates were computed using measured surface velocities and ice thicknesses derived from airborne radio-echo sounding. The volume of ice calved between 1960 and 1999 was estimated to be 20.5 ± 4.7 km3 (or 0.57 km3 a–1). Approximately 89% of this loss occurred along the eastern margin. The largest single source is Belcher Glacier, which accounts for ~50% of the total amount of ice calved.

Supraglacial ponds play a key role in absorbing atmospheric energy and directing it to the ice of debris-covered glaciers, but the spatial and temporal distribution of these features is not well documented. We analyse 172 Landsat TM/ETM+ scenes for the period 1999–2013 to identify thawed supraglacial ponds for the debris-covered tongues of five glaciers in the Langtang Valley of Nepal. We apply an advanced atmospheric correction routine (Landcor/6S) and use band ratio and image morphological techniques to identify ponds and validate our results with 2.5 m Cartosat-1 observations. We then characterize the spatial, seasonal and interannual patterns of ponds. We find high variability in pond incidence between glaciers (May–October means of 0.08–1.69% of debris area), with ponds most frequent in zones of low surface gradient and velocity. The ponds show pronounced seasonality, appearing in the pre-monsoon as snow melts, peaking at the monsoon onset at 2% of debris-covered area, then declining in the post-monsoon as ponds drain or freeze. Ponds are highly recurrent and persistent, with 40.5% of pond locations occurring for multiple years. Rather than a trend in pond cover over the study period, we find high interannual variability for each glacier after controlling for seasonality.

The deformation of the ground surface that is produced by frost heave has motivated almost a century of concerted laboratory, field and theoretical studies. Well before the development of equipment capable of resolving the microscopic films that support liquid transport towards growing ice lenses, early investigators predicted their occurrence and noted their importance. Idealized experiments continue to prompt theoretical advances that have been combined to develop predictive models for the macroscopic frost-heave characteristics that are seen in the field. This contribution highlights steps on the road to our current understanding of the physical interactions that control how ice forms in unconsolidated porous media.

During an ice-tank experiment, samples were taken to study the processes of acquisition and alteration of the gas properties in young first-year sea ice during a complete growth–warming–cooling cycle. The goal was to obtain reference levels for total gas content and concentrations of atmospheric gases (O2, N2, CO2) in the absence of significant biological activity. The range of total gas-content values obtained (3.5–18 mL STP kg−1) was similar to previous measurements or estimates. However, major differences occurred between current and quiet basins, showing the role of the water dynamics at the ice–water interface in controlling bubble nucleation processes. Extremely high CO2 concentrations were observed in all the experiments (up to 57% in volume parts). It is argued that these could have resulted from two unexpected biases in the experimental settings. Concentrations in bubbles nucleated at the interface are controlled by diffusion both from the ice–water interface towards the well-mixed reservoir and between the interface water and the bubble itself. This double kinetic effect results in a transition of the gas composition in the bubbles from values close to solubility in sea water toward values close to atmospheric, as the ice cover builds up.

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.

Spaceborne interferometric synthetic aperture radar (InSAR) techniques for measuring ice flow velocity and topography have developed rapidly over the last decade and a half, revolutionizing the study of ice dynamics. Spaceborne interferometry has contributed to major progress in many areas of glaciological study by: providing the first comprehensive measurements of ice-stream flow velocity over the major outlets of Greenland and Antarctica; revealing that ice-stream and outlet-glacier flow can change rapidly (months to years); improving understanding of several ice-sheet and ice-shelf processes; providing velocity for flux-gate based mass-balance assessment; mapping flow of mountain glaciers; and capturing the geomorphic traces of past ice flow. We review the basic technique development, the measurement characteristics, and the extensive set of results yielded by these measurements.

Annual-layer thickness data, spanning AD 1534–2001, from an ice core from East Rongbuk Col on Qomolangma (Mount Everest, Himalaya) yield an age–depth profile that deviates systematically from a constant accumulation-rate analytical model. The profile clearly shows that the mean accumulation rate has changed every 50–100 years. A numerical model was developed to determine the magnitude of these multi-decadal-scale rates. The model was used to obtain a time series of annual accumulation. The mean annual accumulation rate decreased from ∼0.8 m ice equivalent in the 1500s to ∼0.3 m in the mid-1800s. From ∼1880 to ∼1970 the rate increased. However, it has decreased since ∼1970. Comparison with six other records from the Himalaya and the Tibetan Plateau shows that the changes in accumulation in East Rongbuk Col are broadly consistent with a regional pattern over much of the Plateau. This suggests that there may be an overarching mechanism controlling precipitation and mass balance over this area. However, a record from Dasuopu, only 125 km northwest of Qomolangma and 700 m higher than East Rongbuk Col, shows a maximum in accumulation during the 1800s, a time during which the East Rongbuk Col and Tibetan Plateau ice-core and tree-ring records show a minimum. This asynchroneity may be due to altitudinal or seasonal differences in monsoon versus westerly moisture sources or complex mountain meteorology.

By subtracting surface topographies from 1979, 1994, 2000 and 2008, we measured icethinning rates increasing from 1 ma-1 (1979-94) to >4 ma-1 (2000-08) on the tongue of Mer de Glace, French Alps. The relative contributions of changes in surface mass balance and ice fluxes to this acceleration in the thinning are estimated using field and remote-sensing measurements. Between 1979-94 and 2000-08, surface mass balance diminished by 1.2mw.e.a-1, mainly because of atmospheric warming. Mass-balance changes induced by the growing debris-covered area and the evolving glacier hypsometry compensated each other. Meanwhile, Mer de Glace slowed down and the ice fluxes through two cross sections at 2200 and 2050ma.s.l. decreased by 60%. Between 1979-94 and 2000-08, two-thirds of the increase in the thinning rates was caused by reduced ice fluxes and one- third by rising surface ablation. However, these numbers need to be interpreted cautiously given our inability to respect mass conservation below our upper cross section. An important implication is that large errors would occur if one term of the continuity equation (e.g. surface mass balance) were deduced from the two others (e.g. elevation and ice-flux changes).

Glaciers respond to long-term climate changes and also to the year-to-year fluctuations inherent in a constant climate. Differentiating between these factors is critical for the correct interpretation of past glacier fluctuations and for the correct attribution of current changes. Previous work has established that century-scale, kilometre-scale fluctuations can occur in a constant climate. This study asks two further questions of practical significance: how likely is an excursion of a given magnitude in a given amount of time, and how large a trend in length is statistically significant? A linear model permits analytical answers wherein the dependencies on glacier geometry and climate setting can be clearly understood. The expressions are validated with a flowline glacier model. The likelihood of glacier excursions is well characterized by extreme-value statistics, although probabilities are acutely sensitive to some poorly known glacier properties. Conventional statistical tests can be used for establishing the significance of an observed glacier trend. However, it is important to determine the independent information in the observations which can be effectively estimated from the glacier geometry. Finally, the retreat of glaciers around Mount Baker, Washington State, USA, is consistent with, but not independent proof of, the regional climate warming that is established from the instrumental record.

Observed meteorological data and a high-resolution (5 km) model were used to simulate Greenland ice sheet surface melt extent and trends before the satellite era (1960–79) and during the satellite era through 2010°. The model output was compared with passive microwave satellite observations of melt extent. For 1960–2010 the average simulated melt extent was 15 ± 5%. For the period 1960–72, simulated melt extent decreased by an average of 6%, whereas 1973–2010 had an average increase of 13%, with record melt extent in 2010. The trend in simulated melt extent since 1972 indicated that the melt extent in 2010 averaged twice that in the early 1970s. The maximum and mean melt extents for 2010 were 52% (∼9.5 × 105 km2) and 28% (∼5.2 × 105 km2), respectively, due to higher-than-average winter and summer temperatures and lower-than-average winter precipitation. For 2010, the southwest Greenland melt duration was 41–60 days longer than the 1960–2010 average, while the northeast Greenland melt duration was up to 20 days shorter. From 1960 to 1972 the melting period (with a >10% melt extent) decreased by an average of 3 days a−1. After 1972, the period increased by an average of 2 days a−1, indicating an extended melting period for the ice sheet of about 70 days: 40 and 30 days in spring and autumn, respectively.

The effective viscosity of ice depends upon many factors, including temperature, deviatoric stress, crystal orientation and impurities. A flow law that includes these factors and is simple to implement is a requirement for numerically efficient ice-flow models. The dominant microscale flow mechanism changes as temperature, deviatoric stress or grain-size changes. For both anisotropic and isotropic constitutive relations, this shift in dominant flow mechanism is expressed as a change in the stress exponent. We study the effects of this shift in stress exponent on ice flow using a two-term flow law for isotropic ice. Our stress–strain-rate relationship does not explicitly describe the microscale processes of ice deformation; however, it encompasses a range of deformation behaviors with a simple law. In terrestrial ice, a flow-mechanism shift may occur in low-deviatoric-stress regions near ice divides, resulting in a near-linear constitutive relationship for ice flow. Compared to a non-linear (Glen) divide, a divide dominated by a near-linear flow mechanism has vertical-velocity profiles that are similar at divide and flank sites, internal layers that do not develop a Raymond bump, and a steady-state surface profile that is more rounded near the divide.

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.