Hostname: page-component-7c8c6479df-ph5wq Total loading time: 0 Render date: 2024-03-28T16:47:05.234Z Has data issue: false hasContentIssue false

Distributed shear of subglacial till due to Coulomb slip

Published online by Cambridge University Press:  08 September 2017

Neal R. Iverson
Affiliation:
Department of Geological and Atmospheric Sciences, Iowa State University, Ames, Iowa 50011, U.S.A.
Richard M. Iverson
Affiliation:
U.S. Geological Survey, 5400 MacArthur Boulevard, Vancouver, Washington 98661, U.S.A.
Rights & Permissions [Opens in a new window]

Abstract

In most models of the flow of glaciers on till beds, it has been assumed that till behaves as a viscoplastic fluid, despite contradictory evidence from laboratory studies. In accord with this assumption, displacement profiles measured in subglacial till have been fitted with viscoplastic models by estimating the stress distribution. Here we present a model that illustrates how observed displacement profiles can result from till deformation resisted solely by Coulomb friction. Motion in the till bed is assumed to be driven by brief departures from static equilibrium caused by fluctuations in effective normal stress. These fluctuations result from chains of particles that support intergranular forces that are higher than average and that form and fail at various depths in the bed during shearing. Newton’s second law is used to calculate displacements along slip planes and the depth to which deformation extends in the bed. Consequent displacement profiles are convex upward, similar to those measured by Boulton and colleagues at Breidamerkurjökull, Iceland. The model results, when considered together with the long-term and widespread empirical support for Coulomb models in soils engineering, indicate that efforts to fit viscoplastic flow models to till displacement profiles may be misguided.

Type
Research Article
Copyright
Copyright © International Glaciological Society 2001

Introduction

Deformation of unlithified sediment beneath wet-based ice masses may contribute significantly to the motion of modern ice streams (e.g. Reference Alley, Blankenship, Bentley and RooneyAlley and others, 1986) and surging glaciers (e.g. Reference Clarke, Collins and ThompsonClarke and others, 1984). Bed deformation may have also enabled fast flow of Pleistocene ice masses, perhaps contributing to the 100 000 year climate cycle, to extreme millennial climate variability and to the low aspect ratios of ice sheets during the Last Glacial Maximum (Reference Clark, Alley and PollardClark and others, 1999). There is therefore widespread interest in modeling accurately the flow of glaciers that rest on sediment.

In most models, the till that commonly underlies glaciers has been assumed to behave as a viscoplastic fluid, such that at shear stresses exceeding its Coulomb yield strength, it shears at rates proportional to the excess shear stress (Reference Boulton and HindmarshBoulton and Hindmarsh, 1987; Reference AlleyAlley, 1989; Reference Walder and FowlerWalder and Fowler, 1994; Reference BoultonBoulton, 1996; Reference Clark, Licciardi, MacAyeal and JensonClark and others, 1996; Reference Jenson, MacAyeal, Clark, Ho and VelaJenson and others, 1996; Reference HindmarshHindmarsh, 1998a–c; Reference Licciardi, Clark, Jenson and MacAyealLicciardi and others, 1998). Support for this assumption has come from the field measurements of Boulton and co-workers at Breidamerkurjökull, Iceland (Reference Boulton and HindmarshBoulton and Hindmarsh, 1987; Reference Boulton and DobbieBoulton and Dobbie, 1998). By inserting segmented rods and cable extensometers in till beneath the glacier margin, they measured displacement as a function of depth and found that most glacier motion resulted from pervasive shearing of the upper 0.5–1.0 m of till. Rates of shear strain increased upward toward the glacier sole. These displacement profiles have been interpreted as evidence of linear or slightly non-linear viscoplastic behavior (stress exponent < 2.5) (Reference AlleyAlley, 1989; Reference Boulton and DobbieBoulton and Dobbie, 1998).

As was first noted by Reference KambKamb (1991), however, most experimental data from soil mechanics (e.g. Reference MitchellMitchell, 1993, p.343) and geophysics (e.g. Reference Biegel, Sammis and DieterichBiegel and others, 1989) contradict such inferences. These data indicate that strengths of slowly shearing granular materials do not increase significantly with deformation rate, contrary to the viscoplastic behavior inferred from displacement profiles. Laboratory experiments with various testing devices on a variety of tills confirm that the strength of till is insensitive to its deformation rate (Reference Tika, Vaughan and LemosTika and others, 1996; Reference Iverson, Hooyer and BakerIverson and others, 1998; Reference Tulaczyk, Kamb and EngelhardtTulaczyk and others, 2000). These studies indicate that till is best idealized as a Coulomb material, such that its strength is independent of deformation rate and linearly dependent on effective normal stress. This is not a surprising conclusion: soils engineers have characterized soil deformation successfully with Coulomb models for many decades (e.g. Reference Terzaghi, Peck and MesriTerzaghi and others, 1996), and recent studies of till deformation beneath modern glaciers have yielded results that are difficult to explain with viscoplastic models (e.g. Reference Hooke, Hanson, Iverson, Jansson and FischerHooke and others, 1997; Reference Truffer, Harrison and EchelmeyerTruffer and others, 2000).

Despite the clear motivation for Coulomb models, viscous deformation resistance continues to be assumed in models of subglacial till deformation (Reference NgNg, 2000; Reference Thorsteinsson and RaymondThorsteinsson and Raymond, 2000). This, in part, may reflect the lack of an explicit and simple analytical description of how distributed till deformation can occur in subglacial till without viscous deformation resistance.

