Chemical constituents trapped within glacial ice provide a unique record of climate, as well as repositories for biological material such as pollen grains, fungal spores, viruses, bacteria and dissolved organic carbon. Past research suggests that the veins of polycrystalline ice may provide a liquid microenvironment for active microbial metabolism fueled by concentrated impurities in the veins. Despite these claims, no direct measurements of impurity concentration in ice veins have been made. Using micro-Raman spectroscopy, we show that sulfate and nitrate concentrations in the veins of glacial ice from Greenland (Greenland Ice Sheet Project 2) and Antarctic (Newall Glacier and a Dominion Range glacier) core samples were 104 and 105 times greater than the concentrations measured in melted (bulk) core water. Methanesulfonate was not found in the veins, consistent with its presence as particulate matter within the ice. The measured vein concentration of molecular anions implies a highly acidic (pH < 3) vein environment with high ionic strength (mM-M). We estimate that the vein volume provides 16.7 and 576 km3 of habitable space within the Greenland and Antarctic ice sheets, respectively, which could support the metabolism of organisms that are capable of growing in cold, high ionic strength solutions with low pH.

Interferometric synthetic aperture radar (InSAR) observations of ice-shelf flow contain ocean-tide and atmospheric-pressure signals. A model-based correction can be applied, but this method is limited by its dependency upon model accuracy, which in remote regions can be uncertain. Here we describe a method to determine two-dimensional ice-shelf flow vectors independently of model predictions of tide and atmospheric pressure, by stacking conventional and multiple aperture InSAR (MAI) observations of the Dotson Ice Shelf, West Antarctica. In this way we synthesize a longer observation period, which enhances long-period (flow) displacement signals, relative to rapidly varying (tide and atmospheric pressure) signals and noise. We estimate the error associated with each component of the velocity field to be ~22 ma-1, which could be further reduced if more images were available to stack. With the upcoming launch of several satellite missions, offering the prospect of regular short-repeat SAR acquisitions, this study demonstrates that stacking can improve estimates of ice-shelf flow velocity.

Given the cross-sectional geometry of a valley glacier, effects of lateral drag can be parameterized in flowline models through the introduction of Nye shape factors. Lateral drag also arises due to lateral variability in bed topography and basal flow, which induce horizontal shear stress and differential ice motion. For glaciers with various geometric and basal conditions, we compare three-dimensional Stokes solutions to flowline model solutions to examine both sources of lateral drag. We calculate associated correction factors that help flowline models to capture the effects of lateraldrag. Such parameterizations provide improved simulations of the dynamics of narrow, channelized, fast-flowing glacial systems. We present an example application for Athabasca Glacier in the Canadian Rocky Mountains.

We propose an infinite-dimensional adjoint-based inexact Gauss-Newton method for the solution of inverse problems governed by Stokes models of ice sheet flow with nonlinear rheology and sliding law. The method is applied to infer the basal sliding coefficient and the rheological exponent parameter fields from surface velocities. The inverse problem is formulated as a nonlinear least-squares optimization problem whose cost functional is the misfit between surface velocity observations and model predictions. A Tikhonov regularization term is added to the cost functional to render the problem well-posed and account for observational error. Our findings show that the inexact Newton method is significantly more efficient than the nonlinear conjugate gradient method and that the number of Stokes solutions required to solve the inverse problem is insensitive to the number of inversion parameters. The results also show that the reconstructions of the basal sliding coefficient converge to the exact sliding coefficient as the observation error (here, the noise added to synthetic observations) decreases, and that a nonlinear rheology makes the reconstruction of the basal sliding coefficient more difficult. For the inversion of the rheology exponent field, we find that horizontally constant or smoothly varying parameter fields can be reconstructed satisfactorily from noisy observations.

We present measurements with inclinometer/magnetometer probes inserted in glacier ice, and discuss the underlying mathematical theory to extract velocity-gradient components from these data. Assumptions concerning the ice flow field must be made to reduce the number of unknowns and to close the system of equations obtained from the theory. With the incompressibility assumption and the first-order plane-strain approximation, the evolution equation for the tilt of the main sensor axis can be solved exactly and the obtained function can be fitted to measured data to obtain optimal shear and normal strain components. Daily variations superposed on the tilt evolution reflect variations in sliding and perhaps partial elastic recovery of deformations.

