Introduction

Calving of glacier cliffs at tidewater termini is the major ablation process, in some cases leading to sustained, rapid retreat and possibly catastrophic collapse of a glacier. In the simplest sense, calving represents the primary negative component of the ice flux at a glacier terminus. The relative rates of calving and ice flow, the primary positive component of terminus ice flux, essentially control the advance and retreat of the ice front. In contrast to ice velocity, however, which can be measured directly on a number of time-scales, calving rate is a calculated rather than measured parameter. Determining the factors controlling calving is consequently reduced to a statistical correlation exercise (e.g. Reference Brown, Meier and PostBrown and others, 1982).

In the past two decades a significant amount of observational data has been collected to constrain ice flux of the terminus (e.g. Reference Brown, Meier and PostBrown and others, 1982; Reference Krimmel and VaughnKrimmel and Vaughn, 1987; Reference Meier and PostMeier and Post, 1987; Reference MeierMeier and others, 1994). The most complete dataset comes from Columbia Glacier, Alaska, and demonstrates that flow of tidewater glaciers responds to both short- and long-term forcing from a variety of sources. Rain and wind can cause increases in ice-flow velocity for periods of hours to clays, and tidal fluctuations in water depth result in systematic, daily variations in glacier velocity.

Brown and others' (1982) statistical analysis of Columbia Glacier data showed a strong linear correlation between calving rate and water depth at the glacier terminus. Despite the good correlation, however, a physical basis for this relationship was lacking until Reference HughesHughes (1992) presented a theoretical calving model that not only yielded a linear relation between water depth and calving rate, but also suggested additional calving-rate dependencies on ice buoyancy and crevasse spacing. Reference Van der VeenVan der Veen (1996) called these relationships into question and demonstrated equally strong correlations between the 2 year running-mean caking rate, mean ice thickness, longitudinal stretching rate and the 2 year running-mean ice speed. Reference Van der VeenVan der Veen (1996) further suggested that a universally applicable caking relation may not exist, and that different calving relations may hold for glaciers with stable vs retreating termini. In the well-documented retreat from the Little Ice Age maximum in Glacier Bay, Alaska, a relationship has been found between average rate of terminus retreat and calving rate. However, the relationship does not hold over short periods (Reference Powell and MolniaPowell, 1983) and is not necessarily related to water depth; and other glaciological factors such as down-wasting, fast flow (e.g. Meier and others, 1985a) and relict crevassing (Reference PowellPowell, 1988) may also be important factors.

Boundary Conditions at the Terminus

The model we invoke follows that suggested by Reference PostPost (1975), Reference Meier and PostMeier and Post (1987), Reference MayoMayo (1988) and Reference AlleyAlley (1991). In this model, a glacier builds a morainal bank at the terminus while in a retreat position. During advance, the glacier pushes the bank ahead of it, with bank size continually increasing due to the addition of sediments. Figure 1 shows a generalized longitudinal cross-section through a grounded glacier and the tidewater terminus with an associated push-morainal bank. In our model, sediment is continually added to the bank from supraglacial, englacial and sub-glacial sources, with a minor additional component provided by scraping-off at the morainal toe (e.g. Hunter and others, 1996a). Sedimentation rates increase exponentially toward the terminus, tending to maintain the overall wedge shape of the morainal bank (Reference Cowan, Powell, Anderson and AshleyCowan and Powell, 1991). Sediment is removed from the bank when the glacier overrides its push-morainal bank. Submarine slumping and mass wasting redistribute sediment in the forebank area, but are not significant mechanisms for sediment removal.