The objective of this investigation is to demonstrate that displacement profiles like those observed in till beneath Breidamerkurjökull may arise from deformation resisted only by Coulomb friction and that, therefore, a viscoplastic till rheology cannot be inferred confidently from such profiles. Our analysis also indicates that neither dilatant hardening (Reference Iverson, Hooyer and BakerIverson and others, 1998) nor pore-pressure diffusion into the bed from the ice–till interface (Reference Tulaczyk, Mickelson and AttigTulaczyk, 1999), two proposed mechanisms for distributing deformation in Coulomb models, is necessary to account for the observed distribution of displacement.

Background

Yielding will occur in a granular material when the shear traction applied to a planar element in the material equals the Coulomb yield strength, proportional to the effective normal stress on the element. Everywhere the material is near this condition, called limiting equilibrium, it will be on the verge of deforming irreversibly. As many others have noted (e.g. Reference Boulton and JonesBoulton and Jones, 1979; Reference Alley, Blankenship, Bentley and RooneyAlley and others, 1986), for a subglacial till layer to be weakened sufficiently to approach this limiting state, till pore-water pressure must be a large fraction of the total normal stress.

Laboratory experiments, motivated by diverse problems in physics, geophysics and engineering, provide insight into how deformation occurs as limiting equilibrium is approached. These studies highlight the stress heterogeneity that characterizes granular materials and the temporal variability of intergranular stresses during deformation. Experiments with particles in static and slowly shearing packs indicate that stresses are concentrated along “force chains” which coalesce and form complex networks that include many particles (Fig. 1a) (e.g. Reference Dantu and GlanvilleDantu, 1957; Reference Drescher and de Josselin de JongDrescher and de Jong, 1972; Reference Allersma, Vermeer and LugerAllersma, 1982; Reference LiuLiu and others, 1995; Reference Painter, Tennakoon, Behringer, Herrmann, Houi and LudingPainter and others, 1998; Reference Veje, Howell, Behringer, Schollann, Luding, Herrmann, Herrmann, Houi and LudingVeje and others, 1998; Reference Howell, Behringer and VejeHowell and others, 1999). During shearing, these chains reorganize continuously and cause intergranular stresses to fluctuate significantly (Fig. 1b) (e.g. Reference Painter, Tennakoon, Behringer, Herrmann, Houi and LudingPainter and others, 1998; Reference Howell, Behringer and VejeHowell and others, 1999). Such fluctuations have been measured normal to the shearing direction in ring-shear experiments with simulated and natural tills (Fig. 1c) (Reference Iverson, Hooyer and HookeIverson and others, 1996, Reference Iverson, Baker and Hooyer1997; Reference HookeHooke, 1998, fig. 7.16). Not surprisingly, these spatial and temporal variations in intergranular stresses give rise to instantaneous velocity fields that are highly heterogeneous, with groups of particles pushed forward sporadically along surfaces of concentrated particle slip or rolling (Reference Drescher and de Josselin de JongDrescher and de Jong, 1972; Reference Misra, Herrmann, Houi and LudingMisra, 1998).

Fig. 1. (a) Force chains in 2900 photoelastic (birefringent under strain) disks sheared within a Couette apparatus. The outside diameter of the sample chamber is 0.38 m. Disks are 7 and 9 mm in diameter. Lighter disks are under larger forces. The rate of shearing at the specimen center line is 2 mm s−1 (adapted from Reference Howell, Behringer and VejeHowell and others, 1999). (b) G 2, an optical measure of the force on disks proportional to the number of birefringent bands they display, as a function of displacement during an experiment like that shown in (a). Over the measured range of G 2, its value is nearly linearly proportional to the force on disks. The value plotted is an average over one-tenth of the samplers circumference (adapted from Reference Howell, Behringer and VejeHowell and others, 1999). (c) Local stresses measured normal to the shearing direction at the upper boundary of a till specimen during deformation in a ring- shear device. Forces are measured over circular platens with diameters of 18 and 36 mm. The rate of shearing is 0.01 mm s−1 (adapted from Reference Iverson, Baker and HooyerIverson and others, 1997).

Results of distinct-element numerical models, in which equations of motion (forces and moments) are solved simultaneously for every particle assuming appropriate particle-contact rules (e.g. Reference Cundall and StrackCundall and Strack, 1979), reinforce these laboratory results. Simulations of slow deformation with such models indicate that force chains always develop with the associated tendency for grains to move and rotate sporadically in coherent clusters (Reference Cundall, Drescher, Strack, Vermeer and LugerCundall and others, 1982; Reference Misra, Chang, Babic, Liang and MisraMisra, 1997; Reference Williams, Rege, Chang, Babic, Liang and MisraWilliams and Rege, 1997). In especially relevant simulations of fault-gouge deformation in a shear zone, Reference Morgan and BoettcherMorgan and Boettcher (1999) applied a distinct-element model to grain-size distributions similar to those of till and found that force chains formed and failed continuously, resulting in local stress fluctuations. Shearing was caused by the failure of force chains and focused along planes that became parallel to the shearing direction at low strains (0.1). Although strain accrued through sporadic slip along individual planes, the cumulative effect of these slip events (shear strain = 2.0) was a relatively uniform distribution of finite strain.

Hypothesis

We explore the hypothesis that local stress fluctuations like those observed in experiments cause bulk deformation of till subglacially (e.g. Reference Boulton and HindmarshBoulton and Hindmarsh, 1987). Mean shear and effective normal stresses on a till layer control whether it is near limiting equilibrium and on the verge of deforming. However, particle movement at a particular depth is driven by brief reductions in effective normal stress caused by the failure of force chains. Consequent force imbalances on groups of particles push these groups forward episodically along slip planes at various depths, resulting in bulk deformation of the bed over time.

