Hostname: page-component-cd9895bd7-lnqnp Total loading time: 0 Render date: 2024-12-22T06:09:35.894Z Has data issue: false hasContentIssue false

Sediment, glaciohydraulic supercooling, and fast glacier flow

Published online by Cambridge University Press:  14 September 2017

Richard B. Alley
Affiliation:
EMSEnvironment Institute and Department of Geosciences, The Pennsylvania State University, University Park, PA 16802-7501, U.S.A. E-mail: ralley@essc.psu.edu
Daniel E. Lawson
Affiliation:
U.S. Army Cold Regions Research and Engineering Laboratory, 72 Lyme Road, Hanover, NH 03755-1290, U.S.A.
Edward B. Evenson
Affiliation:
Department of Earth and Environmental Sciences, Lehigh University, Bethlehem, PA 18015, U.S.A.
Grahame J. Larson
Affiliation:
Department of Geological Sciences, Michigan State University, East Lansing, MI 48824, U.S.A.
Rights & Permissions [Opens in a new window]

Abstract

Glaciers often advance over proglacial sediments, which then may enhance basal motion. For glaciers with abundant meltwater, thermodynamic considerations indicate that the sediment–ice contact in the direction of ice flow tends toward an angle opposed to and somewhat steeper than the surface slope (by slightly more than 50%). A simple model based on this hypothesis yields the extent of over-ridden sediments as a function of sediment thickness and strength, a result that may be useful in guiding additional fieldwork for hypothesis testing. Sediment-floored as well as rock-floored overdeepenings are common features along glacier flow paths and are expected based on theories of glacier erosion, entrainment, transport and deposition.

Type
Research Article
Copyright
Copyright © International Glaciological Society 2003

Introduction

Fast glacier flow is favored by subglacial sediment, which can bury bedrock roughness to increase sliding (Reference WeertmanWeertman, 1964), deform to increase ice velocity (e.g. Reference Blankenship, Rooney, Bentley and AlleyBlankenship and others, 1987; Reference AlleyAlley and others, 1989; Reference Clark, Alley and PollardClark and others, 1999; Reference Tulaczyk, Kamb and EngelhardtTulaczyk and others, 2000; Reference Boulton, Dobbie and ZatsepinBoulton and others, 2001) and constrict basal channels to raise basal water pressure, speeding sliding and deformation (Reference AlleyAlley,1989;Reference Walder and FowlerWalder and Fowler, 1994). Subglacial bedrock erosion usually is not expected to produce thick, continuous till sheets (Reference Cuffey and AlleyCuffey and Alley, 1996; Reference Alley, Cuffey, Evenson, Strasser, Lawson and LarsonAlley and others, 1997; cf. Reference Tulaczyk, Scherer and ClarkTulaczyk and others, 2001); thus, over-riding of proglacial sediments is important.

Here, we hypothesize that beneath glaciers with abundant meltwater flow, the slope of the up-glacier side of over-ridden sediments evolves towards an angle controlled by the thermodynamics of the system and the strength of the sediments. We show that this hypothesis is consistent with available data on glacial sedimentary processes, and leads to a simple model that can be used to guide additional hypothesis tests. If this hypothesis is correct, soft subglacial sediments will be extensive owing to the shallow angle taken by their up-glacier sides.

Proglacial Sediment Accumulations

Several processes combine to cause advancing glaciers to encounter thick proglacial sediment accumulations. Webegin with a brief review of these processes; workers familiar with the glacial sedimentary system may wish to skip this section and the next one.

Reported glacial erosion rates range from essentially zero to probably the fastest sustained values observed for any erosive processes on Earth (e.g. Reference Hallet, Hunter and BogenHallet and others, 1996). Erosion of subfreezing beds (Reference Cuffey, Conway, Hallet, Gades and RaymondCuffey and others, 1999) or thawed beds beneath cold ice typically is slow, but the onset of surface meltwater drainage to the bed greatly increases the number and vigor of subglacial erosive processes (e.g. Reference LawsonLawson, 1986; Reference IversonIverson, 1991; Reference HalletHallet, 1996; Reference Alley, Cuffey, Evenson, Strasser, Lawson and LarsonAlley and others, 1997). Despite millions or tens of millions of years of erosion, the melt-water-poor Antarctic ice sheet still rests on upper-crustal rocks in many places (e.g. Reference ten Brink, Hackney, Bannister, Stern and Makovskyten Brink and others, 1997), indicating orders-of-magnitude slower erosion than for, say, Alaskan tidewater glaciers with their abundant meltwater (Reference Hallet, Hunter and BogenHallet and others, 1996).

Valley glaciers and many outlet glaciers characteristically deepen valley bottoms below levels that would have been produced fluvially, giving longitudinal profiles that are less steep and side-walls that are steeper than for fluvial systems (e.g. Reference Brocklehurst and WhippleBrocklehurst and Whipple, 2002), as shown by, for example, Norwegian fjords, Yosemite Valley, California, U.S.A., and the Laurentian Great Lakes (Reference Larson and SchaetzlLarson and Schaetzl, 2001). Some workers may call any such valley “overdeepened”; however, here we restrict the term “overdeepening” to a glacier bed depression in which the elevation of the down-glacier side (the “adverse slope”) increases in the direction of ice flow. We include among overdeepenings those depressions closed by the up-glacier sides of morainal banks or shoals of tidewater glaciers.

For those glaciers with drainage of surface meltwater to the bed, Reference HookeHooke (1991; also see Reference Alley, Strasser, Lawson, Evenson, Larson, Mickelson and AttigAlley andothers, 1999) argued that overdeepenings are favored by the processes of bedrock erosion. As reviewed below, even faster subglacial erosion of sediments also favors overdeepenings. We surmise that overdeepenings are not so much“accidents” as “steady-state” features toward which the glacial landscape evolves.

