A linear two-level atmospheric model is employed to study the influence of ice-sheet topography on atmospheric stationary waves. In particular, the stationary-wave-induced temperature anomaly is considered locally over a single ice-sheet topography, which is computed using the plastic approximation. It is found that stationary waves induce a local cooling which increases linearly with the ice volume for ice sheets of horizontal extents smaller than ∼1400 km. Beyond this horizontal scale, the dependence of stationary-wave-induced cooling on the ice volume becomes gradually weaker. For a certain ice-sheet size, and for small changes of the surface zonal wind, it is further shown that the strength of the local stationary-wave-induced cooling is proportional to the basic state meridional temperature gradient multiplied by the vertical stratification in the atmosphere. These results are of importance for the nature of the feedback between ice sheets and stationary waves, and may also serve as a basis for parameterizing this feedback in ice-sheet model simulations (e.g. through the Pleistocene glacial/interglacial cycles).

Digital optical televiewing (OPTV) of hot-water-drilled boreholes is evaluated as a technique for the investigation of englacial ice and debris structures on the basis of six boreholes drilled in the terminus region of midre Lovénbreen, Svalbard. The resulting OPTV logs successfully reveal several visually distinctive englacial ice properties and deformation structures (e.g. oblique englacial fractures imaged here for the first time). Combining these OPTV logs with surface mapping has resulted in the identification of eight separate structural elements, several of which can be interpolated onto 3-D grids at a node spacing of 1 m vertically and 10 m horizontally. Basally derived englacial sediment layers are also found to be intercalated with primary stratification, elevated into near-vertical planes around a central fold axis by large-scale lateral folding. The analysis also allows supraglacial longitudinal debris ridges to be subclassified into two types: a previously described (type-I) form, which are the exposed fold axes of large-scale lateral folds, and a new (type-II) form experiencing secondary deformation by small-scale horizontal folding in association with vertical displacements across arcuate shear planes in response to longitudinally compressive stresses near the glacier terminus. Although using boreholes to investigate glacier structure is limited (e.g. by parallelism with vertical planes), applying OPTV to multiple boreholes at midre Lovénbreen has successfully revealed a range of 3-D structural elements at high spatial resolution. As such, interpolating between multiple OPTV logs overcomes many of the problems associated with interpretations made solely on the basis of surface-based structural mapping, and combining the two techniques represents a powerful glaciological tool.

Ice cores are remarkably faithful recorders of past climate, providing multiply duplicated reconstructions with small and quantifiable uncertainties. Ice core reconstructions in general do not rely on assumed quantitative time-invariance of empirical calibrations between climate and sedimentary characteristics, but instead rely on assuming little more than the constancy of physical law over time. The history of some of the discoveries that allow these ice-core climatic reconstructions is instructive for students, citizens and policymakers. Much important ice-core science was published in the Journal of Glaciology and the Annals of Glaciology, with scientific impacts much larger than indicated by mere citation counts.

Using a Monte Carlo (MC) method, we determine the accumulation-rate profile along a flowband, the influx of ice into the upstream end of the flowband and the age of an internal layer. The data comprise the depth profile of the internal layer, a few velocity measurements at the surface and the average accumulation at one location. The data in our example were collected at Taylor Mouth, a flank site off Taylor Dome, Antarctica. We present three alternative formulations of this inverse problem. Depending on the formulation used, this particular inverse problem can have up to four solutions, each corresponding to a different spatial accumulation-rate pattern. This study demonstrates the ability of a MC method to find several solutions to this inverse problem, and how to use a Metropolis algorithm to determine the probability distribution of each of these different solutions. The only disadvantage of the MC method is that it is computationally more expensive than other inverse methods, such as the Gradient method.

The field observations of G.S. Boulton stimulated widespread interest in deformable beds. Shear resistance of till in its critical state is insensitive to strain rate and increases linearly with effective pressure. During unsteady deformation, pseudo-viscous shear resistance can be caused by dilation of consolidated tills and resultant pore-pressure decline. This effect is probably uncommon, however, because susceptible tills of low hydraulic diffusivity are also those least likely to consolidate significantly during effective-pressure transients. Stick-slip motion at Whillans Ice Stream, Antarctica, indicates that its basal till must weaken during rapid slip and strengthen during longer periods of slower slip. Recurrence intervals for rapid-slip episodes there (6-18hours) indicate that till-strength variations, if driven by changes in pore pressure either related or unrelated to basal freezing, are focused in the uppermost several centimeters of the bed. Ploughing of grains at the bed surface and associated excess pore pressures in adjacent till can account for rate-weakening during rapid slip, with pore-pressure decay causing strengthening between slip episodes. By promoting shallow, sluggish subglacial water flow and low effective pressure, soft beds may help sustain themselves by slowing their own transport. Soft-bed shear resistance, kinematics and continuity are problems rooted in subglacial hydrology.