The objective of our analysis is to calculate the distribution of motion with depth in the bed that results from these local departures from static equilibrium. Rather than consider a viscoplastic till rheology, we assume simple Coulomb friction, in which rate-independent slip between grains depends linearly on effective stress. Non-synchronous till displacements along hypothetical slip planes are calculated using conservation of momentum, which also dictates the maximum depth to which deformation extends in the bed. Integrating till displacements from this depth to the ice–bed interface yields the bed-normal displacement profile.

Our general approach is similar to that of distinct-element models in that, rather than apply a continuum model for flow, we solve the equation of motion for discrete portions of the till while applying an appropriate friction relation for slip. Although our simple analytical model falls well short of the realism achieved using numerically intensive alternatives like the distinct-element method, it illustrates with a minimum of complexity an important concept in interpreting the subglacial deformation of till: that fluid-like deformation profiles need not be a consequence of a fluid rheology.

Analysis

Consider a thick till bed that slopes at an angle β (Fig. 2). The depth, y, below the ice–till interface is measured normal to the bed surface. Motion in the bed is assumed to occur along a series of bed-parallel slip planes separated by a uniform distance, δ, much smaller than the bed thickness. Although local slip surfaces in granular media are sometimes inclined at an angle to the macroscopic shearing direction, experimental evidence (Reference Mandl, de Jong and MalthaMandl and others, 1977; Reference Logan, Dengo, Higgs, Wang, Evans and WongLogan and others, 1992; Reference Beeler, Tullis, Blanpied and WeeksBeeler and others 1996) indicates that at sufficiently large strains slip surfaces tend to align with the shearing direction, consistent with the results of distinct-element models (Reference Morgan and BoettcherMorgan and Boettcher, 1999). Thus, idealizing the bed as essentially a deck of cards, as we do here, is a reasonable abstraction. The till is assumed to be water-saturated with a constant bulk density, ρ t. Possible ramifications of bulk density varying temporally are discussed later in the paper.

Fig 2. Model bed geometry.

As in distinct-element models of granular shearing, we calculate displacements using Newton’s second law. In the most general case, slip in the bed may occur along neighboring slip planes concurrently. Thus, a slice of till bounded by slip planes may have a frame of reference that is accelerated relative to that of the till slice beneath it. For this case, if V r is the bed-parallel velocity of a till slice relative to the slice immediately underlying it and t is time, the equation of motion for slip is

(1)

where F is the net bed-parallel force acting on till and ice above the slip plane, m is the mass of till and ice above the slip plane, and V is the velocity of the till slice beneath the slip plane relative to a fixed reference frame. The second term on the righthand side of this equation results from the accelerated reference frame (Reference AharoniAharoni, 1972, p. 192). This term can be eliminated, however, if we assume that slip events at various depths occur at discrete times. Distinct-element modeling of fault-gouge deformation by Reference Morgan and BoettcherMorgan and Boettcher (1999) indicates that, although more than one slip surface can be active at any time in a shear zone, there is a general tendency for slip to alternate from one plane to another. Similar behavior has been assumed in models of subglacial till deformation (Reference Tulaczyk, Mickelson and AttigTulaczyk, 1999). For non-synchronous slip events, till below an active slip plane will always be at rest, and Equation (1) reduces to the standard form of Newton’s second law:

(2)

A net driving force F on till and ice above a particular slip plane exists when the Coulomb shear strength becomes less than the downslope weight per unit area of the overlying till and ice. If σ e is the effective normal stress on planes at depth y, and φ is the Coulomb friction angle of the till, then the Coulomb shear strength is σ e tan φ. The till is considered to be deformed sufficiently to be in its critical state, and thus φ is a constant and cohesion is negligible (Reference SkemptonSkempton, 1985). If W d is the downslope component of the weight per unit area of the overlying till and ice at depth y, the net force per unit area that causes motion is W d =− σ e tan φ. If H is the ice thickness, the mass per unit area of the till and overlying ice at depth y is ρ t(y + CH), where C = ρ/ρ t and ρ is the density of ice. Incorporating these relations in Equation (2) yields

(3)

V r can only be positive because up-glacier motion does not occur if W d < σ e tan φ.

To characterize σ e and W d as functions of depth, we first emphasize that σ e has a fluctuating component due to the formation and failure of force chains. Motion is expected to occur when force chains fail (e.g. Reference HookeHooke, 1998, p. 109; Reference Morgan and BoettcherMorgan and Boettcher, 1999). Thus, we define σ e such that it includes a transient perturbation, . Such perturbations are expected to occur randomly with depth and to have a wide range of magnitudes (Reference Painter, Tennakoon, Behringer, Herrmann, Houi and LudingPainter and others, 1998; Reference Veje, Howell, Behringer, Schollann, Luding, Herrmann, Herrmann, Houi and LudingVeje and others, 1998). Rather than attempt to characterize these perturbations statistically, which would add realism but sacrifice clarity, we make the simplest possible assumption: that perturbations, although not concurrent, are of equal magnitude at all depths. If is taken to be positive when the effective normal stress is increased, then the effective stress in the till on any plane normal to y during a perturbation is:

(4)

where g is the gravitational acceleration, σ 0 is the total normal stress at the ice-till interface, and λ is the ambient pore-pressure ratio, which equals the fraction of the total normal stress supported by pore-water pressure in the absence of an effective-stress perturbation. The value of λ can vary from 0 to 1 and accounts in a general way for the full range of possible pore-water pressures, given that the steady-state pore pressure is not necessarily hydrostatic. The downslope component of the weight per unit area of ice and till at depth y is

(5)

