Field-scale spatial variability of saturated hydraulic conductivity on a recently constructed artificial ecosystem

Abstract

Saturated hydraulic conductivity ( K s) influences water storage and movement, and is a key parameter of water and solute transport models. Systematic field evaluation of K s and its spatial variability for recently constructed artificial ecosystems is still lacking. The objectives of the present study were; (1) to determine saturated hydraulic conductivity of an artificial ecosystem using field methods (Philip–Dunne, and Guelph permeameters), and compare their results to the constant-head laboratory method; (2) to evaluate the spatial variability of K s using univariate and geostatistical analyses, and (3) to evaluate the ability of five pedotransfer functions to predict K s. The results showed that K s varied significantly (p < 0.05) among methods, probably reflecting differences in scales of measurement, flow geometry, assumptions in computation routines and inherent disturbances during sampling. Mean K s values were very high for all methods (38.6–77.9 m day − 1 ), exceeding values for natural sandy soils by several orders of magnitude. The high K s values and low coefficients of variation (26–44%) were comparable to that of well-sorted unconsolidated marine sands. Geostatistical analysis revealed a spatial structure in surface K s data described by a spherical model with a correlation range of 8 m. The resulting kriged map of surface K s showed alternating bands of high and low values, consistent with surface structures created by wheel tracks of construction equipment. Vertical K s was also spatially structured, with a short correlation range of 40 cm, presumably indicative of layering caused by post-construction mobilization and deposition of fine particles. K s was linearly and negatively correlated with dry soil bulk density (ρ b) (r 2 = 0.73), and to a lesser extent silt plus clay percentage ( Si + C) (r 2 = 0.21). Combining both ρ b and Si + C significantly (p < 0.05) improved the relationship and gave the best predictor of K s (r 2 = 0.76). However, evaluation of five PTFs developed for natural soils showed that they all underestimated K s by an order of magnitude, suggesting that application of water balance simulation models based on such PTFs to the present study site may constitute a bias in model outputs. Overall, the study demonstrated the influence of material handling, construction procedures and post-construction processes on the magnitude and spatial variability of K s on a recently constructed artificial ecosystem. These unique hydraulic properties may have profound impacts on soil moisture storage, plant water relations and water balance fluxes on artificial ecosystems, particularly where such landforms are intended to restore pre-disturbance ecological and hydrological functions.