A free-drift sea-ice model for advection is described which includes an interactive wind-driven ocean for closure. A reduced system of equations is solved economically by a simple iteration on the water stress. The performance of the model is examined through a sensitivity study considering ice thickness, Ekman-layer scaling, wind speed, and drag coefficients. A case study is also presented where the model is driven by measured winds and the resulting drift rate compared to measured ice-drift rate for a three-day period during March 1981 at about 80 km inside the boundary of the open pack ice in the Bering Sea.

The advective model is shown to be sensitive to certain assumptions. Increasing the scaling parameter A for the Ekman depth in the ocean model from 0.3 to 0.4 causes a 10 to 15% reduction in ice speed but only a slight decrease in rotation angle (α) with respect to the wind. Modeled α is strongly a function of ice thickness, while speed is not very sensitive to thickness. Ice speed is sensitive to assumptions about drag coefficients for the upper (CA) and lower (CW) surfaces of the ice. Specifying CA and the ratio of CA to CW are important to the calculations.

A three-dimensional perfectly plastic ice-sheet model, developed for determining the surface elevations and the flow pattern of an ice sheet with given bottom topography and ice-margin positions, is applied to the reconstruction of the glacial ice covers of Greenland and the Canadian Arctic islands. In the northern regions, two different reconstructions have been performed with ice margins along the present 600 and 200 m sea-depth contours, respectively. In central Greenland, the ice margin is considered to be at the outermost ice-margin deposits on the coastal shelf to the west, and at the present 200 m sea-depth contour to the east.

The main conclusions to be drawn from the reconstructions are: (1). The flow pattern of the glacial ice cover of Greenland shows a great resemblance to the present one, the central ice divide being displaced less than 50 km from its present position and being no more than 200 m higher than today. (2). The main ice divide of the ice sheet covering the Canadian Arctic islands (the Innuitian ice sheet) was located over the highlands of eastern Ellesmere Island with local domes positioned over the present ice caps, indicating that even the deep ice of Wisconsin age in these ice caps is of local origin. This is also the case for the Devon Island ice cap. (3). Even in the not very likely case of a rather extensive glacial ice cover in north-west Greenland, the ice-flow pattern upstream of the Camp Century deep drill site would not have changed radically compared to the present flow pattern. Thus it is concluded that even advanced ice margins in late-Wisconsin time could at most have resulted in an elevation of the deposition site of the late-Wisconsin ice at Camp Century 600 m higher than at present. The consequences of this conclusion are discussed.

A good correlation was found between the ablation of snow and degree day index (cumulative values of positive daily mean air temperature) during the summer of 1978 on the Yukikabe snow patch in the Daisetsu mountains, central Hokkaido. The volume change of the snow patch in the ablation season of any year can hence be estimated from air temperature using this relationship. Each of the heat-balance terms controlling the ablation is evaluated separately by using empirical equations and assumed values for meteorological parameters at the snow patch. Triangular diagrams are constructed in order to illustrate the relative contributions of sensible heat, latent heat, and net radiation, the main three heat sources. A higher contribution from sensible and latent heat is found for the snow patches of Japan than for many glaciers and ice caps elsewhere. This may be due to higher mid-summer air temperatures than in other glaciated parts of the world.

Radiocarbon dates on 16 pine stumps and trunks revealed the position of the alpine pine-forest limit for the biostratigraphic (i.e. pollen)-defined Atlantic period (8.2 to 5.0 ka BP) in the coastal mountains of south-west Norway. Based on the present topographic mean summer temperature gradients (mean temperature from June to September), and the difference in elevation between the alpine pine-forest limits of the present and the Atlantic periods, mean summer temperature differences between these periods in the forest limit zone are calculated as 0.7°C to the west and 1.0°C in the inner valley districts to the east of the watershed. On the eastern side the mean summer temperature difference becomes greater towards lower altitudes. The alpine pine-forest limits of both the present and the Atlantic periods are characterized by a marked west-to-east inclination. This indicates pronounced westerly winds during the Atlantic period, similar to those today, characterized by more cloudy weather, higher precipitation and less sunshine to the west than to the east of the watershed.

