Solution (source code)

= Solution

In the shallow, equal-density approximation, melting supplies water at rate $v_m$ per unit wetted area. Combining this source with the pressure-driven <volume flux> gives the <melting-source gravity-current equation>,
$$
\boxed{h_t+q_x=v_m,\qquad q=-Dh^2h_x,\qquad
h_t=D(h^2h_x)_x+v_m}
$$
inside the wet region, with no melting source ahead of the front. In the PDF's signed convention this is $h_t+\gamma(h^2h_x)_x=v_m$. If ice and water densities are distinguished, the meltwater-volume source is $s_m=(\rho_i/\rho_w)v_m$; the same equations then use $s_m$. Geometrical excavation of ice has rate $v_m$, so a density difference requires the confining ice/water-volume mechanics to be treated consistently. The equal-density shallow model identifies these two volumes.

Let $x_-(t)<x_+(t)$ bound the current, and $\mathcal A(t)=\int_{x_-}^{x_+}h\,dx$ be its volume per span. The moving-boundary integral of continuity is
$$
\boxed{\dot{\mathcal A}=s_m(x_+-x_-)-q(x_+)+q(x_-)
+h(x_+)\dot x_+-h(x_-)\dot x_-.}
$$
This exhibits both the melting contribution and any volume swept out by an advancing boundary. For a dry zero-height nose at each end, with no flux through those ends, the boundary terms vanish and
$$
\boxed{\dot{\mathcal A}=s_m\mathcal L(t),\qquad
\ddot{\mathcal A}=s_m\dot{\mathcal L}(t),\qquad \mathcal L=x_+-x_-,}
$$
where the second relation uses constant $s_m$. For a symmetric pulse $\mathcal L=2L$, these become $\dot{\mathcal A}=2s_mL$ and $\ddot{\mathcal A}=2s_m\dot L$; on a reflecting half-line they are $s_mL$ and $s_m\dot L$. Multiplication by $B_y$ gives the three-dimensional water volume.

The <meltwater production over an advancing footprint> can equivalently be expressed as $\mathcal A_m(t)=s_m\int_0^t\mathcal L(t')\,dt'$. At an individual position melting begins only at its wetting time, so new area contributes no finite melt thickness at the instant of arrival. Advancement enlarges the area subsequently melting and hence increases the volume-production rate. The fixed-volume similarity from part (b) cannot simply be reused after adding this source.