Late Neogene Sirius Group strata from Tillite Spur and Quartz Hills in the Reedy Glacier area, Antarctica, demonstrate the variability in Sirius Group facies and contrasts Sirius Group strata deposited at high and low paleo-elevation, respectively. The Tillite Spur and Quartz Hills Formations (Pliocene) are formally defined here.The Tillite Spur Formation type section crops out on the edge of the Wisconsin Plateau overlooking Tillite Spur. It comprises 32m of alternating coarse gray conglomerate and muddy olive-brown diamictites. The Quartz Hills Formation type section crops out above the western margin of Reedy Glacier in a pre-existing cirque towards the southern end of the Quartz Hills. It comprises c.100m of alternating massive diamictites and rhythmically interbedded sandstone and laminated mudstones which were deposited close to sea level and subsequently rapidly uplifted (>500 m Myr−1) to their present elevation at c. 1500 m. Three orders of paleoclimatic variability are recorded in the Sirius Group strata from Reedy Valley: (1) recycled marine microfloras in glacial diamictites indicate intervals of marine incursion into the Antarctic cratonic interior co-occurring with reductions in the East Antarctic ice sheet; (2) an advancing and retreating paleo-Reedy Glacier deposited a glacial/interglacial sequence alternating on a 10-100 kyr scale; 3) Centimeter and millimeter stratification in strata of the Quartz Hills Formation record annual kyr scale variability.

The set of force equations and stress strain-rate relations for ice masses can be solved with the method of lines and shooting the stress-free conditions at the free surface. Single- and multiple-shooting schemes with fixed point or Newton iterations for solving the stress equations including the deviatoric stress gradients are described and their performances arc discussed. The single-shooting Newton iteration proved to be the fastest seheme for typical valley glaciers, although its horizontal grid limitation is restrictive. Grid resolution can be improved substantially with a multiple-shooting scheme but computation time and storage requirements increase substantially. The Newton iteration allows the handling of mixed basal boundary conditions, partly basal velocity and partly basal shear traction being prescribed. A stick slip free gravity flow illustrates the performance of the scheme.

Numerical methods are used to examine the interaction between the spatial distribution of the basal shear traction and the corresponding basal velocity for an inclined slab geometry. In our improved treatment, we reject the common assumption that basal velocity is a simple function of local variables in favour of a non-local treatment that includes normal deviatoric stress and takes basal velocity to be an integrated response to spatially varying influences. Computationally, one must either iterate the basal velocity with a friction parameterization that relates basal shear traction to basal velocity or, alternatively, prescribe the basal shear traction that results from bed decoupling and substrate déformation.

The average of basal shear traction over the entire bed of the ice mass is invariant under changes in sliding distribution and thus constitutes a useful reference; any local relative reduction of traction leads to basal movement, either sliding over the bed or moving with a deforming subglacial layer. The local stress réduction is accompanied by a concentration of traction up-and down-glacier of the moving base. Growth, decay and possible migration of basal stress concentrations may be closely related to short-lived sliding events and to surges.

The annual layering in the Greenland and Antarctic ice sheets revealed by the b 0 record becomes less distinct with depth because of self-diffusion. However, the calculated diffusion rates are loo slow to explain die observations. It is suggested that the presence of veins of liquid water increases the effective diffusion constant by a factor of about 20.

The synthetic-aperture radar-interferometry technique is used to detect the migration of the limit of tidal flexing, or hinge line, of the floating ice longue of Peter-mann Gletscher, a major outlet glacier of north Greenland. The hinge line is detected automatically from differential interferograms using a model-fitting technique based on an elastic-beam theory. The statistical noise of the model fit is less than 3 mm, and the hinge line is mapped with a precision of 30 m. Following automatic registration of multidate image data to a precision of 5m, hinge-line migration is subsequently detected with a precision of 40 m in the horizontal plane across the glacier width. The results show that the hinge line of Petermann Gletscher migrates back and forth with tide by ±70m, in excellent agreement with the migration calculated from ocean tides predicted by a tidal model combined with the glacier surface and basal slope measured by an ice-sounding radar. Superimposed on the short-term hinge-line migration due to tide, we detect a hinge-line retreat of 270 m in 3.87 years which varies across the glacier width by ± 120 m. The retreat suggests glacier thinning at a rate of 78 + 35 cm ice a−1. Coincidentally, an analysis of ice-volume fluxes indicates that the hinge-line ice flux of Petermann Gletscher exceeds its balance flux by 0.88±1km3 ice a−1, which in turn implies glacier thinning at 83+95cm ice a−1 in the glacier lower reaches. Both methods therefore suggest that Petermann Gletscher is currently losing mass to the ocean.

