Two medium-depth ice cores were retrieved from Berkner Island by a joint project between the Alfred-Wegener-Institut and the British Antarctic Survey in the 1994/95 field season. A 151m deep core from the northern dome (Reinwarthhöhe) of Berkner Island spans 700 years, while a 181 m deep core from the southern dome (Thyssenhöhe) spans approximately 1200 years. Both cores display clear seasonal cycles in electrical conductivity measurements, allowing dating by annual-layer counting and the calculation of accumulation profiles. Stable-isotope measurements (both δ18O and δD), together with the accumulation data, allow us to estimate changes in climate for most of the past millennium: the data show multi-decadal variability around a generally stable long-termmean. In addition, a full suite of major chemistry measurements is available to define the history of aerosol deposition at these sites: again, there is little evidence that the chemistry of the sites has changed over the past six centuries. Finally, we suggest that the southern dome, with an ice thickness of 950 m, is an ideal site from which to gain a climate history of the late stages of the last glacial and the deglaciation for comparison with the records from the deep Antarctic ice cores, and with other intermediate-depth cores such as Taylor Dome and Siple Dome.

Bakaninbreen is a polythermal glacier in southern Spitsbergen, Svalbard, which surged between 1985 and 1995. For 9 days in spring 1987, when the surge front was travelling at ∼2.5-3.0 m d–1, three single-component geophones and two accelerometers were deployed in a T-shaped array immediately downstream of the surge front to record seismic emissions. The events were characterized by their waveforms and spectral content. At least three different categories have been identified: impulsive P- and S-waveforms, surface P- and S-wave trains, and harmonic (75-130 Hz) events. We interpret the impulsive events to originate at the base of the glacier, at or downstream of the surge front; the surface P- and S-wave trains from near-surface brittle fracture associated with the surface expression of the surge front itself; and the harmonic events from deep sources that involve resonance in a water-filled fracture, associated with the base of the surge front. We believe the basal events are related to the activation of stagnant ice downstream of the surge front, which allows water to access the bed and provides the mechanism for its propagation.

The major river systems of India, i.e. the Indus, Ganga and Brahmaputra river systems originating in the Himalayan region, are considered the lifeline of the Indian subcontinent. The main sources maintaining the flow of the Himalayan rivers are snow/glacial melt runoff, rainfall runoff and base flow. The Beas River originates from Beas Kund Glacier in the Himalayan region and flows down to join the Sutlej River, which is a tributary of the Indus River system. In the present study two approaches, namely hydrologic modelling and isotope analysis, have been applied to estimate the contribution of snow and glacier melt. Samples of streamflow, rainfall and snow for isotopic analysis were collected daily from April to September and weekly from October to March during 2010 and 2011. The isotope analysis of samples reveals that the snow/glacier melt contribution to the Beas River at Manali is 50% of the total flow during these 2 years. Snowmelt runoff modelling has been carried out using the SNOWMOD model, and the snow/glacier melt runoff contribution is calculated to be 52% of the total flow during the same period. These findings indicate that the results obtained from the two approaches are similar.

In the past, ionic analyses of deep ice cores tended to consist of a few widely spaced measurements that indicated general trends in concentration. the ion-chromatographic methods widely used provide well-validated individual data, but are time-consuming. the development of continuous flow analysis (CFA) methods has allowed very rapid, high-resolution data to be collected in the field for a wide range of ions. In the European Project for Ice Coring in Antarctica (EPICA) deep ice-core drilling at Dome C, many ions have been measured at high resolution, and several have been analyzed by more than one method. the full range of ions has been measured in five different laboratories by ion chromatography (IC), at resolutions of 2.5–10 cm. In the field, CFA was used to measure the ions Na+, Ca2+, nitrate and ammonium. Additionally, a new semi-continuous in situ IC method, fast ion chromatography (FIC), was used to analyze sulphate, nitrate and chloride. Some data are now available to 788 m depth. In this paper we compare the data obtained by the three methods, and show that the rapid methods (CFA and FIC) give an excellent indication of trends in ionic data. Differences between the data from the different methods do occur, and in some cases these are genuine, being due to differences in speciation in the methods. We conclude that the best system for most deep ice-core analysis is a rapid system of CFA and FIC, along with in situ meltwater collection for analysis of other ions by IC, but that material should be kept aside for a regular check on analytical quality and for more detailed analysis of some sections.

