Ice streams provide major drainage pathways for the Antarctic ice sheet. The stress distribution and style of flow in such ice streams produce elastic and rheological anisotropy, which informs ice-flow modelling as to how ice masses respond to external changes such as global warming. Here we analyse elastic anisotropy in Rutford Ice Stream, West Antarctica, using observations of shear-wave splitting from three-component icequake seismograms to characterize ice deformation via crystal-preferred orientation. Over 110 high-quality measurements are made on 41 events recorded at five stations deployed temporarily near the ice-stream grounding line. To the best of our knowledge, this is the first well-documented observation of shear-wave splitting from Antarctic icequakes. The magnitude of the splitting ranges from 2 to 80 ms and suggests a maximum of 6% shear-wave splitting. The fast shear-wave polarization direction is roughly perpendicular to ice-flow direction. We consider three mechanisms for ice anisotropy: a cluster model (vertical transversely isotropic (VTI) model); a girdle model (horizontal transversely isotropic (HTI) model); and crack-induced anisotropy (HTI model). Based on the data, we can rule out a VTI mechanism as the sole cause of anisotropy – an HTI component is needed, which may be due to ice crystal a-axis alignment in the direction of flow or the alignment of cracks or ice films in the plane perpendicular to the flow direction. The results suggest a combination of mechanisms may be at play, which represent vertical variations in the symmetry of ice crystal anisotropy in an ice stream, as predicted by ice fabric models.

Glaciers on King George Island, Antarctica, have shown retreat and surface lowering in recent decades, concurrent with increasing air temperatures. A large portion of the glacier perimeter is ocean-terminating, suggesting possible large mass losses due to calving and submarine melting. Here we estimate the ice discharge into the ocean for the King George Island ice cap. L-band synthetic aperture radar images covering the time-span January 2008 to January 2011 over King George Island are processed using an intensity-tracking algorithm to obtain surface velocity measurements. Pixel offsets from 40 pairs of radar images are analysed and inverted to estimate a weighted average surface velocity field. Ice thicknesses are derived from simple principles of ice flow mechanics using the computed surface velocity fields and in situ thickness data. The maximum ice surface speeds reach >225 m a-1, and the total ice discharge for the analysed flux gates of King George Island is estimated to be 0.720 ± 0.428 Gt a−1, corresponding to a specific mass loss of 0.64 ± 0.38 m w.e. a-1 over the area of the entire ice cap (1127 km2).

To investigate the spatio-temporal characteristics of sea-ice resource, we used sea-ice volume to measure the amount of sea-ice resource in the Bohai Sea, China. The sea-ice area was extracted from Advanced Very High Resolution Radiometer (AVHRR) remote-sensing images using the zonal threshold method. The sea-ice thickness was estimated using a sea-ice model based on shortwave radiation theory and field measurements. The spatio-temporal characteristics of sea-ice volume were then analysed using GIS technology. The results indicate that the Bohai Sea experienced two sea-ice volume peaks in winter 2009/10. The largest sea-ice volume was in Liaodong Bay (∼80.26% of the entire sea-ice volume of the Bohai Sea). Bohai Bay had the second largest ice volume, and Laizhou Bay the smallest. The relationship between sea-ice volume and distance from shore is essentially exponential. The proportion of total sea-ice volume that is 0–10 km from shore is ∼42.43%, whereas the proportion that is 100–110 km from shore is only 0.002%.

We have investigated the speed of compressional waves in a polythermal glacier by, first, predicting them from a simple three-phase (ice, air, water) model derived from a published ground-penetrating radar study, and then comparing them with field data from four orthogonally orientated walkaway vertical seismic profiles (VSPs) acquired in an 80 m deep borehole drilled in the ablation area of Storglaciären, northern Sweden. The model predicts that the P-wave speed increases gradually with depth from 3700ms–1 at the surface to 3760ms–1 at 80m depth, and this change is almost wholly caused by a reduction in air content from 3% at the surface to <0.5% at depth. Changes in P-wave speed due to water content variations are small (<10 ms–1); the model’s seismic cold–temperate transition surface (CTS) is characterized by a 0.3% decrease downwards in P-wave speed (about ten times smaller than the radar CTS). This lack of sensitivity, and the small contrast at the CTS, makes seismically derived water content estimation very challenging. Nevertheless, for down-going direct-wave first arrivals for zero- and near-offset VSP shots, we find that the model-predicted travel times and field observations agree to within 0.2 ms, i.e. less than the observational uncertainties.

