We present a study of the development of a thermo-erosional notch at the waterline and its influence on calving of Hansbreen, a medium-sized grounded tidewater glacier in southern Svalbard. The study is based on the results of undercut notch melt modelling, based on measurements of sea-water properties, repeated terrestrial laser scans and analysis of time-lapse camera images. There is a strong correlation between observed calving activity and modelled melt rate of the undercut notch. Measured depths of the undercut reach 4 m and vary greatly over time. The calving activity of Hansbreen was significantly lower in 2011 than in 2012, due to the persistent presence of the ice pack in Hornsund fjord, which cooled the sea surface and suppressed the wave action. Calving on Hansbreen is controlled by a local imbalance of forces at the front, due to thermo-erosional undercutting at the sea waterline. Calving activity is therefore sensitive to changes in sea-water temperature and wave height. It may be expected that calving rates will rise with increased advection of warm oceanic water to the Arctic.

We present a set of new volume scaling relationships specific to Svalbard glaciers, derived from a sample of 60 volume–area pairs. Glacier volumes are computed from ground-penetrating radar (GPR)-retrieved ice thickness measurements, which have been compiled from different sources for this study. The most precise scaling models, in terms of lowest cross-validation errors, are obtained using a multivariate approach where, in addition to glacier area, glacier length and elevation range are also used as predictors. Using this multivariate scaling approach, together with the Randolph Glacier Inventory V3.2 for Svalbard and Jan Mayen, we obtain a regional volume estimate of 6700 ± 835 km3, or 17 ± 2 mm of sea-level equivalent (SLE). This result lies in the mid- to low range of recently published estimates, which show values as varied as 13 and 24 mm SLE. We assess the sensitivity of the scaling exponents to glacier characteristics such as size, aspect ratio and average slope, and find that the volume of steep-slope and cirque-type glaciers is not very sensitive to changes in glacier area.

The subglacial chemical weathering environment is largely controlled by low temperatures and the presence of freshly comminuted minerals with a high surface area. These characteristics are believed to promote dissolution processes that give rise to low silica and high Ca2+ fluxes emanating from glacierized basins. We test an alternative hypothesis, that mineral precipitation reactions in the subglacial environment play an equally important role in controlling the water chemistry in glacierized basins. We analyze borehole and proglacial water chemistry from a subarctic polythermal glacier, complemented by mineral XRD analysis of suspended sediment, till and bedrock samples. In conjunction with a thermodynamic analysis of the water and mineral chemistry, we use reaction-path modelling to study the chemical enrichment of water through the glacier system. We find that the high pH of the subglacial environment is conducive to secondary mineral precipitation, and that it is not possible to balance the water chemistry using dissolution reactions alone. We show that low silica can be explained by standard weathering reactions without having to invoke mineral-leaching reactions. Our results suggest that subglacial weathering intensity may be significantly underestimated if the production of secondary minerals is not considered.

Runoff estimation in high-altitude glacierized basins is an important issue on the Tibetan Plateau. To investigate glacier mass balance, runoff and water balance in the Qugaqie basin and Zhadang sub-basin in the southern Tibetan Plateau, two glacier models and three snow models were integrated into the spatially distributed hydrological model JAMS/J2K. The results showed that the temperature index method simulated glacier runoff better than the degree-day factor method. The simulated glacier melt volume in the Qugaqie basin in 2006, 2007 and 2008 contributed 58%, 50% and 41%, respectively, to its total runoff. In the Zhadang basin, the glacier melt volume contributed 78% and 66% to its runoff during 2007 and 2008, respectively. Compared with the observation results, the simulated glacier mass balance showed similar variations with slightly higher values, indicating an underestimation of glacier melt volume. The water balance simulation in the upstream areas (705–874 mm) was comparable to that in the downstream areas (1051–1502 mm) and generally lower than the observed results. In both basins, the glacier mass-balance simulation was relatively accurate in the melt season compared to the other seasons.

