Changes in ice-sheet surface elevation are caused by a combination of ice-dynamic imbalance, ablation, temporal variations in accumulation rate, firn compaction and underlying bedrock motion. Thus, deriving the rate of ice-sheet mass change from measured surface elevation change requires information on the rate of firn compaction and bedrock motion, which do not involve changes in mass, and requires an appropriate firn density to associate with elevation changes induced by recent accumulation rate variability. We use a 25 year record of surface temperature and a parameterization for accumulation change as a function of temperature to drive a firn compaction model. We apply this formulation to ICESat measurements of surface elevation change at three locations on the Greenland ice sheet in order to separate the accumulation-driven changes from the ice-dynamic/ablation-driven changes, and thus to derive the corresponding mass change. Our calculated densities for the accumulation-driven changes range from 410 to 610 kgm–3, which along with 900 kgm–3 for the dynamic/ablation-driven changes gives average densities ranging from 680 to 790 kgm–3. We show that using an average (or ‘effective’) density to convert elevation change to mass change is not valid where the accumulation and the dynamic elevation changes are of opposite sign.

We analyze the global sea-level budget since 1850. Good estimates of sea-level contributions from glaciers and small ice caps, the Greenland ice sheet and thermosteric sea level are available over this period, though considerable scope for controversy remains in all. Attempting to close the sea-level budget by adding the components results in a residual displaying a likely significant trend of ~0.37mma–1 from 1955 to 2005, which can, however, be reasonably closed using estimated melting from unsurveyed high-latitude small glaciers and ice caps. The sea-level budget from 1850 is estimated using modeled thermosteric sea level and inferences from a small number of mountain glaciers. This longer-term budget has a residual component that displays a rising trend likely associated with the end of the Little Ice Age, with much decadal-scale variability that is probably associated with variability in the global water cycle, ENSO and long-term volcanic impacts.

In situ measurements of water depth and optical properties of a melt pond in East Greenland were collected to verify a previously developed algorithm for determining supraglacial water depths using satellite imagery. That algorithm made five simplifying assumptions which we have tested using the in situ data and laboratory analysis of water samples. We conclude that three assumptions were justified, one was not and the remaining one (substrate homogeneity) requires further study, but probably has a minor effect on the retrieved water depths and volumes. Measured water depths of 0.2–3.0m agree well with those derived from a satellite image using the algorithm. Numerically modeled depths also agree well with those from the satellite image. This new analysis demonstrates the validity of our algorithm as a means for determining meltwater volumes in supraglacial ponds and lakes.

In this paper we review the interaction of El Niño Southern Oscillation (ENSO) variability and warming trends recorded in ice-core records from high-altitude tropical glaciers, discuss the implications of the warming trends for the glaciers and consider the societal implications of glacier retreat. ENSO has strong impacts on meteorological phenomena that directly or indirectly affect most regions on the planet and their populations. Many tropical ice fields have provided continuous annually resolved proxy records of climatic and environmental variability preserved in measurable parameters, especially oxygen and hydrogen isotopic ratios (δ18O, δD) and the net mass balance (accumulation). These records present an opportunity to examine the nature of tropical climate variability in greater detail and to extract new information on linkages between rising temperatures on tropical glaciers and equatorial Pacific sea surface temperatures in critical ENSO indicator regions. The long-term climate records from a collection of high-altitude tropical ice cores provide the longer-term context essential for assessing the significance of the magnitude and rate of current climate changes that are in large measure driving glacier retreat. The well-documented ice loss on Quelccaya in the Peruvian Andes, Naimona’nyi in the Himalaya, Kilimanjaro in eastern Africa and the ice fields near Puncak Jaya in Papua, Indonesia, presents a grim future for low-latitude glaciers. The ongoing melting of these ice fields (response) is consistent with model predictions for a vertical amplification of temperature in the tropics (driver) and has serious implications for the people who live in these areas.

