Moraine dams can be inherently unstable, but effective assessment strategies remain poorly identified. We integrate electrical resistivity tomography (ERT) with electrical self-potential (SP) and lake-level measurements to investigate the structure of, and hydrological processes within, a moraine-dam complex adjacent to Miage glacier, Italy. This complex separates two meltwater lakes characterized by an efficient subterranean hydraulic connection. Our ERT data reveal a continuous free surface within the complex, whose morphology reflects the topography of the moraine complex akin to unconfined groundwater aquifers. SP data were corrected for spatial changes in the thickness of the upper unsaturated layer using principles of electrography. The residual streaming-potential map is consistent with Darcian flow of lake waters through the moraine complex, characterized by a negative-to-positive potential change from <-30 mV to >+70 mV. These electrical signatures are consistent with those generated by water seepage through earth dams in various non-glacial settings. Integrated electrical geophysical methods thus provided an inexpensive and unobtrusive evaluation of the hydrological properties of and processes within the moraine complex. Since spatio-temporal patterns of subsurface water flow critically affect the strengths of moraine sediments, such methods promise to be powerful in assessing the long-term stability of moraine-dammed glacial lakes.

Due to remoteness and high altitude, only a few ground-based glacier change studies are available in high-mountain areas in the Himalaya. However, digital elevation models based on remotely sensed data (RS-DEMs) provide feasible opportunities to evaluate how fast Himalayan glaciers are changing. Here we compute elevation changes in glacier surface (total area 183.3 km2) in the Khumbu region, Nepal Himalaya, for the period 1992-2008 using multitemporal RS-DEMs and a map-derived DEM calibrated with differential GPS survey data in 2007. Elevation change is calculated by generating a weighted least-squares linear regression model. Our method enables us to provide the distribution of uncertainty of the elevation change. Debris-covered areas show large lowering rates. The spatial distribution of elevation change shows that the different wastage features of the debris-covered glaciers depend on their scale, slope and the existence of glacial lakes. The elevation changes of glaciers in the eastern Khumbu region are in line with previous studies. The regional average mass balance of -0.40 ± 0.25 m w.e.a-1 for the period 1992-2008 is consistent with a global value of about -0.55 m w.e. a-1 for the period 1996-2005.

A 6 year application of an enhanced temperature-index mass-balance model to Careser and La Mare glaciers, Eastern Italian Alps, is presented. The two glaciers exhibit very different characteristics, and a comprehensive dataset of distributed mass-balance measurements was used to test the model performance. The model was run using meteorological data acquired outside the glaciers. The work was focused on two main aspects: (1) the development of a morphological redistribution procedure for snow, and (2) the comparison of three different melt algorithms proposed in the literature. The results show that the simple method proposed for snow redistribution can greatly improve simulation of winter balance, and further improvements would be achievable by collecting data on inaccessible and high-altitude areas. All three melt formulations displayed a good skill level and very similar results in modelling the mass-balance distribution over glacier areas, with slightly better results from a multiplicative algorithm in capturing the vertical balance gradient. The simulation errors are related to aspect and elevation, and tend to be spatially aggregated. Some assumptions concerning the spatial and temporal distribution of air temperature and incoming solar radiation, although reasonable and widely used in the literature, may be responsible for this aggregation. Hence, there is a need to further investigate the processes that regulate the distribution of melt energy, and that appear to control the current deglaciation phase in this area.

Near ice shelves around Antarctica the ocean becomes supercooled and has been observed to carry small suspended ice crystals. Our measurements demonstrate that these small crystals are persistently present in the water column beneath the winter fast ice, and when incorporated in sea ice they reduce the mean grain size of the sea-ice cover. By midwinter, larger ice crystals below the ice/water interface are observed to form a porous sub-ice platelet layer with an ice volume fraction of 0.25 ± 0.06. The magnitude and direction of the oceanic heat flux varied between (5 ± 6) Wm-2 (upwards) and (-15 ± 10) Wm-2 (downwards) in May, but by September it settled between (-6 ± 2) and (-11 ± 2) W m-2. The negative values imply that the ocean acts as a heat sink which is responsible for the growth of 12% of the ice thickness between June and September. This oceanic contribution should not be ignored in models of Antarctic sea-ice thickness close to an ice shelf.

