Hostname: page-component-78c5997874-v9fdk Total loading time: 0 Render date: 2024-10-31T23:34:23.058Z Has data issue: false hasContentIssue false

On the pulling power of ice streams

Published online by Cambridge University Press:  20 January 2017

T. Hughes*
Affiliation:
Department of Geological Sciences and Institute for Quaternary Studies, University of Maine, Orono, Maine 04469-0110, U.S.A.
Rights & Permissions [Opens in a new window]

Abstract

Gravity wants to pull an ice sheet to the center of the Earth, but cannot because the Earth’s crust is in the way, so ice is pushed out sideways instead. Or is it? The ice sheet “sees” nothing preventing it from spreading out except air, which is much less massive than ice. Therefore, does not ice rush forward to fill this relative vacuum; does not the relative vacuum suck ice into it, because Nature abhors a vacuum? If so, the ice sheet is not only pulled downward by gravity, it is also pulled outward by the relative vacuum. This pulling outward will be most rapid where the ice sheet encounters least resistance. The least resistance exists along the bed of ice streams, where ice-bed coupling is reduced by a basal water layer, especially if the ice stream becomes afloat and the floating part is relatively unconfined around its perimeter and unpinned to the sea floor. Ice streams are therefore fast currents of ice that develop near the margins of an ice sheet where these conditions exist. Because of these conditions, ice streams pull ice out of ice sheets and have pulling power equal to the longitudinal gravitational pulling force multiplied by the ice-stream velocity. These boundary conditions beneath and beyond ice streams can be quantified by a basal buoyancy factor that provides a life-cycle classification of ice streams into inception, growth, mature, declining and terminal stages, during which ice streams disintegrate the ice sheet. Surface profiles of ice streams are diagnostic of the stage in a life cycle and, hence, of the vitality of the ice sheet.

Type
Research Article
Copyright
Copyright © International Glaciological Society 1992

Introduction

The best definition of an ice stream still may be the one Henri Bader gave three decades ago: an “ice stream” is something akin to a mountain glacier consisting of a broad accumulation basin and a narrower valley glacier; but a mountain glacier is laterally hemmed in by rock slopes, while the ice stream is contained by slower-moving surrounding ice. The edges of the ice stream arc often crevassed, and the surface tends to be concave as the ice is “funnelled” down. Many of-the large outlet glaciers in Greenland, particularly in the south, are the narrow outlets of large ice streams which reach back many scores of miles into the ice sheet. The largest ones flowing into Disko Bay have corresponding depressions in the glacier floor. An interesting question is whether the ice stream makes the depression or vice versa; and an interesting hypothesis is that the ice stream, once started, is self-perpetuating because its ice mass, warmed up by heat of internal friction, has a lower viscosity than the surrounding ice (Reference BaderBader, 1961).

Insights into these ideas may be found in the pulling power of ice streams. As defined here, pulling power is postulated as the ability of an ice stream to reach deep into an ice sheet and pull out the ice. Formally, pulling power Px at some position on an ice stream is the product of the pulling force Fx and ice velocity u x at that point, both measured horizontally along direction x of flow:

(1)

where u x is averaged through an ice column of height h in a flow band of width w. For steady-state flow, ux is the equilibrium mass-balance velocity u.

Analyzing the pulling power of ice streams involved both old and new concepts. These concepts will now be identified.

Grcuntationai force, F z = Fg. The gravitational force is the product of the mass of ice (a scalar) multiplied by gravity acceleration (a vector directed toward Earth’s center), so F; is a vector vertically downward.

Pulling force., Fx = FP. The pulling force is the product of lithostatic pressure in ice (a scalar) multiplied by the vertical area across which the pressure is relieved (a vector whose length is proportional to the area and whose direction is away from the ice mass), so Fp is a vector directed horizontally downstream.

Pushing force, Fx = F B. The pushing force is the product of lithostatic pressure in ice (a scalar) multiplied by the vertical area across which the pressure is exerted (a vector whose length is proportional to the area and whose direction is toward the ice mass), so F P is a vector directed horizontally downstream.

Braking force, Fx = FB. Braking forces are either horizontal traction forces or horizontal kinematic forces. A traction braking force is the product of shear stresses σxy and σxz acting against side and basal surfaces of an ice stream, where shear traction exists, and the areas of these traction surfaces. A kinematic braking force is the product of longitudinal deviator stress σ′xx that causes kinematic straining, and the cross-sectional area affected by kinematic straining. The analyses of sheet flow by Orowan (British Glaciological Society, 1949) and of shelf flow by Wecrtman (1957a) were classic examples of using traction and kinematic horizontal braking forces, respectively.

Horizontal force balance, Fp-FB = 0. The horizontal force balance requires dynamic equilibrium between horizontal force Fp, whether pulling or pushing, both of which require gravity acceleration (which appears in the lithostatic pressure), and braking forces FR, whether traction forces or kinematic forces, both of which require horizontal ice velocity (which causes straining).

Pulling stress, The pulling stress is an axial deviator stress in the direction of horizontal ice flow. It represents a deficiency of lithostatic pressure in the horizontal x-direction of ice flow, creating a relative vacuum in that direction that is filled by inrushing ice. Hence, strain rate is tensile and proportional to pulling stress σ′xx, where σσxx is the average difference between the reduced lithostatic stress xx and the full lithostatic stress σxx, both of which are compressive, so that σ′xx is tensile.

Pulling power, . pulling power is the product of the horizontal pulling force Fp and mean ice velocity . Therefore, pulling power requires both gravity acceleration (which causes the pulling force) and ice velocity (which causes the braking force). Pulling power is a measure of the vitality of an ice stream, which changes during the life cycle of an ice stream as gravitational potential energy (represented by ice elevation) is converted into kinetic energy (represented by ice velocity) as stream flow moves the ice mass horizontally to lower elevations.

Basal buoyancy factor, The basal buoyancy factor quantifies ice-bed coupling and ice-shelf buttressing, as the two sources of traction and kinematic braking forces that resist the pulling force. Ice-bed coupling is caused by the fraction of basal lithostatic pressure in ice of height h and density ρ1 that is supported by the bed, instead of by the basal hydrostatic pressure for water of depth d and density ρw Full ice-bed coupling (ø= 0) occurs when d = 0. No ice-bed coupling (ø = øG) occurs when d = (ρ1w)h. Full basal buoyancy (ø = 1) occurs when an ice stream is not coupled to the bed and not buttressed by an ice shelf øG = 1). No basal buoyancy = 0) occurs when an ice stream is either fully coupled to the bed (d = 0), or fully buttressed by an ice shelf grounded around its entire perimeter (øG = 0), or both. The pulling force is maximized for full basal buoyancy and is zero for no basal buoyancy. Ice bed coupling produces traction forces and ice-shelf buttressing produces kinematic forces, both of which constitute the braking force that balances the pulling force.

Background for the concept of pulling power will be presented, followed by an elaboration on the distinction between pushing and pulling, and an elaboration upon ice-bed coupling and ice-shelf buttressing. Then the horizontal force balance is presented for sheet flow, shelf flow and stream flow. Pulling power is then related to four ice-stream surface profiles, (1) with minimal basal buoyancy, (2) with basal buoyancy continuously decreasing upstream, (3) with equilibrium basal buoyancy, and (4) with maximum basal buoyancy. These conditions of basal buoyancy become the basis for postulating a life-cycle classification for ice streams consisting of inception, growth, mature, declining and terminal stages, with the possiblity for rejuvenation at any stage after inception. A way to model the ice-stream life cycle is presented, and driving the life cycle by pulling power is discussed. Next, pulling power is postulated as the mechanism for disintegrating marine ice sheets by its role in down-drawing ice sheets, discerping ice ridges between ice streams, controling the erosive power of ice streams, regulating ice-shelf buttressing, widening ice streams and causing ice streams to surge. Finally, ways to initiate stream flow are discussed, and results for marine ice streams are generalized to include terrestrial ice streams.

Background

Work related to the pulling power of ice streams dates at least from the old dispute between Reference LliboutryLliboutry (1958) and Reference NyeNye (1958) as to why transverse crevasses open when a glacier enters a narrow valley. Lliboutry (Reference Lliboutry1958, fig. 4) maintained “To open a crevasse we need a strain, not a stress.” His analog was a plastic slab being stretched between rolls, such that transverse compression by the rolls causes longitudinal extension that opens pre-existing transverse cracks. Nye (Reference Nye1958, fig. 1) countered, “I say a stress.” His analog was a plastic slab cut into blocks and compressed between frictionless parallel plates, the cuts being the pre-existing transverse cracks. Transverse compressive stress exerted by the plates causes longitudinal extensional strain in the slab, but docs not open the cracks. “If, however, we apply a (longitudinal) tensile stress … the whole thing comes apart.” Reference NyeNye (1958) had proposed pulling power as a disintegration mechanism. His analog was for a valley glacier, but it has significance for ice streams in view of the observation by Reference BaderBader (1961) that “An ice stream is something akin to a mountain glacier (that enters) a narrow valley …”

Another indication that ice streams pull ice out of ice sheets was the analysis by Reference WeertmanWeertman (1963) of outlet glaciers draining an ice sheet fringed by mountain ranges. Although his analysis was for convex surface profiles, if it is applied to the concave surface of a marine ice stream on a horizontal bed, the longitudinal strain rate is a maximum across the surface-inflection line at the head of the ice stream and a minimum across the basal grounding line at the foot of the ice stream, so that basal frictional heat is conducted to the surface faster at the head than at the foot of the ice stream. Chances for a frozen bed are therefore greater at the head, with stream-flow sliding on a thawed bed pulling sheet flow anchored to a frozen bed.

The ice-stream analysis by Reference WeertmanWeertman (1974) involves pulling. He considered an ice stream as the transition zone between an ice sheet and an ice shelf in which the down-slope increase in a longitudinal tension stress allows a corresponding reduction of basal shear stress, and thereby produces a concave profile. Hence, concave “necking” of the ice stream is produced by a longitudinal tensile stress, just as arc the transverse crevasses analyzed by Reference NyeNye (1958). A longitudinal tensile stress requires a longitudinal pulling force.

The concave surface profile of West Antarctic ice streams can be interpreted as a “necking” phenomenon that was causing the maximum-slope surface-inflection line at the heads of ice streams to retreat, ultimately collapsing the West Antarctic ice sheet (Reference HughesHughes, 1973). Necking has been observed directly by Reference Whillans, Bolzan and ShabtaieWhillans and others (1987) in longitudinal velocity gradients at the head of Ice Stream B (Fig. 1). The conclusion that necking is caused by pulling also emerges in using slip-line theory to represent the transition from sheet flow to stream flow to shelf flow as analogous to plastic material (ice) being successively drawn through a die, compressed between plates and indented by a punch. Representing the sheet flow to stream-flow transition as being analogous to pushing ice through a die (Reference HughesHughes, 1977, Fig. 17) did not explain the great transverse crevasses arcing around the head of Byrd Glacier (Fig. 1) in the zone of converging flow. These crevasses would form if Byrd Glacier was pulling ice out of East Antarctica, and the analog was to ice being pulled through the die (Reference Stuiver, Denton, Hughes, Fas-took, Denton and HughesStuiver and others, 1981, Fig. 8-7). This observation was the basis for an early attempt theoretically to quantify the concave profile of ice streams (Reference Hughes, Denton and HughesHughes, 1981a), and to employ retreating ice streams as the vehicle for collapse of the West Antarctic ice sheet in a computer simulation (Reference Stuiver, Denton, Hughes, Fas-took, Denton and HughesStuiver and others, 1981, Fig. 8-7). It can be argued that stream flow also pulls flanking ice into the ice stream, so that the proper analog with a plastic material between parallel plates is that extending plastic flow actively pulls the plates together, rather than extending flow being a passive response to the plates being pushed together.

Fig. 1 Location map for selected Antarctic ice streams.

Pushing versus pulling

Figure 2 illustrates the distinction between pushing and pulling. Figure 2a shows horizontal spreading of an ice sheet as resulting from vertical lowering caused by a downward pulling force Fz equal to the mass of the ice sheet multiplied by gravity acceleration g. When the ice sheet is in steady-state equilibrium, the surface profile is unchanged because accumulation rates match the lowering rate and ablation rates match the spreading rate. Vertical pulling by gravity causes horizontal pushing. As seen in Figure 2b, if the ice sheet beyond distance x from the ice margin is removed, it can be replaced by a horizontal pushing force which is resisted by basal traction force t0wx over distance ι in a flow band of constant width w, where ρ1 is ice density, h is ice height at position x and T0 is basal shear stress over distance x. Note that Fx is directed toward the remaining ice mass, so that by definition Fx = Fp is a pushing force. If To is constant, equating the pushing force with the traction force gives This is the convex surface profile of an ice sheet predicted by plasticity theory, in which T 0 is the yield stress.