This paper contains statistical analyses of parameters to characterize starting zones of destructive avalanches which have resulted from clear-cut logging in British Columbia, Canada. Data from 76 avalanche sites in the Coast Mountains (western British Columbia) and the Columbia Mountains (eastern British Columbia) are analyzed. The parameters include a selection which characterize snow supply (related to potential avalanche frequency), avalanche magnitude and those which are known to affect avalanche formation including terrain features, vegetation density, vegetation height and ground surface roughness. The results provide the data framework for possibly preventing future disasters by altering logging plans.

Area changes at 39 of the widest Greenland marine-terminating glacier outlets are measured in consecutive annual end-of-melt-season Moderate Resolution Imaging Spectroradiometer (MODIS) scenes spanning ten annual intervals (2000–10). The rates of cumulative area change for glaciers and ice shelves are well represented by linear least-squares fits, R = –0.99 and R = –0.94, with average rates of –70 and –65km2 a–1, respectively. Collectively, during this decade, the 39 glaciers lost a cumulative area of 1368 km2. More than three-quarters of the total area change occurred north of 72˚N. The largest 11-year area change for a single glacier during the survey period is the 311km2 loss at Humboldt Glacier. The largest annual change for a single glacier was extreme compared with the others, where Petermann glacier retreated 17km between 3 and 5 August 2010. For the 10 year sample, on average, the count of glaciers retreating is twice that advancing. A larger distinction is evident considering area change, with the ratio of retreat and advance, on average, nine times the gain. For glaciers with ice shelves, we find no year with collective area gain. The area change data from this study are posted at: http://bprc.osu.edu/~jbox/data/GAC/

In the framework of the International Trans-Antarctic Scientific Expedition (ITASE) programme, France and Italy carried out a traverse along one west–east and two north–south transects in East Antarctica from November 2001 to January 2002. Eighteen shallow snow–firn cores were drilled, and surface snow samples were collected every 5km along the traverse. Firn temperatures were measured in boreholes down to 30 m. The cores were analyzed for β radioactivity to obtain snow accumulation-rate data. The surface snow samples were analyzed for δ18O to correlate isotopic values with borehole temperatures. Multiple regression analysis shows a global near-dry-adiabatic lapse rate and a latitudinal lapse rate of 1.05˚C(˚ lat. S)–1, in the Dome C drainage area. Analysis of firn temperatures reveals a super-adiabatic lapse rate along the ice divide between Talos Dome and the Southern Ocean coast, and in some sectors along the ice divide between the Astrolabe Basin and D59. Snow accumulation rates and firn temperatures show warmer temperatures and higher accumulation values close to the ice divides extending from Talos Dome and Dome C to the Southern Ocean. The spatial pattern of data is linked with a katabatic-wind-source basin and moisture-source region.

An annually dated ice core recovered from South Pole (2850 m a.s.l.) in 1995, that covers the period 1487–1992, was analyzed for the marine biogenic sulfur species methanesulfonate (MS). Empirical orthogonal function analysis is used to calibrate the high-resolution MS series with associated environmental series for the period of overlap (1973–92). Utilizing this calibration we present a ~500 year long proxy record of the polar expression of the El Niño–Southern Oscillation (ENSO) and southeastern Pacific sea-ice extent variations. These records reveal short-term periods of increased (1800–50, 1900–40) and decreased sea-ice extent (1550–1610, 1660–1710, 1760–1800). In general, increased (decreased) sea-ice extent is associated with a higher (lower) frequency of El Niño events.

