We conducted two-dimensional continuous multi-offset georadar surveys on Bench Glacier, south-central Alaska, USA, to measure the distribution of englacial water. We acquired data with a multichannel 25 MHz radar system using transmitter–receiver offsets ranging from 5 to 150 m. We towed the radar system at 5–10 kmh–1 with a snow machine with transmitter/receiver positions established by geodetic-grade kinematic differentially corrected GPS (nominal 0.5 m trace spacing). For radar velocity analyses, we employed reflection tomography in the pre-stack depth-migrated domain to attain an estimated 2% velocity uncertainty when averaged over three to five wavelengths. We estimated water content from the velocity structure using the complex refractive index method equation and use a three-phase model (ice, water, air) that accounts for compression of air bubbles as a function of depth. Our analysis produced laterally continuous profiles of glacier water content over several kilometers. These profiles show a laterally variable, stratified velocity structure with a low-water-content (~0–0.5%) shallow layer (~20–30 m) underlain by high-water-content (1–2.5%) ice.

Radar profiles of bed echo intensity can survey conditions at the ice–bed interface and test for the presence or absence of water. However, extracting information about basal conditions from bed echo intensities requires an estimate of the attenuation loss through the ice. We used the relationship between bed echo intensities from constant-offset radar data and ice thickness to estimate depth-averaged attenuation rates at several locations on and near Kamb Ice Stream (KIS), West Antarctica. We found values varying from 29 dBkm–1 at Siple Dome to 15 dBkm–1 in the main trunk region of KIS, in agreement with a previous measurement and models. Using these attenuation-rate values, we calculated the relative bed reflectivity throughout our KIS surveys and found that most of the bed in the trunk has high basal reflectivities, similar to those obtained in the location of boreholes that found water at the bed. Areas of lower bed reflectivity are limited to the sticky spot, where a borehole found a dry bed, and along the margins of KIS. These results support previous hypotheses that the basal conditions at locations like the sticky spot on KIS control its stagnation and possible reactivation.

We present the results of several radio-echo sounding surveys carried out on Johnsons and Hurd Glaciers, Livingston Island, Antarctica, between the 1999/2000 and 2004/05 austral summer campaigns, which included both radar profiling and common-midpoint measurements with low (20– 25 MHz)- and high (200 MHz)-frequency radars. The latter have allowed us to estimate the radio-wave velocity in ice and firn and the corresponding water contents in temperate ice, which vary between 0 and 1.6% depending on the zone. Maximum ice thickness is ~200 m, with a mean value of 93.6±2.5 m. Total ice volume is 0.968±0.026 km3, for an area of 10.34±0.03km2. The subglacial relief of Johnsons Glacier is quite smooth, while that of Hurd Glacier shows numerous overdeepenings and peaks. The radar records suggest that Hurd Glacier has a polythermal structure, contrary to the usual assumption that glaciers in Livingston Island are temperate. This is also supported by other dynamical and geomorphological evidence.

Accumulation rates and their spatio-temporal variability are important boundary conditions for ice-flow models. The depths of radar-detected internal layers can be used to infer the spatial variability of accumulation rates. Here we infer accumulation rates from three radar layers (26, 35 and 41 ka old) in the Vostok Subglacial Lake region using two methods: (1) the local-layer approximation (LLA) and (2) a combination of steady-state flowband modeling and formal inverse methods. The LLA assumes that the strain-rate history of a particle traveling through the ice sheet can be approximated by the vertical strain-rate profile at the current position of the particle, which we further assume is uniform. The flowband model, however, can account for upstream strain-rate gradients. We use the LLA to map accumulation rates over a 150 km × 350 km area, and we apply the flowband model along four flowbands. The LLA accumulation-rate map shows higher values in the northwestern corner of our study area and lower values near the downstream shoreline of the lake. These features are also present but less distinct in the flowband accumulation-rate profiles. The LLA-inferred accumulation-rate patterns over the three time periods are all similar, suggesting that the regional pattern did not change significantly between the start of the Holocene and the last ~20 ka of the last Glacial Period. However, the accumulation-rate profiles inferred from the flowband model suggest changes during that period of up to 1 cma–1 or ~50% of the inferred values.

The 2006/07 and 2007/08 US-ITASE traverses from Taylor Dome to South Pole in East Antarctica provided opportunities to survey the subglacial and englacial environments using 3 MHz and 200MHz radar. We present first results of these new ground-based radar data. A prominent basal deformation layer indicates different ice-flow regimes for the northern and southern halves of the Byrd Glacier drainage. Buried dune stratigraphy that appears to be related to the megadunes towards the west occurs at depths of up to 1500 m. At least two new water-filled subglacial lakes were discovered, while two recently drained lakes identified from repeat ICESat surface elevation surveys appear to be devoid of water.

