Abstract: The Richard’s equation was used to numerically model 1D flow within unsaturated topsoil and solved using a two-step Crank-Nicolson scheme. Under constant water supply, the Darcy velocity obtained at the upper boundary, the ...