During the Japanese-Swedish Antarctic traverse expedition of 2007/08, we measured the surface snow density at 46 locations along the 2800 km long route from Syowa station to Wasa station in East Antarctica. The mean snow density for the upper 1 (or 0.5) m layer varied from 333 to 439 kg m-3 over a region spanning an elevation range of 365-3800 ma.s.l. The density variations were associated with the elevation of the sampling sites; the density decreased as the elevation increased, moving from the coastal region inland. However, the density was relatively insensitive to the change in elevation along the ridge on the Antarctic plateau between Dome F and Kohnen stations. Because surface wind is weak in this region, irrespective of elevation, the wind speed was suggested to play a key role in the near-surface densification. The results of multiple regression performed on the density using meteorological variables were significantly improved by the inclusion of wind speed as a predictor. The regression analysis yielded a linear dependence between the density and the wind speed, with a coefficient of 13.5 kg m-3 (m s-1)-1. This relationship is nearly three times stronger than a value previously computed from a dataset available in Antarctica. Our data indicate that the wind speed is more important to estimates of the surface snow density in Antarctica than has been previously assumed.

Glacier surface mass balance evaluated over the actual glacier geometry depends not only on climatic variations, but also on the dynamic adjustment of glacier geometry. Therefore, it has been proposed that reference-surface balances calculated over a constant glacier hypsometry are better suited for climatic interpretation. Here we present a comparison of 82 year modelled time series (1926-2008) of conventional and reference-surface balance for 36 Swiss glaciers. Over this time period the investigated glaciers have lost 22% of their area, and ice surface elevation close to the current glacier terminus has decreased by 78 m on average. Conventional balance in the last decade, at −0.91 mw.e.a-1, is 0.14 m w.e. a-1 less negative than the reference-surface balance. About half of the negative (stabilizing) feedback on mass balance due to glacier terminus retreat is compensated by more negative mass balances due to surface lowering. Short-term climatic variability is clearly reflected in the conventional mass-balance series; however, the magnitude of the long-term negative trend is underestimated compared to that found in the reference-surface balance series. Both conventional and reference-surface specific balances show large spatial variability among the 36 glaciers.

A laboratory device was built to measure the forces that ice exerts on a 0.05 m diameter rigid plastic sphere in two different configurations: in contact with a flat bed or isolated from the bed. Measurements indicated that bed-normal contact forces were 1.8 times larger than drag forces due to creeping flow past a slippery sphere isolated from the bed. Measurements of forces as a function of the bed-normal ice velocity, estimations of the ice viscosity parameter and observations of markers in the ice indicate ice is Newtonian with a viscosity of ~1.3x 1011 Pas. Newtonian behavior is expected due to small and transient stresses. A model of regelation indicates that it had a negligible (<5%) influence on forces. Water pressure in the cavity beneath the sphere in contact with the bed had a likewise negligible influence on contact forces. When no cavity is present, drag forces can be correctly estimated using Stokes’s law (Newtonian viscosity) for a slippery sphere. The same law with a bed-enhancement factor of 1.8 is appropriate for estimating bed-normal contact forces. These results reinforce previous laboratory measurements and theories but provide no support for explanations of high debris/bed friction or rates of abrasion that depend on high contact forces.

The viability of employing sediment plumes emanating from outlets along the western margin of the Greenland ice sheet as indicators of runoff is assessed. An automated sediment plume quantification system based on daily 250 m Moderate Resolution Imaging Spectroradiometer (MODIS) band 1 reflectance imagery is developed. Coherent plumes are identified using spectral thresholds and polygon tracing. Validation employs imagery quality-control procedures and manual verification of plume areas. Outlets at land-terminating margins with wide and straight fjord geometries deliver the most accurate and consistent results. Plume area observations are also possible at marine-terminating margins with relatively static fronts and low proximal sea-ice concentrations. Variability in plume area is examined with reference to Special Satellite Microwave Imager (SSM/I)-derived daily melt extent at the hydrologic catchment scale. At annual timescales, plume areas tend to co-vary with surface melt extent, indicating that more mass is lost by runoff during years of extensive melting. Some synchronicity in plume areas from different catchments is apparent. At seasonal and daily timescales, plumes from individual outlets primarily relate to catchment-specific melting.