A number of processes in the ocean must be modeled properly in order to produce valid estimates of oceanic heat transport, sea-surface temperature, and sea-ice extent in climate studies. These include: wind-driven turbulent mixing and water transport in the surface layer, internal vertical mixing due to several small-scale mechanisms, horizontal and vertical exchanges by mesoscale eddies, mixing along isopycnals, large-scale transport by currents, deep convection in polar regions, and boundary exchanges with atmosphere, ice, and land. Techniques to model these processes are described. Prospects are given for parameterizing the effects of phenomena that cannot be resolved in climate studies, particularly mesoscale eddies. Past simulations of the ocean in climate studies are reviewed. A modeling strategy is outlined for an improved treatment of the ocean, consistent with the computational power soon to be available.

Simple shear tests on the bottom 17 m of basal ice from Camp Century, Greenland, were carried out in order to study the flow behavior near the bottom of an ice sheet and its implications for ice-sheet modeling. The ice core was recovered in 1966. Our experimental results show that the basal ice tested, (which contained alternating bands of dirty and clean ice) has the highest strain-rate ever reported for polycrystalline ice under simple shear. The enhancement factors obtained are interpreted in terms of fabric, ageing, and impurity. Horizontal velocity profiles are calculated using data reported previously for Camp Century and Dye 3 stations. Various depth-age relationships are compared with these data. The higher than expected shear strain-rates measured on samples of near-bottom ice from Camp Century may very well exist at other locations. If such high shear strain-rates are more prevalent than presently thought, they could have an important bearing on ages calculated by physical or mathematical models of ice sheets.

As in polar ice, 18O variations of precipitation are recorded in carbonate sediments formed in lakes (lake marl). We have analyzed many late-glacial profiles from Europe, There are strong l8O variations which coincide with well-known pollen zone boundaries and which indicate that abrupt, drastic climatic changes occurred in the late glacial period. These events are a major warming around 13 ka BP (pollen zone boundary Oldest Dryas/Bölling) and a marked cold phase between about 10.8 and 10 ka BP (Younger Dryas). Comparison of the δ18O records of European lake sediments and of Greenland ice cores reveal a striking similarity which indicates that climatic changes in the late glacial and early postglacial were parallel in Greenland and in Europe. First results from North American lake-sediment profiles do not exhibit similar δ18O variations. This pattern of climatic changes was probably caused by retreating and readvancing polar water in the high-latitude North Atlantic Ocean, as discussed by Ruddiman and McIntyre (1981).

The effect of the oceanic heat capacity on global temperature response to external forcing is studied by considering a CO2-induced warming. The heat flux into the ocean is calculated using a box-diffusion model, earlier established for simulating the global carbon cycle. The calculated transient warming due to CO2 produced by fossil fuel until 1980 is 0.20 to 0. 26K, or about 50% of the value for radiative equilibrium, corresponding to a delay of 16 to 24 years. Analyses of CO2 in ice cores suggest a lower pre-industrial concentration, e.g. 265 ppm, than previously assumed (c. 290 ppm). The transient temperature increase until 1980 calculated for this case is about twice that for a case which starts from 297 ppm. Comparison with observations since 1880 shows that a long-term warming trend might be due to CO2, but the residual scatter unexplained by CO2 is still large. Finally the different thermal behaviour of oceans and continents is accounted for in a schematic way, considering heat exchange between them. The results indicate stronger temperature variations over the continents than over the sea.

Analyses of atmospheric air extracted from air bubbles in polar glacier ice provide a measure of past CO2 concentrations and their temporal variations. Earlier measurements have shown that the atmospheric CO2 concentration was significantly lower during the late Wisconsin stage and that the change from low to high CO2 concentration occurred in a relatively short time period.

We measured the CO2 concentration of ice samples from the deep ice core recently drilled at Dye 3 (south Greenland). The core section investigated represents ice from about 40 to 30 ka BP. The air extracted from the ice samples shows large variations of CO2 concentration which are correlated with δ18O values of the ice samples. A probable explanation of the results is that corresponding changes of the atmospheric CO2 concentration occurred during that glacial period. These changes could have had a significant influence on the climate. Possible reasons for the variations of the atmospheric CO2 concentration are discussed.

