We introduce a method to construct an ice thickness chart using two satellite datasets and a one-dimensional high-resolution thermodynamic snow and sea-ice model (HIGHTSI). Thin-ice thickness up to 40 cm is retrieved using MODIS ice surface temperature data. For thicker ice, thickness is retrieved from a combination of SAR image and background level ice thickness chart Hi provided by HIGHTSI. Because Hi is inherently static, we utilize passive microwave data based sea-ice concentration charts to introduce mesoscale ice dynamics into Hi. In order to create the thickness chart, the value of modified Hi is scaled by a factor which depends on the magnitude of backscattering coefficient σ0. The proposed method is effective in cold winter conditions. We calculated a series of ice thickness charts covering the Kara Sea, Arctic Russia, during winter 2008/09. The ice thickness charts were validated against Russian ice charts. Our method gave realistic results in the thickness range 30–120 cm and was uncertain for the detection of thinner-ice areas. Due to the limitations in our reference dataset, more validation work is needed to establish the accuracy in more detail. Our thickness charts can be used for operational purposes and in climate studies.

Antarctic precipitation estimations derived from several new sources are examined in comparison to results found previously. The availability of analyzed atmospheric datasets has been a significant and beneficial tool for atmospheric and climate research for a broad range of research interests. This is particularly true for the polar regions, where the observational arrays are sparsely distributed. in high southern latitudes, a comprehensive assimilation of all available observations, including satellite data, is necessary for an accurate depiction of the atmospheric circulation. Recent st udies have found the operational analyses of the European Centre for Medium-range Weather Forecasts to be superior to those of other weather-forecasting centers in depicting the large-scale atmospheric circulation patterns over Antarctica. “Re-analysis” programs at major weather-forecasting centers have produced atmospheric numerical analyses using a “frozen” data-assimilation system. These projects have also derived precipitation and evaporation fields using an ensemble of short-term forecasts. From these new sources, Antarctic Ρ - E (precipitation minus evaporation/sublimation) is compared and evaluated against the long-term glaciological synthesis, as well as results from previous studies. The comparisons indicate significant regional disagreements exist between P — E from the re-analysis forecasts and the glaciological data. For the ensemble forecasting method, the continental-average evaporation is the largest area of uncertainty and differs by an order of magnitude between the rc-analysis datasets. This finding supports the use of the atmospheric moisture budget for determining P — E collectively in atmospheric diagnostic studies for Antarctica.

Cirques in the Rassepautasjtjåkka massif currently lack glaciers and the geomorphology indicates that no glaciers occupied the cirques during the Holocene. The current climatic conditions in the cirques can be assessed using available climatic data; air temperature at Rassepautasjtjåkka, summer and winter balances of adjacent glaciers, and general precipitation patterns in northern Sweden. The data suggest that either a significant change in precipitation and wind regime or a moderate change in temperature is required to initiate a cirque glacier in the massif. Formation of a wet-based erosive glacier requires warmer winters with higher accumulation rates, equivalent to a more maritime influence in the area. Studies of current atmospheric circulation suggest that strong west- east circulation, associated with a northerly position of the polar front, is favourable for increased accumulation. Using typical erosion rates from present glaciers, we see that ~ 10% of the last 3 Myr may be required for forming the Rassepautasjtjåkka cirques. This is a significant portion of time since most of the glacial cycles are spent in states of interglacials, maximum glaciation or mountain-based glaciation. Marine sediments from the Norwegian Sea provide indications of minor glaciations back to ~12.6 Myr and, hence, cirque-formation periods are not restricted to Quaternary. Thus, it is possible that many cirque forms have a much longer history than previously recognized.

