Estimation of hydraulic conductivity and its uncertainty from grain-size data using GLUE and artificial neural networks

peer reviewedaudience: researcher, professionalVarious approaches exist to relate saturated hydraulic conductivity (Ks) to grain-size data. Most methods use a single grain-size parameter and hence omit the information encompassed by the entire grain-size distribution. This study compares two data-driven modelling methods, i.e.multiple linear regression and artificial neural networks, that use the entire grain-size distribution data as input for Ks prediction. Besides the predictive capacity of the methods, the uncertainty associated with the model predictions is also evaluated, since such information is important for stochastic groundwater flow and contaminant transport modelling. Artificial neural networks (ANNs) are combined with a generalized likelihood uncertainty estimation (GLUE) approach to predict Ks from grain-size data. The resulting GLUE-ANN hydraulic conductivity predictions and associated uncertainty estimates are compared with those obtained from the multiple linear regression models by a leave-one-out cross-validation. The GLUE-ANN ensemble prediction proved to be slightly better than multiple linear regression. The prediction uncertainty, however, was reduced by half an order of magnitude on average, and decreased at most by an order of magnitude. This demonstrates that the proposed method outperforms classical data-driven modelling techniques. Moreover, a comparison with methods from literature demonstrates the importance of site specific calibration. The dataset used for this purpose originates mainly from unconsolidated sandy sediments of the Neogene aquifer, northern Belgium. The proposed predictive models are developed for 173 grain-size -Ks pairs. Finally, an application with the optimized models is presented for a borehole lacking Ks data

Estimating Infiltration Parameters from Basic Soil Properties

Infiltration data were collected on two rectangular grids with 25 sampling points each. Both experimental grids were located in tropical rain forest (Guyana), the first in an Arenosol area and the second in a Ferralsol field. Four different infiltration models were evaluated based on their performance in describing the infiltration data. The model parameters were estimated using non‐linear optimization techniques. The infiltration behaviour in the Ferralsol was equally well described by the equations of Philip, Green–Ampt, Kostiakov and Horton. For the Arenosol, the equations of Philip, Green–Ampt and Horton were significantly better than the Kostiakov model. Basic soil properties such as textural composition (percentage sand, silt and clay), organic carbon content, dry bulk density, porosity, initial soil water content and root content were also determined for each sampling point of the two grids. The fitted infiltration parameters were then estimated based on other soil properties using multiple regression. Prior to the regression analysis, all predictor variables were transformed to normality. The regression analysis was performed using two information levels. The first information level contained only three texture fractions for the Ferralsol (sand, silt and clay) and four fractions for the Arenosol (coarse, medium and fine sand, and silt and clay). At the first information level the regression models explained up to 60% of the variability of some of the infiltration parameters for the Ferralsol field plot. At the second information level the complete textural analysis was used (nine fractions for the Ferralsol and six for the Arenosol). At the second information level a principal components analysis (PCA) was performed prior to the regression analysis to overcome the problem of multicollinearity among the predictor variables. Regression analysis was then carried out using the orthogonally transformed soil properties as the independent variables. Results for the Ferralsol data show that the parameters of the Green–Ampt and Kostiakov model were estimated relatively accurately (maximum R2 = 0.76). For the Arenosol, use of the second information level together with PCA produced regression models with an R2 value ranging from 0.38 to 0.68. For the Ferralsol, most of the variance was explained by the root content and organic matter content. In the Arenosol plot, the fractions medium and fine sand explained most of the observed variance