Ice sheet dynamics

Ice sheet dynamics describe the motion within large bodies of ice, such those currently on Greenland and Antarctica. Ice motion is dominated by the movement of glaciers, whose gravity-driven activity is controlled by two main variable factors: the temperature and strength of their bases. A number of processes alter these two factors, resulting in cyclic surges of activity interspersed with longer periods of inactivity, on both hourly and wict|centennial time scales.

Flow dynamics

Ice behaves as a brittle solid until its thickness exceeds about 50 meters (160 ft), after which the pressure permits plastic flow, where the glacier deforms under its own weight by creep. This creep occurs as a function of the weight of ice, the slope it is flowing down, the temperature, and the internal properties of the ice, but can be modelled quite accurately. Glaciers may also move by basal sliding, where the base of the glacier is lubricated by meltwater, allowing the glacier to slide over the terrain on which it sits. Meltwater may be provided by pressure-induced melting, friction or geothermal heat.

The top 50 meters of the glacier form the fracture zone, where ice moves as a single unit. Cracks form as the glacier moves over irregular terrain, which may penetrate the full depth of the fracture zone.

Glacial bottom processes

Most of the important processes controlling glacial motion occur in the ice-bed contact—even though it is only a few metres thick.doi ref|10.1146/annurev.earth.33.092203.122621] Glaciers will move by sliding when the basal shear stress drops below the shear resulting from the glacier's weight.clarifyme

:τ_D = ρgh sin α:where τ_D is the driving stress, and α the angle of repose.

:τ_B is the basal shear stress, a function of bed temperature and softness.

:τ_F, the shear stress, is the lower of τ_B and τ_D. It controls the rate of plastic flow, as per the figure (inset, right).

For a given glacier, the two variables are τ_D, which varies with h, the depth of the glacier, and τ_B, the basal shear stress.clarifyme

Basal shear stress

The basal shear stress is a function of three factors: the bed's temperature, roughness and softness.

Whether a bed is hard or soft depends on the porosity and pore pressure; higher porosity decreases the sediment strength (thus increases the shear stress τ_B). If the sediment strength falls far below τ_D, movement of the glacier will be accommodated by motion in the sediments, as opposed to sliding.Porosity may vary through a range of methods.
*Movement of the overlying glacier may cause the bed to undergo wict|dilatancy; the resulting shape change reorganises blocks. This reorganises closely packed blocks (a little like neatly folded, tightly packed clothes in a suitcase) into a messy jumble (just as clothes never fit back in when thrown in in a disordered fashion). This increases the porosity. Unless water is added, this will necessarily reduce the pore pressure (as the pore fluids have more space to occupy).
*Pressure may cause compaction and consolidation of underlying sediments. Since water is relatively incompressible, this is easier when the pore space is filled with vapour; any water must be removed to permit compression. In soils, this is an irreversible process.
*Sediment degradation by abrasion and fracture decreases the size of particles, which tends to decrease pore space, although the motion of the particles may disorder the sediment, with the opposite effect. These processes also generate heat, whose importance will be discussed later.

A soft bed, with high porosity and low pore fluid pressure, allows the glacier to move by sediment sliding: the base of the glacier may even remain frozen to the bed, where the underlying sediment slips underneath it like a tube of toothpaste. A hard bed cannot deform in this way; therefore the only way for hard-based glaciers to move is by basal sliding, where meltwater forms between the ice and the bed itself.Summarised from cite book|author=Boulton, Geoffrey S.|year=2006|chapter=Glaciers and their coupling with hydraulic and sedimentary processes|title=Glacier Science and Environmental Change|editpr=Peter G. Knight ]

Bed softness may vary in space or time, and changes dramatically from glacier to glacier. An important factor is the underlying geology; glacial speeds tend to differ more when they change bedrock than when the gradient changes.

As well as affecting the sediment stress, fluid pressure (p_w) can affect the friction between the glacier and the bed. High fluid pressure provides a buoyancy force upwards on the glacier, reducing the friction at its base. The fluid pressure is compared to the ice overburden pressure, p_i, given by ρgh. Under fast-flowing ice streams, these two pressures will be approximately equal, with an effective pressure (p_i - p_w) of 30 kPa; i.e. all of the weight of the ice is supported by the underlying water, and the glacier is afloat.

