Influence of precipitation infiltration recharge on hydrological processes of the karst aquifer system and adjacent river

Abstract

Understanding the hydrological processes of the karst aquifer system (KAS) and adjacent river can further accurately evaluate karst water storage. However, the mechanism of the KAS and adjacent river hydrological processes under precipitation infiltration recharge conditions, the micro-flow behavior of groundwater in the KAS, and changes in media saturation remain unclear. This study coupled the KAS and adjacent river through the Volume of Fluid (VOF), simulated the flow behavior of groundwater in different medium of the KAS and coupled the free flow process of the adjacent river by using the multiphase Darcy-Brinkman-Stokes (DBS) equation, applied the van Genuchten-Mualem and van Genuchten-modified Mualem equations to characterize the seepage processes in soil. The results indicate that the saturation of different porous media controlled by recharge intensity directly influences the flow behavior and flow paths of the porous media. The main groundwater flow path under low precipitation intensity (b=3) is primarily the lower porous medium, while under moderate precipitation intensity (3<b≤8), the main flow path of groundwater is the karst conduit. During heavy precipitation (8<b≤12), due to the limited capacity of karst conduit, they cannot bear the large amount of groundwater supplied from upper porous medium. The peak conduit saturation in laminar flow is 0.3 lower than that in turbulent flows. Under laminar flow conditions, the conduit drainage leads to a rapid decrease in saturation, making it easier for porous media to recharge to the conduit and discharge towards the adjacent river. The turbulent state conduit maintains a high saturation recession, but the discharge is consistently lower than laminar flow. The variation process line of average flow velocity and interfacial flow velocity are obtained using the DBS method and N-S simulation method respectively. There is a certain discrepancy between the simulation results of the two methods since the N-S simulation scheme dose not consider the water retention characteristics of the porous medium.