Observations of slab-avalanche releases in alpine terrain have led to the hypothesis that rocky outcrops can influence the spatial distributions of temperature and heal flow in dry alpine snow covers and Unis control the local distribution of depth hoar. We investigate the effects of terrain on crystal growth by using a two-dimensional finite-clemcnl model of heal How coupled with a model of crystal growth from vapor. We used the model to examine the influence of snow properties, terrain geometry and snow depth on this phenomenon. The efleel is stronger in the early winter than in the late winter, because the rock has then had time to cool. In all cases, we found that depth-hoar growth occurs preferentially over the rock. This suggests that snow-pit investigations made over soil can he misleading if rocky outcrops are present.

The time evolution of a subglacial tunnel cross-section is examined usine a two-dimensional finite-element ice-flow model coupled to an idealized drainage system. Simulations are driven by physically based calculations of surface water-input variations at Slorgiaciaren, Sweden. Highlights of the model are its ability to handle unsteady conditions and irregular tunnel shapes. Agreement between modelled water pressure and borehole water levels is good. The following conclusions are reached: (i) Tunnels adapt to fluctuating inflow on time-scales of days. Storms, during which effective pressure ranges from 0 to 0.9 MPa, cause significant adjustments but daily fluctuations due solely to melt-water inflow are minor, (ii) Open-channel flow may become commonplace late in the ablation season, (iii) Initial tunnel shape influences subsequent tunnel evolution and seasonal water-pressure variation. Over the course of a summer, tunnels retain some of their initial shape, though in all experiments the width-to-height ratio increased with time, (iv) Tunnel contraction forms broad low tunnels. However, (v) given two tunnels of equal initial area, the higher narrower one expands more rapidly. Thus, more semi-circular tunnels may capture How from broader neighbours early in the summer.

As part of the winter environment in middle- and high-latitude regions, the interactions between wind, vegetation, topography and snowfall produce snow covers of non-uniform depth and snow water-equivalent distribution. A physically based numerical snow-transport model (SnowTran-3D) is developed and used to simulate this three-dimensional snow-depth evolution over topographically variable terrain. The mass-transport model includes processes related to vegetation snow-holding capacity, topographic modification of wind speeds, snow-cover shear strength, wind-induced surface-shear stress, snow transport resulting from saltation and suspension, snow accumulation and erosion, and sublimation of the blowing and drifting snow. The model simulates the cold-season evolution of snow-depth distribution when forced with inputs of vegetation type and topography, and atmospheric foreings of air temperature, humidity, wind speed and direction, and precipitation. Model outputs include the spatial and temporal evolution of snow depth resulting from variations in precipitation, saltation and suspension transport, and sublimation. Using 4 years of snow-depth distribution observations from the foothills north of the Brooks Range in Arctic Alaska, the model is found to simulate closely the observed snow-depth distribution patterns and the interannual variability.

Force variations on a "ploughmeter" and fluctuations in subglacial water pressure have been measured in the same borehole at Storglaciaren, Sweden, to investigate hydraulic properties of the basal till layer. A strong inverse correlation of the pressure and force records, in conjunction with a significant lime lag between the two signals, suggests that pore-water pressures directly affect the strength of the till. Variations in sub-glacial water pressure result in potential gradients across the water till interlace at the bottom of the borehole that drive pressure waves downwards through the till layer when the borehole water level is high and back upwards when the water level is low. Analysis of the propagation velocity of this pressure wave indicates that the hydraulic diffusivity of Storglaciaren till is in the range 1.9−3.6 x 10−6m2s−1,in good agreement with estimates obtained in the laboratory. Hydraulic conductivity values associated with these difrusivities are between 10−9 and 10−8ms−1 and thus are well within the range of values for other glacial tills.

The relationships between accumulation and microwave emissivity in the dry-snow regions of the Greenland ice sheet are studied with special focus on hoar formation. With the use of a discrete-ordinate radiative-transfer model, the relative contributions of hoar and accumulation to microwave omission are examined. The results show that the link between microwave emission and accumulation rates in the dry-snow area of the Greenland ice sheet is significant; however, the emissivity is also dependent on the extent of hoar development. As a result, accumulation estimates based on passive microwave observations in the low-accumulation dry-snow areas will require successful parameterization of hoar formation. Furthermore, the increasing Tb, trends that were observed during the 1980s, which have been hypothesized to be related to accumulation changes, may be strongly related to variations in hoar development. Other possible explanations for these observed trends, such as slight surface melt or atmospheric conditions, are also discussed.