Ice-thickness measurements critical for flood and mudflow hazard studies are very sparse on Cascade Range (North America) volcanoes. Helicopter electromagnetic (HEM) data collected to detect hydrothermal alteration are used to determine ice thickness over portions of Mount Baker and Mount Adams volcanoes. A laterally continuous inversion method provides good estimates of ice <100 m thick over water-saturated and altered regions where the resistivity of the basement is <200 Ωm. For areas with ice overlying fresh, resistive rocks with small resistivity contrasts between ice and rock, ice thickness is not well resolved. The ice thicknesses derived from HEM data are consistent with the previous drillhole data from Mount Adams and radar data from both volcanoes, with mean thicknesses of 57 m for Mount Adams and 68 m for Mount Baker. The thickest ice on Mount Baker rests on the gentle lower slopes whereas the thickest ice at Mount Adams lies on the flat summit. Ice volume calculations suggest that Mount Baker contains ∽710 × 106 m3 of ice in the HEM survey area, with a crude estimate of ∽1800 × 106m3 for the entire volcano. Ice volume on Mount Adams is 65 × 106m3 in parts of the HEM survey area and ∽200 × 106m3 overall.

We used ground-penetrating radar (GPR), GPS and glaciochemistry to evaluate melt regimes and ice depths, important variables for mass-balance and ice-volume studies, of Upper Yentna Glacier, Upper Kahiltna Glacier and the Mount Hunter ice divide, Alaska. We show the wet, percolation and dry snow zones located below ~2700ma.s.l., at ~2700 to 3900ma.s.l. and above 3900ma.s.l., respectively. We successfully imaged glacier ice depths upwards of 480 m using 40-100 MHz GPR frequencies. This depth is nearly double previous depth measurements reached using mid-frequency GPR systems on temperate glaciers. Few Holocene-length climate records are available in Alaska, hence we also assess stratigraphy and flow dynamics at each study site as a potential ice-core location. Ice layers in shallow firn cores and attenuated glaciochemical signals or lacking strata in GPR profiles collected on Upper Yentna Glacier suggest that regions below 2800ma.s.l. are inappropriate for paleoclimate studies because of chemical diffusion, through melt. Flow complexities on Kahiltna Glacier preclude ice-core climate studies. Minimal signs of melt or deformation, and depth-age model estimates suggesting ~4815 years of ice on the Mount Hunter ice divide (3912ma.s.l.) make it a suitable Holocene-age ice-core location.

A possible cause for accelerated thinning and break-up of floating marine ice shelves is warming of the water in the cavity below the ice shelf. Accurate bathymetry beneath large ice shelves is crucial for developing models of the ocean circulation in the sub-ice cavities. A grid of free-air gravity data over the floating Larsen C ice shelf collected during the IceBridge 2009 Antarctic campaign was utilized to develop the first bathymetry model of the underlying continental shelf. Independent control on the continental shelf geologic structures from marine surveys was used to constrain the inversion. Depths on the continental shelf beneath the ice shelf estimated from the inversion generally range from about 350 to 650 m, but vary from <300 to >1000 m. Localized overdeepenings, 20-30 km long and 900-1000 m deep, are located in inlets just seaward of the grounding line. Submarine valleys extending seaward from the overdeepenings coalesce into two broad troughs that extend to the seaward limit of the ice shelf and appear to extend to the edge of the continental shelf. The troughs are generally at a depth of 550-700 m although the southernmost mapped trough deepens to over 1000 m near the edge of the ice shelf just south of 68° S. The combination of the newly determined bathymetry with published ice-draft determinations based on laser altimetry and radar data defines the geometry of the water-filled cavity. These newly imaged troughs provide a conduit for water to traverse the continental shelf and interact with the overlying Larsen C ice shelf and the grounding lines of the outlet glaciers.

