An Explicit Scheme to Represent the Bidirectional Hydrologic Exchanges Between the Vadose Zone, Phreatic Aquifer, and River

Understanding groundwater storage variations as a result of effects from both the vadose zone and the river is of critical importance in the hydraulically connected surface-subsurface system. In this study, we present a novel Bidirectional Exchange Scheme in Surface and Subsurface (BE3S) that represents bidirectional exchanges between the vadose zone, phreatic aquifer, and river. The approach enables explicit representations of each flow regime while conserving mass, and successfully yields solutions of the coupled system for multiple temporal resolutions. We test the scheme by comparing the BE3S-derived outputs against other models and corresponding observational data. We apply the scheme to simulate hydrologic states in a reach of the Brazos River in Southeast Texas, such as soil moisture content, groundwater level, and river stage as well as net subsurface discharge fluxes. Good agreements between the simulated and observed data for all the components show the suitability of the proposed scheme in modeling the bidirectional flows and exchanges. We also assess how the bidirectional hydrologic exchanges are affected by adjacent flow domains. We find that vertical percolation is significantly affected by unsaturated soil thickness, resulting in the spatial variability of vertical percolation across sloping topography. Hilltop-to-valley convergent groundwater flow is also found to impede vertical percolation in river valleys due to shallow water table while facilitating percolation in the hilltops. The capability of the presented scheme that accounts for the topographically driven lateral groundwater flow and drainage dynamics provides the potential to enhance the representation of the surface-subsurface system in Earth System Models (ESMs).