Newtonian grounding-line flux Created 2026-09-28 Updated 2026-11-05
A Newtonian fluid model of an ice sheet does not determine its grounding line volume flux per unit width from thickness alone without a basal law and matching conditions. Let be flotation thickness, the ice kinematic viscosity, gravity and the reduced gravity. A local cubic scaling can be illustrated by patching two explicit approximations on a horizontal bed.
For no-slip Newtonian shallow-ice flux, . In a locally steady region with negligible accumulation, mass conservation gives and . Imposing the unbuttressed Newtonian grounding-line stress condition at the same matching section then yields
This is the result of the stated elementary patch, not a resolved grounding-zone solution. It shows why a bare numerical coefficient cannot be inferred from Newtonian rheology alone.
Basal sliding, longitudinal stresses and contact constraints require a transition calculation; they can change the flux-thickness law as well as its coefficient. A systematic Stokes flow treatment matches the outer shearing/sliding sheet to an extensional boundary layer and supplies effective flotation and flux conditions.[ref] The patch above should therefore be used as a transparent closure example, rather than as a universal predictive grounding-line law.
The Newtonian shallow-ice flux on the flat bed is
so mass conservation gives
At the ice divide, symmetry gives . At the terminus, and the moving-front mass balance applies; in a steady state . At the snowline, and both and are continuous, hence is continuous because .
In a steady state, in the accumulation region and in the ablation region. Therefore
Flux continuity gives . Integrating the flux law gives the piecewise shape
Matching yields
The largest surface slope is of order , so the shallow-ice approximation requires