Ice sheet dynamics describe the motion within large bodies of ice

Ice sheet

An ice sheet is a mass of glacier ice that covers surrounding terrain and is greater than 50,000 km² , thus also known as continental glacier...

, such those currently on Greenland

Greenland

Greenland is an autonomous country within the Kingdom of Denmark, located between the Arctic and Atlantic Oceans, east of the Canadian Arctic Archipelago. Though physiographically a part of the continent of North America, Greenland has been politically and culturally associated with Europe for...

 and Antarctica. Ice motion is dominated by the movement of glacier

Glacier

A glacier is a large persistent body of ice that forms where the accumulation of snow exceeds its ablation over many years, often centuries. At least 0.1 km² in area and 50 m thick, but often much larger, a glacier slowly deforms and flows due to stresses induced by its weight...

s, 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 time scales.

Flow dynamics

The main cause of flow within glaciers can be attributed to an increase in the surface slope, brought upon by an imbalance between the amounts of accumulation vs. ablation. This imbalance increases the shear stress

Shear stress

A shear stress, denoted \tau\, , is defined as the component of stress coplanar with a material cross section. Shear stress arises from the force vector component parallel to the cross section...

 on a glacier until it begins to flow. The flow velocity and deformation will increase as the equilibrium line between these two processes is approached, but are also affected by the slope of the ice, the ice thickness and temperature.

When the amount of strain (deformation) is proportional to the stress being applied, ice will act as an elastic solid. Ice will not flow until it has reached a thickness of 30 meters (98 ft), but after 50 meters (164 ft), small amounts of stress can result in a large amount of strain, causing the deformation to become a plastic flow

Plasticity (physics)

In physics and materials science, plasticity describes the deformation of a material undergoing non-reversible changes of shape in response to applied forces. For example, a solid piece of metal being bent or pounded into a new shape displays plasticity as permanent changes occur within the...

 rather than elastic. At this point the glacier will begin to deform under its own weight and flow across the landscape. According to the Glen–Nye Flow law, the relationship between stress and strain, and thus the rate of internal flow, can be modeled as follows:


where: = shear strain (flow) rate = stress = a constant between 2–4 (typically 3 for most glaciers) that increases with lower temperature = a temperature-dependent constant

The lowest velocities are near the base of the glacier and along valley sides where friction acts against flow, causing the most deformation. Velocity increases inward toward the center line and upward, as the amount of deformation decreases. The highest flow velocities are found at the surface, representing the sum of the velocities of all the layers below.

Glaciers may also move by basal sliding

Basal sliding

Basal sliding is the act of a glacier sliding over the bed before it due to meltwater under the ice acting as a lubricant. This movement very much depends on the temperature of the area, the slope of the glacier, the bed's sediment size, the amount of meltwater from the glacier, and the glacier's...

, 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 produced by pressure-induced melting, friction or geothermal heat. The more variable the amount of melting at surface of the glacier, the faster the ice will flow.

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 meters thick. Glaciers will move by sliding when the basal shear stress drops below the shear resulting from the glacier's weight.

τ_D = ρgh sin α

where τ_D is the driving stress, and α the ice surface slope in radians.

τ_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.

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 ; 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.

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 lake

Supraglacial lake

A supraglacial lake is any pond of liquid water on the top of a glacier. Although these pools are , they may reach kilometers in diameter and be several meters deep...

s represent another possible supply of liquid water to the base of glaciers, so they 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

Surge (glacier)

Glacial surges are short-lived events where a glacier can advance substantially, moving at velocities up to 100 times faster than normal. Surging glaciers are clustered around a few areas. High concentrations of surging glaciers can be found in Svalbard, Canadian Arctic islands, Alaska and Iceland...

. Water that reaches the bed of a glacier may freeze there, increasing the thickness of the glacier by pushing it up from below.

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 effects

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 s, which protrude above ice sheets, barely erode at all - erosion has been estimated as 5 m per 1.2 million years! This explains, for example, the deep profile of fjord

Fjord

