We present novel results of pollen analyses performed on a 10 m firn core retrieved from Alto dell’Ortles glacier (3840 m a.s.l.), Eastern Italian Alps, in 2009. The objective was to identify and quantify pollen grains retained in the ice to detect annual and interannual variations in the pollen spectra, thus enabling construction of an accurate pollen-based timescale. Up to now, this has been achieved by pollen diagram interpretation. Here we present a statistical approach developed to extract the seasonal/annual signal contained in the pollen spectra of an ice core. The method is based on principal component analyses of pollen assemblages obtained by high-level taxonomical identification. We apply this approach to the Ortles samples, demonstrating that seasonal and yearly variations of the pollen spectra are easily detectable and provide valuable information that can help improve the chronological model of the firn core. This approach can potentially be used for deeper cores as well as other types of archives (e.g. varved sediments), allowing faster, more objective estimation of yearly and seasonal variations than with classical percentage pollen diagrams.

In this paper we undertake a quantitative analysis of the dynamic process by which ice underneath a dry porous debris layer melts. We show that the incorporation of debris-layer airflow into a theoretical model of glacial melting can capture the empirically observed features of the so-called Østrem curve (a plot of the melt rate as a function of debris depth). Specifically, we show that the turning point in the Østrem curve can be caused by two distinct mechanisms: the increase in the proportion of ice that is debris-covered and/or a reduction in the evaporative heat flux as the debris layer thickens. This second effect causes an increased melt rate because the reduction in (latent) energy used for evaporation increases the amount of energy available for melting. Our model provides an explicit prediction for the melt rate and the temperature distribution within the debris layer, and provides insight into the relative importance of the two effects responsible for the maximum in the Østrem curve. We use the data of Nicholson and Benn (2006) to show that our model is consistent with existing empirical measurements.

Flowing snow is a cohesive granular material. Snow temperature and moisture content control the strength of the cohesive bonding between granules and therefore the outcome of granular interactions. Strong, cohesive interactions reduce the free mechanical energy in the avalanche core and therefore play a significant role in defining the avalanche flow regime. We introduce cohesion into avalanche dynamics model calculations by (1) treating cohesion as an additional internal binding energy that must be overcome to expand the avalanche flow volume, (2) modifying the Coulomb stress function to account for the increase in shear because of cohesive interactions and (3) increasing the activation energy to control the onset of avalanche fluidization. The modified shear stress function is based on force measurements in chute experiments with flowing snow. Example calculations are performed on ideal and real terrain to demonstrate how snow cohesion modifies avalanche flow and runout behaviour.

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

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.

Ice cliffs have been identified as a reason for higher ablation rates on debris-covered glaciers than are implied by the insulation effects of the debris. This study aims to improve our understanding of cliff backwasting, and the role of radiative fluxes in particular. An energy-balance model is forced with new data gathered in May and October 2013 on Lirung Glacier, Nepalese Himalaya. Observations show substantial variability in melt between cliffs, between locations on any cliff and between seasons. Using a high-resolution digital elevation model we calculate longwave fluxes incident to the cliff from surrounding terrain and include the effect of local shading on shortwave radiation. This is an advance over previous studies, that made simplified assumptions on cliff geometry and radiative fluxes. Measured melt rates varied between 3.25 and 8.6 cm d−1 in May and 0.18 and 1.34 cm d−1 in October. Model results reproduce the strong variability in space and time, suggesting considerable differences in radiative fluxes over one cliff. In October the model fails to reproduce stake readings, probably due to the lack of a refreezing component. Disregarding local topography can lead to overestimation of melt at the point scale by up to ∼9%.

We apply cross-correlation to Pléiades satellite images to generate accurate, high-resolution monthly surface velocity maps of Monte Tronador glaciers between March and June 2012. Measured surface displacements cover periods as short as 19 days, with a precision of ∼0.58 m (11 m a−1). These glaciers follow a radial flow pattern, with maximum surface speeds of ∼390 m a−1 associated with steep icefalls. The lower reaches of the debris-covered tongues of Verde and Casa Pangue glaciers are almost stagnant, whereas Ventisquero Negro, another debris-covered glacier, shows acceleration at the front due to calving into a proglacial lake. Low-elevation debris-covered glacier tongues show increasing velocities at the beginning of the accumulation season, whereas higher-elevation, clean-ice tongues reduce their speed during this period. This contrasting behavior is probably in response to an increase in water input to the subglacial system from winter rainfall events at low elevations and a decrease in meltwater production at higher elevations. These sequential velocity maps can help to identify the controls on glacier surface velocity, aid in the delimitation of ice divides and could also contribute to more realistic calibration of ice-flux-mass–balance models in this glacierized area.

Cryo-Raman tomography allows us, for the first time, to determine accurate morphologies and volumes of natural air hydrates in Antarctic ice cores. The measurements show complex growth structures that are not accounted for by the available models of hydrate formation.

