Because the glacier snowline is easy to identify on optical satellite images and because in certain conditions it can be used as an indicator of the equilibrium line, it may be a relevant parameter for the study of the relationships between climate and glaciers. Although several studies have shown that the snowline altitude (SLA) at the end of the hydrological year is a good indicator of the equilibrium-line altitude (ELA) for mid-latitude glaciers, such a relationship remains conjectural for tropical glaciers. Indeed, unlike in mid-latitudes, tropical climate conditions result in a distinct seasonality of accumulation/ablation processes. We examine this relationship using direct field ELA and mass-balance measurements made on Glaciar Zongo, Bolivia (~16° S), vand Glaciar Artesonraju, Peru (~9° S), and the SLA retrieved from satellite images acquired in the past two decades. We show that on glaciers in the outer tropics: (1) ablation is reduced during the dry season in austral winter (May-August), the SLA does not change much, and satellite images acquired between May and August could be used to compute the SLA; and (2) the highest SLA detected on a number of satellite images acquired during the dry season provides a good estimate of the annual ELA. However, as snowfall events can occur during the dry season, the SLA detected on satellite images tends to underestimate the ELA. Thus, we recommend validating the SLA computed from satellite images with field data collected on a benchmark glacier before measuring the SLA on other glaciers in the same mountain range for which no field data are available. This study is a major step towards extending the measurement of glacier parameters (ELA and mass balance) at the scale of a whole mountain range in the outer tropics to better document the relationships between climate and glaciers.

The subaqueous margins of calving glaciers have the potential to make significant contributions to glacier mass loss. However, to date, very little is known about the morphology and development of subaqueous margins. A unique combination of sub-bottom profile and bathymetric data collected between 2008 and 2010 in proglacial lakes at Mueller, Hooker and Tasman glaciers in New Zealand’s Southern Alps reveal subaqueous ice ramps extending up to 510 m from the terminus of each glacier. Ice ramp surfaces are undulating and covered with a thick layer (up to 10 m) of unsorted sediment derived from supraglacial and englacial debris, lateral moraines and deltaic deposits. A cyclic calving pattern, relatively stable lake level and the debris cover appear to control the development and maintenance of these ice ramps. High subaerial retreat rates generally correspond to high subaqueous calving rates, although the highest subaerial retreat rates are not associated with the largest ice ramp. Debris mantling the subaqueous ice ramp surfaces insulates the ice from melting and also reduces buoyant forces acting on the terminus. Comparisons with previous studies show that the ice ramps evolve over time with changes in glacier dynamics and water-body properties.

Modelling the hydrology of the Greenland ice sheet, including the filling and drainage of supraglacial lakes, requires melt inputs generated at high spatial and temporal resolution. Here we apply a high spatial (100 m) and temporal (1 hour) mass-balance model to a 450 km2 subset of the Paakitsoq region, West Greenland. The model is calibrated by adjusting the values for parameters of fresh snow density, threshold temperature for solid/liquid precipitation and elevation-dependent precipitation gradient to minimize the error between modelled output and surface height and albedo measurements from three Greenland Climate Network stations for the mass-balance years 2000/01 and 2004/05. Bestfit parameter values are consistent between the two years at 400 kg m-3, 2°C and +14% (100 m)-1, respectively. Model performance is evaluated, first, by comparing modelled snow and ice distribution with that derived from Landsat-7 ETM+ satellite imagery using normalized-difference snow index classification and supervised image thresholding; and second, by comparing modelled albedo with that retrieved from the MODIS sensor M0D10A1 product. Calculation of mass-balance components indicates that 6% of surface meltwater and rainwater refreezes in the snowpack and does not become runoff, such that refreezing accounts for 31% of the net accumulation.

We used Little Ice Age (LIA) trimlines and moraines to assess changes in South American glaciers over the last ~140 years. We determined the extent and length of 640 glaciers during the LIA (~AD1870) and 626 glaciers (the remainder having entirely disappeared) in 1986, 2001 and 2011. The calculated reduction in glacierized area between the LIA and 2011 is 4131 km2 (15.4%), with 660 km2 (14.2%) being lost from the Northern Patagonia Icefield (NPI), 1643km2 (11.4%) from the Southern Patagonia Icefield (SPI) and 306 km2 (14.4%) from Cordillera Darwin. Latitude, size and terminal environment (calving or land-terminating) exert the greatest control on rates of shrinkage. Small, northerly, land-terminating glaciers shrank fastest. Annual rates of area loss increased dramatically after 2001 for mountain glaciers north of 52° S and the large icefields, with the NPI and SPI now shrinking at 9.4km2a-1 (0.23% a-1) and 20.5 km2a-1 (0.15% a-1) respectively. The shrinkage of glaciers between 52° S and 54° S accelerated after 1986, and rates of shrinkage from 1986 to 2011 remained steady. Icefield outlet glaciers, isolated glaciers and ice caps south of 54° S shrank faster from 1986 to 2001 than they did from 2001 to 2011.