Glacial production of steep-walled valleys with nearly flat to overdeepened floors has numerous well-known implications for the glacial and glacial–geomorphicsystems. In steep mountains with bedrock subject to landslides, the hypothesis of perfectly plastic hill-slopes (Reference BurbankBurbank and others, 1996) implies that any downwasting or retreat of glacier ice in excess of sediment replacement will tend to destabilize oversteepened slopes. The resulting landslides can extend upslope to reduce the elevations of ridge-lines and peaks (Reference Meigs and SauberMeigs and Sauber, 2000), producing large sediment fluxes (e.g. Reference Warren and KirkbrideWarren and Kirkbride, 1998; Reference Meigs and SauberMeigs and Sauber, 2000). Glacial retreat may lower base level for streams entering a deglaciated valley, releasing more sediment into the proglacial environment (Reference Meigs and SauberMeigs and Sauber, 2000; personal communication from A. Meigs, 2002). Such behavior may contribute to the observation of Reference Hallet, Hunter and BogenHallet and others (1996) that recent sediment yields from mountain glaciers often have exceeded longer-term averages (which themselves are higher than for equivalent fluvial systems) (cf. Reference Meigs and SauberMeigs and Sauber, 2000; Reference Koppes and HalletKoppes and Hallet, 2002). This implies that sediment fluxes may be higher during glacial retreat than during advance.

The reduced or reversed long-profileslopes of deglaciated valleys often serve as sediment traps, reducing sediment flux to regions farther downstream (e.g. Reference Carter, Neil and McCaveCarter and others, 2000) at times when sediment production may be especially fast upstream of the traps. Resulting sediment accumulations in front of retreating glaciers can reach hundreds of meters or more in thickness (e.g. Reference BarnesBarnes, 1987; Reference Nolan, Motyka, Echelmeyer and TrabantNolan and others, 1995; Reference BjörnssonBjörnsson, 1996;Reference Warren and KirkbrideWarren and Kirkbride, 1998). Sediment may accumulate on newly established fluvial surfaces, and may be finer-grained and thus better at lubricating subsequent glacier advances if accumulation is in fjords following marine incursion or in lakes formed in overdeepenings. Isolation of ice in an overdeepening by downwasting or backstepping of active flow, followed by sediment burial slowing ablation (e.g. Reference ChinnChinn, 1996), migration of the locus of proglacial-stream activity, and subsequent formation of a kettle lake, may contribute to trapping of finer-grained sediment.

Glaciers have rather limited ability to push material in front of them, based on models and observations (e.g.Reference Mickelson, Clayton, Fullerton, Borns and WrightMickelson and others, 1983; Reference BennettBennett, 2001). A stress concentration in ice from interaction with an obstacle occurs over distances comparable to the obstacle size (e.g. Reference WeertmanWeertman, 1964; Reference Kamb and EchelmeyerKamb and Echelmeyer, 1986), so a sediment pile larger than a few tens of meters with typical sediment strength of tenths of a bar would cause the ice stress to exceed the few-bar level at which rapid yielding or failure occurs. Somewhat more extensive pushing may be possible if sediment occurs in thin sheets or overlies a very weak decollement. A range of push-moraine types is reported in the literature (e.g. Reference BennettBennett, 2001), with many mechanisms of formation (and some disagreement on those mechanisms), but ice-pushing sufficient to keep many natural proglacial sediment accumulations in front of a glacier is not expected (e.g. Reference Mickelson, Clayton, Fullerton, Borns and WrightMickelson and others, 1983; Reference BennettBennett, 2001).

Barring certain eventualities (such as formation of highly competent ice-marginal streams), the inability of glaciers to bulldoze large quantities of sediment means that advancing glaciers initially over-ride proglacial materials (e.g. Reference Mickelson, Clayton, Fullerton, Borns and WrightMickelson and others, 1983; Reference Nolan, Motyka, Echelmeyer and TrabantNolan and others, 1995). Over-riding enables subglacial transport processes (subglacial streams, deforming glacier beds, and entrainment of sediment into the ice; reviewed by Reference Alley, Cuffey, Evenson, Strasser, Lawson and LarsonAlley and others, 1997), excavating sediment rapidly. Sediment from subglacial streams and deforming beds is delivered back to the proglacial environment (Reference LawsonLawson,1979), and debris and meltwater from accreted basal ice form hummocky morainal complexes (Reference LawsonLawson,1979).

This recycling of sediment in front of glaciers can lead to very large sediment fluxes; for example, Reference Powell, Dowdeswell and ScoursePowell (1990) reported local sedimentation rates of 30 ma–1 at Riggs Glacier, Alaska, U. S.A. This in turn can lead to filling oflakes or fjords, aggradation of proglacial streams, and an increasingly thick sediment package into which the ice advances.

Supercooling Controls on Moraine Slopes

In previous works (Reference AlleyAlley, 1997; Reference Alley, Strasser, Lawson, Evenson, Larson, Mickelson and AttigAlley and others, 1999) we suggested that supercooling of rising subglacial waters provides an important control on recycling of the over-ridden sediment. Here we develop this idea further.

Many glaciers experience rapid surface melting, and channel the water to their beds and along the beds to the ice front (e.g. Reference Rothlisberger, Lang, Gurnell and ClarkRothlisberger and Lang, 1987; Reference PatersonPaterson, 1994). Typically, basal water flow is driven primarily by the potential gradient arising from the ice-surface slope, from higher-pressure regions beneath the ice to lower-pressure zones at the ice front. (Equations are summarized by Reference Alley, Lawson, Evenson, Strasser and LarsonAlley and others (1998).) Close contact with ice maintains the flowing water near the pressure-melting temperature, which increases with decreasing pressure by approximately 0.001 °C per m reduction in ice thickness. The water flow is nearly non-accelerating, so the work done is dissipated as heat. In many situations, this heat is sufficient to warm the water along the pressure-melting curve and to cause additional melting of adjacent ice; however, when basal water flows up the sufficiently steep adverse slope of an overdeepening, the increase in the potential energy from elevation reduces the heat generation below that needed to maintain the water on the phase boundary. Then, if other heat sources (e.g. geothermal heat) are small, the water will supercool and grow ice (Reference Lawson, Strasser, Evenson, Alley, Larson and ArconeLawson and others, 1998), with the latent heat of freezing warming the water along the pressure-melting curve.

