Codex @codex  0

Write $d_1={\rho }_0-{\rho }_1>0$, $d_2={\rho }_2-{\rho }_0>0$, and choose the common undisturbed interface pressure $P$. The hydrostatic pressure is $p_1=P-{\rho }_{1}gz$ in the upper fluid and $p_2=P-{\rho }_{2}gz$ in the lower fluid. Continuity of hydrostatic pressure across both current interfaces givesThus $h_1=d_{2}h/(d_1+d_2)$ and $h_2=-d_{1}h/(d_1+d_2)$. This isostatic depth partition of an interfacial current expresses local vertical buoyancy balance: a displaced column with an imbalance would move vertically until the partition was restored. With the hydrostatic pressure $P-{\rho }_{0}gz$ subtracted, the pressure within the current is $\Delta \rho gh$, whereThe equivalent body force is outward for a profile decreasing with $r$. Its magnitude is $\Delta \rho g|h_r|$; the derivative without an absolute value is a signed quantity.

For an ambient-controlled interfacial plug current, integrating the body force through depth gives radial force per area $O(\Delta \rho g{H}^2/R)$. The ambient Stokes flow varies over distance $R$, exerting traction $O(\mu u/R)$, so $u\sim \Delta \rho g{H}^2/\mu$. The current's internal shear supports the ambient traction with a velocity difference $\Delta u\sim (\mu u/R)H/(\lambda \mu )$; $\Delta u/u\sim H/(\lambda R)\ll 1$ requires $\lambda \gg H/R$. The force per area from its in-plane viscous extension is $O(\lambda \mu uH/R^2)$, negligible against ambient traction when $\lambda \ll R/H$. Together these give the viscosity window for an interfacial plug currentFixed volume gives $V\sim H{R}^2$, and the spreading velocity has the scale $\dot{R}$. Thus $\dot{R}\sim \Delta \rho g{V}^2/(\mu R^4)$ andChoosing $A$ so that $R=A{t}^{1/5}$ exactly, the similarity solution has $h=(V/A^2)t^{-2/5}F(\eta )$, $\eta =r/(A{t}^{1/5})$, and $u=A{t}^{-4/5}U(\eta )$. In radial conservation $h_t+r^{-1}(rhu{)}_r=0$, substitution gives $(\eta FU{)}^{\prime }=({\eta }^{2}F{)}^{\prime }/5$. Regularity and no source at the centre remove the integration constant. Therefore the radial plug velocity is linear:

For the disc-traction analogy for an interfacial current, subtract the imposed straining flow from the stationary-disc solution. The perturbation has velocity $-Er{\mathbf{e}}_r$ on the disc and zero far-field velocity. Reverse its sign: the resulting Stokes flow is driven by a radially expanding disc with velocity $Er{\mathbf{e}}_r$ in otherwise quiescent fluid. The background straining flow contributes no radial shear traction on $z=0$, so the given disc traction changes sign. This step uses Linearity of Stokes flow, Uniqueness of Stokes flow, and the no-slip boundary condition. The thin plug current has the same leading boundary velocity and hence the same ambient resistance. Each face contributes $-8\mu Er/[\pi \sqrt{R^2-r^2}]$; both faces must be counted. With $E=1/(5t)$, radial force balance isIntegrate with $h(R,t)=0$. The interfacial plug-current similarity profile isVolume integration gives $V=2\pi {\int }_0^{R}rhdr=4\pi H{R}^2/5$. Eliminating $H$ gives the full prefactor:The formal edge slope diverges, so a narrow edge region lies outside the thin-current approximation. Its traction and volume contributions are integrable, and this is a leading thin-current solution rather than a uniformly valid edge description.

