The ability to perform microbial detection and characterization in-field at extreme environments, rather than on returned samples, has the potential to improve the efficiency, relevance and quantity of data from field campaigns. To date, few examples of this approach have been reported. Therefore, we demonstrate that the approach is feasible in subglacial environments by deploying four techniques for microbial detection: real-time polymerase chain reaction; microscopic fluorescence cell counts, adenosine triphosphate bioluminescence assay and recombinant Factor C assay (to detect lipopolysaccharide). Each technique was applied to 12 subglacial ice samples, 12 meltwater samples and two snow samples from Engabreen, Northern Norway. Using this multi-technique approach, the detected biomarker levels were as expected, being highest in debris-rich subglacial ice, moderate in glacial meltwater and low in clean ice (debris-poor) and snow. Principal component analysis was applied to the resulting dataset and could be performed in-field to rapidly aid the allocation of resources for further sample analysis. We anticipate that in-field data collection will allow for multiple rounds of sampling, analysis, interpretation and refinement within a single field campaign, resulting in the collection of larger and more appropriate datasets, ultimately with more efficient science return.

We use satellite radar interferometry to investigate changes in the location of the Petermann Glacier grounding line between 1992 and 2011. The grounding line location was identified in 17 quadruple-difference interferograms produced from European Remote Sensing (ERS)-1/2 data – the most extensive time series assembled at any ice stream to date. There is close agreement (20.6 cm) between vertical displacement of the floating ice shelf and relative tide amplitudes simulated by the Arctic Ocean Dynamics-based Tide Model 5 (AODTM-5) Arctic tide model. Over the 19 a period, the groundling line position varied by 470 m, on average, with a maximum range of 7.0 km observed on the north-east margin of the ice stream. Although the mean range (2.8 km) and variability (320 m) of the grounding line position is considerably lower if the unusually variable north-east sector is not considered, our observations demonstrate that large, isolated movements cannot be precluded, thus sparse temporal records should be analysed with care. The grounding line migration observed on Petermann Glacier is not significantly correlated with time (R2 = 0.22) despite reported ice shelf thinning and episodes of large iceberg calving, which suggests that unlike other ice streams, on the south-west margin of the Greenland ice sheet, Petermann Glacier is dynamically stable.

Supraglacial deposits are known for their influence on glacier ablation. The magnitude of this influence depends on the thickness and the type of the deposited material. The effects of thin layers of atmospheric black carbon and of thick moraine debris have been intensively studied. Studies related to regional-scale deposits of volcanic tephra with thicknesses varying between millimetres and metres and thus over several orders of magnitude are scarce. We present results of a field experiment in which we investigated the influence of supraglacial deposits of tephra from Grímsvötn volcano on bare-ice ablation at Svínafelsjökull, Iceland. We observed that the effective thickness at which ablation is maximized ranges from 1.0 to 2.0 mm. At ~10 mm a critical thickness is reached where sub-tephra ablation equals bare-ice ablation. We calibrated two empirical ablation models and a semi-physics-based ablation model that all account for varying tephra-layer thicknesses. A comparison of the three models indicates that for tephra deposits in the lower-millimetre scale the temperature/radiation-index model performs best, but that a semi-physics-based approach could be expected to yield superior results for tephra deposits of the order of decimetres.

In order to better understand the densification of polar firn, firn cores from the three sites within ~10 km of Dome Fuji, Antarctica, were investigated using surrogates of density: dielectric permittivities εv and εh at microwave frequencies with electrical fields in the vertical and horizontal planes respectively. Dielectric anisotropy Δε (=εv − εh) was then examined as a surrogate of the anisotropic geometry of firn. We find that layered densification is explained as a result of complex effects of two phenomena that commonly occur at the three sites. Basically, layers with initially smaller density and smaller geometrical anisotropy deform preferentially throughout the densification process due to textural effects. Second, layers having a higher concentration of Cl− ions deform preferentially during a limited period from the near surface depths until smoothing out of layered Cl− ions by diffusion. We hypothesize that Cl− ions dissociated from sea salts soften firn due to modulation of dislocation movement. Moreover, firn differs markedly across the three sites in terms of strength of geometrical anisotropy, mean rate of densification and density fluctuation. We hypothesize that these differences are caused by textural effects resulting from differences in depositional conditions within various spatial scales.

