The sensitivity of glaciers to climatic change is key information in assessing the response and sea-level implications of projected future warming. New Zealand glaciers are important globally as an example of how maritime glaciers will contribute to sea-level rise. A spatially distributed energy-balance model is applied to Brewster Glacier, New Zealand, in order to calculate glacier mass balance, run-off and sensitivity to climate change. The model successfully simulates four annual mass-balance cycles. Close to half (52%) of the energy available for melt on the glacier is supplied by turbulent heat fluxes, with radiation less important, except during the winter. Model sensitivity to temperature change is one of the largest reported on Earth, at −2.0 m w.e. a−1 °C−1. In contrast, a 50% change in precipitation is required to offset the mass-balance change resulting from a 1 °C temperature change. Meltwater runoff sensitivity is also very high, increasing 60% with a 1°C warming. The extreme sensitivity of mass balance to temperature change suggests that significant ice loss will occur with even moderate climate warming.

At a site on the ice sheet adjacent to the Jakobshavn ice stream in West Greenland, ice deformation rates and temperatures have been measured in boreholes to the bedrock at 830 m depth. Enhanced deformation rates were recorded just below the Holocene–Wisconsin transition at 680 m depth. A 31 m layer of temperate ice and the temperature minimum of −22°C at 520 m depth were detected. The good agreement of these data with results of a two-dimensional thermomechanically coupled flow model implies that the model input is adequate. Discrepancies between modelled and measured temperature profiles on a flowline at the ice-stream centre have been attributed to effects not accounted for by the model. We have suggested that the convergent three-dimensional flow leads to a vertical extension of the basal ice entering the stream. A thick basal layer of temperate and Wisconsin ice would explain the fast flow of this ice stream. As a test of this hypothesis, the new core-borehole conductivity (CBC) method has been used to compare conductivity sequences from the ice stream to those of the adjacent ice sheet. The correlation thus inferred suggests that the lowest 270 m of the ice sheet correspond to the lowermost 1700 m of the stream, and, consequently, that the lower part of the ice stream has experienced a very large vertical extension.

When a mass balance is computed using an outdated map, that computation does not reveal the actual mass change. But older maps often must be used for practical reasons. We present a method by which, with a few additional measurements each year, a mass balance computed with an outdated map can be transformed into an actual mass change. This is done by taking into account the influence of changes in areal extent and changes in the surface elevation of the glacier since the map was made. This method is applied to South Cascade Glacier, Washington, U.S.A., as an example. The computed cumulative mass balance from 1970 to 1997 would have been 16% too negative if the 1970 map had not been updated. While the actual volume change of a glacier is relevant to hydrological studies, the change that would have occurred on a constant (or static) surface is more relevant to certain glacier dynamics problems and most climate problems. We term this the reference-surface balance and propose that such a balance, which deliberately omits the influence of changes in area and surface elevation, is better correlated to climatic variations than the conventional one, which incorporates those influences.

Mass changes of the Antarctic ice sheet impact sea-level rise as climate changes, but recent rates have been uncertain. Ice, Cloud and land Elevation Satellite (ICESat) data (2003–08) show mass gains from snow accumulation exceeded discharge losses by 82 ± 25 Gt a−1, reducing global sea-level rise by 0.23 mm a−1. European Remote-sensing Satellite (ERS) data (1992–2001) give a similar gain of 112 61 Gt a−1. Gains of 136 Gt a−1 in East Antarctica (EA) and 72 Gt a−1 in four drainage systems (WA2) in West Antarctic (WA) exceed losses of 97 Gt a−1 from three coastal drainage systems (WA1) and 29 Gt a−1 from the Antarctic Peninsula (AP). EA dynamic thickening of 147 Gt a−1 is a continuing response to increased accumulation (>50%) since the early Holocene. Recent accumulation loss of 11 Gt a−1 in EA indicates thickening is not from contemporaneous snowfall increases. Similarly, the WA2 gain is mainly (60 Gt a−1) dynamic thickening. In WA1 and the AP, increased losses of 66 ± 16 Gt a−1 from increased dynamic thinning from accelerating glaciers are 50% offset by greater WA snowfall. The decadal increase in dynamic thinning in WA1 and the AP is approximately one-third of the long-term dynamic thickening in EA and WA2, which should buffer additional dynamic thinning for decades.

