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).
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
The sediment wedge under the glacier extends up-glacier to xmax,with
Taking 7 ≈ 1.5, the extent of the sediment wedge can be calculated from assumed sediment strengths and thicknesses as shown in Figure 3.
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.