Glacier basal motion generates diurnal to multi-annual fluctuations in glacier velocity and mass flux. Understanding these fluctuations is important for prediction of future sea-level rise and for gaining insight into glacier physics and erosion. Here, we derive glacier velocity through cross-correlation of WorldView satellite imagery to document the evolution of ice surface velocity on Kennicott Glacier, Alaska, over the 2013 melt season. The summer speedup is spatially uniform over a ~12 km2 area, over which the spring velocity varies significantly. Velocity increases by 1.4-fold to tenfold across the study domain, with larger values where spring velocities are low. To investigate the cross-glacier distribution of basal motion required to explain the observed surface speedup, we employ a two-dimensional cross-sectional glacier flow model. We find the model is insensitive to the spatial distribution of basal slip because stress gradient ice coupling diffuses the surface expression of the basal velocity field. While the temporal evolution of the subglacial hydrologic system is critical for predicting a glacier's response to meltwater inputs, our work suggests that glacier and ice-sheet models do not require a detailed representation of subglacial hydrology to accurately capture the spatial pattern of glacier speedup.

Lake-ice properties at Lake Abashiri, Hokkaido, Japan, were examined using field observations and a 1-D thermodynamic model to clarify formation processes at mid-latitudes subject to significant snowfall as well as moderate air temperature. At all lake sites examined, the ice comprised two distinct layers: a snow ice (SI) layer on top and a congelation ice (CI) layer below. The SI layer occupied as much as 29–73% of the total ice thickness, a much greater fraction than that reported for lakes at Arctic high latitudes. In the model, the CI growth rate was estimated using the traditional heat budget method, while the SI growth rate was calculated assuming the excessive snowfall from the isostatic balance is converted to SI by a snow compression rate (β) with the surface melting rate added when the surface heat budget becomes positive. By tuning the value of β to the observational results of SI thickness, the model outcome successfully reproduced the observational thicknesses of CI and SI, and the break-up date of the lake. Essentially, the model findings show how snow and its formation into SI reduce, by about half, the seasonal variability of total ice thickness.

In this paper, we reanalyze the glacier mass balance on Glaciar Antisana 15α over the 1995–2012 period. Annual glacier mass balances were quantified on the basis of monthly glaciological measurements using an adaptation of Lliboutry's statistical approach. The geodetic mass balance was computed between 1997 and 2009 giving a cumulative balance of −1.39 ± 1.97 m w.e. and a slightly negative adjusted annual glaciological mass balance (−0.12 ± 0.16 m w.e. a−1). Despite a careful analysis of uncertainties, we found a large discrepancy between the cumulative glaciological and the geodetic mass balances over the common period, of 4.66 m w.e. This discrepancy can mainly be explained by underestimated net accumulation in the glacier upper reaches, which could be due to the peculiar climate conditions of the equatorial zone with year round accumulation, thereby preventing clear identification of annual layers. An increase of ~70% in measured rates of net accumulation would be needed to balance the glaciological and geodetic mass balances; a hypothesis confirmed by estimated ice flux in the vicinity of the ELA. Consequently, the vertical gradient of precipitation may be higher than previously estimated and the accumulation processes (including the role of frost deposition) need to be carefully analyzed.

Results of numerical simulations of co-axial deformation of pure ice up to high-strain, combining full-field modelling with recrystallisation are presented. Grain size and lattice preferred orientation analysis and comparisons between simulations at different strain-rates show how recrystallisation has a major effect on the microstructure, developing larger and equi-dimensional grains, but a relatively minor effect on the development of a preferred orientation of c-axes. Although c-axis distributions do not vary much, recrystallisation appears to have a distinct effect on the relative activities of slip systems, activating the pyramidal slip system and affecting the distribution of a-axes. The simulations reveal that the survival probability of individual grains is strongly related to the initial grain size, but only weakly dependent on hard or soft orientations with respect to the flow field. Dynamic recrystallisation reduces initial hardening, which is followed by a steady state characteristic of pure-shear deformation.