Near-surface air temperature, typically measured at a height of 2 m, is the most important control on the energy exchange and the melt rate at a snow or ice surface. It is distributed in a simplistic manner in most glacier melt models by using constant linear lapse rates, which poorly represent the actual spatial and temporal variability of air temperature. In this paper, we test a simple thermodynamic model proposed by Greuell and Böhm in 1998 as an alternative, using a new dataset of air temperature measurements from along the flowline of Haut Glacier d’Arolla, Switzerland. The unmodified model performs little better than assuming a constant linear lapse rate. When modified to allow the ratio of the boundary layer height to the bulk heat transfer coefficient to vary along the flowline, the model matches measured air temperatures better, and a further reduction of the root-mean-square error is obtained, although there is still considerable scope for improvement. The modified model is shown to perform best under conditions favourable to the development of katabatic winds – few clouds, positive ambient air temperature, limited influence of synoptic or valley winds and a long fetch – but its performance is poor under cloudy conditions.

Snow and ice thickness in the coastal Kara Sea, Russian Arctic, were investigated by applying the thermodynamic sea-ice model HIGHTSI. The external forcing was based on two numerical weather prediction (NWP) models: the High Resolution Limited Area Model (HIRLAM) and the European Centre for Medium-Range Weather Forecasts (ECMWF) model. A number of model experiments were carried out applying different snow parameterization schemes. The modelled ice thickness was compared with in situ measurements and the modelled snow thickness was compared with the NASA Advanced Microwave Scanning Radiometer for Earth Observing System (AMSR-E) snow thickness. The HIRLAM and ECMWF model results agreed with each other on air temperature and wind. The NWP model precipitation forecasts caught up the synoptic-scale snowfall events, but the magnitude was liable to errors. The ice growth was modelled reasonably well applying HIGHTSI either with a simple parameterization for snow thickness or with the HIRLAM or ECMWF model precipitation as input. For the latter, however, an adjustment of snow accumulation in early winter was necessary to avoid excessive accumulation and consequent underestimation of ice thickness. Applying effective snow heat conductivity improved the modelled ice thickness. The HIGHTSI-modelled snow thickness had a seasonal evolution similar to that of the AMSR-E snow thickness. New field data are urgently needed to validate NWP and ice models and remote-sensing products for snow and sea ice in the Kara Sea.

The Patagonia Icefields are characterized by a large number of outlet glaciers calving into lakes and the ocean. In contrast to the recent intensive research activities on tidewater glaciers in other regions, very few observations have been made on calving glaciers in Patagonia. We analysed satellite images of Glaciar Upsala, the third largest freshwater calving glacier in the Southern Patagonia Icefield, to investigate changes in its front position, ice velocity and surface elevation from 2000 to 2011. Our analyses revealed a clear transition from a relatively stable phase to a rapidly retreating and fast-flowing condition in 2008. The glacier front receded by 2.9 km, and the ice velocity increased by 20–50%, over the 2008–11 period. We also found that the ice surface lowered at a rate of up to 39 m a−1 from 2006 to 2010. This magnitude and the rate of changes in the glacier front position, ice velocity and surface elevation are greater than previously reported for Glaciar Upsala, and comparable to recent observations of large tidewater glaciers in Greenland. Our data illustrate details of a rapidly retreating calving glacier in Patagonia that have been scarcely reported despite their importance to the mass budget of the Patagonia Icefields.

The environments underneath ice sheets are of high scientific interest. Wireless sensors offer the prospect of sustained, distributed remote sensing in the subglacial environment. Typically, wireless sensor networks use radio-frequency (RF) electromagnetic communications, but these are highly attenuated in wet environments. In such environments, acoustic communications may be more power-efficient. Here we review the literature on acoustic and RF attenuation through ice and other relevant media, and present the results of new experiments on acoustic attenuation in glacial ice. Link budgets for communications from a range of subglacial environments show that acoustic communications are a viable strategy for transmission through water and ice where RF is too highly attenuated to be detected. Acoustic communication at 30 kHz is predicted to be possible through 1 km of glacial ice, using a 1 W transmitter. Such a strategy may be appropriate for shallow ice-stream environments around the Antarctic and Greenland ice sheet margins.