Fig. 2 Producingstream /loti) from sheet flow, (a) In sheet flow, gravitational pulling force FZ deforms the ice sheet from the solid profile to the clashed profile, (b) Inner dashed part of the ice sheet can be replaced by horizontal pushing force Fx. (c) Outer dashed part of the ice sheet can be replaced by horizontal pulling force Fx. (dj Horizontal pulling converts the connex surface of sheet flouj into the concave surface of stream flow, with stream flow becoming a floating ice shelf in terminus 1 or a grounded ice lobe in terminus 2.

Horizontal pulling can be induced by removing the ice sheet over distance x, so that becomes a pulling force resisted by basal traction force over distance L-x where L is the length of the flow band from the ice divide to the margin, as shown in Figure 2c. Note that Fx is directed away from the remaining ice mass, so that by definition Fx = Fp is a pulling force. Pulling extends only to the ice divide because To reverses sign across the ice divide. Equating the forces for constant T 0 gives as the distance over which Fx is exerted for a given h and To in Figure 2c. With Fx a pulling force instead of a pushing force, the ice surface over distance L-x is quickly pulled down by Fx as ice is pulled forward by Fx. A similar result is obtained by not removing the ice sheet over distance x but, instead, having TO be large over L-x and small over x, which produces the concave surface of stream flow in Figure 2d. An ice stream can end as a floating ice shelf, beneath which T 0 is zero (shown as terminal region 1), or as a grounded ice lobe, beneath which To increases to the ice margin (shown as terminal region 2). A marine ice stream supplies an ice shelf. A terrestrial ice stream supplies an ice lobe. Marine ice streams are the subject of this investigation because predictions can be tested by field studies on the numerous present-day marine ice streams. In Figure 2d, pulling force Fp at distance x from the ice-stream grounding line or the ice-lobe terminus will be net horizontal force where ρw is water density and d is the depth of water that would produce the hydrostatic pressure of basal meltwater at x.

As presented in Figure 2, horizontal force Fx is a pushing force for sheet flow and a pulling force for stream flow, both produced by gravity, and both inducing traction forces that resist gravitational motion. Hence, Fx will be used to represent gravitational forcing, with Fx = Fp causing pushing or pulling, and Fn will represent braking forces caused by stresses that resist pushing or pulling motion.

Ice-bed coupling and ice-shelf buttressing

As presented in Figure 2, stream flow is generated from sheet flow by reducing ice bed coupling toward the terminus of a flow band. If that is the case, marine ice streams are the major dynamic links connecting a marine ice sheet with its floating ice shelves, and stream flow is transitional between sheet flow and shelf flow. However, stream flow cannot develop from sheet flow if reduced ice-bed coupling is completely offset by increased ice-shelf buttressing. To illustrate this, consider an ice cube in a pan. If water is added to the pan, an ice cube of height h will begin to float in water of depth d when the lithostatic pressure of ice and the hydrostatic pressure of water give the same overburden pressure σ0 on the bottom of the pan:

(2)

where ρ1 is ice density, σw is water density and g is gravity acceleration. Ice-bed coupling is nil when Equation (2) is satisfied. Now imagine that the ice cube is a large tabular iceberg and the pan is a larger lake. Elimination of ice-bed coupling allows pulling force Fp to spread the iceberg radially, in the manner analyzed by Reference WeertmanWeertman (1957a), until the iceberg grounds along the shoreline of the lake, thereby producing a braking force FB that completely cancels FP. Hence, ice-bed coupling is removed in the center of the lake by a gain of basal buoyancy, and ice-shelf buttressing is attained around the perimeter of the lake by a loss of basal buoyancy.

The above considerations suggest that the pulling power of a marine ice stream is controled by basal buoyancy. Pulling power is increased when ice-bed coupling beneath an ice stream is decreased by basal buoyancy. Pulling power is decreased when ice-shelf buttressing beyond an ice stream is increased by a loss of basal buoyancy around the grounded parts of the ice shelf. Both ice-bed coupling and ice-shelf buttressing can therefore be quantified by a basal buoyancy factor ø that is defined as follows:

(3)

where basal water would rise to a depth d in ice of height h in imaginary temperate boreholes along length Ls of an ice stream, as shown in Figure 3, and øc is ø at the, ice-stream grounding line of a floating ice shelf. In determining ø from Equation (3), variations in d are controled by variations in ice-bed coupling, and variations in ø are controled by variations in ice-shelf buttressing.

Fig. 3 Pulling forcea defined us net horizontal lithostatic and hydrostatic forces acting in a marine flow band of an ice sheet. Horizontal lithostatic and hydrostatic forces pressing against imaginary ice columns are shown for sheet flow (right), stream flow (center) and shelf flow (left), with the. flow band having floating length Lp, grounded length L and streaming length Ls- Averaged for the ice column, basal water rises to depth d in imag-iiiary temperate boreholes drilled- through the flow band. Four cases analyzed in the text for the ice stream are no buoyancy along Ls (case I), full buoyancy at the basal grounding line decreasing to no buoyancy at the surface-inflection line (case II), constant buoyancy along Ls (case III) and full buoyancy along Ls (case IV).

Since basal water can exist as a thin film at the ice-rock interface, can flow in a network of subglacial channels, can be ponded in subglacial lakes and can occupy interstitial pores in subglacial till or sediments, depth d of basal water standing in imaginary temperate boreholes, as employed in Equation (3), is a statistical averaging of all these forms of basal water across width w of the ice stream at a distance x from its grounding line. In particular, it will be assumed that d = 0 where basal water is a thin film. This assumption is reached by considering a subglacial lake in which bedrock hilltops penetrate temperate basal ice. A thin water film covers the hilltops. If the lake is now drained, the bedrock hills would act like pillars supporting the ceiling of basal ice. Therefore, even though a water film exists on the hilltops, the overburden of ice is supported by bedrock, not basal buoyancy, because bedrock is continuous whereas the water film exists only at the ice-rock interface on hilltops. This fact is incorporated into Equation (3) by setting ci = 0 in imaginary boreholes drilled down to the hilltops, even though the hydrostatic pressure of water in both the film and the lake equals the same lithostatic overburden pressure of ice before the lake was drained, and lake water would rise to a depth d = (ρ1/ ρW)h in a borehole drilled above the lake. Another reason for letting d = 0 when the water layer is a film over rugged bedrock is that d has meaning for pulling power only insofar as ice-bed coupling resists stream flow. Hence, d = 0 for all undrowned bedrock projections because these resist stream flow by retaining ice-bed coupling that retards horizontal ice motion. This is the essence of the Reference WeertmanWeertman (1957b) sliding theory, in which basal water pressure is not involved, only bed roughness which provides traction that retards sliding. This understanding of d can be generalized in Equation (3) by specifying that d = 0 gives ø = 0 for sheet flow with no basal buoyancy, d = (ρ\ρW)h gives ø= øG for shelf flow with full basal buoyancy, and intermediate values of d give 0 < ø < øG for stream flow with partial basal buoyancy. Ice-shelf buttressing of an ice stream is determined by øG in Equation (3), such that øG 0 when the ice shelf is grounded around its entire perimeter, øG 1 when the ice shelf is grounded only along the ice-stream grounding line, and 0 < øG < 1 when the ice shelf occupies an embayment or is pinned by islands and shoals. In Figure 1, examples of øG = 0 are ice above subglacial lakes near the interior ice dome of Wilkes Land (Reference Oswald and RobinOswald and Robin, 1973), and the “pseudo ice shelves” of Ice Stream C (Reference Robin, Swithinbank and SmithRobin and others, 1970); examples of øG= 1 are the floating ice tongues of Thwaites, David, Ninnis and Mcrtz Glaciers; and examples of 0 < øG < 1 are the Ross and Ronne Ice Shelves.

Refer to Figure 3, for which d changes along an ice stream of length Ls and øG = 1 for an ice shelf of length LF. Over length L-Ls, sheet flow with ø = 0 exists when the bed is frozen, or thawed, but with only a basal water film. Over length Ls, stream flow with ø = 0 exists when basal water is only thick enough to drown bedrock projections up to the controling obstacle size in the Reference WeertmanWeertman (1975b) theory of basal sliding, and ø = 1 when basal water saturates basal till or sediments so they are unable to support a basal shear stress; otherwise 0 < ø < 1 in the ice stream, depending on how d varies along Ls, where d is statistically averaged for basal water conditions across a given width w of the ice stream.

Ice-bed coupling and ice-shelf buttressing are quantified by ø to provide a means for assessing the pulling power of ice streams using Equation (1). The plan to be followed consists of:

  • 1. Presenting a force balance for stream flow that reduces to force balances for sheet flow when ø = 0 and shelf flow when ø = 1.

  • 2. Deriving ice-stream profiles for a range of ø values and computing pulling power along these profiles for steady-state equilibrium.

  • 3. Classifying ice streams in a hypothetical scheme in which the pulling-power curve is diagnostic of the position of an ice stream in its life cycle.

  • 4. Applying the concept of pulling power to disintegration of marine ice sheets, notably in West Antarctica.

The goal of this plan is to provide numerical models of ice sheets with a parameter ø that quantifies ice-bed coupling and ice-shelf buttressing for the ice streams that drain the ice sheet. Rapid collapse of marine ice sheets can then be simulated by allowing ø to represent increases in ice-stream pulling power resulting from uncoupling and unbuttressing.

Force Balances Along a Flow Band

Figure 3 shows an ice sheet flow band of constant width on a horizontal bed. The flow band begins as sheet flow (right) and ends as shelf flow (left), with stream flow (center), being a transition from sheet flow to shelf flow. Constant width and a horizontal bed are specified so the transitional nature of stream flow can be analyzed in its simplest form, without complications introduced by laterally converging or diverging flow and variable bed topography. The origin of rectilinear coordinates is the ice-shelf grounding line, with x horizontal and positive toward the ice divide and z vertical and positive upward. The flow band has floating length Lp, and an ice stream of length Ls over grounded length L, with ha being ice height at the grounding line separating shelf flow from stream flow and hG being ice height at the inflection line separating stream flow from sheet flow.

Force balances are computed for ice columns in the sheet-flow, stream-flow and shelf-flow parts of the flow band. Assumptions are (1) dynamic equilibrium exists, (2) longitudinal gradients in stretching stress are small, and (3) the T-term can be neglected (Reference PatersonPaterson, 1981, p. 164). Assumption (1) is justified if pulling force Fy is added to braking force FB such that Fp + FB = 0. Assumptions (2) and (3) arc reasonable for a flow band of constant width on a horizontal bed (Reference MuszynskiMuszynski, 1987; Reference Muszynski and BirchfieldMuszynski and Birchfield, 1987; Reference Van der Veen, Van der Veen and OcrlemansVan der Veen, 1987). The pulling force in all cases is merely the expression of Newton’s second law, Fp being gravity acceleration acting on the mass of the ice column, with Fp in the horizontal x-dircction being an average differential pressure exerted on faces of the ice column that are normal to the x-direction. The braking force in all cases arises from motion or deformation of the ice column, with FR consisting of the sum of stresses resisting the motion or deformation multiplied by the area of the ice column upon which these stresses act. None of this is new, but it will now be reviewed for sheet flow and shelf flow, before analyzing stream flow as being transitional between these two.

Force balance for sheet flow

Pushing force Fp for an ice sheet on a horizontal bed is obtained from the horizontal lithostatic forces shown in Figure 4 (right). A vertical ice column has down-slope height h and upslopc height h + Δh along incremental length Δx over which flow-band width w is constant. The pushing force per unit width is the difference between the areas of force triangles 1 and 2, defined as the mean lithostatic pressures acting on the upslopc and downslope sides of the ice column, respectively. Therefore, for an ice-shect flow band of width w:

(4)

where the term containing (Δh)2 can be ignored. Note that Fp is a force gradient, because it increases with ice-thickness change Δh. Lithostatic pressure pushes against area wΔh to produce Fp.

Fig. 4 Horizontal forces on ice columns having Constant width and no side traction. Horizontal gravitational forces per unit width are the areas of triangles 1, 2 and 4, and rectangle 3. Traction force is basal shear stress To times the basal area of the ice column. Pulling force is tensile deviator stress 2′σxx times the transverse cross-sectional area on the landward side of the ice column. Following Reference Whillans, Van der Veen and OcrlemansWhillans (1987), the role of basal uiater pressure can be emphasized by distinguishing lithostatic pressure in triangles 1 and 2, and rectangle 3 from hydrostatic pressure in triangle 4.

Pushing gravitational motion induces creep controled by deviator shear stress σxz, where T0 = σxz at z = 0 is the basal shear stress. Since resistance to creep is provided by basal shear stress T 0 acting on basal area wΔx, and the braking force is:

(5)

where σxz is expressed in terms of the strain rate of ice in simple shear for the creep law

(6)

in which T is the effective shear stress, n is a viscoplastic exponent, A is an ice-hardness parameter, A averaged through h is Ā. and T = σxz= for simple shear.

When the ice column is in dynamic equilibrium, Fp = −FB and (Reference NyeNye, 1957):

(7)

where is often called the “driving stress” for an ice sheet. Actually, T 0 is a braking stress that resists ice motion. It takes on the role of a driving stress only when Fp is equated with FB, as only Fp expresses the product of ice mass and gravity acceleration.

Force balance for shelf flow

Pulling force Fp for an ice shelf having only longitudinal extension and no side traction is obtained from the horizontal lithostatic and hydrostatic forces shown in Figure 4 (left). A vertical ice column has constant height h and width w over incremental length Δx and floats at depth d in water. Ice on the side facing the calving front can be replaced by an equivalent mass of water without changing the nature of this boundary, whereas ice facing the grounding line cannot. Hence, horizontal forces on opposite sides arc unbalanced, with a net seaward pulling force. The pulling force per unit width is the difference between the areas of force triangles 1 and 4, defined as the mean lithostatic and hydrostatic pressures acting on the sides of the ice column facing the grounding line and the calving front, respectively, of the ice shelf. Therefore, for an ice-shelf flow band of width w:

(8)

where d = h(ρI/ ρ\W) satisfies the buoyancy requirement for floating ice, ø = 1 in Equation (3).

Pulling gravitational motion induces creep controled by deviator tensile stress σ′xx, where

and σxxzz is the deviation of longitudinal stress σxx from overburden lithostatic stress in the ice column. Both σxx and σzz- are zero at the ice surface, so resistance to creep is provided by average longitudinal tensile stress acting at distance h-z below the surface on transverse cross-sectional area wh, which is constant over length Δx. Braking force FB therefore results from σxx, which produces kinematic strain rate as follows:
(9)

where σ’xx is expressed in terms of the strain rate of ice in longitudinal tension for the creep law:

(10)

in which is constant through h over length Δx and T = σ′xx for constant w.

When the ice column is in dynamic equilibrium, FP = − FB and (Reference WeertmanWeertman, 1957a):

(11)

where for constant w and volume conservation, and is the basal lithostatic and hydrostatic pressure.

Force balance for stream flow

Pulling force Fp for an ice stream having longitudinal extension and both basal and side traction is obtained from the horizontal lithostatic and hydrostatic forces shown in Figure 4 (center). A vertical ice column has downslope height h and upslope height h + Δh along incremental length Δx, over which basal water would rise in temperate boreholes to an average depth d above the bed. The average d is based upon the fractions of the bed in basal area wΔx having no buoyancy (ø = 0) and full buoyancy (ø = øG). The pulling force per unit width is the area of force triangle 1 minus the combined areas of force triangle 2, force rectangle 3 and force triangle 4. Therefore, for an ice stream of width w:

(12)

where ø is introduced by using Equation (3). Note that Equation (12) for stream flow reduces to Equation (4) for sheet flow when ø = 0 and to Equation (8) for shelf flow wheø = 1 and Δh = 0 along length Δx.

In Equation (12), downslope pulling by the first term is partly offset by upslopc pushing by the second, third and fourth terms to give a net forward pulling. The first term is the mean lithostatic pressure exerted on column area w(h + Ali.) to produce the horizontal pulling force per unit width given by the area of triangle 1 in Figure 4. The second term is the mean lithostatic pressure above height d(ρW\ρI) exerted on column area w[h]-d(ρW/ρI)] above this height to produce the horizontal pushing force per unit width given by the area of triangle 2 in Figure 4. The third term is the lithostatic pressure at height d(ρW\ρI) exerted on column area wd(ρw/ρi) be low this height to produce the horizontal pushing force per unit width given by the area of rectangle 3 in Figure 4. The fourth term is the mean hydrostatic pressure below height d exerted on column area wd below this height to produce the horizontal pushing force per unit width given by the area of triangle 4 in Figure 4.

Notice that the horizontal hydrostatic force per unit width exerted by water, the area of triangle 4 in Figure 4, is shown only on the seaward side of ice columns for which d > 0 in Figure 3. This is because a marine ice stream and its floating ice-shelf extension interface with water beyond the seaward wide of the ice column, but not the landward side, where only lithostatic forces exist. Hence, any gradient ∂d/∂x along length Ls of stream flow, as shown in Figure 3, need not be considered in balancing horizontal forces on the ice columns.

