Effect of Sedimentation on Ice-Sheet Grounding-Line Stability

See allHide authors and affiliations

Science  30 Mar 2007:
Vol. 315, Issue 5820, pp. 1838-1841
DOI: 10.1126/science.1138396


Sedimentation filling space beneath ice shelves helps to stabilize ice sheets against grounding-line retreat in response to a rise in relative sea level of at least several meters. Recent Antarctic changes thus cannot be attributed to sea-level rise, strengthening earlier interpretations that warming has driven ice-sheet mass loss. Large sea-level rise, such as the ≈100-meter rise at the end of the last ice age, may overwhelm the stabilizing feedback from sedimentation, but smaller sea-level changes are unlikely to have synchronized the behavior of ice sheets in the past.

Ice sheets both cause and respond to sea-level changes. Large, rapid ice-sheet fluctuations have occurred in the past (1) and may recur, but the relative importance of different causes remains uncertain. Here we show that sedimentation beneath ice shelves at grounding lines (where ice begins to float) provides substantial stability against ice-sheet advance or retreat during sea-level changes of up to a few meters or more. Such forcing from sea-level change occurs on the same time scale as does sedimentary stabilization, pointing to other environmental controls as being especially important to ice sheets.

In 1905, R. F. Scott reported evidence of the geologically recent shrinkage of the Antarctic Ice Sheet, despite persistently cold conditions (2). Scott suggested that more snowfall in a warmer past explained the larger ice sheet at that time. Subsequent evidence showing the (near) synchrony of northern and southern ice retreat after the last ice age and that Antarctica had experienced larger ice sheets when the southern climate was colder and drier disproved Scott's hypothesis.

Many workers (3) instead came to accept Penck's argument that sea-level rise from the melting of northern ice sheets had driven the Antarctic ice to retreat (4). Ice-sheet synchronization by smaller sea-level changes linked to Heinrich events or other millennial events has also been frequently proposed (1). However, simple sea-level control is inconsistent with sparse data indicating that, after an initial retreat from the outer shelf (5), additional Antarctic shrinkage was delayed until near the end of the northern-driven deglacial sea-level rise (6). Sea-level control is also inconsistent with the observed slowing of Antarctic ice-stream motion in response to rising tide (7, 8).

Sedimentary stabilization of nonfloating tide-water glaciers is well known (9): Tall, steep-sided sediment bodies deposited between ice and water reduce iceberg calving. Ice flowing from large ice sheets more typically forms floating ice shelves, with sediment being deposited beneath the gradually sloping ice-shelf base. Our model results and the observations of (10) show that such sedimentation also serves to stabilize the grounding line and thus to prevent ice-sheet shrinkage in response to a small (<≈10 m) sea-level rise. Penck's hypothesis (4) may still be valid for ice-sheet response to a larger sea-level rise (≈100 m); however, our results, together with recent evidence that ice shelves respond sensitively to ocean temperature changes and quickly propagate the response inland (11), point to the greater importance of other environmental variables, especially sub–ice-shelf temperature.

Going from the grounded West Antarctic Ice Sheet to the floating Ross Ice Shelf, the upper surface typically drops ≈10 to 25 m in a ramp extending over a few kilometers, which is much steeper than the slope upglacier or downglacier (12). Where data are available, this surface-elevation drop is not caused by thinning of the ice, as might occur from strong basal melting near the grounding line, nor from the effects of basal crevassing or near-surface density changes (12). Instead, the surface-elevation ramp results mainly from a corresponding ramp in basal topography. Thus, the onset of flotation is not caused primarily by a downglacier thinning of ice to the flotation criterion, but more by a downglacier drop in the elevation of the glacier bed. This contributes to the observed stability of the grounding-line position over many decades (12) and to the ice-stream slowdown in response to rising tide (7, 8).

