Parasite inversion for determining the coefficients and time-validity of Philip's two-term infiltration equation

Many different equations have been proposed to describe quantitatively one-dimensional soil water infiltration. The unknown coefficients of these equations characterize soil hydraulic properties and may be estimated from a n record, {t similar to i,I similar to i}i=1n, of cumulative infiltration measurements using curve fitting techniques. The two-term infiltration equation, I(t)=St+cKst, of Philip has been widely used to describe measured infiltration data. This function enjoys a solid mathematical-physical underpinning and admits a closed-form solution for the soil sorptivity, S [L T-1/2], and multiple, c [-], of the saturated hydraulic conductivity, K-s [L T-1]. However, Philip's two-term equation has a limited time validity, t(valid) [T], and thus cumulative infiltration data, I similar to(t similar to), beyond t=tvalid will corrupt the estimates of S and K-s. This paper introduces a novel method for estimating S, c, K-s, and t(valid) of Philip's two-term infiltration equation. This method, coined parasite inversion, use as vehicle Parlange's three-parameter infiltration equation. As prerequisite to our method, we present as secondary contribution an exact, robust and efficient numerical solution of Parlange's infiltration equation. This solution admits Bayesian parameter estimation with the DiffeRential Evolution Adaptive Metropolis (DREAM) algorithm and yields as byproduct the marginal distribution of Parlange's beta parameter. We evaluate our method for 12 USDA soil types using synthetic infiltration data simulated with HYDRUS-1D. An excellent match is observed between the inferred values of S and K-s and their "true" values known beforehand. Furthermore, our estimates of c and t(valid) correlate well with soil texture, corroborate linearity of the c(beta) relationship for 0 <= t <= tvalid, and fall within reported ranges. A cumulative vertical infiltration of about 2.5 cm may serve as guideline for the time-validity of Philip's two-term infiltration equation.