Calculations based on these principles (Reference Rothlisberger, Lang, Gurnell and ClarkRothlisberger and Lang, 1987; Reference HookeHooke, 1991; Reference Alley, Lawson, Evenson, Strasser and LarsonAlley and others, 1998) show that if bed slope and surface slope have opposite signs, water flow will still be directed toward the ice front provided the magnitude of the bed slope is less than about 11 times that of the surface slope (depending in detail on the densities of water and ice). Supercooling is possible if the magnitude of the adverse bed slope is more than about 1.2–1.7 times that of the surface slope; a range is quoted because of uncertainty about the degree to which gas exchange between water and some unspecified reservoir will affect the air saturation of the water and thus the pressure dependence of the melting temperature. For the large water fluxes typical of temperate glaciers, the geothermal and sliding heats usually will be sufficiently small to allow supercooling and freeze-on to occur over appropriate bed slopes, assuming the water remains in intimate contact with ice and so does not reach the lowest points of overdeepenings at a temperature much warmer than the local pressure-melting temperature (Reference Alley, Lawson, Evenson, Strasser and LarsonAlley and others, 1998).

Data from Matanuska Glacier, Alaska, and from other glaciers (e.g. Reference Lawson, Strasser, Evenson, Alley, Larson and ArconeLawson and others, 1998; Reference Evenson, Mickelson and AttigEvenson and others, 1999; Reference RobertsRoberts and others, 2002) show that supercooling does occur beneath glaciers having sufficient drainage of melt-water along a sufficiently steep overdeepening, and that this supercooling causes accretion of debris-rich basal ice. Much evidence additionally shows that this ice growth can largely or completely plug subglacial channels (e.g. Reference Hooke and PohjolaHooke and Pohjola, 1994). At Matanuska Glacier, marginal exposures reveal examples of water channels partially or completely plugged by ice that grew into them (Reference Evenson, Mickelson and AttigEvenson and others, 1999).

On Matanuska Glacier, dye added to a large surface stream entering a moulin emerged from a vigorous vent feeding a proglacial river, but traveled at a rate consistent with flow through a water system with openings averaging just over 10 mm in diameter (Reference Lawson, Strasser, Evenson, Alley, Larson and ArconeLawson and others, 1998). Turbid water discharging from vigorous vents in contact with sediment along the base of Matanuska Glacier almost completely lacks bed load although carrying abundant silt (Reference PearcePearce and others, 2003), consistent with ice having plugged subglacial streams and greatly reduced their competence. Sediment may creep into and affect subglacial drainages as discussed byReference Walder and FowlerWalder and Fowler (1994), but is not modeled to reduce either bed load or water-flow velocities to the very low values observed for Matanuska Glacier. Subglacial streams also may be closed if ice flow carries them sufficiently rapidly against obstacles in the bed (Reference WalderWalder,1982); however, the soft-sediment bed of Matanuska Glacier (Reference Lawson, Strasser, Evenson, Alley, Larson and ArconeLawson and others, 1998) and the unremarkable velocity of the glacier (≈0.2 m d– 1 surface velocity across one overdeepening; Reference Ensminger, Evenson, Alley, Larson, Lawson, Strasser, Mickelson and AttigEnsminger and others, 1999) argue against flow-pinching of channels being important there. Thus, subglacial water channels of Matanuska Glacierprobablyhavebeenpluggedbyicegrowth.

Overdeepenings in Sediment: Hypothesis

An overdeepening canexist whether the adverse slope is bedrock, unconsolidated sediment or both. Drainage of water through subglacial sediments, or in low-gradient channels incised into moraine shoals, would reduce or eliminate supercooling; however, the common observation of turbid water discharging near tidewater-glacier calving fronts, and of streams exiting at the base of ice ending on land, indicates that significant basal water often does flow along, rather than through, glacier beds, aided by the tendency for glaciers to produce low-permeability subglacial diamictons.

Basal sediment transport can occur in basal ice or under the ice in meltwater streams or deforming beds. As reviewed by Reference Hunter, Powell and LawsonHunter and others (1996a, Reference Hunter, Powell and Lawsonb) and Reference Alley, Cuffey, Evenson, Strasser, Lawson and LarsonAlley and others (1997), most of the transport beneath tidewater and terrestrial glaciers with abundant drainage of surface meltwater to subglacial channels will be in those channels even over rather soft tills (Reference Walder and FowlerWalder and Fowler, 1994), unless those channels are clogged by ice growth. The importance of stream transport is indicated, for example, by the commonly quite large sediment fluxes from many glaciers in comparison to the volume of their moraines (e.g. Reference Bell and LaineBell and Laine, 1985; Reference LawsonLawson, 1993) (although basal ice produced from accretion of subglacial water has sufficiently high debris concentrations to produce large ice-marginal moraines; Reference LawsonLawson, 1979; Reference Lawson, Strasser, Evenson, Alley, Larson and ArconeLawson and others, 1998).

Owing to the approximately cubic increase of sediment transport with water flux in a conduit (e.g. Reference Alley, Cuffey, Evenson, Strasser, Lawson and LarsonAlley and others, 1997), a distributed subglacial water system has only a small fraction of the sediment-transport capacity of a channelized system carrying the same total water flux. Any development of supercooling that plugs channels thus will greatly reduce the sediment flux. Even with Reference Walder and FowlerWalder and Fowler (1994) canals in which water pressure increases with water flux, ice growth from supercooling will tend to plug regions of preferred water flow, thus spreading the water flow and reducing its ability to erode and transport sediments.