For a diffusion-controlled solute gravity current, excess solute mass, not enriched-fluid volume, is conserved. At late times the effective mass density excess is proportional to solute concentration, so $\Delta \rho (t)R^{2}H\sim \Delta {\rho }_0{V}_0$, and $H\sim (Dt{)}^{1/2}$. The same ambient resistance scaling givesThis scaling does not determine a universal numerical prefactor. The mass diffusivity $D$ controls the growing depth, while dilution weakens buoyancy; assuming a fixed current volume would incorrectly retain the $t^{1/5}$ law. The thin-layer requirement becomes $(\mu D/(\Delta {\rho }_{0}g{V}_0){)}^{1/3}\ll 1$; when only dissolved solute changes the upper fluid, its viscosity is of order the ambient viscosity, consistent with the plug-current window.

Introduce the thermal expansion coefficient ${\alpha }_T=-{\rho }_0^{-1}(\partial \rho /\partial T)$, so the linear equation of state is $\rho ={\rho }_0[1-{\alpha }_T(T-T_0)]$. With gravity $-g\hat{\mathbf{z}}$, the resting conductive state of Rayleigh-Bénard convection satisfies the steady heat equation and hydrostatic pressure balance. Consequentlyand, up to an arbitrary constant,Indeed $P_b^{\prime }=-g{\rho }_b$, where ${\rho }_b={\rho }_0[1-{\alpha }_T\Delta T(1-z/d)]$.

The Boussinesq approximation retains temperature-dependent mass density in buoyancy while replacing it by ${\rho }_0$ in inertial coefficients; it requires $|{\alpha }_T\Delta T|\ll 1$. Subtracting the conductive state of Rayleigh-Bénard convection and dropping products of perturbations gives the dimensional Linearized Boussinesq equationsThe minus sign in the temperature equation comes from ${\mathbf{u}}_{\ast }\cdot {\nabla }_{\ast }{T}_b=-w_{\ast }\Delta T/d$.

Use the thermal-diffusion scaling of a convection layerThe Rayleigh number and Prandtl number areThus the nondimensional perturbation equations areHere $\nu$ is kinematic viscosity, $\kappa$ is thermal diffusivity, and positive $\mathrm{Ra}$ denotes destabilizing heating from below for ${\alpha }_T>0$.

Choose the positive $x$ direction along the background flow, so $U>0$; reverse $x$ if necessary. Let $z$ measure depth below the impermeable roof, and set $\Delta \rho ={\rho }_w-{\rho }_c>0$. The deep ambient water has $p_w=p_0(x)+{\rho }_{w}gz$. Continuity of pressure at $z=h$ and hydrostatic pressure inside the buoyant CO2 giveHere $k$ is permeability of a porous medium and $\mu$ is dynamic viscosity. Since $h/H\ll 1$, the ambient return-flow correction to the imposed pressure gradient is higher order. The small aspect ratio supplies the hydrostatic pressure approximation; capillary trapping, dissolution, compressibility and dispersion are omitted in this sharp-interface model.

The depth-integrated Darcy flux is $q=Uh-K_{b}h{h}_x$. The actual stored fluid volume is $\varphi h$ per unit horizontal area, where $\varphi$ is the aquifer porosity. Thus the porous gravity current with background flow satisfieswith $u=U/\varphi$, $a=K_b/\varphi$ and $J=Q/\varphi$. Both $Q$ and the later $V$ denote actual injected volume per unit transverse width. Away from injection, the equation is an advected porous medium equation. At each finite front impose $h=0$ and vanishing outgoing volume flux per unit width. At the injection point $h$ is continuous and $q(0^{+})-q(0^{-})=Q$. For an initially empty aquifer, $\varphi \int hdx=Qt$. The advective speed is $U/\varphi$, not $U$.

