Soil Property and Subsurface Heterogeneity Control on Groundwater Recharge of Vadose Zone Injection Wells

Abstract

Vadose zone well (VZW) injection is an effective method of managed aquifer recharge (MAR), and it plays an important role in semiarid and arid regions. Accurate estimation of the recharge of VZW injection remains a challenge in the field of hydrology due to the nonlinear nature of the process in the vadose zone. To improve knowledge of VZW injection, a finite-element numerical model was developed by COMSOL Multiphysics to evaluate the soil property and subsurface heterogeneity control on the recharge of VZWs. Several numerical experiments were conducted to estimate the characteristic arrival time and cumulative recharge volume of injected water for different subsurface conditions. Simulation results indicate that soil properties have a profound influence on the recharge of VZWs and saturated hydraulic conductivity plays a more important role than the Van Genuchten –Mualem soil constitutive model parameters α and n in the influences on the recharge. Subsurface heterogeneity has large impacts on the recharge of VZWs as well. The presence of a low permeable layer in the aquifer hinders the infiltration of injected water and reduces recharge rate Qwt and cumulative recharge volume Vwt. These influences are affected by depth Dc, lateral extension Rc, thickness Bc, and hydraulic conductivity Ks of the low permeable layer. Whether the well screen cuts through the low permeable layer is a factor that influences the effects of these parameters in two forms. When the well screen cuts through the low permeable layer, Qwt and Vwt decrease substantially as Dc increases, and the change of Bc has a greater impact on recharge than Ks; when the low permeable layer is below the well screen, the change of Dc only has slight influence on recharge, and the cases with the same hydraulic conductance value of Ks/Bc have similar Qwt and Vwt. For the low permeable lens, when the well screen cuts through it, the increase of Rc leads to slight decrease of Qwt and Vwt, which are more obvious in the later stage of VZW injection. When the low permeable lens is below the well screen, the increase of Rc makes Qwt decrease in the early stage of VZW injection, and then Qwt for different values of Rc reach similar asymptotic values in the later stage. Understanding the influences of soil properties and subsurface heterogeneity on the recharge of VZWs provides physically based guidance for the design and management of VZWs.