A number of factors can affect bed temperature, which is intimately associated with basal meltwater. The melting point of water decreases under pressure, meaning that water melts at a lower temperature under thicker glaciers. This acts as a "double whammy", because thicker glaciers have a lower heat conductance, meaning that the basal temperature is also likely to be higher.

Bed temperature tends to vary in a cyclic fashion. A cool bed has a high strength, reducing the speed of the glacier. This increases the rate of accumulation, since newly fallen snow is not transported away. Consequently, the glacier thickens, with three consequences: firstly, the bed is better insulated, allowing greater retention of geothermal heat. Secondly, the increased pressure can facilitate melting. Most importantly, τ_D is increased. These factors will combine to accelerate the glacier. As friction increases with the square of velocity, faster motion will greatly increase frictional heating, with ensuing melting - which causes a positive feedback, increasing ice speed to a faster flow rate still: west Antarctic glaciers are known to reach velocities of up to a kilometre per year.Eventually, the ice will be surging fast enough that it begins to thin, as accumulation cannot keep up with the transport. This thinning will increase the conductive heat loss, slowing the glacier and causing freezing. This freezing will slow the glacier further, often until it is stationary, whence the cycle can begin again.

Supraglacial lakes represent another possible supply of liquid water to the base of glaciers, so can play an important role in accelerating glacial motion.Lakes of a diameter greater than ~300 m are capable of creating a fluid-filled crevasse to the glacier/bed interface.When these crevasses form, the entirety of the lake's (relatively warm) contents can reach the base of the glacier in as little as 2-18 hours - lubricating the bed and causing the glacier to surge.cite conference
author = Krawczynski, M.J.
coauthors = Behn, M.D.; Das, S.B.; Joughin, I.
year = 2007
title = Constraints on melt-water flux through the West Greenland ice-sheet: modeling of hydro-fracture drainage of supraglacial lakes
url = http://www.agu.org/cgi-bin/wais?jj=C41B-0474
accessdate = 2008-03-04
booktitle=Eos Trans. AGU,
volume=88(52)
pages= Fall Meet. Suppl., Abstract C41B-0474]

Finally, bed roughness can act to slow glacial motion. The roughness of the bed is a measure of how many boulders and obstacles protrude into the overlying ice. Ice flows around these obstacles by melting under the high pressure on their lee sides; the resultant meltwater is then forced down a steep pressure gradient into the cavity arising in their stoss, where it re-freezes. Cavitation on the stoss side increases this pressure gradient, which assists flow.

Erosional effectsanchor|erosion

Because ice can flow faster where it is thicker, the rate of glacier-induced erosion is directly proportional to the thickness of overlying ice. Consequently pre-glacial low hollows will be deepened and pre-existing topography will be amplified by glacial action, while wict|nunataks, which protrude above ice sheets, barely erode at all - erosion has been estimated as 5 m per 1.2 million years!cite journal | author = | year = | doi = 10.1038/ngeo201| title = Fjord insertion into continental margins driven by topographic steering of ice| journal = Nature Geoscience
*Non-techincal summary: cite journal | author = Kleman | year = 2008 | doi = 10.1038/ngeo210 | title = Geomorphology: Where glaciers cut deep | journal = Nature Geoscience ] This explains, for example, the deep profile of fjords, which can reach a kilometer in depth as ice is topographically steered into them. Being the principal conduits for draining ice sheets, fjords' extension inland increases the rate of ice sheet thinning. It also makes the ice sheets more sensitive to changes in climate and the ocean.

Pipe and sheet flow

The flow of water under the glacial surface can have a large effect on the motion of the glacier itself. Subglacial lakes contain significant amounts of water, which can move fast: cubic kilometres can be transported between lakes over the course of a couple of years.cite doi|10.1126/science.1136897]

This motion is thought to occur in two main modes: pipe flow involves liquid water moving through pipe-like conduits, like a sub-glacial river; sheet flow involves motion of water in a thin layer. A switch between the two flow conditions may be associated with surging behaviour. Indeed, the loss of sub-glacial water supply has been linked with the shut-down of ice movement in the Kamb ice stream. The subglacial motion of water is expressed in the surface topography of ice sheets, which slump down into vacated subglacial lakes.

Boundary conditions