Automated classification of sea-ice types in Synthetic Aperture Radar (SAR) imagery is complicated by the class-dependent decrease of backscatter intensity with Incidence Angle (IA). In the log-domain, this decrease is approximately linear over the typical range of space-borne SAR instruments. A global correction does not consider that different surface types show different rates of decrease in backscatter intensity. Here, we introduce a supervised classification algorithm that directly incorporates the surface-type dependent effect of IA. We replace the constant mean vector of a Gaussian probability density function in a Bayesian classifier with a linearly variable mean. During training, the classifier first retrieves the slope and intercept of the linear function describing the mean value and then calculates the covariance matrix as the mean squared deviation relative to this function. The IA dependence is no longer treated as an image property but as a class property. Based on training and validation data selected from overlapping SAR and optical images, we evaluate the proposed method in several case studies and compare to other classification algorithms for which a global IA correction is applied during pre-processing. Our results show that the inclusion of the per-class IA sensitivity can significantly improve the performance of the classifier.

The Norwegian Antarctic Research Expedition 1978-79 landed on 24 tabular icebergs and flew over many others in the South Atlantic and the Weddell Sea between latitudes 54 and 76°S. Data were obtained on surface mass balance, stratigraphy, density, 10 m temperatures, crevassing, distribution, and age. Ice thicknesses were measured by airborne radio echo-sounding.

All icebergs had experienced surface melting. However, on icebergs south of 66°S., the annual surface melting was only a few centimetres of water equivalent. The average surface mass balance was near zero. Typically the 10 m temperature had increased from about -20° C at the time of calving to -10° C. Only icebergs that had moved northwards from the continent into the west wind drift had snow temperatures close to 0° C. Internal temperatures are increased mainly by the refreezing of percolating melt and rain water. This increased the densities of the upper layers by 100 to 150 kg m-3 above those of nearby ice shelves.

All icebergs measured by radio echo-sounding showed variations in thickness of about 20% of the mean thickness. Nearly all had a convex profile across the short axis and were tilted. An average thickness/freeboard curve indicates that icebergs less than 225 m thick will have permeable layers below sea-level. The ratio of freeboard to thickness varied from 0.21 for a 100 m thick berg to 0.14 for a 350 m thick iceberg.

All icebergs showed systematic surface crevassing parallel with their sides, the crevasse intensity decreasing with distance from the edge. Icebergs with their smallest dimension greater than 400 m usually had a central zone with little crevassing. Grounded icebergs showed severe crevassing, and could not thereafter survive long periods in open water. Bottom crevasses were not detected.

Ice-penetrating radar and modelling data are presented suggesting the presence of a zone of temperate ice, water ponding or saturated sediment beneath the tongue of Taylor Glacier, Dry Valleys, Antarctica. The proposed subglacial zone lies 3–6 km up-glacier of the terminus and is 400– 1000m across. The zone coincides with an extensive topographic overdeepening to 80m below sea level. High values of residual bed reflective power across this zone compared to other regions and the margins of the glacier require a high dielectric contrast between the ice and the bed and are strongly indicative of the presence of basal water or saturated sediment. Analysis of the hydraulic equipotential surface also indicates strong convergence into this zone of subglacial water flow paths. However, thermodynamic modelling reveals that basal temperatures in this region could not exceed –7˚C relative to the pressure-melting point. Such a result is at odds with the radar observations unless the subglacial water is a hypersaline brine.

To assess the safety against failure of rock slopes in cold regions, such as high mountain areas, where stability is potentially maintained by ice in rock discontinuities, the shear strength of ice-filled rock joints was investigated in a series of direct shear-box tests. To permit control and repeatability, the experiments were conducted using simulated rock specimens. These were cast in the laboratory using high-strength concrete. Laboratory measurements showed that at a constant rate of shearing, the interface shear strength between ice and a joint surface of repeatable roughness is a function of both temperature and normal stress.