Measurements made by an underwater glider deployed near the Ross Ice Shelf were used to identify the presence of Ice Shelf Water (ISW), which is defined as seawater with its potential temperature lower than its surface freezing point temperature. Properties logged by the glider included in situ temperature, electrical conductivity, pressure, GPS location at surfacings and time. For most of the first 30 recorded dives of its deployment, evidence suggests the glider was prevented from surfacing due to being under the ice shelf. For dives under the ice shelf, farthest from the ice shelf front, ISW layers of varying thicknesses and depth locations were observed; between 2 m thick (centred at 231 m depth) to >93 m thick (centred at >360 m). For dives under the ice shelf, close to the ice shelf front, either no ISW was observed or ISW layers were centred at shallower depths (116–127 m). Thicker ISW layers (e.g. up to 250 m thickness centred at 421 m) were observed for some glider dives in open water in front of the Ross Ice Shelf. No in situ supercooling (water colder than the pressure-dependent freezing point temperature) was observed.

In 1987 an ice core to the bedrock at a depth of 85.6 m was drilled at the top of Høghetta ice dome in northern Spitsbergen. Chronology of the ice core was examined by tritium and 14C methods showing time gap at about 50 m depth. The age of three bottom ice samples was determined as 4150–5670 year B.P. by 14C method done for frozen bacteria colonies and a frozen petal. This chronology and negative bottom temperature of −9.4°C suggest that glaciers in Spitsbergen shrank considerably during the hypsithermal. The pH of melt-water samples lower than 5.0 corresponds well to large northern hemispheric volcanic eruptions during the last 300 years. Increase of acidity from 30 m depth to the surface may reflect the spread of air pollution to the Arctic during the past 200 years. On the basis of ice-core analyses on electrical conductivity, pH, chemical composition and air bubble pattern, climate and environment in Spitsbergen during the last 6000 years are discussed.

A new ocean wave/sea-ice interaction model is proposed that simulates how a directional wave spectrum evolves as it travels through an arbitrary finite array of circular ice floes, where wave/ ice dynamics are entirely governed by wave-scattering effects. The model is applied to characterize the wave reflection and transmission properties of a strip of ice floes, such as an ice edge band. A method is devised to extract the reflected and transmitted directional wave spectra produced by the array. The method builds upon an integral mapping from polar to Cartesian coordinates of the scattered wave components. Sensitivity tests are conducted for a row of floes randomly perturbed from a regular arrangement. Results for random arrays are generated using ensemble averaging. A realistic ice edge band is then reconstructed from field experiment data. Simulations show good qualitative agreement with the data in terms of transmitted wave energy and directional spreading. In particular, it is observed that short waves become isotropic quickly after penetrating the ice field.

The characteristics of summer energy budgets in the ablation zones during the summer at the surface of two glaciers in the Antarctic Peninsula are investigated and compared to the findings of previous studies. The study areas are located on King George Island (62° S) and in Marguerite Bay (68° S). The summer energy balance was computed from automatic weather station data. The results reveal that turbulent heat fluxes dominate over radiation balance in Marguerite Bay whereas on King George Island ablation is driven by net radiation. The summer energy balance on King George Island reflects the more maritime subpolar climate along the northwest tip of the peninsula in contrast to a more continentally toned polar climate in Marguerite Bay and areas further south. The terms of the energy balances are partitioned very differently but ad-vection from northerly directions causes the highest summer snowmelt rates at both study sites. It is concluded that sensitivity studies should consider not only the mean air-temperature increase, but also possible changes in other climate parameters.

Blowing snow was produced artificially in a cold wind-tunnel, and various measurements were conducted including particle diameters, concentrations, saltation lengths heat transport and electric charge. The mean diameter of blowing snow particles decreased only slightly with increasing height; in the saltation layer, standard deviation was large and velocities were scattered in a wide range, suggesting the complex dynamic process on taking-off. The mean saltation length ranged from a few cm to 40 cm increasing with wind velocity.

When wind blew without snow drifting, the static air pressure on the snow surface was smaller at higher levels, the vertical pressure gradient being negative. The pressure gradient became positive when blowing snow was initiated eg +9.6 Pa/m at 11.2 m/s and -8.3 °C. The magnitude of à downward force acting on a saltating snow partice caused by the pressure gradient was not large enough to explain the downward acceleration found from photographic analyses of particle trajectories.

Blowing snow particles were charged negatively the magnitude of charge increased with lowering temperature. Increase in vertical heat transfer was found in blowing snow by measuring the temperature of the air at various levels; the increase is reflected on that in the apparent turbulent diffusion coefficient.