The geometry of glacier termini against the backs of morainal banks appears to vary in known cases. Examples in Alaska at present termini (Reference Cowan, Powell, Anderson and AshleyPowell, 1991; Reference Hunter, Powell and SmithHunter and others, 1996a) and from seismic reflection records of recently abandoned banks (Reference Seramur, Powell, Carlson, Cowan and DaviesSeramur and others, 1997; Reference Cai, Powell, Cowan and CarlsonCai and others, 1997), as well as those from Pleistocene morainal bank exposures (e.g. Reference Hunter, Powell and SmithHunter and others, 1996b), appear to show a range in geometries from a steep, near-vertical face to a ramp-type form. In some eases there even appear to be buried ice bodies in the bank. However, the variability in form is not well documented, and causes of the variability are not yet understood. In this first attempt, we take the simplest case for an idealized model, that of a vertical contact between glacial ice and the morainal bank. This represents an end-member category and provides for the maximum back-stress case of the bank on the glacier. As bank backs become more ramp-like, the back-stress effect should decrease, and those cases require further modeling. Using the simplest geometry for the ice bank interface, we note that depending on the relative motion of the ice front over time, three general sets of boundary conditions are possible at the glacier terminus.

Advancing glacier

An advancing glacier actively pushes a bank of sediment in front of it, analogous to the pile of sediment pushed in front of a bulldozer. When comprised of uniform, non-cohesive materials, the sediment bank will uniformly taper away from the terminus in the bank-front section as shown in the idealized model in Figure 2. As discussed by Reference Davis, Suppe and DahlenDavis and others (1983) and Reference DahlenDahlen (1984), such "critically tapered wedges" assume a taper determined by several factors, including the frictional characteristics of the constituent particles and the sliding conditions at the base of the wedge. For our model, the actual dimensions of the bank-front wedge are also important, as the proportions of the wedge will be related to the volume of sediment in the whole bank. As sediment accumulates at the terminus, the morainal bank will grow self-similarly (e.g. Turcotte, 1992) at a constant taper (Reference Davis, Suppe and DahlenDavis and others, 1983). Although the type and facies distribution of sediments in the wedge may be significant in determining the geometry of the bank-front wedge, we simplistically treat the bank as a uniform mass of non-cohesive sediment.

Fig. 1. Schematic diagram of the natural, physical system at the grounded tidewater terminus of a temperate glacier

Fig. 2. Idealized model of a grounded tidewater terminus during a prolonged phase of glacier advance.

Retreating glacier

A retreating glacier creates a gap between the cliff and the bank crest. During retreat, the bank may assume different geometries depending on the relative rates of retreat and sedimentation (Fig. 3). If sedimentation cannot keep pace with retreat, the glacier may retreat completely from the bank, leaving an asymmetric form with a gentle foreslope and a backslope near the angle of repose (Fig. 3c). Observations of morainal bank geometry suggest this may be the most common case (Reference Cowan, Powell, Anderson and AshleyPowell, 1991; Reference Hunter, Powell and SmithHunter and others, 1996a).

Fig. 3. Idealized models of morainal bank geometry during prolonged glacier retreat.

Terminus Restraining Forces