We extend the application of control methods to a comprehensive three-dimensional thermomechanical ice-sheet model, SICOPOLIS (SImulation COde for POLythermal Ice Sheets). Lagrange multipliers, i.e. sensitivities, are computed with an exact, efficient adjoint model that has been generated from SICOPOLIS by rigorous application of automatic differentiation. The case study uses the adjoint model to determine the sensitivity of the total Greenland ice volume to various control variables over a 100 year period. The control space has of the order 1.2 × 106 elements, consisting of spatial fields of basal flow parameters, surface and basal forcings and initial conditions. Reliability of the adjoint model was tested through finite-difference perturbation calculations for various control variables and perturbation regions, ascertaining quantitative inferences of the adjoint model. As well as confirming qualitative aspects of ice-sheet sensitivities (e.g. expected regional variations), we detect regions where model sensitivities are seemingly unexpected or counter-intuitive, albeit ‘real’ in the sense of actual model behavior. An example is inferred regions where sensitivities of ice-sheet volume to basal sliding coefficient are positive, i.e. where a local increase in basal sliding parameter increases the ice-sheet volume. Similarly, positive (generally negative) ice temperature sensitivities in certain parts of the ice sheet are found, the detection of which seems highly unlikely if only conventional perturbation experiments had been used. The object of this paper is largely a proof of concept. Available adjoint-code generation tools now open up a variety of novel model applications, notably with regard to sensitivity and uncertainty analyses and ice-sheet state estimation or data assimilation.

A 2500 year integration has been performed with a global coupled atmospheric-sea-ice-ocean model of intermediate complexity with the main objective of studying the climate variability in polar regions on decadal time-scales and longer. The atmospheric component is the ECBILT model, a spectral T21 three-level quasi-geostrophic model that includes a representation of horizontal and vertical heat transfers as well as of the hydrological cycle. ECBILT is coupled to the CLIO model, which consists of a primitive-equation free-surface ocean general circulation model and a dynamic-thermodynamic sea-ice model. Comparison of model results with observations shows that the ECBILT-CLIO model is able to reproduce reasonably well the climate of the high northern latitudes. The dominant mode of coupled variability between the atmospheric circulation and sea-ice cover in the simulation consists of an annular mode for geopotential height at 800 hPa and of a dipole between the Barents and Labrador Seas for the sea-ice concentration which are similar to observed patterns of variability. In addition, the simulation displays strong decadal variability in the sea-ice volume, with a significant peak at about 18 years. Positive volume anomalies are caused by (1) a decrease in ice export through Fram Strait associated with more anticyclonic winds at high latitudes, (2) modifications in the freezing/melting rates in the Arctic due to lower air temperature and higher surface albedo, and (3) a weaker heat flux at the ice base in the Barents and Kara seas caused by a lower inflow of warm Atlantic water. Opposite anomalies occur during the volume-decrease phase of the oscillation.

At least in conventional hydrostatic ice-sheet models, the numerical error associated with grounding line dynamics can be reduced by modifications to the discretization scheme. These involve altering the integration formulae for the basal traction and/or driving stress close to the grounding line and exhibit lower – if still first-order – error in the MISMIP3d experiments. MISMIP3d may not represent the variety of real ice streams, in that it lacks strong lateral stresses, and imposes a large basal traction at the grounding line. We study resolution sensitivity in the context of extreme forcing simulations of the entire Antarctic ice sheet, using the BISICLES adaptive mesh ice-sheet model with two schemes: the original treatment, and a scheme, which modifies the discretization of the basal traction. The second scheme does indeed improve accuracy – by around a factor of two – for a given mesh spacing, but $\lesssim 1$ km resolution is still necessary. For example, in coarser resolution simulations Thwaites Glacier retreats so slowly that other ice streams divert its trunk. In contrast, with $\lesssim 1$ km meshes, the same glacier retreats far more quickly and triggers the final phase of West Antarctic collapse a century before any such diversion can take place.

Since 1994 the Alfred Wegener Institute (AWI) has operated an airborne radio-echo sounding system for remote-sensing studies of the polar ice caps in Antarctica and in Greenland. It is used to map ice thicknesses and internal layernigs of glaciers, ice sheets and ice shelves, and is capable of penetrating ice thicknesses of up to 4 km. The system was designed and built by AWI in cooperation with Aerodata Flugmeßtechnik GmbH, Technische Umversitat Hamburg-Harburg and the Deutsches Zentrum fur Luft- und Raumfahrt e.V. The system uses state-of-the-art techniques, and results in high vertical (5 m) as well as along-track (3.25 m) resolution. The radar signal is a 150 MHz burst with a duration of 60 or 600 ns. The peak power is 1.6 kW, and the system sensitivity is 190 dB. The short backfire principle has been adopted and optimized for antennae used on Polar2, a Dormer 228-100 aircraft, resulting in an antenna gain of 14 dB each. Digital data recording allows further processing. The quality of the recorded data can be monitored on screen and as online analogue plots during the flight.

