We present a new set of values for the spectral extinction coefficients, Kλ, for the interior of first-year (FY) and multi-year (MY) Arctic sea ice during the summer melt season measured during SHEBA (Surface Heat Budget of the Arctic Ocean program) and at Barrow, Alaska, USA. Results for FY ice are consistent with previously reported values, and differences can be understood in terms of variations in the concentration of biological and suspended particulate material. The values for the interior of MY ice are lower than previously reported for both bare and ponded ice. For bare MY ice the new Kλ values predict a substantial increase in the solar radiation transmitted through the ice into the upper mixed layer. Ponded MY ice is only slightly more transparent than previously reported, and FY ice values are generally consistent with previously reported values. Assuming an asymmetry parameter of 0.94, the extinction coefficients are consistent with a volume-scattering coefficient of 77 m–1 that is constant from 400 to at least 720 nm.

The snow/ice albedo was studied during a 4 week field experiment over first-year sea ice in the Gulf of Bothnia, Baltic Sea, in spring 2004. Observations were made on radiative fluxes, cloud cover, wind, air temperature and humidity, as well as snow/ice temperature, thickness, density and grain size. The albedo variation during the observation period was large: the daily mean albedo ranged from 0.79 over a new snow cover to 0.30 over bare, melting ice. The evolution of the albedo was related to the surface properties, but existing parameterizations based on Arctic data did not explain the observations well. The snow thickness was found to be the most critical factor affecting the albedo. A new parameterization was derived for the albedo dependence on snow thickness, to be applied over the Baltic Sea in spring, when periods of melting and freezing alternate but the ice is still relatively thick (about 0.6 m). The diurnal cycle of solar radiation was large, and the snow/ice metamorphism due to the melting during daylight and refreezing during the night caused a diurnal albedo cycle with a maximum in the early morning and a minimum in the afternoon, with an albedo difference up to 0.14 between the two.

Measurements of the incident solar radiation taken during the Antarctic Remote Ice Sensing Experiment (ARISE) aboard the R/V Aurora Australis in the Southern Ocean and springtime Antarctic ice pack are analyzed together with all-sky cloud imagery to determine the incident shortwave cloud radiative forcing at the surface. For most solar zenith angles (Z<82˚) in this dataset, the primary shortwave cloud effect is to induce cooling of the surface; as the sun approaches the horizon, however, the shortwave effects become negligible or even positive. The clear-sky atmospheric transmissivity over the length of the cruise is 0.91, a value comparable to those derived from measurements taken at various locations in the Arctic during daylight periods. Although the presence of clouds has a great effect on the surface heat budget and provides a negative shortwave feedback that may stabilize the polar atmosphere, the effect on the photosynthetically active radiation available to ice algae is relatively small in comparison to the effects of even small amounts of snow on sea ice.

The shortwave albedo is a major component in determining the surface energy balance and thus the evolution of the spring melt cycle. As the melt commences, the ice is partitioned into multiple surface types ranging from highly reflective white ice to absorptive blue ice. The reflectance from these surfaces shows significant spatial and temporal variability. Spectral albedo measurements were made at six different sites encompassing these two surface types, from 19 March to 3 May 2005, on 1.5 m thick landfast sea ice in southwestern Hudson Bay, Canada (58˚ N). Furthermore, the broadband albedo and the surface energy balance were continuously recorded at a nearby site during the 1 month period. Rapid changes in the albedo were found to relate to typical subarctic climate conditions, i.e. frequent incursions of southerly air, resulting snow and rain events and the generally high maximum solar insolation levels. Subsequently, diurnal variations in snow surface temperature were evident, often causing daytime melting and night-time refreezing resulting in the formation of ice lenses and superimposed ice. After rain events and extensive melting, the snowpack was transformed throughout into melt/freeze metamorphosed snow and superimposed ice. The integrated (350–1050 nm) albedo varied between 0.52 and 0.95 at the blue-ice sites, while it varied between 0.73 and 0.91 at white-ice sites. Variability on the order of ±10% in the white-ice broadband albedo resulted from the diurnal freeze–thaw cycle, but also synoptic weather events, such as snowfall and rain events, could rapidly change the surface conditions.