Samples were collected from a snow pit and shallow firn core near Kahiltna Pass (2970 m a.s.l.), Denali National Park, Alaska, USA, in May 2008. The record spans autumn 2003 to spring 2008 and reveals clusters of ice layers interpreted as summertime intervals of above-freezing temperatures. High correlation coefficients (0.75–1.00) between annual ice-layer thickness and regional summertime station temperatures for 4 years (n = 4) indicate ice-layer thickness is a good proxy for mean and extreme summertime temperatures across Alaska, at least over the short period of record. A Rex-block (aka high-over-low) pattern, a downstream trough over Hudson Bay, Canada, and an upstream trough over eastern Siberia occurred during the three melting events that lasted at least 2 weeks. About half of all shorter melting events were associated with a cut-off low traversing the Gulf of Alaska. We hypothesize that a surface-to-bedrock core extracted from this location would provide a high-quality record of summer temperature and atmospheric blocking variability for the last several hundred years.

We use ice-penetrating radar data across grounding lines of Siple Dome and Roosevelt Island, Antarctica, to measure the spatial pattern, magnitude and duration of sub-ice-shelf melting at these locations. Stratigraphic layers across the grounding line show, in places, a large-amplitude downwarp at, or slightly downstream of, the grounding line due to sub-ice-shelf basal melting. Localized downwarping indicates that melting is transient; melt rates, or the grounding line position, have changed within a few hundred years in order to produce the observed stratigraphy. Elsewhere, no melt-related stratigraphic signature is preserved. In part, heterogeneity in the amount of sub-ice-shelf melt is due to regional circulation patterns in the sub-shelf cavity, but local (on the order of tens of kilometers) heterogeneity in the melt pattern may reflect small differences in the shape of the ice-shelf base at the grounding line. We find that all of the grounding lines crossed have been in place for at most ∼400 years.

The ice streams feeding the Ross Ice Shelf, Antarctica, have large tidally modulated (sinusoidal and stick–slip) flow, but the interaction with the ice shelf is poorly understood. We show that the flow of the Ross Ice Shelf front, up to ∼650 km from the ice streams, exhibits smooth, sinusoidal motions corresponding to tidal modulation. These observations suggest a possible linking of the ice shelf with the ice streams to form a unified system that responds to small perturbations in stresses associated with ocean tides. If this is the case, the presence of the sinusoidal motion but the absence of stick–slip motion suggests there is damping of very high-frequency signals. The dissimilar signatures of the motions observed in the ice streams and at the front of the ice shelf present challenges to model development aimed at understanding the dynamics of coupled ice-stream/ice-shelf flow and the movement of ice across grounding lines.

A new vibration apparatus for measuring the shear strength of snow has been designed and fabricated. The force applied to a snow block is calculated using Newton’s second law. Results from this apparatus concerning the dependence of the shear strength on snow density, overburden load and strain rate are in reasonable agreement with those obtained from the work of previous researchers. Snow densities ranged from 160 to 320 kg m−3. The overburden load and strain rate ranged from 1.95 × 10−1 to 7.79 × 10−1 kPa and 2.9 × 10−4 to 9.1 × 10−3 s−1 respectively.

The northwestern sector of the Amery Ice Shelf, East Antarctica, has a layered structure, due to the presence of both meteoric ice and a marine ice layer resulting from sub-shelf freezing processes. Crystal orientation fabric and grain-size data are presented for ice cores obtained from two boreholes ˜70 km apart on approximately the same flowline. Multiple-maxima crystal orientation fabrics and large mean grain sizes in the meteoric ice are indicative of stress relaxation and subsequent grain growth in ice that has flowed into the Amery Ice Shelf. Strongly anisotropic single-maximum crystal orientation fabrics and rectangular textures near the base of the ˜200 m thick marine ice layer suggest accretion occurs by the accumulation of frazil ice platelets. Crystal orientation fabrics in older marine ice exhibit vertical large circle girdle patterns, influenced by the complex stress configurations that exist towards the margins of the ice shelf. Post-accumulation grain growth and fabric development in the marine ice layer are restricted by a high concentration of brine and insoluble particulate inclusions. Differences in the meteoric and marine ice crystallography are indicative of the contrasting rheological properties of these layers, which must be considered in relation to large-scale ice-shelf dynamics.