Within the percolation and soaked facies of the Greenland ice sheet, the relationship between radar-derived internal reflection horizons and the layered structure of the firn column is unclear. We conducted two small-scale ground-penetrating radar (GPR) surveys in conjunction with 10 m firn cores that we collected within the percolation and soaked facies of the Greenland ice sheet. The two surveys were separated by a distance of ~50 km and ~340m of elevation leading to ~40 days of difference in the duration of average annual melt. At the higher site (~1997ma.s.l.), which receives less melt, we found that internal reflection horizons identified in GPR data were largely laterally continuous over the grid; however, stratigraphic layers identified in cores could not be traced between cores over any distance from 1.5 to 14.0 m. Thus, we found no correlation between firn core stratigraphy observed directly and radar-derived internal reflection horizons. At the lower site (~1660ma.s.l.), which receives more melt, we found massive ice layers >0.5m thick and stratigraphic boundaries that span >15m horizontally. Some ice layers and stratigraphic boundaries correlate well with internal reflection horizons that are laterally continuous over the area of the radar grid. Internal reflection horizons identified at ~1997ma.s.l. are likely annual isochrones, but the reflection horizons identified at ~1660ma.s.l. are likely multi-annual features. We find that mapping accumulation rates over long distances by tying core stratigraphy to radar horizons may lead to ambiguous results because: (1) there is no stratigraphic correlation between firn cores at the 1997 m location; and (2) the reflection horizons at the 1660m location are multi-annual features.

A full-stress, thermomechanically coupled, numerical model is used to explore the interaction between basal thermal conditions and motion of a terrestrially terminating section of the west Greenland ice sheet. The model domain is a two-dimensional flowline profile extending from the ice divide to the margin. We use data-assimilation techniques based on the adjoint model in order to optimize the basal traction field, minimizing the difference between modeled and observed surface velocities. We monitor the sensitivity of the frozen/melted boundary (FMB) to changes in prescribed geothermal heat flux and sliding speed by applying perturbations to each of these parameters. The FMB shows sensitivity to the prescribed geothermal heat flux below an upper threshold where a maximum portion of the bed is already melted. The position of the FMB is insensitive to perturbations applied to the basal traction field. This insensitivity is due to the short distances over which longitudinal stresses act in an ice sheet.

Potential energy released from the capsize of ice-shelf fragments (icebergs) is the immediate driver of the brief explosive phase of ice-shelf disintegration along the Antarctic Peninsula (e.g. the Larsen A, Larsen B and Wilkins ice shelves). The majority of this energy powers the rapidly expanding plume of ice-shelf fragments that expands outward into the open ocean; a smaller fraction of this energy goes into surface gravity waves and other dynamic interactions between ice and water that can sustain the continued fragmentation and break-up of the original ice shelf. As an initial approach to the investigation of ice-shelf fragment capsize in ice-shelf collapse, we develop a simple conceptual model involving ideal rectangular icebergs, initially in unstable or metastable orientations, which are assembled into a tightly packed mass that subsequently disassembles via massed capsize. Computations based on this conceptual model display phenomenological similarity to aspects of real ice-shelf collapse. A promising result of the conceptual model presented here is a description of how iceberg aspect ratio and its statistical variance, the two parameters related to ice-shelf fracture patterns, influence the enabling conditions to be satisfied by slow-acting processes (e.g. environmentally driven melting) that facilitate ice-shelf disintegration.

Assessing the significance of current glacier loss on Kilimanjaro, Tanzania, demands a well-constrained temporal perspective. That context is provided by direct measurements, ancillary observations of the ice fields and the analyses of the ice cores collected from them. Ice retreat mechanisms observed there today are consistent with the preservation of the oldest ice, ~11.7 ka, in the central deepest part of the Northern Ice Field (NIF). This ice-core derived paleoclimate history published by Thompson and others (2002) is further confirmed by more recent paleoclimate records from tropical East Africa. Mounting evidence suggests that the (anticipated) loss of the entire NIF will be unprecedented within the past 10 000 years. New evidence bears directly on the mechanisms driving the current ice loss. Measurements made in 2000 on the NIF document that air temperature at 0.5 and 1.5 m above the surface remained below 5°C, while a surface temperature of 0.0°C was sustained for up to 8 hours d-1 under clear conditions, consistent with observations of melting on all Kilimanjaro summit ice fields. The linear relationship between oxygen and hydrogen isotopic ratios for all six ice cores drilled in 2000 lies very close to the global meteoric waterline and does not support sublimation (evaporation) as a major driver of ice loss today or in the past on Kilimanjaro.