The principal engineering problems requiring knowledge concerning the properties of ice are prediction of the maximum force that ice might exert on a structure and determination of the load that can be placed safely on an ice cover. The properties of ice relevant to these problems are discussed. Particular attention is given to the dependence of Young's modulus on the frequency of the stress, strain-rate dependence of the yield and fracture strengths, and the ductile to brittle transition. It is shown that the strain-rate dependence of these properties is consistent with information on the stress dependence of dislocation velocity on the basal plane. The tensile and compressive strengths are also shown to be consistent with current theories of crack initiation and propagation.

Most theories of glacier movement and subglacial erosion have assumed that glaciers rest on rigid bedrock surfaces. Whilst this is probably correct for much of the bed area of most modern glaciers, deformable sediments do occur beneath them and formed a substantial area of the beds of large ice sheets during glacial periods. Observations and theories are presented and reviewed about the processes of glacier erosion of rock and unlithified sediment beds both when they are frozen and unfrozen.

Erosional bedrock landforms, such as roches moutonnées, indicate two principal subglacial erosional processes, plucking and abrasion. Where supraglacially derived debris is unimportant, plucking provides the tools which abrade the bed, and must be a quantitatively more important process than abrasion, though more localized. Where plucking is suppressed, erosion rates must be slow. Subglacial measurements of abrasion rates beneath a temperate glacier are used to test an earlier abrasional theory (Boulton, [C1974]). The form of the predicted abrasion-rate curve for changing ice velocity and pressure is verified. This theory successfully simulates two-dimensional erosional bedforms. Subglacial observations demonstrate how flow basal ice around the flanks of bedrock obstacles causes streaming of debris to occur. It is suggested that this streaming process is primarily responsible for the longitudinally lineated form of large-scale surfaces typical of glacially eroded bedrock.

Plucking and abrasion also occur beneath cold ice, though at slow rates, and are probably restricted to places where the ice thickness is small.

Where the glacier bed is composed of unlithified sediment, subglacial measurements show that deformation can produce very large discharges of subglacial material, which makes this a potential agent of very rapid subglacial landform production. The heterogeneity of subglacial sediment leads to spatially variable rates of deformation, and it is suggested that relatively stronger parts of the sediment body may form the nuclei for drumlin and mega-flute formation.

Whereas unlithified unfrozen sediment deforms beneath the glacier rather than being incorporated within it, ice-cemented subglacial sediments can behave like bedrock, because of their relative rigidity, and are readily plucked and incorporated englacially. They may also deform beneath the glacier.

Glacier surfaces support unique microbial food webs dominated by organic and inorganic debris called ‘cryoconite’. Observations from Longyearbreen, Spitsbergen, show how these aggregate particles can develop an internal structure following the cementation of mineral grains (mostly quartz and dolomite) by filamentous microorganisms. Measurements of carbon and dissolved O2 show that these microorganisms, mostly cyanobacteria, promote significant rates of photosynthesis (average 17 μgC g−1 d−1) which assist aggregate growth by increasing the biomass and producing glue-like extracellular polymeric substances. The primary production takes place not only upon the surface of the aggregates but also just beneath, due to the translucence of the quartz particles. However, since total photosynthesis is matched by respiration (average 19 μgC g−1 d−1), primary production does not contribute directly to cryoconite accumulation upon the glacier surface. The microorganisms therefore influence the surface albedo most by cementing dark particles and organic debris together, rather than simply growing over it. Time-lapse photographs show that cryoconite is likely to reside upon the glacier for years as a result of this aggregation. These observations therefore show that a better understanding of the relationship between supraglacial debris and ablation upon glaciers requires an appreciation of the biological processes that take place during summer.