Effects of damage on elastic properties were studied in columnar-grained specimens of freshwater and saline ice, subjected, at −10°C, to varying levels of inelastic strain. The ice was compressed uniaxially at constant strain rates up to 0.20 strain, which caused localized recrystallization and imparted damage in the form of non-propagating cracks. Damage was quantified in terms of dimensionless crack density, which, along with recrystallized area fraction, was determined from thin sections. The change in porosity due to stress-induced cracks served as another indicator of damage. Elastic properties were derived using P-wave and S-wave ultrasonic transmission velocities measured in across-column directions through the damaged ice, either parallel (x1) or perpendicular (x2) to the initial loading direction. In general, as damage increased with greater strain, the ice became more compliant and (particularly freshwater ice) more anisotropic. Furthermore, with increasing strain rate, the magnitude of these effects and crack density tended to increase, in contrast to the recrystallized area fraction, which tended to decrease. We observed compliance to correspond closely with porosity and with dimensionless crack density, for strains up to 0.10. At greater levels of strain these correspondences became less clear due, in part, to the different character of the damage.

Short-term glacier velocity variations typically occur when a water input is accommodated by an increase in the subglacial water pressure. Although these velocity variations have been well documented on many glaciers, few studies have considered them on glaciers where heavy rain and glacier melt occur year-round. This study investigates the relationship between water inputs and glacier velocity on Franz Josef Glacier, New Zealand. We installed six GNSS stations across the lower glacier during austral summer 2010/11 and one station during summer 2012/13. Glacier velocity remained elevated at all stations for ∼7 days following large rain events. During diurnal melt events, we find velocity variations in the early afternoon (12:00–16:00) at 600 m a.s.l. and in the late evening (20:00–01:00) at 400 m a.s.l. We hypothesize that the late-evening velocity variations occurred as an upstream region of high subglacial water pressures and accelerated ice motion propagated downstream. This mechanism may also explain the increased longitudinal compression and transverse extension across the lower glacier during speed-up events. Our results indicate that the subglacial drainage system likely decreases in efficiency upstream and that the water input variability can still cause short-term velocity variations despite the large year-round water inputs.

Terrestrial radar interferometry (TRI) is a new technique for studying ice motion and volume change of glaciers. TRI is especially useful for temporally and spatially dense measurements of highly dynamic glacial termini. We conducted a TRI survey of Breiðamerkurjökull, a marine-terminating glacier in Iceland, imaging its terminus near the end of the melt season in 2011, 2012 and 2013. The ice velocities were as high as 5 m d−1, with the fastest velocities near the calving front. Retreat of the glacier over the 3 year observation period was accompanied by strong embayment formation. Iceberg tracking with the radar shows high current velocities near the embayment, probably indicating strong meltwater outflow and mixing with relatively warm lagoon water.

Iceberg calving and basal melting are the two primary mass loss processes from the Antarctic ice sheet, accounting for approximately equal amounts of mass loss. Basal melting under ice shelves has been increasingly well constrained in recent work, but changes in iceberg calving rates remain poorly quantified. Here we examine the processes that precede iceberg calving, and focus on initiation and propagation of ice-shelf rifts. Using satellite imagery from the Moderate Resolution Imaging Spectroradiometer (MODIS) and the Multi-angle Imaging Spectroradiometer (MISR), we monitored five active rifts on the Amery Ice Shelf, Antarctica, from 2002 to 2014. We found a strong seasonal component: propagation rates were highest during (austral) summer and nearly zero during winter. We found substantial variability in summer propagation rates, but found no evidence that the variability was correlated with large-scale environmental drivers, such as atmospheric temperature, winds or sea-ice concentration. We did find a positive correlation between large propagation events and the arrival of tsunamis in the region. The variability appears to be related to visible structural boundaries within the ice shelf, e.g. suture zones or crevasse fields. This suggests that a complete understanding of rift propagation and iceberg calving needs to consider local heterogeneities within an ice shelf.