A complete glacier inventory of the Himalaya and Karakoram (H-K) has been created by merging records from the Chinese Glacier Inventory, several regional inventories produced by the International Centre for Integrated Mountain Development, Kathmandu, Nepal, and partial inventories from the Geological Survey of India. The only remaining gap, the Indian part of Kashmir, has been filled by a reconnaissance inventory based on Soviet military maps at 1 : 200 000 scale representing the late 1970s. It contains records for 3526 glaciers covering 9584 km2. The new H-K inventory contains records and outlines for 20 812 glaciers covering 43 178 km2. The extent of ice in the Karakoram is slightly less than in the Himalaya, but the Karakoram glaciers are on average twice as thick (~160m as against ~80 m). A glacier-by-glacier analysis, relying on estimates of mass balance for the entire mountain range and on an extension of the often-used volume–area scaling relation, suggests that up to about one-fifth of the glaciers present in 1985 may have disappeared already. If mass loss were to remain constant at the average rate for 1975–2008, from 3000 to 13 000 more glaciers might disappear by 2035. If mass loss were to continue to accelerate as inferred for 1985–2008, only a few thousand to a few hundred glaciers might remain in 2035. Total area and total mass would each decrease by about one-half (constant-rate assumption) or three-quarters (constant-trend assumption). These projections, which are uncertain and neglect some possibly important mitigating controls, such as variable extents of debris cover and the feedback due to retreat to higher elevations, demonstrate the need for more complete analyses to inform public perceptions of, and policy decisions relating to, the health of H-K glaciers.

Ice surface altimetry from ICESat-1 and NASA aircraft altimeter overflights spanning 2002–09 indicate that a region of lower Crane Glacier, Antarctic Peninsula, shows an unusual temporal pattern of elevation loss: a period of very rapid drawdown (~91ma–1 between September 2004 and September 2005) bounded by periods of large but more moderate rates (23ma–1 until September 2004; 12ma–1 after September 2005). The region of increased drawdown is ~4.5 km ×2.2 km based on satellite (ASTER and SPOT-5) stereo-image digital elevation model (DEM) differencing spanning the event. In a later differential DEM the anomalous drawdown feature is not seen. Bathymetry in Crane Glacier fjord reveals a series of flat-lying, formerly subglacial deeps interpreted as lake sediment basins. We conclude that the elevation-change feature resulted from drainage of a small, deep subglacial lake. We infer that the drainage event was induced by hydraulic forcing of subglacial water past a downstream obstruction. However, only a fraction of Crane Glacier’s increase in flow speed that occurred near the time of lake drainage (derived from image feature tracking) appears to be directly attributable to the event; instead, retreat of the ice front off a subglacial ridge 6 km downstream of the lake is likely the dominant cause of renewed fast flow and more negative mass balance in the subsequent 4 years.

In this study we aim to connect glacier extent in 1998 with general climatic conditions, and glacier changes between 1969 and 1998 with climate change in the Austrian Alps. The investigations are based on two complete glacier inventories, a homogenized gridded precipitation dataset and European Centre for Medium-Range Weather Forecasts re-analysis (ERA-40) data of air temperature at different levels. A relationship between median glacier elevation, minimal elevation, the general elevation of the surrounding mountains and mean climatic values was found. In the Austrian Alps, the existence of glaciers at low elevations can only be maintained with above-average accumulation or strong dynamic ice supply. For debris-free glaciers, we found a limit of ~2080ma.s.l., where mean summer temperatures (June–August) exceed 4°C. Glacier changes from 1969 to 1998 are strongly negative both in relative area and in mean thickness. There is a weak and regionally varying negative trend in precipitation over this period. A spatially consistent sequence of positive temperature anomalies in the early 1980s and after 1990 offers an explanation for the retreat. The study shows that the observed spatial variability of glacier changes is connected more strongly to the topographic differences than to a regionally different climate change signal.

Area changes at 39 of the widest Greenland marine-terminating glacier outlets are measured in consecutive annual end-of-melt-season Moderate Resolution Imaging Spectroradiometer (MODIS) scenes spanning ten annual intervals (2000–10). The rates of cumulative area change for glaciers and ice shelves are well represented by linear least-squares fits, R = –0.99 and R = –0.94, with average rates of –70 and –65km2 a–1, respectively. Collectively, during this decade, the 39 glaciers lost a cumulative area of 1368 km2. More than three-quarters of the total area change occurred north of 72˚N. The largest 11-year area change for a single glacier during the survey period is the 311km2 loss at Humboldt Glacier. The largest annual change for a single glacier was extreme compared with the others, where Petermann glacier retreated 17km between 3 and 5 August 2010. For the 10 year sample, on average, the count of glaciers retreating is twice that advancing. A larger distinction is evident considering area change, with the ratio of retreat and advance, on average, nine times the gain. For glaciers with ice shelves, we find no year with collective area gain. The area change data from this study are posted at: http://bprc.osu.edu/~jbox/data/GAC/