Detailed ionic analyses of Dyer Plateau snow show that major soluble impurities in snow consist of sodium (Na+), chloride (Cl−), nitrate (NO3−), sulfate (SO42−), and acidity (H+). The ratios of Na+ to Cl− concentrations are close to that of sea water, indicating little or no fractionation of sea-salt aerosols. The analyses of core sections from three sites along a 10 km transect show that local spatial variation of snow chemistry in this area is minimal and that temporal (decadal, inter-annual and sub-annual) variations in snow chemistry are very well preserved.

Anion analyses of the upper 181 m section of two 235 m ice cores yield a data set of 485 years (1505-1989) of annual snow accumulation and fluxes of Cl−, NO3−, and non-sea-salt (nss) SO42−. No significant long-term trends are observed in any of the anion fluxes. This is consistent with other Antarctic ice-core records showing no significant anthropogenic atmospheric pollution in the high southern latitudes. Linear regression analysis shows that Cl− flux is independent of snow-accumulation rate. Significant positive correlations are found between accumulation rate and both NO3− flux and background nss-SO42− flux. These results suggest that dry deposition is primarily responsible for air-to-ground Cl− flux while wet deposition dominates the NO3− and nss-SO42− flux (≥90% and ≥75%, respectively). The nss-S042− fluxes provide a chronology of explosive volcanic emissions reaching the Antarctic region for the past 485 years.

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.

The response of glaciers to changing climate is explored with an atmosphere/glacier hierarchical modeling approach, in which global simulations are downscaled with an Arctic MM5 regional model which provides temperature and precipitation inputs to a glacier mass-balance model. The mass balances of Hubbard and Bering Glaciers, south-central Alaska, USA, are simulated for October 1994–September 2004. The comparisons of the mass-balance simulations using dynamically-downscaled vs observed temperature and precipitation data are in reasonably good agreement, when calibration is used to minimize systematic biases in the MM5 downscalings. The responses of the Hubbard (a large tidewater glacier) and Bering (a large surge-type glacier) mass balances to the future climate scenario CCSM3 A1B, a ‘middle-of-the-road’ future climate in which fossil and non-fossil fuels are assumed to be used in balance, are also investigated for the period October 2010–September 2018. Hubbard and Bering Glaciers are projected to have increased accumulation, particularly on the upper glaciers, and greater ablation, particularly on the lower glaciers. The annual net balance for the entire Bering Glacier is projected to be significantly more negative, on average (–2.0ma–1w.e., compared to –1.3ma–1w.e. during the hindcast), and for the entire Hubbard Glacier somewhat less positive (0.3ma–1w.e. compared to 0.4 ma–1w.e. during the hindcast). The Hubbard Glacier mass balances include an estimated iceberg calving flux of 6.5 km3 a–1, which is assumed to remain constant.

We explore spatial aliasing of non-Gaussian distributions of sea-ice thickness. Using a heuristic model and >1000 measurements, we show how different instrument footprint sizes and shapes can cluster thickness distributions into artificial modes, thereby distorting frequency distribution, making it difficult to compare and communicate information across spatial scales. This problem has not been dealt with systematically in sea ice until now, largely because it appears to incur no significant change in integrated thickness which often serves as a volume proxy. Concomitantly, demands are increasing for thickness distribution as a resource for modeling, monitoring and forecasting air–sea fluxes and growing human infrastructure needs in a changing polar environment. New demands include the characterization of uncertainties both regionally and seasonally for spaceborne, airborne, in situ and underwater measurements. To serve these growing needs, we quantify the impact of spatial aliasing by computing resolution error (Er) over a range of horizontal scales (x) from 5 to 500 m. Results are summarized through a power law (Er= bxm) with distinct exponents (m) from 0.3 to 0.5 using example mathematical functions including Gaussian, inverse linear and running mean filters. Recommendations and visualizations are provided to encourage discussion, new data acquisitions, analysis methods and metadata formats.

