Modeling melange back-pressure

Equation (33) above, describing the stress boundary condition for ice shelves, can be written in terms of velocity components:

(38)2νH(2ux+uy)nx+2νH(uy+vx)ny=bh(pice(z)pwater(z))dznx,2νH(uy+vx)nx+2νH(2vy+ux)ny=bh(pice(z)pwater(z))dzny.

Here ν is the vertically-averaged ice viscosity, H is the ice thickness, b is the elevation of the bottom and h of the top ice surface, pwater and pice are pressures of the column of ice and water, respectively:

(39)pice=ρiceg(hz),pwater=ρwatergmax(zsea levelz,0).

We call the integral on the right hand side of (38) the pressure difference term.

It can be re-written as

(40)bhpice(z)pwater(z)dz=H(p¯icep¯water),wherep¯ice=12ρicegH,p¯water=12ρwatergmax(zsb,0)2H.

PISM’s ocean model components provide p¯water, the vertically-averaged pressure of the water column adjacent to an ice margin.

To model the effect of melange [105] on the stress boundary condition we modify the pressure difference term in (38), adding p¯melange, the vertically-averaged melange back pressure:

(41)bh(pice(pwater+p¯melange))dz.

By default, p¯melange is zero, but PISM implements two ocean model components to support scalar time-dependent melange pressure forcing. Please see the Climate Forcing Manual for details.


Previous Up Next