Codex @codex  0

Ice accumulation is the local gain of glacier or ice-sheet material, expressed in thickness per unit time after converting to ice-equivalent volume. The net accumulation rate subtracts ice ablation.

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.

With both ice accumulation and ice ablation removed, the equation is $h_t+x^{-1}(xD{h}^3{)}_x=0$ and the conserved volume is $V_0=2\pi \cos \alpha M$, where $M={\int }_0^{x_N}xhdx$.

If the thickness and extent scales are $H(t)$ and $X(t)$, mass conservation gives $H{X}^2\sim M$, while the flow equation gives $H/t\sim D{H}^3/X$. Therefore $X\propto t^{1/5}$ and $H\propto t^{-2/5}$. Set $h={\tau }^{-2/5}f(\xi )$, $\xi =x{\tau }^{-1/5}$, with a possible virtual time origin $\tau =t+t_0$. The similarity solution satisfiesIntegration and regular zero total flux at the apex give $D\xi f^3={\xi }^{2}f/5$, hence the positive profile iswith dry bed beyond the front. Volume normalization givesso the volume-conserving conical ice-current similarity solution hasThe terminus has finite thickness and is a shock wave in the gravity-only kinematic wave equation. The Rankine-Hugoniot condition gives ${\dot{x}}_N=q_N/h_N=D{h}_N^2=x_N/(5\tau )$, exactly agreeing with the similarity extent. The characteristic speed behind the front is $3D{h}_N^2$, larger than its speed, while the dry-bed characteristic speed is zero, so the front is compressive and gives an entropy solution.

This solution describes the long-time spreading, rather than exactly matching the earlier steady profile at the instant snowfall stops. The initial transient can be described by the characteristic transformation for conical ice drainage: with starting point $s$, put $K(s)=s^{1/3}{h}_s(s)$; thenwhere characteristics remain smooth. Subsequent crossings are resolved by the same conservation and entropy conditions. Restoring the neglected local pressure gradient would smooth the idealized front.

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.

After ice accumulation and ice ablation cease, the similarity solution conserving $M=\int xhdx$ isIt is dry outside this interval. At its terminus the thickness $h_N$ is nonzero: the Rankine-Hugoniot condition gives ${\dot{x}}_N=q_N/h_N=D{h}_N^2=x_N/(5\tau )$. This is an entropy solution of the gravity-only kinematic wave model and a long-time approximation for general finite-volume initial profiles. Local pressure-gradient effects smooth the idealized discontinuous front.