Two methods were used to locate the grounding line on the Rutford Ice Stream. The first method determined where the glacier was floating in hydrostatic equilibrium, while the second method measured the flexing close to the grounding line due to ocean tides. The ratio of surface elevation to ice thickness of the glacier goes through the hydrostatic equilibrium value 1 to 2 km downstream of where tidal flexing was recorded. This behaviour can be explained if the upward pressure of the sea at the base of the ice is augmented by a vertical shear-stress gradient within the glacier to overcome its weight. Simple el as tie-beam theory matches the flexure profile data if a modified elastic modulus and effective thickness are used. Tiltmeters can he used to monitor the position of the grounding line if the geometry of the flexing region can be defined.

A dynamic-thermodynamic sea-ice model (Hibler 1979) is used to simulate northern hemisphere sea ice for a 20-year period, 1961 to 1980. The model is driven by daily atmospheric grids of sea-level pressure (geo-strophic wind) and by temperatures derived from the Russian surface temperature data set. Among the modifications to earlier formulations are the inclusion of snow cover and a multilevel ice-thickness distribution in the thermodynamic computations.

The time series of the simulated anomalies show relatively large amounts of ice during the early 1960s and middle 1970s, and relatively small amounts during the late 1960s and early 1970s. The fluctuations of ice mass, both in the entire domain and in individual regions, are more persistent than are the fluctuations of ice-covered area. The ice dynamics tend to introduce more high-frequency variability into the regional (and total) amounts of ice mass. The simulated annual ice export from the Arctic Basin into the East Greenland Sea varies interannually by factors of 3 to 4.

Very small (ppm) amounts of soil dust in snow can significantly reduce snow albedo and thereby affect the snow-surface energy budget. Ice cores from Greenland show enhanced dust concentrations in ice from the last glacial maximum, in amounts capable of causing measurable effects on snow albedo. This enhanced dust is probably due in part to the expanded desert areas at that time.

Volcanic ash layers visible in the Byrd station core reduced the snow albedo in West Antarctica when they were on the surface. The ash is unlikely to have had a long-term effect on albedo because of the episodic nature of volcanic eruptions.

Very large amounts of dust on snow can inhibit snow-melt by insulating the snow. A debris cover probably slowed the melting of parts of the North American ice sheet during its most recent decay phase.

Snow in the Arctic Ocean is presently suffering large-scale contamination by carbon soot from anthropogenic sources. Preliminary estimates indicate that soot concentrations in Arctic snow are sufficient to reduce snow albedo measurably.

This paper presents a review of the ice-sheet mechanics behind the modeling work of Pleistocene ice sheets driven by astronomical (Milankovitch) forcing.

Detailed studies of the last 20 km of the flow-line leading to the core hole at Dye 3 Greenland, provide a description of ice flow over and around basal hills. The surface pattern is very simple. Velocity vectors are nearly parallel to one another and the largest variations in velocity are speed changes along the direction of flow. The surface elevation is stepped and the speed is faster than average where the surface slope is steepest. These positions correspond to basal highs, and the surface velocity increases as expected, based on the decrease in ice thickness, which indicates that most of the ice thickness must vary in velocity as does surface ice. Further support for this comes from the form of an internal radio-reflecting layer, which, in general, has the same shape as the bed but with much reduced relief. The damping of the relief is the same both along and across the flowline, suggesting that lateral velocity fluctuations are not important and that flow around and between obstacles is not well developed at the surface or at depth. At two sites, however, the internal layer does not match the bed and at one of these there must be important third-dimensional flow at depth.

The repetitive synoptic ice data obtainable by satellite sensing provide a means of studying the time-dependent behavior of both sea ice and ice sheets on climatic time scales. Examples of sea-ice parameters which may be measured are extent, concentration, and multiyear fraction; and examples of ice-sheet/ice-shelf parameters are surface elevation, ice-front position, extent and duration of summer melting, and ice accumulation rates. Desired snow-cover parameters include extent and snow depth or water-equivalent depth. The unique ability of satellites to measure such ice parameters and the characteristics of the consequent data sets significantly influence the structure of ice models that can be successfully used with the data. Ice data sets recently acquired by satellite sensing are described. The past decade of sea-ice data provides a detailed description of the interannual variability of sea ice on a regional and seasonal basis. Because of the longer time scales involved in ice-sheet variations, a comparable record of ongoing ice-sheet variations has not yet been established, but important baseline data sets are being developed.