As an inhomogeneous mixture of pure ice, brine, air and solid salts the physical properties of sea ice depend on its highly temperature-dependent microstructure. Understanding the microstructure and the way it responds to variations in temperature and salinity is crucial in developing an improved understanding of the interaction between sea ice and the environment. However, measurements monitoring the internal structure of sea ice are difficult to obtain without disturbing its natural state. We have recently developed an application of cross-borehole d.c. resistivity tomography to make in situ measurements that resolve the anisotropic resistivity structure of first-year sea ice. We present results from measurements made in 2008 off Barrow, Alaska, and in 2009 off Ross Island, Antarctica. the sea ice in these two regions forms in different environments: at Barrow, relatively quiescent conditions typically lead to a predominance of columnar ice, while more turbulent conditions and underwater ice formation in McMurdo Sound tend to produce a larger component of frazil or platelet ice. Interpretation of the resistivity measurements carried out in association with temperature and salinity data collected simultaneously allows both observation of the temporal evolution of the ice structure and, in the case of the Antarctic measurements, the identification of different ice types.

There are some reports on the measurement of the charge-to-mass ratio of blowing-snow particles, but there are few studies concerned with individual snow-particle charge. We measured the charge-to-mass ratios using snow particles selected according to size, and discussed individual charges. Experiments were conducted in a cryogenic wind tunnel. Charge-to-mass ratios measured in our experiment were all negative and their absolute values tended to increase with a decrease in particle diameter. Individual snow-particle charges were calculated from the average of particle diameter distributions. The charges were all approximated by the power function of diameter at each temperature. Assuming that the coefficient of these approximations is a function of air temperature, we could easily predict the individual snow-particle charge.

The spatial distribution and the local amount of snow in mountainous regions strongly depend on the spatial characteristics of snowfall, snow deposition and snow redistribution. Uniform altitudinal gradients can only represent a part of these influences but are without alternative for use in larger-scale models. How well altitudinal gradients represent the true snow distribution has not been assessed. We analyse altitudinal characteristics of snow stored in two high-alpine catchments in Switzerland. Peak winter snow depths were monitored using high-resolution airborne laser scanning technology. These snow depths were transferred to snow water equivalent by applying simple density estimations. From these data, altitudinal gradients were calculated for the total catchment areas and for selected subareas characterized by different accumulation patterns. These gradients were then compared with gradients resulting from automated snow depth measurements obtained from several snow stations on different height levels located in the catchments, and with estimations from climatological precipitation gradients. The analysis showed that neither precipitation gradients nor flat-field stations estimate catchment-wide snow amounts accurately. While the climatological gradient showed different trends for different areas and years, the snow stations tended to overestimate mean snow amounts.

We perform a linear-stability analysis of a system of ordinary differential equations describing the motion of dense flowing avalanches. The depth-averaged equations relate the mean translational velocity of the avalanche to the production of turbulent kinetic energy associated with the random motion of the snow granules. Flow friction is described by an extended Voellmy model, where both the Coulomb and turbulent resistance are related to the random kinetic energy. We identify two snow avalanche flow regimes: (1) an unstable frictional flow regime characterized by low random kinetic energy production and (2) a stable, collisional flow regime at the avalanche front. Flow-regime transitions are governed by the production of random kinetic energy, which is controlled by mass. The unstable regime is characterized by a saddle point which the avalanche encounters at the onset of motion and deposition. Depending on the release mass and/or mass growth, the avalanche head will either fluidize, leading to far-reaching avalanches, or starve and die out. At the tail of the avalanche this saddle point also controls the commencement of deposition and therefore the stopping behaviour of avalanches. We discuss the conditions for flow-regime transitions and the enhanced mobility of snow avalanches.

