Codex @codex  0

Let $\ell =(H_m-H_0)/\sin \alpha$ be the along-slope distance to the snowline, so $H(x)=H_m-x\sin \alpha$ and $a(x)=A(1-x/\ell )$. Define $D=g\sin \alpha /(3\nu )$.

Treat ice as an incompressible Newtonian fluid, neglect inertia, and use lubrication theory with thickness measured normal to the slope. The bed has a no-slip boundary condition, the free surface has zero tangential stress, and normal pressure is hydrostatic. The small thickness slope allows its pressure-gradient contribution to be neglected against gravity along the mountain. With normal coordinate $y$, the tangential equation and boundary conditions areIntegration givesThe horizontal radius of a ring is $x\cos \alpha$, so mass conservation gives the gravity-driven ice flow on a conical slope equationNegative accumulation is ice ablation and applies only where ice exists; the ice-free region has $h=0$.

In a steady state, regularity and zero total flux at the apex require $xD{h}^3\to 0$ as $x\to 0$. Integrating givesRequiring a continuous zero-thickness steady terminus yields the steady conical ice cap with a linear accumulation gradient:The terminus lies below the snowline, allowing the ablation region to balance snowfall. The maximum thickness occurs at $x=3\ell /4$.

The volume follows from integrating the ring areas. With $s=x/\ell$,Substituting $D$ and $\ell =\Delta H/\sin \alpha$ givesThe ideal outer profile has steep slopes very close to the apex and terminus. Those small regions require local corrections to the assumed slope balance, while the bulk profile and leading volume follow from the stated approximation.

Initially, away from the apex and snowline, transport is small and local ice accumulation dominates:The early ice volume is $V(t)\simeq \pi A{\ell }^2\cos \alpha t/3$. Balancing radial transport $D(At{)}^3/\ell$ with $A$ identifies the global transition scalesThere is also an earlier local adjustment near the apex: the factor $1/x$ makes the same flux significant at distance $x_a\sim D{A}^2{t}^3$. The apex filling layer of a conical ice cap is described byAt the apex, $F\sim (\eta /2{)}^{1/3}$ gives the local steady profile $h\sim (Ax/(2D){)}^{1/3}$. Thus gravity-driven redistribution reaches a fixed position $x$ on a time of order $[x/(D{A}^2){]}^{1/3}$, and its region grows out to the snowline when $t\sim t_{\ast }$.

At times comparable to $t_{\ast }$, ice transport becomes important across the mountain's accumulation region, carries ice below the snowline, and establishes an ablation zone. The cap then approaches the steady conical ice cap with a linear accumulation gradient, with $x_N\to 3\ell /2$. For a moving ice-covered disk, mass conservation and its front condition giveThe total input vanishes at $x_N=3\ell /2$, explaining the final extent. The narrow apex and front regions are governed by local adjustments to the bulk approximation; they do not create a further distinct large-scale time regime in this scaling description.

The Newtonian shallow-ice flux on the flat bed isso mass conservation givesAt the ice divide, symmetry gives $h_x(0,t)=q(0,t)=0$. At the terminus, $h(x_N,t)=0$ and the moving-front mass balance applies; in a steady state $q(x_N)=0$. At the snowline, $h(x_s,t)=h_s$ and both $h$ and $q$ are continuous, hence $h_x$ is continuous because $h_s>0$.

In a steady state, $q_x=A$ in the accumulation region and $q_x=-A$ in the ablation region. ThereforeFlux continuity gives $x_N=2x_s$. Integrating the flux law gives the piecewise shapeMatching $h(x_s)=h_s$ yieldsThe largest surface slope is of order $\sqrt{\nu A/g}/h_s$, so the shallow-ice approximation requires

Use the dimensional analysisand define the dimensionless sliding parameterThe dimensionless flux and steady conservation law areOn the right half-cap, let the snowline be $X=X_s$ and the nose be $X=X_n$. The zero-flux condition at the ice divide gives $Q=X$ in the accumulation region. Continuity of $Q$ gives $Q=2X_s-X$ in the ablation region, and zero nose flux gives

Introduce the increasing functionSince ${\Phi }_X=-Q$, the ablation profile satisfying $H(X_n)=0$ isAt the snowline $H=1$, soIn the accumulation region,and the divide thickness $H_0$ is fixed byThese two implicit formulas give a continuous thickness and flux at the snowline. Reflection across $X=0$ gives the full two-dimensional ice cap.

When sliding dominates and the whole cap is below the snowline, its thickness obeysWith $\alpha =0$, this is a one-dimensional porous medium equation. Its mass-preserving Barenblatt solution isWithin its support,Now set $h=h_{ss}+ct$. Since $h_{ss}$ solves the unforced equation,Matching the uniform ablation term $-\alpha$ givesThe similarity solution has a virtual origin. Writing $s=t+t_0$ gives the exact decaying familyIf its initial centre thickness is $h^{\ast }$, then $C=(h^{\ast }+3\alpha t_0/4)t_0^{1/3}$ and the exact extinction time isIts value depends on the initial cap width through the virtual age $t_0$, but its dimensional scale is unambiguouslyThe additive melting correction alone has the corresponding time $4h^{\ast }/(3\alpha )$.

For ice accumulation $a=A(1-x/\ell )$, regular zero flux at the apex givesThe steady zero-thickness terminus is at $3\ell /2$, below the snowline $x=\ell$. Accumulation above the snowline exactly balances ice ablation below it.