Prediction of unsaturated flow and water backfill during infiltration in layered soils

We develop a new analytical infiltration model to determine water flow dynamics around layer interfaces during infiltration process in layered soils. The model mainly involves the analytical solutions to quadratic equations to determine the flux rates around the interfaces. Active water content profile behind the wetting front is developed based on the solution of steady state flow to dynamically update active parameters in sharp wetting front infiltration equations and to predict unsaturated flow in coarse layers before the front reaches an impeding fine layer. The effect of water backfill to saturate the coarse layers after the wetting front encounters the impeding fine layer is analytically expressed based on the active water content profiles. Comparison to the numerical solutions of the Richards equation shows that the new model can well capture water dynamics in relation to the arrangement of soil layers. The steady state active water content profile can be used to predict the saturation state of all layers when the wetting front first passes through these layers during the unsteady infiltration process. Water backfill effect may occur when the unsaturated wetting front encounters a fine layer underlying a coarse layer. Sensitivity analysis shows that saturated hydraulic conductivity is the parameter dictating the occurrence of unsaturated flow and water backfill and can be used to represent the coarseness of soil layers. Water backfill effect occurs in coarse layers between upper and lower fine layers when the lower layer is not significantly coarser than the upper layer.