Driving stress

Contents

Driving stress#

Ice-sheet flow is driven by gravity. Specifically it is a quantity called the driving stress that drives flow:

τd=ρgHHx,

where ρ is the density of ice ( 9.17 kg m3), g is gravitational acceleration ( ~9.8 m s2), H is the ice thickness, x is the horizontal coordinate oriented in teh ice flow direction.

Derivation.#

Consider a column of ice that extends from the ice-sheet bed to the surface and has length δx in the x direction and length δy in the y directon perpendicula to flow. We will assume that the bed is flat. The pressure at each elevation z is given by

P=ρg(Hz),

where P is the pressure, and z is the vertical coordinate (increasing upward). The equation above comes from integrating the hydrostatic pressure equation

dPdx=ρgz.

(todo: replace this figure, add figure caption)

image

To calculate the force acting on one side of the ice column (blue arrow), with thickness H, we can consider that force = pressure × area. However, ppressure p depends on z, so we need to integrate over the thickness of the ice, rather than just multiply by the area. We will call the force on the left-hand side of the column F1:

F1=δy0HP(z)dz
F1=δy0Hρg(Hz)dz=δyρg[Hzz22]=δyρg[H2H22]
F1=12δyρgH2

Note that in this calculation, we have assumed constant ρ and g, and that air pressure is negligible compared to ice pressure.

The more relevant quantity to consider here is not just the force acting on one side of the ice column, but rather the difference between the forces on either side. If we imagine that the ice column is part of an ice sheet with a surface slope of dHdx, and has thickness x in the x-direction, then the force acting on the opposite side of the ice column can be calculated as follows:

F2=12δyρgH2+δxddx(12δyρgH2)
ΔF=F1F2=12δyρgH2(12δyρgH2+δxddx(12δyρgH2))
ΔF=12δxδyρgH2x

We can now invoke the chain rule to rewrite this as

ΔF=12δxδyρg(2HHx)=δxδyρgHHx

This is the stres imbalance integrated over the is ice thickness and across a width δy perpendicular to flow. We define the driving stress as the force imbalance per unit area of the bed:

τd=ΔFA=ΔFδxδy
τd=ρgHHx

This makes sense intuitively; there is more driving stress where the ice is thicker (higher H) and steeper (higher Hx).

We can generalize this to two dimensions as follows:

τd=ρgH(Hxi^+Hyj^)=ρgHH

The driving stress field will point in the direction in which the slope changes the fastest.

We can also generalize this to account for cases when the be topography is not flat:

τd=ρgHhx

where h is the ice surface height.