Available data indicate that the topographic step in the glacier bed was formed by the recent deposition of primarily subglacially transported sediment. As shown in (10), a sediment wedge occurs just upglacier of the grounding line of Whillans Ice Stream, and this wedge is remarkably similar to the numerous till-dominated deposits formed at the grounding line of the extended ice sheet during its retreat from the continental-shelf edge as the last ice age ended (13). Observations are available from beneath the grounded portions of two active ice streams (Whillans and Bindschadler) feeding the Ross Ice Shelf, and both showed deformation of soft tills (14); although the deformation is probably discontinuous in space and time (14), notable sediment transport results. The geology near the grounding zone of the Ross Ice Shelf has not been mapped in detail, but available data indicate that the ice flows along fault-controlled basins containing poorly consolidated Tertiary and perhaps Quaternary sediments (15). It is highly unlikely that the sedimentary deposit discovered by Anandakrishnan et al. (10) is a localized feature and that the topographic step is caused by an unknown transverse geologic structure extending across the entire Ross Ice Shelf, except where surveyed by Anandakrishnan et al. (10). Instead, sedimentary control of most or all of the grounding line is highly likely, with the available data supporting wedge formation within the troughs of modern and past ice streams (13).

Using three different ice-flow models [labeled Dt (Dupont), Pk (Parizek), and Pd (Pollard) for their lead developers (11, 16, 17)] including longitudinal as well as basal shear stresses, we have assessed the effect of sub–ice-shelf sedimentation on an ice-stream/ice-shelf system. The models were run for a coupled ice-stream/ice-shelf configuration approximating Whillans Ice Stream and the adjacent Ross Ice Shelf across the grounding-zone wedge. In some runs, ice-shelf buttressing was prescribed, offsetting much of the spreading tendency at the grounding line. After reaching a steady state with a flat bed, the model ice-stream/ice-shelf system was perturbed by instantaneously adding a wedge of sediment similar to the 31-m-high wedge observed by Anandakrishnan et al. (10), with the same basal friction as beneath the preexisting ice stream (18). The evolution to a new steady state was assessed while the geometry of the sediment wedge was held constant. The response of steady-state configurations to instantaneous sea-level change was simulated with and without grounding-zone wedges.

The frictional drag from increased ice/sediment contact with the addition of the wedge slows and thickens the ice above, which in turn causes the geometric effect of the wedge to become important and to further restrict the ice flow (Fig. 1). The resulting grounding-line advance beyond the crest of the wedge causes flotation to occur where the bed falls away, matching observations at the modern grounding line, whereas without the wedge, flotation occurs where the ice thins sufficiently.

Fig. 1.

Effect of grounding-line deposition (Pd model). Dashed lines show the steady ice-stream profile in the absence of a grounding-line wedge, and solid lines show the wedge and the corresponding steady ice-stream profile. The wedge causes the ice to thicken, the grounding line to advance past the wedge crest, and an inflection point to form in the upper surface at the upglacier end of the wedge, which might serve to increase water storage there.

In all of our models, sea-level rise in the absence of a wedge causes the grounding line to retreat from its initial position because the reduction in basal friction from the flotation of previously grounded ice is more important than the increase in back pressure from the deeper water. The wedge causes the ice above the wedge crest to thicken to well above the flotation level; small instantaneous sea-level rise then causes very small grounding-line retreat to a new position that still lies beyond the crest of the wedge and well beyond the no-wedge position. Larger sea-level rise (≈5 m or more) causes the ice to float free of the wedge and to reach the same steady state as in the no-wedge case (Figs. 1 and 2).

Fig. 2.