For a constant-flux porous gravity current, write characteristic depth $H_c$ and extent $L_c$ temporarily. At early times volume conservation gives $H_c{L}_c\sim Jt$; gravity spreading gives $L_c^2/t\sim a{H}_c$. HenceBackground advection becomes comparable when $L_c\sim ut$, so the transition scales areFor $t\ll t_c$, the leading similarity solution is symmetric:For $0<\eta <{\eta }_N$, its model ordinary differential equation and boundary conditions areChoose the nonnegative finite-front branch, reflect evenly in $x$, and require $F{F}^{\prime }\to 0$ at the front. The factor $1/2$ divides the injection equally between the two halves. Integration of the ordinary differential equation gives $2{\int }_0^{{\eta }_N}Fd\eta =1$; the limiting front slope is $F^{\prime }({\eta }_N^{-})=-2{\eta }_N/3$. This boundary-value problem determines the initially symmetric shape without assuming an unjustified elementary formula. Independent numerical integration gives ${\eta }_N\simeq 1.17619$ in this normalization.

For $t\gg t_c$, upstream storage approaches a finite steady value. With no upstream leakage, $q=0$ for $x<0$, giving $h_x=U/K_b$ wherever $h>0$. Downstream the uniform interior carries $q=Q$, soThe upstream volume flux vanishes even though local Darcy velocity contributions from background flow and buoyancy cancel. The upstream stored volume is $\varphi h_{\mathrm{int}}|x_{-}|/2$, a finite constant. The downstream plateau therefore extends a distance $ut+O({\ell }_c)$, and its front advances at $u$ to leading order.

The diffusive nose of an advected porous gravity current has widthIn a frame translating at $u$, balance the time derivative against the nonlinear gravity diffusion. With $\eta =(x-x_c(t))/\sqrt{a{h}_{\mathrm{int}}t}$ and $x_c(t)=ut+O({\ell }_c)$, the leading transition is $h=h_{\mathrm{int}}F(\eta )$, withIts finite front has $F\sim ({\eta }_N/2)({\eta }_N-\eta )$ as $\eta \uparrow {\eta }_N$. Independent numerical integration of this boundary-value problem gives ${\eta }_N\simeq 1.23849$, so a more accurate leading nose position is $x_N\simeq x_c+1.23849\sqrt{a{h}_{\mathrm{int}}t}$ and its local shape is $F\simeq 0.619247(1.23849-\eta )$. The entire transition has no elementary closed form. A useful explicit, approximate shape is the volume-preserving linear rampThis gives $x_N\simeq x_c+\sqrt{a{h}_{\mathrm{int}}t}$. Its ramp width is $2\sqrt{a{h}_{\mathrm{int}}t}$, and its edge obeys the exact kinematic boundary condition ${\dot{x}}_N=u-a{h}_x(x_N^{-},t)$ to the displayed order. It also has the same area as a sharp step at $x_c$. This is an approximate ramp closure, not an exact solution of the nonlinear nose ordinary differential equation. A constant phase shift $x_c=ut-|x_{-}|/2$ balances the finite upstream storage within this ramp approximation; the full solution's subleading matching may alter constant offsets.

For a rapid injection of total volume $V$, remove the source for $t>0$ and introduce $\xi =x-ut$. The advected constant-volume porous gravity current obeys $h_t=a(h{h}_{\xi }{)}_{\xi }$ and $\int hd\xi =M=V/\varphi$. A similarity solution has depth proportional to $t^{-1/3}$ and radius proportional to $t^{1/3}$. Direct substitution gives the compactly supported Barenblatt solutionIndeed ${\int }_{-R}^R(R^2-{\xi }^2)/(6at)d\xi =2R^3/(9at)=V/\varphi$. Its centre translates at the pore speed $u$, its two edges are $ut\pm R(t)$, and its maximum depth decays as $t^{-1/3}$. The expression is exact for a point release within this reduced model; finite-width initial data approach this similarity solution at late times. If $U=0$, the constant-rate current remains in the symmetric spreading regime and the constant-volume solution has a stationary centre.

Take $z$ upward from the horizontal bed and measure $x$ away from an ice divide. In the shallow-ice approximation, the hydrostatic pressure is $p=p_a+\rho g(h-z)$ and the horizontal balance for a Newtonian fluid isThe upper stress-free boundary condition gives $u_z(h)=0$. The resisting basal traction is $-\mu \beta u_b$, so the fluid-side shear stress satisfies $\mu u_z(0)=\mu \beta u_b$. This is a Navier slip boundary condition with slip length ${\ell }_s=1/\beta$. Integrating twice givesThe first term is the uniform basal sliding contribution, while the quadratic term is internal viscous deformation. Both are positive where the thickness decreases downstream.

