A rigid mast, equipped with five multi-turn potentiometers, was bolted to the bedrock floor of a frontal cavity at Glacier de Tsanfleuron, Switzerland. Each potentiometer was linked by a thin cord to an anchor emplaced in the adjacent ice wall, enabling basal ice motion to be measured continually over a 6 day period in the late summer. Results indicate that basal ice velocity increased quasi-linearly away from the glacier bed, rising from 10.66 to 11.82 mm d−1 between 25 and ∼ 265 mm above the ice– bed interface. Extrapolation of this gradient indicates that ∼ 10.55 mm d−1 may have been accommodated by pure slip between the ice and the glacier bed. Basal ice motion is temporally variable, generally being characterized by tens of minutes to hours of little or no motion interspersed with rapid-motion “events” lasting for between < 120 s (a single measurement interval) and 360 s. These motion events were responsible for speeds of up to > 400 mm d−1 over individual measurement intervals. The magnitude–frequency distribution of these events is consistent with a pattern of infrequent and large “slips” initiated at the ice–bed interface, that are manifested as more frequent and smaller motion events some tens of cm above that interface.

The freezing of sea water to the base of an ice shelf can give rise to large patches of accumulated ice, a phenomenon known as marine ice. In this study a numerical method is presented for calculating the thickness of the marine-ice layer using an ice- shelf-ocean model. The present-day modeling paradigm of ice-shelf–ocean interaction usually involves the fixed specification of the ice-shelf geometry while the ocean circulation in the cavity beneath the ice shelf evolves freely. This approach relies on several assumptions, such as steady-state ice-shelf thickness and ice-shelf flow fields, in order to make reasonable quantitative estimates of the thermodynamic exchange processes occurring at the ice-shelf base. This paper discusses the impact of these and other assumptions on the estimation of the thickness of the marine-ice layer. Model simulation results are presented for an idealized ice-shelf–ocean configuration as a demonstration of the feasibility of the numerical method. A sensitivity analysis is given so as to quantify the relative uncertainty in the marine-ice thickness that arises from uncertainties in the model input parameters, these being principally the ice-shelf flow field, the basal accumulation rate and the ice-shelf thickness field.

The transition from inland- to streaming-style ice flow near to and upstream from the onset to Ice Stream D, West Antarctica, is investigated using the force-balance technique. Basal drag provides the majority of the flow resistance over the study area but is substantially modified by non-local stress gradients. Lateral drag increases with distance downstream, balancing ∼50–100% of the driving stress at the onset. Longitudinal stress gradients (LSG) are also found to be significant, an observation that distinguishes ice flow in this region from the inland- and streaming-flow regimes that bound it, in which LSG are usually negligible. LSG decrease the spatial variability in basal drag and sliding speed and increase the area of the bed over which frictional melting occurs. Overall, LSG decrease the resistive influence of basal stress concentrations and increase the spatial uniformity of basal sliding. These observations suggest that streaming flow develops as an integrated response to the physical interaction between the ice and its bed over an extended region upstream from the onset, rather than being solely due to changes in basal characteristics at the onset. An implication is that non-steady-flow behavior upstream from the onset may ultimately propagate downstream and result in non-steady behavior at the onset.

A major discontinuity in the variation of δ18O (δD) with altitude in high mountains was first seen in data from Mount Logan, YukonTerritory, Canada (Holdsworth and others, 1991). The profile of δ vs altitude revealed three well-defined regions: (1) a lower, monotonic, fractionation sequence below ∼3 km; (2) a middle layer, typically 1–2 km thick, within which δ values are nearly constant or stepped with altitude, and (3) part of another fractionation sequence in the “quasi-geostrophic flow region” above ∼5.3 km. The middle region was inferred to be a “mixed layer”, combining moisture from regions (1) and (3). This type of structure is now seen to occur on other high-altitude mountains, including Cerro Aconcagua, Argentina, where observations reach almost 7 km. The new observations confirm the general occurrence of a multi-layered atmosphere during precipitation at high-altitude glacier sites. This structure is linked to synoptic-scale polar cyclones, where the middle layer is identified as being the signature of the warm-front zone. These results have implications for the common practice of using a specific, spatially derived, isotopic thermometer in the time domain for the paleoclimatic interpretation of high-altitude ice-core δ records.