Amplitude-versus-angle (AVA) analysis of a seismic reflection line, imaged 13 km from Russell Glacier terminus, near the western margin of the Greenland ice sheet (GrIS), suggests the presence of sediment at the bed. The analysis was complicated by the lack of identifiable multiples in the data due to a highly irregular and crevassed ice surface, rendering deeper seismic returns noisy. A modified technique for AVA processing of glacial seismic data using forward modelling with primary reflection amplitudes and simulated multiple amplitudes is presented here. Our analysis demonstrates that AVA analysis can be applied to areas with noisy seismic returns and indicates that sediment underlies the seismic study site. Our data are inconsistent with the common assumption that the GrIS is underlain only by hard bedrock, but consistent with the presence of subglacial sediment with porosity between 30% and 40%. As analysis and modelling of ice-sheet dynamics requires a sound knowledge of the underlying basal materials, subglacial sediment should be taken into account when considering ice dynamics in this region of the GrIS.

In this paper, a molecular-dynamics sea-ice model is used to study contact and force networks in fragmented sea ice, composed of separate floes with power-law size distribution. The momentum equations for individual floes, taking into account floe/floe collisions (with Hertzian contact mechanics), are formulated in a way suitable for a computationally efficient numerical algorithm, allowing simulation of systems of thousands of floes. The simulations are performed for a number of scenarios: pure convergence without wind, through a jamming phase transition; constant wind at a constant ice concentration; and an idealized marginal ice zone. An analysis of the statistical properties of the contact and force networks reveals a highly localized, intermittent character of internal stress in the ice, as well as the role of the size-dependent response of floes to the forcing in formation of spatial patterns of internal stress at lower ice concentrations. The results provide a valuable starting point for formulating improved rheology models for fragmented sea ice.

We estimate temperature sensitivity of observed mass balance for glaciers in the European Alps to compare with our earlier estimates from modelling. We treat 1961-90 as our reference period and compare mass balance for 1991-2010 with this reference period. There are eight Alpine glaciers with long enough records. The mean mass balance for the eight glaciers is close to zero for the reference period, very negative for 1991-2010 and highly negative in the exceptionally warm year 2003. Mass-balance changes are compared with 1991-2010 temperature anomalies to give an average temperature sensitivity of approximately -0.7 ± 0.2 m w.e. a-1 K-1 for the eight glaciers. There is considerable variation between individual glaciers, with higher sensitivity for western glaciers and lower sensitivity for eastern glaciers. Temperature sensitivity for the eight glaciers is strongly correlated with the standard deviation of observed mass balance, which from earlier work is known to be correlated with winter balance. The mean temperature sensitivity for observed mass balance is in reasonable agreement with values from earlier models, but changes in glacier hypsometry during 1991-2010 may have further increased observed mass loss.

We report the results of a comparative analysis focusing on grain size, mineralogical composition and spectral reflectance values (400-2500 nm) of cryoconite samples collected from Jakobshavn Isbræ, West Greenland, and Canada Glacier, McMurdo Dry Valleys, Antarctica. The samples from the Greenland site were composed of small particles clumped into larger rounded agglomerates, while those from the site in Antarctica contained fragments of different sizes and shapes. Mineralogical analysis indicates that the samples from Jakobshavn Isbræ contained a higher percentage of quartz and albite, whereas those from Canada Glacier contained a higher percentage of amphibole, augite and biotite. Spectral measurements confirmed the primary role of organic material in reducing the reflectance over the measured spectrum. The reflectance of the samples from the Antarctic site remained low after the removal of organic matter because of the higher concentration of minerals with low reflectance. The reflectance of dried cryoconite samples in the visible region was relatively low (e.g. between ∼0.1 and ∼0.4) favouring increased absorbed solar radiation. Despite high reflectance values in the shortwave infrared region, the effect of the presence of cryoconite is negligible at infrared wavelengths where ice reflectance is low.

