 Solution

The water mass conservation equation inherited from the earlier model is $ϵh_t+(hu{)}_x=0$. It is essential: without it the printed displayed equations are not closed. Take the slow bed branch $ϵ=0$, and perturb the uniform state by $h=1+\hat{h}{e}^{\lambda t+ikx}$, $u=1+U{e}^{\lambda t+ikx}$, $c=1+C{e}^{\lambda t+ikx}$, $s=S{e}^{\lambda t+ikx}$ and $\eta =1+\hat{\eta }{e}^{\lambda t+ikx}$. For $k\ne 0$, mass conservation gives $\hat{h}=-U$ and $\hat{\eta }=S-U$.

Put $m=E^{\ast \prime }(1)+1$; the extra $1$ comes from linearizing $c/u$. The linearized momentum and suspended-sediment equations areEliminating $U,C$ gives the dispersion relationAway from $F=1$, for fixed nonzero $k$ and negligible $\delta$,Thus growth requires $m/(F^2-1)>0$. For $F>1$, any $E^{\ast \prime }(1)>-1$ gives an unstable branch; for $F<1$ growth would require $E^{\ast \prime }(1)<-1$. A physically increasing erosion function therefore gives the supercritical instability in this simplified model, rather than proving growth of ordinary subcritical dunes. Retaining drag giveswhich shows its stabilizing correction and the need to reconsider the limit near critical flow.

The negative bed-to-surface ratio for subcritical flow means bed crests lie beneath surface troughs: the dune phase relation. Its positive sign for supercritical flow gives aligned crests: an antidune. In the unstable supercritical case the imaginary part gives upstream migration, since $-\mathrm{Im}\lambda /k=-m/[(F^2-1)(1+k^2)]$. The figure shows these two phase relations. They identify the bedforms but do not by themselves establish growth. The qualitative limitation of the suspended-sediment mechanism is also discussed in the official ^[ref].