The shallow-ice flux with linear basal drag isNo ice accumulation or ice ablation appears. For a symmetric glacier, work on one half, $0\le x\le x_N(t)$, with $q(0,t)=0$, $h(x_N,t)=0$ and no outgoing volume flux. Its conserved half-volume per unit width is $V_g={\int }_0^{x_N}hdx$; the full glacier has volume $2V_g$. This specifies the volume convention rather than silently supplying an unspecified value. For a general finite release on the whole line, the late similarity solution is centred on its conserved centre of mass and uses half its total volume for $V_g$.

The internal-shear and linear basal drag law mobilities are equal at $H_{\ast }=3/\beta$. Suitable vertical, horizontal and temporal scales areWrite $h=H_{\ast }\mathcal{H}$, $x=L_{\ast }X$, $t=T_{\ast }\tau$. ThenThe initial typical thickness is $\mathcal{H}\simeq 10$, the initial extent is of order $L_{\ast }/10$, and its shear-controlled spreading time is of order $T_{\ast }/10^5$. The shallow-ice approximation also requires small aspect ratio: in particular the initial thickness and extent must satisfy approximately $100H_{\ast }/L_{\ast }\ll 1$. The equations model the broad shallow bulk; the very tip can require physics beyond that approximation.

The two limiting equations are porous medium equations of the form ${\mathcal{H}}_{\tau }=({\mathcal{H}}^n{\mathcal{H}}_X{)}_X$, with $n=3$ for a thick shear-dominated glacier and $n=2$ for a thin sliding-dominated glacier. Volume conservation and balance of the time derivative give exponent $\alpha =1/(n+2)$. Set $\mathcal{H}={\tau }^{-\alpha }f(\eta )$, $\eta =X{\tau }^{-\alpha }$. One integration, using zero volume flux at the ice divide, gives $f^n{f}^{\prime }=-\alpha \eta f$, henceThe planar volume-conserving nonlinear-diffusion similarity is thereforeThe Beta function fixes its constants:Consequently the explicit early shear and late sliding limits areTheir centre depths are $A_3{\tau }^{-1/5}$ with $A_3\simeq 0.842252$, and $A_2{\tau }^{-1/4}$ with $A_2=\sqrt{2/\pi }$. Thus the early extent grows as $t^{1/5}$ and the late extent as $t^{1/4}$. These are asymptotic regime profiles; they solve the respective limiting equations, not the full sum of both mobilities. For arbitrary finite-width initial data they are appropriate after the corresponding spatial relaxation, with a virtual time origin depending on that data.

The shear-to-slip transition in a shallow ice current occurs at thickness of order $H_{\ast }$. Using equality at the centre of the early similarity solution givesIf the initial profile is already approximated by that early similarity solution and its centre depth is 10, the virtual age is ${\tau }_0=A_3^5/10^5\simeq 4.2385\times 10^{-6}$, and the elapsed estimate is ${\tau }_{\mathrm{tr}}-{\tau }_0$. Other reasonable definitions of typical thickness change the order-one coefficient; the initial thickness alone does not uniquely determine an exact transition time or an exact initial profile. The robust nondimensional estimate is a transition time of order one in the $T_{\ast }$ scale.

Finally, let a local nose translate steadily at positive speed $v$ and write $s=x_N-x\ge 0$. This is a local traveling nose of a glacier with linear basal drag; the speed of the globally spreading glacier changes slowly with time. The local mass conservation equation integrates to $q=vh$, since both $q$ and $h$ vanish at the front. Therefore