The field observations of G.S. Boulton stimulated widespread interest in deformable beds. Shear resistance of till in its critical state is insensitive to strain rate and increases linearly with effective pressure. During unsteady deformation, pseudo-viscous shear resistance can be caused by dilation of consolidated tills and resultant pore-pressure decline. This effect is probably uncommon, however, because susceptible tills of low hydraulic diffusivity are also those least likely to consolidate significantly during effective-pressure transients. Stick-slip motion at Whillans Ice Stream, Antarctica, indicates that its basal till must weaken during rapid slip and strengthen during longer periods of slower slip. Recurrence intervals for rapid-slip episodes there (6-18hours) indicate that till-strength variations, if driven by changes in pore pressure either related or unrelated to basal freezing, are focused in the uppermost several centimeters of the bed. Ploughing of grains at the bed surface and associated excess pore pressures in adjacent till can account for rate-weakening during rapid slip, with pore-pressure decay causing strengthening between slip episodes. By promoting shallow, sluggish subglacial water flow and low effective pressure, soft beds may help sustain themselves by slowing their own transport. Soft-bed shear resistance, kinematics and continuity are problems rooted in subglacial hydrology.

An understanding of the reflection and transmission of light by snow is important for snow thermodynamics, hydrology, ecology and remote sensing. Snow has an intricate microstructure replete with ice/air interfaces that scatter light. Spectral observations of light reflection and transmission, from 400 to 1000 nm, were made in temperate snowpacks, under cold and under melting conditions. The optical observations were made using a dual-detector spectroradiometer. One detector was placed above the snow surface to monitor the incident and reflected solar irradiance, and the second detector was placed at the base of snow cover to measure downwelling irradiance. The optical measurements were supplemented by a physical characterization of the snow, including depth, density and an estimate of grain size. In general, transmitted light levels were low and showed a strong spectral dependence, with maximum values between 450 and 550 nm. For example, a 10 cm thick snow layer reduced visible transmission (500 nm) to about 5% of the incident irradiance, and infrared transmission (800 nm) to less than 1%. Extinction coefficients were in the range 3–30 m−1, and tended to decrease slightly as the snow aged and increase as snow density increased.

The behavior of marine-terminating ice sheets, such as the West Antarctic ice sheet, is of interest due to the possibility of rapid grounding-line retreat and consequent catastrophic loss of ice. Critical to modeling this behavior is a choice of basal rheology, where the most popular approach is to relate the ice-sheet velocity to a power-law function of basal stress. Recent experiments, however, suggest that near-grounding line tills exhibit Coulomb friction behavior. Here we address how Coulomb conditions modify ice-sheet profiles and stability criteria. The basal rheology necessarily transitions to Coulomb friction near the grounding line, due to low effective stresses, leading to changes in ice-sheet properties within a narrow boundary layer. Ice-sheet profiles ‘taper off’ towards a flatter upper surface, compared with the power-law case, and basal stresses vanish at the grounding line, consistent with observations. In the Coulomb case, the grounding-line ice flux also depends more strongly on flotation ice thickness, which implies that ice sheets are more sensitive to climate perturbations. Furthermore, with Coulomb friction, the ice sheet grounds stably in shallower water than with a power-law rheology. This implies that smaller perturbations are required to push the grounding line into regions of negative bed slope, where it would become unstable. These results have important implications for ice-sheet stability in a warming climate.

Discharge measurements, climate observations and hydrochemical samples gathered monthly (1998/99) in the Yanamarey and Uruashraju glacier-fed catchments of the Cordillera Blanca, Peru, permit an analysis of the glacier meltwater contribution to stream-flow. These glacier catchments feed the Río Santa, which discharges into the Pacific Ocean. Based on a water-balance computation, glacier melt contributes an estimated 35% of the average discharge from the catchments. For comparison, a volumetric end-member mixing model of oxygen isotopes shows glacier melt contributes 30–45% to the total annual discharge. Based on stream geochemistry, discharge from the Yanamarey glacier catchment provides 30% of the annual volume discharged from the Querococha watershed, which is <10% glacierized. By analogy, the larger Río Santa watershed, also <10% glacierized, receives at least 12% of its annual discharge from melting glacier ice. Tributary watersheds to the Río Santa with larger fractions of glacier cover have less variable runoff and enhanced discharge, demonstrating that the glaciers effectively buffer stream discharge seasonally. With continued glacier melting, stream-flow will likely become more variable, and there will be less dry-season runoff.

