Antarctic ice-elevation maps are presented, based on 4 000 measurement points, analysed from traverses made by constant-density balloons. These balloons were launched during the Tropical Wind, Energy Conversion and Reference Level Experiment (TWERLE). The TWERLE data yielded daily maps of the elevation of the 150 mbar pressure surface in the southern hemisphere. These maps serve as a reference level from which the altimeter readings are subtracted to yield the ice elevation.

The accuracy of the ice-elevation data is estimated at better than 60 m. The accuracy of the contouring depends on the density of measurements in the particular area. Two contouring techniques were used. One is based on generating elevation data at fixed grid-points from arbitrarily distributed measurements, using Cressman's method. The second is a surface spline interpolation technique by Duchon.

The elevation data are of special value on the inland plateau areas. There, the typical distance between balloon measurements is sufficiently small, relative to the typical distance between contour lines at 100 m intervals. Furthermore, the inland areas are the least accessible to other mapping techniques.

A new deep drill has been developed within the framework of the European Programme for Ice Coring in Antarctica (EPICA). Several versions of the EPICA drill exist. The version used at Dome Concordia (75˚06'1" S, 123˚23'71" E) was operated with a new electronic control system developed by the Ente per le Nuove tecnologie, l’Energia e l’Ambiente (ENEA) Research Center in Brasimone, Italy. This electronic control system was used for the first time during the 1997/98 Antarctic summer season.

Due to the nonlinear nature of the ice interaction, sea-ice build-up against coasts and structures is a complex process. This build-up signifycantly affects mesoscale (10 to 100 km) ice motions over typical forecast time scales of several days. To examine the ramifications of assuming a non-linear ice interaction in ice forecast models, we have carried out a series of idealized simulations employing a viscous plastic sea-ice rheology (Hibler 1979). These simulations employ constant wind fields at a grid resolution of 18.5 km and allow the ice to build up and strengthen. With the plastic ice interaction the ice build-up is found to take place by means of a ridging front. Depending on the nature of the strength-thickness coupling, this build-up is accompanied by kinematic wave propagation effects. The nonlinear interaction can also result in fluctuating velocities in certain locations, even though the forcing is fixed. The build-up results are found to be consistent with the analytic solution of a one-dimensional rigid plastic model.

In the period from 1982 to 1984, the Silesian University Expeditions investigated the annual cycle of the Hans Glacier velocity and front fluctuations. They also analysed the factors influencing these processes (J.Jania, L.Kolondra, E. Bukowska-Jania 1983).

Surrounding the Hans Glacier frontal zone, permanent tripods had been installed (by cementation in monolithic rocks) for the photo-theodolite, to establish three stereo-photogrammetric bases. While two of the three bases were used for surveying the fluctuations of the glacier front along the entire width (ca. 2.5 km), the remaining one was installed to record the velocity of the glacier by the time-parallactic method. Photogrammetric pictures were reiterated approximately every 10 days during August in 1982, 1983 and 1984. During the polar winters of 1982/83 and 1983/84, the oscillations of the glacier front were recorded at one base only. Pictures were taken once a month (also by moonlight). Using a Gornik-type seismograph, natural micro-tremors coming from the glacier were recorded continually, The seismograph works at the nearby Polish Polar Station, which operates a meteorological station.

The application of permanent metal tripods with an auto-centering disc made it possible to take successive pictures at the same external orientation of the camera on one hand, while, on the other hand, improving the convergent photographs (a similar approach was reported by U. Voigt 1966). On the glacier surface, ground points of control were signalled with a Maltese cross. The investigators made use of natural reference points, i.e. some characteristic features of the glacier surface. Maximum errors of the photogrammetric survey were mxy = ±0.3 m, mz = ±0.1 m.

The results of glacier tongue velocity measurements (ca 50 ma−1), as well as the results of measuring the fluctuations of the glacier front position enabled the rate of calving to be calculated. Thus, the calving velocity amounts to ca 100 ma−1 and the mass loss at the contact with sea water approaches ca 20% of the annual mass loss due to ablation.

The calving speed and the velocity of the glacier undergo variations in different periods of the year and the maxima of the processes do not overlap. While the glacier velocity reaches its maximum value at the beginning of the summer season (July), maximum calving speed is recorded in autumn (September-October). However, there may appear a shift in the time at which these maxima occur. It depends on the meteorological conditions and the thermal state of the sea in the given year. The effect produced by the two “antagonistic” glacial processes is the change in position of the glacier front in the sea. Its maximum and minimum extension appears by the end of July and in October, respectively. The amplitude of the Hans Glacier front fluctuations, measured for the period of August, 1982 to August, 1983, amounted to 60 ma−1 on the centre line. The results of photogrammetric surveying by C. Lipert, from 1957 to 1959, have shown that the maximum changes in the extension of the glacier front amount to 250 ma−1. These fluctuations display regularities similar to those reported for the Columbia Glacier, Alaska by C.S. Brown, M.F. Meier and A. Post (1982).

Analyses of micro-tremors coming from the glacier involved their variability in scale and frequency from one day to the next and throughout the year, as well as photogrammetric survey. Attempts were also made to find the englacial source of those micro-tremors. Source location was attempted in the summer of 1980, by using three geophones situated in the frontal part of the glacier. Thus, the majority of the micro-tremors owe their origin to the zone situated at a distance of 200 to 300 m from the front line and not to the ice cliff, as had been expected (A. Cichowicz, personal communication). It is worth noting that there exists an overlap of the annual distribution of the frequency of micro-tremor occurrence with the curve of glacier velocity variations. This enabled the investigator to determine the glacier dynamics by tremor recording.

The amount of snow cover in the Lake Biwa catchment area has changed significantly from 1.9 × 1012kg in 1974 to 0.3 × 1012kg in 1979. The amount of snow cover was less in the warmer winter and when the amount of snow cover was less than 1012 kg, the average amount of snow cover in recent years, the lowest dissolved-oxygen concentration in the deep layer rapidly decreases. Climatic warming will significantly decrease the amount of snow cover in Lake Biwa catchment area and the dissolved-oxygen concentration in the deep layer of the lake, which may further enhance eutrophication of the lake.

Large scale melting phenomena like meltwater drainage channels and meltwater accumulation basins or frozen “lakes” have been surveyed on the land ice mass in Dronning Maud Land, Antarctica. These melting features were also detected in a Landsat Thematic Mapper (TM) image recorded on 12 February 1990. Image processing techniques such as principal component (PC) analysis, band ratioing and histogram-equalizing are carried out to emphasize the melting phenomena. Interestingly, a histogram-equalized single TM Band 5 image appears to be a good discriminator, leaving the melting phenomena as bright areas in the TM scene. The largest frozen “lake” is close to 1km wide while some of the drainage channels stretch more than 5 km. There is also a comparison of satellite-derived and in situ reflectance and temperature. The atmospheric corrected satellite data underestimate surface snow albedo by about 16% while TM measurements of surface temperature are 2° to 6°C lower than that measured on the ground.

Trends in climate affecting the West Antarctic ice sheet may be detected first in the Antarctic Peninsula region. Although the area contains the most comprehensive weather records for any part of Antarctica, reliable snow-accumulation data are lacking.

Mainly as a result of the large snow-accumulation rate in the region (typically in the range 4.0–10.0 kg m−2 a−1), stratigraphie evidence of climate derived from ice cores can be resolved in much greater detail than is possible over most of the continent. Ice cores have been drilled at two sites, representing the extremes of climate type encountered in the region. A 133 m core has been obtained from Dolleman Island (70°35.2′S, 60°55.5′W) to represent the continental-type climate of the Weddell coast region, and an 87 m core has been obtained from the Palmer Land plateau (74°01′S, 70°38′W) to represent the more maritime regime of the west coast and central areas. Replicated cores were obtained at both sites in order to assess the contribution of local noise factors to the climatic signal preserved in the cores. Climatic trends during the period 1938–86 have been assessed on the basis of stable-isotope analysis of the top 47 m of the Palmer Land core and of the top 32 m of the Dolleman Island core.

A statistical analysis of derived profiles of mean annual δ18O and accumulation rate indicates that the local noise factors at these sites are sufficiently small that data averaged over periods as short as 5 years should reveal climatic shifts at the level of 0.2% and 5% respectively. These changes are much smaller than trends that have actually occurred during the past 50 years.

The most notable trend over the past 30 years is an increase of more than 30% in the snow-accumulation rate that has occurred in parallel with an overall temperature increase of 0.06°C/a during the same period. Increases of similar magnitude can be inferred from studies in East Antarctica, and may be related to a significant increase in precipitation rate that has been documented recently at mid-to high-latitude stations in the Northern Hemisphere. The finding may have relevance to studies of the possible consequences of a CO2-induced climate change. More extensive accumulation time series are now required from Antarctica, if satisfactory models of the long-term balance of the ice sheet are to be derived.

Bias-corrected suspended sediment rating equations are developed for the data from six separate ablation seasons in the Glacier de Tsidjiore Nouve basin. These equations differ significantly from one another. Suspended loads as predicted by equations used beyond their year of origin fall in the range of 34-278% of the true load. Mean absolute-percentage errors (MAPE) based on forecasts computed from each year’s model applied to the other 5 years of data range from 35 to 81%, with a mean of 52%. The equivalent MAPE for a lumped multi-year model is 38%, for an autocorrelation-correcting generalized least-squares (GLS) regression model, is 15%, and for a transfer function (TF) developed from a single 25 d period, only 5%. Simple structure-transfer functions appear to be more robust to temporal transfer than the usual type of rating equation, and offer the possibility of greatly improving upon the poor performance of OLS rating equations in providing estimates of suspended loads for periods beyond frames of origin.

The purpose of this Study is to propose a new constitutive law for pack ice, which is not only capable of Simulating the in-plane Shear and out-of-plane uniaxial compression, but also capable of avoiding overestimating divergence during Shear. The pack ice is treated as a two-dimensional granular plastic, obeying Coulomb’s friction law with a maximum principal Stress limit. During the out-of-plane uniaxial compression process the flow rule is normal, while during the in-plane Shear process the flow rule is coaxial with a linearly varying dilatancy angle describing the divergence. The Strength parameterizations are based on thickness and compactness of the pack ice; weakening or hardening can be achieved through advection and redistribution.

Topographic sketch mapping of the whole ablation area of the Khumbu Glacier, East Nepal Himalaya, is performed, using a simple, stereo-photogrammetric method. This map shows that the surface morphology can be classified into 11 morphologic characteristics. Depending on their distribution and combination, the ablation area can be divided into four morphologic areas. Detailed maps, on a scale of 1:1000–2500, of these four areas indicate that the distribution and combination of these 11 morphologic characteristics result from thickness of the debris cover, supraglacial streams and ponds, and glacier dynamics. The irregularity in the ablation area of the Khumbu Glacier can be considered to be a consequence of the mass balance between rate of ice charge from upstream and irregular distribution of ablation rate, depending on debris-cover characteristics in situ.

Formation mechanisms of bent icicles and thin spikes were studied by observations of natural icicles and by artificially growing icicles in a cold laboratory, and the results are discussed with reference to the general growth process of icicles. It is shown that icicles can be bent by the displacement of their roots, change in water supply rates, and wind. Spikes were observed to grow by the freezing of a water column captured at the centre of icicles after their length and diameter growth stopped.

Traditional models of Glaciol isostasy, derived predominantly for studying the response of the earth to the retreat of Northern Hemisphere ice sheets since the last ICC age, use earth models which assume constant lithospheric thickness. For Antarctica, where the continent below the ice sheet is two separate land masses differing in geological form, the assumption thai a uniform lithosphere can explain isostatic bchavioitr is questionable. Here, a method to calculate the glacio-isostatic adjustment of the continent with a laterally varying lithospheric thickness IS presented. The method is then used in a time-dependent ice-sheet model to model the isostatic response of the continent when the ice sheet passes through a Glaciol/interglacial transition. Various relations bet ween the “crustal thickness” beneath Antarctica, derived from seismic data, and the “lithospheric thickness” estimate used in glacio-isostatic calculations are assumed in a sensitivity study. Using the simplest relationships between crustal and lithospheric thickness, the greatest sensitivity of the ice sheet to the crustal structure of Antarctica is not in the interior but at coastal locations, particularly near the major ice shelves. During periods of ice-sheet advance the grounding-line migration of the shelves varies according to the depression of the earth peripheral to the ice sheet. The peripheral depression depends on the regional elasticity of the lithosphere which is controlled by the lithospheric thickness. Therefore the capacity for the ice sheet to advance varies with the regional thickness of the crust.

Dynamic ice-sheet models are used to assess the contribution of mass loss from the Greenland ice sheet to sea-level rise. Mass transfer from ice sheet to ocean is in a large part through outlet glaciers. Bed topography plays an important role in ice dynamics, since the acceleration from the slow-moving inland ice to an ice stream is in many cases caused by the existence of a subglacial trough or trough system. Problems are that most subglacial troughs are features of a scale not resolved in most ice-sheet models and that radar measurements of subglacial topography do not always reach the bottoms of narrow troughs. The trough-system algorithm introduced here employs mathematical morphology and algebraic topology to correctly represent subscale features in a topographic generalization, so the effects of troughs on ice flow are retained in ice-dynamic models. The algorithm is applied to derive a spatial elevation model of Greenland subglacial topography, integrating recently collected radar measurements (CReSIS data) of the Jakobshavn Isbræ, Helheim, Kangerdlussuaq and Petermann glacier regions. The resultant JakHelKanPet digital elevation model has been applied in dynamic ice-sheet modeling and sea-level-rise assessment.

A coupled ice-ocean model is used to study the sensitivity of the Beaufort Sea climatology to representation of sub-grid-scale eddies; to hypothetically not present and double Mackenzie River discharge; and to approximate climate warming specified through a surface air-temperature increase of 3° C. The eddy effect is considered in two ways: as eddy interaction with sea-floor topography yielding a driving force ("neptune" parameterization) and as eddy diffusion and viscosity. The model with neptune parameterization reproduces surface layer circulation, as well as the bathymetrically steered Beaufort Undercurrent, while the model with usual damping does not simulate the Beaufort Undercurrent. The absence of the strong boundary Beaufort Undercurrent affects the thermohaline structure of the Beaufort Sea which becomes less consistent with observational data. The increase of the Mackenzie River discharge causes more northward transport of sea ice, resulting in sea-ice thinning in Mackenzie Bay, while the absence of the Mackenzie River discharge induces southward sea-ice drift and sea-ice thickening in Mackenzie Bay. The sensitivity study of surface air-temperature warming shows a shrinkage of sea ice by 6% in area and 15% in volume, causing the freshening and warming of the surface ocean layer. The sensitivity studies of river discharge and surface air temperature use the neptune parameterization.

The equation describing the surface evolution of a large ice sheet is examined on the basis of a scale analysis applied to Antarctic conditions. Changes in the surface elevation are mainly driven by mass-balance fluctuations which approximately follow global atmospheric temperature variations. The essential spatial non-uniformity of the accumulation rate and the resultant difference between central and coastal regions in reaction time-scales are taken into account. The dynamic interaction of the time-lagging interior with the quasi-stationary margin is described. As a result, a simplified model is deduced to simulate the surface-elevation variations in the central parts of the Antarctic ice sheet caused by mass-balance perturbations corresponding to the main Milankovich cycles with the periods of 19–100 kyears. Special computational tests are performed to validate the model through intercomparison with the predictions obtained with a two-dimensional thermomechanical model. The sensitivity of the model to physical factors (represented by dimensionless tuning parameters) is discussed. Climatically controlled variations of the ice-sheet thickness in the vicinity of Vostok Station during the past 200 kyears are estimated.

Ice sheets consist of several disjoint regions, each with physically distinct behavior. Parts are cold; others are temperate and either partly or completely saturated. Varying dust content, impurities, debris content, etc., may affect the ice flow. Usually, these regions are separated by material or non-material surfaces or boundaries. We use mixture concepts, involving the balances of mass, momentum, energy and entropy. When applied to regular domain points these concepts and appropriate constitutive postulates yield the field equations for the evolution of the constituents. When formulated in terms of singular surfaces, boundary and transition conditions emerge. Our presentation takes the form of an extended abstract of work that is presently under consideration (Hutter and Engelhardt 1988; Hutter and Engelhardt, in preparation).

Blue Glacier was subject to 5 to 10% increase in thickness during the period 1958 to 1978. A surface velocity increase of 15 to 50% has accompanied this change in ice thickness. A non-uniform distribution of thickening over the glacier produced a general decrease in surface slope. The flow law parameters in a power-law relation for ice can be obtained from the resulting increase in velocity.

In 1957–58 the surface velocities and surface elevation along several transverse profiles were determined by others. During the period 1977 to 1980 similar measurements were made of surface velocity, elevation, and ice thickness along profiles spatially equivalent to those of the 1957–58 survey. Using these results a detailed study of the flow response has been made. Preliminary results indicate that the controlling surface slope a is one that is averaged over approximately ten times the glacier thickness and also indicate a relatively high stress exponent n = 4 to 5.

Blue Glacier flows in a complex channel and any detailed interpretation of the flow response requires the inclusion of this geometry in a model. For this reason multi-dimensional finite-element methods using a non-linear constitutive relation were employed. Various degrees of model sophistication allowed for the specification of different parameters governing the flow of the glacier. Variations in the flow law parameters and the curvature of the channel were found to have significant effects on the velocity and stress fields, as well as on the characteristics of the flow response to the thickness perturbation. Fully three-dimensional models of the flow of Blue Glacier and its response to the thickness increase were developed. Important effects such as those arising from longitudinal thickness and surface slope variations were included in three-dimensional modeling. The three-dimensional effects were seen to play an important role in flow of a rheologically non-linear fluid such as ice.

Intense deformation within an extensive englacial debris layer appears to be a dominant flow mechanism for glaciers on Mount St Helens, Washington. Observations of this phenomenon have been made on Shoestring Glacier and more recently on other glaciers on the volcano. Transverse velocity profiles on Shoestring Glacier have shown a marked plug flow behavior occurring between narrow marginal zones of crevassing. Both longitudinal and transverse velocity profiles on Shoestring Glacier were measured during the year before and the two years after the cataclysmic eruption of the volcano on 18 May 1980.

The geometry and regime of every glacier on Mount St Helens dramatically changed during the eruption. The entire accumulation area of Shoestring Glacier was removed by a volcanic landslide and blast, and the glacier surface was deeply incised through the action of pyroclastic flows and mudflows. Surface velocity measurements were supplemented by vertical profiles taken along the newly exposed ice cliffs. Large velocity gradients seen in the transverse, longitudinal and vertical profiles coincided with the mapped location of a recently exposed englacial debris layer. Further detailed measurements on the vertical profiles show that the intense shear occurs by deformation internal to the debris layer, A study of the properties of debris-laden ice was undertaken to determine the mechanism for localization of deformation to the debris layer and for the associated response of the glacier to a sudden change in stress.