We therefore expect that a too-steep adverse slope of an overdeepening will become a sediment trap, capturing material delivered by streams from up-glacier as those streams lose their ability to transport that sediment owing to ice plugging the channels. (Trapping will preferentially affect coarser clasts because sufficiently fine material can be carried in the distributed water flow, with interesting implications for sorting of subglacial sediments (Reference PearcePearce and others, 2003).) Such sedimentation in an overdeepening shoulddecrease the magnitude of the adverse slope until supercooling is reduced enough that sediment continuity can be maintained in streams. The alternative, that till deformation, basal freeze-on or other processes will remove sediment supplied by streams and thus allow an adverse slope that maintains a large amount of supercooling, is possible, but the very large sediment fluxes in some subglacial streams suggest that non-stream sediment removal will not be sufficient (e.g. Reference Hunter, Powell and LawsonHunter and others, 1996a, b; Reference Alley, Cuffey, Evenson, Strasser, Lawson and LarsonAlley and others, 1997).

The highly variable sediment load observed at different times for a given stage in many proglacial streams (e.g. Reference LawsonLawson, 1993) shows that subglacial streams often have unfilled sediment-transport capacity. We expect that subglacial streams of a bedrock-bedded glacier overrunning a sediment deposit with a sufficiently gradual slope so that supercooling is not approached will tend to erode that sediment deposit to fill this transport capacity, starting with the sediment first encountered by the streams on the up-glacier side of the deposit. This in turn would lead to steepening of the up-glacier side of the remaining deposit.

Steepening may be favored further by any subglacial till deformation (e.g. Reference Tulaczyk, Scherer and ClarkTulaczyk and others, 2001). Where sediment occurs in discontinuous patches on bedrock, flux by subglacial deformation will be less than in an equivalent situation with continuous sediment. If subglacial deformation occurs and there is a down-glacier transition from a glacier bed with exposed bedrock to one completely covered by unconsolidated sediment, the down-glacier increase in sediment flux will cause sediment loss near the bedrock-to-sediment transition.

Thus, data and theory indicate that glaciers commonly advance into and overrun proglacial sediments and then recycle them. We hypothesize that where abundant surface meltwater is present and reaches the bed, the up-glacier sides of the overrun sediments tend to that angle (slightly more than ∼1.5 times steeper than and opposed to the surface slope) capable of causing some supercooling of subglacial waters, reducing sediment transport somewhat to balance the enhanced sediment flux from till deformation and freeze-on accretion of the overrun sediments.

Overdeepenings in Sediment: Hypothesis Test

Available data support the importance of glaciohydraulic supercooling controlling the up-glacier slope of over-ridden sediment. Several modern glaciers (Reference Lawson, Strasser, Evenson, Alley, Larson and ArconeLawson and others, 1998; Reference Evenson, Mickelson and AttigEvenson and others, 1999; Reference RobertsRoberts and others, 2002; Reference Spedding and EvansSpedding and Evans, 2002) have shown supercooling of subglacial water emerging from overdeepenings. Sediment floors are likely for some of these, and known for Matanuska Glacier (Reference Arcone, Lawson and DelaneyArcone and others, 1995; Reference Lawson, Strasser, Evenson, Alley, Larson and ArconeLawson and others, 1998) and Kvíarjökull, Iceland (Reference Spedding and EvansSpedding and Evans, 2002). At Matanuska Glacier, the occurrence of a sediment-floored overdeepening with basal-ice accretion (Reference Lawson, Strasser, Evenson, Alley, Larson and ArconeLawson and others, 1998), exposures showing conduits plugged by accreted ice (Reference Evenson, Mickelson and AttigEvenson and others, 1999), dye-trace evidence of distributed drainage (Reference Lawson, Strasser, Evenson, Alley, Larson and ArconeLawson and others, 1998), and lack of bed load in basal discharge (Reference PearcePearce and others, 2003) provide much support for the hypothesis, although we as yet lack evidence for time evolution of the sediment slope in the overdeepening.

Muir Glacier, Alaska, experienced a rapid retreat, then stabilized and built a moraine shoal above sea level, primarily from subglacial fluvial discharge forming a fan delta along much of the ice margin. Recent data show that the up-glacier slope of this shoal is about twice as steep as the ice–air surface slope (Fig. 1); this basal slope is sufficient to cause considerable supercooling (cf. Reference Arcone, Lawson, Moran, Delaney, Noon, Stickley and LongstaffArcone and others, 2000). During the time the terminus was becoming terrestrial, turbid water frequently was observed emerging on the surface of the glacier, forming debris-laden ice around vents after short-lived water discharge, with vent sites shifting rapidly over time. In addition, overthrusting of ice above the moraine shoal’s up-glacier side was accompanied by water carrying gravelly sands being forced along thrust planes, indicating the presence of high basal-water pressures. This behavior is fully consistent with our understanding of the effects of supercooling plugging subglacial low-pressure channels to create high-pressure water capable of emerging on the surface (Reference Hooke and PohjolaHooke and Pohjola, 1994), and with supercooling during the upward flow of that water creating the debris-laden ice that grew in vents (Reference Evenson, Mickelson and AttigEvenson and others, 1999; Reference Ensminger, Alley, Evenson, Lawson and LarsonEnsminger and others, 2001).

Fig. 1. Data and interpreted ground-penetrating radar profile of the lower 500 m of Muir Glacier, obtained approximately to ice flow, left to right, in 1996. Interpreted contacts are shown between englacial ice and debris-rich basal ice, and between ice and subglacial sediment. The basal slope is steep enough relative to the surface to cause supercooling of water flowing ice–sediment interface. Modified from Reference Arcone, Lawson, Moran, Delaney, Noon, Stickley and LongstaffArcone and others (2000).