Rutford Ice Stream, West Antarctica, drains a catchment of >45 000 km2 into the Ronne Ice Shelf through a 26 km wide, 2.4 km deep subglacial trough adjacent to the Ellsworth Mountains. Forty-two per cent of its catchment boundary is common with Pine Island Glacier, where rapid change to ice dynamics is currently underway. These changes may eventually affect adjacent catchments such as Rutford. Radar sounding data were acquired over the Rutford ice-drainage basin that show the internal structure. In particular, distinctive reflector groups were identified that mark the boundaries between four different flow elements. The flow-margin reflector groups include curvilinear events that cross-cut isochrones and are therefore likely to be post-depositional. These reflections may arise from crystal orientation fabrics generated by localized strain in a flow margin. One of the sectors of the ice-drainage basin supplies the largest share (38%) of the ice volume flux through the main trunk of Rutford Ice Stream. This sector may be preferentially affected by continuing surface lowering in the Pine Island Glacier catchment.

We present the results of low-frequency (20 MHz) radio-echo sounding (RES) carried out in December 2000 and December 2006 on the main ice divides of Livingston Island, South Shetland Islands (SSI), Antarctica, and Bowles Plateau, Antarctica, respectively, as well as high-frequency (200 MHz) RES on the latter, aimed at determining the ice thickness, internal structure and subglacial relief. Typical ice thickness along the main ice divides is ~150 m, reaching maxima of ~200 m. On Bowles Plateau the ice is much thicker, with an average of 265 m and maxima of ~500 m. The bed below the main ice divides is above sea level, while part of the outlet glaciers from Bowles Plateau lies significantly below sea level, down to –120 m. The strong scattering of the radio waves in the areas under study constitutes further evidence that the ice in the accumulation area of the ice masses of the SSI is temperate. Typical thickness of the firn layer in Bowles Plateau is 20–35 m, similar to that found in King George ice cap. A strong internal reflector within the firn layer, interpreted as a tephra layer from the 1970 eruption at Deception Island, has allowed a rough estimate of the specific mass balances for Bowles Plateau within 0.20–0.40ma–1w.e., as average values for the period 1970–2006.

We present 1 and 100 MHz ground-based radar data from the onset region of Rutford Ice Stream, West Antarctica, which indicate the form and internal structure of isochrones. In the flow-parallel lines, modelled isochrone patterns reproduce the gross pattern of the imaged near-surface layers, assuming steady-state flow velocity from GPS records and the current accumulation rate for the last 200 years. We interpret this as indicating overall stability in flow in the onset region of Rutford Ice Stream throughout this period. However, in the cross-flow lines some local variability in accumulation is seen in areas close to the ice-stream margin where a number of tributaries converge towards the ice-stream onset zone. Episodic surface lowering events are observed followed by rapid fill episodes. The fill events indicate deposition towards the northwest, most likely generated by storm winds, which blow at an oblique angle to ice flow. More problematic is explaining the generation of episodic surface lowering in this area. We speculate this may be due to: changing ice-flow direction in the complex tributary area of the onset zone; a change in basal sediments or sedimentary landforms; a change in basal melt rates or water supply; or episodic lake drainage events in the fjord systems of the Ellsworth Subglacial Highlands. The study highlights the difficulty of assessing flow stability in the complex onset regions of West Antarctic ice streams.

We present the thermal distribution in the confluence area of Gorner- and Grenzgletscher, Valais, Switzerland. The area was mapped by ice-penetrating radar at 1–5 and 40 MHz. The higher-frequency data reveal a thick surface layer of low backscatter in the center of the Grenzgletscher branch. Based on datasets of borehole-temperature measurements and flow velocity, we interpret this as a thick layer of cold ice, advected from the accumulation region of Grenzgletscher. Along seven profiles the base of the low-backscatter zone can be found at a maximum depths between approximately 100 and 200 m. Laterally, the layer extends some 400 m, ∼1/3 of the width of the Grenzgletscher branch. The lower boundary of the low-backscatter zone is systematically higher than the cold–temperate transition surface found in the boreholes. This discrepancy is attributed to the direct sensitivity of radar backscatter to liquid-water inclusions, rather than to the temperature distributions as observed in boreholes. We present the current state of the cold layer and discuss its influence on other glacier characteristics.