Vapor transport in snow and firn plays a key role in post-depositional changes of many reactive species, in mass and energy balance of large snow-covered areas, and in metamorphic changes of the crystal structure of the snow itself. While conventional estimates of vapor transport and sublimation rate are based on diffusion of water vapor from the snow surface to the atmosphere, ventilation (airflow through interstitial pore spaces) in snow and firn can affect the top several meters of firn, yet the effects of ventilation on sublimation rates in firn have not been previously investigated. In this paper we present the first calculations of sublimation rates due to advective/ diffusive airflow and vapor transport in the snow and firn. Airflow velocities, vapor transport and sublimation rates are calculated using a two-dimensional finite-element model. the airflow patterns induce regions of sublimation and regions of condensation within the snow and firn. Because small surface sastrugi move in time as the snow surface profile changes due to wind redeposition, the subsurface flow and regions of condensation and sublimation will also likely change on a time-scale of days to weeks. However, if the roughness features are such that they move very little over time, it is likely that the regions of condensation and sublimation will have a noticeable effect on the microstructure of the firn. While the highest mass-transfer rates occur near the surface, the depth of vapor transport and phase change depends on firn properties, wavelength of the snow surface roughness, and temperature.

Observations during February and March 1980 of structures in 66 separate floes in Weddell Sea pack ice show widespread occurrence of frazil ice in amounts not previously reported in sea ice of comparable age and thickness in the Arctic. It is estimated that as much as 50% of the total ice production in the Weddell Sea is generated as frazil. Average floe salinities also appear higher than those of their Arctic counterparts. Comparative studies of fast ice at 28 locations in McMurdo Sound show this ice to be composed almost entirely of congelation ice that exhibits crystalline textures and orientations that are similar to those observed in Arctic fast ice. However, average fast-ice salinities in McMurdo Sound are higher than those reported for Arctic fast ice of comparable age and thickness.

King George Island is the largest of the South Shetland Islands, close to the tip of the Antarctic Peninsula. The annual mean temperature on the island has increased by 1°C during the past three decades, and the ice cap that covers the majority of the island is sensitive to climatic change. We present data from two field campaigns (1997 and 2007): 700 km of global positioning system (GPS) and ground-penetrating radar (GPR) profiles were collected on Arctowski Icefield and on the adjacent central part. The data were analysed to determine the surface and bed topography and the thermal regime of the ice. Average ice thickness is 250 m and maximum thickness is 420 m. The GPR profiles show isochrones throughout the ice cap which depict the uparching of Raymond bumps beneath or close to the ice divides. A water table from percolation of meltwater in the snowpack shows the firn-ice boundary at ∼ 3 5 m depth. The firn layer may be temperate due to the release of latent heat. In the area below 400ma.s.l., backscatter by water inclusions is abundant for ice depths below the water table. We interpret this as evidence for temperate ice. Scatter decreases significantly above 400 m. Ice temperatures below the water table in this part of the ice cap are subject to further field and modelling investigations.

We investigate whether the recent slow-down of Whillans Ice Stream (WIS), West Antarctica, may lead to its complete stoppage in the near future, using a numerical model. Basal resistance to ice-stream motion is represented by a continuous till layer whose strength changes in response to basal melting and freezing. We implement a basal drainage system, which acts to hinder ice-stream stoppage through supply of extra water to those parts of the bed that are experiencing freezing. The ice module is a standard flowline model with parameterized key out-of-plane effects (Raymond, 1996). The most important result of our modeling effort is that we never obtained a slow-down that would not be followed by a complete stoppage within <100 years. WIS slow-down can be avoided in the model if large basal water-input rates are assumed, to satisfy basal freezing. In comparison, the tested perturbations in WIS width and mass balance had a relatively small effect on the tendency of the simulated ice stream to slow down. These results underscore the need for more quantitative constraints on the efficiency of sub-ice-stream water drainage. We conjecture that the present-day slow-down of WIS will evolve to shut-down in the next few decades, unless an addition of basal water prevents freeze-on-driven bed strengthening.