Table 1 Comparison of saturation values of ψ for WMS model: ψ^sat = ψ₂ (12) and (7), ψ^sat = ψ₃ (13), and ψ^sat = ψ₄, where ψ₄ is estimated from the solution of (8)-(11) and (7) at t = 10 000 s

k ₁

k ₂

ψ ₂

ψ ₃

ψ ₄

k₃ = 0

k₃ = 0.000033