Although satellite data are useful for obtaining ice-thickness distribution for perennial sea ice or in stable thin-sea-ice areas, they are still largely an unresolved issue for the seasonal ice zone (SIZ). We address this problem using L-band synthetic aperture radar (SAR). In the SIZ, ice-thickness growth is closely related to deformation, so surface roughness is expected to correlate with ice thickness. L-band SAR, suitable for detecting such surface roughness, is a promising tool for obtaining thickness distribution. This idea was supported by an airborne polarimetric and interferometric SAR (Pi-SAR) validation. To extend this result to spaceborne L-band SAR with coarser resolution, we conducted in situ measurements of ice thickness and surface roughness in February 2008 in the southern Sea of Okhotsk with an icebreaker in coordination with the Advanced Land Observing Satellite (ALOS)/Phased Array-type L-band SAR (PALSAR) orbit. A helicopter-borne laser profiler was used to improve the estimation of surface roughness. It was found that backscatter coefficients (HH) correlated well with ice thickness (R = 0.86) and surface roughness (R = 0.70), which confirms the possibility of determining ice-thickness distribution in the SIZ. the interannual variation of PALSAR-derived ice-thickness distribution in the southern Sea of Okhotsk is also discussed.

Melt ponds are regularly observed on the surface of Arctic sea ice in late spring and summer. They strongly reduce the surface albedo and accelerate the decay of Actic sea ice. Until now, only a few studies have looked at the spatial extent of melt ponds on Arctic sea ice. Knowledge of the melt-pond distribution on the entire Arctic sea ice would provide a solid basis for the parameterization of melt ponds in existing sea-ice models. Due to the different spectral properties of snow, ice and water, a multispectral sensor such as Landsat 7 ETM+ is generally applicable for the analysis of distribution. an additional advantage of the ETM+ sensor is the very high spatial resolution (30 m). an algorithm based on a principal component analysis (PCA) of two spectral channels has been developed in order to determine the melt-pond fraction. PCA allows differentiation of melt ponds and other surface types such as snow, ice or water. Spectral bands 1 and 4 with central wavelengths at 480 and 770 nm, respectively, are used as they represent the differences in the spectral albedo of melt ponds. A Landsat 7 ETM+ scene from 19 July 2001 was analysed using PCA. the melt-pond fraction determined by the PCA method yields a different spatial distribution of the ponded areas from that developed by others. A MODIS subset from the same date and area is also analysed. the classification of MODIS data results in a higher melt-pond fraction than both Landsat classifications.

The seasonal duration of sea-ice cover in the region immediately north of Hokkaido in the Sea of Okhotsk was estimated from daily sea-ice concentration maps (1995–2009) produced by the Japanese Coast Guard. The results show that this period is generally well represented by the proxy record of sea-ice duration in the region of Utoro in the Shiretoko peninsula. Comparison of the duration of seasonal sea-ice cover along the coast of Hokkaido with meteorological conditions over the period of sea-ice mapping by the Japanese Coast Guard indicates that all sea ice will disappear from the coast of Hokkaido once air temperature is increased by 2–4˚C over present values by the effects of climate change.

Snow’s interaction with the environment is an important area in environmental science, particularly as climatic conditions change. We consider the local influence at the slope or smaller scale. Near-surface properties drive the dynamic interface with the atmosphere and with surrounding terrain. While accounting for topography, complex interactions involving energy and mass transfer at the snow surface are considered, using a computer simulation (RadTherm/RT). Digital elevation maps are used to numerically fabricate terrain features and vegetation, while applying appropriate thermal properties to specified terrain types. Conduction, convection radiation and phase change (for dry snow) are considered. Particularly relevant to this study are longwave infrared and shortwave radiation, which in the model account for shadowing, multiple reflective and emissive contributions. An example of a north-facing 30˚, snow-covered clearing bounded by trees is examined using measured meteorological conditions. Applying the same weather conditions, the model is used to examine the difference if the trees are assumed bare or covered with snow. Results indicate that, for the conditions considered, when trees are covered with snow, the open slope is cooler and the snow mass loss is less. Spatial variability across the slope is also noted. Differences are largely due to topographic radiation exchange.