This study calculates area, volume and elevation changes of two glaciers on James Ross Island, Antarctica, during the period 1979-2006. Davies Dome is a small ice cap. Whisky Glacier is a valley glacier. Ground-penetrating radar surveys indicate ice thickness, which was used for calculations of the bed topography and volume of both glaciers. Maximum measured ice thicknesses of Davies Dome and Whisky Glacier are 83 ± 2 and 157 ± 2 m, respectively. Between 1979 and 2006, the area of the ice cap decreased from 6.23 ± 0.05 km2 to 4.94 ± 0.01 km2 (-20.7%), while the area of the valley glacier reduced from 2.69 ± 0.02 km2 to 2.40 ± 0.01 km2 (-10.6%). Over the same period the volume of the ice cap and valley glacier reduced from 0.23 ± 0.03 km3 to 0.16 ± 0.02 km3 (-30.4%) and from 0.27 ± 0.02 km3 to 0.24 ± 0.01 km3 (-10.6%), respectively. The mean surface elevation decreased by 8.5±2.8 and 10.1 ±2.8m. The average areal (~0.048-0.011 km2a-1) and volumetric (~0.003−0.001 km3 a-1) changes are higher than the majority of other estimates from Antarctic Peninsula glaciers.

Fluxes of atmospheric CO2 have been reported over sea ice during winter and spring. These fluxes are partly driven by the gradient of the CO2 concentration between sea ice and the atmosphere. We present a new non-destructive method to measure the pCO2 of bulk sea ice at its in situ temperature. This method is based on an equilibration procedure between sea ice and a standard gas of known CO2 concentration. The concentration is measured by gas chromatography with a precision of 5%. Tests were performed on artificial standard sea ice and confirmed the reproducibility of the technique in the range of precision of the gas chromatograph. To test the accuracy of this method, the first profiles of pCO2 measured in bulk sea ice are reported and compared with direct in situ measurements of brine pCO2 over depth-integrated intervals.

Laboratory creep deformation experiments have been conducted on initially isotropic laboratory-made samples of polycrystalline ice. Steady-state tertiary creep rates, , were determined at strains exceeding 10% in either uniaxial-compression or simple-shear experiments. Isotropic minimum strain rates, , determined at ˜1 % strain, provide a reference for comparing the relative magnitude of tertiary creep rates in shear and compression through the use of strain-rate enhancement factors, E, defined as the ratio of corresponding tertiary and isotropic minimum creep rates, i.e. . The magnitude of strain-rate enhancement in simple shear was found to exceed that in uniaxial compression by a constant factor of 2.3. Results of experiments conducted at octahedral shear stresses of to = 0.040.80 MPa indicate a creep power-law stress exponent of n = 3 for isotropic minimum creep rates and n = 3.5 for tertiary creep rates. The difference in stress exponents for minimum and tertiary creep regimes can be interpreted as a t 0 stress-dependent level of strain-rate enhancement, i.e. .The implications of these results for deformation in complex multicomponent stress configurations and at stresses below those used in the current experiments are discussed.

The evolution of ice dynamics through surges of four tidewater-terminating glaciers in northwest Svalbard is investigated by remote sensing. A 20 year time series of glacier surface flow speeds and frontal positions is presented covering the recent surges of Monacobreen, Comfortlessbreen, Blomstrandbreen and Fjortende Julibreen. Surface flow speeds were derived using feature tracking between pairs of ERS SAR and ALOS PALSAR images, while frontal positions were taken from the same imagery, as well as more frequent but lower-spatial-resolution Envisat Wide Swath Mode images. During all four surges, increased ice flow caused the tidewater margin to advance while the calving flux was initially reduced to near zero due to compressive stresses limiting crevasse propagation. As ice speed decreased, the terminus continued to advance, until the glacier’s speed had returned to its pre-surge flow rate. Only at this time did the terminus start to retreat and peak iceberg calving flux was established. We conclude that terminus advance closely tracks glacier speed-up, that there is little mass loss through calving during the most active phase of the surge, and that seasonal cycles of terminus positions diminish during the active surge phase.

