Two non-invasive optical Schlieren methods have been adapted to visualize brine channel development and convective processes in experimentally grown sea ice obtained when a NaCl aqueous solution is cooled from above in a quasi-two-dimensional Hele–Shaw cell. The two different visualization methods, i.e. traditional and synthetic Schlieren optical imaging, produce high spatial resolution images of transport processes during ice growth, without any external perturbation. These images allow observations of the flow dynamics simultaneously within the ice layer, around the ice/water interface, and in the liquid water layer, revealing connections between the processes occurring within the two phases. Results from these methods show that desalination of the growing ice layer occurs by two concurrent, yet independent, mechanisms: (1) boundary layer convection persisting throughout the ice growth period, with short fingers present just below the ice/water interface, and (2) gravity-driven drainage from the brine channels producing deep penetrating convective streamers, which appear after a given time from the beginning of ice growth. The improved visualization and qualitative characterization of these processes show that Schlieren optical methods have exciting potential applications for future study of convective processes during sea-ice growth.

Fractured systems developing at the base of temperate glaciers play a significant role in the transport of subglacial water. Their influence on dynamics and hydraulics of glaciers and the seasonal alterations they experience during the hydrological year, remain however, poorly understood, as they lay beyond the reach of most conventional glaciological techniques. In this work, we investigate the seasonal evolution of glacier seismicity in the ablation zone of Rhonegletscher (Switzerland) and focus on basal fracture processes. To do this, we use seismic measurements performed between June 2012 and July 2013 and complementary measurements of surface motion and changes in melt water input inferred from degree day modeling. Results of a spectrogram analysis highlight predominant seismic emissions until end of October 2012 and from mid April 2013 in line with measured surface velocities, and indicate a rapid adaptation of the basal drainage system at initiation and end of the winter period. The analysis of icequakes emitted from two active regions shows that the basal events, due to opening and closing of cracks within fracture networks, persist after winter initiation, and points out the year-long preservation of the associated sources. In this regard, we deduce that basal crack networks aseismically evolve during the winter period.

We offer the first sub-seasonal view of glacial age archives from the Siple Dome-A (SDMA) ice core using the ultra-high resolution capabilities of a newly developed laser ablation-inductively coupled plasma-mass spectrometry (LA-ICP-MS; 121 µm sampling resolution) system capable of conducting multi-element glaciochemical analysis. Our ultra-high resolution data demonstrates that: (1) the SDMA ice core record can be annually dated based on seasonality in chemical inputs at a depth not previously possible using previous glaciochemical sampling methods, (2) winter accumulation at the SD site was greater than summer accumulation during the three late glacial periods selected (~15.3, 17.3, 21.4 Ka ago) in this study and (3) resulting annual layer thicknesses results show greater variability than the current SD ice core depth/age model (Brook and others, 2005), possibly due to depositional effects such as wind scouring and/or decadal variability in snow accumulation that is not captured by the resolution of the current depth/age model.

Antarctic ice shelves are buttressed by numerous pinning points attaching to the otherwise freely-floating ice from below. Some of these kilometric-scale grounded features are unresolved in Antarctic-wide datasets of ice thickness and bathymetry, hampering ice flow models to fully capture dynamics at the grounding line and upstream. We investigate the role of an 8.7 km2 pinning point at the front of the Roi Baudouin Ice Shelf, East Antarctica. Using ERS interferometry and ALOS-PALSAR speckle tracking, we derive, on a 125 m grid spacing, surface velocities deviating by −5.2 ± 4.5 m a−1 from 37 on-site global navigation satellite systems-derived velocities. We find no evidence for ice flow changes on decadal time scales and we show that ice on the pinning point virtually stagnates, deviating the ice stream and causing enhanced horizontal shearing upstream. Using the BISICLES ice-flow model, we invert for basal friction and ice rigidity with three input scenarios of ice velocity and geometry. We show that inversion results are the most sensitive to the presence/absence of the pinning point in the bathymetry; surface velocities at the pinning point are of secondary importance. Undersampling of pinning points results in erroneous ice-shelf properties in models initialised by control methods. This may impact prognostic modelling for ice-sheet evolution in the case of unpinning.

Using a field emission scanning electron microscope and energy-dispersive X-ray spectrometer (FESEM−EDS), individual insoluble dust particles in ice cores recovered from the northern (Dunde), western (Muztagata) and south-eastern (Palong-Zangbo) Tibetan Plateau were analysed in order to reveal the mineralogical and morphological characteristics of the analogue of long-range transported Asian dust. The results reveal that the dust particles are mainly composed of quartz (17–36%, number abundance), clay (37–48%) and feldspar (12–18%). Illite and chlorite are the dominant clay species, while kaolinite is rarely observed. For the three sites, regional differences in the mineral assemblages are significant, in particular the abundance of quartz, chlorite and muscovite species, reflecting the regional provenance of these dust particles and the climatic regime in their source areas. Oxide ratios in clays indicate different weathering strengths of the particles in their source regions, with a higher K2O/(SiO2 + Al2O3) in the western and northern Tibetan Plateau. The individual particles have modes of 1.48 and 1.53 in the aspect ratio and circularity distribution respectively. No significant relationship between aspect ratio and circularity was found. Quartz and feldspar particles have a narrow aspect ratio distribution. Muscovite particles have the coarsest grain size, while chlorite particles have the finest.

