The main outlets of karst systems are springs, the hydrographs of which are largely affected by flow processes in the unsaturated zone. These processes differ between the epikarst and transmission zone on the one hand and the matrix and conduit on the other hand. However, numerical models rarely consider the unsaturated zone, let alone distinguishing its subsystems. Likewise, few models represent conduits through a second medium, and even fewer do this explicitly with discrete features. This paper focuses on the interest of hybrid models that take into account both unsaturated subsystems and discrete conduits to simulate the reservoir-scale response, especially the outlet hydrograph. In a synthetic karst aquifer model, we performed simulations for several parameter sets and showed the ability of hybrid models to simulate the overall response of complex karst aquifers. Varying parameters affect the pathway distribution and transit times, which results in a large variety of hydrograph shapes. We propose a classification of hydrographs and selected characteristics, which proves useful for analysing the results. The relationships between model parameters and hydrograph characteristics are not all linear; some of them have local extrema or threshold limits. The numerous simulations help to assess the sensitivity of hydrograph characteristics to the different parameters and, conversely, to identify the key parameters which can be manipulated to enhance the modelling of field cases.
Combining geostatistics and simulations of flow and transport to characterize contamination within the unsaturated zone