Active seismic imaging of glaciers and ice sheets is important for constraining inputs to climate models, such as englacial ice fabric and the nature of the basal interface. However, acquiring high-quality seismic data is time-consuming and resource-intensive. Using traditional single-element geophones requires ideal weather conditions (e.g. light winds) and excellent source coupling. In addition, deploying and retrieving these geophones is slow and cumbersome. We have developed a four-element ‘georod’ that enhances signal levels by 20–30dB in a variety of conditions, including blowing snow and poorly coupled source detonations. The long, slender design of these georods makes them easy to deploy and retrieve, allowing researchers to acquire greater line-kilometers of seismic data during field campaigns that are commonly time-constrained.

Solar radiation penetrating the ice is one of the most important factors that determine the functioning of lake ecosystem in late winter. Parameterization of the attenuation of solar radiation in the snow-ice sheet is an essential tool in the study of the light regime of ice-covered lakes. The optical properties of the snow-ice sheet in Vendyurskoe lake, northwestern Russia, are investigated on the basis of long-term field observations (1995–2012). The four-layer approach (snow, white ice, slush and congelation ice) is used to study the attenuation of the downwelling planar irradiance in the snow-ice sheet. The bulk attenuation coefficients for four layers (18.8 m–1 for snow, 6 m−1 for white ice, 3.5 m−1 for slush and 2.1 m−1 for congelation ice) are calculated by the quasi-Newton method. A comparison of observed and calculated values of the irradiance beneath the ice shows that the determined coefficients adequately describe the attenuation of the downwelling irradiance by snow-ice cover.

We identified ikaite crystals (CaCO3·6H2O) and examined their shape and size distribution in first-year Arctic pack ice, overlying snow and slush layers during the spring melt onset north of Svalbard. Additional measurements of total alkalinity (TA) were made for melted snow and sea-ice samples. Ikaite crystals were mainly found in the bottom of the snowpack, in slush and the surface layers of the sea ice where the temperature was generally lower and salinity higher than in the ice below. Image analysis showed that ikaite crystals were characterized by a roughly elliptical shape and a maximum caliper diameter of 201.0±115.9 μm (n = 918). Since the ice-melting season had already started, ikaite crystals may already have begun to dissolve, which might explain the lack of a relationship between ikaite crystal size and sea-ice parameters (temperature, salinity, and thickness of snow and ice). Comparisons of salinity and TA profiles for melted ice samples suggest that the precipitation/dissolution of ikaite crystals occurred at the top of the sea ice and the bottom of the snowpack during ice formation/melting processes.

The Bering Glacier–Bagley Icefield system in Alaska is currently surging (2011). Large-scale elevation changes and small-scale elevation-change characteristics are investigated to understand surge progression, especially mass transport from the pre-surge reservoir area to the receiving area and propagation of the kinematic surge wave as manifested in heavy crevassing characteristic of rapid, brittle deformation. This analysis is based on airborne laser altimeter data collected over Bering Glacier in September 2011. Results include the following: (1) Maximal crevasse depth is 60 m, reached in a rift that separates two deformation domains, indicative of two different flow regimes. Otherwise surge crevasse depth reaches 20–30 m. (2) Characteristic parameters of structural provinces are derived by application of geostatistical classification. Parameters include significance and spacing of crevasses, surface roughness and crevasse-edge curvature (indicative of crevasse age). A classification based on these parameters serves to objectively discriminate structural provinces, indicative of surge progression down-glacier and up-glacier. (3) Elevation changes from 2011 and 2010 altimetry show 40–70 m surface lowering in the reservoir area in lower central Bering Glacier and 20–40m thickening near the front in Tashalich arm. Combining elevation changes with results of crevasse profilometry and pattern analysis, the rapid progression of the surge can be mathematically–physically reconstructed.

A passive seismology experiment was conducted across the main overdeepening of Storglaciären in the Tarfala valley, northern Sweden, to investigate the spatial and temporal distribution of basal microseismic waveforms in relation to known dynamics of this small polythermal sub-arctic glacier. The high ablation rate made it difficult to keep geophones buried and well coupled to the glacier during the experiment and reduced the number of days of good-quality data collection. The characterization of typical and atypical waveforms showed that the dominant waveforms were from near-surface events such as crevassing. Waveforms resembling basal microseismic signals were very rare, and seldom observed on more than two seismic stations simultaneously. The analysis of waveforms, amplitudes and particle motions suggested a near-field origin for most events. Even though basal sliding is known to occur in the overdeepening, no convincing examples of basal waveforms were detected, suggesting basal microseismic signals are rare or difficult to detect beneath polythermal glaciers like Storglaciären. We discuss the reasons for failing to locate basal signals, consider the origin of common waveforms and make recommendations for setting up passive seismology experiments on glaciers with high ablation rates.