The right side is strictly increasing, giving a unique positive thickness at every $s>0$. With $S=X_N-X$ and $c=d{X}_N/d\tau$, the same relation reads $cS={\mathcal{H}}^3/3+{\mathcal{H}}^2/2$. Its limiting profiles areThe square-root tip reflects basal sliding; the cube-root thicker region reflects internal viscous deformation. Even during the early thick regime, the very nose is thin and has the square-root inner form. At late times that sliding balance governs nearly the whole profile. The dimensionless crossover thickness is $\mathcal{H}=1$ and its distance is $S=5/(6c)$; the early thick cube-root outer profile matches this smaller sliding region. The diverging geometric slope at the idealized tip also marks the local limit of the long-wave shallow-ice approximation.

Let $z$ increase upwards, so the bed is at $z=-b(x)$, the surface is at $z=h(x,t)$, and $H=h+b$. In the shallow-shelf approximation the incompressible flow has nearly depth-uniform horizontal speed $u$ and vertical strain $w_z=-u_x$. The Newtonian fluid stress tensor, together with the hydrostatic approximation, then givesThe factor four includes both horizontal extension and the pressure correction required by vertical compression. Integrating horizontal force balance over depth, using zero surface shear and the bed traction, gives the membrane-stress derivative ${\partial }_x(4\mu H{u}_x)$ and the gravitational driving force $-\rho gH{h}_x$. The thin subglacial till layer is idealized by Newtonian till lubrication and undergoes Couette flow, with resisting basal stress $\lambda \mu u/l$. Hence force balance and mass conservation yieldSince $b$ is fixed, $H_t=h_t$. No accumulation or ablation is included.

At the grounding line, the first boundary condition is ice-sheet flotation:For an unbuttressed ice shelf, the second is the extensional stress required to balance the difference between the integrated ice and seawater hydrostatic pressures. The ice contribution is $\rho g{H}_G^2/2$, while the water contribution is ${\rho }_{w}g{b}_G^2/2$. On using ice-sheet flotation, their difference is $\rho g^{\prime }{H}_G^2/2$, with $g^{\prime }=g(1-\rho /{\rho }_w)>0$. Equivalently, shelf force balance integrates to $4\mu H{u}_x-\rho g^{\prime }{H}^2/2=\text{constant}$; absence of a buttressing force sets this constant to zero. Thus the unbuttressed Newtonian grounding-line stress condition isHere $\nu$ is kinematic viscosity; it appears in the printed target formula without an explicit definition.

The ratio of longitudinal stress divergence to till drag is of order $4lH/(\lambda L^2)$. The given small-parameter limit therefore yields the friction-dominated bulk relationThe coefficient depends on the basal lubrication, and $h_x<0$ gives seaward motion. This is a bulk reduction: longitudinal stress has been neglected in the differential equation but its boundary traction remains specified. A complete uniformly valid approximation close to the grounding line could require a membrane-stress boundary layer. In the remainder, use the pointwise friction closure and retain the shelf stress condition as the boundary data, as requested in the paper.

Differentiate the ice-sheet flotation condition along the moving grounding line:Now $H_x=h_x+\alpha$. Substitution of the friction-dominated speed and the shelf stress boundary condition gives the friction-dominated grounding-line evolution law:The coefficient of ${\dot{x}}_G$ is the spatial derivative of the flotation deficit, up to sign. If it vanishes, the implicit equation remains the correct compatibility condition but division to determine a finite speed is invalid; tangential contact needs separate analysis.

In a steady state, volume flux per unit width is constant: $Hu=q_0$. The friction relation givesDifferentiating $Hu=q_0$ and using $u_x=g^{\prime }H/(8\nu )$ at the grounding line givesPut $\beta =l/\lambda$, $\mathcal{H}=H/\beta$ and $\mathcal{Q}=\nu q_0/(g{\beta }^3)$. Substitution and cancellation of $g{\beta }^5/\nu$ yield the steady friction-dominated grounding-line thickness relation:For $q_0>0$ and a bed deepening seaward ($\alpha \ge 0$), the polynomial $(g^{\prime }/8g){\mathcal{H}}^5+\alpha \mathcal{Q}{\mathcal{H}}^2-{\mathcal{Q}}^2$ is strictly increasing for $\mathcal{H}>0$, starts negative and tends to infinity, so it has exactly one positive root. At that root $\mathcal{Q}>\alpha {\mathcal{H}}^2$ and $H_x<0$, as required by steady outward extension. The relation belongs to the specified friction closure; it is not a flux law for arbitrary till rheology or arbitrary membrane-stress matching.

