Inverse estimation of hydraulic parameters of soils with rock fragments

The application of inverse modeling to determine the hydraulic properties of layered soils with rock fragments (RF > 2 mm) is associated with complex challenges such as selecting suitable hydraulic functions and defining a strategy to optimize many parameters. These two issues were addressed in this study by performing field drainage experiments in layered soil profiles (three layers) with different particle size distributions and different percentages of RF. Water contents (& theta;) and pressure heads (h) were monitored at different soil depths during drainage experiments (h & GE;-200 cm). Data of & theta; and h for a drier range (h <-5000 cm) were determined with a dew point potentiometer on disturbed soil samples. The unimodal and bimodal van Genuchten-Mualem functions were evaluated, and their parameters were optimized using the Hydrus-1D Inverse solution. The RETC curve fitting software was used to estimate the second set of parameters of the bimodal function. The parameters for the uni-and bimodal models (15 and 27 parameters, respectively) were repeatedly optimized by fixing some pa-rameters while estimating others. The unimodal van Genuchten-Mualem function provided a satisfactory fit of & theta; and h measurements only when either drainage or dew point potentiometer measurements were used. On the other hand, the bimodal van Genuchten-Mualem function provided a satisfactory fit of & theta; and h measurements when both sets of data (drainage and dew point potentiometer measurements) were used simultaneously. Furthermore, the optimized parameters were physically consistent except for a few high saturated hydraulic conductivity values. Therefore, bimodal functions should be considered to represent the hydraulic properties of soils with RF. The problem of many calibrated parameters when using bimodal functions can be managed using sequential optimization, which allowed us to estimate 27 parameters for layered soils with RF successfully.

Automated summary

This study addressed the complex challenges associated with the application of inverse modeling in determining the hydraulic properties of layered soils with rock fragments (RF > 2 mm). Field drainage experiments were conducted in layered soil profiles with different particle size distributions and percentages of rock fragments to optimize the parameters of hydraulic functions. The bimodal van Genuchten-Mualem function provided a satisfactory fit of water contents and pressure heads measurements when both drainage and dew point potentiometer measurements were used simultaneously.