Automated digital cameras were installed in May–June 2007 beside major West Greenland marine-terminating glaciers as part of the Extreme Ice Survey (EIS). EIS cameras began imaging the lowest 4 km2 of the glacier at hourly intervals throughout sunlit periods of the year. This study presents the development of techniques for quantifying glacier velocity from a single camera perspective. A Multi-Image/Multi-Chip matching procedure yields higher matching skill than conventional matching, and facilitates false-match rejection via a clustering scheme. The matching of motionless on-land features facilitates compensating camera motion. Ray projection to a known terrain elevation allows the assigning of scale to convert pixel displacements to velocity units. With the 10.2-megapixel camera system, velocities on relatively fast glaciers can be resolved at distances up to ∼4 km. At a distance of 2 km, a demonstrated precision of ∼0.5 pixels yields a ∼0.5 m footprint size. Daily velocities indicate progressive multi-day velocity accelerations associated with calving. Deceleration trends are associated with the regrowth of resistive stress after major calving. The higher observation frequency available to terrestrial photogrammetry indicates higher observed intra-seasonal velocity range than observable by the at-best weekly satellite snapshots.

Air bubbles in ice cores play an essential role in climate research, not only because they contain samples of the palaeoatmosphere, but also because their shape, size and distribution provide information about the past firn structure and the embedding of climate records into deep ice cores. In this context, we present profiles of average bubble size and bubble number for the entire EDML (Antarctica) core and the top 600 m of the EDC (Antarctica) core, and distributions of bubble sizes from selected depths. The data are generated with an image-processing framework which automatically extracts position, orientation, size and shape of an elliptical approximation of each bubble from thick-section micrographs, without user interaction. The presented software framework allows for registration of overlapping photomicrographs to yield accurate locations of bubble-like features. A comparison is made between the bubble parameterizations in the EDML and EDC cores and data published on the Vostok (Antarctica) ice core. The porosity at the firn/ice transition is inferred to lie between 8.62% and 10.48% for the EDC core and between 10.56% and 12.61 % for the EDML core.

Analyses of shallow cores obtained at the European Project for Ice Coring in Antarctica (EPICA) drilling site Kohnen station (75°00′ S, 00°04′ E; 2892 m a.s.l.) on the plateau of Dronning Maud Land reveal the presence of conserved snow dunes in the firn. In situ observations during three dune formation events in the 2005/06 austral summer at Kohnen station show that these periods were characterized by a phase of 2 or 3 days with snowdrift prior to dune formation which only occurred during high wind speeds of >10 m s-1 at 2 m height caused by the influence of a low-pressure system. The dune surface coverage after a formation event varied between 5% and 15%, with a typical dune size of (4 ± 2) m × (8 ± 3) m, a maximum height of 0.2 ± 0.1 m and a periodicity length of about 30 m. The mean density within a snow dune varied between 380 and 500 kg m-3, whereas the mean density at the surrounding surface was 330 ± 5 kgm-3. The firn cores covering a time-span of 22 ± 2 years reveal that approximately three to eight events per year occurred, during which snow dunes had been formed and were preserved in the firn.

Over the period of modern satellite observations, Arctic sea-ice extent at the end of the melt season (September) has declined at a rate of >11% per decade, and there is evidence that the rate of decline has accelerated during the last decade. While climate models project further decreases in sea- ice mass and extent through the 21st century, the model ensemble mean trend over the period of instrumental records is smaller than observed. Possible reasons for the apparent discrepancy between observations and model simulations include observational uncertainties, vigorous unforced climate variability in the high latitudes, and limitations and shortcomings of the models stemming in particular from gaps in understanding physical process. The economic significance of a seasonally sea-ice-free future Arctic, the increased connectivity of a warmer Arctic with changes in global climate, and large uncertainties in magnitude and timing of these impacts make the problem of rapid sea-ice loss in the Arctic a grand challenge of climate science. Meaningful prediction/projection of the Arctic sea-ice conditions for the coming decades and beyond requires determining priorities for observations and model development, evaluation of the ability of climate models to reproduce the observed sea-ice behavior as a part of the broader climate system, improved attribution of the causes of Arctic sea-ice change, and improved understanding of the predictability of sea-ice conditions on seasonal through centennial timescales in the wider context of the polar climate predictability.