The ablation areas of debris-covered glaciers typically consist of a complex mosaic of surface features with contrasting processes and rates of mass loss. This greatly complicates glacier response to climate change, and increases the uncertainty of predictive models. In this paper we present a series of high-resolution DEMs and repeat lake bathymetric surveys on Ngozumpa Glacier, Nepal, to study processes and patterns of mass loss on a Himalayan debris-covered glacier in unprecedented detail. Most mass loss occurs by melt below supraglacial debris, and melt and calving of ice cliffs (backwasting). Although ice cliffs cover only ~5% of the area of the lower tongue, they account for 40% of the ablation. The surface debris layer is subject to frequent re-distribution by slope processes, resulting in large spatial and temporal differences in debris-layer thickness, enhancing or inhibiting local ablation rates and encouraging continuous topographic inversion. A moraine-dammed lake on the lower glacier tongue (Spillway Lake) underwent a period of rapid expansion from 2001 to 2009, but later experienced a reduction of area and volume as a result of lake level lowering and sediment redistribution. Rapid lake growth will likely resume in the near future, and may eventually become up to 7 km long.

Using a slightly modified form of the Spring–Hutter equations, glacial outburst floods are simulated from three classic sites, “Hazard Lake”, Yukon, Canada, Summit Lake, British Columbia, Canada, and Grímsvötn, Iceland, in order to calibrate the hydraulic roughness associated with subglacial conduits. Previous work has suggested that the Manning roughness of the conduits is remarkably high, but the new calibration yields substantially lower values that are representative of those for natural streams and rivers. The discrepancy can be traced to a poor assumption about the effectiveness of heat transfer at the conduit walls. The simulations reveal behaviour that cannot be inferred from simplified theories: (1) During flood onset, water pressure over much of the conduit can exceed the confining pressure of surrounding ice. (2) Local values of fluid potential gradient can differ substantially from the value averaged over the length of the conduit, contradicting an assumption of simple theories. (3) As the flood progresses, the location of flow constrictions that effectively control the flood magnitude can jump rapidly over large distances. (4) Predicted water temperature at the conduit outlet exceeds that suggested by measurements of exit water temperature.

Since the 1970s, the sudden, rapid collapse of 20% of ice shelves on the Antarctic Peninsula has led to large-scale thinning and acceleration of its tributary glaciers. The leading hypothesis for the collapse of most of these ice shelves is the process of hydrofracturing, whereby a water-filled crevasse is opened by the hydrostatic pressure acting at the crevasse tip. This process has been linked to observed atmospheric warming through the increased supply of meltwater. Importantly, the low-density firn layer near the ice-shelf surface, providing a porous medium in which meltwater can percolate and refreeze, has to be filled in with refrozen meltwater first, before hydrofracturing can occur at all. Here we build upon this notion of firn air depletion as a precursor of ice-shelf collapse, by using a firn model to show that pore space was depleted in the firn layer on former ice shelves, which enabled their collapse due to hydrofracturing. Two climate scenario runs with the same model indicate that during the 21st century most Antarctic Peninsula ice shelves, and some minor ice shelves elsewhere, are more likely to become susceptible to collapse following firn air depletion. If warming continues into the 22nd century, similar depletion will become widespread on ice shelves around East Antarctica. Our model further suggests that a projected increase in snowfall will protect the Ross and Filchner–Ronne Ice Shelves from hydrofracturing in the coming two centuries.

Models are proposed for channelized and distributed flow of meltwater at the base of an ice sheet. The volumes of both channel and distributed systems evolve according to a competition between processes that open drainage space (e.g. sliding over bedrock, melting of the ice) and processes that close it (e.g. viscous creep of the ice due to a positive effective pressure). Channels are generally predicted to have lower water pressure and therefore capture water from the surrounding regions of distributed flow. There is a natural length scale associated with the distributed system that determines the width of the bed from which water can be drawn into a channel. It is suggested that this determines the spacing between major channels and that this may be reflected in the spacing of eskers. A more permeable distributed system results in more widely spaced, and therefore larger, channels. Calculations of the flow into the head of a channel reveal that there is a critical discharge necessary for it to form, and provide a criterion for where channels can exist.

