State Key Laboratory of Water Resources and Hydropower Engineering Science, Wuhan University, Wuhan 430072, P. R. China

A coupled unsaturated-saturated water flow numerical model was developed. The water flow in the unsaturated zone is considered the one-dimensional vertical flow, which changes in the horizontal direction according to the groundwater table and the atmospheric boundary conditions. The groundwater flow is treated as the three-dimensional water flow. The recharge flux to groundwater from soil water is considered the bottom flux for the numerical simulation in the unsaturated zone, and the upper flux for the groundwater simulation. It connects and unites the two separated water flow systems. The soil water equation is solved based on the assumed groundwater table and the subsequent predicted recharge flux. Then, the groundwater equation is solved with the predicted recharge flux as the upper boundary condition. Iteration continues until the discrepancy between the assumed and calculated groundwater nodal heads have a certain accuracy. Illustrative examples with different water flow scenarios regarding the Dirichlet boundary condition, the Neumann boundary condition, the atmospheric boundary condition, and the source or sink term were calculated by the coupled model. The results are compared with those of other models, including Hydrus-1D, SWMS-2D, and FEFLOW, which demonstrate that the coupled model is effective and accurate and can significantly reduce the computational time for the large number of nodes in saturated-unsaturated water flow simulation.

The work was supported by the National Basic Research Program of China (Grant No. 2010CB428802), the Scholarship Award for Excellent Doctoral Students granted by Ministry of Education, the Fundamental Research Funds for the Central Universities, and the Ph. D. Candidates’ Self-research Program of Wuhan University in 2008.

Carrera, J., and Medina, A. 1999. A discussion on the calibration of regional groundwater models. Feyen, J., and Wiyo, K., eds., Proceedings of International Workshop of EurAgEng’s Field of Interest on Soil and Water, Modelling of Transport Processes in Soils at Various Scales in Time and Space, 629-640.  Wageningen: Wageningen Pers.
Chen, Z. Q., Govindaraju, R. S., and Kavvas, M. L. 1994. Spatial averaging of unsaturated flow equations under infiltration conditions over areally heterogeneous fields: 1. Development of models. Water Resources Research, 30(2), 523-534. [doi:10.1029/93WR02885]
Diersch, H. J. G., and Kolditz, O. 1998. Coupled groundwater flow and transport: 2. Thermohaline and 3D convection systems. Advances in Water Resources, 21(5), 401-425. [doi:10.1016/S0309-1708(97)  00003-1]
Facchi, A., Ortuani, B., Maggi, D., and Gandolfi, C. 2004. Coupled SVAT: Groundwater model for water resources simulation in irrigated alluvial plains. Environmental Modelling and Software, 19(11), 1053-1063. [doi:10.1016/j.envsoft.2003.11.008]
Havard, P. L., Prasher, S. O., Bonnell, R. B., and Madani, A. 1995. LINKFLOW, a water flow computer model for water table management, Part I: Model development. Transactions of the ASABE, 38(2), 481-488.
Lin, L., Yang, J. Z., Zhang, B., and Zhu, Y. 2010. A simplified numerical model of 3D groundwater and solute transport at large scale area. Journal of Hydrodynamics, 22(3), 319-328. [doi:10.1016/S1001-6058(09) 60061-5]
Mantoglou, A. 1992. A theoretical approach for modeling unsaturated flow in spatially variable soils: Effective flow models in finite domains and nonstationarity. Water Resources Research, 28(1), 251-267. [doi: 10.1029/91WR02232]
McDonald, M. G., and Harbaugh, A. W. 1988. A modular three-dimensional finite-difference groundwater flow model. Techniques of Water-resources Investigations. Denver: U.S. Geological Survey.
Niswonger, R. G., Prudic, D. E., and Regan, R. S. 2006. Documentation of the unsaturated-zone flow (UZF1) package for modeling unsaturated flow between the land surface and the water table with MODFLOW-2005. Techniques and Methods, 6-A19. Denver: U.S. Geological Survey.
Pikul, M. F., Street, R. L., and Remson, I. 1974. A numerical model based on coupled one-dimensional Richards and Boussinesq equations. Water Resources Research, 10(2), 295-302. [doi:10.1029/ WR010i002p00295]
Vachaud, G., and Vauclin, M. 1975. Comments on ‘A numerical model based on coupled one dimensional Richards and Boussinesq equations’ by Mary F. Pikul, Robert L. Street, and Irwin Remson. Water Resources Research, 11(3), 506-509. [doi:10.1029/WR011i003p00510]
Vogel, T., Huang, K., Zhang, R., and van Genuchten, M. Th. 1996. The Hydrus Code for Simulating One-Dimensional Water Flow, Solute Transport, and Heat Movement in Variably-satruated Media. Riverside: U.S. Salinity Laboratory, Agriculture Research Service, U.S. Department of Agriculture.
Zhu, Y., Yang, J. Z., and Tong, J. X. 2010. A simplified 3D numerical model for groundwater flow in declining aquifer-aquitard system. Journal of Sichuan University (Engineering Science Edition), 42(6), 43-50. (in Chinese)