This study aims to demonstrate that deep ice cores can be synchronized using internal horizons in the ice between the drill sites revealed by airborne radio-echo sounding (RES) over a distance of >1000km, despite significant variations in glaciological parameters, such as accumulation rate between the sites. In 2002/03 a profile between the Kohnen station and Dome Fuji deep ice-core drill sites, Antarctica, was completed using airborne RES. The survey reveals several continuous internal horizons in the RES section over a length of 1217 km. The layers allow direct comparison of the deep ice cores drilled at the two stations. In particular, the counterpart of a visible layer observed in the Kohnen station (EDML) ice core at 1054 m depth has been identified in the Dome Fuji ice core at 575 m depth using internal RES horizons. Thus the two ice cores can be synchronized, i.e. the ice at 1560 m depth (at the bottom of the 2003 EDML drilling) is ∼49ka old according to the Dome Fuji age/depth scale, using the traced internal layers presented in this study.

We have conducted a series of experiments designed to investigate elastic and electromagnetic (EM) velocity anisotropy associated with a preferentially aligned fracture system on a temperate valley glacier in south-central Alaska, USA. Measurements include a three-dimensional compressional wave (P-wave) seismic reflection survey conducted over a 300 m x 300 m survey patch, with uniform source grid and static checkerboard receiver pattern. Additionally, we acquired a multi-azimuth, multi-offset, polarimetric ground-penetrating radar (GPR) reflection experiment in a wagon-wheel geometry with 94° of azimuthal coverage. Results show azimuthal variation in the P-wave normal-moveout velocity of >3% (3765 and 3630 m s–1 in the fast and slow directions respectively) and difference of nearly 5% between the fast (0.164 m ns–1) and slow (0.156 m ns–1) EM velocities. Fracture orientations estimated from the GPR and seismic velocity data are consistent and indicate a preferred fracture orientation that is 30-45° oblique to glacier flow; these measurements agree with borehole observations. Anisotropic analysis of the polarimetric data gives a single volumetric water content estimate of 0.73 ±0.11%. We conclude that meaningful estimates of physical properties in glaciers based on EM or seismic velocity measurements require collecting data such that the presence of anisotropy can be evaluated and an anisotropic analysis employed when necessary.

The ongoing retreat of mountain glaciers necessitates the development of future scenarios of glacier runoff. These scenarios are not only governed by future climate scenarios influencing glacier mass balance but also by the glacier volumes, which are subject to melt. Ground-penetrating radar (GPR) is a valuable tool for measuring the thickness of mountain glaciers, although ground-based measurements are labour-intensive, so not all glaciers can be surveyed. This study presents the results of GPR surveys on 64 Alpine glaciers, carried out between 1995 and 2010. The glacier areas range from 0.001 to 18.4 km2, and their ice thickness was surveyed with an average density of 36 points km-2. The point measurements were extrapolated manually to derive volume maps. The mean ice thickness varies between 10 and 92 m; the maximum ice thickness is about three times the mean thickness. According to the glacier state recorded in the second glacier inventory, the 64 glaciers cover an area of 223.3 ± 3.6 km2, with a mean thickness of 50 ± 3 m and a glacier volume of 11.9 ± 1.1 km3. The mean maximum ice thickness is 119 ± 5m.

Explosive seismic reflection data from Halvfarryggen, a 910 m thick local ice dome of the Antarctic ice sheet, show numerous laterally continuous reflections within the ice between 300 and 870 m depth. We compare the quality of data obtained with explosive sources with that obtained using a vibroseis source for detecting englacial reflections with a snowstreamer, and investigate the origin of englacial reflections. We find vibroseis in combination with a snowstreamer is ten times more productive than explosive seismics. However, englacial reflections are more clearly visible with explosives, which have a broader bandwidth signature, than the vibroseis, which is band-limited at the high-frequency end to 100 Hz. Only the strongest and deepest englacial reflection is detected with vibroseis. We interpret the majority of englacial reflections to originate from changes in the crystal orientation fabric in closely spaced layers, less than the vibro-seismic tuning thickness of 13.5 m. Phase analysis of the lowermost englacial reflector, 40 m above the bed, indicates a sharp increase in seismic wave speed. We interpret this reflector as a transition to a vertical single-maximum fabric. Our findings support current results from anisotropic ice-flow models, that crystal fabric is highly anisotropic at ice domes, both laterally and vertically.

Two seismic surveys were carried out on the high-altitude glacier saddle, Colle Gnifetti, Monte Rosa, Italy/Switzerland. Explosive and vibroseismic sources were tested to explore the best way to generate seismic waves to deduce shallow and intermediate properties (<100 m) of firn and ice. The explosive source (SISSY) excites strong surface and diving waves, degrading data quality for processing; no englacial reflections besides the noisy bed reflector are visible. However, the strong diving waves are analyzed to derive the density distribution of the firn pack, yielding results similar to a nearby ice core. The vibrator source (ElViS), used in both P- and SH-wave modes, produces detectable laterally coherent reflections within the firn and ice column. We compare these with ice-core and radar data. The SH-wave data are particularly useful in providing detailed, high-resolution information on firn and ice stratigraphy. Our analyses demonstrate the potential of seismic methods to determine physical properties of firn and ice, particularly density and potentially also crystal-orientation fabric.

As part of ongoing work to obtain a reliable estimate of the total ice volume of Svalbard glaciers and their potential contribution to sea-level rise, we present here volume calculations, with detailed error estimates, for ten glaciers on western Nordenskiöld Land, central Spitsbergen, Svalbard. The volume estimates are based upon a dense net of GPR-retrieved ice thickness data collected over several field campaigns spanning the period 1999-2012. The total area and volume of the ensemble are 116.06 ± 4.53 km2 and 10.439 ±0.373 km3, respectively, while the individual areas, volumes and average ice thickness lie within 2.6-50.4 km2, 0.08-5.54 km3 and 29-108 m, respectively. Volume/area scaling relationships overestimate the total volume of these glaciers by up to 35% with respect to our calculation. On the basis of the pattern of scattering in the radargrams, we also analyse the hydrothermal structure of these glaciers. Nine of the ten are polythermal, while only one is entirely cold.

P-band interferometric synthetic aperture radar (InSAR) data at 5 m resolution from Kahiltna Glacier, the largest glacier in the Alaska Range, Alaska, USA, show pronounced spatial variation in penetration depth, δP. We obtained δP by differencing X- and P-band digital elevation models. δP varied significantly over the glacier, but it was possible to distinguish representative zones. In the accumulation area, δP decreased with decreasing elevation from 18 ± 3 m in the percolation zone to 10 ± 4 m in the wet snow zone. In the central portion of the ablation area, a location free of debris and crevasses, we identified a zone of very high δP (34 ± 4 m) which decreased at lower elevations (23 ± 3 m in bare ice and 5-10m in debris-covered ice). We observe that the spatial configuration of δP is consistent with the expected thermal regime of each zone: δP is high in areas where cold firn/ice likely occurs (i.e. percolation zone and upper ablation area) and low in areas where temperate surface firn/ice likely exists (wet snow zone and lower ablation area). We suggest that the very high δP observed in the upper ablation area is due to the presence of a cold surface layer.