Self-organizing maps (SOMs) provide a powerful, non-linear technique to optimally summarize a complex geophysical dataset using a user-selected number of ‘icons’ or SOM states, allowing rapid identification of preferred patterns, predictability of transitions, rates of transitions, and hysteresis in cycles. The use of SOMs is demonstrated here through application to a 24 year dataset (1973–96) of monthly Antarctic sea-ice edge positions. Variability in sea-ice extent, concentration and other physical characteristics is an important component of the Earth’s dynamic climate system, particularly in the Southern Hemisphere where annual changes in sea-ice extent (temporarily) double the size of the Antarctic cryosphere. SOM-based patterns concisely capture the spatial and temporal variability in these data, including the annual progression of expansion and retreat, a general eastward propagation of anomalies during the winter, and sub-annual variability in the rate of change in extent at different times of the year (e.g. retreat in January is faster than in November). There is also often a general seasonal hysteresis, i.e. monthly anomalies during cooling follow a different spatial path than during warming.

An experimental investigation was performed on the kinetic friction coefficient of laboratory-grown, columnar saline ice sliding against itself. Tests were performed on a dual-opposing load apparatus specially manufactured for attachment to an MTS testing system. The mean kinetic friction coefficient, μ, was measured for sliding velocities from 10−6 to 5 × 10−2 m s−1 at temperatures from —3° to —40°C under a contact pressure of about 20 kPa. The ice specimens were oriented with grain columns perpendicular to the sliding interface. At -3°C and at —10°C, three distinct regions were observed: from 10−6 to about 10−5ms−1, μwas nearly constant at 0.5; at velocities from 10−5 to 10−3 m s−1, μ began to drop rapidly to about 0.1; and, above 10−3 m s−1, μ began to level off at ~0.05. The velocity at which μ began to decline increased with decreasing temperature. At temperatures below —10°C, μ increased from ~0.5 at v =10−6ms−1 to a peak value of ~0.7 near a velocity of 5 × 10−5ms−1 and then fell rapidly to about 0.1 at 10−2ms−1. In general, μ increased with decreasing temperature and sliding velocity.

The processes involved in the evolution of vertical profiles of Mg2+, Ca2+ and microparticle concentrations, as well as their seasonal variation in surface snow, were studied by weekly sampling from September 2003 to September 2004 of a snow pit on Ürümqi glacier No. 1, eastern Tien Shan, China. The development of the microparticle and Mg2+ and Ca2+ stratigraphy in the snow pit is closely related to the physical development of the snow–firn pack. The sampling site is located at 4130 ma.s.l. in the percolation zone of the glacier, and in addition to the effects of sublimation and wind erosion, melting plays a crucial role in both the physical and chemical evolution processes. During the winter, soluble aerosol concentrations in the surface layers are altered slightly by sublimation and wind erosion, and the concentrations are further modified as the wet season begins in late April. In contrast, soluble aerosol stratigraphy in the deeper layers remains relatively unchanged through the winter. In early summer, as melting occurs in the upper part of the snow–firn pack, meltwater carries chemical species to different depths in the underlying snow–firn layers, such that at the end of the ablation season, all of the surface cations might be leached out from the upper layers. In addition, the possible source of calcium and magnesium is discussed in this paper.

Ice algae are a key component in polar marine food webs and have an active role in large-scale biogeochemical cycles. They remain extremely under-sampled due to the coarse nature of traditional point sampling methods compounded by the general logistical limitations of surveying in polar regions. This study provides a first assessment of hyperspectral imaging as an under-ice remote-sensing method to capture sea-ice algae biomass spatial variability at the ice/water interface. Ice-algal cultures were inoculated in a unique inverted sea-ice simulation tank at increasing concentrations over designated cylinder enclosures and sparsely across the ice/water interface. Hyperspectral images of the sea ice were acquired with a pushbroom sensor attaining 0.9 mm square pixel spatial resolution for three different spectral resolutions (1.7, 3.4, 6.7 nm). Image analysis revealed biomass distribution matching the inoculated chlorophyll a concentrations within each cylinder. While spectral resolutions >6 nm hindered biomass differentiation, 1.7 and 3.4 nm were able to resolve spatial variation in ice algal biomass implying a coherent sensor selection. The inverted ice tank provided a suitable sea-ice analogue platform for testing key parameters of the methodology. The results highlight the potential of hyperspectral imaging to capture sea-ice algal biomass variability at unprecedented scales in a non-invasive way.

