Let be the pore fraction and put . In the long-wave approximation, the pore pressure is hydrostatic, , so Darcy's law gives the horizontal volume flux per unit width
Local mass conservation therefore gives the Boussinesq equation for an unconfined aquifer
The river and drainage-divide conditions are
The long-wave approximation makes independent of at leading order. Mass conservation and symmetry about then give
For an incompressible Newtonian fluid,
The leading normal-stress balance on either surface is
Since , it follows that
For a slice of length , the axial forces are the integrated normal stresses on its vertical ends, ambient pressure on the varying end height, and the horizontal components of surface tension on its two sloping faces. Expanding their difference to first order in cancels the uniform terms and the lower-order capillary terms, leaving
The kinematic condition on is . Substituting gives the second required relation,
Let measure height above the sponge. Under the long-wave approximation, hydrostatic pressure gives , and the horizontal balance for viscous fluid flow is
The no-slip boundary condition and stress-free boundary condition give
Integrating this lubrication theory profile gives the horizontal volume flux per unit span
Local mass conservation includes the downward loss found in part i:
Consequently the required nonlinear partial differential equation is
If the imposed inlet flux is and the moving front is , sufficient boundary conditions are
For a current released onto a dry substrate one also takes the initial condition away from the source. The front position is part of this moving-boundary problem.
Let be the local water-film thickness. To first order in the interface amplitudes,
The lubrication theory flux down the vertical surface, with a no-slip boundary condition at the ice and a stress-free boundary condition at the water-air interface, is
For the unperturbed film , so
The linearized Young–Laplace equation gives the capillary-pressure perturbation
Because the prescribed volume flux is uniform, its first-order perturbation must vanish. Linearization of the flux law gives
Thus, with ,
and hence
The complex amplitude ratio records the phase shift caused by surface tension in the long-wave approximation.