Modeling the Matrix-Conduit Exchanges in Both the Epikarst and the Transmission Zone of Karst Systems

Abstract

Usual conceptual models of karst hydrodynamics highlight the important role of unsaturated subsystems in recharge repartition. However, few of them have been compared with scarce suitable physically-based numerical models. Hybrid models that couple single continuum medium with discrete features promise an improved consideration of karst specificities. Here we evaluate their capability to properly reproduce interactions between a vertical conduit and the surrounding unsaturated matrix. We simulate the response of such a configuration to a single recharge event for various sets of parameters. We show the ability of hybrid models to reproduce the most significant behaviors described in the literature, i.e., transient storage and distribution of recharge, flow concentration towards conduits in the epikarst, and matrix-conduit exchanges varying in time and space. In addition to the explicit conduits, simulating variably saturated flows with the Richards equation and distinguishing the epikarst and the transmission zone are key elements to reproduce most processes. The contrasts between subsystems necessary to observe desired behaviors have been quantified. They are reinforced by the varying matrix saturation that causes realistic competition between matrix and explicit conduits. The study also highlights the need to deepen knowledge of the scaled medium properties we need to know to apply such models to actual cases.