Hydrology and Earth System Sciences www.hydrol-earth-syst-sci.net 1027-5606 1607-7938 13 8 2009 10.5194/hess-13-1503-2009 http://www.hydrol-earth-syst-sci.net/13/1503/2009/ http://www.hydrol-earth-syst-sci.net/13/1503/2009/hess-13-1503-2009.html http://www.hydrol-earth-syst-sci.net/13/1503/2009/hess-13-1503-2009.pdf 1503 1518 2009-08-20 Simulation and validation of concentrated subsurface lateral flow paths in an agricultural landscape Q. Zhu H. S. Lin henrylin@psu.edu The Pennsylvania State University, Department of Crop and Soil Sciences, 116 Agricultural Sciences and Industry Building, University Park, PA 16802, USA The importance of soil water flow paths to the transport of nutrients and contaminants has long been recognized. However, effective means of detecting concentrated subsurface flow paths in a large landscape are still lacking. The flow direction and accumulation algorithm based on single-direction flow algorithm (D8) in GIS hydrologic modeling is a cost-effective way to simulate potential concentrated flow paths over a large area once relevant data are collected. This study tested the D8 algorithm for simulating concentrated lateral flow paths at three interfaces in soil profiles in a 19.5-ha agricultural landscape in central Pennsylvania, USA. These interfaces were (1) the interface between surface plowed layers of Ap1 and Ap2 horizons, (2) the interface with subsoil water-restricting clay layer where clay content increased to over 40%, and (3) the soil-bedrock interface. The simulated flow paths were validated through soil hydrologic monitoring, geophysical surveys, and observable soil morphological features. The results confirmed that concentrated subsurface lateral flow occurred at the interfaces with the clay layer and the underlying bedrock. At these two interfaces, the soils on the simulated flow paths were closer to saturation and showed more temporally unstable moisture dynamics than those off the simulated flow paths. Apparent electrical conductivity in the soil on the simulated flow paths was elevated and temporally unstable as compared to those outside the simulated paths. The soil cores collected from the simulated flow paths showed significantly higher Mn content at these interfaces than those away from the simulated paths. These results suggest that (1) the D8 algorithm is useful in simulating possible concentrated subsurface lateral flow paths if used with appropriate threshold value of contributing area and sufficiently detailed digital elevation model (DEM); (2) repeated electromagnetic surveys can reflect the temporal change of soil water storage and thus is a useful indicator of possible subsurface flow path over a large area; and (3) observable Mn distribution in soil profiles can be used as a simple indicator of water flow paths in soils and over the landscape; however, it does require sufficient soil sampling (by excavation or augering) to possibly infer landscape-scale subsurface flow paths. In areas where subsurface interface topography varies similarly with surface topography, surface DEM can be used to simulate potential subsurface lateral flow path reasonably so the cost associated with obtaining depth to subsurface water-restricting layer can be minimized. Auerswald, K., Simon, S., and Stanjek, H.: Influence of soil properties on electrical conductivity under humid water regimes, Soil. Sci., 166, 382–390, 2001. Bakhsh, A. and Kanwar, R. S.: Soil and landscape attributes interpret subsurface drainage clusters, Aust. J. Soil. Res., 46, 735–744, 2008. Birkhead, A. L., Heritage, G. L., White, H., and van Niekerk, A. W.: Ground-penetrating radar as a tool for mapping the phreatic surface, bedrock profile, and alluvial stratigraphy in the Sabie River, Kruger National Park, J. Soil Water Cons., 51, 234–241, 1996. Bogaart, P. W. and Troch, P. A.: Curvature distribution within hillslopes and catchments and its effect on the hydrological response, Hydrol. Earth Syst. Sci., 10, 925–936, 2006. Burns, D. A., Hooper, R. P., McDonnell, J. J., Freer, J. E., Kendall, C., and Beven, K.: Base cation concentrations in subsurface flow from a forested hillslope: The role of flushing frequency, Water Resour. Res., 34, 3535–3544, 1998. Buttle, J. M. and McDonald, D. J.: Coupled vertical and lateral preferential flow on a forested slope, Water Resour. Res., 38, 1060, doi:10.1029/2001WR000773, 2002. Cassel, D. K., Afyuni, M. M., and Robarge, W. P.: Manganese distribution and patterns of soil wetting and depletion in a piedmont hillslope, Soil Sci. Soc. Am. J., 66, 939–947, 2002. Chen, S. K., Liu, C. W., and Huang, H. C.: Analysis of water movement in paddy rice fields (II) simulation studies, J. Hydrol., 268, 259–271, 2002. Corwin, D. L. and Lesch, S. M.: Apparent soil electrical conductivity measurements in agriculture, Comput. Electron. Agr., 46, 11–43, 2005. De Lannoy, G. J. M., Verhoest, N. E. C., Houser, P. R., Gish, T. J., and van Meirvenne, M.: Spatial and temporal characteristics of soil moisture in an intensively monitored agricultural field (OPE3), J. Hydrol., 331, 719–730, 2006. Elliot, J. A., Cessna, A. J., Best, K. B., Nicholaichuk, W., and Tollefson, L. C.: Leaching and preferential flow of clopyralid under irrigation: field observations and simulation modeling, J. Environ. Qual., 27, 124–131, 1998. Erskine, R. H., Green, T. R., Ramirez, J. A., and MacDonald, L. H.: Comparison of gridbased algorithms for computing upslope contributing area, Water Resour. Res., 42, W09416, doi:10.1029/2005WR004648, 2006. Fairfield, J. and Leymarie, P.: Drainage networks from grid digital elevation models, Water Resour. Res., 27, 709–717, 1991. Fiori, A., Romanelli, M., Cavalli, D. J., and Russo, D.: Numerical experiments of streamflow generation in steep catchments, J. Hydrol., 339, 183–192, 2007. Freer, J., McDonnell, J. J, Beven, K. J., Brammer, D., Burns, D. A, Hooper, R. P., and Kendall, C.: Topographic controls on subsurface storm flow at the hillslope-scale for two hydrologically distinct small catchments, Hydrol. Processes., 11, 1347–1352, 1997. Freer, J., McDonnell, J. J, Beven, K. J., Peters, N. E., Burns, D. A, Hooper, R. P., Aulenbach, B., and Kendall, C.: The role of bedrock topography on subsurface storm flow, Water Resour. Res., 38, 1269, doi:10.1029/2001WR000872, 2002. Gish, T. J., Dulaney, W. P., Kung, K.-J. S., Daughtry, C. S. T., Doolittle, J. A., and Miller, P. T.: Evaluating use of ground-penetrating radar for identifying subsurface flow pathways, Soil Sci. Soc. Am. J., 66, 1620–1629, 2002. Gish, T. J., Walthall, C. L., Daughtry, C. S. T., and Kung, K.-J. S.: Using soil moisture and spatial yield patterns to identify subsurface flow pathways, J. Environ. Qual., 34, 274–286, 2005. Guo, J. H., Liang, X., and Leung, L. R.: A new multiscale flow network generation scheme for land surface models, Geophys. Res. Lett., 32, L23502, doi:10.1029/2004GL021381, 2004. Haga, H., Matsumoto, Y., Matsutani, J., Fujita, M., Nishida, K., and Sakamoto, Y.: Flow paths, rainfall properties, and antecedent soil moisture controlling lags to peak discharge in a granitic unchanneled catchment, Water Resour. Res., 41, W12410, doi:10.1029/2005WR004236, 2005. Haria, A. H., Johnson, A. C., Bell, J. P., and Batchelor, C. H.: Water-movement and isoproturon behaviour in a drained heavy clay soil. 1. Preferential flow processes, J. Hydrol., 163, 201–216, 1994. Heppell, C. M., Burt, T. P., and Williams, R. J.: Variations in the hydrology of an underdrained clay hillslope, J. Hydrol., 227, 236–256, 2000. Jones, R.: Algorithms for using a DEM for mapping catchment areas of stream sediment samples, Comput. Geosci., 28, 1051–1060, 2002. Ju, S. H. and Kung, K.-J. S.: Finite element simulation of funnel flow and overall flow property induced by multiple soil layers, J. Environ. Qual., 22, 432–442, 1993. Kenny, F., Matthews, B., and Todd, K.: Routing overland flow through sinks and flats in the interpolated raster terrain surfaces, Comput. Geosci., 34, 1417–1430, 2008. Kettler, T. A., Doran, J. W., and Gilbert, T. L.: Simplified method for soil particle-size determination to accompany soil-quality analyses, Soil Sci. Soc. Am. J., 65, 849–852, 2001. Kirkby, M. J.: Thresholds and instability in stream head hollows: a model of magnitude and frequency for wash processes, in: Process Models and Theoretical Geomorphology, edited by: Kirkby, M. J., Wiley, Chichester, UK, 295–314, 1994. Kitahara, H., Terajima, T., and Nakai, Y.: Ratio of pipe flow to throughflow discharge, J. Jpn. Forest. Soc., 76, 10–17, 1994. Krovchenko, A. N.: Influence of spatial structure on accuracy of interpolation methods, Soil Sci. Soc. Am. J., 67, 1564–1571, 2003. Kravchenko, A. N. and Robertson, G. P.: Can topographical and yield data substantially improve total soil carbon mapping by regression kriging? Agron. J., 99, 12–17, 2007. Kung, K.-J. S.: Preferential flow in a sandy vadose zone: 1. Field observation, Geoderma, 46, 30 51–71, 1990. Kung, K.-J. S.: Laboratory observation of the funnel flow mechanism and its influence on solute transport, J. Environ. Qual., 22, 91–102, 1993. Lin, H. S., Bouma, J., Wilding, L. P., Richardson, J.L., Kutilek, M., and Nielsen, D. R.: Advances in hydropedology, Adv. Agron., 85, 1–89, 2005. Lin, H. S., Kogelmann, W., Walker, C., and Bruns, M. A.: Soil moisture patterns in a forested catchment: A hydropedological perspective, Geoderma, 131, 345–368, 2006. Lin, H. S., Brook, E., McDaniel, P., and Boll, J.:Hydropedology and Surface/Subsurface Runoff Processes, in: Encyclopedia of Hydrologic Sciences.edited by: Anderson, M. G., John Wiley & Sons, Ltd., doi:10.1002/0470848944.hsa306, 1–25, 2008. Maidment, D.: Arc Hydro: GIS for Water Resources, ESRI, Redlands, CA, USA, 55–86, 2002. Marks, D., Dozier, J., and Frew, J.: Automated basin delineation from digital elevation data, GeoProcessing, $2$, 299–311, 1984. McDaniel, P. A. and Buol, S. W.: Manganese distributions in acid soils of the North Carolina Piedmont, Soil Sci. Soc. Am. J., 55, 152–158, 1991. McDaniel, P. A., Regan, M. P., Brooks, E., Boll, J., Bamdt, S., Falen, A., Young, S. K., and Hammel, J. E.: Linking fragipans, perched water tables, and catchment-scale hydrological processes, Catena, 73, 166–173, 2008. Noguchi, S., Tsuboyama, Y., Sidle, R. C., and Hosoda, I.: Morphological characteristics of macropores and the distribution of preferential flow paths in a forested slope segment, Soil Sci. Soc. Am. J., 63, 1413–1423, 1999. O'Callaghan, J. F. and Mark, D. M.: The extraction of drainage networks from digital elevation data, CVGIP, 28, 323–344. 1984. Orlandini, S., Moretti, G., Franchini, M., Aldighieri, B., and Testa, B.: Path-based methods for the determination of nondispersive drainage directions in grid-based digital elevation models, Water Resour. Res., 39, 1144, doi:10.1029/2002WR001639, 2003. Paik, K.: Global search algorithm for nondispersive flow path extraction, J. Geophys. Res., 113, F04001, doi:10.1029/2007JF000964, 2008. Palkovics, W. E. and Petersen, G. W.: Contribution of lateral soil-water movement above a fragipan to streamflow, Soil Sci. Soc. Am. J., 41, 394–400, 1977. Patrick, W. H. and Henderson, R. E.: Reduction and reoxidation cycles of Manganese and Iron in flooded soil and in water solution, Soil Sci. Soc. Am. J., 45, 855–859, 1981. Perillo, C. A., Gupta, S. C., and Moncrief, J. F.: Prevalence and initiation of preferential flow paths in a sandy loam with argillic horizon, Geoderma, 89, 307–331, 1999. Quinn, P. F., Beven, K. J., Chevallier, P., and Planchon, O.: The prediction of hillslope flowpaths for distributed modelling using digital terrain models, Hydrol. Proc., 5, 59–80, 1991. Reedy, R. C. and Scanlon, B. R.: Soil water content monitoring using electromagnetic induction, J. Geotech. Geoenviron. Eng., 129, 1028–1039, 2003. Rhoades, J. D., Raats, P. A. C., and Prather, R. S.: Effects of liquid-phase electrical conductivity water content and surface conductivity on bulk soil electrical conductivity, Soil Sci. Soc. Am. J., 40, 651–665, 1976. Sander, T. and Gerke, H. H.: Preferential flow patterns in paddy fields using a dye tracer, Vadose Zone. J., 6, 105–115, 2007. Seibert, J. and McGlynn, B. L.: A new triangular multiple flow direction algorithm for computing upslope areas from gridded digital elevation models, Water Resour. Res., 43, W04501, doi:10.1029/2006WR005128, 2007. Schäuble, H., Marinoni, O., and Hinderer, M.: A GIS-based method to calculate flow accumulation by considering dams and their specific operation time, Comput. Geosci., 34, 635–646, 2008. Sherlock, M. D. and McDonnell, J. J.: A new too for hillslope hydrologists: spatial distributed groundwater level and soilwater content measured using electromagnetic induction, Hydrol. Proc., 17, 1965–1977, 2003. Sidle, R. C., Noguchi, S., Tsuboyama, Y., and Laursen, K.: A conceptual model of preferential flow systems in forested hillslopes: evidence of self-organization, Hydrol. Proc., 15, 1675–1962, 2001. Soil Survey Division Staff.: Soil Survey Manual, U.S. Department of Agriculture Handbook No. 18, US Government Printing Office, Washington, DC, USA, 1993. Sørensen, R., Zinko, U., and Seibert, J.: On the calculation of the topographic wetness index: evaluation of different methods based on field observations, Hydrol. Earth Syst. Sci, 10, 101–112, 2006. Tarboton, D. G.: A new method for the determination of flow directions and upslope areas in grid digital elevation models, Water Resour. Res., 33, 309–319, 1997. Thompson, J. A., Pena-Yewtukhiw, E. M., and Grove, J. H.: Soil-landscape modelling across a physiographic region: Topographic patterns and model transportability, Geoderma, 133, 57–70, 2006. Tsukamoto Y. and Ohta, T.: Runoff process on a steep forested slope, J. Hydrol., 102, 165–178, 1988. Vachaud, G., De Silans Passerat, A., Balabanis, P., and Vauclin, M.: Temporal stability of spatial measured soil water probability density function, Soil Sci. Soc. Am. J., 49, 822–827, 1985. Van Genuchten, M. T.: A closed form equation for predicting the hydraulic conductivity of unsaturated soils, Soil Sci. Soc. Am. J., 44, 892–898, 1980. Vogt, J. V., Colombo, R., and Bertolo, F.: Deriving drainage networks and catchment boundaries: a new methodology combining digital elevation data and environmental characteristics, Geomorphology, 53, 281–298, 2003. Walker, C. and Lin, H. S.: Soil property changes after four decades of wastewater irrigation: a landscape perspective, Catena, 73, 63–74, 2008. Wu, S., Li, J., and Huang, G. H.: A study on DEM-derived primary topographic attributes for hydrologic applications: sensitivity to elevation data resolution, Appl. Geogr., 28, 210–223,~2008. Yaalon, D. H., Jungreis, C., and Koyumjisky, H.: Distribution and reorganization of manganese in three catenas of Mediterranean soils, Geoderma, 7, 71–78, 1972. Zhu, Q. and Lin, H. S.: Combining sample size, spatial structure, and auxiliary variables to determine optimal kriging in contrasting landscapes, Ecol. Model., in review, 2009. Zhu, Q., Lin, H. S., and Doolittle, J.: Repeated electromagnetic induction surveys for understanding landscape soil and water dynamics, Soil Sci. Soc. Am. J. in review, 2009. Zinko, U., Seibert, J., Dynesius, M., and Nilsson, C.: Plant species numbers predicted by topography-based groundwater flow index, Ecosystems, 8, 430–441, 2005.