where W 0 and ρ t gy sin β are the contributions due to the weights of the ice and till, respectively. Substitution of Equations (4) and (5) into Equation (3) yields

(6)

To condense Equation (6), we define the following parameters that describe the state of stress on potential slip planes:

(7a)

(7b)

(7c)

where S 0 is the difference between the downslope component of the weight of the ice per unit area and the available shear strength of the bed at y = 0 in the absence of an effective-stress perturbation, S′ is the shear-strength perturbation caused by the effective-stress perturbation, and a is a constant (0 < α < 1) that characterizes the steady-state stress due to the weight of the till. Since we are assuming that the bed is at rest unless perturbed sufficiently, the shear strength of the bed exceeds the downslope weight per unit area of ice, and thus S 0 ≤ 0. Only a reduction in effective normal stress, corresponding to a negative value of , can cause slip, and thus we consider only the case in which S′ ≥ 0. Using Equations (7a–c) allows the equation of motion (Equation (6)) to be written in a compact form:

(8)

If the shear-strength perturbation, S″, is sufficiently large (S″ + S 0 > ρ t gαy)), W d (Equation (5)) exceeds resisting stresses at depth y, slip occurs there in the till, and deformation will extend to a depth in the bed, y 0, that depends on the magnitude of the numerator of Equation (8). Note that at this depth V r will always equal zero, so this depth can be determined by setting the lefthand side of Equation (8) equal to zero and solving for y to obtain

(9)

Deformation decreases to zero at depth because the bed shear strength, controlled by intergranular friction, increases more rapidly with depth than W d. At depth y 0, friction is sufficiently large that W d is balanced exactly by the till shear strength, and there is no motion.

Next we determine the relative velocity between neighboring till slices for a given shear-strength perturbation, S′, by solving Equation (8). Two time periods are relevant: the period during the shear-strength perturbation when the till slice (and overlying till and ice) are accelerating and the period immediately after the shear-strength perturbation when the till slice is decelerating but still moving owing to the inertia of the till and ice. For the duration of the perturbation, V r is obtained by integrating Equation (8) and applying an initial condition which states that motion commences from a state of rest. The integration yields

(10)

where T is the duration of the perturbation. During the second stage of motion, S′ = 0 because the perturbation has ceased. Integrating Equation (8) with S′ = 0 and applying Equation (10) as an initial condition at t = T yields

(11)

where t stop is the time required for till at a particular depth y to come to rest. This time is determined by setting V r = 0 in Equation (11) and solving for t to obtain

(12)

Note that values of t stop diminish with depth, indicating that motion events caused by effective-stress perturbations near the glacier sole persist longest. This persistence indicates that effective-stress perturbations of a given magnitude induce larger force imbalances near the glacier sole.

The bed-normal profile of displacement can be obtained by first determining the relative displacement, X r(y), between till slices. To do so, we integrate Equation (10) from 0 to T and Equation (11) from T to t stop, and add the two results to obtain

(13)

Writing this equation in terms of t stop (Equation (12)) yields a simpler expression,

(14)

which illustrates that at large depths where t stop approaches T, X r(y) approaches zero. At depth y 0, T = t stop, and X r(y) = 0.

To calculate the total displacement, X(y), the relative displacements on slip surfaces between y 0 and y must be summed. The number of slip surfaces is given by the range of depth over which the summation is performed, y 0y, divided by the distance between slip surfaces, δ. Thus, keeping in mind that δ is assumed to be very small, the summation is expressed by the integral,

(15)

where the negative sign accounts for the negative (upward) direction of the integration. Substituting Equation (13) into this equation and integrating results in a lengthy expression for the desired displacement profile (see Appendix).

For the important special case in which the downslope component of the weight per unit area of the glacier (W 0) and the bed shear strength are balanced sufficiently at the ice–till interface to assume S 0 = 0, Equation (13) can be written conveniently in terms of y 0 (Equation (9) ) to obtain

(16)

Substituting this equation into Equation (15) and integrating yields the displacement profile for this special case:

(17)

Equation (17) is similar to the general expression presented in the Appendix, but that expression contains a more complex logarithmic term and cannot be written as concisely in terms of y 0. Note that displacements indicated by Equations (16) and (17) become infinite as y → 0, so the model results are meaningful only at depths yδ/2. Displacements at y = 0 depend on interactions between basal ice and till particles that are not included in the model.

Equation (17) indicates that total displacements depend on the square of the duration of the shear-strength perturbation, T 2, and depend non-linearly on the magnitude of the shear-strength perturbation, S′, through dependence on y 0. Displacements are largest if the slip-plane spacing δ is small because δ is inversely related to the number of slip planes, and the total displacement of a given till slice increases with the number of slip planes beneath it. Note also that the normalized total displacement, X(y)/y 0, scales with the dimension-less parameter αgT 2/2δ.

Parameter Choices and Results

The objective now is to explore whether profiles generated by the model are consistent with those observed at Breidamerkurjökull, Iceland. Reference Boulton and DobbieBoulton and Dobbie (1998) presented the most complete data on the distribution of till motion beneath that glacier (Fig. 3a). Through a borehole about 1 km from the glacier margin, these authors placed a displacement marker at the level of the glacier sole, and four other markers at various depths in the underlying till. The markers were connected with wires to potentiometer-gauged spools that were held in the borehole at the base of the ice. As the bed sheared, the spools were carried away from the markers, and the length of wire paid out was monitored for 17 days to obtain the distribution of displacement in the bed over that period. Unfortunately, certain details of the experiment were not discussed, such as how the spools were held fixed at the bottom of the borehole and how the depth of marker emplacement in the bed was measured. In addition, the trajectories of the wires between the markers and the spools were unknown as a function of time. Owing to this significant uncertainty, others have restricted use of this technique to measurements in the bed only very near the glacier sole (Reference Blake, Fischer and ClarkeBlake and others, 1994; Reference Engelhardt and KambEngelhardt and Kamb, 1998). Despite these uncertainties, however, inferred profiles were qualitatively similar to those measured directly by excavation nearer to the glacier margin (Reference Boulton and HindmarshBoulton and Hindmarsh, 1987), with shear strain increasing upward toward the glacier sole. For the purposes of this paper, therefore, the displacement profiles of Figure 3a will be taken at face value.