To investigate recent glacier changes in the Himalayan region, we carried out GPS and ground-penetrating radar (GPR) measurements at Yala Glacier, a benchmark glacier in Nepal. Glacier surface elevation and ice thickness were surveyed along a 1.5 km profile from the glacier top to the terminus. Ice flow velocity was measured at five locations by surveying stakes for either 1 year or 4 day periods. Obtained surface elevation and ice velocity were compared with those measured in 1982 and 1996. The mean ice thickness along the radar profile was 36 m in 2009 and the ice has been thinning at rates of-0.69 ±0.25 and -0.75 ± 0.24 m a-1 during the periods 1982-96 and 1996-2009, respectively. The thinning rate increases down-glacier, reaching a magnitude up to -1.8 m a-1 near the terminus from 1996 to 2009. The ice velocity has reduced by >70% from 1982 to 2009 in the lower half of the glacier. By assuming a constant driving stress over the glacier, the total ice volume in 2009 was estimated as 0.061 km3. Our results indicate that Yala Glacier has lost ∼40% of its ice volume over the last 27 years and that the rate of the mass loss has accelerated over the last decade.

Surface elevation data for sea ice in the northwesternty - Weddell Sea, Antarctica, collected by a helicopter-borne laser altimeter during the Winter Weddell Outflow Study 2006, were used to estimate the form drag on pressure ridges and its contribution to the total wind drag, and the air-ice drag coefficient at a reference height of 10 m under neutral stability conditions (C dn(10)). This was achieved by partitioning the total wind drag into two components: form drag on pressure ridges and skin drag over rough sea-ice surfaces. The results reveal that for the compacted ice field, the contribution of form drag on pressure ridges to the total wind drag increases with increasing ridging intensity Ri (where Ri is the ratio of mean ridge height to spacing), while the contribution decreases with increasing roughness length. There is also an increasing trend in the air-ice drag coefficient C dn(10) as ridging intensity R i increases. However, as roughness length increases, C dn(10) increases at lower ridging intensities (Ri < 0.023) but decreases at lower ridging intensities (0.023 < Ri < 0.05). These opposing trends are mainly caused by the dominance of the form drag on pressure ridges and skin drag over rough ice surfaces. Generally, the form drag becomes dominant only when the ridging intensity is sufficiently large, while the skin drag is the dominant component at relatively larger ridging intensities. These results imply that a large value of C dn(10) is caused not only by the form drag on pressure ridges, but also by the skin drag over rough ice surfaces. Additionally, the estimated drag coefficients are consistent with reported measurements in the northwestern Weddell Sea, further demonstrating the feasibility of the drag partition model.

Deriving glacier outlines from satellite data has become increasingly popular in the past decade. In particular when glacier outlines are used as a base for change assessment, it is important to know how accurate they are. Calculating the accuracy correctly is challenging, as appropriate reference data (e.g. from higher-resolution sensors) are seldom available. Moreover, after the required manual correction of the raw outlines (e.g. for debris cover), such a comparison would only reveal the accuracy of the analyst rather than of the algorithm applied. Here we compare outlines for clean and debris-covered glaciers, as derived from single and multiple digitizing by different or the same analysts on very high- (1 m) and medium-resolution (30 m) remote-sensing data, against each other and to glacier outlines derived from automated classification of Landsat Thematic Mapper data. Results show a high variability in the interpretation of debris-covered glacier parts, largely independent of the spatial resolution (area differences were up to 30%), and an overall good agreement for clean ice with sufficient contrast to the surrounding terrain (differences ∼5%). The differences of the automatically derived outlines from a reference value are as small as the standard deviation of the manual digitizations from several analysts. Based on these results, we conclude that automated mapping of clean ice is preferable to manual digitization and recommend using the latter method only for required corrections of incorrectly mapped glacier parts (e.g. debris cover, shadow).