Differential spaceborne radar interferometry observations of West Bagley Icefield are used to measure surface velocity and topography. Bagley Icefield is the accumulation area fur Bering Glacier which surged in two phases from spring 1993 through summer 1905. The observations presented are based on data collected during the winter of 1992, prior to the surge, and during winter 1994 while the surge was in full progress. Both observation intervals correspond to 3 day repeat orbit phases of the ERS-I C-band SAR. This paper gives an overview of the algorithms used to derive surface-velocity vector fields and topography for valley glaciers from SAR images. The resulting high-resolution velocity data clearly show West Bagley Icefield accelerating from its quiescent pre-surge velocity by a factor of 2.7 in response to the Bering Glacier surge. Persistence of lnterfero-metric phase coherence and the relatively moderate degree of acceleration on the western arm of Bagley Icefield suggest that the velocity increase may have been caused by increased longitudinal stress gradients resulting from coupling to the surging main trunk of Bering Glacier.

Debris-laden ice accretes to the base of Matanuska Glacier, Alaska, U.S.A., from water that supercools while flowing in a distributed drainage system tip the adverse slope of an overdeepening. Frazil ice grows in the water column and forms aggregates, while other ice grows on the glacier sole or on substrate materials. Sediment is trapped by this growing ice, forming stratified debris-laden basal ice. Growth rates of >0.l ma−1 of debris-rich basal ice are possible. The large sediment fluxes that this mechanism allows may have implications for interpretation of the widespread deposits from ice that flowed through other overdeepenings, including Heinrich events and the till sheets south of the Laurentian Great Lakes.