We present the first snow/ice chemistry and ice radar results ever collected from South Georgia as part of an initial reconnaissance with the ultimate goal of assessing the feasibility of a South Georgia ice core to reconstruct past climate in the South Atlantic. South Georgia is well situated to capture a record of past atmospheric chemical composition over the South Atlantic and of past variability in the position and intensity of the austral westerlies. The question is how well preserved an ice core record can be recovered from a region experiencing accelerated melting? The results presented in this paper offer only a preliminary step in determining the feasibility of future deep ice coring on South Georgia. However, this initial reconnaissance does provide some basic information including: the chemistry of the atmosphere over South Georgia relative to other Southern Hemisphere ice coring sites; the potential for preservation of ‘annual layers’ in old ice on the island; a possible age for deep ice in the region; and an estimate of glacier health in the lower elevation regions of the island.

Glaciers in Greenland are changing rapidly. To better understand these changes, we have produced a series of seven synthetic aperture radar (SAR) backscatter mosaics for seven winters during the period from 2000 to 2013. Six of the mosaics were created using RADARSAT Fine-Beam data and the seventh used ALOS-PALSAR Fine-Beam Single-Polarization data. The RADARSAT mosaics are radiometrically calibrated and capture changes in the backscatter coefficient related to melt and other events, particularly the strong melting in the summer of 2012. Comparison of features in the ascending-orbit ALOS mosaic and the descending-orbit RADARSAT mosaics indicate that in areas of smooth to moderate topography their locations are consistent to within a few tens of meters. The locations of features identifiable in the RADARSAT mosaics, which were collected with the same imaging parameters, generally agree to within better than the 20 m posting of the data. With such geometric accuracy, these data establish a record of change in Greenland for the early part of the 21st century, thus providing a baseline that can be compared with new radar and optical datasets.

The mass flux from marine terminating glacier outlets from the Greenland ice sheet is a significant mass budget term. Yet, glacier calving sensitivity to climate parameters is not well understood. In this study, the front area changes of 42 marine terminating Greenland outlet glaciers were measured annually from 1999 to 2013 using optical satellite imagery. Thirty-five glaciers retreated, while seven remained stable. The average front position change was −240 m a−1 and the total rate of front area change was −114 km2 a−1. The yearly average area changes of all measured glaciers are compared with interannual variations in sea surface temperature (SST), sea-ice concentration (SIC), surface melt, Northern Hemisphere (NH) surface air temperatures and the North Atlantic Oscillation (NAO) index. We find: (1) northwestern glaciers area change have a significant correlation with surface melt and NH land temperature variations; (2) southeastern glaciers have a significant correlation with local SST, SIC and NH temperatures; and (3) southwestern glaciers area change correlate significantly with local SST and NAO index. We conclude that a climate sensitivity signal emerges considering a population of glaciers. Further, we find a significant correlation between the date of minimum extent and glacier latitude. Area change data are available from http://GAC.geus.dk.

I explore the tidewater glacier cycle with a 1-D, depth- and width-integrated flow model that includes a mass-flux calving parameterization. The parameterization is developed from mass continuity arguments and relates the calving rate to the terminus velocity and the terminus balance velocity. The model demonstrates variable sensitivity to climate. From an advanced, stable configuration, a small warming of the climate triggers a rapid retreat that causes large-scale drawdown and is enhanced by positive glacier-dynamic feedbacks. Eventually, the terminus retreats out of deep water and the terminus velocity decreases, resulting in reduced drawdown and the potential for restabilization. Terminus readvance can be initiated by cooling the climate. Terminus advance into deep water is difficult to sustain, however, due to negative feedbacks between glacier dynamics and surface mass balance. Despite uncertainty in the precise form of the parameterization, the model provides a simple explanation of the tidewater glacier cycle and can be used to evaluate the response of tidewater glaciers to climate variability. It also highlights the importance of improving parameterizations of calving rates and of incorporating sediment dynamics into tidewater glacier models.

The glacier system covering Europe's highest mountain, Elbrus, has exhibited an accelerated retreat since 1980. Some studies have related this retreat to a significant summer temperature increasing trend. Relief- and aspect-related parameters for the glacierized area have an important impact on glacier changes. In this paper, the changes in glacier area are identified, quantified and correlated with relief parameters for the period 1985–2007. Spatial analysis was performed using the GLAM-CD (Glacier Mapper – Change Detector) algorithm. The input data for this algorithm were Landsat 5 images, the Aster Global Digital Elevation and the glacier outlines from the GLIMS project (Global Land Ice Measurements from Space). Regression analyses between glacier area losses and relief-related parameters indicate a significant positive relation with the altitude and a significant negative relation with the glacier surface area. In this context, we used a correlated component regression to model these relations. The model explains >50% of the total variation.

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.

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.

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

