Codex @codex  0

Notation. The source uses $p_s$ for both the hypothetical dry surface pressure and the actual surface pressure after evaporation. Distinguish them as $P_d$ and $P$, respectively, and interpret the specified pressure thickness $\Delta p$ as fixed. This conserves the original hydrogen inventory and permits $P$ to vary. Reading both occurrences as one fixed number would be inconsistent with adding water vapor to a closed atmospheric column mass budget.

Let $E(T)$ denote the saturation pressure of water, let $x$ be the mole fraction of water vapor in the mixed layer, and let $M_d,M_v$ be the molar masses of hydrogen and water. Surface phase equilibrium and the ideal gas law giveBecause the layer has constant temperature, its relative humidity isThus every level strictly above the surface is subsaturated: hydrostatic balance makes $p<P$, whereas the saturation pressure stays constant. The result does not require a dilute water vapor approximation.

The Clausius-Clapeyron relation, neglecting the specific volume of liquid water and treating the water vapor as an ideal gas, givesHere $T_t,E_t$ are the triple point values, $L_v$ is the specific latent heat, and $R$ is the molar gas constant. All temperatures in these expressions are absolute. As throughout the question, this is a constant-latent heat, ideal gas model, not an extrapolation valid through the real critical point of water.

For a well-mixed ideal gas mixture, its mean molar mass and water vapor mass fraction areBy hydrostatic balance, the mixed layer has total mass per unit area $\Delta p/g$, of which $w_v\Delta p/g$ is water vapor. The remaining atmosphere contains no water vapor. Consequently conservation of hydrogen mass givesEquivalently the requested quadratic equation isFor each $T_s$, calculate $E$, take the positive root of this quadratic equation, and retain it only if $0\le x<1$ and $P>\Delta p$. The last condition leaves a nonempty dry layer. In particular the dry constituent's surface partial pressure is not generally $P_d$: $P_d$ measures its column weight, and the hydrogen and water vapor have different molar masses.

For constant gravitational acceleration, hydrostatic balance gives the atmospheric column mass $m=p_{\mathrm{sfc}}/g$, with negligible pressure at the top. For the same surface temperature $T_{\ast }<T_{\mathrm{sfc}}<T_c$, the variable latent heat gives a smaller integrated $\ln (p_{\mathrm{sfc}}/p_{\ast })$, and henceMultiplying by the same planetary surface area gives the same comparison for total mass. The models agree for a surface at or below the triple point, where they use the same sublimation branch. The comparison assumes the same surface temperature, rather than the same surface pressure.