Speckle tracking of repeat RADARSAT-2 fine-beam imagery acquired over 24 day periods in March 2009 allowed the creation of updated surface motion maps for the entire Devon Ice Cap, Canada. Error analysis indicates that speckle tracking can determine ice motion to an accuracy of ~5 ma-1. Comparisons with earlier velocity maps from the mid-1990s and 2000 reveal velocity patterns that largely agree with flow regimes described previously. However, motion determined along East5 Glacier indicates an increase in surface velocities between the studies. Additionally, Southeast2 Glacier has significantly accelerated over the past decade, with velocities greater in 2009 than in the early 1990s along almost the entire length of the glacier. This is likely indicative of a surge. Present-day total mass loss from Devon Ice Cap due to iceberg calving is calculated as 0.40 ± 0.09 Gta-1, similar to that reported by Burgess and others (2005), with Belcher Glacier accounting for ~42% of the entire loss.

Unconformable firn stratigraphy exists throughout a 650 km long radar profile that we recorded down-flow of megadune fields in the Byrd Glacier (East Antarctica) catchment. Profile segments reveal cosets of prograding bedding sequences up to 90 m thick and with lateral, along-crest dimensions up to tens of kilometers. We profiled them in oblique section and nearly parallel to the prevailing wind. The prograding snow accumulates on broad, low windward slopes located above ice-bed depressions, which implies long-term slope stability. The apparent subglacial control implies that the accumulation progrades in balance with ice velocity, which we measured at ~30 ma”1. The sequences prograde over intensely modified and recrystallized wind-glaze firn, visible in the profiles as unstratified layers and zones up to several tens of meters thick. The intense recrystallization eliminates density stratification, and the altered layers appear to thicken into a connected network. Modeling of coset formation using wind and ice flow reproduces their dimensions and morphology. However, accumulation rates well above current regional estimates and existing data for megadunes are required because of the measured ice speed and required slope stability. The consistent unconformable strata along our traverse show that coset and recrystallized morphology extend far beyond the megadune fields.

The nature and behaviour of sediment beneath glaciers influences how they flow and respond to changing environmental conditions. The difficulty of accessing the bed of current glaciers is a key constraint to studying the processes involved. This paper explores an alternative approach by relating sediments under the beds of former mid-latitude ice sheets to changing ice behaviour during a glacial cycle. The paper focuses on the partly marine-based Pleistocene British-Irish ice sheet in the Clyde basin, Scotland. A three-dimensional computation of subsurface glacial sediment distribution is derived from 1260 borehole logs. Sediment distribution is linked to an empirically based reconstruction of ice-sheet evolution, permitting identification of distinctive phases of sedimentation. Maximum sediment mobilization and till deposition (~0.04ma_1) occurred during ice advance into the basin from adjacent uplands. Transport distances were generally short. Subglacial processes were influenced locally by the relative stiffness of pre-existing sediments, the permeability of the sub-till lithology, and topography; the resulting mean till thickness is 7.7 m with a high standard deviation of 7.0 m. In places, focused till deposition sealed pre-existing permeable substrates, promoting lower effective pressures. Sediment remobilization by meltwater was a key process as ice margins retreated through the, basin.

In order to apportion the dust sources of mountain glaciers in western China, the Sr-Nd isotopic compositions of insoluble particles were determined in snow samples collected from 13 sites. The combined plot of 87Sr/86Sr and εNd(0) demonstrates a distinctive geographic pattern over western China, which can be classified into three regions from north to south. Samples from the Altai mountains show the lowest 87Sr/86Sr ratio and the highest εNd(0) value, similar to the data of deserts in the north of China such as the Gurbantunggut desert. Samples from the southern Tibetan Plateau (TP) and Himalaya show the highest 87Sr/86Sr and lowest εNd(0) values, resembling the local and regional dust sources found in the southern TP and Himalaya-India region. Samples from the Tien Shan and northern Tibetan Plateau exhibit intermediate 87Sr/86Sr and εNd(0) values, similar to the data reported for the northern margin of the TP (NM_TP). However, three sampling sites, JMYZ (Jiemayangzong) located in the Himalaya and ZD (Zadang) and YL (Yulong) located in the southeast TP, presented distinctive Sr-Nd isotopic signatures typical of the NM_TP, suggesting potential long-range and high-altitude dust transport across the TP.