Fig. 3. (a) Daily displacement of anchors in the bed of Breidamerkurjökull over a 17 day period, as reconstructed by Reference Boulton and DobbieBoulton and Dobbie (1998, fig 13). (b) Daily displacement indicated by the model over 17 days, assuming that limiting equilibrium was approached each day during one period of sufficient duration for a perturbation in effective normal stress to occur once at each depth in the bed. Parameter values used in the model are listed in Table 1.

Many of the parameters of the model are known from field observations at Breidamerkurjökull and from ancillary laboratory measurements (Reference Boulton and HindmarshBoulton and Hindmarsh, 1987; Reference Boulton and DobbieBoulton and Dobbie, 1998). These are H, σ 0 , ρ t , C, β and φ (Table 1). Other parameters, such as λ and W 0, can be estimated. Although pore-water pressure, recorded by sensors in the displacement markers, fluctuated through a range of about 200 kPa (Reference Boulton and DobbieBoulton and Dobbie, 1998), averaging over the period of record indicates that λ = 0.9 is a good approximation. The value of W 0 is specified such that it equals the resisting stress at the bed surface in the absence of effective-stress perturbations. This is a convenient but not essential assumption that satisfies the special case noted earlier in which S 0 = 0, and allows use of Equation (17) in computing the velocity profile. Thus, using values from Table 1 and Equation (7a), W 0 = σ 0(1 − λ) tan φ = 58.7 kPa, a reasonable value. An alternative would be to calculate W 0 from the glacier thickness and surface slope, but that approach would be suspect owing to unknown longitudinal stress gradients near the glacier margin.

Table 1. Model parameters used to generate Figure 3b

Other important parameters in the model — the spacing of slip surfaces, δ, the shear-strength perturbation, S′, and its duration, T — can be estimated only within broad limits. The value of δ must be sufficiently large to approximate the movement of many grains in coherent groups (e.g. Reference Drescher and de Josselin de JongDrescher and de Jong, 1972; Reference Misra, Chang, Babic, Liang and MisraMisra, 1997) but sufficiently small to result in time-integrated strain that is pervasive macroscopically. Values of S′, driven by perturbations in effective normal stress as given by Equation (7b), are expected to be small for slow deformation. Moreover, distinct-element modeling of fault-gouge deformation indicates that cushioning of large particles by adjacent smaller particles results in short, discontinuous force chains that keep stress fluctuations small (Reference Morgan and BoettcherMorgan and Boettcher, 1999). A reasonable assumption, therefore, is that the value of S′ is a small fraction of the mean shear strength of the till layer. Values of T are likely to be small since individual force chains remain intact for only small increments of strain (Reference Morgan and BoettcherMorgan and Boettcher, 1999).

Given that values of δ, S′ and T are not well known, we now explore whether displacement profiles, calculated with Equation (17), can be fitted successfully to those from Breidamerkurjökull (Fig. 3a) by adjusting the values of these parameters through reasonable ranges. The value of S′ controls the depth to which deformation will extend after sufficient time has elapsed for perturbations in effective normal stress to occur at all depths in the bed (Equation (9)). Thus, the value of S′ is adjusted to match the observed depth of deformation. Then T is adjusted until the total displacement of the uppermost till slice agrees with the observed displacement there. The value of δ is taken to be 0.01 m.

Figure 3b depicts displacement in the bed indicated by the model if this fitting procedure is used. We assume that limiting equilibrium was approached for one period each day at Breidamerkurjökull, consistent with the roughly diurnal variation in pore-water pressure measured there (Reference Boulton and DobbieBoulton and Dobbie, 1998, fig. 11), and that during each period a perturbation in effective normal stress occurred once at every depth in the bed. Fitted values of S’ and T are 3.1 kPa and 0.16 s, respectively. Somewhat larger values of T are obtained if a larger value of δ is chosen, but since displacements vary with T 2, T cannot be much larger without producing displacements that are unrealistically large. Conversely, the fitted value of T is smaller if a smaller value of δ is used.

The model results agree with the observed convex-upward displacement profile. This displacement pattern is independent of the fitted values for S′, T and δ and is therefore a fundamental feature of the model. The profile is convex upward because friction along slip surfaces increases with depth in the bed more rapidly than the downslope component of the weight per unit area of till and ice. Thus, a given perturbation in effective normal stress applied to each depth in the bed causes larger force imbalances and relative displacements between till slices near the glacier sole than at depth.

Discussion

The model presented here should be viewed as illustrative but not predictive. It is based on a fundamental principle, conservation of momentum, and a well-established empirical rule, Coulomb friction. It illustrates that the observed distribution of strain at Breidamerkurjökull arises naturally from these principles if combined with the fact that brief, local departures from static equilibrium instigate motion in the bed. On the other hand, the model is not a useful predictive tool at this time. It contains several parameters that are difficult to measure or estimate and that we have adjusted to bring the model results into quantitative agreement with observations. In addition, as discussed hereinafter, certain subglacial complexities are neglected.