All radar power interpretations require a correction for attenuative losses. Moreover, radar attenuation is a proxy for ice-column properties, such as temperature and chemistry. Prior studies use either paired thermodynamic and conductivity models or the radar data themselves to calculate attenuation, but there is no standard method to do so; and, before now, there has been no robust methodological comparison. Here, we develop a framework meant to guide the implementation of empirical attenuation methods based on survey design and regional glaciological conditions. We divide the methods into the three main groups: (1) those that infer attenuation from a single reflector across many traces; (2) those that infer attenuation from multiple reflectors within one trace; and (3) those that infer attenuation by contrasting the measured power from primary and secondary reflections. To assess our framework, we introduce a new ground-based radar survey from South Pole Lake, comparing selected empirical methods to the expected attenuation from a temperature- and chemistry-dependent Arrhenius model. Based on the small surveyed area, lack of a sufficient calibration surface and low reflector relief, the attenuation methods that use multiple reflectors are most suitable at South Pole Lake.

The high-salinity surface layer of young sea ice was subjected to field and laboratory experiments. Artificial pools, in which young ice was formed, were opened within a fast-ice sheet in the Saroma lagoon, Hokkaido, in February of 1983 and 1984. The salinity of 1 mm thick surface layer of the young ice was observed as high as 42.4‰, which exceeds the seawater salinity of 31‰. The surface salinity increased with rising surface temperature. When a load was placed on the fast ice near the pool, seeped brine of salinity 72.5‰ was observed on the surface of the young ice; and when the load was removed, the brine disappeared. Meanwhile, brine permeabilities, both upward and downward, were measured in the laboratory, Both permeabilities decreased logarithmically with lowering surface temperature. A remarkable anisotropy was observed: the upward permeability was greater than downward, and the ratio of upward to downward premeability increased with lowering surface temperature from 5 at -3 °C to 33 at -5°C. Upward and downward permeabilities in ms-1 were respectively 1x10-4 and 2x10-5 at -3°C, 2x10-5 and 6x10-7 at -5°C, and at -10°C upward permeability was 3x10-7.

Comparison of the distribution of seasonal variations in surface elevation derived from a firn-densification–elevation model with observed variations derived from ERS-1/-2 satellite radar altimetry shows close similarity in the patterns of the amplitude of the variations over the North Greenland ice sheet. The amplitudes of the seasonal variations decrease from west to east and from south to north, determined by the accumulation rate and the surface-temperature distribution pattern. Several methods of estimating the amplitude of the seasonal variation in the observations are compared, including the use of a three-frequency sinusoidal function derived from the modeled seasonal variation that is asymmetric. The resulting correlation coefficient between the observed amplitude, estimated with the three-frequency function, and the modeled amplitude is 0.66 and the slope is 0.7. Residual differences may be caused by interannual variability in accumulation and temperature and other approximations in the model.

Uniaxial compression tests were performed on samples of the Greenland Ice Gore Project (GRIP) deep ice core, both in the field and later in a cold-room laboratory, in order to understand the ice-flow behavior of large ice sheets. Experiments were conducted under conditions of constant strain rate (type A) and constant load (type B). Fifty-four uniaxial-compression test specimens from 1327-2922 m were selected. Each test specimen (25 mm x 25 mm x 90 mm) was prepared with its uniaxial stress axis inclined 45° from the core axis in order to examine the flow behavior of strong single-maximum ice-core samples with basal planes parallel to the horizontal plane of the ice sheet. The ice-flow enhancement factors show a gradual increase with depth down to approximately 2000 m. These results can be interpreted in terms of an increase in the fourth-order Schmid factor. Below 2000 m depth, the flow-enhancement factor increases to about 20-30 with a relatively high variability When the Schmid factor was > 0.46, the enhancement factor obtained was higher than expected from the .-axis concentrations measured. The higher values of flow-enhancement factor were obtained from specimens with a cloudy band structure. It was revealed that cloudy bands affect ice-deformation processes, but the details remain unclear.