Annual and winter (December–April) sea-ice area and extent are calculated for the Greenland Sea (GS) and Barents Sea (BS) from daily ice concentrations obtained from space-borne microwave radiometry for 1979–2003. The ice extent decreases significantly, particularly during winter, by 65 000 km2 (decade)–1 in the GS and by 72 000 km2 (decade)–1 in the BS. Ice-extent fractions (of these total extents) occupied by ice of five different ice-concentration ranges are calculated and analyzed. Changes in these fractions are again significant and most pronounced during winter. In the GS, the fraction of close to very compact ice (65–95%) decreases by 17 000 km2 (decade)–1 and the fraction of very compact ice (>95%) increases by 29 000 km2 (decade)–1, corresponding to a loss of 19% and a gain of 58% relative to the 25 year mean, respectively. In the BS, the fraction of close to compact ice (65–85%) increases by 26 000km2 (decade)–1 and the fraction with compact to very compact ice (>85%) decreases by 66 000 km2 (decade)–1, corresponding to a gain of 30% and a loss of 67% relative to the 25 year mean, respectively. The changing surface wind pattern analyzed from ERA-40 data favours this increasing (decreasing) ice compactness in the GS (BS).

To understand the role of wind conditions on the summertime surface ocean system, the ocean, ice and atmosphere in the central Arctic Ocean were observed using two drifting buoys, one in 2002 under stormy conditions, one in 2003 under calm conditions. Although the ice concentration near the North Pole was the same in 2002 and 2003 during early summer, the heat used in bottom melting in 2003 was only about half of that in 2002. To obtain the total heat input into the upper ocean, heat used in lateral melting was additionally derived from a time series of ice concentration in 2002. Assuming the same heat input into the upper ocean, the heat used in lateral and bottom melting was estimated and compared between the years. It is thought that the warm fresh water embedded within the ice cover was mixed downward during the frequently stormy mid-summer of 2002, enhancing bottom melting. By contrast, the warm water in 2003 tended to be used for lateral melting due to the relatively calm conditions, suggesting that a continuously weak wind favours ice-cover decrease during summer. A simple calculation of the ice-cover evolution reveals that the difference in ice concentration during August between 2002 and 2003 reached 10%, which is consistent with the satellite-derived ice concentration.

Measurements of under-ice turbulence were performed using an acoustic three-dimensional current meter with an attached fast-repetition temperature–conductivity sensor at two coastal areas in the Baltic Sea during two winters. Observations covered both the ice-growth and spring-melt periods. The objective of these measurements was to obtain knowledge of under-ice turbulence and oceanic heat and salt fluxes to and from the ice in the coastal fast-ice region using eddy correlation techniques. The maximum observed daily average heat flux was 1 Wm–2, and the maximum for 10 min periods was an order of magnitude larger. Under-ice turbulence was much smaller than that recorded in the oceans and in coastal regions with tide. These results provide better knowledge of under-ice turbulence and heat-flux variations and are useful for the future development of a Baltic Sea ice salinity model.

The reflection, absorption and transmission of ultraviolet light by a sea-ice cover strongly impacts primary productivity, higher trophic components of the food web, and humans. Measurements of the incident irradiance at 305, 320, 340 and 380 nm and of the photosynthetically active radiation were made from April through September 1998 as part of the SHEBA (Surface Heat Budget of the Arctic Ocean program) field experiment in the Arctic Ocean. In addition, observations of snow depth and ice thickness were made at more than 100 sites encompassing a comprehensive range of conditions. The thickness observations were combined with a radiative transfer model to compute a time series of the ultraviolet light transmitted by the ice cover from April through September. Peak values of incident ultraviolet irradiance occurred in mid-June. Peak transmittance was later in the summer at the end of the melt season when the snow cover had completely melted, the ice had thinned and pond coverage was extensive. The fraction of the incident ultraviolet irradiance transmitted through the ice increased by several orders of magnitude as the melt season progressed. Ultraviolet transmittance was approximately a factor of ten greater for melt ponds than bare ice. Climate change has the potential to alter the amplitude and timing of the annual albedo cycle of sea ice. If the onset of melt occurs at increasingly earlier dates, ultraviolet transmittance will be significantly enhanced, with potentially deleterious biological impacts.

Snow depth on sea ice plays a critical role in the heat exchange between ocean and atmosphere because of its thermal insulation property. Furthermore, a heavy snow load on the relatively thin Southern Ocean sea-ice cover submerges the ice floes below sea level, causing snow-to-ice conversion. Snowfall is also an important freshwater source into the weakly stratified ocean. Snow-depth on sea-ice information can be used as an indirect measure of solid precipitation. Satellite passive microwave data are used to investigate the interannual and regional variability of the snow cover on sea ice. In this study we make use of 12 years (1992–2003) of Special Sensor Microwave/Imager (SSM/I) radiances to calculate average monthly snow depth on the Antarctic sea-ice cover. For the Antarctic sea-ice region as a whole, we find that September snow depth and sea-ice area are negatively correlated, which is not the case for individual regions. An analysis of the snow depth around Antarctica was undertaken. The results show an overall increase in snow depth for each of the five Antarctic sectors and the region as a whole, but only the Indian Ocean sector and the entire Southern Ocean show a statistically significant increase. There is a partial eastward propagation of maximum snow depths, which may be related to the Antarctic Circumpolar Wave. The overall trend and the variability of regional snow-depth distributions are also in agreement with cyclone density.

