Identifying Covariates to Assess the Spatial Variability of Saturated Soil Hydraulic Conductivity Using Robust Cokriging at the Watershed Scale

Abstract

The mapping of saturated soil hydraulic conductivity (K Sat) is essential to understanding soil water dynamics and is a sensitive input in hydrological modeling. The objectives of this study were to provide a reference for the selection of soil hydrology and other environmental attributes that can be used as covariates for estimating K Sat and to compare the efficiency of univariate ordinary kriging versus ordinary robust cokriging, using selected soil hydrology and environmental attributes. Data sets were obtained from a sample grid of 179 points established in the Ellert creek watershed (ECW), located in Rio Grande do Sul state, Southern Brazil. K Sat , macroporosity, microporosity, total porosity, and bulk density were determined from soil sampled at each point. Data of land use and elevation were also applied. All data sets were firstly submitted to classical statistics. Boxplot graphics were constructed to evaluate the relationship between K Sat and land uses. Spearman coefficient of correlation between K Sat and the other attributes was also assessed. For the assortment of covariates, cluster analysis was applied. Classical and robust estimators were applied to calculate the auto and cross-semivariograms and hereafter the ordinary kriging and cokriging. The Spearman coefficient showed some inconsistencies among the applied variables, suggesting that the multivariate method was more appropriate. All cross-semivariograms, except for land use, showed results with better accuracy than the auto-semivariograms. From the methods applied, the best estimates of K Sat were obtained using the robust cokriging method, using macroporosity and soil bulk density as covariates.

Hydraulic Conductivity in a Soil Cultivated with Wheat-Rapeseed Rotation Under Two Tillage Systems

“…Villarreal et al (2020) concluded the same, but indicate that in soils with low clay content the differences between tillage systems disappear. Alletto and Coquet (2009) found a gradual decrease of K s through time, and in general lower values in the wheel track than outside, depending on variations of coarse porosity and bulk density, as was observed by Soares et al (2020) at a large-scale approach. Wahl et al (2004) found a correlation of the infiltration values with the macroporosity depending on the dynamic of CT in time, being the initial mechanical loosening and subsequent soil settlement the main factors during the season.…”