Journal of the Croatian Geological Survey and the Croatian Geological Society heterogeneous hydraulic properties manifested by rapid groundwater flows concentrated in conduits and disolutionally widened fractures, usually sparsely distributed, while the storage of groundwater is mainly related to surrounding fractured rock matrix and surficial epikarst (or subcutaneous) zone (GOLDS-CHEIDER et al., 2007). Because of its complexity, karst aquifer research necessitates an interdisciplinary approach (DREYBRODT, 1988; COVING-TON et al., 2011). In addition to conventional methods, the use of environmental stable isotope tracers supplemented by statistical analysis, contributes to a better understanding of karst hydrogeology (CRISS et al., 2007; PALMER, 2010; GRIMMEISEN et al., 2017). Stable isotopes of water (2 H and 18 O) are analysed for several reasons, e.g. to explore recharge processes and reservoir mixing (LONG & PUTNAM, 2004; MALOSZEWSKI et al., 2002), to explore vadose water mixing, residence times and storage in non-conduit parts of karst systems (PALMER, 2010; MO-ERMAN et al., 2014; DÓŠA et al., 2011), to estimate mean recharge area elevations (BARBIERI et al., 2005; CHARIDEH & RAHMAN, 2007), etc. Detailed knowledge of the various processes within a karst hydrological system is crucial for the sustainable management of water resources, their allocation, and protection against pollution (MANCE et al., 2018 a, b). The black-box approach is often used