Glacier change has been recognized as an important climate variable due to its sensitive response to climate change. Although there are a large number of glaciers distributed over the southeastern Qinghai–Tibetan Plateau, the region is poorly represented in glacier databases due to seasonal snow cover and frequent cloud cover. Here, we present an improved glacier inventory for this region by combining Landsat observations acquired over 2011–13 (Landsat 8/OLI and Landsat TM/ETM+), coherence images from Advanced Land Observing Satellite Phased Array type L-band Synthetic Aperture Radar images and the Shuttle Radar Topography Mission (SRTM) DEM. We present a semi-automated scheme for integrating observations from multi-temporal Landsat scenes to mitigate cloud obscuration. Further, the clean-ice observations, together with coherence information, slope constraints, vegetation cover and water classification information extracted from the Landsat scenes, are integrated to determine the debris-covered glacier area. After manual editing, we derive a new glacier inventory containing 6892 glaciers >0.02 km2, covering a total area of 6566 ± 197 km2. This new glacier inventory indicates gross overestimation in glacier area (over 30%) in previously published glacier inventories, and reveals various spatial characteristics of glaciers in the region. Our inventory can be used as a baseline dataset for future studies including glacier change assessment.

We investigate an internal surge of Karayaylak Glacier, which was reported by the media in May 2015. To differentiate the May 2015 glacier surge from other glacier advances, we surveyed changes in velocity, crevasses and glacier area using Landsat 8 OLI L1T, ZY-1-02C and Gaofen-1 images from October 2014 to July 2015. The velocity, measured by automatic feature extraction and tracking during the active phase, was 10–100 times the velocity during the quiescent phase, with a maximum of (20.2 ± 0.9) m d−1 (mean ± standard error) from 8 to 15 May 2015 in the west branch of the glacier. The surge initiation and termination took place from 13 April to 16 June 2015. Ice in the west branch (length, 7 km; area, 6.8 km2) of Karayaylak Glacier accelerated down to the east branch, leading to the development of crevasses and ice covering an additional 0.1 km2 of summer pasture on the northwestern side. However, we detected no advance of the glacier's terminus.

To explore links between glacier dynamics, sediment yields and the accumulation of glacial sediments in a temperate setting, we use extensive glaciological observations for Columbia Glacier, Alaska, and new oceanographic data from the fjord exposed during its retreat. High-resolution seismic data indicate that 3.2 × 108 m3 of sediment has accumulated in Columbia Fjord over the past three decades, which corresponds to ~5 mm a−1 of erosion averaged over the glaciated area. We develop a general model to infer the sediment-flux history from the glacier that is compatible with the observed retreat history, and the thickness and architecture of the fjord sediment deposits. Results reveal a fivefold increase in sediment flux from 1997 to 2000, which is not correlated with concurrent changes in ice flux or retreat rate. We suggest the flux increase resulted from an increase in the sediment transport capacity of the subglacial hydraulic system due to the retreat-related steepening of the glacier surface over a known subglacial deep basin. Because variations in subglacial sediment storage can impact glacial sediment flux, in addition to changes in climate, erosion rate and glacier dynamics, the interpretation of climatic changes based on the sediment record is more complex than generally assumed.