In November 2004 we used an acoustic Doppler velocimeter (ADV) to measure the three components of velocity in Langmuir circulation (LC) cells when frazil was forming in a lake. LC circulation was indicated by windrows of slush ice on the water surface. We also collected samples of lake water, ice interstitial water and ice to determine their sediment concentrations. The ADV record showed rotating currents in the cross-wind/vertical plane indicative of LC. Downward-directed velocities were large enough to entrain frazil into the water column for 27% of the 97 min observation time at the ADV location. Sediment concentrations in the ice and interstitial water samples were greater than concentrations in water-column samples collected in upwelling zones of LC cells. We conclude that suspension freezing was forming particle-laden ice as frazil trapped below LC convergence zones grew in downwelled, supercooled water. The rotating roll vortices of Langmuir circulation play an important role in forming sediment-laden ice.

Based on biophysical ice-core data collected in the landfast ice off Barrow, Alaska, USA, in 2002 and 2003, a one-dimensional ice–ocean ecosystem model was developed to determine the factors controlling the bottom-ice algal community. The data and model results revealed a three-stage ice-algal bloom: (1) onset and early slow growth stage before mid-March, when growth is limited by light; (2) fast growth stage with increased light and sufficient nutrients; and (3) decline stage after late May as ice algae are flushed out of the ice bottom. Stages 2 and 3 are either separated by a transition period as in 2002 or directly connected by ice melting as in 2003, when in situ light and nutrient enrichment experiments showed only light limitations. The modeled net primary production of ice algae (NPPAi) from March to June is 1.2 and 1.7 g Cm–2 for 2002 and 2003, respectively, within the range of previous observations. Model sensitivity studies found that overall NPPAi increased almost proportionally to the initial nutrient concentrations in the water column. A phytoplankton bloom (if it occurs as in 2002) would compete with ice algae for nutrients and lead to reduced NPPAi. About 45% of the NPPAi was exported to the shallow benthos.

A polynya appears regularly in Storfjorden on the east side of the Svalbard archipelago. It is mainly forced by offshore winds and contributes around 10% of the brine water produced on Arctic shelves. We have applied a regional ocean model (ROMS), including a sea-ice model, on a fine grid (2 km) to simulate a full year of sea-ice growth and decay starting on 1 August 1999. This allows us to reproduce some key processes of the polynya opening and closing events during January–April 2000. The polynya remains open as long as the offshore winds exist, and reaches a width along the direction of the wind of 10–20 km. We suggest using a mean sea-ice thickness of <0.3m as the polynya criterion, as our simulations show varying strength in the horizontal gradients in sea-ice concentration and thickness. Results show a general freeze-up during December and January, with a mean polynya area during February–April within the fjord of 33 by 50 km, being 13% of the total fjord area. Some model results including sea-ice cover and drift speed can be partially validated using satellite imagery and field data, but in general are new measurements from the polynya interior needed to further improve the modelling of solid- and grease-ice processes.

Horizontal variation of landfast sea-ice properties was studied in the Gulf of Bothnia, Baltic Sea, during March 2004. In order to estimate their variability among and within different spatial levels, 72 ice cores were sampled on five spatial scales (with spacings of 10 cm, 2.5 m, 25 m, 250m and 2.5 km) using a hierarchical sampling design. Entire cores were melted, and bulk-ice salinity, concentrations of chlorophyll a (Chl a), phaeophytin (Phaeo), dissolved nitrate plus nitrite (DIN) as well as dissolved organic carbon (DOC) and nitrogen (DON) were determined. All sampling sites were covered by a 5.5–23 cm thick layer of snow. Ice thicknesses of cores varied from 26 to 58 cm, with bulk-ice salinities ranging between 0.2 and 0.7 as is typical for Baltic Sea ice. Observed values for Chl a (range: 0.8–6.0 mg Chl a L–1; median: 2.9 mg Chl a L–1) and DOC (range: 37–397 μM; median: 95 μM) were comparable to values reported by previous sea-ice studies from the Baltic Sea. Analysis of variance among different spatial levels revealed significant differences on the 2.5km scale for ice thickness, DOC and Phaeo (with the latter two being positively correlated with ice thickness). For salinity and Chl a, the 250 m scale was found to be the largest scale where significant differences could be detected, while snow depth only varied significantly on the 25 m scale. Variability on the 2.5 m scale contributed significantly to the total variation for ice thickness, salinity, Chl a and DIN. In the case of DON, none of the investigated levels exhibited variation that was significantly different from the considerable amount of variation found between replicate cores. Results from a principal component analysis suggest that ice thickness is one of the main elements structuring the investigated ice habitat on a large scale, while snow depth, nutrients and salinity seem to be of secondary importance.