Surface roughness enhances the net ablation rate associated with direct solar radiation relative to smooth surfaces, because roughness allows solar energy reflected from one part of the surface to be absorbed by another part. In this study we examine the feedback between solar-radiation-driven ablation and growth of surface roughness on the Greenland ice sheet, using a numerical model of radiative transfer. Our experiments extend previous work by examining: (1) the effects of diurnal and seasonal variation of solar zenith angle and azimuth relative to incipient roughness features, (2) the evolution of roughness geometry in response to radiatively driven ablation and (3) the relative solar energy collection efficiencies of various roughness geometries and geographic locations and orientations. A notable result of this examination is that the time evolution of the aspect ratio of surface features under solar-driven ablation collapses onto a roughly universal curve that depends only on latitude, not the detailed shape of the feature. The total enhancement of surface melt relative to a smooth surface over a full ablation season varies with this ratio, and this dependence suggests a way to parameterize roughness effects in large-scale models that cannot treat individual roughness features. Overall, our model results suggest that surface roughness at the latitudes spanned by the Greenland ice sheet tends to dissipate as the ablation season progresses.

Investigations into glacial changes, including understanding variations in the rates of glacial volume and surface-elevation changes, have increased over the past decade. This study uses historical glacier elevation data in the form of topographic maps from 1970 and a digital elevation model from the mid-1980s to calculate surface-elevation and volume changes for Ha-Iltzuk Icefield, southwest British Columbia, Canada. Results indicate that the icefield thinned at an average rate of 0.76±0.25 ma–1 during this period. A previous study of Ha-Iltzuk Icefield also using the geodetic method found a thinning rate of 1.0±0.20ma–1 between the mid-1980s and 1999, indicating a slight increase in the amount of icefield thinning. Within the ablation zone, thinning increased with decreasing elevation at a rate of 1.9±0.68 ma–1 km–1 between these two periods (1970 to mid-1980s versus mid-1980s to 1999). Analysis of meteorological data suggests that increases in both temperature and rainfall, as well as decreases in snowfall, likely contributed to the increased thinning rate.

The East Antarctic ice sheet (EAIS) is Earth’s largest reservoir of fresh water and has the potential to raise sea level by ~50 m. A significant amount of the ice sheet’s mass is discharged by outlet glaciers draining through the Transantarctic Mountains, the balance characteristics of which are largely unknown. Here the mass balance is estimated for four glaciers draining ice from the EAIS through the Transantarctic Mountains into the Ross Sea embayment: David, Mulock, Byrd and Nimrod glaciers. Remote-sensing observations are used to map changes in ice flow and surface elevation, and ultimately to compute the mass balance of each glacier using the input–output method and three separate estimates for accumulation rate. Results computed using this method indicate small positive balances for David (2.41±1.31 Gt a–1), Mulock (1.91±0.84 Gt a–1) and Nimrod (0.88±0.39 Gt a–1) glaciers, and a large positive imbalance for Byrd Glacier (21.67±4.04 Gt a–1). This large imbalance for Byrd Glacier is inconsistent with other observations, and is likely due to an overestimation of accumulation rates across large regions of the interior catchment.