Observations of the motion and basal conditions of Worthington Glacier, Alaska, U.S.A., during late-winter and spring melt seasons revealed no evidence of a relationship between water pressure and sliding velocity. Measurements included borehole water levels (used as a proxy for basal water pressure), surface velocity, englacial deformation, sliding velocity, and time-lapse videography of subglacial water flow and bed characteristics. The boreholes were spaced 10–15 m apart; six were instrumented in 1997, and five in 1998. In late winter, the water-pressure field showed spatially synchronous fluctuations with a diurnal cycle. The glacier’s motion was relatively slow and non-cyclic. In spring, the motion was characterized by rapid, diurnally varying sliding. The basal water pressure displayed no diurnal signal, but showed high-magnitude fluctuations and often strong gradients between holes. This transition in character of the basal water-pressure field may represent a seasonal evolution of the drainage system from linked cavities to a network of isolated patches and conduits. These changes occurred as the glacier was undergoing a seasonal-velocity peak. The apparent lack of correlation between sliding velocity and water pressure suggests that local-scale water pressure does not directly control sliding during late winter or early in the melt season.

Avertical numerical model has been developed that simulates tank experiments of sediment entrainment into sea ice. Physical processes considered were: turbulent vertical diffusion of heat, salt, sediment, frazil ice and their aggregates; differential growth of frazil-ice crystals; secondary nucleation of crystals; and aggregation between sediment and ice. The model approximated the real size distribution of frazil ice and sediment using five classes of each. Frazil crystals (25 μm to 1.5 cm) were modelled as discs with a constant thickness of their diameter. Each class had a constant rise velocity based on the density of ice and drag forces. Sediment grains (1–600 μm) were modelled as constant density spheres, with corresponding sinking velocities. The vertical diffusion was set constant for experiments based on calculated turbulent rms velocities and dissipation rates from current data. The balance between the rise/sinking velocities and the constant vertical diffusion is an important feature of the model. The efficiency of the modeled entrainment process was estimated through α, an aggregation factor. Values for α are in the range 〈0.0003, 0.1〉, but average values are often close to 0.01. Entrainment increases with increasing sediment concentration and turbulence of the water, and heat flux to the air.

During the 1999–2000 Italian Expedition, an airborne radar survey was performed along 12 transects across Lake Vostok, Antarctica, and its western and eastern margins. Ice thickness, subglacial elevation and the precise location of lake boundaries were determined. Radar data confirm the geometry derived from previous surveys, but with some slight differences. We measured a length of up to 260 km, a maximum width of 81 km and an area of roughly 14 000 km2. Along the major axis, from north to south, the ice thickness varies from 3800 to 4250 m, with a decreasing gradient. From west to east the ice thickness is fairly constant, except for two narrow strips located on the western and eastern margins, where it increases with high thickening rate. Over the lake the surface elevation increases from 3476 m a.s.l. (south) to 3525 (north), with a decreasing gradient, while the lake surface elevation decreases from −315 to −750 m a.s.l., with a decreasing gradient (absolute value). The ice-surface and lake-ceiling slopes suggest that the lake is in a state of hydrostatic equilibrium.

The interaction between ice sheets and the rest of the climate system at long time-scales is not well understood, and studies of the ice ages typically employ simplified parameterizations of the climate forcing on an ice sheet. It is important therefore to understand how an ice sheet responds to climate forcing, and whether the reduced approaches used in modeling studies are capable of providing robust and realistic answers. This work focuses on the accumulation distribution, and in particular considers what features of the accumulation pattern are necessary to model the steady-state response of an ice sheet. We examine the response of a model of the Greenland ice sheet to a variety of accumulation distributions, both observational datasets and simplified parameterizations. The predicted shape of the ice sheet is found to be quite insensitive to changes in the accumulation. The model only differs significantly from the observed ice sheet for a spatially uniform accumulation rate, and the most important factor for the successful simulation of the ice sheet’s shape is that the accumulation decreases with height according to the ability of the atmosphere to hold moisture. However, the internal ice dynamics strongly reflects the influence of the atmospheric circulation on the accumulation distribution.

