Quantification of seepage in a multi-layered disconnected river-aquifer system

Abstract

Quantification of seepage in disconnected river-aquifer systems is significant for local water management and groundwater pollution control, especially in areas with water shortage or contamination. The vadose zone under riverbeds usually exhibits a multi-layered structure, particularly when paved with low permeability liners. To evaluate the impact of engineering solutions to seepage under such conditions, we proposed an approach by combining GIS and the minimum flux in saturation layer (MFSL) method. MFSL can calculate the stable seepage rate by assessing the dominant low permeability layers (including but not limited to the liners) in multi-layered disconnected river-aquifer systems. We used MFSL to calculate local seepage rate, and used GIS to extend the results to a regional scale. The reliability of MFSL is discussed by comparing the results with the double ring infiltration test, the numerical simulation by HYDRUS, and the methods of stream package in MODFLOW, including its improved form. A case study was conducted in the Yongding River with river-aquifer seepage calculated under various conditions, including different river water levels (i.e., under the designated water level, drought stage level, flood stage level and flood inundation level) and with/without low permeability liners (i.e., ecological membranes or geomembrane). Results showed that low permeability liners could greatly reduce the seepage rate. However, if an unlined inundation area exists, the seepage rate may increase greatly. The results indicated that the proposed method was reliable and convenient for calculating long-term, large area seepage in disconnected river-aquifer systems especially those paved with liners.

Keywords

infiltration / low permeability liner / multi-layered porous media / river-aquifer seepage

References

[1] Swain E D, Howie  B, Dixon J . Description and field analysis of a coupled ground-water/surface-water flow model (MODFLOW/BRANCH) with modifications for structures and wetlands in southern Dade County. Tallahassee: U.S. Geological Survey, 1996

[2] Wilson L G, De Cook  K J. Field observations on changes in the subsurface water regime during influent seepage in the Santa Cruz River. Water Resources Research, 1968, 4(6): 1219–1234

[3] Jonnalagadda S B ,  Mhere G . Water quality of the Odzi River in the eastern highlands of Zimbabwe. Water Research, 2001, 35(10): 2371–2376

[4] Brunner P, Cook  P G, Simmons  C T. Disconnected surface water and groundwater: from theory to practice. Ground Water, 2011, 49(4): 460–467

[5] McDonald M G, Harbaugh  A W. A modular three-dimensional finite difference groundwater model. United States Geological Survey, Openfile Report No. 6, 1996

[6] Fok Y S. One-Dimensional infiltration into layered soils. Journal of the Irrigation and Drainage Division, 1970, 96(2): 121–129

[7] Bouwer H. Theory of seepage from open channels.In: Chow V T, eds. Advances in hydroscience. Vol. 5. New York: New York Academic Press, 1969, 121–172

[8] Rovey C E K . Numerical model of flow in stream aquifer system. Fort Collins: Colorado State University,1975

[9] Osman Y Z, Bruen  M P. Modelling stream–aquifer seepage in an alluvial aquifer: an improved loosing-stream package for MODFLOW. Journal of Hydrology, 2002, 264(1–4): 69–86

[10] Crosbie R S, Taylor  A R, Davis  A C, Lamontagne  S, Munday T . Evaluation of infiltration from losing-disconnected rivers using a geophysical characterisation of the riverbed and a simplified infiltration model. Journal of Hydrology, 2014, 508(1): 102–113

[11] Wang C, Mao  X, Zhao B . Experiments and simulation on infiltration into layered soil column with sand interlyaer under ponding condition. Transactions of the Chinese Society of Agricultural Engineering, 2010, 26(11): 61–67 (in Chinese)

[12] Šimunek J, Šejna  M, Saito H ,  Sakai M ,  Van-Genuchten T . The HYDRUS-1D software package for simulating the one-dimensional movement of water, heat, and multiple solutes in variably-saturated media. California: Department of environmental sciences, University of California Riverside, 2009

[13] Šimůnek J ,  van Genuchten M T ,  Šejna M . Development and applications of the HYDRUS and STANMOD software packages and related codes. Vadose Zone Journal, 2008, 7(2): 587–600

[14] Trefry M G, Muffels  C. Feflow: a finite-element ground water flow and transport modeling tool. Ground Water, 2007, 45(5): 525–528

[15] Wang C, Mao  X, Hatano R . Modeling ponded infiltration in fine textured soils with coarse interlayer. Soil Science Society of America Journal, 2014, 78(3): 745–753

[16] Mao X M, Shang  S H. Method of minimum flux in saturation layer for calculating stable infiltration through layered soil. Journal of Hydraulic Engineering, 2010, 7(41): 810–817 (in Chinese)

[17] Beijing Water Authority. Promote the construction of Yongding River ecological corridor and establish the foundation of water bank economic development in Mentougou District. Beijing Water, 2010, (3): 74 (In Chinese)

[18] Jiang B, Wong  C P, Lu  F, Ouyang Z Y ,  Wang Y J . Drivers of drying on the Yongding River in Beijing. Journal of Hydrology, 2014, 519(1): 69–79

[19] Hu H Z. Study of reactive transport of reclaimed water through vadose zone and its impact on the groundwater environment with prevention measures of seepage in lakes and rivers. Dissertation for the Doctoral Degree. Beijing: China Agricultural University, 2014 (in Chinese)Lin J, Yang Q, Yang Q F, Liu J, Chen Z, Dou W, Jiang Y, Liu Z, Zhao W, Ni T. Impacts of reclaimed water reuse in the typical river channel of the Yongding River on the groundwater and model establishment. Research report. Beijing: Beijing institute of hydrogeology and engineering geology. 2011 (in Chinese)

[20] Huang C Y, Xu  J M. Soil Science. Beijing: China Agriculture Press, 2010 (in Chinese)

[21] Hillel D, Baker  R S. A descriptive theory of fingering during infiltration into layered soils. Soil Science, 1988, 146(146): 207–217

[22] van Genuchten M T . A closed-form equation for predicting the hydraulic conductivity of unsaturated soils. Soil Society of America Journal, 1980, 44(44): 892–898

[23] Mualem Y. A new model for predicting the hydraulic conductivity of unsaturated porous media. Water Resources Research, 1976, 12(3): 513–522

[24] van Genuchten M T . A closed-form equation for predicting the hydraulic conductivity of unsaturated soils 1. Soil Science Society of America Journal, 1980, 44(5): 892–898

[25] Hao Q C, Shao  J L, Cui  Y L, Xie  Z H. Applicability of artificial recharge of groundwater in the Yongding River alluvial fan in Beijing through numerical simulation. Journal of Earth Science, 2014, 25(3): 575–586

[26] Pan L, Warrick  A W, Wierenga  P J. Finite element methods for modeling water flow in variably saturated porous media: numerical oscillation and mass-distributed schemes. Water Resources Research, 1996, 32(6): 1883–1889

Acknowledgements

This work was supported by the National Natural Science Foundation of China (51379207, 51321001) and Beijing Municipal Science and Technology Project (D090409004009004).

Compliance with ethics guidelines

Jiang Li, Haizhu Hu, Xiaomin Mao, and Yi Liu declare that they have no conflict of interest or financial conflicts to disclose.

This article does not contain any studies with human or animal subjects performed by any of the authors.

RIGHTS & PERMISSIONS

The Author(s) 2017. Published by Higher Education Press. This is an open access article under the CC BY license (http://creativecommons.org/licenses/by/4.0)