The summer extent of the Arctic sea-ice cover has decreased in recent decades and there have been alterations in the timing and duration of the summer melt season. These changes in ice conditions have affected the partitioning of solar radiation in the Arctic atmosphere–ice–ocean system. the impact of sea-ice changes on solar partitioning is examined on a pan-Arctic scale using a 25 km × 25 km Equal-Area Scalable Earth Grid for the years 1979–2007. Daily values of incident solar irradiance are obtained from NCEP reanalysis products adjusted by ERA-40, and ice concentrations are determined from passive microwave satellite data. the albedo of the ice is parameterized by a five-stage process that includes dry snow, melting snow, melt pond formation, melt pond evolution, and freeze-up. the timing of these stages is governed by the onset dates of summer melt and fall freeze-up, which are determined from satellite observations. Trends of solar heat input to the ice were mixed, with increases due to longer melt seasons and decreases due to reduced ice concentration. Results indicate a general trend of increasing solar heat input to the Arctic ice–ocean system due to declines in albedo induced by decreases in ice concentration and longer melt seasons. the evolution of sea-ice albedo, and hence the total solar heating of the ice–ocean system, is more sensitive to the date of melt onset than the date of fall freeze-up. the largest increases in total annual solar heat input from 1979 to 2007, averaging as much as 4%a–1, occurred in the Chukchi Sea region. the contribution of solar heat to the ocean is increasing faster than the contribution to the ice due to the loss of sea ice.

Satellite radar altimetry provides data to monitor winter Arctic sea-ice thickness variability on interannual, basin-wide scales. When using this technique an assumption is made that the peak of the radar return originates from the snow/ice interface. This has been shown to be true in the laboratory for cold, dry snow as is the case on Arctic sea ice during winter. However, this assumption has not been tested in the field. We use data from an airborne normal-incidence Ku-band radar altimeter and in situ field measurements, collected during the CryoSat Validation Experiment (CryoVEx) Bay of Bothnia, 2006 and 2008 field campaigns, to determine the dominant scattering surface for Arctic snow-covered sea ice. In 2006, when the snow temperatures were close to freezing, the dominant scattering surface in 25% of the radar returns appeared closer to the snow/ice interface than the air/snow interface. However, in 2008, when temperatures were lower, the dominant scattering surface appeared closer to the snow/ice interface than the air/snow interface in 80% of the returns.

In Antarctica, blowing snow accounts for a major component of the surface mass balance near the coast. Measurements of precipitation and blowing snow are scarce, and therefore collected data would allow testing of numerical models of mass flux over this region. A present weather station (PWS), Biral VPF730, was set up on the coast at Cap Prud’homme station, 5 km from Dumont d’Urville (DDU), principally to quantify precipitation. Since we expected to be able to determine blowing-snow fluxes from the PWS data, we tested this device first on our experimental site, the Lac Blanc pass. An empirical calibration was made with a snow particle counter. Although the physics of the phenomenon was not well captured, the flux outputs were better than those from FlowCapts. The first data from Antarctica were reanalyzed. The new calibration seems to be accurate for estimating the high blowing-snow flux with an interrogation of the precipitation effects.

We present results from a laboratory tank study of ice growing from saline water in a wave field, focusing on the transition from a predominantly frazil/grease-ice cover to a pancake-ice cover. Combining surface temperature observations with direct and indirect determinations of ice salinity and solid fraction, we describe the evolution of frazil- and pancake-ice area fraction, salinity and solid-ice volume fraction over the course of 1 day. In the investigated stage of transition, frazil ice surrounding the pancakes was found to have rather constant properties: a surface temperature 0.4–0.6 K below the freezing point of the underlying sea water, a salinity of 24–26 g kg–1 and a solid volume fraction of 0.25–0.29. the average salinity of young pancake ice, estimated from heat and salt budgets, decreased from 18 to 15 g kg–1, while average solid volume fractions increased from 0.6 to 0.7. the transition from frazil to pancake is estimated to take place when solid fractions reach 0.37–0.40 and surface temperatures are 0.7–0.9 K below the freezing point. We find that, for proper modelling of the grease–pancake system, it is important to distinguish between a surface and volumetric fraction of pancakes and grease ice.

Following an International Geophysical Year project, we conducted meteorological observations during 2004–07 around the Suntar–Khayata range in eastern Siberia, where a strong temperature inversion exists throughout the winter. The temperature on the flat plain around Oymyakon (~700ma.s.l.) was ~20°C lower than that in a glaciated area located at ~2000ma.s.l. The inversion remained stable from October to April due to the Siberian high. Snowfall was limited to the beginning and end of winter. The stable conditions prevented atmospheric disturbances and inhibited snowfall during midwinter. From 1945 to 2003, glaciers in the Suntar–Khayata range retreated, with an area reduction of 19.3%. To assess this retreat, we estimated the response of the glaciers to climate change. According to US National Centers for Environmental Prediction (NCEP) data, the temperature in this region increased by ~1.9°C over 60 years. By calculating snow accumulation and ablation, the sensitivity of the equilibrium-line altitude (ELA) to the temperature shift was evaluated. We estimated snow precipitation based on precipitation at <0°C and ablation using the degree-day method. By these estimates, the ELA of Glacier No. 31, assumed 2350 m at present, could rise ~150m if temperature rises an average of 1°C. Furthermore, a 1.8°C temperature rise could cause the ELA to rise to 2600ma.s.l., removing the accumulation zone. With no accumulation zone, the glacier body would decrease, roughly halving in volume after ~400 years.