An important tenet of the model is that over brief intervals and at specific depths in the bed, driving stresses exceed resisting stresses as force chains rearrange and local normal stresses fluctuate. Although such departures from quasi-static stress equilibrium are usually not considered in glaciology, they drive the slip events that cause the slow deformation of granular materials. Indeed, this concept is essential in distinct-element models of slow granular flow, in which net forces and torques on individual particles are computed and used to calculate particle displacements and rotations from Newton’s second law (e.g. Reference Cundall, Drescher, Strack, Vermeer and LugerCundall and others, 1982; Reference Morgan and BoettcherMorgan and Boettcher, 1999).

An idealization of the model is that perturbations in effective normal stress are assumed to occur at discrete times with equal frequency and magnitude at all depths. A more realistic model, which would require a numerical approach, would apply perturbations randomly with depth and time and a distribution of perturbation magnitudes consistent with that observed experimentally (e.g. Reference Veje, Howell, Behringer, Schollann, Luding, Herrmann, Herrmann, Houi and LudingVeje and others, 1998). There is no obvious reason to believe, however, that the displacement profile in that case would differ significantly from the concave-down profile obtained with our simpler model. If there were independent evidence that the character of force chains varied systematically with depth, causing associated variations in the magnitude or frequency of stress fluctuations, a more complicated model would be justified. Although such a situation is possible and might be caused, for example, by a progressive change in grain-size distribution with depth, it is probably not typical.

Another neglected complexity is that during a shearing episode, till porosity may vary, contrary to the constant bulk density assumed in the model. Local variations in till porosity are likely as force chains form and fail causing particle rearrangement. In addition to affecting the local friction angle of the till, a change in porosity may cause a change in pore-water pressure, owing to the generally small hydraulic diffusivity of till. If, for example, porosity increases during slip, pore-water pressure may be reduced, thereby strengthening slip surfaces and reducing the displacement expected from a given perturbation in effective normal stress. Reference Iverson, Hooyer and BakerIverson and others (1998) suggested that this process, called dilatant hardening (e.g. Reference ReynoldsReynolds, 1885), might distribute strain in till due to strengthening of active slip surfaces and their consequent migration. Although this process may well occur and help distribute deformation, our model illustrates that even in the absence of this effect, distributed deformation will occur in a Coulomb till.

Reference Tulaczyk, Mickelson and AttigTulaczyk (1999) and Reference Tulaczyk, Kamb and EngelhardtTulaczyk and others (2000) emphasized that water-pressure fluctuations at the ice–till interface and consequent pore-pressure diffusion into the bed will help distribute strain in a basal till layer. Pore-water pressure variations will decrease in amplitude with depth and arrive later at depth than at shallower levels. These pore-pressure variations cause temporal variations in effective normal stress with depth that are hypothesized to drive vertical migration of bed-parallel shear zones in the bed. Like dilatant hardening, this process may be important, but our model illustrates that it is not necessary to distribute deformation.

Conclusion

Despite its limitations, our model of subglacial till deformation demonstrates that displacement profiles like those observed at Breidamerkurjökull do not contradict the results of laboratory tests that indicate till behaves as a Coulomb material. This is an important conclusion since Coulomb behavior has widespread empirical support and has for many decades been used successfully to explain diverse phenomena associated with soil deformation. Coulomb behavior is, indeed, the appropriate null hypothesis for analyzing till deformation, so a viscoplastic till rheology should not be invoked if Coulomb models can account for observations. Our Coulomb model of the bed can account for observations at Breidamerkurjökull. Thus, fitting viscoplastic flow relations to displacement profiles measured in till there or in till beneath other glaciers may yield little insight.

Acknowledgements

This work was supported in part by the U.S. National Science Foundation, grant 0PP-9725360. We thank G. K. C. Clarke, D. Cohen, M. Truffer and J. S. Walder for helpful criticism and W. Harrison, in particular, for motivating important improvements.

Appendix

To obtain the bed-normal displacement profile for S 0 ≤ 0, Equation (13) is substituted into Equation (15) to give

(A1)

or

(A2)

Integrating yields the displacement profile:

(A3)

where

References