That stream flow is transitional between sheet flow and shelf flow is seen in Figure 4. Moving d(ρW/ρI) downward increases Δh by causing triangle 2 to encroach upon rectangle 3 and triangle 4 until sheet flow occurs atd(ρW/ρI) = 0, reducing Equation (12) to Equation (4) at ø = 0. Moving d(ρW/ρI) upward decreases Δh by causing triangle 4 to encroach upon rectangle 3 and triangle 2 until shelf flow occurs at d(ρW/ρI) = h, reducing Equation (12) to Equation (8) at ø = 1 and Δh = 0.

Pulling gravitational motion induces creep controled by deviator tensile stress σ′xx along Δx and deviator shear stresses σxy and σxyts along w and h, respectively, with maximum values σxy = Ts rs for side traction and σxz = To for basal traction in the ice stream. Multiplying these deviator stresses by the areas against which they act gives braking force FR. Taking as the mean ice height along Δx and where σxx − σzz is the deviation of longitudinal stress σxx from lithostatic stress in the Reference RobinRobin (1967) analysis:

(13)

where is the change in σ′xx along σx and the term containing can be ignored.

When the ice column is in dynamic equilibrium,

and the following expression for Δh/Δx is obtained:
(14)

where Plgh and:

(15)

Equation 14) and (15) apply to stream flow, which begins with sheet flow at the head and ends as shelf flow at the foot of a marine ice stream. Stream flow becomes laminar sheet flow when in which case Equation (14) can be combined with Equation (6) to produce Equation (7). Stream flow becomes longitudinal shelf flow when

, in which case Equation (14) can be combined with Equation (10) to produce Equation (11).

The unknown variation of ø along length Ls of stream flow prevents a single simple expression relating Fp to strain rates, such as Equation (7) for sheet flow and Equation (11) for shelf flow. Nevertheless, the pulling power of an ice stream can be appreciated by comparing for sheet flow at the head of an ice stream with for shelf flow at the foot of the ice stream using these equations. In Equation (7), h = 2 km and Δh/Δx = 0.005 are typical at an ice-stream surface inflection line, where the maximum value of at Z= 0 exists for sheet flow. In Equation (11), 1 km is typical at an ice-stream basal grounding line, where the maximum value of exists for unconfined shelf flow. Taking pI = 0.92 Mg m-3, ρw = 1.02 Mg m-3, and n = 3, Equation (7) and (11) then give

= 15.6 if Ā is the same at both locations. Hence, for maximum shelf flow is an order of magnitude greater than for maximum sheet flow, and the ability of an ice stream to reach far into an ice sheet and pull out ice, its pulling power, is measured by the length Ls of stream flow and how effectively at the grounding line is transmitted along Ls.

The one-dimensional force balance considered here as pushing or pulling horizontal gravitational forces, which produce sheet, stream and shelf flow resisted by braking forces, has much in common with the two-dimensional force balance of driving forces and resistive forces analyzed by Reference Whillans, Van der Veen and OcrlemansWhillans (1987) and Reference Van der Veen and WhillansVan der Veen and Whillans (1989). For example, Reference Whillans, Van der Veen and OcrlemansWhillans (1987) began by stating “Glaciers flow in response to the interaction of gravity, which pulls the glacier forward, and resistive forces, which hold it back.” His Table 1 compares the conventional division of stresses into spherical and deviatoric components, on the one hand, with his lithostatic and resistive components, on the other hand, where his lithostatic force is gravitational and gradients in it produce his driving stress for glacial ice. In his Figure 2, Reference Whillans, Van der Veen and OcrlemansWhillans (1987) showed a nearly exponential decrease of his driving stress from the ice divide of West Antarctica to the calving front of the Ross Ice Shelf, along a flowline containing Ice Stream B. This same depiction of stream flow as being transitional from sheet flow to shelf flow is shown in Figure 3 as resulting from a general increase in basal buoyancy from the surface-inflection line to the basal grounding line of an ice stream. Specific changes in basal buoyancy, as measured by the change in d along Ls in Figure 3 (bottom), constitute four cases in a spectrum of ø variations that are possible for ice streams. These four cases will be analyzed in detail.

Table. 1 Variations of basal shear stress Towith distance L along flow bands of length L and constant width for sheet flow over frozen (T0 =TM) and thawed (T0 = TM) beds for constant surface accumulation and ablation rates a and b separated by an equilibrium line at distance x = E from the margin of the ice sheet (Reference HughesHughes, 1981b)*

Pulling Power and Surface Profiles

Pulling force Fp given by Equation (12) applies for stream flow but it does not describe the pulling force adequately because the ΔAh term depends on incremental length Δx along the ice stream, such that this term vanishes as Δx approaches zero, leaving only the ø term. Hence, the actual pulling force Fx is the limit of Fp as and is:

(16)

where pulling stress 2σ’xx, given by Equation (15), acts on cross-sectional area wh. The G term in the vertically integrated force balance (Patcrson, 1981, p. 164) is approximately the gradient of Fx. Hence, the pulling force can be understood in relation to vertically integrated momentum conservation. Mean ice velocity u can be computed to a good approximation by the following expression for conservation of ice flux and steady-state equilibrium:

(17)

where x = 0 at the grounding line of a marine ice stream, x = Ls at the head of the ice stream, where us is the mean ice velocity, and rates of accumulation and ablation along Ls are a and b, respectively. Entering Equation (16) and (17) into Equation (1) gives the pulling power:

(18)

Solutions of Equation (18) require knowing how h varies along h over length Ls of the ice stream. This requires knowing the ice-stream profile.

Surface profiles and the resulting pulling power will be computed for four kinds of marine ice streams:

  • Case I. Stream flow with minimal basal buoyancy, so that ø = 0 along Ls.

  • Case II. Stream flow with decreasing basal buoyancy, such that along Ls.

  • Case III. Stream flow with equilibrium basal buoyancy, such that along Ls.

  • Casc IV. Stream flow with maximum basal buoyancy, so that ø = øG along Ls.

Subscripts S and G refer to the surface-inflection line and to the basal grounding line, the respective beginning and end of the ice stream.

Cases I through IV belong to a spectrum of ø variations along length Ls of stream flow, with I and IV being end members of the spectrum, as shown in Figure 5. At this early stage in our understanding of ice streams, it is adventurous to speculate on bed conditions that may produce the ø spectrum in Figure 5, or other versions of a ø spectrum, but an attempt to depict possible bed conditions is offered in Figure 6. In all cases, bedrock is depicted as consisting of a series of ridges and troughs elongated along the axis of stream flow, such as mapped seismically beneath Rut ford Ice Stream (Reference Doake, Frolich, Mantripp, Smith and VaughanDoake and others, 1987) and Ice Streams A, Β and C (Reference Shabtaie and BentleyShabtaie and Bentley, 1988) in West Antarctica (see Fig. 1).

Fig. 5 A possible spectrum of ø variations along x for ice streams of length Ls. The four cases of basal buoyancy shown in Figure 3 (bottom) as average water depths d representative of average basal water pressure are located in this spectrum as shown.

Fig. 6 Four possible cases of basal water configurations shown as average d variations in Figure 3 (bottom) and average ø variations in Figure 5. Ice-stream beds are shown in plan view (top) and transverse cross-section (bottom), with thin basal water films for which ø = 0 shown in white, thick basal water layers for which ϕ=0 shown in black, and water-saturated till or sediment for which ø = øG shown as dotted. The vertical is exag-gerated.

  • Case I. shows the network of narrow water channels along trough bottoms, for which ø = øG in channels requires a downslope pressure-gradient current, and ø ø≈=0 in a water film between water channels, so that ø ≈ 0 averaged across width w all along length Ls of the ice stream. This situation exists when decreasing basal water (bottom freezing) and smoother bedrock topography occur in the direction of ice flow.

  • Case II. shows the water channels widening in the direction of ice flow, such that averaged across w, ø = 0 where stream flow begins at x = Ls and ø = øG where stream flow ends at x = 0. This widening can be caused if increasing basal water (bottom melting) and smoother bedrock topography occur in the direction of ice flow.

  • Case III. shows wide water channels separated by narrow ridge crests along the entire length Ls of the ice stream, such that ø ≈ (hG/h)øG averaged across w along Ls. This condition exists when basal water depths nearly match bedrock topographic variations in the direction of ice flow, and no net current exists in the water. It differs from case I only in that more basal water is present all along Ls, so that hydrostatic equilibrium is maintained in the water channels (channel water pressure where h > hG is the same as at h hG, SO there is no pressure-gradient current).

  • Case IV. shows a pad of water-soaked till or sediments between basal ice and bedrock, with deformation and permeability of the pad being such that ø = øG across w and along LQ (local basal water pressure is that required for the buoyancy of local ice thickness h, whereas in case III it is that required for the buoyancy of ice thickness fiQ at the grounding line). This situation can exist when permafrost having a high ice fraction becomes thawed, and it requires thawing at x = Ls to produce new melt-water as fast as pad deformation and permeability allow basal water to be discharged at x = 0. Case IV would then require the surface-inflection line of the ice stream to retreat into the ice sheet, as it is a thermal boundary (the melting isotherm), and this retreat collapses the ice sheet.

The amount of basal water along Ls depends on basal melting and freezing rates in cases I through IV, and in other cases that could be postulated. In cases I, II and III, the basic assumption is that the ice overburden is mostly supported by basal water pressure above water channels and by bedrock between water channels (even though the water film between channels has the overburden ice pressure). These cases require ice-stream applications of subglacial hydrology theories by Reference WeertmanWeertman (1986), Reference KambKamb (1987) and Reference LliboutryLliboutry (1987), among others. Case IV can be tested by theories of ice streams on deformable and permeable beds developed by Reference Alley, Blankenship and RooncyAlley and others (1987), Reference Linglc, Brown, Van der Veen and OcrlemansLinglc and Brown (1987), Alley (Reference Alley1989a, Reference Alleyb) and Reference MacAyeal, Van der Veen and OcrlemansMacAyeal (1989), for example.

Figure 7 shows an alternative to Figure 6 for presenting ø variations in terms of basal water distribution. Bedrock hills and hollows determine the distribution instead of bedrock ridges and troughs aligned in the direction of the ice stream. In case I, lakes exist only in the bottoms of isolated hollows, which act as slippery spots beneath the ice stream and a thin water film exists elsewhere, so that ø = 0 along Ls. In case II, basal water begins as lakes in isolated hollows and ends surrounding isolated hilltops along Ls, with ø = øG for ice above isolated lakes, ø = (hg/h)øG for ice surrounding isolated hills, and ø = 0 for ice in remaining areas where a thin water film allows bedrock to support the ice overburden. In case III, the thin water film exists only on isolated hilltops, which act as sticky spots where ø = 0, but for the ice stream in general, ø = (hg/h)øG along Ls. In case IV, sediments beneath the ice stream have been eroded away and redeposited as till at the grounding line, leaving basal water in a bedrock trough dammed by a till sill, so ø = øG along Ls for ice above the trough. Case IV produces fore-deepened troughs in the inter-island channels, straits and fiords of continental shelves that had been glaciated by marine ice sheets. In these cases, øG = 0 for ice above isolated bedrock lakes, but øG > 0 if the till sill is permeable.