A similar high-pressure basal water system in a region likely experiencing supercooling was observed at Columbia Glacier, Alaska. Before its catastrophic retreat, Columbia Glacier terminated on a steep moraine shoal of its own deposits (Reference Meier and PostMeier and Post, 1987). Boring into the terminal overdeepening then (Reference Humphrey, Kamb, Fahnestock and EngelhardtHumphrey and others, 1993) allowed turbid water to rise 180m fromthe bed along the borehole andexit to the side into an englacial drainage system. Had this water continued along the bed to the ice front, the slopes would have caused supercooling, whereas the most direct englacial path from the point where the water exited the borehole to the ice front would not have involved supercooling (Reference Humphrey, Kamb, Fahnestock and EngelhardtHumphrey and others, 1993). Such a basal water system operating at higher potential than an englacial system is consistent with closure of low-pressure basal channels by ice growth (Reference Hooke and PohjolaHooke and Pohjola, 1994). (Note that the observations of upwelling water were maintained for less than 2 days, perhaps not long enough to allow supercooling in the upwelling water to plug the borehole. Note also that closure of low-pressure channels is not uniquely a sign of supercooling, but might arise by ice flow carrying basal conduits against bedrock obstacles (Reference WalderWalder, 1982) or by creep of soft sediments into channels (Reference AlleyAlley,1989; Reference Walder and FowlerWalder and Fowler, 1994), as noted above.)

Taku Glacier, Alaska, overrode proglacial sediments initially with a low surface slope, and then eroded them at rates of meters per year while building a moraine shoal above sea level (Reference Nolan, Motyka, Echelmeyer and TrabantNolan and others, 1995). Data from Reference Nolan, Motyka, Echelmeyer and TrabantNolan and others (1995, fig. 3 assuming nadir reflections) indicate an average slope of over-ridden sediments of about –1.4 times the surface slope in the reach 400–1600m from the terminus. This approximates the supercooling criterion (–1.2 to –1.7 times the surface slope); slope fluctuations along the flowline fall above as well as below the supercooling criterion. Consistent with this, the glacier exhibits about 2 m of basal ice, with characteristics similar to those of basal ice grown from supercooled water at Matanuska Glacier (personal communication from R. Motyka, 2002).

We still do not have full documentation of glaciers overriding both too-steep and too-flat proglacial sediments, causing interactions between subglacial fluvial transport and supercooling that produce slopes causing some supercooling. Collecting such data is likely to be difficult, and is made more so by scarcity of advancing glaciers today, by loss of numerous overdeepened glaciers during retreat from Ice Age and Little Ice Age maxima, and by non-steady effects associated with recent retreat. Nonetheless, agreement of physical insight with available data from Matanuska, Muir, Taku, Columbia and other glaciers supports our hypothesis that subglacial sediment slopes with abundant meltwater fluxes do tend to that angle causing some supercooling.

We next develop this hypothesis further using a simple perfect-plasticity model (cf. Reference NyeNye, 1951) to aid in additional hypothesistesting. This is a steady-state model. Taku Glacier eroded over-ridden sediments by meters per year (Reference Nolan, Motyka, Echelmeyer and TrabantNolan and others, 1995), indicating an adjustment time on the order of 100 years for hundreds of meters of sediment. It is likely that for many glaciers, post-Little Ice Age adjustment is approaching, but has not reached, steady state; thus, a steady-state model may be instructive, but non-steady modeling would be useful in the future.

Perfect-Plasticity Model of Overridden Sediments

Much work (e.g. Reference Tulaczyk, Kamb and EngelhardtTulaczyk and others, 2000) indicates that a Coulomb-plastic model is appropriate for subglacial sediment deformation. Regardless of the dominant processes controlling sediment flux, if subglacial sediment is continuous or nearly so and can be characterized by a yield strength, then the flow and form of the glacier will be influenced by that yield strength.

Assume that a glacier is over-riding sediment of uniform thickness hbmax having plastic yield strength τ . Assume that the system has reached the steady-state configuration shown in Figure 2, with up-glacier horizontal distance x, ice-surface elevation above the origin at the terminus hs, and depth of the overdeepened bed hb as shown. The surface slope is σs = ∂hs/∂x, and the bed slope is α b = -∂hb /∂ x = 7σs, where 7 falls in the range 1.2–1.7 depending on the air-saturation state of the water. Allow the ice thickness to go to zero at the ice front as shown. Let the basal shear stress be the usual pghσs. (Modification for side drag is made by replacing the quantity p g by pgs wherever the quantity appears in this derivation, in which the shape factor s (Reference NyeNye, 1952) is 0.5 for a semicircular cross-section and 1 for an infinitely wide glacier, and thicknesses are understood to be measured along the center line for ice flow.) Here h = hs + h b as shown in Figure 2, with p the density of ice and g the acceleration of gravity. This yields parabolic profiles with

(1)

(2)

Fig. 2. Geometry used in perfect-plasticity model of subglacial sediment.

The sediment wedge under the glacier extends up-glacier to xmax,with

(3)

Taking 7 ≈ 1.5, the extent of the sediment wedge can be calculated from assumed sediment strengths and thicknesses as shown in Figure 3.

Fig. 3. Lateral extent of over-ridden subglacial sediments as a function of sediment thickness and strength, calculated using perfect-plasticity model.

For stiff, thin sediments (strength – = 1bar = 105 Pa, thickness hbmax = 10 m), the calculated extent is quite small (xmax = 5 m), and certainly falls in the complex ice-marginal zone where these calculations are not especially meaningful, but the result suggests minimal subglacial sediment. For thick and soft sediment (n = 0.01 bar = 103 Pa; h b max = 1000 m), xmax = 5000 km of overrun sediment, sufficiently long that such uniform preglacial conditions are unlikely to have applied, and if they did, a glacier responding to ice-age cycles probably would not have approached steady state; nonetheless, extensive subglacial sediment is suggested. For a wide range of intermediate values, physically plausible yet important lengths of sediment-floored glacier are obtained (e.g. τ = 0.1 bar = 104 Pa, hbmax = 100 m, xmax = 5 km sediment wedge).

Discussion