Thinning rates for the debris-covered Gangotri Glacier and its tributary glaciers during the period 1968–2014, length variation and area vacated at the snout from 1965 to 2015, and seasonal variation of ice-surface velocity for the last two decades have been investigated in this study. It was found that the mass loss of Gangotri and its tributary glaciers was slightly less than those reported for other debris-covered glaciers in the Himalayan regions. The average velocity during 2006–14 decreased by ~6.7% as compared with that during 1993–2006. The debris-covered area of the main trunk of Gangotri Glacier increased significantly from 1965 until 2015 with the maximum rate of increase (0.8 ± 0.2 km2 a−1) during 2006–15. The retreat (~9.0 ± 3.5 m a−1) was less in recent years (2006–2015) but the down-wasting (0.34 ± 0.2 m a−1) in the same period (2006–2014) was higher than that (0.20 ± 0.1 m a−1) during 1968–2006. The study reinforced the established fact that the glacier length change is a delayed response to climate change and, in addition, is affected by debris cover, whereas glacier mass balance is a more direct and immediate response. Therefore, it is recommended to study the glacier mass balance and not only the glacier extent, to conclude about a glacier's response to climate change.

Balls of frost (yukimarimo) up to 12 cm in diameter and weighing up to 14 g were observed forming in light winds in the temperature range, −70 to −60 °C at Dome Charlie, Antarctica, during the 2014 winter. Their density ranged from 15 to 60 kg m−3. Meteorological conditions during four formation periods are tabulated. Their formation seems to be due to a clumping of crystals caused by electrostatic attraction, followed by the collection of hoarfrost crystals while rolling along the snow surface.

Rapid changes in Helheim Glacier and other Greenland outlet glaciers since 2000 are well-known, but knowledge on earlier decades is fragmentary. Here we exploit the satellite image archives to produce and analyze a monthly-to-seasonal record of Helheim Glacier front position, 1980–2011. Statistical analysis identifies decadal periods with abrupt changes in variability and mean. The record also reveals evidence of volatile advance/retreat behavior in the 1980s. In one of several cases of large-amplitude subannual changes, the glacier front ‘surged’ forward in 1984/85, advancing ~6 km within a few months – surpassing its Little Ice Age maximum position – and afterward retreated ~5 km within a few weeks. These findings challenge the prevailing view of front position stability in the decades before the multi-year retreat in the early 2000s. Cold conditions including rigid ice mélange appear to be a factor in the high-amplitude seasonal advances in the 1980s. However the magnitude and abruptness of the changes in the record cannot be explained solely as a climatic response, such that glacio-dynamics must be invoked. Further, the volatile advance/retreat behavior in the cold 1980s resulted in increased dynamic ice loss, complicating the interpretation of increased calving activity as a response to warming.

The presence of water-filled crevasses is known to increase the penetration depth of crevasses and this has been hypothesized to play an important role controlling iceberg calving rate. Here, we develop a continuum-damage-based poro-mechanics formulation that enables the simulation of water-filled basal and surface crevasse propagation. The formulation incorporates a scalar isotropic damage variable into a Maxwell-type viscoelastic constitutive model for glacial ice, and the effect of the water pressure on fracture propagation using the concept of effective solid stress. We illustrate the model by simulating quasi-static hydrofracture in idealized rectangular slabs of ice in contact with the ocean. Our results indicate that water-filled basal crevasses only propagate when the water pressure is sufficiently large, and that the interaction between simultaneously propagating water-filled surface and basal crevasses can have a mutually positive influence leading to deeper crevasse propagation, which can critically affect glacial stability. Therefore, this study supports the hypothesis that hydraulic fracture is a plausible mechanism for the accelerated breakdown of glaciers.

The structural evolution of snow under metamorphism is one of the key challenges in snow modeling. The main driving forces for metamorphism are curvature differences and temperature gradients, inducing water vapor transport and corresponding crystal growth, which is detectable by the motion of the ice/air interface. To provide quantitative means for a microscopic validation of metamorphism models, a VTK-based image analysis method is developed to track the ice/air interface in time-lapse μCT experiments to measure local interface velocities under both, isothermal and temperature gradient conditions. Using estimates of local temperatures from microstructure-based finite element simulations, a quantitative comparison of measured interface velocities with theoretical expressions is facilitated. For isothermal metamorphism, the data are compared with a kinetics and a diffusion limited growth law. In both cases the data are largely scattered but consistently show a mean curvature dependency of the interface velocity. For temperature gradient metamorphism, we confirm that the main contribution stems from the temperature gradient induced vapor flux, accompanied by effects of mean curvature as a secondary process. The scatter and uncertainties are discussed in view of the present theoretical understanding, the experimental setup and complications such as mechanical deformations.