We conducted a pit study in July 2009 at the NEEM (North Greenland Eemian Ice Drilling) deep ice-coring site in northwest Greenland. To examine the seasonal variations of snow chemistry and characteristics of the drill site, we collected snow/firn samples from the wall of a 2 m deep pit at intervals of 0.03 m and analyzed them for electric conductivity, pH, Cl–, NO3–, SO42–, CH3SO3– (MSA), Na+, K+, Mg2+, Ca2+ and stable isotopes of water (δ18O and δD). Pronounced seasonal variations in the stable isotopes of water were observed, which indicated that the snow had accumulated regularly during the past 4 years. Concentrations of Na+, Cl– and Mg2+, which largely originate from sea salt, peaked in winter to early spring, while Ca2+, which mainly originates from mineral dust, peaked in late winter to spring, slightly later than Na+, Cl– and Mg2+. Concentrations of NO3– showed double peaks, one in summer and the other in winter to spring, whereas those of SO42– peaked in winter to spring. The winter-to-spring concentrations of NO3– and SO42– seem to have been strongly influenced by anthropogenic inputs. Concentrations of MSA showed double peaks, one in spring and the other in late summer to autumn. Our study confirms that the NEEM deep ice core can be absolutely dated to a certain depth by counting annual layers, using the seasonal variations of stable isotopes of water and those of ions. We calculated the annual surface mass balance for the years 2006–08. The mean annual balance was 176 mm w.e., and the balances for winter-to-summer and summer-to-winter halves of the year were 98 and 78 mm, respectively. Snow deposition during the winter-to-summer half of the year was greater than that during the summer-to-winter half by 10–20mm for all three years covered by this study.

Springtime measurements of CO2 exchange over seasonal sea ice in the Canadian Arctic Archipelago using eddy covariance show that CO2 was generally released to the atmosphere during the cold (ice surface temperatures less than about –6˚C) early part of the season, but was absorbed from the atmosphere as warming advanced. Hourly maximum efflux and uptake rates approached 1.0 and –3.0 μmol m–2 s–1, respectively. These CO2 flux rates are far greater than previously reported over sea ice and are comparable in magnitude to exchanges observed within other systems (terrestrial and marine). Uptake generally occurred for wind speeds in excess of 6 m s–1 and corresponded to local maxima in temperature at the snow–ice interface and net radiation. Efflux, on the other hand, occurred under weaker wind speeds and periods of local minima in temperature and net radiation. the wind speeds associated with uptake are above a critical threshold for drifting and blowing snow, suggesting that ventilation of the snowpack and turbulent exchange with the brine-wetted grains are an important part of the process. Both the uptake and release fluxes may be at least partially driven by the temperature sensitivity of the carbonate system speciation in the brine-wetted snow base and upper sea ice. the period of maximum springtime CO2 uptake occurred as the sea-ice permeability increased, passing a critical threshold allowing vertical brine movement throughout the sea-ice sheet. At this point, atmospheric CO2 would have been available to the under-ice sea-water carbonate system, with ramifications for carbon cycling in sea-ice-dominated polar waters.