Annual glacier mass balances are reconstructed for 96% of the Austrian glacier-covered area (451 of 470 km2) between 1969 and 1998. The volume change derived from two complete glacier inventories (1969 and 1998) serves as the boundary condition that is aimed to be reproduced. ERA-40 reanalysis data as well as a gridded precipitation dataset (HISTALP) are used to drive a positive degree-day (PDD) model. The results are verified with four independent long-term mass-balance series. The spatial and vertical distribution of the tuning parameters is altered in order to reproduce the measured mean annual surface mass balances of selected glaciers, and a strong correlation is found between the median elevation of a glacier and the degree-day factor (DDF) at this elevation. This result implies that the lower a glacier’s median elevation is, the less melt occurs at a given elevation and temperature. We attribute this to the fact that lower-altitude glaciers are generally those with more accumulation, which leads to later exposure of bare ice and a longer period of high-albedo snow cover. A further improvement of the model was achieved by making DDF a function of time as well as space. The results indicate that mean DDFs generally increase for a given date over a sequence of consecutive negative mass-balance years, which probably reflects the reduction in albedo related to that. Finally, the major drivers of the observed mass-balance evolution are investigated: summer PDD sums correlate significantly better with the observed mass-balance changes than annual PDD sums or precipitation do. This implies that annual mass balances in the study area are governed by summer temperatures.

Glacier inventories provide the baseline data to perform climate-change impact assessment on a regional scale in a consistent and spatially representative manner. In particular, a more accurate calculation of the current and future contribution to global sea-level rise from heavily glacierized regions such as Alaska is much needed. We present a new glacier inventory for a large part of western Alaska (including Kenai Peninsula and the Tordrillo, Chigmit and Chugach mountains), derived from nine Landsat Thematic Mapper scenes acquired between 2005 and 2009 using well-established automated glacier-mapping techniques (band ratio). Because many glaciers are covered by optically thick debris or volcanic ash and partly calve intowater, outlineswere manually edited in these wrongly classified regions during post-processing. In total we mapped ~8830 glaciers (>0.02 km2) with a total area of ~16 250 km2. Large parts of the area (47%) are covered by a few (31) large (>100 km2) glaciers, while glaciers less than 1 km2 constitute only 7.5% of the total area but 86% of the total number. We found a strong dependence of mean glacier elevation on distance from the ocean and only a weak one on aspect. Glacier area changes were calculated for a subset of 347 selected glaciers by comparison with the Digital Line Graph outlines from the US Geological Survey. The overall shrinkage was ~23% between 1948–57 and 2005–09.

Meltwater from glaciers in the European Alps plays an important role in hydropower production, and future glacier development is thus of economic interest. However, an up-to-date and alpine-wide inventory for accurate assessment of glacier changes or modelling of future glacier development has not hitherto been available. Here we present a new alpine-wide inventory (covering Austria, France, Italy and Switzerland) derived from ten Landsat Thematic Mapper (TM) scenes acquired within 7 weeks in 2003. Combined with the globally available digital elevation model from the Shuttle Radar Topography Mission, topographic inventory parameters were derived for each of the 3770 mapped glaciers, covering 2050 km2. The area-class frequency distribution is very similar in all countries, and a mean northerly aspect (NW, N, NE) is clearly favoured (arithmetic counting). Mean glacier elevation is ~2900 m, with a small dependence on aspect. The total area loss since the previous glacier inventory (acquired around 1970±15 years) is roughly one-third, yielding a current area loss rate of ~2%a–1. Digital overlay of the outlines from the latest Austrian glacier inventory revealed differences in the interpretation of glacier extents that prohibit change assessment. A comparison of TM-derived outlines with manually digitized extents on a high-resolution IKONOS image returned 1.5% smaller glaciers with TM.

Pronounced changes in glacier mass and length were observed for the monitored glaciers in the Jostedalsbreen region, Norway, since the last glacier inventories were compiled in the 1960s and 1980s. However, the current overall extent of the glaciers in the region is not well known. To obtain this information, we have compiled a new inventory from two mosaicked Landsat Thematic Mapper (TM) scenes acquired in 2006 that have excellent snow conditions for glacier mapping, the first suitable scenes for this purpose after 22 years of imaging with TM. Drainage divides and topographic inventory parameters were derived from a 25 m national digital elevation model for 1450 glaciers. By digitizing glacier outlines from 1 : 50 000 scale topographic maps of 1966, we calculated changes in glacier area for ~300 glaciers. Cumulative length changes for the 1997–2006 period were derived from an additional TM scene and compared with field measurements for nine glaciers. Overall, we find a 9% area loss since 1966, with a clear dependence on glacier size; however, seasonal snow in 1966 in some regions made area determination challenging. The satellite-derived length changes confirmed the observed high spatial variability and were in good agreement with field data (±1 pixel), apart from glacier tongues in cast shadow. The new inventory will be freely available from the Global Land Ice Measurements from Space (GLIMS) glacier database.