Data-model comparisons of isotopic and hydrological variabilities of the karstic vadose zone above Villars Cave, SW-France based on 20 years' monitoring record

Because of the great heterogeneity of the limestone fissures, conduits, and reservoirs in karst areas, the understanding of the infiltration in caves and karst terrains is complicated. For example, the time delays observed between rainfall and the various flows into caves, such as shaft flow and seepage flow, exhibit significant variability. To better understand the infiltration above a typical shallow cave, Villars Cave, South-west France, we used both a very long monitoring (>20 years) isotopic data set of rainfall and dripping rates at two different cave gallery levels (14 m and 36 m below the surface, respectively) and a conceptual model (KarstFor) that mix water flows and isotopes. In order to determine the most appropriate reservoir and flow settings, we used one time-lapse electrical resistivity tomography (ERT) image made above the cave, which were incorporated into the KarstFor model. The modeled variations in storage water and infiltration flows are in a good agreement with the observed dripping discharge rates, indicating seasonal variability and various time delays between water excess ((Rainfall-Evapotranspiration) > 0) and drip rates in the different cave levels. The fact that the model strikingly simulates the small, but constant, isotopic difference of drip water delta O-18 (delta O-18(d)) values between the upper and the lower galleries of Villars Cave, not only reinforces reliability in the model itself but also reveals the impact of multiple infiltration routes and karst reservoir dynamics on the delta O-18(d). On an inter-annual scale, the observed delta O-18(d) values follow the variations of the water excess and of the cave temperature, which, in the period analyzed, show a marked trend, possibly induced by the climate variability patterns above the cave. Overall, these results are not only important to better understand the infiltration in karstic zones, but also to better interpret speleothem isotopic records widely used in paleo-climate studies.