Measurements of natural and artificial radioisotopes (32Si, 210Pb and 137Cs) and oxygen isotopes (δ18O) have been carried out on surface snow and ice, shallow snow pits and an ice core collected from Dokriani Bamak glacier, central Himalaya, to study the dynamics of glacier ice and short-term climatic changes. Based on the 32Si and 210Pb activities in the meltwaters, the age of the snout ice is 400 years and the flow rate of ice along the glacier length is ∼14 m a−1. The specific activity of 137Cs, corresponding to 1963 fallout, in the surface ice at the equilibrium line yields a flow rate of 32 m a−1, a factor of two higher than that derived for the snout ice. The depth variation of 137Cs concentration in a shallow ice core yields a mean accumulation rate of 0.43 m a−1 for the glacier ice over the past decade. The δ18O of snout ice (−13.4‰) is significantly depleted compared to the average value of −9.2‰ in the shallow ice core, indicating that cooler climatic conditions prevailed around AD 1600. Based on the oxygen isotopic ratios in the shallow pits, an “altitude effect” of 0.9‰ per 100 m in δ18O variation is documented for this glacier.

In this study we investigate the ability of the Regional Atmospheric Climate MOdel (RACMO) to represent temporal variations of Antarctic accumulation using weekly stake array measurements from Neumayer Station. The model uses European Centre for Medium-Range Weather Forecasts re-analysis data to force the atmospheric variables at the lateral boundaries of the model domain. Accumulation is defined as precipitation minus sublimation. Generally the model represents reasonably well the synoptic situations that lead to precipitation. The amounts of accumulation are, however, usually lower in the model than in the measurements. It cannot be determined whether the model underestimates precipitation or whether this effect is due to the redistribution of snow by the wind, which is not taken into account in the model but affects the accumulation at the measuring site. Significant differences between model and measurements also occur in cases of net ablation due to wind erosion or when accumulation at the stakes was due to snowdrift from the southwest without precipitation being observed.

We use simple numerical and analytical models of ice flow and heat flow to characterize the shape of isochrones and isotherms beneath moving ice divides. Both nonlinear ice flow and reduced accumulation (wind scouring) at a divide can cause reduced downward flow in a region about one ice thickness wide under a divide. Greater downward velocities on the flanks cause isochrones and isotherms to become arched at depth. The magnitudes and shapes of these arches depend on the history of divide position. Arch amplitudes decrease by approximately e−1 for each increase in migration rate of 3–5 times the accumulation rate, the arches become asymmetric, with steeper leading edges and more gentle trailing edges, and the arch apex lags behind the divide. Isochrone and isotherm shapes can be used to infer past divide motions. In advection-dominated ice sheets, isochrone shapes record a longer history of divide position than do isotherm shapes. The opposite is true for diffusion-dominated ice sheets, in which a spatial array of ice-temperature measurements might extend the recorded history of divide position.

There are several cases in which large overthrusts and sub-horizontal faults appear to have occurred in temperate or cold glaciers. As a contribution to solving the problem of their origin, the stress field when there is an abrupt change in the bottom boundary conditions is determined, assuming ice to be isotropic, third-power-law viscous. Deviatoric stresses vary with the distance r to the singularity as r−1/4, and strain rates as r−3/4. They are computed numerically to a multiplicative factor, which is determined by the conditions at a large distance, but not computed here. Although the apparent viscosity varies as r1/2 times a function of the polar angle, the stress field around the singularity is not essentially different from that obtained assuming a constant viscosity. Some considerations on the apparition of faults follow, but at the scale of the microrelief the adopted model becomes oversimplified.

Glacier de Saint-Sorlin, French Alps, left terminal moraines at 1.3, 2.9 and 3.7 km ahead of the present terminus. According to proxy data and to historical maps, these were formed in the 19th, 18th and 17th centuries, respectively. A plateau at 2700–2625 m was then surrounded by ice but never became an accumulation area. This fact shows that the equilibrium-line altitude (ELA) on the glacier never dropped below 2300 m. The following simple models apply sufficiently to yield reliable estimations of past ELA: (1) a uniform and constant vertical gradient of the mass balance, down to the terminus; and (2) a plane bed, with a slope of 8.5° and a uniform width. Then in a steady situation the accumulation–area ratio is 1/2. Compared to the mean for 1956–72, at the onset of the Little Ice Age the balances were higher by 3.75 m ice a−1, and the ELA was 400 m lower. Correlations between 1956–72 balances and meteorological data suggest that during the melting season the 0°C isotherm was about 800 m lower, while the winter precipitation at low altitudes did not change. These correlations may have been different in the past, but an equal lowering of the ELA and of the 0°C isotherm, as assumed by several authors, seems excluded.