Since their discovery, Antarctic subglacial lakes have become of great interest to the science community. It is hypothesized that they may hold unique forms of biological life and that they hold detailed sedimentary records of past climate change. According to the latest inventory, a total of 387 subglacial lakes have been identified in Antarctica (Wright and Siegert, 2011). However, exploration using scientific probes has yet to be performed. We propose a generic, formal approach to manage the operational risk of deploying probes during clean access to subglacial lake exploration. A representation of the entire probe deployment process is captured in a Markov chain. The transition from one state to the next depends on several factors, including reliability of components and processes. We use fault trees to quantify the probability of failure of the complex processes that must take place to facilitate the transition from one state to another. Therefore, the formal framework consists of integrating a Markov chain, fault trees, component and subsystem reliability data and expert judgment. To illustrate its application we describe how the approach can be used to address a series of what-if scenarios, using the intended Ellsworth Subglacial Lake probe deployment as a case study.

Microbially mediated carbon fluxes on the surface of the Greenland ice sheet (GrIS) were recently quantified by Hodson and others (2010) using measurements of the surface coverage of debris (cryoconite) and rates of biological production associated with debris near the ice-sheet margin. We present updated models that do not assume the same spatial uniformity in key parameters employed by Hodson and others (2010) because they make use of biomass distribution and biological production data from a 79 km transect of the GrIS. Further, the models presented here also include for the first time biomass associated with both cryoconite holes and surficial algae. The predicted annual carbon flux for a small (1600 km2) section of ice surrounding the field transect is about four times that estimated using spatially uniform biomass and production in this area. When surficial algae are included, the model predicts about 11 times more carbon fixation via photosynthesis per year than the cryoconite-only models. We therefore suggest that supraglacial carbon fluxes from the GrIS have previously been underestimated by more than an order of magnitude and that the hitherto overlooked surficial algal ecosystem can be the most crucial contributor. The GrIS is shown to be in a relatively stable state of net autotrophy according to our model and so a strong link between bare-ice area and total carbon fluxes is evident. The implication is a biomass feedback to surface albedo and enhanced ablation as a result. Climate predictions for the year 2100 show that greater carbon fixation could also result from climate warming.

Chemical constituents trapped within glacial ice provide a unique record of climate, as well as repositories for biological material such as pollen grains, fungal spores, viruses, bacteria and dissolved organic carbon. Past research suggests that the veins of polycrystalline ice may provide a liquid microenvironment for active microbial metabolism fueled by concentrated impurities in the veins. Despite these claims, no direct measurements of impurity concentration in ice veins have been made. Using micro-Raman spectroscopy, we show that sulfate and nitrate concentrations in the veins of glacial ice from Greenland (Greenland Ice Sheet Project 2) and Antarctic (Newall Glacier and a Dominion Range glacier) core samples were 104 and 105 times greater than the concentrations measured in melted (bulk) core water. Methanesulfonate was not found in the veins, consistent with its presence as particulate matter within the ice. The measured vein concentration of molecular anions implies a highly acidic (pH < 3) vein environment with high ionic strength (mM-M). We estimate that the vein volume provides 16.7 and 576 km3 of habitable space within the Greenland and Antarctic ice sheets, respectively, which could support the metabolism of organisms that are capable of growing in cold, high ionic strength solutions with low pH.

Given the cross-sectional geometry of a valley glacier, effects of lateral drag can be parameterized in flowline models through the introduction of Nye shape factors. Lateral drag also arises due to lateral variability in bed topography and basal flow, which induce horizontal shear stress and differential ice motion. For glaciers with various geometric and basal conditions, we compare three-dimensional Stokes solutions to flowline model solutions to examine both sources of lateral drag. We calculate associated correction factors that help flowline models to capture the effects of lateraldrag. Such parameterizations provide improved simulations of the dynamics of narrow, channelized, fast-flowing glacial systems. We present an example application for Athabasca Glacier in the Canadian Rocky Mountains.