The interface between an ice stream and the ocean is a significant control of the rate of flow. Ice shelves - thick layers of ice floating on the sea - can stabilise the glaciers that feed them. These tend to have accumulation on their tops, may experience melting on their bases, and calve icebergs at their periphery. The catastrophic collapse of the Larsen B ice shelf in the space of three weeks during February 2002 yielded some unexpected observations. The glaciers that had fed the ice sheet (Crane, Jorum, Green, Hektoria - see image) increased substantially in velocity. This cannot have been due to seasonal variability, as glaciers flowing into the remnants of the ice shelf (Flask, Leppard) did not accelerate.cite doi|10.1029/2004GL020670]

Ice shelves exert a dominant control in Antarctica, but are less important in Greenland, where the ice sheet meets the sea in fjords. Here, melting is the dominant ice removal process, resulting in predominant mass loss occurring towards the edges of the ice sheet, where icebergs are calved in the fjords and surface meltwater runs into the ocean.

Tidal effects are also important; the influence of a 1 m tidal oscillation can be felt as much as 100 km from the sea. On an hour-to-hour basis, surges of ice motion can be modulated by tidal activity. During larger spring tides, an ice stream will remain almost stationary for hours at a time, before a surge of around a foot in under an hour, just after the peak high tide; a stationary period then takes hold until another surge towards the middle or end of the falling tide.doi ref|10.1126/science.1087231] cite doi|10.1029/2002GL016329] At neap tides, this interaction is less pronounced, without tides surges would occur more randomly, approximately every 12 hours.

Ice shelves are also sensitive to basal melting. In Antarctica, this is driven by heat fed to the shelf by the circumpolar deep water current, which is 3 °C above the ice's melting point.cite doi|10.1029/2006GL028154]

As well as heat, the sea can also exchange salt with the oceans. The effect of latent heat, resulting from melting of ice or freezing of sea water, also has a role to play. The effects of these, and variability in snowfall and base sea level combined, account for around 80 mm a^-1 variability in ice shelf thickness. In contrast, the ice shelf in Pine Island Bay is thinning at 5.5 m a^-1. Clearly the influence of tides, combined with processes affecting the strength of the glacier's bed, are far more significant than these oceanic effects.

Long term changes

Over long time scales, certain patterns may emerge in the flux of ice sheets. For instance, Heinrich events are thought to represent the disintegration of the Laurentide ice sheet, and the resulting influx of armardas of icebergs into the North Atlantic. These events appear to have a 7,000 year peridoicity; if this is the case, it is likely that internal ice sheet "binge-purge" cycles are responsible for the observed effects.

Effects of climate change on ice sheet dynamics

The implications of the current climate change on ice sheets are difficult to constrain. It is clear that increasing temperatures are resulting in reduced ice volumes globally.Sections 4.5 and 4.6 of IPCC4/wg1/4] (Due to increased precipitation, the mass of parts of the Antarctic ice sheet may currently be increasing, but the total mass balance is unclear.)

Since the surging nature of ice sheet motion is a relatively recent discovery, and is still a long way from being entirely understood, no models have yet made a comprehensive evaluation of the effects of climate change. However, it is clear that climate change will act to destabilise ice sheets by a number of mechanisms.

Rising sea levels will reduce the stability of ice shelves, which have a key role in reducing glacial motion. Some Antarctic ice shelves are currently thinning by tens of metres per year, and the collapse of the Larsen B shelf was preceded by thinning of just 1 metre per year. Further, increased ocean temperatures of 1 °C may lead to up to 10 m a^-1 of basal melting. Ice shelves are always stable under mean annual temperatures of −9 °C, but never stable above −5 °C; this places regional warming of 1.5 °C, as preceded the collapse of Larsen B, in context.

Increasing global temperatures take around 10,000 years to directly propagate through the ice before they influence bed temperatures, but may have an effect through increased surfacal melting, producing more supraglacial lakes, which may feed warm water to glacial bases and facilitate glacial motion.

Also, in areas of increased precipitation, such as Antarctica, the addition of mass will increase rate of glacial motion, hence the turnover in the ice sheet.

Observations, while currently limited in scope, do agree with these predictions of an increasing rate of ice loss from both Greenland and Antarctica.

A possible positive feedback may result from shrinking ice caps, in volcanically active Iceland at least. Isostatic rebound may lead to increased volcanic activity, causing basal warming - and, through co2 release, further climate change. [cite doi|10.1029/2008GL033510]

References