While numerous maps of Greenland ice flow velocity exist, most have gaps in coverage and/or accuracy is limited. We processed a large volume of synthetic aperture radar and Landsat 8 imagery collected between 1995 and 2015 to produce a nearly complete map of ice flow velocity for Greenland at a far greater accuracy than most prior products. We evaluated the accuracy of this map by comparing it with a variety of measured and estimated velocities. For the slow-moving interior of the ice sheet, where estimates are determined from interferometric phase, the errors are ~2 m a−1 or better. For coastal areas, where estimates are determined entirely from speckle- or feature-tracking methods, errors are 2–3 m a−1, which is in good agreement with the estimated formal errors. Especially for the slow-moving majority of the ice sheet, this map provides an important source of data for numerous types of glaciological studies.

We have used ground-based radio-echo sounding (RES) profiles to reveal the spatial distribution of basal and internal ice properties across Siple Dome, West Antarctica, and under the dormant ice streams on its flanks. The RES-detected bed-reflection power, corrected for the effects of instrumentation and ice-thickness variation, is nearly constant across Siple Dome at a value suggesting spatially homogeneous basal properties of ice frozen to bedrock. Till, if present under the dome, must be thin (<0.1 m). The high basal reflectivity measured under now dormant “Siple Ice Stream” (SIS) and Ice Stream C suggests that they are underlain by either a thin (<0.05 m) water layer or a thick (>1 m) thawed or frozen till layer. The evidence that the dormant SIS is not frozen directly to underlying bedrock (but is separated by a water or till layer) is a further indication that it was once an active ice stream, and suggests that streaming motion may have ceased before the basal layer was frozen. The absence of a thick till layer beneath Siple Dome is consistent with its apparent stability as an inter-ice-stream ridge in the past and may suggest that it will remain as a stable limitation of ice-stream width in the future.

Observations show that the Greenland ice sheet has been losing mass at an increasing rate over the past few decades, which makes it a major contributor to sea-level rise. Here we use a three-dimensional higher-order ice-flow model, adaptive mesh refinement and inverse methods to accurately reproduce the present-day ice flow of the Greenland ice sheet. We investigate the effect of the ice thermal regime on (1) basal sliding inversion and (2) projections over the next 100 years. We show that steady-state temperatures based on present-day conditions allow a reasonable representation of the thermal regime and that both basal conditions and century-scale projections are weakly sensitive to small changes in the initial temperature field, compared with changes in atmospheric conditions or basal sliding. We conclude that although more englacial temperature measurements should be acquired to validate the models, and a better estimation of geothermal heat flux is needed, it is reasonable to use steady-state temperature profiles for short-term projections, as external forcings remain the main drivers of the changes occurring in Greenland.

Equations are developed that can be used to determine the amount of gas present in sea ice from measurements of the bulk ice density, salinity, and temperature in the temperature range of −2 to −30 °C. Conversely these relationships can be used to give the density of sea ice as a function of its temperature and salinity, considering both the presence of gas and of solid salts in the ice. Equations are also given that allow the calculation of the gas and brine volumes in the ice at temperatures other than that at which the bulk density was determined.

Flood and debris-flow hazards at Grubengletscher near Saas Balen in the Saas Valley, Valais, Swiss Alps, result from the formation and growth of several lakes at the glacier margin and within the surrounding permafrost. In order to prevent damage related to such hazards, systematic investigations were carried out and practical measures taken. The evolution of the polythermal glacier, the creeping permafrost within the large adjacent rock glacier and the development of the various periglacial lakes were monitored and documented for the last 25 years by photogrammetric analysis of annually flown high-resolution aerial photographs. Seismic refraction, d.c. resistivity and gravimetry soundings were performed together with hydrological tracer experiments to determine the structure and stability of a moraine dam at a proglacial lake. The results indicate a maximum moraine thickness of > 100 m; extremely high porosity and even ground caverns near the surface may have resulted from degradation of sub- and periglacial permafrost following 19th/20th-century retreat of the partially cold glacier tongue. The safety and retention capacity of the proglacial lake were enhanced by deepening and reinforcing the outlet structure on top of the moraine complex. The water level of an ice-dammed lake was lowered and a thermokarst lake artficially drained.