Grounding-line stability in the Pk (solid lines) and Dt (dashed lines) models. Ice flows from left to right. In the absence of a wedge, the grounding line is at 1.1 km on this scale. The Pk and Dt experiments are not identical because of effects of the underlying model physics on the ice-shelf and wedge geometry, but the experiments are roughly comparable, and the results are clearly qualitatively identical and quantitatively similar. The ability of the grounding line to retreat is greatly limited by the constant-thickness boundary condition at the upglacier (left) end of the model domain; additional experiments with the same boundary condition applied farther away from the grounding line show much greater sensitivity to perturbations when the grounding line is upglacier of the wedge, and they emphasize the stabilizing influence of the wedge. (A) Response to sea-level rise. Adding the wedge moves the grounding line to 10.1 km, and increasing the stiffness of the wedge (the frictional resistance for a specified sliding velocity) by an order of magnitude moves the grounding line to 10.5 km. Sea-level rise causes retreat in all cases. However, the wedge causes the grounding-line retreat to be small up to some limit, beyond which the behavior is identical to the no-wedge case. The stiffer wedge requires a higher sea-level rise to reach that limit. (B) Responsetoa reduction in buttressing. Initial buttressing opposes 69% of the spreading tendency of the ice shelf at the grounding line. A 6% reduction means that buttressing opposes 0.94 × 69% = 65% of the spreading tendency. A reduction of ≈25% is required to cause retreat from the stiffer wedge.

Larger wedges have greater ice thickness above flotation levels at the wedge crest and so require larger sea-level rise to force the grounding line upglacier of the wedge crest. However, as shown in Fig. 3, the dependence of the resistance on the wedge height or length is steeper for smaller wedges than for larger ones. Furthermore, if the wedge volume grows at a constant rate, the wedge height and length will increase more rapidly at first and more slowly later. Both effects combine to produce a more rapid increase in resistance at early stages of growth than at later stages (19).

Fig. 3.

Effect of wedge height on resistance to grounding-line retreat (blue) and on the time required for wedge deposition (red) in the Pd model. In this case, simulations were run for wedges having an upglacier slope of 0.001, a downglacier slope of 0.01, and heights of 0, 4, 8, 12, and 16 m with corresponding lengths of 0, 4, 8, 12, and 16 km, respectively. The thickness of ice above the wedge crest, which provides the resistance to grounding-line retreat, increases with wedge size; however, the increase in ice height above the flotation level is slower than the increase in wedge height, so earlier increments of wedge growth are more effective at stabilizing the grounding line, especially as successive height increments take longer to form if the sediment supply is constant. (The time depends on the sediment flux and so is plotted in relative units.) In comparison, the effect of doubling the stiffness of an 8-km wedge is also indicated. Similar experiments in the Dt model yield comparable trends but somewhat smaller height above the flotation level.

Stiffer wedges (those with greater frictional resistance for a given basal velocity) also require larger sea-level rise to force the grounding line upglacier of the wedge crest (Fig. 2). We expect that at least portions of a wedge will often be stiffer than the surrounding no-wedge regions. The steepened ice/air surface slope that develops in response to friction from a wedge (Figs. 1 and 4) will steepen the subglacial hydrologic-potential gradient and speed the flow of subglacial water, reducing water storage and pressure and thus reducing the lubrication of ice motion by sliding or till deformation (20). Consideration of our numerical experiments with stiffer wedges and of the additional stabilizing effect that would occur from sedimentation during sea-level rise at ordinary rates suggests to us that a sea-level rise of >≈10 m may be required to force retreat from the wedge deposited beneath Whillans Ice Stream over the past millennium or so (10).

Fig. 4.

Time-dependent response to instantaneous insertion of a wedge (Dt model). The final (Fnl.) geometries of the surface and bed are shown at the top, and the ice/air surface elevation (in meters) as a function of time and distance is color-coded. Ice thickness was held constant at the upglacier end (0 km, left). Most of the response to the insertion of the wedge was completed in ≈20 years.