Time-varying elevations near the calving front of Jakobshavn Isbræ, Greenland were observed with a terrestrial radar interferometer (TRI) in June 2015. An ice block with surface dimensions of 1370 m × 290 m calved on 10 June. TRI-generated time series show that ice elevation near the calving front began to increase 65 h prior to the event, and can be fit with a simple block rotation model. We hypothesize that subsurface melting at the base of the floating terminus breaks the gravity-buoyancy equilibrium, leading to slow subsidence and rotation of the block, and its eventual failure.

Recent studies have suggested a physical link between Ca++ ions and an increase in the ductility or ‘softening’ of polycrystalline ice. In order to investigate the potential effects of Ca++ on deformation, we created sets of both undoped and CaSO4-doped specimens of polycrystalline ice for testing in uniaxial tension or compression. Deformation tests in tension were carried out under a constant load at an initial stress of 0.75 MPa and a temperature of −6°C. Compression tests were carried out at −10 and −20°C at constant strain rates of 1×10−4 s−1, 1 × 10−5 s−1 and 1 × 10−6 s−1 and taken to 5% strain. Our results show that CaSO4 increases the strength of polycrystalline ice at higher strain rates and lower temperatures, an effect that decreases with decreasing strain rate and higher temperatures. A microstructural analysis of the post-test compression specimens reveals mean grain diameters much larger in the CaSO4-doped specimens tested at the lowest applied strain rate of 1 × 10−6 s−1. Precipitates were found to have formed along grain boundaries in some doped specimens and evidence of intergranular fracture was observed in all specimens tested at 1 × 10−4 and 1 × 10−5 s−1. In tension-tested specimens, there was no difference in the mean grain diameter between doped and undoped specimens at 25% strain.

The mechanical properties of glacier beds play a fundamental role in regulating the sensitivity of glaciers to environmental forcing across a wide range of timescales. Glaciers are commonly underlain by deformable till whose mechanical properties and influence on ice flow are not well understood but are critical for reliable projections of future glacier states. Using synoptic-scale observations of glacier motion in different seasons to constrain numerical ice flow models, we study the mechanics of the bed beneath Hofsjökull, a land-terminating ice cap in central Iceland. Our results indicate that the bed deforms plastically and weakens following incipient summertime surface melt. Combining the inferred basal shear traction fields with a Coulomb-plastic bed model, we estimate the spatially distributed effective basal water pressure and show that changes in basal water pressure and glacier accelerations are non-local and non-linear. These results motivate an idealized physical model relating mean basal water pressure and basal slip rate wherein the sensitivity of glacier flow to changes in basal water pressure is inversely related to the ice surface slope.

Canada's Queen Elizabeth Islands (QEI) contain ~14% of the world's glacier and ice-cap area. Sparse in-situ measurements indicate that interannual variability in glacier surface mass balance in this region is driven primarily by variations in summer melt, and that the annual surface mass balance of four index glaciers has become increasingly negative since 2007. Here, we use a 16-a record of satellite-derived mean summer (June–August) land surface temperatures (LST) from NASA's Moderate Resolution Imaging Spectroradiometer to investigate large-scale spatial and temporal variability in the duration and intensity of summer melt across glaciated surfaces in the QEI from 2000 to 2015. During this period, QEI mean summer glacier surface temperatures increased at an average rate of 0.06 ± 0.04°C a−1, for a total of nearly 1°C. Most of this increase occurred between 2005 and 2012, when mean summer near-surface (2 m) and upper-air (700 hPa) temperatures were 1.0–1.2°C higher than the 1948–2015 mean. There is a strong correlation between the glacier LST and 700 hPa air temperature records (r> 0.8). The period 2005–12, when mean summer LSTs were anomalously high, was likely the warmest period in the region since at least 1948.