Global warming is causing an apparent rapid retreat of many glaciers worldwide. In addition to mass-balance investigation, the determination and monitoring of total glacial ice volume and ice-thickness distribution are important parameters for understanding the interactions between climate and the complex glacier system. Because of spatially irregular and sparse datasets, scaling of volume and ice-thickness distribution is often a challenging problem. This study focuses on two small (<2 km2) temperate glaciers in the Hohe Tauern (Eastern Alps) region of central Austria. The period 2003–04 saw the first use of ground-penetrating radar (GPR) to determine the total ice volume and ice-thickness distribution of the two glaciers. A centre frequency of 20 MHz was used in point measuring mode. Despite variable data quality, bedrock reflections up to depths of >100m were identified in the data. The acquired GPR data are irregularly distributed and the spatial density is too low to calculate reasonable bedrock topography with standard interpolation approaches. Thus one main focus of this study was to develop an appropriate interpolation technique. Eventually, kriging technique and a glacial mechanically based interpolation parameter were used. Mean calculated ice thicknesses for the two investigated glaciers are 40–50 m, with a maximum of 150–165 m. No direct validation data are available, so different considerations support the computed bedrock topography.

During three summer field seasons (2003, 2005 and 2006) we carried out radio-echo sounding measurements with a 5MHz (central frequency) ice-penetrating radar on McCall Glacier, Arctic Alaska, USA, along the central flowline and 17 cross-profiles. Two-way travel time was, after migration, converted to ice thickness, which, in combination with a recent digital elevation model of the surface of the glaciated area, resulted in a detailed map of the bed topography. This reveals a complex basal topography in the confluence area of the different glacial cirques. The pattern of subglacial water flow following the hydraulic potential gradient was calculated for the whole glacier area and shows a confluence of subglacial water downstream from the confluence of the glacier cirques. From the ice-thickness map the total ice volume was estimated as slightly less than 0.5 km3. Bed reflection power (BRP) was determined for the glacier after correction for ice-thickness dependence. Results reveal a clear relationship between the BRP pattern and basal sliding anomalies along the central flowline.

We have used ground-penetrating radar (GPR) to observe englacial structural control upon the development of an esker formed during a high-magnitude outburst flood (jökulhlaup). The surge-type Skeiðarárjökull, an outlet glacier of the Vatnajökull ice cap, Iceland, is a frequent source of jökulhlaups. The rising-stage waters of the November 1996 jökulhlaup travelled through a dense network of interconnected fractures that perforated the margin of the glacier. Subsequent discharge focused upon a small number of conduit outlets. Recent ice-marginal retreat has exposed a large englacial esker associated with one of these outlets. We investigated structural controls on esker genesis in April 2006, by collecting >2.5km of GPR profiles on the glacier surface up-glacier of where the esker ridge has been exposed by meltout. In lines closest to the exposed esker ridge, we interpret areas of englacial horizons up to ~30m wide and ~10–15m high as an up-glacier continuation of the esker sediments. High-amplitude, dipping horizons define the base of esker materials across many lines. Similar dipping surfaces deeper in the profiles suggest that: (1) the dipping surfaces beneath the esker are englacial tephera bands; (2) floodwaters were initially discharged along structurally controlled englacial surfaces (tephra bands); (3) the rapid increase in discharge resulted in hydrofracturing; (4) establishment of preferential flow paths resulted in conduit development along the tephra bands due to localized excavation of surrounding glacier ice; and (5) sedimentation took place within the new accommodation space to form the englacial structure melting out to produce the esker.

Ice-thickness measurements are needed to calculate fluxes through fast-flowing outlet glaciers in Greenland, Alaska, Patagonia and Antarctica. However, relatively high attenuation of radio waves by dielectric absorption and volume scattering from englacial water hampers detection of the bed through warm deep ice. In the past we have had success measuring ice thickness of temperate glaciers using a ground-based monopulse radar system operating at low frequencies (2 MHz). Here we adapt the same system to operate from an airplane. Test flights over Bering Glacier, Alaska, USA, detected the bed through ice up to 1250m thick. Flights across the Seward–Malaspina Glacier system, Alaska, resolved the ice thickness of Malaspina Glacier, but strong hyperbolic-shaped returns obscured the bed echo through the Seward throat. It is likely that this clutter in the signal was caused by off-nadir returns from chaotic surface crevasses that are ubiquitous in the throat region.

Visible and infrared satellite images reveal numerous lineations on the Siple Coast region of West Antarctica. We used 5 MHz ice-penetrating radar to probe the interior and the bed of the ice sheet beneath a lineation at the boundary between Engelhardt Ice Ridge and flat-ice terrain to the south of the Kamb Ice Stream (KIS) outlet. Results show curved reflectors that emerge from the bed beneath 600 m thick ice. The tops of the reflectors extend about 100m into the ice above the bed, where they become almost horizontal. Apparent reflectivity of the horizontal section is about 20 dB less than that of the bed. We conclude that the likely cause of such strong reflection is sea water that was accreted into basal crevasses when the flat-ice terrain was floating. Internal layers are warped downward just downslope from the basal reflectors. It is thought that the downwarping was caused by localized basal melting in the past. The spatial pattern of downwarping suggests that localized basal melting was stronger on the north side than on the south side of KIS; apparently ice/ocean interactions on the two sides of KIS were different.