We report on the radiative transfer process and optical properties of sea ice in the thermal infrared (TIR) band, presenting two new linear kernel driver models (Relative Emissivity Distribution Function, REDF) that describe TIR emission characteristics of smooth and rough ice. In order to test the models and determine the necessary coefficients, in situ measurements from the Bohai Sea were carried out during the 2011/12 and 2012/13 boreal winters. The results show that the relative emissivity of smooth sea ice decreases along with increasing viewing zenith angle, and the shape of the relative emissivity curve is similar to that of an ideal plane. Affected by parameters such as roughness and surface temperature distribution, the anisotropy of relative emissivity of sea ice with a high degree of roughness is stronger relative to the cosine emitter. The model coefficients were also obtained using a robust regression method based on the measured data. The presented models are more practical than the numerical radiative transfer model and can be used for multi-angular TIR remote sensing.

The ApRES (autonomous phase-sensitive radio-echo sounder) instrument is a robust, lightweight and relatively inexpensive radar that has been designed to allow long-term, unattended monitoring of ice-shelf and ice-sheet thinning. We describe the instrument and demonstrate its capabilities and limitations by presenting results from three trial campaigns conducted in different Antarctic settings. Two campaigns were ice sheet-based – Pine Island Glacier and Dome C – and one was conducted on the Ross Ice Shelf. The ice-shelf site demonstrates the ability of the instrument to collect a time series of basal melt rates; the two grounded ice applications show the potential to recover profiles of vertical strain rate and also demonstrate some of the limitations of the present system.

We infer the horizontal velocity fields of the ice caps Langjökull and Hofsjökull, central Iceland, using repeat-pass interferometric synthetic aperture radar (InSAR). NASA’s uninhabited aerial vehicle synthetic aperture radar (UAVSAR) acquired airborne InSAR data from multiple vantage points during the early melt season in June 2012. We develop a Bayesian approach for inferring three-dimensional velocity fields from multiple InSAR acquisitions. The horizontal components generally agree with available GPS measurements wherever ice motion is well constrained by InSAR observations. We provide evidence that changes in volumetric moisture content near the glacier surface induce phase offsets that obfuscate the vertical component of the surface velocity fields, an effect that could manifest itself on any glacier that experiences surface melt. Spatial patterns in the InSAR-derived horizontal speeds are broadly consistent with the results of a simple viscous flow model, and the directionality of the InSAR-derived horizontal flow field is nearly everywhere consistent with the ice surface gradient. Significant differences between the InSAR-derived horizontal speed and the speed predicted by the viscous flow model suggest that basal slip accounts for more than half the observed outlet glacier flow.

Due to solute impurities and freezing-point depression in polycrystalline ice, a complicated and dynamic network of liquid water forms within the solid ice matrix at the boundaries between ice crystal grains. Impurity concentrations, temperature and pressure influence this network structure and impact physical, transport and rheological properties of ice. However, the nature of this internal network structure is not fully understood. Here we utilize nuclear magnetic resonance (NMR) measurements of diffusion and magnetic relaxation to study the geometry and interconnectivity of the liquid-filled network in laboratory ice, formed from a 7 g L−1 NaCl solution, and its evolution due to recrystallization processes. Additionally, we apply these NMR measurements to observe the impact on ice microstructure of an ice-binding protein (IBP) excreted by the V3519-10 organism (Flavobacteriaceae family) isolated from the Vostok ice core in Antarctica. Recrystallization inhibition was observed as a function of IBP concentration. This work demonstrates the utility of advanced NMR techniques for applications to ice microstructure and has broader implications for understanding geophysical properties of cryospheric systems.

