Friction-dominated grounding-line evolution law ($[\alpha {\rho }_w/\rho -H_x]{\dot{x}}_G$)

For a bed depth $b=\alpha x$, a thin subglacial till layer of thickness $l$ is idealized by Newtonian till lubrication, with dynamic viscosity $\lambda \mu$ and basal drag $\lambda \mu u/l$. Neglecting longitudinal stress divergence in the bulk yields $u=-(l/\lambda )(g/\nu )H{h}_x$. Retaining the unbuttressed Newtonian grounding-line stress condition as boundary data and differentiating ice-sheet flotation along the moving grounding line givesIndeed, mass conservation gives $H_t=-u{H}_x-H{u}_x$, $H_x=h_x+\alpha$, and ice-sheet flotation gives $H_t=(\alpha {\rho }_w/\rho -H_x){\dot{x}}_G$. This is a formal local friction closure, not a uniformly valid removal of membrane stress near every grounding line. If the coefficient of ${\dot{x}}_G$ vanishes, it is an implicit compatibility condition and cannot be divided to obtain a finite speed.

With constant volume flux per unit width $Hu=q_0>0$, define $\beta =l/\lambda$, $\mathcal{H}=H/\beta$, and $\mathcal{Q}=\nu q_0/(g{\beta }^3)$. The friction closure gives $H_x=\alpha -q_0\nu /(\beta g{H}^2)$. Differentiating constant flux and applying the unbuttressed Newtonian grounding-line stress condition gives $q_0{H}_x/H+g^{\prime }{H}^2/(8\nu )=0$. Substitution proves the displayed relation. For $\alpha \ge 0$, $g^{\prime }>0$ and $\mathcal{Q}>0$, the polynomial $(g^{\prime }/8g){\mathcal{H}}^5+\alpha \mathcal{Q}{\mathcal{H}}^2-{\mathcal{Q}}^2$ increases strictly on the positive half-line, begins negative and tends to infinity; hence it has exactly one positive root. Its root satisfies $\mathcal{Q}>\alpha {\mathcal{H}}^2$, so the thickness decreases seaward at that steady grounding line.