Our analysis considers only (he forces at the glacier terminus that resist flow of the glacier. We use the simple model in Figures 4 and 5 to formulate the problem, and discuss three principal restraining forces (see Appendix for symbols used): a horizontal normal force due to the water (F _w), a horizontal normal force due to the morainal bank (F _M), and a horizontal shear force associated with frictional sliding along the base of the forebank wedge F>_BS Because different combinations of restraining forces operate under each set of boundary conditions, we first discuss the forces common to all sets of boundary conditions, and then discuss restraining forces associated with specific cases. We do not discuss the forces associated with long-term stationary termini, as these conditions are potentially rare. In all cases, we assume the morainal bank sediments are saturated with water at hydrostatic pressure.

Fig. 4. Geometry, coordinate system and important parameters in the model. Coordinate system is anchored to the toe of the glacier and moves with it during advance or retreat.

Fig. 5. Restraining forces associated with grounding-line push-morainal banks.

Hydrostatic pressure

The hydrostatic fluid pressure (P_f at any vertical position in the system is given by

(1)

where ρ_w is the density of the water, g is the acceleration of gravity, d_max is the water depth beyond the bank foreslope, and the equation is valid for all z P d _max. At the vertical ice-water interface, where the water depth is D the hydrostatic restraining force F _W> per unit width ai the glacier terminus is given by

(2)

This hydrostatic restraining force is operative during all stages of glacier advance and retreat, as long as the hank remains submerged.

Shear force from basal sliding

As a glacier advances it pushes a bank of sediment in front of it. As the forebank wedge is transported, it overrides finer-grained, more distal sediments on the sea floor. Some of these sediments are scraped from the sea floor and accreted to the toe or underside of the sedimentary forebank wedge. Those sediments not accreted serve as a weak detachment layer beneath the forebank wedge. Evidence of such detachment layers, including high vertical-strain gradients and the presence of narrow, subhorizontal shear zones, is seen in many subareally developed push moraines (e.g. Van der Wateren, 1986,1994). Any time a push moraine advances, a shear force related to basal sliding of the moraine operates to resist glacier flow.

The basal sliding shear force per unit width of the morainal bank is given by

(3)

where τ _b,(x) is the distribution of shear stress along the base of a sediment wedge of horizontal length L (Fig. 4). For a uniform, non-cohesive wedge of sediment,

(4)

where μ _b is the coefficient of internal friction (Reference Jaeger and CookJaeger and Cook, 1979) of the sediments comprising the basal detachment layer, and σ_zz ′(x, z = 0) is the distribution of normal effective stress along the base of the forebank wedge. To fully describe F _BS we must explicitly know this normal stress distribution along the base of the wedge, a quantity we now derive.

From the schematic vertical stress profile in Figure 6, it is clear that the total vertical stress at any point in the bank consists of two components. Consequently, vertical effective stress at any point in the bank can be written as

(5)

Seeing from Figure 4 that the upper surface of the forebank wedge is defined by

(6)

the distribution of vertical stress along the top of the bank is

(7)

and anywhere in the wedge, the vertical stress due only to the weight of the sediment is

(8)

where p _s is the mean density of the sediments in the bank. Substituting Equations (1), (7) and (8) into Equation (5) yields the equation for vertical effective stress anywhere in the bank:

(9)

Substituting z = 0, simplifying and rearranging yields the equation for vertical effective stress, at the base of the morainal bank,

(10)

The shear force acting along the base of the bank is now found by substituting Equation (10) into Equations (4) and (3), and evaluating the integral to obtain:

(11)

Fig. 6. Variation of total vertical stress (σ_zz) along any vertical profile (i.e. x is constant) through the morainal bank. Total vertical stress at any point in the forebank sediment wedge is the sum of the weight of the water column acting on top of the wedge (σ_zz ^w(x)),plus the weight of the sediment column (σ_zz ^s(x, z)) above the point in question.

Horizontal normal force during advance

Reference Davis, Suppe and DahlenDavis and others (1983) demonstrated that a non-cohesive mass of material will naturally assume a prismatic wedge geometry when pushed by a rigid boundary. Similar conclusions were drawn by Reference ElliottElliott (1976) and Reference ChappleChapple (1978) for cohesive, ductile materials. Although derived to describe the geometry and mechanics of thin-skinned fold and thrust belts and accretionary prisms along subducting plate boundaries, these "critical wedge" models are also applicable to moraines pushed before an advancing glacier (e.g. Van der Wateren, 1994). The following relies extensively on the work of Reference Davis, Suppe and DahlenDavis and others (1983) and Reference DahlenDahlen (1984), to which the reader is referred for details of non-cohesive critical-wedge mechanics.

According to Reference Davis, Suppe and DahlenDavis and others (1983), a "critically tapered" wedge is the "thinnest body that ran be thrust over its basal décollement without any internal deformation". When new material is added to a forebank wedge, by either sedimentation, underplating or accretion at the toe, the wedge "deforms internally while sliding in order to accommodate the influx and to maintain its critical taper" (Reference Davis, Suppe and DahlenDavis and others, 1983). The sediment in the morainal bank is therefore considered to be in a "critical state", because it is continually in a state of failure everywhere internally and along the base.

The critical state of stress in a bank composed of uniform, cohesionless material is illustrated by the Mohr circle in Figure 7. Shear and normal effective stress along the wedge base are given respectively by τT _b(x) and σ _zz′{x,z = 0). These two stresses are related by the criterion for frictional sliding along the wedge base: a linear Coulomb shear failure envelope with slope tan (oslash; _b = μb, where μb is the coefficient of friction along the base of the forebank wedge. Shear and normal effective stresses along failure planes in the wedge interior, τT _i{x) and σ _ni′(x, z),are related by a second linear Coulomb shear failure criterion with slope tan (oslash; _i = μ _i where μ _i is the coefficient of friction of sediment comprising the wedge.

Because of frictional sliding along the wedge base, principal stress orientations in the forebank wedge are inclined at an angle ψ to the x-z coordinate system describing the wedge (Fig. 8). For uniform, cohesionless materials, this mis-orientation angle is constant throughout the wedge (Reference DahlenDahlen, 1984), and is related to the wedge taper, material density and fluid pressure ratio in the wedge (λΛ = P _f/σ_zz Reference Hubbert and RubeyHubbert and Rubey, 1959) by

(12)

where α′ is Reference DahlenDahlen's (1984) "modified slope angle" defined by

(13)

and α is the forebank wedge taper as shown in Figure 4. The graph in Figure 9 illustrates the relation between sediment density and α′ for wedge tapers of 15° and 20°, and the typical upper and lower bounds to morainal bank foreslope angles measured by Hunter 1994) for three Alaskan glaciers. Figure 10 illustrates the relation between ψ and ø _i for various sediment densities and critical taper angles of 15° and 20°.

Fig. 7. Slate of effective normal and shear stress in an actively advancing, critically tapered forebank wedge of uniform, cohesionless sediment. After Reference Davis, Suppe and DahlenDavis and others (1983).

Fig. 8. Geometry of stresses in a critically tapered forebank wedge of sediment. Principal stresses are inclined at an angle ψ to the wedge base because of frictional sliding. For the simple case of a planar, horizontal base, wedge taper is defined by the angle α.

In the following analysis, we assume a sediment density of 2.2 g cm⁻³, a wedge taper of 15° hydrostatic fluid pressure (λ = 0.465), internal friction angle of 40° and a water density of 1 g cm⁻³. In a morainal bank where these assumptions hold, the stress misorientation angle is ~5° (Fig. 10), and failure within the wedge will be on normal and thrust faults with orientations shown in Figure 8. Although deformation both within the forebank wedge and at the wedge base may be complex, we need only know the distribution of normal stress along the back of the morainal bank in order to determine the horizontal, normal restraining force operating during glacier advance.

Fig. 9. Variation of modified slope angle as a function of bank sediment density.

Fig. 10. Variation of stress misorientation angle as a function of function angle for wedge tapers of 15° and 20° and various sediment densities.

To determine the horizontal normal force associated with pushing of the morainal bank, we utilize the Mohr circle in Figure 7, and note that for a critically tapered wedge,

(14)

After expanding the horizontal effective stress and rearranging, we find that total horizontal stress at any position within the bank is given by

(15)

The horizontal normal force per unit width along the back of an advancing critically tapered forebank wedge (Fig. 4) is therefore found from

(16)

where σ_xx (x, z) is given by Equation (15).

Setting x = 0 and substituting Equations (1), (4) and (9) into Equation (16), we arrive at a complete expression for the horizontal normal force per unit width of the morainal bank during glacier advance:

(17)

Simplifying, rearranging, and evaluating this integral yields:

(18)

Horizontal normal force during retreat

When a glacier retreats, only one bank-related restraining force acts at the glacier terminus. This restraining force is a horizontal normal force (F _M Fig. 5) derived from the wall of sediment laterally supported by the ice cliff. As noted in Figure 3, a variety of bank geometries can develop as a glacier retreats. The normal force associated with each of these geometries is intimately associated with the distribution of horizontal total stress along the ice—sediment interface, such that the horizontal force per unit width of the glacier terminus is given by Equation (16), where σ_xx (x=0,z) is now a new function, yet to be determined.

Because a full analysis of the horizontal stress distribution associated with each of the wide array of possible bank geometries is beyond the scope of this paper, we instead present a solution for σ_xx (x=0,z)when sedimentation rate does not keep pace with glacier retreat i Fig. 3c). We examine this case because it represents a minimum bound to the morainal-bank-derived horizontal normal force acting at the terminus of a retreating glacier.

For a non-cohesive material resting at the angle of repose, the ratio of vertical to horizontal effective stresses is controlled by the frictional properties of the material particles such that

(19)

where μ is the coefficient of friction of the material, and is equal to tan ø where ø is the angle of repose of the material (Reference Zoback and HealyZoback and Healy, 1984). Substituting C for the constant in brackets in Equation (19) above, expanding the horizontal effective stress term σ_xx ′(x, z) = σ_xx (x, z)— P_f (z) and rearranging, yields an equation for total horizontal stress anywhere in a non-cohesive pile of sediment resting at the angle of repose:

(20)

Substituting Equations (9) and (1) into Equation (20), and setting x = 0, we find the distribution of horizontal stress along the ice-morainal-bank interface at the glacier terminus, the grounding line, during retreat to be

(21)

Because we are now only interested in the height of the bank at the grounding line (H; Fig. 4), which is no longer simply related to the dimensions of the entire morainal bank, we can substitute H for L tan α to obtain:

(22)

The horizontal restraining force per unit width of the glacier terminus during retreat is now found by integrating the horizontal stress distribution derived in Equation (22), and is given by:

(23)

Comparison of Restraining Forces

Figure 11 illustrates the behavior of the four restraining forces at glacier termini derived in Equations (2), (11), (18) and (23). Figure 12 illustrates the combined effect of these forces during advance and retreat of a marine-ending glacier with a tidewater terminus. These results suggest that glacial advance is an inherently stable process, where restraining forces continually increase as the morainal bank builds. Rapid glacier advance is therefore likely to be associated, not with a change in restraining forces at the terminus, but rather with some change in restraining forces developed elsewhere in the system (e.g. the base; Reference KambKamb, 1987). In contrast, glacier retreat is an inherently unstable process because terminus restraining forces continually decrease as the glacier retreats from a morainal bank, especially a large terminal bank. By the time the moraine has built to 20-30% of the maximum water depth, the restraining force associated with water depth at the terminus is relatively insignificant, having been surpassed by horizontal normal forces due to the bank.

Fig. 11. Restraining forces acting at a tidewater terminus. Plot shows the force per unit width due to the hydrostatic pressure in the water column (F_W),basal sliding of the bank (F_BS)normal stress along the bank-ice interface during advance (F_MA) and normal stress along the bank-ice interface during retreat (F_MR). Note that F_W decreases as a morainal bank builds, while all other restraining forces increase with increasing height of the bank (H in Fig. 4). Calculation done for a bank that builds from 0 to 100 m height in an area where the maximum water depth is 100m (i.e. d_max = 100). Basal sliding force calculated for a basal coefficient, of friction, μ_b = tan (20°), and normal force during retreat calculated for a coefficient of internal friction, μ_i = tan (40°).

Fig. 12. Restraining forces operating during advance (F_a) and retreat (F_r) of a tidewater terminus. Retreat forces are the sum of F_W, a decreasing function of morainal bank height at the ice cliff {H), and F_Mr, an increasing function of H.

Excluding F _W, all of the restraining forces are dependent on the material properties of the bank constituent sediments, particularly the coefficient of friction (i.e. angle of repose). Although we do not discuss the effects in detail, we note that decreasing basal friction reduces F _BS and F _Mr, whereas decreasing internal friction reduces F _Ma. Frictional properties of unconsolidated sediments are at least related to average particle size, uniformity of particle size, particle composition and angularity of particles (Reference Lambe and WhitmanLambe and Whitman, 1969; Reference Dackombe and GardinerDackombe and Gardiner, 1983), suggesting that the restraining forces associated with any specific glacier will be highly dependent on the sedimentological characteristics of a terminal morainal bank. Full analysis of these dependencies will require a sedimentologically accurate and detailed model of facies distribution and evolution in a bank, and is beyond the scope of this paper.

We note that if there is a relationship between calving rate and water depth, this relationship will be most apparent in the absence of other restraining forces at the terminus. When glaciers are grounded along significant morainal banks, the influence of water in restraint at the terminus will be clouded, and may be insignificant compared to other operative forces.

Glacier Advance, Retreat and Stability

This analysis was designed to evaluate whether calving is a cause of glacier retreat (e.g. Reference Brown, Meier and PostBrown and others, 1982; Reference Meier and PostMeier and Post, 1987; Reference MeierMeier and others, 1994), or merely an effect of increased glacier flow leading to glacier thinning, and its subsequent approach to buoyancy (e.g. Reference Van der VeenVan der Veen, 1996). The inverse statistical correlation between calving rate and water depth at the terminus is presently interpreted in the context of Reference HughesHughes (1992) calving model, where calving is driven by bending creep of the glacier terminus, and bending creep is suppressed by greater water depths at the terminus. As noted by Reference Powell and MolniaPowell (1983, Reference Powell and Molnia1988) and Reference Van der VeenVan der Veen (1996), however, this interpretation cannot explain the sustained rapid retreat of Alaskan glaciers, which have continued to retreat, despite encountering a variety of water depths at the terminus. Our results suggest why temperate glaciers with grounded tidewater termini may exhibit such rapid retreats.

Seasonal terminus fluctuations have been described for Alaskan tidewater termini for many years, but the best-documented has been that of Columbia Glacier. Reference Meier, Rasmussen and MillerMeier and others (1985b) quantitatively showed that the Columbia Glacier terminus regularly advances between October/November (late autumn/winter) and May/June (spring/early summer) and retreats for the rest of the year. They also showed glacier speed increases synchronously with retreat rate and decreases synchronously with advance rate of the terminus. These regular, small-scale, seasonal fluctuations in glacier motion occur even during the present phase of unstable retreat. Our analysis indicates that bank back-stress may be an important factor in such seasonal terminus fluctuations.

Due partly to the influence of terminus restraining forces (i.e. F _MA, F _M, F _BS; Fig. 4), glacier flow, calving rate and terminus ice flux decrease and increase cyclically, perhaps resulting in, or at least amplifying, seasonal advance-retreat oscillations. Beginning in late winter and continuing through spring, flow velocities are greatest, with the glacier advancing, thinning, and actively building a morainal bank (Fig. 13a). Steepening of the morainal-bank fronts during winter and spring has been attributed to these advances (Hunter and others, 1996a), As the bank builds, restraining forces continually increase (Fig. 13d), which assists in eventually halting net forward motion, and enhancing the onset of glacier thickening behind a large bank through summer and early fall (Fig. 13b). The seasonal glaciological cause of switching this process is not well constrained. However, perhaps a response to previous thinning and an approach toward buoyancy forces an increase in calving subsequently leading to glacial retreat. By late fall and through early winter, the glacier has retreated partly off the terminal morainal bank (Fig. 13c), restraining forces are at a minimum (Fig. 13d) and a rapid acceleration of flow begins. Because the glacier is now at its thickest yearly average, gravitational driving stresses are large, buoyancy is greatly reduced and calving decreases (Fig. 13c).

The cyclic variation in restraining forces at the terminus has several implications for the three general conditions of long-term glacial movement: advance, retreat and quasi-static. A quasi-static terminus position is possible if short-term advance and retreat distances are similar. Long-term, steady-state advance occurs if the distance covered by short-term advance cycles, on average, is greater than the average distance covered by short-term retreat cycles (Fig. 14a). Long-term, steady-state retreat occurs if the distance covered by short-term advance cycles, on average, is less than the average distance covered by short-term retreat cycles. Note that long-term advance is a steady-state, oscillatory process because restraining forces are continually increasing as the glacier advances. Unstable, long-term, rapid advance cannot occur in this situation. This does not apply to surging, however, which is a short-term response to a unique set of conditions wherein ice flux is so high that losses due to calving are overwhelmed. As noted previously, surging behavior is most likely brought about by a change in conditions at the glacier bed, rather than the terminus. Comparatively, although long-term unstable advance is not likely, unstable, long-term retreat is possible, and probably more likely than steady-state retreat, for reasons we now discuss.

Processes which result in glacial retreat, namely, thinning, buoyancy and calving, also result in decreased restraint of the terminus, thereby triggering advance. Therefore, retreat tends to be naturally halted, and steady-state short-term oscillatory behavior is preserved. However, this behavior is possible only in the presence of significant moraine-related restraining forces. These forces not only resist advance, but ultimately induce thickening, reduce buoyancy and therefore inhibit processes that cause retreat. If, during short-term advance, the bank-related restraining forces are not increased to some threshold value (e.g. the bank is not built to some threshold height), the glacier will be initially thinner and more buoyant at the onset of the next short-term retreat cycle. Should this occur, the following seasonal cycle of retreat may reduce restraining forces below some critical level, allowing rapid retreat to begin (Fig. 14b).

This critical level of restraint at the terminus represents a significant instability in the short-term advance-retreat cycle that may be triggered by one or more larger, successive retreat cycles. During such a cycle, the glacier may retreat significantly from the terminal morainal bank, and does not re-establish considerable contact with the bank during the following advance cycle. Unless sedimentation rates are very high, short-term advances behind the terminal morainal bank do not produce mechanically significant morainal banks (i.e. the banks are much smaller), resulting in continually decreasing restraint at the terminus at the onset of each new seasonal advance cycle. In this situation, the glacier never slows and thickens appreciably, and enters a period of rapid, sustained retreat. Because significant restraint at the terminus is not generated during advance behind the terminal morainal bank, steady-state retreat is unlikely.

Although the specifies of individual glacial systems will vary, our analysis generally suggests that quasi-static or steady-state advancing glaciers require restraining forces at the terminus of some threshold magnitude (i.e. terminal morainal banks of some threshold size). We are not presently-capable of determining this threshold level of restraining force at the terminus, but note that it will at least be related to bank height, sedimentation rates, sediment types and the evolution of morainal bank geometry. Glaciers with smaller restraint at the terminus will naturally tend toward unstable retreat. Moreover, rapid retreat can be triggered in any system by one or more successive, large, seasonal retreats in which contact with the terminal morainal bank is lost. The cause of such larger short-term retreats is unknown.

Conclusions

We have quantified the restraining forces at the tidewater terminus of a temperate glacier building a bank of sediment at the grounding line. Restraining horizontal normal forces arise from the hydrostatic water pressure at the submerged glacier terminus, and from the support of the sediments in the bank by the ice cliff. Restraining shear forces operate when the glacier is actively pushing the morainal bank in front of it. Although the simple model we present suffers from inaccuracies in the depiction of bank and terminus geometry, and does not include aspects of sediment distribution or sedimentation rate (e.g. if the glacier flows up the back-side of the bank, resistive forces of the bank on the glacier decrease), several results are noteworthy.

Fig. 13. Variation of flow rate, glacier thickness, morainal bank geometry and calving rate during a seasonal cycle of advance and retreat. (d) shows the idealized, schematic, cyclic variation in restraining forces at the terminus during stages of the cycle illustrated in (a), (b) and (c). Note that the linear, vertical segments of the cycle in (d) are not constrained, and that they are most likely linear, but slightly inclined to the right. For this paper, this deviation from vertical is slight enough to be disregarded.

Fig. 14. Schematic variation in short-term restraining forces at the terminus during long-term, steady-state advance (a) and unstable retreat (b).

During advance, restraining forces increase at an increasing rate. Glacial advance is thus a stable process, governed by a negative feedback, eventually leading to a complete cessation of advance. During retreat, restraining forces decrease at a decreasing rate. This behavior suggests retreat is a positive feedback, and an inherently unstable process. The instability in short-term seasonal advance-retreat cycle is moderated by the presence of a morainal bank, which supplies a significant restraining force to the glacier terminus. Sustained, rapid retreat can only occur if bank-related restraining forces diminish below some threshold value during the retreat phase of a normal oscillatory cycle.

Our results are consistent with the idea that calving is a response to changing glacier conditions, not a cause of increased glacier retreat. The controls on glacial advance and retreat are ultimately tied to glacier flow, thinning and restraining forces at the terminus. Rapid How leads to thinning, which in turn leads to increased calving. Flow rates are closely linked to restraining forces at the terminus. which we show are intimately associated with morainal banks.