In 2005 an airborne survey was carried out from a Twin Otter aircraft at Pâkitsup Akuliarusersua (Paakitsoq) near Ilulissat in West Greenland. The survey aimed to measure ice thickness with a 60 MHz coherent radar and surface elevation with a scanning laser altimeter. Positioning information came from multiple on-board differential GPS units and an inertial navigation system. The region surveyed covers >80km along the ice margin and has a total area of ~2700km2 with varying density of measurements: the between-track distance was ~1 km near the margin, increasing to ~3km away from the margin. Regional high-resolution maps of basal topography under the Greenland ice sheet are useful for resolving important glaciological and hydrological questions and for enhancing related process studies, such as the influence of basal meltwater on ice dynamics. The ice-sheet margin in this region is also currently under consideration for hydropower development and has a long and continuing history of glaciological investigations, lately with emphasis on the connection between surface meltwater formation and surface velocity of the ice sheet. Here we present a new regional map of the surface and basal topography of the ice-sheet margin and discuss some of the implications for reported observations at Swiss Camp.

We used internal ice layers from a radio-echo sounding profile between the Kohnen and Dome Fuji deep drilling sites to infer the spatio-temporal pattern of accumulation rate in this sector of Dronning Maud Land, East Antarctica. Continuous internal reflection horizons can be traced to about half of the ice thickness and have a maximum age of approximately 72.7 ka BP. To infer palaeo-accumulation rates from the dated layers, we derived the thinning functions from a flow calculation with a high-resolution higher-order model of Dronning Maud Land embedded into a three-dimensional thermomechanical model of the Antarctic ice sheet. The method takes into account complex ice-flow dynamics and advection effects that cannot be dealt with using traditional local approaches. We selected seven time intervals over which we determine the average accumulation rate and average surface temperature at the place and time of origin of the layer particles. Our results show lower accumulation rates along eastern parts of the profile for the late Holocene (0–5 ka BP) than are shown by existing maps, which had no surface control points. During the last glacial period we find a substantially lower accumulation rate than predicted by the usual approach linking palaeo-accumulation rates to the condensation temperature above the surface inversion layer. These findings were used to fine-tune the relation between accumulation rate and temperature.

Englacial horizons deeper than 100m are absent within 100MHz ground-penetrating radar (GPR) surface profiles we recorded on Clark and Commonwealth Glaciers in the Antarctic Dry Valleys region. Both glaciers show continuous bottom horizons to 280 m, with bottom signal-to-noise ratios near 30 dB. Density horizons should fade below 50 m depth because impermeable ice occurred by 36 m. Folding within Commonwealth Glacier could preclude radar strata beneath about 80 m depth, but there is no significant folding within Clark Glacier. Strong sulfate concentrations and contrasts exist in our shallow ice core. However, it appears that high background concentration levels, and possible decreased concentration contrasts with depth placed their corresponding reflection coefficients at the limit of, or below, our system sensitivity by about 77m depth. Further verification of this conclusion awaits processing of our deep-core chemistry profiles.

A finite-difference model for the calculation of radar layer geometries in large ice masses is presented. Balance velocities are used as coefficients in the age equation and in the heat equation. Solution of the heat equation allows prediction of sliding areas and computation of basal melt rates. Vertical distributions of velocity are parameterized using shape functions. These can be set uniformly, or allowed to vary in space according to the distribution of sliding. The vertical coordinate can either be uniformly distributed within the thickness of the ice, or be uniformly distributed within the flux. The finite-difference scheme results in a large set of linear equations. These are solved using a nested factorization preconditioned conjugate gradient scheme. The convergence properties of some other iteration solution schemes are studied. The output is computations of age and temperature assuming steady state, in large ice masses at high resolution. Age calculations are used to generate isochrones which show the best fit to observed layers. Comparisons with analytical solutions are made, and the influence of the order of the finite-difference approximation and the choice of vertical coordinate on solution accuracy is considered.

This paper presents an overview of internal layering across Pine Island Glacier, West Antarctica, as measured from airborne-radar data acquired during a survey conducted by the British Antarctic Survey and the University of Texas in the 2004/05 season. Internal layering is classified according to type (continuous/discontinuous/missing) and the results compared with InSAR velocities. Several areas exhibit disruption of internal layers that is most likely caused by large basal shear stresses. Signs of changes in flow were identified in a few inter-tributary areas, but overall the layering classification and distribution of layers indicate that only minor changes in ice-flow regime have taken place. This is supported by bed-topography data that show the main trunk of the glacier, as well as some of the tributaries, are topographically controlled and located in deep basins.