This simple model suggests a general picture of ice–sediment interactions in terminal regions of advancing or steady glaciers with abundant surface meltwater reaching the glacier bed, in which glaciers over-ride thick proglacial sediments, and supercooling limits the ability of subglacial streams to transport the over-ridden sediments back into the proglacial environment. If this is accurate, several implications follow.

Terminal regions of advancing glaciers are commonly sediment-floored, often for distances that are significant compared to the total glacier length.

Because sediment typically increases the basal velocity (by burying bedrock bumps to smooth the bed for sliding, by allowing additional velocity through sediment deformation, and by creeping into low-pressure water channels and raising basal water pressure (Reference Walder and FowlerWalder and Fowler, 1994) to speed sliding and bed deformation), terminal regions of advancing glaciers often have high velocity and low ice–air surface slopes, and thus sensitive response to climate changes.

As argued by Reference Nolan, Motyka, Echelmeyer and TrabantNolan and others (1995), dynamics of subglacial and proglacial sediment transport then have major effects on the behavior of glaciers. For example, rapid advance over a sediment wedge will be followed by downcutting through subglacial sediments and associated surface lowering.

Sudden increase in the ability of the non-glacial environment to remove proglacial sediment (say, because a tidewater glacier has recycled a moraine shoal to the edge of the continental shelf, where the steeper continental slope increased loss of sediment from the shoal) can expose the front of the ice to environmental influences that were previously walled off by the sediment, increasing melting or calving and perhaps triggering retreat (e.g. Reference AlleyAlley, 1991).

Sediment fluxes from a glacier, hence aggradation of proglacial streams or moraine shoals, often will peak just after rapid glacial advance has over-ridden sediments to produce slopes shallower than required for supercooling and thus to produce maximum subglacial transport, as documented by Reference Nolan, Motyka, Echelmeyer and TrabantNolan and others (1995). Subsequent overdeepening from sediment transport would reach the supercooling threshold and reduce sediment flux.

Glacier fluctuations favor glacier fluctuations. Subglacial processes produce overdeepenings that act as proglacial sediment traps during ice retreat; subsequent readvance over-rides accumulated sediment, favoring the low-elevation, low-slope ice surface and fast ice flow associated with rapid variations in behavior.

These results lead to two additional points, which have been made by many previous workers.

Sediment-floored as well as bedrock-floored overdeepenings, far from being accidents or curiosities, are fundamental elements of the glacier environment that, through supercooling of subglacial water and related sedimentary and hydrologic processes, exert first-order controls on glacier behavior and geomorphic effects.

To really understand glacier flow, one should know glacial and climatic history as well as physics.

Conclusions

Glaciers commonly advance over proglacial sediments that then promote fast glacier flow. We hypothesize that in the presence of abundant meltwater, the ice–sediment contact tends to that angle (slightly greater than ∼ 5 0% steeper than and opposed to the surface slope) causing some supercooling of the water. Available data are consistent with this hypothesis. A simple plasticity model based on this hypothesis yields the extent of the over-ridden sediments as a function of the sediment strength and thickness, and may aid in testing the hypothesis.

Acknowledgements

We thank the U.S. National Science Foundation (including grants OPP 0126187 and 9814774) and the U.S. Army Cold Regions Research and Engineering Laboratory for partial funding, G. Denton and the rest of the U.S. National Oceanic and Atmospheric Administration’s Abrupt Climate Change Panel, T. Chinn, and A. Meigs for stimulating discussions, D. Mickelson, S. Tulaczyk and an anonymous reviewer for helpful suggestions, and our colleagues at Matanuska Glacier for long-standing collaborations.

References