The reflection and transmission of ice-coupled waves under a sea-ice sheet is re-examined in this work. Recent theory can account for inhomogeneities in ice sheets such as the sails of pressure ridges and cracks or leads, but has only examined small numbers of features. Here the scattering coefficients can be obtained for an ice sheet containing a region of arbitrarily varying thickness that more closely approximates a real ice sheet and, consequently, models wave scattering more effectively. We explore the consequences of inverting ridge keels, placing the mass above the ice sheet as a technique by which the keel can be partially accounted for. It is found that the ridge sails themselves can reasonably be neglected. Upward-looking submarine sonar data are used to demonstrate the use of the model, noting that there are choices as to how to interpret such data. The model is tested to establish its sensitivity to these choices.

Safe exploration and transportation of natural energy resources in polar and subpolar seas such as the Sea of Okhotsk and the Arctic Ocean requires an understanding of the characteristics of ice-coupled wave propagation. Using a time-domain solution involving both the boundary element method and the finite-element method, a numerical procedure is developed to analyze the wave properties for arbitrary ice conditions. This is done by applying a distinction index to discrete nodes representing the dynamic boundary conditions. The numerical results agree well with experimental data for different floe lengths, thicknesses and elastic moduli, obtained by using model ice plates. The elastic deformation of the ice floe depends strongly on the flexural rigidity of individual ice plates.

We present a model for wave propagation in the marginal ice zone (MIZ). The model is based on a two-dimensional linear water wave floating-elastic-plate solution in the frequency domain using a matched eigenfunction expansion. The model is fully coherent so the results are dependent on exact geometry, and we show firstly that this dependence can be effectively removed by averaging over random floe lengths. We present a range of predictions from the model which show various curves of transmitted energy as a function of the floe thickness, number of floes and wave period. These curves show that it is possible with this simple model to make predictions about wave intensity with distance into the MIZ. The results predict that the MIZ acts as a low-pass filter and and that the filtering is strongly dependent on floe thickness and number but that it is independent of floe length.

Traditional methods of measuring the propagation of waves originating from ocean swell and other sources have relied on wire strain gauges, accelerometers or tiltmeters. All methods required constant attention to keep in range, while data recovery has demanded that the instrument site be revisited. In this paper, we describe the use of ultra-sensitive tiltmeters and novel re-zeroing techniques to autonomously gather wave data from both polar regions. A key feature of our deployments has been the use of the Iridium satellite communications system as a way of ensuring continuous data recovery and remote control of the instrumentation. Currently four instruments have been successfully reporting from the Arctic Ocean for over 18 months, with two further units deployed in 2005, one in the Weddell Sea, Antarctica, and one additional unit in the Arctic.

The kinetics of ice growth on the prismatic and basal planes is studied by molecular dynamics simulations. The time evolution of two systems has been investigated. In one a slab of ice is initially in contact with supercooled water, while in the second the ice is in contact with a supercooled salt solution. The simulations were done at a temperature below the eutectic temperature, and complete solidification is observed. The total freezing time is longer in the systems with ions than in the systems with pure water. The final state for the salt systems always shows the formation of ion clusters. For the ionic system growing on the prismatic plane, an intermediate metastable state is observed before total solidification. The duration of this metastable state depends on the ability of the system to get all the ions participating in cluster formation. The simulations enable understanding of the mechanisms for ice formation under different solution conditions.

We present new analytical results relating to the growth and evolution of sea ice. It is noteworthy that thin sea ice plays a central role in the surface heat and mass balance of the Arctic Ocean. In order to describe these balances, we analyze highly resolved temperature data taken through the air/sea/ice interface during the transition from an ice-free to an ice-covered Arctic Ocean surface. Our detailed analysis of the field data is based on the classical model of a mushy layer, which is modified in order to obtain analytical solutions in explicit form (so, for example, ice thickness and growth rate, temperature distributions, conductive and latent heat fluxes are determined). Furthermore, we find that the sea-ice growth is not simply a square-root function of time. It depends on the temperature variations in the atmosphere and lies between two square-root functions of time for the maximum and minimum temperatures found during observations. The theory under consideration is in good agreement with observations.