Ice velocities observed in 2005/06 at three GPS stations along the Sermeq Avannarleq flowline, West Greenland, are used to characterize an observed annual velocity cycle. We attempt to reproduce this annual ice velocity cycle using a 1-D ice-flow model with longitudinal stresses coupled to a 1-D hydrology model that governs an empirical basal sliding rule. Seasonal basal sliding velocity is parameterized as a perturbation of prescribed winter sliding velocity that is proportional to the rate of change of glacier water storage. The coupled model reproduces the broad features of the annual basal sliding cycle observed along this flowline, namely a summer speed-up event followed by a fall slowdown event. We also evaluate the hypothesis that the observed annual velocity cycle is due to the annual calving cycle at the terminus. We demonstrate that the ice acceleration due to a catastrophic calving event takes an order of magnitude longer to reach CU/ETH (‘Swiss’) Camp (46 km upstream of the terminus) than is observed. The seasonal acceleration observed at Swiss Camp is therefore unlikely to be the result of velocity perturbations propagated upstream via longitudinal coupling. Instead we interpret this velocity cycle to reflect the local history of glacier water balance.

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.

At the time of its calving in February 2010, Mertz Glacier, East Antarctica, was characterized by a 145 km long, 35 km wide floating tongue. In this paper, we use GPS data from the Collaborative Research into Antarctic Calving and Iceberg Evolution (CRAC-ICE) 2007/08 and 2009/10 field seasons to investigate the dynamics of Mertz Glacier. Two months of data were collected at the end of the 2007/08 field season from two kinematic GPS stations situated on each side of the main rift of the glacier tongue and from rock stations located around the ice tongue during 2008/09. Using Precise Point Positioning with integer ambiguity fixing, we observe that the two GPS stations recorded vibrations of the ice tongue with several dominant periods. We compare these results with a simple elastic model of the ice tongue and find that the natural vibration frequencies are similar to those observed. This information provides a better understanding of their possible effects on rift propagation and hence on the glacier calving processes.

A detailed history of volcanism covering the last 2840 years is reconstructed from the top 100.42 m of a 109.91 m ice core from Dome A (DA2005 ice core), East Antarctica. Using two known volcanic stratigraphic markers, the mean accumulation rate during the period AD 1260-1964 is found to be 23.2 mmw.e. a-1, consistent with the previously reported accumulation rate at Dome A. This mean accumulation rate is used to date the entire core. Volcanic eruptions in the period 840 BC-AD1998 are detected as outstanding sulphate events. Seventy-eight eruptions are identified, with a mean of 2.7 eruptions per century. Comparisons with previous Antarctic ice-core volcanic records are made to assess the quality of this new DA2005 record. In terms of dates for volcanic events, the DA2005 record is in good agreement with previous records in the second millennium ad (ad 1000-1998). A series of volcanic signatures found in both the DA2005 record and several other Antarctic ice-core records in the first millennium ad (ad 1-1000) appear to validate the DA2005 record during this time period. For the older periods, direct comparisons are difficult between the DA2005 record and other Antarctic ice-core records due to the lack of well-dated stratigraphic horizons.

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).

The processes that form intraglacial lakes are poorly understood because of the difficulty in detecting and assessing such hidden lakes. Extensive geophysical surveys were performed between 2007 and 2010 in order to reassess the risk of an outburst flood from Glacier de Tete Rousse, French Alps, where outburst flooding from an intraglacial lake caused 175 fatalities in 1892. Our geophysical survey combined ground-penetrating radar measurements and nuclear magnetic resonance imaging. We found a subglacial water-filled reservoir with a volume of 55 000 m3. Artesian outpourings occurred when the subglacial cavity was reached by two borehole drillings, indicating that the hydrostatic water pressure exceeded the ice pressure at the bottom of the cavity. On the basis of these geophysical and glaciological findings, we warned the public authorities in July 2010 of the risk facing the 3000 inhabitants downstream of the glacier. The subglacial reservoir was drained artificially. This example demonstrates how geophysical surveys can be used to detect this type of hazard when it is suspected, in particular when no hydraulic outlet from the snout exists. Numerous ice temperature measurements have shown that the tongue of this glacier is cold-based. This thermal regime could explain the accumulation of water in this glacier.

Complex unconformable englacial stratigraphy, including a segment of distinctive cosets of bed sequences, occurs throughout the thickness of a 3.2 MHz ice-sheet radar profile we acquired across the upper Byrd Glacier (East Antarctica) catchment. Some cosets span >10 km, are >100 m thick and are delineated by distinct horizons. At 40-90 m depth in firn, comparisons between 200 MHz and specially processed 3.2 MHz profiles reveal that the delineating horizons result from density-modified layers produced by decades to millennia of subaerial exposure, as detailed in our related paper (Part I). These comparisons, together with reflected waveforms at depth, also reveal that the modified layers retain their chemical stratification, and therefore the original unconformable surface. Two profile segments show high-amplitude transverse folds spanning much of the ice-sheet thickness. The parallel nature of most of them suggests basal sliding beneath long-term up-ice-flow accumulation zones, which we identify in satellite images as the likely sources for the cosets. The unconformable stratigraphy at depths greater than 2000 m shows that antidunal deposition and intense firn recrystallization zones have persisted for tens of thousands of years in this region of East Antarctica.

