Numerically Predicting Seepage Gradient Forces and Erosion: Sensitivity to Soil Hydraulic Properties

Authors

Garey A. Fox, Oklahoma State University - Main CampusFollow
Derek M. Heeren, University of Nebraska-LincolnFollow
Glenn V. Wilson, USDA National Sedimentation LaboratoryFollow
E. J. Langendoen, USDA-ARS National Sedimentation LaboratoryFollow
Amanda K. Fox, Stillwater, OK
Maria L. Chu-Agor, University of Florida

Date of this Version

2010

Citation

Published in Journal of Hydrology 389 (2010) 354–362. DOI:10.1016/j.jhydrol.2010.06.015

Comments

US government work

Abstract

Research has suggested that streambank seepage can be an important mechanism of bank instability; however, limited information is available on the level of soil characterization necessary to accurately predict seepage gradient forces and erosion. The objective of this research was to quantify the expected range of predicted seepage gradients for various degrees of site characterization. Uncertainty analysis on seepage gradient predictions was performed relative to variability in soil hydraulic properties. A two-dimensional unsaturated/saturated groundwater flow model was used to simulate a homogeneous soil layer for sand and loamy sand soils packed at various bulk densities, ρ_b. A pedotransfer function (ROSETTA), designed to estimate unsaturated hydraulic properties from surrogate soil data (i.e., texture and bulk density), was used to derive the saturated hydraulic conductivity, K_s, and water retention parameters for various levels of site information (i.e., only textural class; percent sand, silt, and clay (%SSC); %SSC and ρ_b; and %SSC, ρ_b, and K_s). Statistical distributions were derived for each soil hydraulic parameter and Monte Carlo simulations were performed to generate distributions of maximum seepage gradient. The deviation in predicted seepage gradient was calculated using assumed baseline conditions. Ranges in predicted soil hydraulic parameters and maximum seepage gradients were considerably reduced when using %SSC as compared to soil texture. Therefore, at a minimum, soil samples should be taken for particle size analysis. For ρ_b between 1450 and 1500 kg m^-3, soil hydraulic parameters could be derived using ROSETTA and inputting %SSC, with little additional benefit provided by measuring ρ_b and/or K_s. When the ρ_b was less than 1450–1500 kg m^-3, inputting ρ_b and/or K_s consistently reduced the magnitude of deviations from the baseline and therefore should be measured from undisturbed soil samples. The opposite was observed for ρ_b greater than 1450–1500 kg m^-3 due to discrepancies between ROSETTA-derived and actual values of soil hydraulic parameters other than K_s. Considerable deviations (i.e., around 20%) were observed in seepage gradients under this scenario. When ROSETTA-derived and actual values of soil hydraulic parameters more closely matched, inputting ρ_b and/or K_s benefitted seepage gradient predictions as deviations in seepage gradients were less than 5% for the sand and loamy sand soils. Therefore, it is vital to quantify all soil hydraulic parameters for high ρ_b soils and textures with a wide range in %SSC.