Additional simulations of grounding-line/wedge interactions reveal a rich range of behavior (Fig. 4). Steady grounding-line positions on the upglacier side of a wedge were not simulated, consistent with the theory that such solutions are unstable (21). Migration of the grounding line across an otherwise horizontal bed is strongly influenced by any wedges encountered (22). The advance of the grounding line to a wedge crest is accelerated, but advance beyond the wedge crest slows relative to the no-wedge case. Similarly, the grounding line lingers on the downglacier side of a wedge during retreat and then races off of the upglacier edge. Under some circumstances, an advancing ice shelf has a more gradual basal slope than the upglacier side of a wedge, causing ice to bridge across water and ground on the wedge top, with the possibility of trapping water subglacially if the transverse geometry is favorable (23).

Our models show that the thickening of overlying ice in response to sedimentation stabilizes the grounding-line position against sea-level rise or buttressing reduction, and we infer that the extra thickness will tend to stabilize against any other environmental perturbation. Even a few centuries of sedimentation typically should provide stability against sudden sea-level rise from an outburst flood or surge of plausible-magnitude from another ice sheet, so sea-level change that is sufficiently large to control grounding-line positions will typically occur through noncatastrophic processes at usual rates (on the order of 1 cm/year or less) over several hundred to thousands of years. An ice-thickness change of ≈1.1 m has the same effect on grounding-line flotation as does a sea-level rise of 1 m, but this effect can be achieved more rapidly through sustained changes in snowfall or temperature. (Although slight warming causes very little surface melting in cold places, such as Whillans Ice Stream, a 1°C change can increase the surface melt rate by many tens of centimeters per year in warmer regions, such as some south-Greenland outlet glaciers.)

Changes in ice-shelf buttressing may be especially effective in overcoming the stabilizing effect of grounding-line sedimentation. As shown in Fig. 2, a 5% reduction in buttressing has about the same influence on the grounding line as does a 5-m rise in sea level for our reference case. A 5-m sea-level rise requires 500 years at the typical rate of the last deglaciation and 2500 years at 20th-century rates. Under some circumstances, a 5% decrease in buttressing will be achieved by an ≈5% decrease in the ice-shelf area generating side drag (24, 25). Because the basal melt rate of ice shelves is observed to increase with water temperature by ∼10 m year–1 °C–1 (26), even a fairly subtle perturbation in water temperature can cause a change of >5% in buttressing over years or decades. As an extreme example, a much larger change occurred in a few weeks during 2002 for the Larsen B Ice Shelf (27).

Grounding-zone sedimentation appears to be widespread to ubiquitous under major ice streams that control ice sheets (28, 29). Any near-stillstand of a grounding line (from an encounter with a preexisting basal high or narrowing of a trough or from time changes in basal lubrication, ice-shelf melting, or other environmental conditions) will be stabilized by sedimentation, except under special circumstances. In turn, because of spatial variability in sedimentation rates and physiographic variations between ice-stream troughs, any grounding-line retreats that are forced by sea-level rise are likely to be asynchronous for different ice sheets or different sectors of a large ice sheet.

Penck's hypothesis that Northern Hemisphere ice-sheet melting drove Antarctic ice-sheet retreat at the end of the ice age (4) may be accurate [also see (30)], with the delay in Antarctic retreat after the onset of sea-level rise being linked to sedimentary dynamics. Although our results indicate that sea-level changes of a few meters are unlikely to substantially affect ice-sheet behavior, 100-m changes should have considerable effects. However, sea level may exert the primary control on the ice sheet only if there is multimillennial stability in the other variables that affect ice sheets more quickly, such as water temperature under ice shelves. Because oceanic temperature probably changed with the deglaciation, perhaps linked to the return of North Atlantic Deep Water to the Southern Ocean (31), we consider it possible that ice-shelf changes contributed to or even dominated grounding-line retreat in some sectors of the ice sheet. Regardless, our results show that synchronous behavior of ice sheets with active sedimentary systems on millennial time scales is unlikely to indicate ice-sheet teleconnections via sea level and instead probably indicates common climatic forcing, which demonstrably can have very large and rapid effects on ice sheets.

Supporting Online Material

Materials and Methods

References and Notes

References and Notes

View Abstract

Stay Connected to Science

Navigate This Article