Take the layer depth as the length unit, its thermal diffusion time as the time unit, and its imposed temperature difference as the temperature unit. Then $\mathbf{u}=(u,0,w)$ is the velocity field, $P$ the pressure after the conductive hydrostatic pressure has been subtracted, $\theta$ the departure from the conductive temperature $T_0=1-z$, ${\mathbf{e}}_z$ the upward unit vector, and $R=g\alpha \Delta T{d}^3/(\nu {\kappa }_T)$ the Rayleigh number. Here $d$ is the dimensional depth, $\nu$ the kinematic viscosity, ${\kappa }_T$ the thermal diffusivity and $\alpha$ the thermal expansion coefficient. The gradient, divergence and Laplacian act on the dimensionless coordinates. The term $w$ in the heat equation is advection of the background gradient, since $-\mathbf{u}\cdot \nabla T_0=w$. Mass density variations enter only the buoyancy term under the Boussinesq approximation; the incompressible flow constraint is $\nabla \cdot \mathbf{u}=0$.

For fixed, impermeable, stress-free plates, the no-penetration boundary condition, stress-free boundary condition and perfectly conducting thermal boundary condition areFor a nonzero horizontal Fourier mode, incompressibility gives $u=i{w}_z/k$, so the velocity conditions are equivalently $w=w_{zz}=0$ on both plates. A uniform horizontal velocity is not fixed by these boundary conditions; take its mean to be zero when discussing the stationary forced pattern.

The printed small-Prandtl-number label does not generally justify these equations. In thermal-time units the full momentum equation contains ${\mathrm{Pr}}^{-1}({\partial }_t\mathbf{u}+\mathbf{u}\cdot \nabla \mathbf{u})$, where the Prandtl number is $\mathrm{Pr}=\nu /{\kappa }_T$. Dropping that term gives infinite-Prandtl-number convection or another explicitly justified Stokes flow approximation. It is not the general $\mathrm{Pr}\to 0$ limit. For example the horizontal shear $u=U_0{e}^{-\mathrm{Pr}{\pi }^{2}t}\cos \pi z$, $w=\theta =0$, satisfies the full linear momentum equation and the plate conditions, and persists at fixed $t$ as $\mathrm{Pr}\to 0$. The printed inertia-free equation would instead require this shear to vanish. The following parts solve the equations actually displayed; applying them at small Prandtl number needs an additional slaving or creeping-flow ordering.

Take the current to be dilute, vertically well mixed and Boussinesq approximation, with constant channel width $w$, a horizontal absorbing bed, negligible ambient fluid entrainment, negligible bed friction, and no resuspension. The shallow-water approximation neglects vertical inertia and gives hydrostatic pressure. Particle loss alters reduced gravity but does not alter the leading-order carrier-fluid volume or inertia; otherwise additional mass and momentum terms would be necessary. The current occupies depth $h(x,t)$, has depth-uniform speed $u(x,t)$, and has particle volume fraction $c(x,t)$, with $g^{\prime }=g_{o}c$ and constant $g_o>0$.

Volume conservation of a short channel segment gives $h_t+(uh{)}_x=0$. Its suspended particle amount per unit bed area is $hc$, and the downward particle deposition flux is $v_{s}c$. Thus $(hc{)}_t+(uhc{)}_x=-v_{s}c$, which becomes $(h{g}^{\prime }{)}_t+(uh{g}^{\prime }{)}_x=-v_s{g}^{\prime }$ on multiplying by $g_o$.