Aharoni, J. 1972. Lecture on mechanics. Oxford, Clarendon.Google Scholar
Allersma, H. G. B. 1982. Photo-elastic stress analysis and strains in simple shear. In Vermeer, P. A. and Luger, H. J., eds. Deformation and failure of granular materials. Brookfield, VT, A. A. Balkema, 345354.Google Scholar
Alley, R. B. 1989. Water-pressure coupling of sliding and bed deformation: II. Velocity–depth profiles. J. Glaciol., 35(119), 119129.Google Scholar
Alley, R. B., Blankenship, D. D., Bentley, C. R. and Rooney, S. T.. 1986. Deformation of till beneath Ice Stream B, West Antarctica. Nature, 322(6074), 5759.CrossRefGoogle Scholar
Beeler, N. M., Tullis, T. E., Blanpied, M. L. and Weeks, J. D.. 1996. Frictional behaviour of large displacement experimental faults. J. Geophys. Res., 101(B4), 86978715.CrossRefGoogle Scholar
Biegel, R. L., Sammis, C. G. and Dieterich, J.. 1989. The frictional properties of simulated gouge having a fractal particle-size distribution. J. Struct. Geol., 11(7), 827846.CrossRefGoogle Scholar
Blake, E. W., Fischer, U. H. and Clarke, G. K. C.. 1994. Direct measurement of sliding at the glacier bed. J. Glaciol., 40(136), 595599.CrossRefGoogle Scholar
Boulton, G. S. 1996. Theory of glacial erosion, transport and deposition as a consequence of subglacial sediment deformation. J. Glaciol, 42(140), 4362.Google Scholar
Boulton, G. S. and Dobbie, K. E.. 1998. Slow flow of granular aggregates: the deformation of sediments beneath glaciers. Philos. Trans. R. Soc. London, 356(1747), 27132745.CrossRefGoogle Scholar
Boulton, G. S. and Hindmarsh, R. C. A.. 1987. Sediment deformation beneath glaciers: rheology and geological consequences. J. Geophys. Res, 92(B9), 90599082.CrossRefGoogle Scholar
Boulton, G. S. and Jones, A. S.. 1979. Stability of temperate ice caps and ice sheets resting on beds of deformable sediment. J. Glaciol, 24(90), 2943.CrossRefGoogle Scholar
Clark, P. U., Licciardi, J. M., MacAyeal, D. R. and Jenson, J. W.. 1996. Numerical reconstruction of a soft-bedded Laurentide ice sheet during the last glacial maximum. Geology, 24(8), 679682.2.3.CO;2>CrossRefGoogle Scholar
Clark, P. U., Alley, R. B. and Pollard, D.. 1999. Northern Hemisphere ice-sheet influences on global climate change. Science, 286(5442), 11041111.Google Scholar
Clarke, G. K. C., Collins, S. G. and Thompson, D. E.. 1984. Flow, thermal structure, and subglacial conditions of a surge-type glacier. Can. J. Earth Sa., 21(2), 232240.CrossRefGoogle Scholar
Cundall, P. A. and Strack, O. D. L.. 1979. A discrete numerical model for granular assemblies. Géotechnique, 29(1), 4765.Google Scholar
Cundall, P. A., Drescher, A. and Strack, O. D. L.. 1982. Numerical experiments on granular assemblies: measurements and observations. In Vermeer, P. A. and Luger, H. J., eds. Deformation and failure of granular materials. Brookfield, VT, A. A. Balkema, 355370.Google Scholar
Dantu, P. 1957. Contribution à l’étude mécanique et géométrique des milieux pulverulents. In Glanville, W. H., ed. Proceedings of the Fourth International Conference on Soil Mechanics and Engineering Foundation. Vol. 1. London, Butterworth Scientific Publications, 144148.Google Scholar
Drescher, A. and de Josselin de Jong, G.. 1972. Photoelastic verification of a mechanical model for the flow of a granular material. J. Mech. Phys. Solids, 20(5), 337351.Google Scholar
Engelhardt, H. and Kamb, B.. 1998. Basal sliding of Ice Stream B, West Antarctica. J. Glaciol., 44(147), 223230.Google Scholar
Hindmarsh, R. C. A. 1998a. Drumlinization and drumlin-forming instabilities: viscous till mechanisms. J. Glaciol., 44(147), 293314.Google Scholar
Hindmarsh, R. C. A. 1998b. Ice-stream surface texture, sticky spots, waves and breathers: the coupled flow of ice, till and water. J. Glaciol., 44(148), 589614.Google Scholar
Hindmarsh, R. C. A. 1998c. The stability of a viscous till sheet coupled with ice flow, considered at wavelengths less than the ice thickness. J. Glaciol., 44(147), 285292.CrossRefGoogle Scholar
Hooke, R.LeB. 1998. Principles of glacier mechanics. Upper Saddle River, NJ, Prentice Hall.Google Scholar
Hooke, R. LeB., Hanson, B., Iverson, N. R., Jansson, P. and Fischer, U. H.. 1997. Rheology of till beneath Storglaciaren, Sweden. J. Glaciol., 43(143), 172179.Google Scholar
Howell, D., Behringer, R. P. and Veje, C.. 1999. Stress fluctuations in a 2D granular Couette experiment: a continuous transition. Phys. Rev. Lett., 82(96), 52415244.Google Scholar
Iverson, N. R., Hooyer, T. S. and Hooke, R. LeB.. 1996. A laboratory study of sediment deformation: stress heterogeneity and grain-size evolution. Ann. Glaciol., 22, 167175.CrossRefGoogle Scholar
Iverson, N. R., Baker, R. W. and Hooyer, T. S.. 1997. A ring-shear device for the study of till deformation: tests on tills with contrasting clay contents. Quat. Sci. Rev., 16(9), 10571066.CrossRefGoogle Scholar
Iverson, N. R., Hooyer, T. S. and Baker, R. W.. 1998. Ring-shear studies of till deformation: Coulomb-plastic behavior and distributed strain in glacier beds. J. Glaciol., 44(148), 634642.Google Scholar
Jenson, J. W., MacAyeal, D. R., Clark, P. U., Ho, C. L. and Vela, J. C.. 1996. Numerical modeling of subglacial sediment deformation: implications for the behavior of the Lake Michigan lobe, Laurentide ice sheet. J Geophys. Res., 101(B4), 87178728.CrossRefGoogle Scholar
Kamb, B. 1991. Rheological nonlinearity and flow instability in the deforming bed mechanism of ice stream motion. J. Geophys. Res., 96(B10), 16,58516,595.CrossRefGoogle Scholar
Licciardi, J. M., Clark, P. U., Jenson, J. W. and MacAyeal, D. R.. 1998. Deglaciation of a soft-bedded Laurentide ice sheet. Quat. Sci. Rev., 17(4–5), 427448.Google Scholar
Liu, C. H. and 6 others. 1995. Force fluctuations in bead packs. Science, 269(5223), 513515.Google Scholar
Logan, J. M., Dengo, C. A., Higgs, N. G. and Wang, Z. Z.. 1992. Fabrics of experimental fault zones: their development and relationship to mechanical behaviour. In Evans, B. and Wong, T.-F., eds. Fault mechanics and transport properties of rocks. London, Academic Press, 3367.Google Scholar
Mandl, G., de Jong, L. N. J. and Maltha, A.. 1977. Shear zones in granular material — an experimental study of their structure and mechanical genesis. Rock Mechanics, 9(2–3), 95144.CrossRefGoogle Scholar
Misra, A. 1997. Micromechanical parameters of particle assemblies: experiments and numerical simulations. In Chang, C. A., Babic, M., Liang, R. Y. and Misra, A., eds. Mechanics of deformation and flow of particulate materials. NewYork, American Societyof Civil Engineers, 174188.Google Scholar
Misra, A. 1998. Particle kinematics in sheared rod assemblies: experimental observations. In Herrmann, H. J., Houi, J. P. and Luding, S., eds. Physics of drygranularmedia. Dordrecht, etc., KluwerAcademic, 261266.Google Scholar
Mitchell, J. K. 1993. Fundamentals of soil behaviour. Second edition. New York, John Wiley and Sons Inc.Google Scholar
Morgan, J. K. and Boettcher, M. S.. 1999. Numerical simulations of granular shear zones using the distinct-element method. I. Shear zone kinematics and the micromechanics of localization. J. Geophys. Res., 104(B2), 27212732.Google Scholar
Ng, F. S. L. 2000. Coupled ice–till deformation near subglacial channels and cavities. J. Glaciol., 46(155), 580598.Google Scholar
Painter, B., Tennakoon, S. and Behringer, R. P.. 1998. Collisions and fluctuations for granular materials. In Herrmann, H. J., Houi, J. P. and Luding, S., eds. Physics of dry granular media. Dordrecht, etc., Kluwer Academic, 217228.Google Scholar
Reynolds, O. 1885. On the dilatancy of media composed of rigid particles in contact, with experimental illustrations. Philos. Mag., 20(127), 469481.Google Scholar
Skempton, A. W. 1985. Residual strength of clays in landslides, folded strata and the laboratory. Geotechnique, 35(1), 318.Google Scholar
Terzaghi, K., Peck, R. B. and Mesri, G.. 1996. Soil mechanics in engineering practice. Third edition. New York, etc., John Wiley and Sons.Google Scholar
Thorsteinsson, T. and Raymond, C. F.. 2000. Sliding versus till deformation in the fast motion of an ice stream over a viscous till. J. Glaciol., 46(155), 633640.CrossRefGoogle Scholar
Tika, T. E., Vaughan, P. R. and Lemos, L. J.. 1996. Fast shearing of pre-existing shear zones in soil. Géotechnique, 46(2), 197233.Google Scholar
Truffer, M., Harrison, W. D. and Echelmeyer, K. A.. 2000. Glacier motion dominated by processes deep in underlying till. J Glaciol., 46(153), 213221.CrossRefGoogle Scholar
Tulaczyk, S. 1999. Ice sliding over weak, fine-grained tills: dependence of ice–till interactions on till granulometry. In Mickelson, D. M. and Attig, J. W., eds. Glacial processes: past and present. Boulder, CO, Geological Society of America, 159177. (Special Paper 337.)Google Scholar
Tulaczyk, S. M., Kamb, B. and Engelhardt, H. F.. 2000. Basal mechanics of Ice Stream B, West Antarctica. I. Till mechanics. J. Geophys. Res., 105(B1), 463481.Google Scholar
Veje, C. T., Howell, D. W., Behringer, R. P., Schollann, S., Luding, S. and Herrmann, H. J.. 1998. Fluctuations and flow for granular shearing: results from experiment and simulation. In Herrmann, H. J., Houi, J. P. and Luding, S., eds. Physics of dry granular media. Dordrecht, etc., KluwerAcademic, 237242.Google Scholar
Walder, J. S. and Fowler, A.. 1994. Channelized subglacial drainage over a deformable bed. J. Glaciol., 40(134), 315.Google Scholar
Williams, J. and Rege, N.. 1997. Granular vortices and shear band formation. In Chang, C. A., Babic, M., Liang, R. Y. and Misra, A., eds. Mechanics of deformation and flow of particulate materials. New York, American Society of Civil Engineers, 6276.Google Scholar
Figure 0

