Integration of 2D lateral groundwater flow into the variable infiltration capacity (VIC) model and effects on simulated fluxes for different grid resolutions and aquifer diffusivities

Better representations of groundwater processes need to be incorporated into large-scale hydrological models to improve simulations of regional- to global-scale hydrology and climate, as well as understanding of feedbacks between the human and natural systems. We incorporated a 2D groundwater flow model into the variable infiltration capacity (VIC) hydrological model code to address its lack of a lateral groundwater flow component. The water table was coupled with the variably saturated VIC soil column allowing bi-directional exchange of water between the aquifer and the soil. We then investigated how variations in aquifer properties and grid resolution affect modelled evapotranspiration (ET), runoff and groundwater recharge. We simulated nine idealised, homogenous aquifers with different combinations of transmissivity, storage coefficient, and three grid resolutions. The magnitude of cell ET, runoff, and recharge significantly depends on water table depth. In turn, the distribution of water table depths varied significantly as grid resolution increased from 1° to 0.05° for the medium and high transmissivity systems, resulting in changes of model-average fluxes of up to 12.3% of mean rainfall. For the low transmissivity aquifer, increasing the grid resolution has a minimal effect as lateral groundwater flow is low, and the VIC grid cells behave as vertical columns. The inclusion of the 2D groundwater model in VIC will enable the future representation of irrigation from groundwater pumping, and the feedbacks between groundwater use and the hydrological cycle.