Global warming is causing widespread glacier retreat, with small glaciers disappearing. We investigate changes in glaciers over the western Lenglongling mountains, located in the northeastern margin of the Tibetan Plateau. Glacier extent over the western Lenglongling mountains is estimated by comparing digitized glacier outlines obtained from aerial photographs and satellite imagery. These results suggest that all 179 glaciers in the western Lenglongling mountains shrunk between 1972 and 2007. The total area loss was ~24.4 km2, accounting for ~28.3% (0.81% a-1) of the glacierized area in 1972. The average area retreat rates differ over different time intervals: they are approximately 0.68, 0.90, 0.77 and 0.56 km2 a-1 over the periods 1972-95,1995-99,1999-2002 and 2002-07, respectively. Based on analysis of meteorological data, glacier shrinkage in the study area can probably be attributed to the increase in air temperature. Furthermore, the smaller glaciers display a higher shrinkage rate than larger glaciers, and glaciers on southwest-facing slopes appear to retreat faster than those on northeast- facing slopes.

In order to account for the effects of debris cover in model scenarios of the response of glaciers to climate change and water resource planning, it is important to know the distribution and thickness of supraglacial debris and to monitor its change over time. Previous attempts to map surface debris thickness using thermal band remote sensing have relied upon time-specific empirical relationships between surface temperature and thickness, limiting their general applicability. In this paper, we develop a physically based model that utilizes Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER) thermal band remotely sensed imagery and is based on a solution of the energy balance at the debris surface. The model is used to estimate debris thickness on Miage glacier, Italy, and is validated using field debris-thickness measurements and a previously published debris-thickness map. The temporal transferability of the model is demonstrated through successful application to a separate ASTER image from a different year using reanalysis meteorological input data. This model has the potential to be used for regional-scale supraglacial debris-thickness mapping and monitoring for debris up to at least 0.50 m thickness, but improved understanding of the spatial patterns of air temperature, aerodynamic roughness length and thermal properties across debris-covered glaciers is needed.

The tongue of Triftgletscher, Switzerland, is particularly susceptible to major break-off events due to its steep slope. In order to detect precursors of such an event, we monitored the local seismic activity and detected 2426 icequakes with sources located in an area ranging between 2050 and 2350ma.s.l. Events triggered by cracks and icefalls were recorded, but no precise distinction between the two sources was possible using duration or frequency criteria. Clusters of seismic activity were located and confirmed by visual observations. We performed a surface motion analysis and found that surface motion was driven significantly by runoff changes at a timescale of 2-3 days. By means of a statistical analysis, a power-law behaviour of the released seismic energy distribution was discerned at certain times during the period investigated. Variations in power-law exponent values indicated that low- and high-energy events predominantly occurred during phases of enhanced and reduced surface motion, respectively. Substantial releases of seismic energy likely to signal the glacier recoupling were detected during phases of decreasing runoff. Clues to potential seismic precursors of break-off events are discussed.

Higher-order models represent a computationally less expensive alternative to the Stokes model for ice-sheet modeling. In this work, we develop linear and quadratic finite-element methods, implemented on parallel architectures, for the three-dimensional first-order model of Dukowicz and others (2010) that is based on the Blatter-Pattyn model, and for the depth-integrated model of Schoof and Hindmarsh (2010). We then apply our computational models to three of the ISMIP-HOM benchmark test cases (Pattyn and others, 2008). We compare results obtained from our models with those obtained using a reliable Stokes computational model, showing that our first-order model implementation produces reliable and accurate solutions for almost all characteristic length scales of the test geometries considered. Good agreement with the reference Stokes solution is also obtained by our depth-integrated model implementation in fast-sliding regimes and for medium to large length scales. We also provide a comprehensive comparison between results obtained from our first-order model implementation and implementations developed by ISMIP-HOM participants; this study shows that our implementation is at least as good as the previous ones. Finally, a comparison between linear and quadratic finite- element approximations is carried out, showing, as expected, the better accuracy of the quadratic finite-element method.