Fig. 1. (a) Force chains in 2900 photoelastic (birefringent under strain) disks sheared within a Couette apparatus. The outside diameter of the sample chamber is 0.38 m. Disks are 7 and 9 mm in diameter. Lighter disks are under larger forces. The rate of shearing at the specimen center line is 2 mm s−1 (adapted from Howell and others, 1999). (b) G2, an optical measure of the force on disks proportional to the number of birefringent bands they display, as a function of displacement during an experiment like that shown in (a). Over the measured range of G2, its value is nearly linearly proportional to the force on disks. The value plotted is an average over one-tenth of the samplers circumference (adapted from Howell and others, 1999). (c) Local stresses measured normal to the shearing direction at the upper boundary of a till specimen during deformation in a ring- shear device. Forces are measured over circular platens with diameters of 18 and 36 mm. The rate of shearing is 0.01 mm s−1 (adapted from Iverson and others, 1997).

Figure 1

Fig 2. Model bed geometry.

Figure 2

Fig. 3. (a) Daily displacement of anchors in the bed of Breidamerkurjökull over a 17 day period, as reconstructed by Boulton and Dobbie (1998, fig 13). (b) Daily displacement indicated by the model over 17 days, assuming that limiting equilibrium was approached each day during one period of sufficient duration for a perturbation in effective normal stress to occur once at each depth in the bed. Parameter values used in the model are listed in Table 1.

Figure 3

Table 1. Model parameters used to generate Figure 3b