Fig. 7 A basal water configuration in plan view and longitudinal cross-section for the ø variation of case I through case IV in Figure 5 that is different from the basal water configuration for these cases shown in Figure. 6. The bed consists of tills and hollows instead of ridges and valleus, with ø = 0 in white areas, ø = (h/hG in connected black areas, and ø = øG m disconnected black areas, for which øG = 0. In case IV, basal water is trapped in a glacially eroded trough dammed by glacially deposited sediment till(dotted) at the grounding line, situationsfound in the inter-island channels and fiords of deglaciated landscapes. The vertical is exaggerated.

A water film over hilltops distributed like islands in a subglacial lake allows the hills to be local “sticky spots” where bedrock bears the ice overburden, and therefore retards ice flow. Applying the Reference WeertmanWeertman (1986) analysis of basal water pressure to ice streams, the sticky spots in Figure 7 would be the upstream side of bedrock hills, where a thin water film allows bedrock to resist stream flow. Whether this situation exists in ice streams is unknown, but Reference Whillans and JohnsenWhillans and Johnsen (1983) observed something like it, possible frozen patches on a thawed bed, in sheet flow near Byrd Station in West Antarctica (see Fig. 1). Reference Vornbcrgcr and WhillansVornberger and Whillans (1986) have postulated sticky spots in Ice Stream B (see Fig. 1) as producing bands of surface crevasses.

Case I. Stream Flow with Minimal Basal Buoyancy

Stream flow with ø = 0 was used to produce ice-stream profiles in the CLIMAP reconstructions of global ice sheets at the last glacial maximum (Reference Hughes, Denton and HughesHughes, 1981a). Ice streams without basal buoyancy may exist as outlet glaciers through fiords in coastal mountains, where ice is often thick and has a high surface slope for ice streams, implying some traction from a bedrock floor that is rough and may be partly frozen.

Considerable attention has been focused on longitudinal deviator stress σ′xx and its downslope gradient in ice streams, notably by Reference WeertmanWeertman (1974), Reference Bindschadler and GoreBindschadlcr and Gore (1982), Reference AlleyAlley (1984), Reference Alley and WhillansAlley and Whillans (1984), Reference Mclnnes and BuddMclnnes and Budd (1984), Van der Veen (Reference Van der Veen, Van der Veen and Ocrlemans1985, Reference Van der Veen and Whillans1987), Reference MuszynskiMuszynski (1987), Reference Muszynski and BirchfieldMuszynski and Birchfield (1987) and Reference Whillans, Van der Veen and OcrlemansWhillans (1987). These studies, notably those of Reference Van der Veen, Van der Veen and OcrlemansVan der Veen (1987) and Reference MuszynskiMuszynski (1987), show that

is relatively unimportant in an ice stream. That could be ignored in ice streams fully buttressed by an ice shelf was assumed by Reference Hughes, Denton and HughesHughes (1981a) in the CLIMAP icc-shcct reconstructions. This assumption was based on the proposition that σ′xx at the ice-shelf grounding line is the buoyancy pulling stress minus a braking back-stress imposed by the ice shelf:
(19)

where hG and σG arc the ice thickness and the back-stress at the grounding line. For full ice-shelf buttressing, Reference Hughes, Denton and HughesHughes (1981a) assumed that:

(20)

Taking (d= 0 in Equation (3) maximizes ice-bed coupling beneath an ice stream and σG given by Equation (20) maximizes ice-shelf buttressing beyond an ice stream. Equation (20) holds exactly only when grounding lines completely enclose an ice shelf, making it a “pscudo ice shelf” above a subglacial lake (Reference Robin, Swithinbank and SmithRobin and others, 1970; Reference Oswald and RobinOswald and Robin, 1973). The ice stream has no pulling power in this case.

Equation (19) can be related to the pulling force, given by Equation (10), by requiring that ø = øG at the grounding line is given by:

(21)

A reliable expression for σG is a major problem in glaciology, because ice shelves have complex shapes and dynamics that impose both form drag and dynamic drag at the grounding line (Reference MacAyeal, Van der Veen and OcrlemansMacAyeal, 1987). However, in the idealized case of an ice shelf having width w and thickness h identical to wG and hG at the grounding line, a balance of forces along floating length LF from the grounding line to the calving front of the ice shelf gives where side shear stress Ts is assumed to be constant along LF (Reference ThomasThomas, 1977). Values of σG for the Ross Ice Shelf in Antarctica have been computed from field data by Reference Thomas and MacAyealThomas and MacAyeal (1982) and by Reference JezekJezek (1984).

The concave surface profile that is characteristic of stream flow can be produced even if ø = øG = 0. As seen in Equation (3), ϕ→0 exists when d→0, a physical condition in which the basal water is thick enough to drown bedrock projections that control the basal sliding velocity, but not the large projections, which therefore continue to support nearly all of the ice overburden. For smaller projections, the ice overburden is supported by the water layer. Hence, the fast sliding velocity of an ice stream can be attained without ubiquitous basal buoyancy, provided there is no icc-shclf buttressing G = 1) and the basal water film drowns rate-controling bedrock projections (Reference WeertmanWeertman, 1986).

In ice streams reconstructed for CLIMAP, the bed was assumed to be rigid, impermeable and rough in the manner specified in the basal sliding theory for sheet flow proposed by Reference WeertmanWeertman (1957b). Basal sliding for stream flow was attained by having the bedrock projections of height Λ, that control the sliding velocity for sheet flow, be progressively drowned in the downslope direction by a basal water layer of thickness λ that increased from λ ≈ 0 at x = Ls to λ = A at x = 0 along the ice stream. The variation of λ with x over length Ls was unknown, but on the assumption that ∂ λ/∂x = 0 at both x = Ls and x = 0, where x = 0 was the grounding line, the following cosine variation was employed:

(22)

The undrowned height of these projections was Λ –λ, and the bed-roughness factor for stream flow would then be:

(23)

where λ was the distance between projections of height λ and λ/λ was a constant bed-roughness factor for sheet flow in the Reference WeertmanWeertman (1957b) theory.

Basal sliding velocity u0 for sheet flow in the Reference WeertmanWeertman (1957b) theory was:

(24)

where

, and B0 was a constant that probably depended on basal water pressure (Reference Budd, Keagc and BlundyBudd and others, 1979). In CLIMAP ice-sheet reconstructions, Equations (C) and (24) were used to compute convex surface flowline profiles characteristic of sheet flow from a formula proposed by Reference NyeNye (1952). For a horizontal bed, the case considered here, his formula was:
(25)

The CLIMAP reconstructions assumed that basal shear stress To had two characteristic variations for sheet flow from the ice divide at x = L to the ice margin at x= 0, one being TF for a frozen bed computed from Equation (6) and the other being TM for a melted bed computed from Equation (24). If the bed consisted of mixed frozen and thawed patches, as may exist in transition zones bet ween zones where the bed is either completely frozen or completely thawed, then:

(26)

where F was the thawed fraction in a given Δx step of Equation (25).

Ablation in the CLIMAP ice sheets was confined to the first step in from the ice-sheet margin, with accumulation being constant over all subsequent steps to the ice divide. This constraint was replaced by Reference HughesHughes (1981b), and Table 1 lists the TM and Tf expression for the mass-balance expression:

(27)

where the equilibrium line is at x = E, accumulation rate a is constant over How-band length L-E, ablation rate b is constant over flow-band length L-E and net mass balance E is zero for a steady-state ice sheet, positive for an advancing ice sheet and negative for a retreating ice sheet.

Stream flow develops from sheet flow over a thawed bed, and E = hE = b = 0 in ice streams buttressed by ice shelves, which was the case for CLIMAP ice sheets. Applying these limitations to TM in Table 1 gives:

(28)

Since λ is raised to the 4 m/(2 m+l) power in Equation (28) for sheet flow in which bedrock projections A are not drowned, progressive drowning for stream flow according to Equation (22) requires that λ –λ replaces A as the undrowned height. Then, from Equation (23) and (28):

(29)

Substituting this expression for To into Equation (25) gives:

(30)

Equation (30) converts the convex sheet-flow profile over length L given by TM in Equation (28) into a concave stream-flow profile over length Ls. This is the concave profile for ø = 0, so that Px = 0 in Equation (18).

Case II. Stream Flow with Decreasing Basal Buoyancy

An expression for Δh/Δx can be derived by using the continuity condition to specify σ’w in Equation (15). For constant w, the continuity equation is:

(31)

where u is the average velocity of the ice column for mass-balance equilibrium. Integrating Equation (31) for constant a over length Ls and noting that u is negative:

(32)

where hG and UG are values of h and u at the grounding line. Equation (31) also yields:

(33)

where is the longitudinal strain-rate and is given by Equation (10), in which T is the effective shear stress and σ′xx is the longitudinal deviator stress. For constant width, the transverse deviator stress σ′yy is zero:

(34)

where is the depth-averaged hydrostatic pressure. Solving Equation (34) for σyy yields:

(35)

Important stress components in t are σxy, σxx, σyy, ayy and σzz. Using Equation (35):

(36)

Using Equation (35) to express σ′xx in terms of σxx and σ′zz

(37)

Equation (10) now becomes:

(38)

Equation (32) and (33) can be solved for AΔh/Δx, and Equation (38) can be substituted for in the resulting expression to give:

(39)

Equation (39) is difficult to evaluate further.

A major simplification of Equation (39) is possible if the braking effect of shear stresses xy and σxz can be represented by a longitudinal braking stress au that reduces pulling by the gravity-driven longitudinal deviator stress given by Equation (15). This implies a longitudinal braking force FR that arises from side and basal traction and increases with distance x upstream from the grounding line as follows:

(40)

where τ S, h and τ 0 are the average values of rs, h and r0over distance x. Equation (39) can then be replaced by:

(41)

There is no mathematical way to derive Equation (41) from Equation (39), but physically it means that the braking effect of shear stresses σxy and σxz in restraining stream flow is mimicked if σ’xx is reduced by σB.

Since ru results from basal traction caused by ice-bed coupling due to the weight of ice above buoyancy height being supported by the bed, the effect of r0 in Equation (41) is included if σ′xx is given by Equation (15). Equation (41) then becomes:

(42)

where ø expresses the degree of ice-bed coupling and ice-shelf buttressing, and is included in the τ S term because the braking force cannot exceed the pulling force.

The variation of Ts along x can be represented as a smooth decrease from a maximum at the head of the ice stream, where converging flow has produced strain-hardened ice having an upper viscophistic yield stress (TV)S at x = Ls, to a minimum at the foot of the ice stream, where laminar flow has produced strain-softened ice having a lower viscophistic yield stress (TV)Q at X = 0. An empirical expression that represents this transition is:

(43)

If strain hardening remains all along Ls, then (Tv)S =(Tv)G=Tv and Ts. If strain softening ends as shear rupture at the grounding line, (Tv)<J = 0 and If shear rupture occurs all along Ls, then (TV)S = (TV)Q = 0 and Ts = 0. This spectrum of Ts behavior obtained from Equation (43) is illustrated by the flow curves in Figure 8, in which increasing strain softening allows a given side shear strain ϵ xy to be produced by decreasing side shear stresses TS.

Fig. 8 Possible flow curves for lateral shear in an ice stream. Strain hardening in the zone of converging flow at the head of an ice stream causes an increase in side-shear stress Ts and transverse shear strain τ s until mscoplastic yie.ld stress TV is reached. Over length Ls of stream flow, Ts is constant if strain hardening and softening rates are in balance (cruve 1), Ts decreases if strain softening dominates (cruve 2 and 3) and rs drops to zero if strain softening leads to shear rupture (curue 3).

The average value of Ts is provided by the relationship:

(44)

Substituting Ts from Equation (43) for Ts in Equation (44) gives:

(45)

for which

at x = 0 and ϵ xy at x=s

Strain softening along x is not necessarily described by Equation (45). It can be related to the minimum value of TS if h is replaced by mean height and side-strain softening controls the side-shear force such that:

(46)

Figure 9 shows the effect of Equation (46) as x increases in successive Δx steps. The mean side area is (h1 + hG)Δx for the first step and (hS + hG)Ls for the last step. The side area along χ is overestimated more with each step if the ice-stream surface is highly concave. Progressive ovcrcstimation of side area offsets keeping Ts at its minimum in a way that links strain softening to the concavity of the ice stream. The sooner strain softening develops in its side-shear zones, the more concave the ice stream becomes. This is reasonable, because the ultimate strain softening is shear rupture, for which TS = 0 and Equation (42) has the most concave profile.

Fig. 9 A method for approximating the side area of an ice stream that uses successive straight lines to represent the concave top surface, with one line for each additional Δx inrecrement added to x. The less concave the surface, the better the approximation.

Case II is a situation in which basal permafrost thaws to produce a layer of water-soaked sediments through which bedrock crops out. An ice stream develops, with some of the ice overburden supported by these bedrock “sticky spots” and the remainder supported by basal pore-water pressure. Hence, To exists and its variation along χ can be represented by modifying Equation (28), in which T0 = TM for a thawed bed, to include the effect of basal buoyancy that arises when basal permafrost thaws. Since the basal permafrost consists of unconsolidated sediments with a high ice fraction, thawing allows the basal water pressure to support the entire ice overburden except at those sticky spots where bedrock projects through the sediments and penetrates the overlying ice. If λ is the height of these bedrock projections and λ is their average separation, and if λis the thickness of thawed water-soaked sediments between these projections, they can be related to basal buoyancy factor ø as follows:

(47)

Equation (47) requires that λ= 0 when ø= 0 for no buoyancy and λ=λ when ø = øG for full buoyancy, with a linear variation in between. Equation (28) can be modified to include this situation so that the following expression for T0 results:

(48)

where in the Reference WeertmanWeertman (1957b) theory of basal sliding. This assumes the large “sticky-spot” bedrock projections that control stream flow have a height λ –λ penetrating ice above thawed permafrost that is comparable to the “controling obstacle height” for sheet flow that Reference WeertmanWeertman (1957b) proposed.

An expression for Δh/Δx in which stream flow begins as sheet flow and ends as shelf flow is the following:

(49)

where Equation (46) is substituted into Equation (42), Equation (48) is substituted into Equation (25) and the two arc added. Since braking vanishes when pulling vanishes, (TV)G = 0 when ø = 0. In Equation (49), sheet flow exists when ø=0 and the second righthand term dominates, whereas shelf flow exists when ø = 0G and that term vanishes. Since braking stops when pulling stops, the limiting ice-stream length is:

(50)

If (Tv)G = 100 kPa, hG= 1000 m, h = 2000 m and w =30 km, Ls = 100 km.

Many ways to decrease ø as x increases along length Ls of stream flow arc possible, but one for which so that in Equation (3) both ø and dø/dx arc smooth functions is:

(51)

Basal water pressure decreases gradually from to x = 0 to σ0=0 at x = Ls- Equation (49) becomes:

(52)

Casc III. Stream Flow with Equilibrium Basal Buoyancy

If bedrock “sticky spots” cropping out in the water-soaked sediments beneath an ice stream are small and widely separated, or are not even present, and if the permeability of the sediments allows the pore-water pressure to be in hydrostatic equilibrium with basal water pressure at the grounding line, then all along Ls and Equation (3) becomes:

(53)

This is the case of constant basal water pressure, for which Equation (49) becomes:

(54)

Case IV. Stream Flow with Maximum Basal Buoyancy

If bedrock docs not crop out in the water-soaked sediments beneath an ice stream, and a till or bedrock sill at its grounding line is too impermeable to pore water in the sediments to establish hydrostatic equilibrium under the ice stream, then pore-water pressure supports all of the ice overburden. In this case, all along Ls, so no ice-bed coupling exists and ø in Equation (3) reduces to its value for ice-shelf buttressing:

(55)

This is the case of maximum basal water pressure, for which Equation (49) reduces to:

(56)

Pulling Power and a Life-Cycle Classification for Ice Streams

Pulling power provides an analytical foundation for classifying marine ice streams in a hierarchy in which loss of gravitational potential energy over time produces an ice-stream history characterized by inception, growth, mature, declining and terminal stages of a life cycle, but allowing for rejuvenation as well. Pulling power decreases as an ice stream ages, and quantifies the vitality of the ice stream at a given stage in its life cycle. The life cycle of an ice stream must be inferred from present-day ice streams and from the glacial geological record of former ice streams. These observations support an hypothesis in which the pulling power of an ice stream is closely related to the degree of ice-bed coupling beneath the ice stream and the degree of ice-shelf buttressing beyond the ice stream. Letting numbers 1, 2, 3, 4 and 5 denote no, weak, moderate, strong and full ice-bed coupling, and letters A, B, C, D and E denote no, weak, moderate, strong and full ice-shelf buttressing. Table 2 presents a life-cycle classification of ice streams in which various combinations of ice-bed coupling and ice-shelf buttressing arc quantified as values of ø. In this classification of an ice-stream life cycle, ø = 1 for the inception stage, ø = 3/4 for growth stages, ø = 1/2 for mature stages, ø = 1/4 for declining stages and ø = 0 for terminal stages.

Table. 2 A life-cycle classification for ice streams*

The fundamental assumption underlying the ice-stream classification in Table 2 is that marine ice streams begin in inter-island channels on the outer edge of broad polar continental shelves or in the straits where large continental embaymcnts open on to polar oceans. These continental shelves and embaymcnts arc sedimentary basins floored by permafrost at the beginning of a glaciation cycle, when sea ice thickens and grounds on the ice-cemented sediments, and continued thickening produces marine ice domes whose terrestrial margins transgress on to land and whose marine margins encroach on to the outer continental shelf (Reference Denton and HughesDenton and Hughes, 1981; Reference HughesHughes, 1982, Reference Hughes1987a; Reference Lindstrom and MacAyealLindstrom and MacAyeal, 1987, Reference Lindstrom and MacAyeal1989; Reference LindstromLindstrom, 1990). The life cycle of marine ice streams follows from this assumption.

Since the mass balance of an ice stream changes during its life cycle, pulling power Px should be computed for the actual mean ice velocity u x , not the velocity u for mass-balance equilibrium given by Equation (17). An equation for u x that is based on the basal buoyancy factor ø and combines mean sliding velocity u 0 with mean creep velocity u c along length Ls of stream flow is:

(57)

where and are given by Equation (11) and (7), respectively. Equation (57) gives u x for sheet flow when ø = 0 and u x for shelf flow when ø = 1 and Ls = 0, so that u S becomes u G , the mean ice velocity across the grounding line and x is negative for floating ice. When given by Equation (57) is multiplied by Fx given by Equation (10), the resulting expression for pulling power Px is a complex function of basal buoyancy factor ø.

From Equation (3), ø is a measure of ice-bed coupling. From Equation (21), ø is also a measure of ice-shelf buttressing. Therefore, a classification of ice streams for ice bed coupling increasing from 1 to 5 and for ice-shelf buttressing increasing from A to E can be based on specifying ø along an ice stream.

Inception Stage of Ice Streams

Marine ice accelerates as it enters inter-island channels and straits, where the bed is already thawed or becomes thawed by the increased basal frictional heat. Since the thawed bed consists of water-saturated sediments, 1A ice streams arc nucleated in these channels and straits because ice-bed coupling is nil and only sea ice exists beyond the ice streams. The Wilkes Land margin of the

East Antarctic ice sheet has the most likely candidates for 1A ice streams, notably Ninnis, Mertz, Tottcn and Vandcrford Glaciers in Figure 1 (Reference HughesHughes, 1987b).

Growth Stage of Ice Streams

Basal frictional heat is concentrated at the heads of channels and straits, where basal traction is a maximum at the basal melting-point isotherm that separates frozen from thawed parts of the bed. This is a highly unstable condition because frictional heat causes the isotherm to retreat rapidly. A maximum-slope surface-inflection line marks the head of these ice streams and exists above the basal isotherm. Hence the ice streams retreat rapidly, drawing down the marine ice dome and thereby raising sea level. Rising sea level floats the grounding line of marine ice streams over the submarine sills that typically exist at the oceanward ends of channels and straits, causing an ice shelf to form as grounding lines retreat. The sills are either bedrock ridges or grounding-line moraines, and they form as the thawed sediments are eroded by the ice streams, causing bedrock to crop out on the bed in places where sediments were thin. This situation produces IB and 2A ice streams, where grounding-line retreat provides increasing ice-shelf buttressing (IB ice streams) and bedrock cropping out provides increasing ice-bed coupling (2A ice streams). Pine Island Glacier and Thwaites Glacier arc respective IB, and 2A or 2B ice streams in the Pine Island Bay polynya of the Amundsen Sea embayment of the marine West Antarctic ice sheet (Reference HughesHughes, 1981b; Lindst rom and Hughes, 1984).

Mature Stage of Ice Streams

Pulling power of marine ice streams increases rapidly as grounding lines retreat downslope into the isostatically depressed marine sedimentary basins and decreases rapidly as grounding lines retreat upslopc out of these basins, since the pulling force is greatest when the grounding line is furthest below sea level. Hence, the mature stage of marine ice streams exists when a floating ice shelf occupies the seaward half of a marine sedimentary basin and a grounded ice dome occupies the landward half, with ice streams drawing down the marine ice dome and eroding thawed marine sediments beneath the dome. This situation produces 1C, 2C, 3C, 3B and 3A ice streams, depending on the degree of ice-bed coupling and ice-shelf buttressing for ice streams at various sites in the embayment. The Ross Sea and Wcddcll Sea embaymcnts of the marine West Antarctic ice sheet contain mature ice streams (Reference HughesHughes, 1977).

Declining Stage of Ice Streams

Retreating grounding lines of marine ice streams eventually stabilize when they lie in water too shallow to sustain a significant pulling force or lie against the head-walls of fiords through coastal mountains, while a floating ice shelf fills much of the marine sedimentary basin beyond the grounding line. Retreating inflection lines at the heads of ice streams can continue to retreat, pulling more ice into the marine basins and drawing down terrestrial ice domes. Marine ice streams become terrestrial ice streams at this stage, and they persist so long as accumulation over terrestrial ice domes replaces ice pulled out by stream flow. This loss of marine character is typical of ID, 2D, 3D, 4D, 4C, 4B and 4A ice streams, and therefore represents the declining stage of marine ice streams. Outlet glaciers through the Transantarc-tic Mountains are probably 3D, 4D and 4C ice streams that drain the largely terrestrial East Antarctic ice sheet and are buttressed by ice shelves floating in the marine sedimentary basins of West Antarctica (Reference SwithinbankSwithinbank, 1963).

Terminal Stage of Ice Streams

When an ice stream is no longer able to deliver ice to the calving front of its ice shelf fast enough to offset the iceberg-calving rate, the ice stream enters its terminal stage. Ice streams are terminal in stages IE, 2E, 3E, 4E, 5E, 5D, 5C, 5B and 5A. A calving bay carves away the ice shelf until its calving front coincides with its grounding line to produce a calving ice wall. A 5E ice stream ideally exists in a fiord floored by rugged bedrock, so that ice-bed coupling is optimized; and supplies an ice shelf grounded on all sides, so that icc-shclf buttressing is optimized. Lambert Glacier, which pulls terrestrial East Antarctic ice into the Amery Ice Shelf, may be a rejuvenated 5E ice stream (Reference AllisonAllison, 1979; Reference HughesHughes, 1987b). Rejuvenation of a 5E ice stream occurs when a calving bay reduces ice-shelf buttressing, converting it into an ice stream in stages moving from 5E to 5A, thereby causing a great increase of ice-stream velocity. This may be the case with Jakobshavns Isbra? that drains much of the west-central Greenland ice sheet into Jakobshavns Isfjord and is the fastest ice stream known (Reference HughesHughes, 1986).

Modeling Ice-Stream Life Cycles

About 90% of the ice in present-day ice sheets covering Antarctica and Greenland is discharged by marine ice streams. If this was true of the marine margins of former ice sheets, and if terrestrial ice streams were prominent along their terrestrial margins, then computer models that simulate the advance and retreat history of ice sheets should include any proposed life cycle for ice streams. The life cycle for ice streams proposed here would apply to deglaciation episodes in the history of a marine ice sheet, if the marine ice-transgression hypothesis explains its glaciation episodes (Reference HughesHughes, 1986). A marine ice sheet is postulated to form when sea ice thickens and grounds in shallow marine cmbayments or on shallow continental shelves fringed by islands. Examples are Hudson Bay, Foxe Basin and inter-island channels of the Queen Elizabeth Islands in North America, the Baltic, Barents and Kara Seas in Eurasia, Byrd Sub-glacial Basin and the Ross and Weddcll Seas in Antarctica. Grounding occurs mainly on sediments in the permafrost condition, so that the marine ice sheet grows on a frozen bed that precludes formation of marine ice streams. When a marine ice sheet becomes thick enough to thaw the basal permafrost, both by depressing the basal melting point and increasing basal-friction heating, ice streams will develop in outer inter-island channels where most marine ice is discharged. This begins the deglaciation history of the marine ice sheet and the life cycle of its ice streams.

Computer models of ice sheets that include the life cycle of marine ice streams need to incorporate three boundary conditions that constrain stream flow: (1) ice-shelf buttressing, (2) ice-bed coupling, and (3) side-shear traction. As defined by Equation (3), the basal buoyancy factor ø has potential for parameterizing the first two of these constraints, and the third constraint is made to depend on ø by way of Equation (46). The resulting expression for Δh/Δ is Equation (49) and it includes the effects of all boundary conditions for specified values of ø.

An empirical expression for ø that produces a ø spectrum close enough to the theoretical ø curves in Figure 5 to justify using it to simulate the life cycle of ice streams in computer models of ice-sheet dynamics is:

(58)

where c is a constant ranging from zero to infinity. Figure 10 is a plot of Equation (58) for comparison with Figure 5, showing that c = ∞ reproduces case I and c = 0 reproduces casc IV, the two end members of the ø spectrum. Case II is approximated by c = 1.0 and case III is approximated by c = 0.1. During the life cycle of an ice stream, c increases from zero to infinity and øG decreases from one to zero as time progresses, with c representing ice bed coupling and øG representing icc-shclf buttressing along the various pathways depicted in Table 2.

Fig. 10 The variation of ice-bed coupling ø/øa with distance x upstream /rom the grounding line of a marine ice stream of length Ls for c values ranging from zero to infinity in Equation (58).

Equation (49) can be integrated numerically to produce the changing profiles of ice streams during anice stream life cycle. If fiowline length L is divided into i steps of constant length Δx over which ice elevation increases by variable height where i is an integer and stream flow occurs over length Ls of L, Equation (49) can then be written:

(59)

Equation (59) can be numerically integrated directly using the Euler method or more accurately using the Runge-Kutta method. Ice-bed coupling is specified for 0 ≤ c ≤ ∞ and ice-shelf buttressing is specified for 0 ≤ øG ≤ 1.

Ice-stream surface profiles constructed on a horizontal bed using Equation (59) are plotted in Figures 11 through 14 for ρI = 960 kg m−3, g = 9.8 ms−2, L = 1000 km, Ls = 400 km, w = 30 km, h = 500 m, uG = 1000 m a−1, m = 2, n = 3, (τv)G = 10kPa, B = 2kPaa½ m½ (Reference Hughes, Denton and HughesHughes, 1981a) and A = for a mean temperature of −15°C through the ice thickness (Reference PatersonPaterson, 1981, Table 3.3).

Table. 3 Flowline. surface slopes above a horizontal bed

Figure 11 shows the effect of ice-bed coupling for the midpoint value of ice-shelf buttressing, øG = ½ Surface profiles arc plotted for c values of 0, 1/10, 1/5, 1/2, 1, 2, 5 and 10, following the c spectrum in Figure 10. The ice-stream surface elevation does not begin to climb along distance x upstream from the grounding line until c = 1/5, and the surface-inflection point separating lower concave parts from upper convex parts begins to migrate toward the grounding line when c > 1, being at 240 km for c = 2, 120 km for c = 5, 80 km for c = 10 and at the grounding line for c = ∞. Hence, the pulling power of the ice stream, as represented by a down-draw that creates a lowered concave surface reaches further and more strongly into the ice sheet as ice-bed coupling, represented by c, decreases.

Fig. 11 Evolution of an ice flowline profile from convex sheet flow to concave stream flow as ice-bed coupling exponent c in Equation(58) decreases from infinity to zero for intermediate ice-shelf buttressing given by øG = 0.5.

Figure 12 shows an instability in the c spectrum when øG = 1 for no icc-shclf buttressing. The first righthand term in Equation (59) is the pulling term, as it contains σ′ xx and it takes charge as c drops from 0.48 to 0.47. For higher values of c, the upper part of the ice stream has a convex profile that first migrates upstream and then migrates downstream as c increases, producing the wholly convex surface of sheet flow at c = ∞. For lower values of c, the lower part of the ice stream has a concave profile that moves toward the grounding line as c decreases, and produces a concave surface for c = 0 that intersects the convex surface for c = ∞ at 93 km from the grounding line and at an ice elevation of 930 m.

Fig. 12 The lengthening reach of a marine ice stream into an ice sheet during the inception and growth stages, shown as flow-line profiles for various c values in Equation (58) when øaG = 1 specifies no ice-shelf buttressing. As c increases, the concave stream-flow profile (excluding dashed parts of surface-elevation curves) migrates toward the ice divide, until substantial down-draw of sheet flow occurring after a surface-inflection instability at 0.47 < c < 0.48. Stream flow steadily reverts to sheet flow for c < 1, causing the flow-line elevation to increase and its profile to become increasingly convex.

Figure 13 shows the effect of ice-shelf buttressing for the midpoint value of ice-bed coupling c = 1. Surface profiles are plotted for ø G values of 0, 0.2, 0.4, 0.6, 0.8 and 1.0. All of these profiles arc concave and they are nearly superimposed for the middle range of 0.2 ≤ ø G ≤ 0.8. The concave profiles are higher for øa = 0, and especially for ø G = 1, because the surface elevation climbs close to the grounding line.

Fig. 13 Surface profiles of marine ice streams for ice-slielf buttressing increasing from, no buttressing (øG = 1) to fall buttressing (øG = 0) for moderate ice-bed coupling (c = 1 ). The lowest profiles arc for moderate ice-shelf buttressing (0.4 < øG < 0.6).

Figure 14 shows the effect of ice-shelf buttressing for less extensive ice-bed coupling represented by c = ½. Values of ø G are 0, 0.20, 0.50, 0.90, 0.95 and 1.00 to bring out the rapid upstream increase of ice elevation from ø G = 0.95 to ø G = 1.00, with a convex part developing as ice-shelf buttressing vanishes. The lowest icestrcam profile is for ø G = 0.50, the midpoint of ice-shelf buttressing.

Fig. 14 Surface profiles of marine ice streams for ice-shelf buttressing increasing from no buttressing (øG = 1) to full buttressing (øG = 0) for low ice-bed coupling (c = 0.5). The concave surface of stream, flow and the surface lowering of sheet flow develop most rapidhj at the onset of ice-shelf buttressing (1.00 < øG < 0.05).

Pulling Power Drives the Ice-Stream Life Cycle

The ice-stream pulling power obtained from Equation (18) is:

(60)

where h is obtained from Equation (59) and disregarding accumulation or ablation rates along the ice stream:

(61)

for conservation of volume flow under equilibrium conditions. Ideally, ux should be given by u x in Equation (57), which also depends on ø, but changes in height and length of an ice stream over time are slow enough to justify using Equation (61). Figures 15 through 18 plot Px versus a: using Equation (60) and various combinations of ice-bed coupling and ice-shelf buttressing. For convenience, Px is plotted as being positive, although ux would make it negative.

The inception stage in the life cycle of an ice stream is illustrated in Figure 15, for which c = 0.1 for minimal ice-bed coupling and ø G decreases from unity to zero as ice-shelf buttressing increases. Inception of stream flow begins without ice-shelf buttressing, and when basal mcltwatcr is produced from thawing permafrost at the head of the ice stream faster than it is discharged across the grounding line, so that the ice stream migrates into the ice sheet. Ideally, this is case IV with ø G 1. Retreat of the head and lengthening of the body of the ice stream arc too rapid for hydrostatic equilibrium to be established beneath the ice stream and for sheet flow to be down-drawn into the ice stream. Figure 15 shows concave stream-flow surface profiles as retreating and lengthening over the range 1.00 ≤ ø G ≤ 0.90 (top), during which pulling power Px increases with distance upstream from the grounding line (bottom). The sheet-flow surface profile, for which ø G = 0, is undisturbed beyond the ice stream. However, down-draw of sheet flow is dramatic over the range from 0.90 ≤ ø G ≤ 0.85 (top), during which Px changes dramatically from increasing to decreasing upstream from the grounding line (bottom). The conclusion to be drawn from this is that a little ice-shelf buttressing goes a long way in regulating both ice-stream retreat and ice-shect down-draw when h is formulated by Equation (59) and Px is formulated by Equation (60). In this formulation, the inception stage becomes the growth stage when ø G decreases from 0.90 to 0.85 and c is nil.

Fig. 15 Pulling power (hiring the life cycle of an ice stream emphasizing the inception stage, with c = 0.1 and øG decreasing from unity to zero. As the ice stream retreats without concomitant down-draw of the ice sheet (top), pulling power is concentrated at the head of the ice stream (bottom). After down-draw begins for øG < 0.90 (top), pulling power becomes greatest toward the foot of the ice stream (bottom).

An example of the growth stage in the ice-stream life cycle when c = 0.5 is illustrated by Figure 16. Pulling power is strongest at the grounding line and remains strong farther up the ice stream than for the other combinations of c and ø G in Figures 16 through 18. Compared to the inception stage, pulling power reaches further into the ice sheet but pulls ice out with weaker grip during the growth stage. This is because hydrostatic equilibrium becomes established beneath the ice stream during the growth stage.

The mature stage in the ice-stream life cycle is represented by the combination c = 1, øa G = 0.75 in Figure 17. Pulling power remains strong at the grounding line because an ice shelf is providing only moderate buttressing. Pulling power decreases rapidly as ø G decreases.

Fig. 17 Pulling power during the life cycle of an ice stream emphasizing the mature stage, with c = 1.0 and øG decreasing from unity to zero. Pulling power decreases steadily upstream from the grounding line, beginning at lower values as øG decreases,

Declining stages in the ice-stream life cycle arc shown when c = 0.5, 1, 2 are combined with 0.75 > ø G > 0.25 in Figures 16 through 18. The ice shelf has become somewhat confined in the embaymcnt and pinned to islands or shoals, so pulling power at the grounding line has been reduced. It reaches up the ice stream with decreasing strength as ice-bed coupling increases.

Fig. 16 Pulling power during the life cycle of an ice stream emphasizing the growth stage, with c = 0.5 and øG decreasing from unity to zero. Pulling power decreases slowly up-stream, from the grounding line, especially for the higher øG values.

Fig. 18 Pulling power during the life, cycle of an ice stream emphasizing the declining stage, with c = 2 and øG decreasing from unity to zero. Pulling power decreases rapidly upstream from the grounding line for larger øG values and is low even at the grounding line for small øG values. The terminal stages in the life cycle of an ice stream are represented by øG = 0 in Figures 15 through 18, for which Px = 0.

Terminal stages in the ice-stream life cycle are approached most closely in Figures 16 through 18 when c = 0.5, 1, 2 are combined with øG < 0.25. The ice shelf is now both confined and pinned in an cmbaymcnt studded with islands and shoals. Pulling power at the grounding line is so weak that the degree of ice- bed coupling upstream is relatively unimportant.

Pulling Power Disintegrates Marine Ice Sheets

The concept of a life cycle for marine ice streams in which pulling power is initially strong and remains strong upstream to the inflection line, as in 1A ice streams, and becomes weaker and more concentrated near the grounding line, to terminate finally as a stagnant 5E ice stream, provides a blueprint for disintegration of marine ice sheets. Several mechanisms linked to pulling power contribute to disintegration.

The first mechanism is down-draw, which measures how far longitudinal pulling power can reach up an ice stream and into the heart of a marine ice sheet. The second mechanism is discerpation of ice ridges between ice streams by the transverse pulling power of these ice streams. The third mechanism is basal erosion caused by the pulling power of ice streams, so that stream flow migrates toward the ice divide of marine ice sheets. The fourth mechanism is creation and destruction of ice rises by pulling power, thereby controling the ability of an ice shelf to buttress marine ice streams. The fifth mechanism is widening ice streams by lateral pulling power, thereby enlarging the ice-stream drainage basin. The sixth mechanism is ice-stream surges caused by intermittent increases in pulling power. These mechanisms will now be examined.

Pulling Power Down-Draws Marine Ice Sheets

Ice streams have surface profiles that reflect how far and how strong pulling power can reach up an ice stream and down-draw the ice sheet. The length of reach is controled by ice-bed coupling. The strength of reach is controled by icc-shclf buttressing.

The longest reach is provided by c = 0 and the greatest strength is provided by øG = 1 in Equation (59) and (60). Pulling power reaches full strength up the entire length L s of the ice stream. This occurs during the inception stage, and produces stage 1A ice streams. However, as shown in Figure 12, the concave part of stream flow is only about 93 km long and it is too early in the life cycle for interior ice to be down-drawn into the ice stream. Hence, the unlowered convex surface of sheet flow extends from there to the ice divide. There is no ice bed coupling beneath the ice stream and no ice-shelf buttressing beyond the ice stream.

Down-draw begins during the growth stage of ice streams, when some sticky spots appear beneath ice streams and a confined or pinned ice shelf begins to form seaward of retreating ice-stream grounding lines. These conditions exist for 0 < c. < 0.5 and 0.75 < øG < 1.00 in Equation (59) and (60). The reach and strength of pulling power are reduced, but the ice sheet is beginning to lower as interior ice is down-drawn into the ice streams. Surface profiles illustrating this situation are shown in Figure 14. Down-draw is greatest during the initial reduction of øG That is, when the ice shelf becomes confined and pinned enough to cause some initial buttressing, weak as it is, as seen by the drastic drop in ice elevation from øG = 1.00 to øG = 0.95.

Down-draw stabilizes during the mature stage of ice streams, when sticky spots and slippery spots beneath the ice stream are in approximately equal abundance and the ice shelf beyond the ice stream floats in a large confining cmbayment, where it is pinned to islands or shoals at several points. This situation is represented by c ≈ 1 and øG ≈ 0.75 in Equation (59) and (60). Surface profiles of mature ice streams are shown in Figure 13. The instability of surface profiles depicted in Figure 12 over the narrow range 0.47 < c < 0.48 when øG = 1 vanishes quickly as øG becomes smaller, and occurs at lower values of øG for lower c values, as seen in Figure 15. Whether this has any physical significance is unclear. Comparing Figures 13 and 14 shows that surface elevations at the head of the ice stream can be higher for the mature stage than for the growth stage in many cases. This reflects the steady reduction of pulling power from the growth stage to the mature stage. With less ice being pulled out through ice streams, interior ice elevations can stabilize and even recover.

Down-draw is absent during the declining stage of ice streams because their pulling power is too weak and too short in reach. The bed beneath the ice streams consists of slippery spots scattered in a sticky matrix and a confined and pinned ice shelf strongly buttresses the ice streams. Equation (59) and (60) produce these conditions for 1 < c < ∞ and 0.75 < øG < 0.25, with the convex surface of sheet flow extending closer to the grounding line as c decreases and øG decreases.

Sheet flow reaches the grounding line at the terminal stage of ice streams, for which c = ∞ and øG = 0 in Equation (59) and (60). With full ice-bed coupling and ice-shelf buttressing, pulling power vanishes because ø = 0 makes Fx = 0 in Equation (60). The ice sheet thickens because ux in Equation (60) is too small to discharge ice as rapidly as it precipitates over the ice sheet.

Pulling Power Discerps Ice Ridges

Ice ridges lie between ice streams. Since transverse profiles are often steeper that longitudinal profiles on ice ridges, flow down the flanks into ice streams is usually faster than flow along the crest into an ice shelf. The longitudinal profile of an ice ridge can be computed from solutions of Equation (39) at the head (x = L s) and the foot (x = 0) of its flanking ice streams, since these are also boundaries of the ice ridge. The longitudinal profile is the ice divide for transverse flow into the flanking ice streams.

Equation (39) has simple solutions at x = 0 and x = L s. At x = 0, σ xy is small compared to and σ xz τ0 = 0. Equation (39) at h = h G is then:

(62)

Reference WeertmanWeertman (1957a) derived an expression for the longitudinal strain rate in freely floating ice having constant width, and if the ice-shelf back-stress σG at the grounding line is included (Reference ThomasThomas, 1977), it becomes:

(63)

Equation (39), making use of Equation (21) and (63), is then:

(64)

At x = L s, let = 0 and σ xy << σ xz = τ0 = τM. Equation (39) is then replaced by Equation (25) with h = h s, because buoyancy vanishes at x = L s

(65)

where τ0 is given by Equation (28).

It is reasonable to assume that progressive ice-bed coupling from x = 0 to x = L s causes Equation (64) to be replaced gradually by Equation (65). One way to accomplish this is to write:

(66)

In Equation (66), τs is evaluated at the grounding line, where in Equation (14), so equating these equations requires that:

(67)

where τS is the side-shear stress between an ice stream and its flanking ice ridge. Equation (GG) is, of course, purely arbitrary. However, it has the virtue of illustrating that the pulling force Fx defined by Equation (16) controls Δhx at x = 0, and that the contribution to Δhx by Fx decreases to zero at x = L s in the same way that ø decreases in Equation (51).

A feature of the concave stream-flow profile given by Equation (66) is that it becomes a convex ridge-flow profile if:

(68)

so that Equation (25) becomes:

(69)

Since τ0 is constant along x in Equation (68), Equation (G9) can be integrated to give the parabolic-surface profile that is predicted by plasticity theory:

(70)

Equation (70) gives a parabolic ice-ridge profile when τ0 is given by Equation (8).

Comparing Equation (67) for τs alongside an ice stream with Equation (68) for τ0 beneath an ice ridge having no basal buoyancy shows that side traction exceeds basal traction for unit traction area of an ice ridge, because:

(71)

where w >> h s The ice ridge will become part of its flanking ice streams when total side traction equals total basal traction for an ice ridge, so that where H is the mean side height of the ice ridge above its bed and W is the width of the ice ridge. Hence, W is the minimum spacing between ice streams and, for equal basal and side traction and

:
(72)

Examples of how ice-stream pulling power disccrps intervening ice ridges are found along the West Antarctic grounding line of the Ross Ice Shelf, particularly the Siplc Coast, where Reference Whillans, Bolzan and ShabtaieWhillans and others (1987) have measured surface velocities and deformation, and Reference Shabtaie, Whillans and BentleyShabtaie and others (1987) have mapped surface and basal topography (sec Fig. 1). Ice ridge AB has a contorted surface similar to that of its flanking Ice Streams A and B, whereas ice ridge BC has a smooth surface similar to Ice Stream C, which flanks it on the north. All three ice streams have surfaces not far above buoyancy along much of their length, with Ice Stream B underlain by deforming sediments (Reference Blankenship, Bentley, Rooney and AlleyBlankenship and others, 1966), and much basal water along Ice Stream C (Reference Robin, Swithinbank and SmithRobin and others, 1970). Ice Streams A and B move rapidly, inland ice augments local accumulation to give ice ridge AB a moderate velocity, inland ice is diverted into Ice Streams B and C by a local ice dome on ice ridge BC, giving it a slow velocity from local accumulation only, and Ice Stream C scarcely moves at all (Reference Whillans, Bolzan and ShabtaieWhillans and others, 1987). From these observations, Ice Streams A and ? are type 1C or 2C, whereas Ice Stream C is type IE or 2E, being almost fully buttressed by the Ross Ice Shelf. Hence, Ice Streams A and B exert strong transverse pulling forces but Ice Stream C does not. Moreover, transverse ice velocity entering these ice streams is much higher from ice ridge AB than from ice ridge BC. Ice ridge AB is therefore discerped by pulling power, whereas ice ridge BC is not, so ice ridge AB has the contorted pulled-apart sur-face of an ice stream, whereas ice ridge BC has a smooth surface.

Pulling Power Measures the Erosive Power of Ice Streams

Subglacial erosion requires a pulling force to loosen basal material and an ice velocity to transport the loosened material. Pulling power, being the product of this force and velocity, is a direct measure of the erosive capacity of an ice stream. Therefore, classifying ice streams according to pulling power provides a framework for answering the question posed by Reference BaderBader (1961) as to whether an ice stream is likely to create its own subglacial depression. If stream flow develops independently of bed topography and if pulling power can migrate upstream, the ice stream will erode a linear basal depression. Flow converging downslope from saddles on the ice divide should be able to produce an ice stream without a trough in the bed. However, pulling power at the head of tin ice stream co-exists with a peak in basal shear stress, so the bed should be strongly eroded beneath the surface-inflection line. Retreat of the inflection line will be at a velocity that equals the erosion rate, and will produce a trough beneath the ice stream. The trough should form most rapidly for a 1A ice stream, because basal erosion is merely a consequence of thawing ice-cemented per-mafrost, and most slowly for a 4D ice stream, because basal erosion requires quarrying bedrock. In short, the greater its pulling power, the more easily an ice stream can erode a trough.

Pulling power also allows an assessment of the hypothesis by Reference BaderBader (1961) that an ice stream is self-perpetuating once it forms, because ice viscosity is lowered by the heat of internal friction. Frictional heat per unit area is generated by basal traction, and is greatest for a 5A ice stream, which has minimal basal buoyancy and icc-shclf buttressing. Frictional heat per unit volume is generated by side traction, which is greatest for a 1A ice stream, which has maximum basal buoyancy and minimum icc-shclf buttressing. Basal and side traction combine not only to reduce effective viscosity along these boundaries by generating easy-glide ice fabrics, but also to produce fiord-like channels by basal and side erosion. Both processes tend to stabilize the ice stream, making it self-perpetuating, and thereby allowing long-term erosion.

An early stage in channel formation might be taking place at the head of Ice Stream B on the Siple Coast of West Antarctica. Reference Whillans, Bolzan and ShabtaieWhillans and others (1987) reported that the head of Ice Stream B is a region where “rafts” of inland ice arc being pulled into the ice stream. This is a region of converging flow, where the bed may be a mosaic of frozen and thawed patches that constitutes a transition from sheet flow over frozen permafrost to stream flow over thawed permafrost. If the “rafts” of relatively undeformcd ice exist over frozen patches of permafrost, disintegration of the converging flow zone may proceed because pulling power is a maximum and can tear out these ice rafts. Peak pulling power plucks!

Pulling Power may Regulate Ice-Shelf Buttressing

Ice shelves are able to buttress ice streams because they exert form drag and dynamic drag (Reference MacAyeal, Van der Veen and OcrlemansMacAyeal, 1987), largely because of ice rises, where ice shelves are locally grounded. Crary Ice Rise is a complex of local ice domes on the Ross Ice Shelf that has been variously described as becoming grounded or ungrounded from a ridge on the Ross Sea floor (Reference MacAyeal, Bindschadlcr, Shabtaie, Stephenson and BentleyMacAyeal and others, 1987). This bedrock ridge seems to continue ice ridge AB separating West Antarctic Ice Streams A and B (Reference BentleyBentley, 1984). It may therefore represent an advanced stage of disintegration of ice ridge AB. However, it also lies immediately downstream from Ice Stream B, which is car rying several “rafts” of relatively thick and undeformed ice on to the Ross Ice Shelf. Reference Whillans, Bolzan and ShabtaieWhillans and others (1987) traced these “rafts” to high-stress regions at the head of Ice Stream B, notably the “Unicorn” where Ice Stream B forks into branches B and B2. This suggests that the Crary Ice Rise complex resulted from the pile-up of ice rafts against the bedrock ridge. In that case, Crary Ice Rise was created by pulling power disintegrating the head of Ice Stream B. Moreover, pulling power may be disintegrating Crary Ice Rise by plucking rafts of ice from its sides and lee end. Hence, pulling power may both create and destroy certain kinds of ice rises. These ice rises are temporary docking sites for ice rafts discharged on to ice shelves by ice streams. Pulling power regulates ice-shelf buttressing for these kinds of ice rises.

Pulling Power Widens Ice Streams

Ice Stream A supplies the southernmost corner of the Ross Ice Shelf, and it forms from the confluence of East Antarctic ice from Reedy Glacier and West Antarctic ice from Horlick Ice Stream. Some East Antarctic ice also enters Horlick Ice Stream from Shimizu Ice Stream and local ice from the Wisconsin Range supplies the head of Ice Stream A, in addition to its branches into East and West Antarctica.

Radio-echo sounding by Reference Shabtaie, Whillans and BentleyShabtaie and others (1987) showed that the lateral shear zones defining the sides of Ice Stream A lie about midway between two channels, a deep narrow channel along the center of a shal-low wide channel. The deep narrow channel is presumably a fiord that turns southward beneath Reedy Glacier into East Antarctica, where it ends at a headwall in the Transantarctic Mountains. The wide shallow channel may have been eroded in ice-cemented sediments that form a mantle of permafrost over the bedrock of West Antarctica. The deep narrow channel would have been cut by East Antarctic ice from Reedy Glacier at a time when the West Antarctic ice sheet was much smaller than today or was absent. The shallow wide channel would have formed when the West Antarctic ice sheet was much larger than today, so that increased discharge from Horlick Ice Stream thawed the permafrost and eroded the sediments.

The fact that Ice Stream A has a width today that lies between the deep and shallow channels suggests that it is narrowing as West Antarctic ice lowers from its elevation at the last glacial maximum. Narrowing would then be a consequence of decreased pulling power as Ice Stream A ages from a mature 3C stage to a terminal 5E stage. Renewed growth of West Antarctic ice or, conversely, disintegration of the Ross Ice Shelf would increase pulling power and widen Ice Stream A. As Equation (72) shows, ice streams widen as they steepen, since W/w decreases when hG/hs decreases.

Pulling Power Causes Ice Streams to Surge

As Reference BaderBader (1961) observed, “An ice stream is something akin to a mountain glacier … but a mountain glacier is laterally hemmed in by rock slopes, while the ice stream is contained by slower moving surrounding ice.” Some mountain glaciers surge periodically, so periodic surges of ice streams should also be possible. Indeed, the Dibble and Dalton Iceberg Tongues grounded for 100 km on the continental shelf beyond unimpressive present-day ice streams in Wilkes Land, East Antarctica, may have formed when these ice streams surged, and Thwaites Glacier in Pine Island Bay, West Antarctica, may now be surging (see Fig. 1). Pulling power first increases and then decreases during a surge cycle, as both surface slope and ice velocity increase and decrease at the head of the mountain glacier or ice stream, where the surge begins, probably in phase with increases and decreases of basal buoyancy (Reference KambKamb and others, 1985).

If basal buoyancy, and therefore pulling power, controls surges of marine ice streams, surges can be triggered by rising sea level when the grounding line retreats faster than the inflection line. Rising sea level that moves the grounding line closer to the inflection line, thereby increasing surface slope, also lifts the ice shelf off its pinning points, thereby decreasing buttressing. Both processes increase the pulling power of the ice stream. As more ice is pulled out, the inflection line must retreat and pulling power extends its reach into the ice sheet. A 5E ice stream is most sensitive to changing sea level because basal buoyancy and ice-shelf buttressing change most dramatically at the grounding line. In contrast, a 1A ice stream is least sensitive to changing sea level because no ice shelf buttresses its grounding line and full buoyancy extends back to its inflection line. There may be a natural evolution from a 1A to a 5E ice stream, provided that a 1A ice stream can erode its thawed basal sediments down to bedrock and its actively retreating inflection line can drag along its grounding line, to leave behind a buttressing ice shelf.

Table 3 compares the surface-slope expressions for sheet flow and ridge flow with those for stream flow having ice bed coupling conditions specified by basal buoyancy factor ø. Inception, growth, mature, declining and terminal stages of an ice-stream surge may mirror these stages for the life cycle of an ice stream by having a similar dependence on ø through time that includes the ø variations in Table 3. The Δh/Δx expression for stream flow for minimum basal buoyancy is different from Equation (30), the Δh/Δx expression for case I, because their derivations were different. However, the Δh/Δx expressions in cascs II, III and IV in Equations (52), (54) and (56) for decreasing equilibrium and maximum basal buoyancy were based on the Δh/Δx expression in Table 3. Note that (τv)G = 0 when ø = 0.

Discussion

Badcr (1961) wondered whether “the ice stream, once started, is self-perpetuating because its mass, warmed up by heat of internal friction, has a lower viscosity than the surrounding ice.” This lower viscosity need not exist in the bulk of the ice stream, only in its side and basal boundary zones, where frictional heating can substantially decouple the ice stream from flanking ice ridges and the bed, and where laminar flow might produce an easy-glide ice fabric. A question not asked by Reference BaderBader (1961) is “How do ice streams start in the first place?” The answer to this question is part of “Glaciology’s grand unsolved problem” posed by Reference WeertmanWeertman (1976).

An attempt at answering this question was to suggest that ice streams were a natural consequence of ther-mal convection in polar ice sheets, as horizontal advection of ice toward the ice-sheet margin increasingly modified vertical convection of ice (Reference HughesHughes, 1976). Being warmed by internal frictional heat, the ice stream would by buoyed upward, thereby decoupling it from the bed and allowing stream flow. A more plausible alternative is that stream flow begins when the cold-ice ceiling collapses into the warm-ice basement of a polar ice sheet, creating the concave longitudinal and transverse profiles that characterize stream flow. This is somewhat analogous to the cold lithosphre descending into the warm aesthenosphcrc of the Earth’s mantle to produce crustal trenches. Polar ice sheets are like a miniature mantle in this respect, the major difference being that advection of ice from interior domes toward the margin compels ice-sheet “trenches” to form near the margin and move toward the margin as ice streams. Being partly collapsed “ceiling” ice, ice streams arc colder than flanking uncollapscd ice, but they ture also thinner than flanking ice so that basal mcltwater is driven toward ice streams by the hydrostatic pressure gradient. It is this continuous supply of mcltwater beneath ice streams that uncouples them from the bed enough to allow stream flow to start. This process would actually prevent thermal convection in ice streams, because the return cycle of convectivo flow, which ordinarily would allow thermal convection cells or rolls to develop, is short-circuited by advective flow, so that ice is re-routed through ice streams to the ice margin, never to return. In a sense, this is aborted thermal convection of the kind that can be initiated, but not sustained (Reference HughesHughes, 1985).

Comparing vertical temperature profiles just above and below the sheet-flow to stream-flow transition, and in the flanking ice ridges, might determine if collapse of ceiling ice takes place. The ratio of cold ceiling ice to warm basement ice should be greatest in the ice stream. Ice Stream B, located in Figure 1, would be suitable for this experiment, since it is already being studied in detail.

Although this analysis of pulling power was focused on marine ice streams, it can be generalized to include terrestrial ice streams. This can be done through the basal buoyancy factor ø. Ice-bed coupling for marine and terrestrial ice streams should have the same range, perhaps with total uncoupling being somewhat less likely for terrestrial ice streams, because drainage of basal mcltwater around the perimenter of their terminal ice lobes is possible. Ice-lobe buttressing is similar to ice-shelf buttressing, where back-stress øG is at the grounding line for an ice shelf and at the minimum-slope inflection line for an ice lobe.

The simplest braking force provided by an ice shelf of constant width w, grounded along both sides of float-ing length L F and pinned by an island or shoal having transverse diameter D, is:

(73)

where h G is grounding-line ice thickness, h is the average ice-shelf thickness, τs is side-shear stress alongside L F and σ c is longitudinal compressive stress pushing against D. The simplest braking force provided by an ice lobe of constant width w and grounded length L G is:

(74)

where hG is the grounded ice thickness where the lobe begins and τ0 is the basal shear stress beneath the lobe.

Solving Equation (73) and (74) for σG gives expressions for σG in Equation (21) that allows øG to be computed for an ice shelf and an ice lobe having these simple configurations. These values of øG can then be used in Equation (3) to compute ø for buttressing of a marine ice stream by an ice shelf or a terrestrial ice stream by an ice lobe, for these simple configurations. In general, there-fore, ice-bed coupling, ice-shelf buttressing and ice-lobe buttressing can be used to compute the basal buoyancy factor ø for any ice stream, marine or terrestrial.

Acknowledgements

For mc, Charles Swithinbank is the Father of Ice Streams, for his pioneering work on Antarctic outlet glaciers supplying the Ross Ice Shelf (Reference SwithinbankSwithinbank, 1963). On one of them, Byrd Glacier, he revealed Swith-inbank’s law to me: “All knowledge is obvious — once you have found it out.” This work is dedicated to him and the pursuit of his law. In my pursuit reported here, I became aware of another law: if something is obvious to me, but I cannot make it obvious to my readers, only one thing is obvious, it isn’t obvious to me either.

I presented the germ of the pulling-power concept at the 1986 AGU Chapman Conference on Fast Glacier Flow. Since then I have become indebted to glaciologists on four continents who helped mc to clarify and quantify the concept. In particular, I thank J. Fastook, who produced and plotted the numerical computations of surface profiles and pulling power for ice streams presented here, and R. Franzosa, for a mathematical analysis of the assumptions upon which these plots are based. These gentlemen will publish their own original work separately. Among editors and referees, G. Clarke, K. Echehneycr, D. MacAyeal and C. Raymond have been particularly helpful. Most importantly of all, I thank B. Hughes for typing and processing seven versions of this manuscript for two journals. This work was funded by the U.S. National Science Fondation (grant DPP-8400886) and the U.S. Department of Energy (Battelle, PNL, subcontract 07112-A-B1).

The accuracy of the references in the text, and in this list is the responsibility of the author, to Moro queries should be addressed.

References

Alley, R. B. 1984. A non-steady ice sheet model incorporating longitudinal stresses. Ohio Stute Univ. Inst. Polar Stud. Rep. 84.Google Scholar
Alley, R. B. 1989a. Water-pressure coupling of sliding and bed deformation: I. Water system. J. Glaciol., 35 (119), 108118.Google Scholar
Alley, R. B. 1989b. Water-pressure coupling of sliding and bed deformation: II. Velocity-depth profiles. J. Glaciol., 35 (119), 119-129.Google Scholar
Alley, R. B. and Whillans, I. M. 1984. Response of the East Antarctic ice sheet to sea-level rise. J. Geophys. Res., 89 (C4), 64876493.Google Scholar
Alley, R. B., Blankenship, D. D., Rooncy, S. T. and CR. Bentley. 1987. Till beneath Ice Stream B. 4. A coupled ice-till flow model. J. Geophys. Res., 92 (B9), 8931 8940.Google Scholar
Allison, I. 1979. The mass budget of the Lambert Glacier drainage basin, Antarctica. J. Glaciol., 22 (87), 223235 CrossRefGoogle Scholar
Bader, H. 1961. The Greenland ice sheet. Hanover, NH, U.S. Army Cold Regions Research and Engineering Laboratory.Google Scholar
Bentley, C. R. 1984. The Ross Ice Shelf Geophysical and Glaciological Survey (RIGGS): introduction and summary of measurements performed. Antarct. Res. Ser., 42, 120.Google Scholar
Bindschadler, R. A. and Gore, R. 1982. A time-dependent ice sheet model: preliminary results. J. Geophys. Res., 87 (C12), 96759685.Google Scholar
Blankenship, D. D., Bentley, C. R., Rooney, S. T. and Alley, R. B. 1986. Seismic measurements reveal a saturated porous layer beneath an active Antarctic ice stream. Nature, 322 (0074), 5457.Google Scholar
British Glaciological Society. 1949. Joint meeting of the British Glaciological Society, the British Rhcologists’ Club and the Institute of Metals. J. Glaciol., 1 (5), 231240.Google Scholar
Budd, W. F., Keagc, P. L. and Blundy, N. A. 1979. Empirical studies of ice sliding. J. Glaciol., 23 (89), 157170.CrossRefGoogle Scholar
Denton, G. H. and Hughes, T. J., eds. 1981. The last great ice sheets. New York, etc., John Wiley and Sons.Google Scholar
Doake, C. S. M., Frolich, R. M., Mantripp, D. R., Smith, A. M. and Vaughan, D. G. 1987. Glaciological studies on Rutford Ice Stream, Antarctica. J. Geophys. Res., 92 (B9), 8951 8960.Google Scholar
Hughes, T. J. 1973. Is the West Antarctic ice sheet disintegrating? J. Geophys. Res., 78 (33), 78847910.CrossRefGoogle Scholar
Hughes, T. J. 1970. The theory of thermal convection in polar ice sheets. J. Glaciol., 16 (74), 4171.CrossRefGoogle Scholar
Hughes, T. J. 1977. West Antarctic ice streams. Rev. Geophys. Space Phys., 15 (1), 146.Google Scholar
Hughes, T. J. 1981a. Numerical reconstruction of paleo-ice sheets. In Denton, G. H. and Hughes, T. J., eds. The last great ice sheets. New York, etc., John Wiley and Sons, 221261.Google Scholar
Hughes, T. J. 1981b. Correspondence. The weak underbelly of the West Antarctic ice sheet. J. Glaciol., 27 (97), 518525.Google Scholar
Hughes, T. J. 1982. Did the West Antarctic ice sheet create the East Antarctic ice sheet? Ann. Glaciol., 3, 138145.Google Scholar
Hughes, T. J. 1985. Thermal convection in ice sheets: we look but do not sec. J. Glaciol., 31 (107), 3948.Google Scholar
Hughes, T. J. 1986. The Jakobshavns effect. Geophys. Res. Lett., 13 (1), 4648.Google Scholar
Hughes, T. J. 1987a. Deluge II and the continent of doom: rising sea level and collapsing Antarctic ice. Boreas, 16 (2), 89100 CrossRefGoogle Scholar
Hughes, T. J. 1987b. The marine ice transgression hypothesis. Geogr. Ann., 69A (2), 237250.CrossRefGoogle Scholar
Jezek, K. C. 1984. A modified theory of bottom crevasses used as a means for measuring the buttressing effect of ice shelves on inland ice sheets. J. Geophys. Res., 89 (B3), 19251931.Google Scholar
Kamb, W. B. 1987. Glacier surge mechanism based on linked cavity configuration of the basal water conduit system. J. Geophys. Res., 92 (B9), 90839100.Google Scholar
Kamb, W. B. and 7 others. 1985. Glacier surge mechanism: 1982-1983 surge of Variegated Glacier, Alaska. Science, 227 (4686), 469479 Google ScholarPubMed
Lindstrom, D. R. 1990. The Eurasian ice sheet: formation and collapse resulting from natural CO9 concentration variations. Paleoceanography, 5 (2), 207228.Google Scholar
Lindstrom, D. R. and Hughes, T. J. 1984. Downdraw of the Pine Island Bay drainage basins of the West Antarctic ice sheet. Antorct. J. U.S., 19 (5), 5658.Google Scholar
Lindstrom, D. R. and MacAyeal, D. R. 1987. Environ-mental constraints on West Antarctic ice-sheet for-mation. J. Glaciol., 33 (115), 346356.CrossRefGoogle Scholar
Lindstrom, D. R. and MacAyeal, D. R. 1989. Scandinavian, Siberian, and Arctic Ocean glaciation: effect of Holoccnc atmospheric CO2 variations. Science, 245 (4918), G28631.CrossRefGoogle ScholarPubMed
Lindstrom, D. R. and Tyler, D. 1984. Preliminary results of Pine Island and Thwaites glaciers study. Antarct. J. U.S., 19 (5), 5355.Google Scholar
Linglc, C. S. and Brown, T. J. 1987. A subglacial aquifer bed model and water pressure dependent basal sliding relationship for a West Antarctic ice stream. In Van der Veen, C. J. and Ocrlemans, J. eds Dynamics of the West Antarctic ice sheet. Dordrecht, etc., D. Reidel Publishing Company, 249285.CrossRefGoogle Scholar
Lliboutry, L. 1958. Glacier mechanics in the perfect plasticity theory. J. Giaciol., 3 (23), 1G2169.Google Scholar
Lliboutry, L. 1987. Realistic, yet simple bottom boundary conditions for glaciers and ice sheets. J. Geophys, Res., 92 (B9), 91019109.Google Scholar
MacAyeal, D. R. 1987. Icc-shclf backpressure: form drag versus dynamic drag. In Van der Veen, C. J. and Ocrlemans, J., eds. Dynamics of the West Antarctic ice sheet. Dordrecht, etc., D. Reidel Publishing Company, 141160 Google Scholar
MacAyeal, D. R. 1989. Large-scale ice flow over a viscous basal sediment: theory and application to Ice Stream B, Antarctica. J. Geophys. Res., 94 (B4), 40714087.Google Scholar
MacAyeal, D. R., Bindschadlcr, R. A., Shabtaie, S., Stephenson, S. and Bentley, C. R. 1987. Force, mass, and energy budgets of the Crary Ice Rise complex, Antarctica, J. Glaciol., 33 (114), 218230.Google Scholar
Mclnnes, B. J. and Budd, W. F. 1984. A cross-sectional model for West Antarctica. Ann. Glaciol., 5, 9599 Google Scholar
Muszynski, I. 1987. The dynamics of coupled marine ice stream-ice shelf systems. (Ph.D. dissertation, Northwestern University, Evanston.)Google Scholar
Muszynski, I. and Birchfield, G. E. 1987. A coupled marine icc-strcam-icc shelf model. J. Glaciol., 33 (113), 315 Google Scholar
Nye, J. F. 1952. A method of calculating the thicknesses of the ice-sheets. Nature, 169 (4300), 529530.Google Scholar
Nye, J. F. 1957. The distribution of stress and velocity in glaciers and ice-sheets. Proe. R. Soc. London, Ser. A, 239 (1216), 113133.Google Scholar
Nye, J. F. 1958. Comments on Professor Lliboutry’s paper. J. Glaciol., 3 (23), 170172.Google Scholar
Oswald, G.K. A. and Robin, G.deQ. 1973. Lakes beneath the Antarctic ice sheet. Nature, 245 (5423), 251254 Google Scholar
Paterson, W. S. B. 1981. The physics of glaciers. Second edition. Oxford, etc., Pergamon Press.Google Scholar
Robin, G.deQ. 19G7. Surface topography of ice sheets. Nature, 215 (5105), 10291032.Google Scholar
Robin, G. deQ., Swithinbank, C. W. M. and Smith, B. M. E. 1970. Radio echo exploration of the Antarctic ice sheet. International Association of Scientific Hydrology Publication 8G (ISAGE), 97115.Google Scholar
Shabtaie, S. and Bentley, C. R. 1988. Ice-thickness map of the West Antarctic ice streams by radar sounding. Ann. Glaciol., 11, 126136.Google Scholar
Shabtaie, S., Whillans, I. M. and Bentley, C. R. 1987. The morphology of ice streams A, B, and C, West Antarctica, and their environs.. J. Geophys. Res., 92 (B9), 88658883.Google Scholar
Stuiver, M., Denton, G. H., Hughes, T. J. and Fas-took, J. L. 1981. History of the marine ice sheet in West Antarctica during the last glaciation: a working hypothesis. In Denton, G.H. and Hughes, T. J., eds. The last great ice sheet. New York, John Wiley and Sons, 319436.Google Scholar
Swithinbank, C. W. M. 1963. Ice movement of valley glaciers flowing into the Ross Ice Shelf, Antarctica. Science, 141 (3580), 523524.Google Scholar
Thomas, R. H. 1977. Calving bay dynamics and ice sheet retreat up the St Lawrence valley systems. Géogr. Phys. Quat., 31 (34), 347356.Google Scholar
Thomas, R, H. and MacAyeal, D. R. 1982. Derived characteristics of the Ross Ice Shelf, Antarctica.. J. Glaciol, 28 (100), 397412 Google Scholar
Van der Veen, C. J. 1987. Longitudinal stresses and basal sliding: a comparative study. In Van der Veen, C. J. and Ocrlemans, J., eds. Dynamics of the West Antarctic ice sheet. Dordrecht, etc., D. Reidel Publishing Company, 223248.Google Scholar
Van der Veen, C. J. and Whillans, I. M. 1989. Force budget: I. Theory and numerical methods. J. Glaciol., 35 (119), 5360.Google Scholar
Vornbcrgcr, P. L. and Whillans, I. M. 1986. Surface features on Ice Stream B, Marie Byrd Land, West Antarc-tica. Ami, Glaciol., 8, 168170 Google Scholar
Weertman, J. 1957a. Deformation of floating ice shelves. J. Glaciol., 3 (21), 3842 Google Scholar
Weertman, J. 1957b. On the sliding of glaciers. J. Glaciol, 3 (21), 3338 Google Scholar
Weertman, J. 1963. Profile and heat balance at the bottom surface of an ice sheet fringed by mountain ranges. International Association of Scientific. Hydrology Publication 61 (General Assembly of Berkeley 1963 Snow and Ice), 245252 Google Scholar
Weertman, J. 1974. Stability of the junction of an ice sheet and an ice shelf. J. Glaciol., 13 (67), 311 Google Scholar
Weertman, J. 1976. Glaciology’s grand unsolved problem. Nature, 260 (5549), 284286 CrossRefGoogle Scholar
Weertman, J. 1986. Basal water and high-pressure basal ice. J. Glaciol., 32 (112), 4554G3.CrossRefGoogle Scholar
Whillans, I. M. 1987. Force budget of ice sheets. In Van der Veen, C. J. and Ocrlemans, J., eds. Dynamics of the West Antarctic ice sheet. Dordrecht, etc., D. Reidel Publishing Company, 173G.CrossRefGoogle Scholar
Whillans, I. M. and Johnsen, S. J. 1983. Longitudinal variations in glacial flow: theory and test using data from the Byrd Station strain network, Antarctica. J. Glaciol., 29 (101), 7897.Google Scholar
Whillans, I. M., Bolzan, J. and Shabtaie, S. 1987. Velocity of ice streams B and C, Antarctica. J. Geophys. Res., 32 (B9), 88958902.Google Scholar
Figure 0

Fig. 1 Location map for selected Antarctic ice streams.

Figure 1

Fig. 2 Producingstream /loti) from sheet flow, (a) In sheet flow, gravitational pulling force FZ deforms the ice sheet from the solid profile to the clashed profile, (b) Inner dashed part of the ice sheet can be replaced by horizontal pushing force Fx. (c) Outer dashed part of the ice sheet can be replaced by horizontal pulling force Fx. (dj Horizontal pulling converts the connex surface of sheet flouj into the concave surface of stream flow, with stream flow becoming a floating ice shelf in terminus 1 or a grounded ice lobe in terminus 2.

Figure 2

Fig. 3 Pulling forcea defined us net horizontal lithostatic and hydrostatic forces acting in a marine flow band of an ice sheet. Horizontal lithostatic and hydrostatic forces pressing against imaginary ice columns are shown for sheet flow (right), stream flow (center) and shelf flow (left), with the. flow band having floating length Lp, grounded length L and streaming length Ls- Averaged for the ice column, basal water rises to depth d in imag-iiiary temperate boreholes drilled- through the flow band. Four cases analyzed in the text for the ice stream are no buoyancy along Ls (case I), full buoyancy at the basal grounding line decreasing to no buoyancy at the surface-inflection line (case II), constant buoyancy along Ls (case III) and full buoyancy along Ls (case IV).

Figure 3

Fig. 4 Horizontal forces on ice columns having Constant width and no side traction. Horizontal gravitational forces per unit width are the areas of triangles 1, 2 and 4, and rectangle 3. Traction force is basal shear stress To times the basal area of the ice column. Pulling force is tensile deviator stress 2′σxx times the transverse cross-sectional area on the landward side of the ice column. Following Whillans (1987), the role of basal uiater pressure can be emphasized by distinguishing lithostatic pressure in triangles 1 and 2, and rectangle 3 from hydrostatic pressure in triangle 4.

Figure 4

Table. 1 Variations of basal shear stress Towith distance L along flow bands of length L and constant width for sheet flow over frozen (T0 =TM) and thawed (T0 = TM) beds for constant surface accumulation and ablation rates a and b separated by an equilibrium line at distance x = E from the margin of the ice sheet (Hughes, 1981b)*

Figure 5

Fig. 5 A possible spectrum of ø variations along x for ice streams of length Ls. The four cases of basal buoyancy shown in Figure 3 (bottom) as average water depths d representative of average basal water pressure are located in this spectrum as shown.

Figure 6

Fig. 6 Four possible cases of basal water configurations shown as average d variations in Figure 3 (bottom) and average ø variations in Figure 5. Ice-stream beds are shown in plan view (top) and transverse cross-section (bottom), with thin basal water films for which ø = 0 shown in white, thick basal water layers for which ϕ=0 shown in black, and water-saturated till or sediment for which ø = øG shown as dotted. The vertical is exag-gerated.

Figure 7

Fig. 7 A basal water configuration in plan view and longitudinal cross-section for the ø variation of case I through case IV in Figure 5 that is different from the basal water configuration for these cases shown in Figure. 6. The bed consists of tills and hollows instead of ridges and valleus, with ø = 0 in white areas, ø = (h/hG in connected black areas, and ø = øG m disconnected black areas, for which øG = 0. In case IV, basal water is trapped in a glacially eroded trough dammed by glacially deposited sediment till(dotted) at the grounding line, situationsfound in the inter-island channels and fiords of deglaciated landscapes. The vertical is exaggerated.

Figure 8

Fig. 8 Possible flow curves for lateral shear in an ice stream. Strain hardening in the zone of converging flow at the head of an ice stream causes an increase in side-shear stress Ts and transverse shear strain τs until mscoplastic yie.ld stress TV is reached. Over length Ls of stream flow, Ts is constant if strain hardening and softening rates are in balance (cruve 1), Ts decreases if strain softening dominates (cruve 2 and 3) and rs drops to zero if strain softening leads to shear rupture (curue 3).

Figure 9

Fig. 9 A method for approximating the side area of an ice stream that uses successive straight lines to represent the concave top surface, with one line for each additional Δx inrecrement added to x. The less concave the surface, the better the approximation.

Figure 10

Table. 2 A life-cycle classification for ice streams*

Figure 11

Fig. 10 The variation of ice-bed coupling ø/øa with distance x upstream /rom the grounding line of a marine ice stream of length Ls for c values ranging from zero to infinity in Equation (58).

Figure 12

Table. 3 Flowline. surface slopes above a horizontal bed

Figure 13

Fig. 11 Evolution of an ice flowline profile from convex sheet flow to concave stream flow as ice-bed coupling exponent c in Equation(58) decreases from infinity to zero for intermediate ice-shelf buttressing given by øG = 0.5.

Figure 14

Fig. 12 The lengthening reach of a marine ice stream into an ice sheet during the inception and growth stages, shown as flow-line profiles for various c values in Equation (58) when øaG = 1 specifies no ice-shelf buttressing. As c increases, the concave stream-flow profile (excluding dashed parts of surface-elevation curves) migrates toward the ice divide, until substantial down-draw of sheet flow occurring after a surface-inflection instability at 0.47 < c < 0.48. Stream flow steadily reverts to sheet flow for c < 1, causing the flow-line elevation to increase and its profile to become increasingly convex.

Figure 15

Fig. 13 Surface profiles of marine ice streams for ice-slielf buttressing increasing from, no buttressing (øG = 1) to fall buttressing (øG = 0) for moderate ice-bed coupling (c = 1 ). The lowest profiles arc for moderate ice-shelf buttressing (0.4 < øG < 0.6).

Figure 16

Fig. 14 Surface profiles of marine ice streams for ice-shelf buttressing increasing from no buttressing (øG = 1) to full buttressing (øG = 0) for low ice-bed coupling (c = 0.5). The concave surface of stream, flow and the surface lowering of sheet flow develop most rapidhj at the onset of ice-shelf buttressing (1.00 < øG < 0.05).

Figure 17

Fig. 15 Pulling power (hiring the life cycle of an ice stream emphasizing the inception stage, with c = 0.1 and øG decreasing from unity to zero. As the ice stream retreats without concomitant down-draw of the ice sheet (top), pulling power is concentrated at the head of the ice stream (bottom). After down-draw begins for øG < 0.90 (top), pulling power becomes greatest toward the foot of the ice stream (bottom).

Figure 18

Fig. 17 Pulling power during the life cycle of an ice stream emphasizing the mature stage, with c = 1.0 and øG decreasing from unity to zero. Pulling power decreases steadily upstream from the grounding line, beginning at lower values as øG decreases,

Figure 19

Fig. 16 Pulling power during the life cycle of an ice stream emphasizing the growth stage, with c = 0.5 and øG decreasing from unity to zero. Pulling power decreases slowly up-stream, from the grounding line, especially for the higher øG values.

Figure 20

Fig. 18 Pulling power during the life, cycle of an ice stream emphasizing the declining stage, with c = 2 and øG decreasing from unity to zero. Pulling power decreases rapidly upstream from the grounding line for larger øG values and is low even at the grounding line for small øG values. The terminal stages in the life cycle of an ice stream are represented by øG = 0 in Figures 15 through 18, for which Px = 0.