Ice-thickness measurements critical for flood and mudflow hazard studies are very sparse on Cascade Range (North America) volcanoes. Helicopter electromagnetic (HEM) data collected to detect hydrothermal alteration are used to determine ice thickness over portions of Mount Baker and Mount Adams volcanoes. A laterally continuous inversion method provides good estimates of ice <100 m thick over water-saturated and altered regions where the resistivity of the basement is <200 Ωm. For areas with ice overlying fresh, resistive rocks with small resistivity contrasts between ice and rock, ice thickness is not well resolved. The ice thicknesses derived from HEM data are consistent with the previous drillhole data from Mount Adams and radar data from both volcanoes, with mean thicknesses of 57 m for Mount Adams and 68 m for Mount Baker. The thickest ice on Mount Baker rests on the gentle lower slopes whereas the thickest ice at Mount Adams lies on the flat summit. Ice volume calculations suggest that Mount Baker contains ∽710 × 106 m3 of ice in the HEM survey area, with a crude estimate of ∽1800 × 106m3 for the entire volcano. Ice volume on Mount Adams is 65 × 106m3 in parts of the HEM survey area and ∽200 × 106m3 overall.

We present a method of measuring ice ablation using an absolute pressure transducer as part of an automatic weather station (AWS) system, which we have installed in 17 locations on the Greenland ice sheet. The pressure transducer assembly is drilled into the ice, enclosed in a hose filled with antifreeze liquid. The pressure signal registered by the transducer is that of the vertical column of liquid over the sensor, which can be translated to depth, and ice ablation rate, knowing the density of the liquid. Measuring at sub-daily timescales, this assembly is well suited to monitoring ice ablation in remote regions, with clear advantages over other, well-established methods. The pressure transducer system has the potential to monitor ice ablation for several years without re-drilling, and the system is suitable for high-ablation areas (>5ma-1). A routine to transform raw measurements into ablation values is presented, including a physically based method to remove air-pressure variability from the signal. The pressure transducer time series is compared to that recorded by a sonic ranger for the climatically hostile setting on the Greenland ice sheet.

Information about glacier volume and ice thickness distribution is essential for many glaciological applications, but direct measurements of ice thickness can be difficult and costly. We present a new method that calculates ice thickness via an estimate of ice flux. We solve the familiar continuity equation between adjacent flowlines, which decreases the computational time required compared to a solution on the whole grid. We test the method on Columbia Glacier, a large tidewater glacier in Alaska, USA, and compare calculated and measured ice thicknesses, with favorable results. This shows the potential of this method for estimating ice thickness distribution of glaciers for which only surface data are available. We find that both the mean thickness and volume of Columbia Glacier were approximately halved over the period 1957–2007, from 281 m to 143 m, and from 294 km3 to 134 km3, respectively. Using bedrock slope and considering how waves of thickness change propagate through the glacier, we conduct a brief analysis of the instability of Columbia Glacier, which leads us to conclude that the rapid portion of the retreat may be nearing an end.

A 141m ice core was recovered from Combatant Col (51.385° N, 125.258° W; 3000ma.s.l.), Mount Waddington, Coast Mountains, British Columbia, Canada. Records of black carbon, dust, lead and water stable isotopes demonstrate that unambiguous seasonality is preserved throughout the core, despite summer surface snowmelt and temperate ice. High accumulation rates at the site (>4 m ice eq. a-1) limit modification of annual stratigraphy by percolation of surface meltwater. The ice-core record spans the period 1973–2010. An annually averaged time series of lead concentrations from the core correlates well with historical records of lead emission from North America, and with ice-core records of lead from the Greenland ice sheet. The depth-age scale for the ice core provides sufficient constraint on the vertical strain to allow estimation of the age of the ice at bedrock. Total ice thickness at Combatant Col is ~250 m; an ice core to bedrock would likely contain ice in excess of 200 years in age. Accumulation at Combatant Col is significantly correlated with both regional precipitation and large-scale geopotential height anomalies.

Virtually no information is available on the response of land-terminating Antarctic Peninsula glaciers to climate change on a centennial timescale. This paper analyses the topography, geomorphology and sedimentology of prominent moraines on James Ross Island, Antarctica, to determine geometric changes and to interpret glacier behaviour. The moraines are very likely due to a late-Holocene phase of advance and featured (1) shearing and thrusting within the snout, (2) shearing and deformation of basal sediment, (3) more supraglacial debris than at present and (4) short distances of sediment transport. Retreat of ~100 m and thinning of 15–20 m has produced a loss of 0.1 km3 of ice. The pattern of surface lowering is asymmetric. These geometrical changes are suggested most simply to be due to a net negative mass balance caused by a drier climate. Comparisons of the moraines with the current glaciological surface structure of the glaciers permits speculation of a transition from a polythermal to a cold-based thermal regime. Small land-terminating glaciers in the northern Antarctic Peninsula region could be cooling despite a warming climate.

