Search

[Display omitted] • Saturated hydraulic conductivity (K s) maps are essential for land use planning. • K s mapping at the catchment scale is very expensive and time-consuming. • Variable selection, machine learning and geostatistics are used for K s mapping. • The spatial variability of K s is addressed using a "double-scale" approach. • The map accounts for both anthropic and natural factors of K s spatial variability. Saturated hydraulic conductivity (K s) mapping is an application tool in many fields of science, technology, and land management. Hydrologically, K s is the main factor controlling aquifer recharge, but it is one of the more difficult hydraulic properties of the soil to predict. This work develops a novel methodology for topsoil K s mapping at the catchment scale, based on auxiliary data that are quick to determine and low-cost. It includes a double-scale sampling of the K s to account for small-scale variability in the spatial geostatistical interpolation. A supervised selection of variables through correlation analysis and hierarchical clustering of variables precedes the development of site-specific pedotransfer functions (PTFs) with machine learning techniques. The approach was applied to a heavily anthropized area subject to flooding on the island of Mallorca, Balearic Islands, Spain. The variable selection process filtered out four predictor variables from the initial pool of fourteen predictors. An artificial neural network (ANN) with one hidden layer and six input variables (latitude, longitude, silt, clay, medium sand, and land use) provided the best K s prediction model. Latitude and longitude coordinates and land use are surrogates for other physical and environmental factors. This ANN was found to be much more accurate than other classical PTFs. The ANN was used to estimate new K s values from a subsequent sampling of model covariates, which allowed doubling the input information for spatial interpolation using ordinary kriging (OK). The consistency of the resulting K s map was evaluated by direct comparison with a geomorphological base map. Overall, the spatial distribution of K s was consistent with the lithological variability as well as with other superimposed anthropic factors. [ABSTRACT FROM AUTHOR]