Terrestrial time-lapse photography offers insight into glacial processes through high spatial and temporal resolution imagery. However, oblique camera views complicate measurement in geographic coordinates, and lead to reliance on specific imaging geometries or simplifying assumptions for calculating parameters such as ice velocity. We develop a novel approach that integrates time-lapse imagery with multitemporal DEMs to derive full three-dimensional coordinates for natural features tracked throughout a monoscopic image sequence. This enables daily independent measurement of horizontal (ice flow) and vertical (ice melt) velocities. By combining two terrestrial laser scanner surveys with a 73 days sequence from Sólheimajökull, Iceland, variations in horizontal ice velocity of ~10% were identified over timescales of ~25 days. An overall decrease of ~3.0 m surface elevation showed asynchronous rate changes with the horizontal velocity variations, demonstrating a temporal disconnect between the processes of ice surface lowering and mechanisms of glacier movement. Our software, ‘Pointcatcher’, is freely available for user-friendly interactive processing of general time-lapse sequences and includes Monte Carlo error analysis and uncertainty in projection onto DEM surfaces. It is particularly suited for analysis of challenging oblique glacial imagery, and we discuss good features to track, both for correction of camera motion and for deriving ice velocities.

C band backscatter parameters contain information about the upper snowpack/firn in the dry snow zone. The wide incidence angle diversity of the Advanced Scatterometer (ASCAT) gives unprecedented characterisation of backscatter anisotropy, revealing the backscatter response to climatic forcing. The A (isotropic component) and M2 (bi-sinusoidal azimuth anisotropy) parameters are investigated here, in conjunction with data from atmospheric and snowpack models, to identify the backscatter response to surface forcing parameters (wind speed and persistence, precipitation, surface temperature, density and grain size). The long-term mean A parameter is successfully recreated with a regression using these drivers, indicating strong links between the A parameter and precipitation on long timescales. While the ASCAT time series is too short to determine which factors drive observed trends, factors influencing the seasonal and short timescale variability are revealed. On these timescales, A strongly responds to the propagation of surface temperature cycles/anomalies downward through the firn, via direct modulation of the dielectric constant. The influence of precipitation on A is small at shorter timescales. The M2 parameter is controlled by wind speed and persistence, through modification of monodirectionally-aligned surface roughness. This variability indicates that throughout much of coastal Antarctica, a microwave ‘snapshot’ is generally not representative of longer-term conditions.

Near-surface air temperature is an important determinant of the surface energy balance of glaciers and is often represented by a constant linear temperature gradients (TGs) in models. Spatio-temporal variability in 2 m air temperature was measured across the debris-covered Miage Glacier, Italy, over an 89 d period during the 2014 ablation season using a network of 19 stations. Air temperature was found to be strongly dependent upon elevation for most stations, even under varying meteorological conditions and at different times of day, and its spatial variability was well explained by a locally derived mean linear TG (MG–TG) of −0.0088°C m−1. However, local temperature depressions occurred over areas of very thin or patchy debris cover. The MG–TG, together with other air TGs, extrapolated from both on- and off-glacier sites, were applied in a distributed energy-balance model. Compared with piecewise air temperature extrapolation from all on-glacier stations, modelled ablation, using the MG–TG, increased by <1%, increasing to >4% using the environmental ‘lapse rate’. Ice melt under thick debris was relatively insensitive to air temperature, while the effects of different temperature extrapolation methods were strongest at high elevation sites of thin and patchy debris cover.

We study the evolution of the Juneau Icefield, one of the largest icefields in North America (>3700 km2), using the Parallel Ice Sheet Model (PISM). We test two climate datasets: 20 km Weather Research and Forecasting Model (WRF) output, and data from the Scenarios Network for Alaska Planning (SNAP), derived from spatial interpolation of observations. Good agreement between simulated and observed surface mass balance was achieved only after substantially adjusting WRF precipitation to account for unresolved orographic effects, while SNAP's climate pattern is incompatible with observations of surface mass balance. Using the WRF data forced with the RCP6.0 emission scenario, the model projects a decrease in ice volume by 58–68% and a 57–63% area loss by 2099 compared with 2010. If the modeled 2070–99 climate is held constant beyond 2099, the icefield is eliminated by 2200. With constant 1971–2010 climate, the icefield stabilizes at 86% of its present-day volume. Experiments started from an ice-free state indicate that steady-state volumes are largely independent of the initial ice volume when forced by identical scenarios of climate stabilization. Despite large projected volume losses, the complex high-mountain topography makes the Juneau Icefield less susceptible to climate warming than low-lying Alaskan icefields.