Three ice type regimes at Ice Station Belgica (ISB), during the 2007 International Polar Year SIMBA (Sea Ice Mass Balance in Antarctica) expedition, were characterized and assessed for elevation, snow depth, ice freeboard and thickness. Analyses of the probability distribution functions showed great potential for satellite-based altimetry for estimating ice thickness. In question is the required altimeter sampling density for reasonably accurate estimation of snow surface elevation given inherent spatial averaging. This study assesses an effort to determine the number of laser altimeter ‘hits’ of the ISB floe, as a representative Antarctic floe of mixed first- and multi-year ice types, for the purpose of statistically recreating the in situ-determined ice-thickness and snow depth distribution based on the fractional coverage of each ice type. Estimates of the fractional coverage and spatial distribution of the ice types, referred to as ice ‘towns’, for the 5 km2 floe were assessed by in situ mapping and photo-visual documentation. Simulated ICESat altimeter tracks, with spot size ~70m and spacing ~170 m, sampled the floe’s towns, generating a buoyancy-derived ice thickness distribution. 115 altimeter hits were required to statistically recreate the regional thickness mean and distribution for a three-town assemblage of mixed first- and multi-year ice, and 85 hits for a two-town assemblage of first-year ice only: equivalent to 19.5 and 14.5 km respectively of continuous altimeter track over a floe region of similar structure. Results have significant implications toward model development of sea-ice sampling performance of the ICESat laser altimeter record as well as maximizing sampling characteristics of satellite/airborne laser and radar altimetry missions for sea-ice thickness.

Wet-snow avalanches are difficult to forecast, as the change from stable to unstable snow conditions occurs rapidly in a wet snowpack, often in response to water production and movement. Snow stratigraphy plays a vital role in determining flux behaviour. Capillary barriers or melt–freeze crusts can impede and divert water horizontally over large areas and thus may act as a failure layer for wet-snow avalanches. We present a comparison of measured and modelled liquid water content, θw, and snow stratigraphy during periods of wet-snow instabilities. Special attention is given to the reproducibility of capillary barriers, ponding of water on melt–freeze crusts and the timing of first wetting and of water arrival at the bottom of the snowpack, because these factors are believed to play a major role in the formation of wet-snow avalanches. In situ measurements were performed in the vicinity of automatic weather stations or close to recent wet-snow avalanches in order to compare them with model results. The simulations are based on two different water flux models incorporated within the 1-D snow-cover model SNOWPACK. The comparison of the two model runs with observed θw and stratigraphy revealed that both water-transport models reproduced the ponding of water on melt–freeze crusts. However, in both models melt–freeze crusts were transformed to normal melt forms earlier than observed in nature, so still existing ponding was not captured by the models. Only one of the models was able to reproduce capillary barriers in agreement with observations. The time of the first wetting at the surface was well predicted, but the simulated arrival time of the wetting front at the bottom of the snowpack differed between the simulations; it was either too early or too late compared with the observation.

During sea-ice formation, the dissolved constituents of water are rejected from ice crystals. the initial fractionation of dissolved organic matter (DOM) in the Baltic Sea ice was studied through two freezing experiments and by sampling natural sea ice. DOM was characterized by the spectral absorption of chromophoric DOM (CDOM) and the parallel factor analysis of DOM fluorescence. Molecular weight measurements of DOM were applied to assess changes in the molecular size distribution of DOM in new sea ice relative to parent sea water. Both in the newly formed artificial and the natural sea ice, CDOM was enriched by 34–39% relative to salinity. the same three identified DOM fluorophores were present both in sea water and ice but enriched by 15–54% relative to salinity in ice. After the incorporation of DOM into ice, the ageing of ice decreased the spectral slope coefficient and the molecular weight of DOM. This study shows that physical processes during freezing of brackish water enrich chromo- and fluorophoric DOM relative to salts in sea ice.

Meteorological data from mountainous areas of Japan have been collected by the National Research Institute for Earth Science and Disaster Prevention (NIED) for almost 20 years. The collected long-period data indicate that neither a notable increase in mean winter temperature nor a reduction in snow depth has occurred in these areas. The maximum snow depth, SDmax, and maximum snow water equivalent, SWEmax, show similar fluctuation trends, although with large year-to-year variations in value and a larger fluctuation range for SWEmax than for SDmax. This result suggests that monitoring of only SDmax in mountainous areas is not sufficient for understanding the quantitative fluctuation of water resources originating from snow. The SDmax fluctuation trends in mountainous areas sometimes differ from those in flatland areas because mountain SDmax depends more on winter precipitation than on mean winter air temperature, whereas the opposite is true for flatlands. In addition, the dependence ratio of SDmax on fluctuations in winter precipitation changes with altitude because the distributions of precipitation with air temperature change with altitude.