D. Mallants et al., EVALUATION OF MULTIMODAL HYDRAULIC FUNCTIONS IN CHARACTERIZING A HETEROGENEOUS FIELD SOIL, Journal of hydrology, 195(1-4), 1997, pp. 172-199
Soil water retention curves are often used to estimate the hydraulic c
onductivity function. Unfortunately, single S-shaped functions cannot
adequately describe water retention curves of structured soil, especia
lly near saturation. The approach of superposition of two or more unim
odal retention functions such as the van Genuchten model was used here
to describe retention data of a macroporous soil. A total of 180 core
s, 0.05 m diameter and 0.051 m long, were sampled along a 31-m-long tr
ansect in three overlying soil horizons. Use of unimodal retention cur
ves leads to an underestimation of observed water contents both near s
aturation and in the midpore range, while an overestimation is found i
n the drier range. Superposition of two unimodal retention curves sign
ificantly improved the estimation over the entire pressure range. Howe
ver, the predictions were still not ideal near saturation, With three
unimodal curves, a perfect fit was obtained from saturation to residua
l water content, Most of the multimodal parameter values were moderate
ly heterogeneous along the transect, with the surface horizon slightly
more heterogeneous than the deeper layers. The coefficient of variati
on (CV) for multimodal parameters was generally in the range of 20 to
70%. Use of the multimodal van Genuchten model with the conductivity e
stimation model of Mualem resulted in conductivities that were general
ly much smaller than those estimated by the classical unimodal van Gen
uchten-Mualem model. A preliminary evaluation of the estimated bimodal
and trimodal unsaturated hydraulic conductivity model was based on a
comparison with independent conductivity measurements using a combinat
ion of crust test, hot-air method, and an unsteady drainage flux exper
iment on large columns. The crust and hot-air data compared best with
the estimated trimodal conductivity function. The unsteady drainage da
ta did not match well with the crust and hot-air data and could not be
described with any of the estimated conductivity functions.