We present a method of measuring ice ablation using an absolute pressure transducer as part of an automatic weather station (AWS) system, which we have installed in 17 locations on the Greenland ice sheet. The pressure transducer assembly is drilled into the ice, enclosed in a hose filled with antifreeze liquid. The pressure signal registered by the transducer is that of the vertical column of liquid over the sensor, which can be translated to depth, and ice ablation rate, knowing the density of the liquid. Measuring at sub-daily timescales, this assembly is well suited to monitoring ice ablation in remote regions, with clear advantages over other, well-established methods. The pressure transducer system has the potential to monitor ice ablation for several years without re-drilling, and the system is suitable for high-ablation areas (>5ma-1). A routine to transform raw measurements into ablation values is presented, including a physically based method to remove air-pressure variability from the signal. The pressure transducer time series is compared to that recorded by a sonic ranger for the climatically hostile setting on the Greenland ice sheet.

Radio-echo soundings were collected on Black Rapids Glacier, Alaska, USA, from mid-May to mid-July 1993 to investigate spring speed-up and summer slowdown including high-speed events associated with three lake drainages. Temporal changes in echo power from all depths were highly correlated, indicating a strong effect from varying amounts of near-surface water. Evaluation of bed reflectivity was corrected for this effect based on the time variation of spatially stable patterns of internal scattering identified using principal component analysis. Hourly time series collected at two fixed locations over the deepest part of two valley cross sections showed no detectable change in bed reflection power (<5%) or phase (<0.05 rad). Reoccupation of fixed locations toward the margins at several-day intervals revealed changes in bed power reflectivity up to 50%, but with no definable relation to lake drainages. Theoretical analyses indicate that changes in reflectivity of <5% from a rock bed constrain basal water thickness changes to centimeter scale or less. Conductive basal till degrades the constraint to decimeter scale or more. Changes in bed reflectivity of 50% indicate probable absence of thick conductive till at such locations, and that the changes were caused by centimeter to decimeter changes in equivalent water thickness.

Dye-trace breakthrough curves (BTCs) that increase in velocity and decrease in dispersivity through a melt season have been interpreted as indicating a switch from a distributed to a conduit subglacial drainage system, but this interpretation has not been validated in glaciers where the drainage system configuration was independently known. To test if processes other than a change in the configuration of the subglacial drainage system could produce similar BTCs, we measured BTCs from a persistent, mapped subglacial conduit beneath Rieperbreen, Svalbard, which lacks a distributed system because it is frozen to its bed. This conduit produced slow and highly dispersed BTCs early in the melt season when meltwater delivery rates were low, and fast and sharply peaked BTCs after the snowpack had retreated past the injection moulin. At Rieperbreen, the seasonal evolution of BTCs was controlled by decreases in conduit roughness as increased rates of meltwater delivery increased the relative submergence depths of rocks on the conduit floor. Because seasonal changes in roughness can produce slow and highly dispersed BTCs, dye-tracing studies may not be capable of uniquely identifying subglacial drainage system configurations. As a result, conduits may form earlier in melt seasons than previously recognized.

Snow slab avalanches primarily release by propagation of shear fractures within thin weak layers under thicker slabs. The weak layer is typically on the order of 1 mm thick and fails in mode II. In some cases, the weak layer is thicker and there may be a need to consider slope-normal deformation as part of the energy condition prior to rapid propagation. In this paper, field measurements from shear fracture initiation and high-speed films are combined to consider the effects of slope-normal deformation on bending of the slab prior to propagation and its relation to the propagation condition. Slab bending is modelled using two limiting cases: (1) a uniformly loaded beam supported by a deforming weak layer, analogous to a Winkler foundation, and (2) a uniformly loaded unsupported cantilevered beam. The experimental and analytical results suggest that slab bending prior to fracture initiation is small or negligible. Two previous approaches to modelling slab avalanche initiation involving slab bending are discussed. Both models proposed strong slab-bending effects prior to initiation, which conflicts with our results. Our field observations and modelling both show that strong bending is a dynamic effect following slope-parallel weak-layer fracture initiation.