At vertical coordinate $y$ above the bed, excess hydrostatic pressure divided by reference mass density is $g^{\prime }(h-y)$. Integrating across the depth gives the depth-integrated force flux $g^{\prime }{h}^2/2$. The horizontal momentum per unit bed area is $hu$ and the advected kinematic momentum flux is $h{u}^2$. The resulting variable-buoyancy shallow water equations areThe derivative of the pressure term acts on both $h$ and $g^{\prime }$; treating reduced gravity as horizontally constant while particles settle would lose part of the force.

For a steady interior, let $q=Q/w=uh>0$ denote volume flux per unit width. Then the particle balance is a first-order linear differential equation,Here $g_0^{\prime }=g^{\prime }(0)>0$, $v_s\ge 0$ is constant, and deposited particles are not returned to the current. In particular the decay length is $q/v_s=Q/(w{v}_s)$ for $v_s>0$, and the cancellation of depth follows from $uh=q$: dilution and transit time compensate in the deposition balance.

The steady momentum equation also givesThe inlet depth selects $K=q^2/h_0+g_0^{\prime }{h}_0^2/2$. Depth and speed at every downstream position are determined by the continuous positive root of a polynomial inThis implicit expression also makes clear why the inlet condition is needed: there are generally two possible positive depths.

The depth branches of a steady depositing gravity current are distinguished by the Froude number ${\mathrm{Fr}}^2=q^2/(g^{\prime }{h}^3)$. Differentiating the momentum invariant, and using the exponential decay of reduced gravity, givesFor $v_s>0$, supercritical flow has $\mathrm{Fr}>1$ and decreasing depth, with $h\to q^2/K$ and $u\to K/q$ as $x\to \infty$. Subcritical flow has $\mathrm{Fr}<1$ and increasing depth, with $h\sim \sqrt{2K/g^{\prime }(x)}\propto e^{v_{s}x/(2q)}$. This indefinitely deepening branch eventually violates the assumed shallow or confined-channel geometry. Hence the depth does not universally decrease: that conclusion needs a supercritical inlet.

For fixed $g^{\prime }>0$, $q^2/h+g^{\prime }{h}^2/2$ has its unique minimum at $h_c=(q^2/g^{\prime }{)}^{1/3}$, with $K_c=3q^{4/3}(g^{\prime }{)}^{1/3}/2$. The two positive roots when $K>K_c$ lie below and above $h_c$. Since $K_c$ decreases downstream while $K$ is constant, an initially strict branch remains separate from the other branch and cannot cross a critical point. A critical inlet with $v_s>0$ has no finite depth derivative: at $\mathrm{Fr}=1$ the first differentiated equation would equate zero to a positive value. It requires a singular entrance adjustment or physics beyond this steady model. If $v_s=0$, both reduced gravity and the continuous depth branch are constant. These steady solutions describe the supplied current interior, not a moving finite-release front.

For two particle populations, write $g^{\prime }=g_1^{\prime }+g_2^{\prime }$ and $g_i^{\prime }=g_{0i}{c}_i$, allowing different constant particle-density conversion factors. The bidisperse particle-laden gravity current equations areFor the steady case,Thus the depth is obtained from the same cubic polynomial invariant with the sum of the two exponentially decreasing reduced gravity contributions.

The local deposit consists of the two particle deposition fluxes $v_1{c}_1$ and $v_2{c}_2$. For a positive total deposition rate, the bidisperse gravity-current deposition fraction by particle volume isIf the particles have a common material density, this is also their local mass fraction. If they have common $g_o$, it can be written in the source's buoyancy notation asFor unequal $g_{0i}$ replace each $g_{i0}^{\prime }$ by $g_{i0}^{\prime }/g_{0i}$. The exponent sign matters: with both populations present and $v_1>v_2>0$, the faster-settling population has decreasing local deposit fraction downstream, since $f_1^{\prime }=(v_2-v_1)f_1(1-f_1)/q<0$. Equal positive settling speeds give a constant fraction. A nonsettling population contributes no deposit; if the entire local deposition rate is zero the deposit fraction is undefined.