Southern Ocean hydrography has undergone substantial changes in recent decades, concurrent with an increase in the rate of Antarctic ice-shelf melting (AISM). We investigate the impact of increasing AISM on hydrography through a twin numerical experiment, with and without AISM, using a global coupled sea-ice/ocean climate model. The difference between these simulations gives a qualitative understanding of the impact of increasing AISM on hydrography. It is found that increasing AISM tends to freshen the surface water, warm the intermediate and deep waters, and freshen and warm the bottom water in the Southern Ocean. Such effects are consistent with the recent observed trends, suggesting that increasing AISM is likely a significant contributor to the changes in the Southern Ocean. Our analyses indicate potential positive feedback between hydrography and AISM that would amplify the effect on both Southern Ocean hydrography and Antarctic ice-shelf loss caused by external factors such as changing Southern Hemisphere winds.

We extend through 2011 an ice-sheet elevation and surface velocity record across three measurement networks established in south-central Greenland by The Ohio State University in 1980/81. Surface parameters are derived from repeat GPS in situ observations, elevations measured by airborne laser altimetry and by the Ice, Cloud and land Elevation Satellite (ICESat). Elevations at the western network steadily rose early in the record by 0.10 ± 0.02 m a-1, but an eastward-progressing thinning trend began in the mid-1990s followed by a ~1m elevation drop at all stations from 2005 to 2011. Measurements weakly suggest a surface velocity increase at the western cluster from 1980 to 2005. At the central network, elevations rose by 0.08 ± 0.02 m a-1 through 2005 and surface speed increased by 0.5–0.7 m a-1. Surface elevations at the central network remained nearly constant thereafter. Thickening occurred at the southern ice divide by 0.05 ± 0.02 m a-1, while east of the divide the ice sheet thinned with increasing rate from the divide, likely because of decreasing accumulation rate trends and drawdown into rapidly retreating coastal glaciers. Our most recent data show that thinning rates are slowing at several sites just east of the divide and that the elevation at the divide continues to increase.

The loss of the floating ice tongue on Jakobshavn Isbrӕ, Greenland, in the early 2000s has been concurrent with a pattern of thinning, retreat and acceleration leading to enhanced contribution to global sea level. These changes on decadal timescales have been well documented. Here we identify how the glacier responds to forcings on shorter timescales, such as from variations in surface melt, the drainage of supraglacial lakes and seasonal fluctuations in terminus position. Ice motion and surface melt were monitored intermittently from 2006 to 2008. Dual-frequency GPS were deployed 20–50 km upstream of the terminus along the glacier center line. Gaps in surface melt measurements were filled using a temperature-index model of ablation driven by surface air temperatures recorded during the same time period. Our results corroborate the premise that the primary factors controlling speeds on Jakobshavn Isbrӕ are terminus position and geometry. We also observe that surface speeds demonstrate a complex relationship with meltwater input: on diurnal timescales, velocities closely match changes in water input; however, on seasonal timescales a longer, more intense melt season was observed to effectively reduce the overall ice flow of the glacier for the whole year

A numerical model for an interacting ice shelf and ocean is presented in which the ice- shelf base exhibits a channelized morphology similar to that observed beneath Petermann Gletscher’s (Greenland) floating ice shelf. Channels are initiated by irregularities in the ice along the grounding line and then enlarged by ocean melting. To a first approximation, spatially variable basal melting seaward of the grounding line acts as a steel-rule die or a stencil, imparting a channelized form to the ice base as it passes by. Ocean circulation in the region of high melt is inertial in the along-channel direction and geostrophically balanced in the transverse direction. Melt rates depend on the wavelength of imposed variations in ice thickness where it enters the shelf, with shorter wavelengths reducing overall melting. Petermann Gletscher’s narrow basal channels may therefore act to preserve the ice shelf against excessive melting. Overall melting in the model increases for a warming of the subsurface water. The same sensitivity holds for very slight cooling, but for cooling of a few tenths of a degree a reorganization of the spatial pattern of melting leads, surprisingly, to catastrophic thinning of the ice shelf 12 km from the grounding line. Subglacial discharge of fresh water along the grounding line increases overall melting. The eventual steady state depends on when discharge is initiated in the transient history of the ice, showing that multiple steady states of the coupled system exist in general.

Firn microstructure properties from six different sites in Greenland and Antarctica are investigated by means of X-ray microcomputer tomography. The optical effective radius is calculated from the specific surface area (SSA) and used as a measure of grain size. It is shown that the recently introduced spherical approximation of firn grains using the effective radius Reff is representative of grain size in the microwave frequency region. The measured profiles show the well-known increase of grain size with depth at all sites, where the increase is largest at near-surface depths. A large variability in grain size on the decimeter-to-centimeter scale as a result of different grain properties of single layers is superimposed on the overall trend at each site. A simple empirical parameterization of grain-size evolution is developed which allows the rapid grain growth in the uppermost layers of the firn to be predicted. The growth is driven by strong seasonal and diurnal temperature gradients. The model can be used to simulate grain-size profiles required by models of firn/microwave interaction (e.g. for retrieval of accumulation rates from satellite microwave sensors) in a more realistic fashion.