Multi-annual records of glacier surface meteorology and energy balance are necessary to resolve glacier–climate interactions but remain sparse, especially in the Southern Hemisphere. To address this, we present a record from the ablation zone of Brewster Glacier, New Zealand, between October 2010 and September 2012. The mean air temperature was 1.2°C at 1760 m a.s.l., with only a moderate temperature difference between the warmest and coldest months (∼8°C). Long-term annual precipitation was estimated to exceed 6000 mm a−1, with the majority of precipitation falling within a few degrees of the freezing level. The main melt season was between November and March (83% of annual ablation), but melt events occurred during all months. Annually, net radiation was positive (a source of energy) and supplied 64% of the melt energy, driven primarily by net shortwave radiation. Net longwave radiation was often positive during cloudy conditions in summer, demonstrating the radiative importance of clouds during melt. Turbulent sensible and latent heat fluxes were directed towards the surface in the summer months, accounting for just over a third of the energy for melt (34%). The energy gain associated with rainfall was small except during heavy events in summer.

Ice-sheet models can be used to forecast ice losses from Antarctica and Greenland, but to fully quantify the risks associated with sea-level rise, probabilistic forecasts are needed. These require estimates of the probability density function (PDF) for various model parameters (e.g. the basal drag coefficient and ice viscosity). To infer such parameters from satellite observations it is common to use inverse methods. Two related approaches are in use: (1) minimization of a cost function that describes the misfit to the observations, often accompanied by explicit or implicit regularization, or (2) use of Bayes’ theorem to update prior assumptions about the probability of parameters. Both approaches have much in common and questions of regularization often map onto implicit choices of prior probabilities that are made explicit in the Bayesian framework. In both approaches questions can arise that seem to demand subjective input. One way to specify prior PDFs more objectively is by deriving transformation group priors that are invariant to symmetries of the problem, and then maximizing relative entropy, subject to any additional constraints. Here we investigate the application of these methods to the derivation of priors for a Bayesian approach to an idealized glaciological inverse problem.

This study shows how a new generation of terrestrial laser scanners can be used to investigate glacier surface ablation and other elements of glacial hydrodynamics at exceptionally high spatial and temporal resolution. The study area is an Alpine valley glacier, Haut Glacier d’Arolla, Switzerland. Here we use an ultra-long-range lidar RIEGL VZ-6000 scanner, having a laser specifically designed for measurement of snow- and ice-cover surfaces. We focus on two timescales: seasonal and daily. Our results show that a near-infrared scanning laser system can provide high-precision elevation change and ablation data from long ranges, and over relatively large sections of the glacier surface. We use it to quantify spatial variations in the patterns of surface melt at the seasonal scale, as controlled by both aspect and differential debris cover. At the daily scale, we quantify the effects of ogive-related differences in ice surface debris content on spatial patterns of ablation. Daily scale measurements point to possible hydraulic jacking of the glacier associated with short-term water pressure rises. This latter demonstration shows that this type of lidar may be used to address subglacial hydrologic questions, in addition to motion and ablation measurements.

To understand the dynamics of ice shelves, a knowledge of their internal and basal structure is very important. As the capacity to perform local surveys is limited, remote sensing provides an opportunity to obtain the relevant information. We must prove, however, that the relevant information can be obtained from remote sensing of the surface. That is the aim of this study. The Jelbart Ice Shelf, Antarctica, exhibits a variety of surface structures appearing as stripe-like features in radar imagery. We performed an airborne geophysical survey across these features and compared the results to TerraSAR-X imagery. We find that the stripe-like structures indicate surface troughs coinciding with the location of basal channels and crevasse-like features, revealed by radio-echo sounding. HH and VV polarizations do not show different magnitude. In surface troughs, the local accumulation rate is larger than at the flat surface. Viscoelastic modelling is used to gain an understanding of the surface undulations and their origin. The surface displacement, computed with a Maxwell model, matches the observed surface reasonably well. Our simulations show that the surface troughs develop over decadal to centennial timescales.

In this study we compare the anisotropic flow relations for polycrystalline ice of Azuma and Goto-Azuma (1996), Thorsteinsson (2002), Placidi and others (2010) and Budd and others (2013). Observations from the Dome Summit South (DSS) ice-coring site at Law Dome, East Antarctica, are used to model the vertical distribution of deviatoric stress components at the borehole site. The flow relations in which the anisotropic rheology is parameterized by a scalar function, so that the strain-rate and deviatoric stress tensor components are collinear, provide simple shear and vertical compression deviatoric stress profiles that are most consistent with laboratory observations of tertiary creep in combined stress configurations. Those flow relations where (1) the anisotropy is derived from the magnitude of applied stresses resolved onto the basal planes of individual grains and (2) the macroscopic deformation is obtained via homogenization of individual grain responses provide stress estimates less consistent with laboratory observations. This is most evident in combined simple shear and vertical compression flow regimes where shear is dominant. Our results highlight the difficulties associated with developing flow relations which incorporate a physically based description of microdeformation processes. In particular, this requires that all relevant microdeformation, recrystallization and recovery processes are adequately parameterized.

We report an in situ measurement of the electric field attenuation length Lα at radio frequencies for the bulk ice at Summit Station, Greenland, made by broadcasting radio-frequency signals vertically through the ice and measuring the relative power in the return ground bounce signal. We find the depth-averaged field attenuation length to be at 75 MHz. While this measurement has clear radioglaciological applications, the radio clarity of the ice also has implications for the detection of ultra-high energy (UHE) astrophysical particles via their radio emission in dielectric media such as ice. Assuming a reliable extrapolation to higher frequencies, the measured attenuation length at Summit Station is comparable to previously measured radio-frequency attenuation lengths at candidate particle detector sites around the world, and strengthens the case for Summit Station as a promising northern site for UHE neutrino detection.