We present a newly developed high-resolution frequency-modulated continuous-wave (FMCW) radar system for sounding the sub-ice topography. The system, working in the frequency range from 200 to 400 MHz, was developed to measure thickness and internal layers in cold ice with a resolution better than 1 m. This system has the potential to measure accumulation rates, an important input for improved knowledge of the mass balance of polar ice sheets. First measurements for the test, calibration and optimization of the new ice sounder were made in December 2005 near the Ellsworth Mountains, specifically in Patriot Hills, West Antarctica, at 80°18′ S, 81°22′ W. The complete radar system was installed on a ski-mounted support frame, and towed by hand across the ice surface. The measurement results show the capability of this system to measure ice thickness up to 1000 m and to define internal layers within the ice body.

The deformation of the ground surface that is produced by frost heave has motivated almost a century of concerted laboratory, field and theoretical studies. Well before the development of equipment capable of resolving the microscopic films that support liquid transport towards growing ice lenses, early investigators predicted their occurrence and noted their importance. Idealized experiments continue to prompt theoretical advances that have been combined to develop predictive models for the macroscopic frost-heave characteristics that are seen in the field. This contribution highlights steps on the road to our current understanding of the physical interactions that control how ice forms in unconsolidated porous media.

Modulation of the flow of Rutford Ice Stream, Antarctica, has been reported previously at semi-diurnal, diurnal, fortnightly and semi-annual periods. A model that includes non-linear response to tidal forcing has been shown to fit closely observations at fortnightly frequencies. Here we examine that model further and test its fit at a larger set of observed frequencies, including the large semi-annual displacement. We show analytically that, when forced by major tidal terms, the model (using a basal shear stress exponent m = 3) predicts several discrete response periods from 4 hours to 0.5 years. We examine a 1.5 year GPS record from Rutford Ice Stream and find that the model, when forced by a numerical tide model, is able to reproduce the reported semi-annual signal. We confirm that about 5% of the mean flow is due solely to the (m = 3) non-linear response to tidally varying basal shear stress. Our best-fitting set of model parameters is similar to those originally reported using a much shorter data record, although with noticeably improved fit, suggesting these parameters are robust. We find that m ≈ 3 fits the data well, but that m ≈ 2 does not. Furthermore, we find that a small variation in flow over the 18.6 year lunar node tide cycle is expected. Fits to semi-diurnal and diurnal terms remain relatively poor, possibly due to viscoelastic effects that are not included in the model and reduced GPS data quality at some discrete periods. For comparison, we predict the response of Bindschadler Ice Stream and Lambert Glacier and show, given identical model parameters, a similar response pattern but with ∼1–2 orders of magnitude smaller variability; these may still be measurable and hence useful in testing the applicability of this model to other locations.

Current trends show a rise in Arctic surface and air temperatures, including over the Greenland ice sheet where rising temperatures will contribute to increased sea-level rise through increased melt. We aim to establish the uncertainties in using satellite-derived surface temperature for measuring Arctic surface temperature, as satellite data are increasingly being used to assess temperature trends. To accomplish this, satellite-derived surface temperature, or land-surface temperature (LST), must be validated and limitations of the satellite data must be assessed quantitatively. During the 2008/09 boreal winter at Summit, Greenland, we employed data from standard US National Oceanic and Atmospheric Administration (NOAA) air-temperature instruments, button-sized temperature sensors called thermochrons and the Moderate Resolution Imaging Spectroradiometer (MODIS) satellite instrument to (1) assess the accuracy and utility of thermochrons in an ice-sheet environment and (2) compare MODIS-derived LSTs with thermochron-derived surface and air temperatures. The thermochron-derived air temperatures were very accurate, within 0.1 ± 0.3°C of the NOAA-derived air temperature, but thermochron-derived surface temperatures were ∼3°C higher than MODIS-derived LSTs. Though surface temperature is largely determined by air temperature, these variables can differ significantly. Furthermore, we show that the winter-time mean air temperature, adjusted to surface temperature, was ∼11°C higher than the winter-time mean MODIS-derived LST. This marked difference occurs largely because satellite-derived LSTs cannot be measured through cloud cover, so caution must be exercised in using time series of satellite LST data to study seasonal temperature trends.