In glacial environments particle-size analysis of moraines provides insights into clast origin, transport history, depositional mechanism and processes of reworking. Traditional methods for grain-size classification are labour-intensive, physically intrusive and are limited to patch-scale (1 m2) observation. We develop emerging, high-resolution ground- and unmanned aerial vehicle-based ‘Structure-from-Motion’ (UAV-SfM) photogrammetry to recover grain-size information across a moraine surface in the Heritage Range, Antarctica. SfM data products were benchmarked against equivalent datasets acquired using terrestrial laser scanning, and were found to be accurate to within 1.7 and 50 mm for patch- and site-scale modelling, respectively. Grain-size distributions were obtained through digital grain classification, or ‘photo-sieving’, of patch-scale SfM orthoimagery. Photo-sieved distributions were accurate to <2 mm compared to control distributions derived from dry-sieving. A relationship between patch-scale median grain size and the standard deviation of local surface elevations was applied to a site-scale UAV-SfM model to facilitate upscaling and the production of a spatially continuous map of the median grain size across a 0.3 km2 area of moraine. This highly automated workflow for site-scale sedimentological characterization eliminates much of the subjectivity associated with traditional methods and forms a sound basis for subsequent glaciological process interpretation and analysis.

Despite the important role of supraglacial debris in ablation, knowledge of debris thickness on Himalayan glaciers is sparse. A recently developed method based on reanalysis data and thermal band satellite imagery has proved to be potentially suitable for debris thickness estimation without the need for detailed field data. In this study, we further develop the method and discuss possibilities and limitations arising from its application to a glacier in the Himalaya with scarce in situ data. Surface temperature patterns are consistent for 13 scenes of Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER) and Landsat 7 imagery and correlate well with incoming shortwave radiation and air temperature. We use an energy-balance approach to subtract these radiation or air temperature effects, in order to estimate debris thickness patterns as a function of surface temperature. Both incoming shortwave and longwave radiation are estimated with reasonable accuracy when applying parameterizations and reanalysis data. However, the model likely underestimates debris thickness, probably due to incorrect representation of vertical debris temperature profiles, the rate of heat storage and turbulent sensible heat flux. Moreover, the uncertainty of the result was found to increase significantly with thicker debris, a promising result since ablation is enhanced by thin debris of 1–2 cm.

Snow avalanche simulation software is a commonly used tool for hazard estimation and mitigation planning. In this study a depth-averaged flow model, combining a simple entrainment and friction relation, is implemented in the software SamosAT. Computational results strongly depend on the simulation input, in particular on the employed model parameters. A long-standing problem is to quantify the influence of these parameters on the simulation results. We present a new multivariate optimization approach for avalanche simulation in three-dimensional terrain. The method takes into account the entire physically relevant range of the two friction parameters (Coulomb friction, turbulent drag) and one entrainment parameter. These three flow model parameters are scrutinized with respect to six optimization variables (runout, matched and exceeded affected area, maximum velocity, average deposition depth and mass growth). The approach is applied to a documented extreme avalanche event, recorded in St Anton, Austria. The final results provide adjusted parameter distributions optimizing the simulation–observation correspondence. At the same time, the degree of parameter–variable correspondence is determined. We show that the specification of optimal values for certain model parameters is near-impossible, if corresponding optimization variables are neglected or unavailable.

In this study, we analyze a large dataset of seismic signals, recorded by station TROLL in Dronning Maud Land, Antarctica. The signals, recorded in April–December 2012, came from sources near the edge of the ice shelves, at distances of 230–500 km from TROLL. The sources, which moved westward with time, could be associated with four large, tabular icebergs, drifting between 15° E and 8° W. Combining the seismological data with information from satellite remote sensing, we find that one-third of the signals can be attributed to individual icebergs. The trajectories of three of the associated icebergs are known through iceberg-tracking databases, whereas the fourth, a fragment of one of the other three, is untracked, and only scarce information is available from satellite imagery. The observed seismic signals exhibit a wide variety of frequency characteristics, from unstructured episodes to occurrences of iceberg harmonic tremor. Although we are not able to determine the exact cause of the signals, we classify them into five classes on a phenomenological basis. This study demonstrates the potential of regional seismic networks for iceberg monitoring as supplementary resources to information obtained with remote-sensing technologies.