Simple theory supports field observations (Lawson and others, 1998 that subGlaciol water flow out of overdeepenings can cause accretion of layered, debris-bearing ice to the bases of glaciers. The large meltwater flux into a temperate glacier at the onset of summer melting can cause rapid water flow through expanded basal cavities or other flow paths. If that flow ascends a sufficiently steep slope out of an overdeepèning, the water will supercool as the pressure-melting point rises, and basal-ice accretion will occur. Diurnal, occasional or annual fluctuations in water discharge will cause variations in accretion rate, debris content of accreted ice or subsequent diagenesis, producing layers. Under appropriate conditions, net accretion of debris-bearing basal ice will allow debris fluxes that are significant in the glacier sediment budget.

Airborne surface elevation profiles of the Harding Ice field, south-central Alaska, were made in 1991 and 1996. Thirteen glaciers were profiled, along with the tipper region of the icefield. The profiles were compared to U.S. Geological Survey topographie maps made in the 1950s, to obtain elevation and volume changes. Comparison of the changes for the different glaciers shows no significant correlation between volume change and the type of glacier or characteristics such as location, aspect, size, slope or terminus changes. Estimated total volume change tor this ~43 year period is about -34 km3, which corresponds to an area-average elevation change of-21 m. The estimated error in this elevation change of 5 m is mainly due to errors in the maps at higher elevations. Our measurements provide an accurate baseline against which future determinations of volume change can be made.

Ablation climate studies were made at two locations in northern Greenland in the summers of 1993 and 1994, respectively. Daily ablation was measured at ten stakes within a small area, and the data were compared with each other to detect gross errors. For example, high standard deviations for data taken on the same day, or low correlations between data series at different stakes, indicate erroneous data. After discarding data for one stake in 1993 and two stakes in 1994, random errors in daily ablation data for individual stakes are ± 5 kg m 2 d−1, which is further reduced to only about db 2 kg m−2d−1 by averaging over eight or nine stakes. Random errors in calculated energy balances using the present ablation data are much lower than found in earlier stuthes in West Greenland where ablation was only measured on three stakes without any attempt to detect gross errors. Aside from day-to-day errors, there are ±10% differences in mean ablation at different stakes, which are probably caused by small-scale variations in surface albedo. Such interstate differences give a ± 10% uncertainty in positive degree-day factors, which are 9.8 ± 0.9 and 5.9 ± 0.6 kg m 2 d −1 deg −1 for the two sites.

This paper addresses the coupling of flows of ice, till and water, and the issue of whether such coupling provides mechanisms for meso-scale (kilometres to tens of kilometres) variability in ice-sheet flow and texture. The question of whether effective pressures at the ice-bed interface are statically or hydraulically controlled is examined in this paper. The answer is scale dependent, and has a significant effect on the relationship between ice surface and basal topography.

The consequences of these considerations on till flow, coupled ice–till flow and coupled ice, till and water flow are examined. An analysis of till-flow kinematics and shock formation is carried out. The linear stability of coupled long-wavelength ice-till flow is analysed, and regions in parameter space where this flow is unstable, with rather small rate constants are found. Upstream-moving ice surface waves are predicted. The linear stablity of coupled ice–till–water flow is examined, where water flow is modelled using a basal flow system with effective-pressure-dependent properties. Again, regions in parameter space where the system is linearly unstable are found, this time with relatively rapid rate constants. The water pressure exhibits “breather” modes.

These analyses assume that there is a substantial basal traction. A problem with models of ice streams wholly restrained at the side is identified: they seem to predict erosion rates which are unfeasibly large.

There appears to be sufficient variability in the ice–till– water system to potentially explain texture in ice-stream surfaces, variations in ice-stream thickness of tens of metres not directly relatable to topography, and waves moving upstream or downstream. Most importantly, the ice-stream–bed system is shown to exhibit meso-scale variability simply by coupling ice flow according to the shallow-ice approximation, till flow according to the hydrostatic thin-till approximation and water flow according to an effective-pressure-dependent hydraulics.

Ice temperature was measured in and around the chaotically crevassed south margin of Ice Stream B, Antarctica, from 1992 to 1994. The temperatures at 30 m depth in the chaotic zone are about 12 K lower than in the adjacent uncrevassed ice, due to the ponding of cold winter air. At depths greater than 150 m, there is clear evidence of internal heating of the ice due to the large shear déformation rate in the marginal zone. Analysis of the depth of cooling below the crevasses and of the internal heating gives two pieces of information. First, over the last half century the lateral shear stress averaged 2.0 x 105 Pa in the top third of the margin and, second, the margin moved outward at an average rate of 7.3 m a−1. These values do not involve any assumptions about the How law of ice. The uncertainties are roughly 20%. The value of lateral shear stress indicates that the most of the drag on the ice stream is along its sides.

Using an active transponder with the ERS-I and ERS-2 radar altimeters, the distance to the satellite was measured at a location close to the GRIP site, Greenland, at an altitude of 3.2 km. The measurement was executed while the transponder was in the “ice-tracking mode”. It includes a bias due to the propagation delay. The location of the transponder was determined using the global positioning system.

The transponder signal was modeled and the distance from the altimeter to the effective reflection point of the transponder was determined. Since the transponder was located within 1 km of the ground tracks, the measurement was corrected for this offset. A correction was also done for the surface slope, resulting in the distance (plus bias) to the closest sub-satellite point on the surface of the (compact) snow.

The transponder signal was then removed from the radar altimeter waveform, enabling the determination of the distance (plus bias from the altimeter to the first reflective surface within the snow. The différence between this distance and that obtained using the transponder was < 2 m. This shows that the surface which gives rise to the first return of the reflection agrees with the surface of the (compact, dry) snow at this high-altitude location. This is an important result to be used when studying ice-cap topography using satellite radar altimetry.

A ring-shear device was used to study the factors that control the ultimate(steady) strength of till at high shear strains.Tests at a steady strain rate and at different stresses normal to the shearing direction yielded ultimate friction angles of 26.3° and 18.6° for tills containing 4% and 30% clay-sized particles, respectively Other tests at steady normal stresses and variable shear-strain rates indicated a tendency for both tills to weaken slightly with increasing strain rate. This weakening may be due to small increases in till porosity.

These results provide no evidence of viscous behavior and suggest that a Coulomb-plastic idealization is reasonable for till deformation. However, viscous behavior has often been suggested on the basis of distributed shear strain observed in subglacial till. We hypothesize that deformation may become distributed in till that is deformed cyclically in response to fluctuations in basal water pressure. During a deformation event, transient dilation of discrete shear zones should cause a reduction in internal pore-water pressure that should strengthen these zones relative to the surrounding till, a process called dilatant hardening. Consequent changes in shear-zone position, when integrated over time, may yield the observed distributed strain.

The non-linearity of the ice-flow law or a local accumulation low over an ice divide can cause isochrones (internal layers) to be shallower under the divide relative to the flanks, forming a “divide bump” in the internal layer pattern. This divide signature is analyzed using ice-flow models and inverse techniques to detect and quantify motion of the Siple Dome ice divide, West Antarctica. The principal feature indicating that migration has occurred is a distinct tilt of the axis of the peaks of the warped internal layers beneath the divide. The calculated migration rate is 0.05-0.50 m a−1 toward Ice Stream D and depends slightly on whether the divide bump is caused by the non-linearity of ice flow or by a local accumulation low. Our calculations also suggest a strong south-north accumulation gradient of 5-10 x 10−6 a−1 in a narrow zone north of the divide. A consequence of divide migration is that pre-Holocene ice is thickest about 0.5 km south of the present divide position. Divide motion indicates that non-steady processes, possibly associated with activity of the bounding ice streams, are affecting the geometry of Siple Dome. The migration rate is sufficiently slow that the divide bump is maintained in the internal layer pattern at all observable depths. This suggests that major asynchronous changes in the elevation or position of the bounding ice streams are unlikely over at least the past 103-104 years.