Geologically, a fjord is a long, narrow inlet with steep sides or cliffs, created in a valley carved by glacial activity.-Formation:A fjord is formed when a glacier cuts a U-shaped valley by abrasion of the surrounding bedrock. Glacial melting is accompanied by rebound of Earth's crust as the ice...

s, 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.

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

Ice calving

Ice calving, also known as glacier calving or iceberg calving, is a form of ice ablation or ice disruption. It is the sudden release and breaking away of a mass of ice from a glacier, iceberg, ice front, ice shelf, or crevasse...

 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.

Ice shelves exert a dominant control in Antarctica, but are less important in Greenland, where the ice sheet meets the sea in fjord

Fjord

Geologically, a fjord is a long, narrow inlet with steep sides or cliffs, created in a valley carved by glacial activity.-Formation:A fjord is formed when a glacier cuts a U-shaped valley by abrasion of the surrounding bedrock. Glacial melting is accompanied by rebound of Earth's crust as the ice...

s. 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. 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.

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⁻¹ variability in ice shelf thickness.

Long term changes

Over long time scales, ice sheet mass balance is governed by the amount of sunlight reaching the earth. This variation in sunlight reaching the earth, or insolation

Insolation

Insolation is a measure of solar radiation energy received on a given surface area in a given time. It is commonly expressed as average irradiance in watts per square meter or kilowatt-hours per square meter per day...

, over geologic time is in turn determined by the angle of the earth to the sun and shape of the Earth's orbit, as it is pulled upon by neighboring planets; these variations occur in predictable patterns called Milankovitch cycles

Milankovitch cycles

Milankovitch theory describes the collective effects of changes in the Earth's movements upon its climate, named after Serbian civil engineer and mathematician Milutin Milanković, who worked on it during First World War internment...

. Milankovitch Cycles dominate climate on the Glacial/Interglacial timescale, but there exist variations in ice sheet extent that are not linked directly with insolation.

For instance, during at least the last 100,000 years, portions of the ice sheet covering much of North America, the Laurentide ice sheet

Laurentide ice sheet

The Laurentide Ice Sheet was a massive sheet of ice that covered hundreds of thousands of square miles, including most of Canada and a large portion of the northern United States, multiple times during Quaternary glacial epochs. It last covered most of northern North America between c. 95,000 and...

 broke apart sending large flotillas of icebergs into the North Atlantic. When these icebergs melted they dropped the boulders and other continental rocks they carried, leaving layers known as ice rafted debris. These so-called Heinrich events, named after their discoverer Hartmut Heinrich

Hartmut Heinrich

Hartmut Heinrich is a German marine geologist and climatologist. Dr. Heinrich is Head of the Physics Department at the Federal Maritime and Hydrographic Agency in Hamburg. He is actively involved in GOOS...

, appear to have a 7,000-10,000 year peridoicity, and occur during cold periods within the last interglacial.

Internal ice sheet "binge-purge" cycles may be responsible for the observed effects, where the ice builds to unstable levels, then a portion of the ice sheet collapses. External factors might also play a role in forcing ice sheets. Dansgaard-Oeschger events are abrupt warmings of the northern hemisphere occurring over the space of perhaps 40 years. While these D-O events occur directly after each Heinrich event, they also occur more frequently - around every 1500 years; from this evidence, paleoclimatologists surmise that the same forcings may drive both Heinrich and D-O events.

Hemispheric Asynchrony in Ice Sheet Behavior has been observed by linking short term spikes of methane in Greenland ice cores and Antarctic ice cores. During Dansgaard-Oeschger events, the northern hemisphere warmed considerably, dramatically increasing the release of methane from wetlands, that were otherwise tundra during glacial times. This methane quickly distributes evenly across the globe, becoming incorporated in Antarctic and Greenland ice. With this tie, paleoclimatologists have been able to say that the ice sheets on Greenland only began to warm after the Antarctic ice sheet had been warming for several thousand years. Why this pattern occurs is still open for debate.

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. (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 metres per year 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 release, further climate change.