Alley, R. B. 1989.Water-pressure coupling of sliding and bed deformation: I. Water system. J. Glaciol., 35(119),108118.Google Scholar
Alley, R. B. 1991. Sedimentary processes may cause fluctuations of tidewater glaciers. Ann. Glaciol., 15, 119124.Google Scholar
Alley, R. B. 1997. Water, sediment and tidewater glaciers; simplistic review and weakly constrained speculations. Byrd Polar Res. Cent. Rep. 0, 5155.Google Scholar
Alley, R. B., Cuffey, K. M., Evenson, E. B., Strasser, J. C., Lawson, D. E. and Larson, G. J.. 1997. How glaciers entrain and transport basal sediment: physical constraints. Quat. Sci. Rev.,16(9), 10171038.CrossRefGoogle Scholar
Alley, R. B., Lawson, D. E., Evenson, E. B., Strasser, J. C. and Larson, G. J.. 1998. Glaciohydraulic supercooling: a freeze-on mechanism to create stratified, debris-rich basal ice. II. Theory. J. Glaciol., 44(148), 563569.CrossRefGoogle Scholar
Alley, R. B., Strasser, J. C., Lawson, D. E., Evenson, E. B. and Larson, G. J.. 1999. Some glaciological and geological implications of basal-ice accretion in an overdeepening. in Mickelson, D. M. and Attig, J.W., eds. Glacial Processes: Past and Present. Boulder, CO, Geological Society of America, 19. (Special Paper 337.)Google Scholar
Arcone, S.A., Lawson, D. E. and Delaney, A.J..1995. Short-pulse radar wavelet recovery and resolution of dielectric contrasts within englacial and basal ice of Matanuska Glacier, Alaska, U.S.A. J. Glaciol, 41(137), 6886.CrossRefGoogle Scholar
Arcone, S. A., Lawson, D. E., Moran, M. and Delaney, A.J. 2000. 12–100-MHz profiles of ice depth and stratigraphy of three temperate glaciers. In Noon, D., Stickley, G. F. and Longstaff, D., eds. GPR 2000, Eighth International Conference On Ground Penetrating Radar, 23–26May 2000, Gold Coast, Australia. Bellingham, WA, International Society of Photo-optical Instrumentation Engineers, 377382. (SPIE Proceedings 4084.)Google Scholar
Barnes, D. F. 1987 Gravity anomaly at a Pleistocene lake bed in NW Alaska interpreted by analogy with Greenland’s Lake Taserssauq and its floating ice tongue. J. Geophys. Res., 92(B9), 89768984.CrossRefGoogle Scholar
Bell, M. and Laine, E. P.. 1985. Erosion of the Laurentide region of North America by glacial and glaciofluvial processes. Quat. Res., 23 (2), 154174 CrossRefGoogle Scholar
Bennett, M. R. 2001. The morphology, structural evolution and significance of push moraines. Earth Sci. Rev., 53(3–4), 197236.Google Scholar
Björnsson, H. 1996. Scales and rates of glacial sediment removal: a 20 km long, 300 m deep trench created beneath Breiðamerkurjokull during the Little Ice Age. Ann. Glaciol., 22,141146.CrossRefGoogle Scholar
Blankenship, D. D., Rooney, C. R. Bentley, S.T. and Alley, R. B. 1987 Till beneath Ice Stream B. 1. Properties derived from seismic travel times. J. Geophys. Res, 92(B9), 89038911.Google Scholar
Boulton, G. S, Dobbie, K. E. and Zatsepin, S.. 2001. Sediment deformation beneath glaciers and its coupling to the subglacial hydraulic system. Quat. Int., 86, 328.CrossRefGoogle Scholar
Brocklehurst, S. H. and Whipple, K. X.. 2002. Glacial erosion and relief production in the eastern Sierra Nevada, California. Geomorphology, 42(1–2)1–24.Google Scholar
Burbank, D. W. and 6 others. 1996. Bedrock incision, rock uplift and threshold hillslopes in the north-western Himalayas. Nature, 379(6565), 505510.Google Scholar
Carter, L., Neil, H. L. and McCave, I. N.. 2000. Glacial to interglacial changes in non-carbonate and carbonate accumulation in the SW Pacific Ocean, New Zealand. Palaeogeogr., Palaeoclimatol, Palaeoecol, 162(3–4) 333356 CrossRefGoogle Scholar
Chinn, T.J. 1996. New Zealand glacier responses to climate change of the past century. N.Z J. Geol. Geophys., 39(3), 415428.CrossRefGoogle Scholar
Clark, P. U., Alley, R. B and Pollard, D.. 1999. Northern Hemisphere ice-sheet influences on global climate change. Science, 286(5442) 11041111.CrossRefGoogle Scholar
Cuffey, K. and Alley, R. B.. 1996. Is erosion by deforming subglacial sediments significant? (Toward till continuity) Ann. Glaciol., 22,1724.CrossRefGoogle Scholar
Cuffey, K. M., Conway, H., Hallet, B, Gades, A. M. and Raymond, C. F. 1999. Interfacial water in polar glaciers and glacier sliding at –17°C. Geophys. Res. Lett., 26(6), 751754.Google Scholar
Ensminger, S. L, Evenson, E. B., Alley, R. B., Larson, G.J., Lawson, D. E. and Strasser, J. C.. 1999. Example of the dependence of ice motion on subglacial drainage system evolution: Matanuska Glacier, Alaska, United States. In Mickelson, D. M. and Attig, J. W, eds. Glacial Processes: Past and Present. Boulder, CO, Geological Society of America, 1122. (Special Paper 337)Google Scholar
Ensminger, S. L., Alley, R. B., Evenson, E. B., Lawson, D. E. and Larson, G. J.. 2001. Basal-crevasse-fill origin of laminated debris bands at Matanuska Glacier, Alaska, U.S.A. J. Glaciol, 47(158), 412422.CrossRefGoogle Scholar
Evenson, E. B. and 6 others. 1999. Field evidence for the recognition of glaciohydraulic supercooling. In Mickelson, D. M. and Attig, J.W., eds. Glacial Processes: Past and Present. Boulder, CO, Geological Society of America, 2335. (Special Paper 337)Google Scholar
Hallet, B. 1996. Glacial quarrying: a simple theoretical model. Ann. Glaciol, 22,18.Google Scholar
Hallet, B, Hunter, L. E. and Bogen, J.. 1996. Rates of erosion and sediment evacuation by glaciers: a review of field data and their implications. Global Planet. Change, 12(1–4) 213235.Google Scholar
Hooke, R. LeB. 1991. Positive feedbacks associated with erosion of glacial cirques and overdeepenings. Geol. Soc. Am. Bull, 103(8), 11041108.2.3.CO;2>CrossRefGoogle Scholar
Hooke, R. LeB. and Pohjola, V A.. 1994. Hydrology of a segment of a glacier situated in an overdeepening, Storglaciaren, Sweden. J. Glaciol, 40(134), 140148.CrossRefGoogle Scholar
Humphrey, N., Kamb, B., Fahnestock, M. and Engelhardt, H.. 1993. Characteristics of the bed of the lower Columbia Glacier, Alaska. J. Geophys. Res.,98(B1), 837846.CrossRefGoogle Scholar
Hunter, L. E, Powell, R. D. and Lawson, D. E. 1996a. Flux of debris transported by ice at three Alaskan tidewater glaciers. J. Glaciol, 42(140), 123135 CrossRefGoogle Scholar
Hunter, L. E, Powell, R. D. and Lawson, D. E. 1996b. Morainal-bank sediment budgets and their influence on the stability of tidewater termini of valley glaciers entering Glacier Bay, Alaska, U.S.A. Ann. Glaciol, 22, 211216.Google Scholar
Iverson, N. R. 1991. Potential effects of subglacial water-pressure fluctuations on quarrying. J. Glaciol., 37(125),2736.CrossRefGoogle Scholar
Kamb, B. and Echelmeyer, K.A.. 1986. Stress-gradient coupling in glacier flow: I. Longitudinal averaging of the influence of ice thickness and surface slope. J. Glaciol., 32(111), 267284.CrossRefGoogle Scholar
Koppes, M. N. and Hallet, B.. 2002. Influence of rapid glacial retreat on the rate of erosion by tidewater glaciers. Geology, 30(1), 4750.Google Scholar
Larson, G. and Schaetzl, R.. 2001. Origin and evolution of the Great Lakes. J. Great Lakes Res., 27(4), 518546.Google Scholar
Lawson, D. E. 1979. Sedimentological analysis of the western terminus region of the Matanuska Glacier, Alaska. Crrel Rep. 79-9.Google Scholar
Lawson, D. E. 1986. Observations on hydraulic and thermal conditions at the bed of Matanuska Glacier, Alaska. In Hydraulic Effects At The Glacier Bed and Related Phenomena, International Workshop, 16–19 September Interlaken, Switzerland. Zürich, ETH. Versuchsanstalt fur Wasserbau, Hydrologie und Glaziologie,6971. (Mitteilungen 90.)Google Scholar
Lawson, D. E. 1993. Glaciohydrologic and glaciohydraulic effects on runoff and sediment yield in glacierized basins. Crrel Monogr. 93-02.Google Scholar
Lawson, D. E., Strasser, J. C., Evenson, E. B., Alley, R. B., Larson, G. J. and Arcone, S. A.. 1998b. Glaciohydraulic supercooling:a freeze-on mechanism to create stratified, debris-rich basal ice. I. Field evidence. J. Glaciol., 44(148), 547562.CrossRefGoogle Scholar
Meier, M. F. and Post, A.. 1987. Fast tidewater glaciers. J. Geophys. Res., 92(B9), 90519058.Google Scholar
Meigs, A. and Sauber, J.. 2000. Southern Alaska as an example of the longterm consequences of mountain building under the influence of glaciers. Quat. Sci. Rev., 19(14–15),15431562.CrossRefGoogle Scholar
Mickelson, D. M., Clayton, L., Fullerton, D. S. and Borns, H.W. Jr. 1983.The Late Wisconsin glacial record of the Laurentide ice sheet in the United States. in Wright, H. E. Jr, ed. Late-Quaternary Environments of The United States. Volume 1. Minneapolis, MN, University of Minnesota Press, 337.Google Scholar
Nolan, M., Motyka, R.J., Echelmeyer, K. and Trabant, D. C.. 1995. Ice-thickness measurements of Taku Glacier, Alaska, U.S.A., and their relevance to its recent behavior. J. Glaciol., 41(139), 541553. (Erratum: 42(141), 1996, p.400.)CrossRefGoogle Scholar
Nye, J.F. 1951. The flow of glaciers and ice-sheets as a problem in plasticity. Proc. R. Soc. London, Ser. A, 207(1091), 554572.Google Scholar
Nye, J. F. 1952. The mechanics of glacier flow. J. Glaciol., 2(12),8293.CrossRefGoogle Scholar
Paterson, W. S. B. 1994 .The Physics of Glaciers. third edition. Oxford, etc., Elsevier.Google Scholar
Pearce, J.T. and 6 others. 2003. Bedload component of glacially discharged sediment: insights from the Matanuska Glacier, Alaska. Geology, 31(1),710.2.0.CO;2>CrossRefGoogle Scholar
Powell, R. D. 1990. Glacimarine processes at grounding-line fans and their growth to ice-contactdeltas. in Dowdeswell, J. A. and Scourse, J. D., eds. Glacimarine Environments: Processes and Sediments. London, Geological Society, 5373. (Special Publication 53.)Google Scholar
Roberts, M. J. and 7 others. 2002. Glaciohydraulic supercooling in Iceland. Geology, 30(5), 439442.Google Scholar
Rothlisberger, H. and Lang, H.. 1987. Glacial hydrology. in Gurnell, A. M. and Clark, M. J., eds. Glacio-Fluvial Sediment Transfer: An Alpine Perspective Chichester, etc., John Wiley and Sons, 207284.Google Scholar
Spedding, N. and Evans, D. J. A.. In press. Sediments and landforms at Kvíarjökull, south-east Iceland: a reappraisal of the glaciated valley landsystem. Sediment. Geol. Google Scholar
ten Brink, U. S., Hackney, R. I., Bannister, S., Stern, T. A. and Makovsky, Y.. 1997. Uplift of the Transantarctic Mountains and the bedrock beneath the East Antarctic ice sheet. J. Geophys. Res.,102(B12), 27,60327,621.Google Scholar
Tulaczyk, S. M., Kamb, B. and Engelhardt, H. F.. 2000. Basal mechanics of Ice Stream B,West Antarctica. I.Till mechanics. J. Geophys. Res.,105(B1), 463481.Google Scholar
Tulaczyk, S. M., Scherer, R. P. and Clark, C. D.. 2001. A ploughing model for the origin of weak tills beneath ice streams: a qualitative treatment. Quat. Int., 86(1), 5970.Google Scholar
Walder, J. S. 1982. Stability of sheet flow of water beneath temperate glaciers and implications for glacier surging. J. Glaciol., 28(99), 273293.Google Scholar
Walder, J. S. and Fowler, A.. 1994. Channelized subglacial drainage over a deformable bed. J. Glaciol., 40(134),315.Google Scholar
Warren, C. R. and Kirkbride, M. P.. 1998. Temperature and bathymetry of ice-contact lakes in Mount Cook National Park, New Zealand. N. Z. J. Geol. Geophys., 41(2),133143.Google Scholar
Weertman, J. 1964 .The theory of glacier sliding. J. Glaciol., 5(39), 287303.Google Scholar
Figure 0

Fig. 1. Data and interpreted ground-penetrating radar profile of the lower 500 m of Muir Glacier, obtained approximately to ice flow, left to right, in 1996. Interpreted contacts are shown between englacial ice and debris-rich basal ice, and between ice and subglacial sediment. The basal slope is steep enough relative to the surface to cause supercooling of water flowing ice–sediment interface. Modified from Arcone and others (2000).

Figure 1

Fig. 2. Geometry used in perfect-plasticity model of subglacial sediment.

Figure 2

Fig. 3. Lateral extent of over-ridden subglacial sediments as a function of sediment thickness and strength, calculated using perfect-plasticity model.