Stake surveys on Glacier de Saint-Sorlin, French Alps, during the period 1957–76 show that annual surface velocities fit a linear vectorial model, with a term depending on the site, another on the year, and an important random component. Strain rates, viscosities and stresses at shallow depth are computed using strain triangles of hectometric (102 m) size. Between 1961/62 and 1972/73 the isotropic point, where streamlines cease to converge, moved downstream about 200 m. This trend may explain increasing velocities. Crevasses appear where annual strain is > 1.2%. Faults limit the effective shear stress at the surface to about 0.38 bar. Ten modes of flow are distinguished, instead of only two for the two-dimensional model (compressive and extensive). The gradients of shallow stresses, which ensure extra driving forces, are computed with another mosaic of triangles of similar size (stress triangles). There are also important extra driving forces at the bottom, which force the flow to deviate from the direction of the steepest surface slope. Two criteria allow elimination of stress triangles where these unknown basal extra driving forces are important. Even so, no sliding law in terms of mean annual values can be obtained. This study shows that the classical perturbation theory, which explains advances and retreats by the arrival of kinematic waves, is unsuitable for glaciers of kilometric size.

The 50 year time series of mass balance on Glacier de Sarennes is one of the longest in the French Alps, and so is often used as a reference for glacier variations in the French Alps. Meteorological data can be used to extend the series backwards in time. Martin (1978) proposed such a reconstruction for the 1882–1977 period. With 50 years of observations, we show that the classical method used by Martin is too dependent on the calibration period. We therefore try to improve the accuracy of this reconstruction using the Vincent and Vallon (1997) method which takes into account the albedo change of the surface during the ablation period (this is called the daily method). This new method appears to be stable in time. Once calibrated, the daily method is applied to reconstruct the 1881–1949 period. The new reconstruction is compared to a volumetric balance between two maps from 1906 and 1981. It appears that both reconstructions (classical and daily) fail to render the trend correctly over a long period of time. The cumulative centred mass balance correlates well (r2 = 0.62) with the hydrological mass-balance series of Aletschgletscher, Switzerland.

This paper describes the physical and chemical properties of the snow- pack on John Evans Glacier, Ellesmere Island, Canadian Arctic Archipelago, and investigates the controls on snowpack solute concentrations and atmospheric deposition. The snowpack contains three layers that are traceable across the whole glacier. These represent fall accumulation that has been metamorphosed to depth hoar, winter accumulation mixed with snow reworked by wind from the underlying depth hoar, and spring accumulation mixed with wind-reworked snow. The seasonal cycle in snow chemistry closely reflects changes in the composition of the atmospheric aerosol at Alert, with some modification of NO3− concentrations by post-depositional processes. Mean water-weighted solute concentrations in the snowpack are largely independent of accumulation, while atmospheric deposition tends to increase with accumulation. This suggests that, for most species, wet deposition is the dominant depositional process throughout the year. However, concentrations of Ca2+ and K+ increase with both accumulation and elevation, implying an enhanced input from dry deposition of soil dust above 800 m elevation. Concentrations of SO42− are inversely related to accumulation, especially in the winter layer, suggesting a significant input from non-precipitating events, such as dry deposition or riming, during this period of very limited snowfall.

A completely automated instrument for ice-fabric analysis has been built, tested and applied. The analysis of a thin section is performed in two steps. First, the section is loaded into the machine and is automatically imaged in about 16 min. Then, at a later time and convenient location, the grains in the section are identified from a stored image (either automatically or by the user), and the c-axis orientation of each grain is calculated at a rate of 1200 grains h−1. Comparison of automated results with manual Rigsby stage results shows excellent agreement. Diagnostics indicate that the accuracy of the new technique is 0.5° in the c-axis orientation. Two specific applications are demonstrated: (1) the ability to measure internal grain orientations that may reveal internal stress, and (2) the ability to measure very thin thin sections. The latter, combined with the ability to zoom in on small regions of a section, will be useful for examining the fabric of fine-grained (< 1 mm size grains) ice.