A model subglacial drainage system, coupled to an ice-dammed reservoir that receives a time-varying meltwater input, is described and analysed. In numerical experiments an ice-marginal lake drains through a subglacial channel, producing periodic floods, and fills with meltwater at a rate dependent on air temperature, which varies seasonally with a peak value of Tm. The analysis reveals regions of Tm parameter space corresponding to ‘mode locking’, where flood repeat time is independent of Tm; resonance, where decreasing Tm counter-intuitively increases flood size; and chaotic dynamics, where flood cycles are sensitive to initial conditions, never repeat and exhibit phase-space mixing. Bifurcations associated with abrupt changes in flood size and timing within the year separate these regions. This is the first time these complex dynamics have been displayed by a glaciohydraulic model and these findings have implications for our understanding of ice-marginal lakes, moulins and subglacial lakes.

Lake-calving Yakutat Glacier in southeast Alaska, USA, is undergoing rapid thinning and terminus retreat. We use a simplified glacier model to evaluate its future mass loss. In a first step we compute glacier-wide mass change with a surface mass-balance model, and add a mass loss component due to ice flux through the calving front. We then use an empirical elevation change curve to adjust for surface elevation change of the glacier and finally use a flotation criterion to account for terminus retreat due to frontal ablation. Surface mass balance is computed on a daily timescale; elevation change and retreat is adjusted on a decadal scale. We use two scenarios to simulate future mass change: (1) keeping the current (2000–10) climate and (2) forcing the model with a projected warming climate. We find that the glacier will disappear in the decade before 2110 or 2070 under constant or warming climates, respectively. For the first few decades, the glacier can maintain its current thinning rates by retreating and associated loss of high-ablating, low-elevation areas. However, once higher elevations have thinned substantially, the glacier can no longer counteract accelerated thinning by retreat and mass loss accelerates, even under constant climate conditions. We find that it would take a substantial cooling of 1.5°C to reverse the ongoing retreat. It is therefore likely that Yakutat Glacier will continue its retreat at an accelerating rate and disappear entirely.

During the 2010/11 boreal winter, a distributed set of backscatter measurements was collected using a ground-based Ku-band (17.2 GHz) scatterometer system at 26 open tundra sites. A standard snow-sampling procedure was completed after each scan to evaluate local variability in snow layering, depth, density and water equivalent (SWE) within the scatterometer field of view. The shallow depths and large basal depth hoar encountered presented an opportunity to evaluate backscatter under a set of previously untested conditions. Strong Ku-band response was found with increasing snow depth and snow water equivalent (SWE). In particular, co-polarized vertical backscatter increased by 0.82 dB for every 1 cm increase in SWE (R2 = 0.62). While the result indicated strong potential for Ku-band retrieval of shallow snow properties, it did not characterize the influence of sub-scan variability. An enhanced snow-sampling procedure was introduced to generate detailed characterizations of stratigraphy within the scatterometer field of view using near-infrared photography along the length of a 5 m trench. Changes in snow properties along the trench were used to discuss variations in the collocated backscatter response. A pair of contrasting observation sites was used to highlight uncertainties in backscatter response related to short length scale spatial variability in the observed tundra environment.

Marine ice accretes at the base of ice shelves, often infilling open structural weaknesses and increasing ice-shelf stability. However, the timing and location of marine ice formation remain poorly understood. This study determines marine ice source water composition and origin by examining marine ice crystal morphology, water isotope and solute chemistry in ice samples collected from the southern McMurdo Ice Shelf (SMIS), Antarctica. The measured co-isotopic record together with the output of a freezing model for frazil crystals indicate a spatio-temporally varying water source of sea water and relatively fresher water, such as melted meteoric or marine ice. This is in agreement with the occurrence of primarily banded and granular ice crystal facies typical for frazil ice crystals that nucleate in a supercooled mixture of water masses. We propose that marine ice exposed at the surface of SMIS, which experiences summer melt, is routed to the ice-shelf base via the tide crack. Here frazil crystals nucleate in a double diffusion mechanism of heat and salt between two water masses at their salinity-dependent freezing point. Recycling of previously formed marine ice facilitates ice-shelf self-sustenance in a warming climate.