Capturing soil-water and groundwater interactions with an iterative feedback coupling scheme: new HYDRUS package for MODFLOW

Abstract

Abstract. Accurately capturing the complex soil-water and groundwater interactions is vital for describing the coupling between subsurface–surface–atmospheric systems in regional-scale models. The nonlinearity of Richards' equation (RE) for water flow, however, introduces numerical complexity to large unsaturated–saturated modeling systems. An alternative is to use quasi-3-D methods with a feedback coupling scheme to practically join sub-models with different properties, such as governing equations, numerical scales, and dimensionalities. In this work, to reduce the nonlinearity in the coupling system, two different forms of RE are switched according to the soil-water content at each numerical node. A rigorous multi-scale water balance analysis is carried out at the phreatic interface to link the soil-water and groundwater models at separated spatial and temporal scales. For problems with dynamic groundwater flow, the nontrivial coupling errors introduced by the saturated lateral fluxes are minimized with a moving-boundary approach. It is shown that the developed iterative feedback coupling scheme results in significant error reduction and is numerically efficient for capturing drastic flow interactions at the water table, especially with dynamic local groundwater flow. The coupling scheme is developed into a new HYDRUS package for MODFLOW, which is applicable to regional-scale problems.