www.geologia-croatica.hr Journal of the Croatian Geological Survey and the Croatian Geological Society 2024 | 77/2 | 145–158 | 6 Figs. | 2 Tabs. | 2 Supplements | Article history: Manuscript recieved: April 18, 2024 Revised manuscript accepted: May 27, 2024 Available online: June 21, 2024 Keywords: karst hydrogeology, major ions, isotopes, time series analysis, emerging organic contaminants, Jadro and Žrnovnica springs 1. INTRODUCTION Karst aquifers are intricate, highly dynamic and heterogene- ous systems with enlarged fractures, an often well-developed conduit network, direct surface-subsurface connection via ponors, and high hydraulic conductivity, leading to preferen- tial flow pathways and significant drinking water resources susceptible to anthropogenic contamination (PADILLA & VESPER, 2018; ŽIVANOVIĆ et al., 2022). Interpreting hydro- chemical spring responses and identifying the origins of chem- ical compounds in karst water, some of which as natural traces allow the understanding of system functioning and structure (PERRIN, 2003; MUDARRA et al., 2011; SHEIKHY NA- RANY et al., 2019), can help to pinpoint contamination sources and pathways (HILLEBRAND et al., 2014). The hy- drochemical research of the karst aquifer in question is funda- mental in characterising the behaviour and origin of emerging organic contaminants (EOCs). In 26 sampling campaigns within the Jadro and Žrnovnica Dinaric karst catchment (Cro- atia), water samples were collected for ion and stable isotope analysis. Notably, in seven of these campaigns, samples were Hydrochemical characterization of a Dinaric karst catchment in relation to emerging organic contaminants Ana Selak1, Jasmina Lukač Reberski1,*, Maja Briški1 and Lorena Selak2 1 Croatian Geological Survey, Department of Hydrogeology and Engineering Geology, Sachsova 2, 10000 Zagreb, Croatia; (*corresponding author: jlukac@hgi-cgs.hr) 2 Aarhus University, Department of Biology, Section for Microbiology, Ny Munkegade 114, 8000 Aarhus C, Denmark doi: 10.4154/gc.2024.10 Abstract The main findings of a hydrochemical investigation conducted within a typical Dinaric karst catchment located in Southern Croatia are outlined. The studied aquifer is drained by the Jadro and Žrnovnica springs, which are important for the regional and local water supplies, respectively. Presumably, there is intercatchment groundwater flow coming from the neigh- bouring Cetina River catchment. Various factors governing aquifer hydrochemistry and their interplay with emerging organic contaminants (EOCs) that were detected at different water resources in ng/L concentrations was assessed. A total of 26 sampling campaigns (Octo- ber 2019 – October 2022) were conducted at two springs, in a river and at a deep borehole, all representative of this complex hydrogeological system. Assessment of major ion con- stituents and saturation indices calculated with PHREEQC revealed the sampled water re- sources are of a Ca-HCO3 type due to the predominant weathering of the carbonate min- eral calcite. Sharp spikes observed in chemographs indicated a highly karstified system with an effective conduit network allowing rapid spring responses to precipitation events. Water resources are of good chemical status, as affirmed by anthropogenic contamination indicators, with nitrates, chlorides and sulphates all below maximum threshold values. Strong positive correlations were found between EOCs concentration, number of detected compounds, and nitrates in the Cetina River, indicating a common origin, most likely waste- water. Identification of persistent EOCs including widely used repellent N,N-diethyl-meta- toluamide (DEET) during base flow conditions and its strong positive correlation with the Ca2+ content in both the Cetina and Jadro samples, suggests potential storage in the epi- karst and aquifer matrix. This coupling of conventional hydrochemical indicators and novel markers of anthropogenic impacts, including EOCs, in vulnerable karst water resources is a crucial advancement in the assessment and management of emerging environmental and potential human health risks. Such an approach is pivotal for the sustainable protection of hydrogeologically intricate sites. specifically taken for EOC analysis. EOCs are typically en- countered in karst water resources at concentrations ranging from ng/L to µg/L, encompassing a myriad of anthropogenic organic compounds (LUKAČ REBERSKI et al., 2023). Some EOCs exert adverse effects on aquatic ecosystems and human health (CIZMAS et al., 2015; SCHRIKS et al., 2010). Numerous uncertainties persist regarding the presence and behaviour of EOCs within large hydrogeological domains (LLAMAS et al., 2022). Due to diverse physico-chemical properties and sources, coupled with the intricacy of Dinaric karst structures and (ground)water pathways, datasets and analyses necessary to comprehend EOCs behaviour can be extensive. Multivariate statistical analyses have a broad application in simplifying large datasets, facilitating characterization of interrelationships among numerous hydrogeological variables, uncovering main hydrogeological processes, and providing insights into temporal and spatial variations (MATIATOS et al., 2014; JEBREEN et al., 2018). These invaluable analytical tools were employed in this study for the reliable interpretation of intricate karst dynamics and mailto:jlukac@hgi-cgs.hr G eo lo gi a C ro at ic a 146 Geologia Croatica 77/2 hydrochemical fluctuations (CAETANO BICALHO, 2012; VASIĆ et al., 2020; ĆUK ĐUROVIĆ et al., 2022). The hydrological and hydrogeological features of the Jadro and Žrnovnica aquifer have been extensively documented in prior studies by JUKIĆ & DENIĆ-JUKIĆ (2008), KAPELJ et al. (2012), BONACCI & ANDRIĆ (2015), KADIĆ et al., (2019), and SELAK et al. (2024). Additionally, the (eco)toxicological aspects related to EOCs occurrence in this karst aquifer were previously explored by SELAK et al. (2022, 2024). This study seeks to enhance our understanding of karst aquifer behaviour under varying hydrological conditions by utilizing major ion and stable water isotopes data. The analysis of stable isotopes of oxygen δ18O and hydrogen δ2H has been effective in discerning groundwater provenance and karst aquifer recharge processes (CLARK & FRITZ, 1997). We posited that certain identified EOCs could correlate with physico-chemical parameters of water resources within the studied system. These relationships may facilitate the identification of potential contamination sources, behaviour, and transport. We pursued the objective of elucidating factors that govern the hydro- chemistry of the Jadro and Žrnovnica springs catchment. This involved monthly sampling at two aforementioned karst springs, a deep borehole, and the adjacent Cetina River from which presumably groundwater inflows into the studied catchment (FRITZ, 1979; JUKIĆ & DENIĆ-JUKIĆ, 2015). 2. MATERIALS AND METHODS 2.1. Study area The immediate catchment of the Jadro and Žrnovnica springs in Southern Dalmatia (Croatia) extends to an area of around 500 km2 of heterogeneous relief defined by interchanging mountains and karst poljes. Both Jadro (35 m a.s.l.) and Žrnov- nica (spring zone 78 to 90 m a.s.l.) emerge at the foot of Mosor Mountain where an impermeable coastal flysch belt meets with the permeable carbonate rocks of the hinterland (for detailed hydrogeological map readers are referred to SELAK et al., 2024). Typical karst, well-expressed geomorphology has a predominant NW-SE Dinaric orientation. The catchment mainly consists of karstified carbonate rocks of Mesozoic and partly Eocene age, characterized by their high permeability, lack of soil cover and consequently absence of surface water- courses (KAPELJ et al., 2012). Jadro with a mean discharge of 9.4 m3/s (data period 2011-2022) is used for the regional water supply of Split and its wider area. Žrnovnica with mean discharge of 1.7 m3/s (data period 2011-2022) supplies the nearby settlement and local agricultural areas. Previous studies denote overlap of the springs’ catchments and ground- water exchange between them (JUKIĆ & DENIĆ-JUKIĆ, 2008; BONACCI & ROJE-BONACCI, 1997; BONACCI & ANDRIĆ, 2015; KADIĆ et al., 2017). Moreover, a tracer test undertaken at the Grabov mlin ponor proved the presence of intercatchment groundwater flows coming from the adjacent Cetina River catchment (Fig. 1) (GEOTEHNIKA, 1975; FRITZ, 1979), thereby expanding the overall catchment area. The diverse relief affects the local climate, with a preva- lence of temperate humid Köppen climate type Cfa, while the mountainous areas of Kozjak and Mosor exhibit a temperate humid climate type Cfb (ŠEGOTA & FILIPČIĆ, 2003). The mean daily air temperature, recorded at Split Marjan meteor- ological station at 122 m a.s.l. (data period 2009-2022), is 17.4°C. Annual precipitation averages 1304 mm (data period 2009-2022) with non-uniform spatial distribution within the catchment. Figure 1. Location of the study area (World Karst Aquifer Map by BGR et al. (2017)). G eologia C roatica 147Selak et al.: Hydrochemical characterization of a Dinaric karst catchment in relation to emerging organic contaminants Both springs respond rapidly (within 24 hours) to rainfall, indicating the presence of a substantial rapid flow component (KADIĆ et al., 2019) primarily generated by rainfall from the Dugopolje area (JUKIĆ et al., 2022). Consequently, rapid al- terations in spring hydrochemistry, accompanied by prompt arrival of potential anthropogenic inputs can be anticipated. 2.2. Sampling and analysis of physico-chemical and hydrochemical parameters In 26 sampling campaigns between October 2019 and October 2022, physico-chemical parameters were measured in situ us- ing a WTW multi-parameter probe at the Jadro and Žrnovnica springs, Cetina River and Gizdavac borehole. Only Jadro was monitored in March 2019. The Covid-19 pandemic disrupted regular monthly surveys in 2019 and 2020. Measurements and sampling were undertaken directly at the springs, while a bucket of water was collected from the midstream of the Cetina River. Before conducting sampling and measurements in the bucket, three volumes of groundwater were pumped from the borehole Gizdavac to ensure a representative sample of the aquifer. Dur- ing our research, among the three deep boreholes used for groundwater monitoring, only the Gizdavac borehole featured a functional pump (installed at the depth of 266 metres). Alkalinity was determined by volumetric titration using 1.6N H2SO4 until reaching a pH 4.5, employing HACH digital titrator Model 16900 with bromocresol green-methyl red as an indicator. The water temperature and electrical conductivity data in the Jadro and Žrnovnica springs, as well as the Cetina River (upstream of Trilj City at Vedrine hydrological station), were recorded using HOBO data loggers U24 from Onset Computer Corporation (Bourne, MA, USA). The data was captured with hourly readings, commencing on 10th March 2021 for the Jadro and Žrnovnica springs, and on 11th March 2021 for the Cetina River. A total of 84 water samples were collected, including 26 from each of the Cetina River, Jadro and Žrnovnica springs, and an additional 6 samples from the Gizdavac borehole. Sam- ples were collected in 200 ml and 100 ml polyethylene bottles for analysis of major anions, cations, and stable water isotopes. For conservation purposes, samples were transported to the laboratory at 4°C and were analysed within 48h at the Hydro- chemical laboratory of the Croatian Geological Survey. Prin- cipal ion composition was determined using ion chromatogra- phy on a DIONEX ICS-6000 DP. Water quality analyses were assessed for cation-anion balance by calculating the relative deviation from charge balance (Δmeq = 100 x (Σmeq+ − Σmeq−)/(Σmeq+ + Σmeq−) < ±5%) (DOMENICO & SCHWARTZ, 1990). Stable water isotopes δ18O and δ2H were measured using isotope-ratio mass spectrometry with a Pic- arro L2130i (Santa Clara, USA) employing Cavity Ring-Down Spectroscopy (CRDS) technology. Isotope measurements were crosschecked against Picarro’s standards, which are pe- riodically verified against International Atomic Energy Agency (IAEA) standards: Vienna Standard Mean Ocean Wa- ter 2 (VSMOW2) and Standard Light Antarctic Precipitation 2 (SLAP2). Measurement precision was ± 0.3 ‰ for δ18Ο and ± 1 ‰ for δ2H. Deuterium excess (d-excess), defined as d = δ2H - 8 x δ18O (DANSGAARD, 1964), was used as an indica- tor of precipitation air mass origin and non-equilibrium con- ditions during evaporation. In order to take into account the infiltration of precipitation into the system, mean isotopic pre- cipitation values were weighted with the amount of precipita- tion using the following equation (YURTSEVER & GAT, 1981): d d w i n i i i P P = × ∑= ∑ 1 (1) where Pi is the monthly precipitation amount (mm), δi is the isotopic composition (‰), and n is the number of months. The Croatian Meteorological and Hydrological Service provided discharge data for the Jadro and Žrnovnica springs for the period 2011-2022 and precipitation data for 2009-2022. Groundwater temperature and electrical conductivity data from the Gizdavac borehole logger (period 2010-2022) were provided by Croatian Waters. They also provided discharge data for the Cetina River, measured at the Trilj žičara station. A Piper plot was used for classifying the predominant water types for the springs, groundwater, and surface water. This plot characterizes water based on major cation and anion content in meq/L, enabling analysis of the chemical composition for each water source (PIPER, 1944; APPELO & POSTMA, 2005). Saturation indices (SI) were calculated to determine the chemical equilibrium between minerals and groundwater using PHREEQC v. 3.7.3.15968 (PARKHURST & APPELO, 1999; USGS, 2021). When SI values are negative, minerals will be dissolved, since the water is undersaturated with respect to them. Conversely, positive SI values indicate oversaturation and mineral precipitation. 2.3. Multivariate statistical analysis Statistical analyses were performed in R (v. 1.4.1717; R CORE TEAM, 2021) using packages stats (R CORE TEAM, 2023), dunn.test (DINNO, 2017), corrplot (WEI & SIMKO, 2021), vegan (OKSANEN et al., 2022), and Hmisc (HARELL, 2023). The Shapiro-Wilk test was performed to investigate whether the hydrogeochemical data is normally distributed (p-value > 0.05). A nonparametric Kruskal-Wallis test with a post hoc Dunn test with Bonferroni correction (HOLLANDER & WOLFE, 1973; OGLE et al., 2023) was used to test for a statistically significant difference between the hydrochemical parameters from different sampling locations. To gain insights into the relationships between the hydrochemical parameters in each sample, a correlation matrix was performed using Spearman’s rank correlation coefficient (SPEARMAN, 1904) on a normalized dataset. Furthermore, the total sum of EOCs and their detection rate (number of detected compounds per sampling location) were correlated with major cations and anions. 4. RESULTS AND DISCUSSION 4.1. Hydrochemical composition of the springs, river and groundwater Temporal variability in water quality inherent to the karstic system of Jadro and Žrnovnica springs was comprehensively monitored for three years. This facilitated the inclusion of the G eo lo gi a C ro at ic a 148 Geologia Croatica 77/2 entire spectrum of hydrological conditions and detailed hydro- chemical characterization of the Dinaric karst aquifer during both dry and wet seasons. Table 1 summarizes the descriptive statistics of eleven physico-chemical parameters (T, pH, EC, Na+, K+, Mg2+, Ca2+, Cl-, NO3 -, SO4 2-, and HCO3 -) ascertained across 26 distinct sampling campaigns and different hydro- logical conditions. Consistent with findings for sizable karst aquifers (PULIDO-BOSCH, 2021), water temperature in both the springs and groundwater exhibited limited variability over the observed period. As expected, the Cetina River displayed a broader temperature range with pronounced daily oscillations due to direct contact with the atmosphere (Table 1). All water constituents in the investigated catchment exhibited a slightly alkaline pH typical of carbonate aquifers (MATIĆ et al., 2012; PATEKAR et al., 2022; MALDINI et al., 2023). Given the relatively low anthropogenic interference in the studied catch- ment (LOBOREC et al., 2015), we postulate that the electrical conductivity oscillations are predominantly governed by the dissolution of karst aquifer rocks. Comparatively, groundwater exhibited somewhat higher mean EC values, ref lecting prolonged residence periods. The Piper plot (Fig. 2) revealed a Ca-HCO3 hydrogeo che- mical facies present in all the examined samples. This facies is indicative of prevailing carbonate weathering processes, aligning with the characteristic hydrochemical profile of Dinaric karst water resources, as substantiated in prior studies (MATIĆ et al., 2012; TERZIĆ et al., 2014; FILIPOVIĆ et al., 2023; MALDINI et al., 2023). In contrast to Ca2+, Mg2+ exhibits relatively low abundance (Table 1, Fig. 1), suggesting the predominance of limestone deposits over dolomite. This observation aligns with the lithological analysis of the Table 1. Main statistical descriptors of the physico-chemical parameters and discharge (m3/s) observed within the study area. Water temperature (°C) and electrical conductivity (EC, µS/cm) data encompass the period from October 2019 to October 2022 for the Jadro, and Žrnovnica springs, and the Cetina river, and from 9.11.2010.-31.10.2022 for the Gizdavac borehole. Commencing in March 2021, logger data for both springs and the river were incorporated into the statistical analysis. pH data is given for the period from October 2019 to October 2022, while discharge data is from January 2011- December 2022 for the Jadro and Žrnovnica springs and from January 2019 – December 2022 for the Cetina river (Trilj žičara hydrological station). Major ion data (mg/L) is shown for the period October 2019 – October 2022. Site Statistics T EC pH Q Na+ K+ Mg2+ Ca2+ Cl– NO3 – SO4 2– HCO3 – Ja dr o sp rin g Min 12.47 259.62 7.25 3.73 2.38 0.15 2.21 70.71 3.69 1.50 3.17 204.50 Max 13.99 433.24 8.02 56.62 21.44 0.95 8.03 79.50 40.89 4.32 22.70 274.50 Mean 13.12 355.36 7.47 9.42 9.43 0.52 5.56 75.55 14.84 2.31 12.47 236.07 Median 13.09 351.45 7.46 6.62 6.39 0.50 5.75 75.62 9.47 2.27 11.61 236.07 Variance 0.18 1573.60 0.03 49.9 36.17 0.03 2.60 4.01 120.00 0.36 31.48 171.41 Žr no vn ic a sp rin g Min 12.32 263.50 7.35 0.27 2.96 0.22 1.49 66.52 4.10 0.88 2.87 192.76 Max 13.38 399.84 8.52 17.00 20.47 1.97 6.77 73.10 39.31 3.89 15.99 256.20 Mean 12.84 320.01 7.79 1.73 8.99 0.72 4.44 70.50 14.29 1.66 9.00 222.36 Median 12.88 312.20 7.76 0.87 6.88 0.57 4.59 71.08 10.49 1.50 8.20 224.48 Variance 0.08 1311.65 0.05 4.96 30.57 0.18 2.34 4.51 104.19 0.45 12.23 176.59 Ce tin a Ri ve r Min 6.89 263.10 7.80 11.10 1.71 0.32 3.32 55.55 2.53 0.52 9.36 176.90 Max 16.90 562.00 8.57 305.00 52.76 2.25 7.37 73.07 85.92 3.26 28.49 267.18 Mean 11.98 362.64 8.25 86.76 20.02 0.83 5.26 65.13 30.11 1.24 15.99 202.40 Median 11.85 334.10 8.29 70.90 12.06 0.80 5.15 64.93 18.80 1.05 14.43 198.86 Variance 8.26 5359.92 0.03 2771.38 304.56 0.19 0.91 10.57 713.30 0.43 23.10 316.55 G iz da va c bo re ho le Min 12.00 329.60 7.26 / 2.36 0.19 2.30 91.92 3.83 2.47 2.59 296.46 Max 12.80 501.60 7.58 / 3.38 0.78 6.43 108.89 9.93 4.22 5.65 318.42 Mean 12.44 399.82 7.38 / 2.89 0.61 5.43 97.96 5.59 3.02 3.73 306.73 Median 12.40 396.17 7.36 / 2.94 0.68 6.24 96.39 4.98 2.72 3.54 307.44 Variance 0.005 0.0003 0.01 / 0.20 0.05 2.64 40.79 5.07 0.42 1.06 65.85 Figure 2. Piper plot presentation of the hydrochemical composition of the springs, river, and groundwater samples within the study area. G eologia C roatica 149Selak et al.: Hydrochemical characterization of a Dinaric karst catchment in relation to emerging organic contaminants Gizdavac borehole, which identified limestone as the sole aquifer rock constituent in that area (KAPELJ et al., 2006). SI values calculated using PHREEQC (Table S3) indicate that most water samples are oversaturated with respect to calcite (Fig. 4). SIcalcite values ranged between -0.02 and 1.2 in the springs, ranged from 0.2 to 0.5 in groundwater, and varied between 0.4 and 1 in the Cetina River. According to WHITE (1997), a SIcalcite of -0.3 marks the kinetic threshold for slowdown in dissolution processes, while a SIcalcite value of +0.5 indicates the onset of precipitation. In contrast to the Jadro spring and groundwater samples, the Cetina River exhibited oversaturation with respect to dolomite, while over half of the Žrnovnica spring samples showed dolomite undersaturation. These carbonate minerals are the prominent constituents of the host rocks within the studied aquifer. Accordingly, river and spring samples exhibited Ca2+ prevalence for cations, with an abundance order generally observed as Ca2+ > Na+ > Mg2+ > K+. Groundwater samples demonstrated a distinctive abundance order of cations, with Ca2+ being the most prevalent, followed by Mg2+ > Na+ > K+. In groundwater and springs, a robust positive correlation (0.69 and 0.67, respectively) between Ca2+ and HCO3 - suggests their common origin from calcite dissolution. Most water samples were enriched with Ca2+ and Mg2+ relative to HCO3 - (Fig. 3a) which confirms how the dissolution of calcite and dolomite by carbonic acid is not the only source of alkaline earths. Excess Ca2+ and Mg2+ charges must be balanced by other major anions through dissolution or ion exchange processes. Figure 3b shows samples aligning with the 1:1 line, suggesting gypsum and anhydrite dissolution alongside calcite and dolomite dissolution. This observation is further supported by gypsum saturation indices averaging around -2.5 for the springs, -2.9 for groundwater, and -2.4 for the Cetina River, along with anhydrite saturation indices averaging around -3 for the springs, -3.3 for groundwater, and -2.8 for the Cetina River (Table S3). According to FRITZ’S hydrogeological study (1979), gypsum and anhydrite are expected near the surface in the northern part of the catchment along a fault delineating contact with Triassic deposits and subsequent strata, north of Muć polje. The average [Mg2+]/[Ca2+] molar ratios observed in all water samples approximated 0.1, as expected in aquifers drain- ing predominantly limestone rocks. Such ratios imply dolo- mite dissolution and calcite precipitation processes (APPELO & POSTMA, 2005). We noted a conspicuous decline in [Mg2+]/[Ca2+] during the hydrograph peaks and higher values in the recession periods. The highest values occurred at the later stages of recessions, consistent with the prevailing inter- pretation of increased mean residence time of karst ground- water during recession periods. Elevated [Mg2+]/[Ca2+] ratios can be indicative of matrix water or diffuse recharge through the epikarst, whereas lower [Mg2+]/[Ca2+] values observed in spring water are associated with conduit recharge (TORAN & REISCH, 2013). The observed lower [Mg2+]/[Ca2+] values, dis- tinct positive and negative fluctuations in chemographs (Fig. 4), and rapid hydrodynamic responses corroborate a highly karstified system with effective conduit network that transmits infiltration to springs under high flow conditions and the ab- sence or presence of only thin overlying layers that fail to dampen the springs’ response to rainfall events. Mobile com- pounds including Cl- and SO4 2-, which tend to accumulate in soil solution due to elevated evapotranspiration processes dur- ing the summer months (GOLDSCHEIDER & DREW, 2007), display a concentration peak with the first significant autumn rainfall events (Fig. 4). The measured levels of Cl-, SO4 2- and NO3 - at all sampling sites remained appreciably below the maximum allowed concentrations (250 mg/L, 250 mg/L, and 50 mg/L, respectively) specified in the Drinking Water Direc- tive (2020/2184). In contrast to the springs and groundwater, the Cetina River displayed a greater degree of variability in chemical composition, particularly with respect to concentrations of Cl-, Na+, HCO3 -, and K+ ions (Table 1, Fig. 4). The strongest posi- tive correlation of R=0.99 was observed between Na+ and Cl- Figure 3. A) Biplots of (Ca2+ + Mg2+) versus (HCO3 −). B) (Ca2+ + Mg2+) versus (HCO3 − + SO4 2−). G eo lo gi a C ro at ic a 150 Geologia Croatica 77/2 ions for the Jadro, Žrnovnica, springs and the Cetina River (Fig. S1). Kruskal-Wallis post hoc Dunn test revealed no sta- tistically significant distinction in terms of Cl- ion content among the observed locations (except for Cetina vs. Gizdavac) (Table S2). This lends support to the presence of intercatch- ment groundwater flow, whereby groundwater from the Cetina catchment traverses towards the Jadro and Žrnovnica springs. Elevated Na+ and Cl- concentrations are usually anticipated within coastal aquifers affected by seawater intrusion or input of airborne salts (LUKAČ REBERSKI et al., 2020; PLAN- Figure 4. Multivariate chemographs showing the temporal variation in major ion and saturation indices for calcite and dolomite observed at the Jadro and Žrnovnica springs, Cetina River and in groundwater from the Gizdavac borehole. Total EOCs concentration per location (bar graph) were reported in SELAK et al. (2022a). X on the x-axis marks the sampling campaign with no EOCs detected above limits of detection. G eologia C roatica 151Selak et al.: Hydrochemical characterization of a Dinaric karst catchment in relation to emerging organic contaminants TAK et al., 2021; PATEKAR et al., 2022; FILIPOVIĆ et al., 2023). In our case, the impermeable Eocene flysch coastal belt hinders Adriatic Sea intrusion into the land, as observed by FRITZ (1979). However, strong positive correlations between SO4 2-, Mg2+, Na+, and Cl- in all water samples, and between those ions and K+ in the Jadro spring and Cetina River may indicate sea spray influence (Fig. S1). Furthermore, the pres- ence of chloride may stem from evaporite dissolution (SCHMIDT et al., 2013) or anthropogenic sources including wastewater and fertilizers (DOGDEN et al., 2017), warranting additional research for comprehensive elucidation. Elevated nitrate concentrations are typically associated with agricultural practices, encompassing the application of synthetic fertilizers (e.g. N/P/K and ammonium sulphate), the utilization of manure in cultivated fields, and septic systems discharges (DOGDEN et al., 2017). The mean NO3 - concen- trations observed within the research area remained below the recorded average in Croatian Dinaric karst groundwater, which stands at approx. 5 mg/L (ONDRASEK et al., 2021). Over the course of six sampling campaigns at the deep bore- hole site, groundwater exhibited NO3 - concentrations some- what higher than those observed at other sampling locations, with mean values of 3.02 mg/L compared to 2.31 mg/L for the Jadro, and 1.66 mg/L for the Žrnovnica springs, and 1.24 mg/L for the Cetina River (Table 1, Fig. 4). LOBOREC et al. (2015) observed the trajectory of increasing nitrate concentrations in the Jadro spring over nearly four decades (1975-2014), substan- tiating the growing anthropogenic pressures exerted upon the catchment area. The presence of nitrates in groundwater may suggest rapid infiltration processes (CELLE-JEANTON et al., 2003), which is also indicated by the occurrence of biodegrad- able EOCs gabapentin identified in groundwater from the Giz- davac borehole (Selak et al., 2024). NO3 - displayed a strong positive correlation with SO4 2- in the Žrnovnica spring (R=0.77) and groundwater from the Gizdavac well (R=0.81) (Fig. S1), as well as with K+ in groundwater (R=0.75). These ions could have similar anthropogenic sources in fertilizers and wastewater (HAN et al., 2016; MASS et al., 2019; TOR- RES-MARTÍNEZ et al., 2020). As previously mentioned SO4 2- may originate from the dissolution of gypsum deposits north of Muć, highlighting the intricate interplay between in- herent hydrogeological processes and potential site-specific pollution sources. Given the absence of concurrent stable ni- trogen isotope analysis alongside nitrate measurement, a defi- nite attribution of nitrate concentrations in groundwater to wastewater or agriculture pollution sources cannot be conclu- sively established. The strategically situated borehole along the regional Muć-Gizdavac-Klis fault, delineates a prevailing groundwater flow direction from the NW to the SE towards the Jadro and Žrnovnica springs (identified through a tracer test from the Jablan ponor with high apparent groundwater flow velocities of 10.58 cm/s and 12.18 cm/s, respectively) (FRITZ, 1979). According to this, there is the potential con- veyance of agricultural pollutants including fertilizers and pesticides from Muć polje towards the sampled groundwater and spring locations within a relatively short time. It should be noted how recent research questions the reliability of the tracer tests conducted, as the electrical conductivity of the Jadro spring does not display the rainfall effect from Muć despite the highest apparent groundwater velocity (JUKIĆ et al., 2022). Further investigation is necessary to deepen our understand- ing of the sources and pathways of anthropogenic contami- nants within the studied catchment. 4.2. Stable water isotopes Table 2 provides a descriptive statistical overview of the stable water isotopes δ18Ο and δ2H, including deuterium excess, for springs, river, groundwater, and precipitation. Cetina River δ18Ο data deviated from a normal distribution (Table S1). The results for the Jadro and Žrnovnica springs corroborates earlier research by KAPELJ (2001; published in BRKIĆ et al., 2020) and POLJAK (2014), aligning with the average stable water isotopic values observed in Dalmatian karst springs, which typically vary between -6.7‰ and -8.7‰ for δ18Ο and -39.7‰ and -58‰ for δ2H (BRKIĆ et al., 2020). The isotopic content of precipitation is slightly depleted compared to the data presented by HUNJAK et al. (2013) for the Split station, (located outside of catchment), that had an average of -5.5‰ for δ18O, -34.7‰ for δ2H and d-excess of 9.2 (for period 2007- 2010). This is to be expected as our rain gauges were positioned in the hinterland and at higher altitudes than the Split station (located near the coast at 40 m a.s.l.; Fig. 1). The average d-excess values of karst groundwater, (springs and borehole) and surface water, correspond to the d-excess values of West Mediterranean precipitation (GAT & CARMI, 1970; CELLE- JEANTON et al., 2001). Precipitation at the Dugopolje rain Table 2. Elementary statistics of isotopic signatures for all sampling points (data period October 2019 – October 2022 for springs, river, and groundwater, and data period March 2021 – October 2022 for precipitation). Sampling location Altitude (m.a.s.l.) δ18O (‰) δ2H (‰) d-excess (‰) min max mean sd min max mean sd mean Jadro 35 –7.82 –7.05 –7.47 0.17 –47.31 –41.44 –44.76 1.43 15 Žrnovnica 90 –8.05 –7.09 –7.56 0.22 –47.58 –40.79 –44.97 1.94 15.53 Cetina ~294 –8.46 –7.20 –8.13 0.3 –52.93 –44.06 –50.1 2.02 14.96 Gizdavac ~356 –7.09 –6.66 –6.94 0.17 –41.85 –39.31 –40.89 0.96 14.62 Rain gauge Dugopolje* 497 –9.09 –1.59 –6.42*** 1.94 –58.45 –6.33 –39.72*** 14.16 11.66*** Rain gauge Mosor** 869 –10.25 –3.43 –7.45*** 1.85 –62.53 –13.69 –43.33*** 13.47 16.28*** Due to insufficient rain amounts collected in the rain gauge, the following data is missing: *no data for July 2021 and 2022; **no data from July to November 2021, and for July 2022. *** amount-weighted mean values. G eo lo gi a C ro at ic a 152 Geologia Croatica 77/2 gauge displayed a weighted-mean d-excess of 11.66‰, which is close to the global mean d-excess of 10‰ associated with Atlantic air masses, but indicates a slight Mediterranean influence. The Mosor rain gauge recorded predominantly higher weighted-mean d-excess values (16.28‰), signifying a stronger Mediterranean influence, but also the altitude effect (the rain gauge is located at 869 m a.s.l. in a mountainous area). Water isotopic content is affected by meteoric processes and orographic effects (BRKIĆ et al., 2020). KAPELJ et al. (2012) determined an altitude gradient variation in δ18O, ranging from -0.2‰ to -0.4‰/100 m. Slightly more negative δ18O and δ2H values for Žrnovnica than Jadro can be attributed to the altitude effect, i.e. recharge with precipitation from higher altitudes. The Cetina River had an even lower average δ18O and δ2H values, which is attributed to a recharge area of higher altitudes and a stronger continental influence. In June 2021, there were notable similarities in the stable isotopic composition, particularly for δ18O and Cl- content (Fig. 4), between both springs and the Cetina River. This observation hints at inter- catchment groundwater flow coming from the Cetina River to the Jadro and Žrnovnica springs. Such inflows from the adjacent Cetina River catchment were identified by JUKIĆ & DENIĆ-JUKIĆ (2015) in the partial auto-correlation function of the discharge signal. This phenomenon is likely more pronounced during dry periods, as indicated by the increase in the minimum flow rates at both springs (KADIĆ et al., 2017). 4.3. Correlation between hydrochemical parameters and emerging organic contaminants Exploring relationships between emerging organic contami- nants and other chemical indicators helps to understand the sources, occurrence, and transport of emerging contaminants in karst aquifers (KATZ & GRIFFIN, 2008; HILLEBRAND et al., 2014; ZEMANN et al., 2015). Owing to a limited quan- tity of samples, the correlation between the hydrochemical pa- rameters and concentrations of a particular EOC was solely feasible for the Cetina River. Readers are referred to the work by Selak et al. (2022b) for a comprehensive inventory of ana- lysed and identified EOCs along with their limits of detection and physico-chemical properties. A very mobile pharmaceuti- cal compound metformin (C4H11N5) was the only EOC de- tected more than twice at a single location (recorded in 4 out of 7 campaigns). Metformin exhibited a statistically significant and strong positive correlation with Ca2+ and NO3 - ions in the Cetina River, with correlation coefficients of R=0.72 and R=0.89, respectively (Fig. 5). Additionally, a strong positive correlation (R=0.87) was observed between the total concentration of EOCs and NO3 - ion in samples collected from the Cetina River (Fig. 6). Nitrates, often indicative of anthropogenic contamination from wastewater or agriculture, were previously associated with particular EOCs, including herbicides (HILLEBRAND et al., 2014), X-ray contrast media (ZEMANN et al., 2015), and carbamazepine (DOUMMAR et al., 2014). Nitrates showed a positive correlation with increasing EOCs detection number (ZEMANN et al., 2015; RICHARDS et al., 2023) and total EOC sum (SCHAIDER et al., 2014). Similar to findings of ZEMANN et al. (2015) and RICHARDS et al. (2023), this research identified a statistically significant correlation (R=0.90) between the EOCs number and NO3 - concentrations in the Cetina River (Fig. 6). These positive correlations validate the EOCs sources map and conceptual model presented by SELAK et al. (2024). In groundwater samples, EOCs detection number and to- tal concentration were strongly negatively correlated with HCO3 - and Ca2+ ions, while a strong positive correlation was observed between the EOCs detection number, the total con- centration and Mg2+ ion (Fig. 6). The prevalence of calcium and bicarbonate ions in groundwater, attributed to the pre- dominance of limestone in the aquifer rocks, suggests that in this tectonically disturbed and highly permeable part of the aquifer with low matrix porosity typical of Dinaric karst aq- uifers, the detected EOCs undergo swift infiltration. This ex- plains the presence of rapidly biodegradable gabapentin in the sample extracted from the Gizdavac borehole. In the Cetina River, we observed strong negative correlations (in descend- ing order of correlation coefficient strength) between the num- ber of EOCs detections, the total concentration sum of EOCs, and SO4 2-, Cl-, and Na+, as well as between EOCs detection rate and Mg2+ and K+ ions, as depicted in Fig. 6. Similar pat- terns were noted in Jadro samples, where the total EOCs con- centration and detection number displayed strong negative correlations with K+, Na+, Cl-, SO4 2-, and Mg2+ (Fig. 6). More- over, Žrnovnica spring samples exhibited strong negative cor- relations between EOCs detection rate and Mg2+, Na+, Cl-, SO4 2-, and K+ ions. Since Cl-, Na+, K+, SO4 2-, and Mg2+ ions can partly originate from sea spray, the negative correlation discerned with the EOCs number and total concentration may suggest the absence of atmospheric input of EOCs in the study Figure 5. Correlation matrix of hydrochemical parameters for Cetina (correlation coefficients order with AOE algorithm - the angular order of the eigenvectors) and EOC metformin. Values marked with X are not significant. G eologia C roatica 153Selak et al.: Hydrochemical characterization of a Dinaric karst catchment in relation to emerging organic contaminants area. Conversely, a strong positive correlation was observed between the EOCs detection number and total concentration sum with respect to Ca2+ content in the Cetina River and EOCs detection number and Ca2+ in the Jadro spring samples. HILLEBRAND et al. (2014) established a robust correlation between calcium ions and the herbicide atrazine, which they attributed to the gradual release of atrazine from the karst rock matrix into the groundwater. Persistent and mobile com- pounds including DEET detected in the Cetina River and Jadro spring, showed limited attenuation within the studied karst aquifer (SELAK et al., 2024). This compound showed elevated concentrations during autumn runoff and persisted even during base flow conditions, suggesting potential stor- age within the epikarst and matrix. 5. CONCLUSIONS This hydrochemical characterization study has provided a comprehensive examination of the Jadro and Žrnovnica karst aquifer in both spatial and temporal dimensions, while analyzing factors influencing the quality of its water resources. The hydrochemical fingerprint of this aquifer is typical for coastal Dinaric aquifers characterised by highly karstified structures, a prevalence of limestone, as indicated by the Ca- HCO3 hydrogeochemical facies, and clear evidence of seaspray influence. Stable isotope signatures confirmed a common catchment area shared by the Jadro and Žrnovnica springs while also revealing intercatchment groundwater f low originating from the Cetina River. Figure 6. Correlation matrix of hydrochemical parameters for the Jadro (A), Žrnovnica (B), Gizdavac (C), and Cetina (D) (correlation coefficients order with AOE algorithm - the angular order of the eigenvectors), EOCs detection number and total concentration sum. Values marked with X are not significant. G eo lo gi a C ro at ic a 154 Geologia Croatica 77/2 Alongside the hydrodynamics, hydrochemical characte­ rization has a pivotal role in comprehending the complexity of karst systems and assessing their susceptibility to EOCs. Sharp spikes in chemographs, alongside rapid hydrodynamic responses are evidence of a highly karstified and inherently vulnerable system. Its well­developed conduit network can swiftly transport EOCs via infiltrated water to springs under high flow conditions, as the absence of or thin overlying “protective” layers fail in their attenuation. While conventional indicators of anthropogenic pollution including NO3 ­, Cl­, and SO4 2­ remain below their respective maximum allowable concentrations in drinking water, detections of EOCs reflect the anthropogenic impact on water resources within the studied catchment. Significant strong positive correlation between metformin and NO3 ­ concentration in the Cetina River, as well as significant strong positive correlation between total EOCs concentration, number and NO3 ­ ion, indicates potential contamination stemming from wastewater or agriculture. Robust significant positive correlation between EOCs detection rate and Ca2+ content observed in the Jadro spring samples suggests the persistence of certain EOCs including DEET within the epikarst and aquifer matrix. To enhance the reliability of multivariate analyses of hydrochemical parameters and EOCs, we recommend that future research should involve a larger sample number and focus on monitoring prioritized EOCs. In the context of this Dinaric karst catchment, we propose directing research efforts towards persistent and mobile EOCs, for example; DEET, 1­H benzotriazole, and metformin. The future long­term monitoring and concurrent utilization of hydrochemical markers and novel tracers including EOCs will offer enhanced clarity in discerning the intricate interplay of dominant factors shaping hydrochemical patterns within karst aquifers subject to various anthropogenic pressures. ACKNOWLEDGEMENT The authors express their gratitude to the Croatian Meteorological and Hydrological Service and Croatian Waters for the provided data. This work was financially supported by Interreg Central Europe boDEREC­CE project co­funded by the European Regional Development Fund, GeoTwinn project funded under H2020­EU.4.b., [grant number 809943] and NONa project funded by National Recovery and Resilience Plan 2021–2026 of the European Union – NextGenerationEU, and monitored by the Ministry of Science and Education of the Republic of Croatia. The field investigations and EOCs analysis were funded and supported by the Croatian Geological Survey, Department of Hydrogeology and Engineering Geology. REFERENCES APPELO, C.A.J., & POSTMA, D. (2005): Geochemistry, groundwater and pollution.– 2nd ed., A.A. Balkema publishers, 672 p. BONACCI, O., & ANDRIĆ, I. (2015): Karst spring catchment: an example from Dinaric karst.– Environ Earth Sci, 74/7, 6211–6223. doi: 10.1007/ s12665­015­4644­8 BONACCI, O. & ROJE-BONACCI, T. (1997): Hidrološki vid određivanja bi­ ološkog minimuma rijeke Žrnovnice. [Hydrological aspects of determi- ning the biological minimum of the Žrnovnica River – in Croatian].– Hrvat Vode, 5, 1997, 339–349. BGR, IAH, KIT, UNESCO (2017): World Karst Aquifer Map, 1:40 000 000. Berlin, Reading, Karlsruhe and Paris. URL: https://gdk.gdi-de.org/geo­ network/srv/api/records/473d851c­4694­4050­a37f­e421170eca8. BRKIĆ, Ž., KUHTA, M., HUNJAK, T., & LARVA, O. (2020): Regional Iso­ topic Signatures of Groundwater in Croatia.– Water, 12/7. doi: 10.3390/ w12071983 CAETANO BICALHO, C., BATIOT­GUILHE, C., SEIDEL, J.L., VAN EXTER, S. & JOURDE, H. (2012): Geochemical evidence of water so­ urce characterization and hydrodynamic responses in a karst aquifer.– J Hydrol, 450–451, 206–218. doi: 10.1016/j.jhydrol.2012.04.059 CELLE­JEANTON, H., EMBLANCH, C., MUDRY, J. & CHARMOILLE, A. (2003): Contribution of time tracers (Mg2+, TOC, δ13CTDIC, NO3 −) to understand the role of the unsaturated zone: A case study – Karst aqu­ ifers in the Doubs valley, eastern France.– Geophys Res Lett, 30/6. doi: 10.1029/2002GL016781 CELLE­JEANTON, H., TRAVI, Y. & BLAVOUX, B. (2001): Isotopic typo­ logy of the precipitation in the Western Mediterranean Region at three different time scales.– Geophys Res Lett, 28/7, 1215–1218. doi: 10.1029/2000GL012407 CIZMAS, L., SHARMA, V.K., GRAY, C.M. & MCDONALD, T.J. (2015): Pharmaceuticals and personal care products in waters: occurrence, toxi­ city, and risk.– Environ Chem Lett, 13/4, 381–394. doi: 10.1007/s10311­ 015­0524­4 CLARK, I., & FRITZ, P. (1997): Environmental Isotopes in Hydrology.– Lewis Publication. ĆUK ĐUROVIĆ, M., PETRIČ, M., JEMCOV, I., MULEC, J., GRUDNIK, Z. M., MAYAUD, C., BLATNIK, M., KOGOVŠEK, B. & RAVBAR, N. (2022): Multivariate Statistical Analysis of Hydrochemical and Micro­ biological Natural Tracers as a Tool for Understanding Karst Hydrodyna­ mics (The Unica Springs, SW Slovenia).– Water Resour Res, 58/11, e2021WR031831. doi: 10.1029/2021WR031831 DANSGAARD, W. (1964): Stable isotopes in precipitation.– Tellus, 16/4, 436–468. doi: 10.1111/j.2153­3490.1964.tb00181.x DINNO, A. (2017): Dunn.test: Dunn’s Test of Multiple Comparisons Using Rank Sums. https://cran.r­project.org/package=dunn.test DIRECTIVE (EU) 2020/2184 of the European Parliament and of the Council of 16 December 2020 on the quality of water intended for human con­ sumption. DODGEN, L.K., KELLY, W.R., PANNO, S.V, TAYLOR, S.J., ARMSTRONG, D.L., WILES, K.N., ZHANG, Y. & ZHENG, W. (2017): Characterizing pharmaceutical, personal care product, and hormone contamination in a karst aquifer of southwestern Illinois, USA, using water quality and stream flow parameters.– Sci Total Environ, 578, 281–289. doi: 10.1016/j. scitotenv.2016.10.103 DOMENICO, P.A. & SCHWARTZ, F.W.(1990). Physical and Chemical Hy­ drogeology. John Wiley and Sons, New York, 824 p. DOUMMAR, J., GEYER, T., BAIERL, M., NÖDLER, K., LICHA, T. & SA­ UTER, M. (2014): Carbamazepine breakthrough as indicator for specific vulnerability of karst springs: Application on the Jeita spring, Lebanon.– Appl Geochemistry, 47, 150–156. doi: 10.1016/j.apgeochem.2014.06.004 FILIPOVIĆ, M., FRANGEN, T., TERZIĆ, J. & LUKAČ REBERSKI, J. (2023): Hydrogeology of a complex karst catchment in Southern Dalma­ tia (Croatia) and Western Herzegovina (Bosnia and Herzegovina).– J Maps, 19/1, 2112775. doi: 10.1080/17445647.2022.2112775 FRITZ, F. (1979): Općina Split. Hidrogeološka studija. [Split municipality. Hydrogeological study – in Croatian].– Hrvatski geološki institut. Fond stručne dokumentacije, 191. GAT, J.R. & CARMI, I. (1970): Evolution of the isotopic composition of atmospheric waters in the Mediterranean Sea area.– J Geophys Res, 75/15, 3039–3048. doi: 10.1029/JC075i015p03039 GEOTEHNIKA (1975): Akumulacija Đale – Pribransko rješenje. Istraživač­ ko projektna dokumentacija. Idejni projekt za brtvljenje akumulacijskog bazena [Đale reservoir – dam design. Research project documentation. Conceptual design for sealing the accumulation basin. – in Croatian].– Fond struč. dok. Elektroprivrede Dalmacije, Split. G eologia C roatica 155Selak et al.: Hydrochemical characterization of a Dinaric karst catchment in relation to emerging organic contaminants GOLDSCHEIDER, N. & DREW, D. (2007): Methods in Karst Hydrogeo­ logy.– IAH: International Contributions to Hydrogeology, 26 (1st ed.). CRC Press, 260 p. doi: 10.1201/9781482266023 HAN, D., SONG, X., & CURRELL, M.J. (2016): Identification of anthropo­ genic and natural inputs of sulfate into a karstic coastal groundwater system in northeast China: evidence from major ions, δ13CDIC and δ34SSO4.– Hydrol Earth Syst Sci, 20/5, 1983–1999. doi: 10.5194/he­ ss­20­1983­2016 HARRELL, JR. F. (2023): Hmisc: Harrell Miscellaneous. R package version 5.1­0, https://CRAN.R­project.org/package=Hmisc. HILLEBRAND, O., NÖDLER, K., GEYER, T., & LICHA, T. (2014): Investi­ gating the dynamics of two herbicides at a karst spring in Germany: Consequences for sustainable raw water management.– Sci Total Envi­ ron, 482–483, 193–200. doi: 10.1016/j.scitotenv.2014.02.117 HOLLANDER, M. & WOLFE, D.A. (1973): Nonparametric Statistical Methods.– John Wiley, New York, New York. HUNJAK, T., LUTZ, H.O. & ROLLER-LUTZ, Z. (2013): Stable isotope com­ position of the meteoric precipitation in Croatia.– Isotopes Environ He­ alth Stud, 49/3, 336–345. doi: 10.1080/10256016.2013.816697 JEBREEN, H., BANNING, A., WOHNLICH, S., NIEDERMAYR, A., GHA­ NEM, M. & WISOTZKY, F. (2018): The Influence of Karst Aquifer Mi­ neralogy and Geochemistry on Groundwater Characteristics: West Bank, Palestine.– Water, 10/12. doi: 10.3390/w10121829 JUKIĆ, D. & DENIĆ-JUKIĆ, V. (2008): Estimating parameters of ground- water recharge model in frequency domain: Karst springs Jadro and Žrnovnica.– Hydrol Process, 22/23, 4532–4542. doi: 10.1002/hyp.7057 JUKIĆ, D. & DENIĆ-JUKIĆ, V. (2015): Investigating relationships between rainfall and karst­spring discharge by higher­order partial correlation functions.– J Hydrol, 530, 24–36. doi: 10.1016/j.jhydrol.2015.09.045 JUKIĆ, D., DENIĆ-JUKIĆ, V. & KADIĆ, A. (2022): Temporal and spatial characterization of sediment transport through a karst aquifer by means of time series analysis.– J Hydrol, 609. doi: 10.1016/j.jhydrol.2022.127753 KADIĆ, A., DENIĆ-JUKIĆ, V. & JUKIĆ, D. (2017): Revealing hydrological relations of adjacent karst springs by partial correlation analysis.– Hydrol Res, 49/3, 616–633. doi: 10.2166/nh.2017.064 KADIĆ, A., DENIĆ-JUKIĆ, V. & JUKIĆ, D. (2019): Analiza meteoroloških i hidroloških odnosa u kršu primjenom parcijalne kros­korelacijske funk­ cije višeg reda. [Analysis of meteorological and hydrological relations in the karst using higher-order partial cross-correlation function – in Croatian].– Hrvatske vode, 27/109, 201­210. KAPELJ, S., KAPELJ, J., BIONDIĆ, R., BIONDIĆ, B., KOVAČ, I., TUŠAR, B., PRELOGOVIĆ, E., MARJANAC, T., ANDRIĆ, M., KOVAČIĆ, D., STRELEC, S. & GAZDEK, M. (2006): Studija upravljanja vodama sliva Jadra i Žrnovnice – Prva faza studijsko istraživačkih radova [Water ma- nagement study of the Jadro and Žrnovnica springs recharge area – first phase – in Croatian]. EVV:1/2005.– Hrvatske vode, Split KAPELJ, S., KAPELJ, J., NOVOSEL, A. & SINGER, D. (2001): Hidrogeo­ loška istraživanja slivnog područja izvora Jadra i Žrnovnice - prelimi­ narni izvještaj, I faza istraživanja [Hydrogeological study of the Jadro and Žrnovnica catchment area – preliminary report, first phase – in Croatian]. Fond str. dok. HGI, br. 42/01, Zagreb. KAPELJ, S., KAPELJ, J. & ŠVONJA, M. (2012): Hidrogeološka obilježja sliva Jadra i Žrnovnice [Hydrogeological characteristics of the Jadro and Žrnovica catchement – in Croatian].– Tusculum, 5/1, 205–216. KATZ, B.G. & GRIFFIN, D.W. (2008): Using chemical and microbiological indicators to track the impacts from the land application of treated mu­ nicipal wastewater and other sources on groundwater quality in a karstic springs basin.– Environ Geol, 55/4, 801–821. doi: 10.1007/s00254­007­ 1033­y LLAMAS, M.I., JIMÉNEZ­GAVILÁN, P., LUQUE­ESPINAR, J.A., BENA­ VENTE­HERRERA, J., CANDELA, L., SANMIGUEL­MARTÍ, M., RAMBLA­NEBOT, J., ARANDA­MARES, J.L., & VADILLO­PÉREZ, I. (2022): Hydrogeological, hydrodynamic and anthropogenic factors affecting the spread of pharmaceuticals and pesticides in water resources of the Granada plain (Spain).– J Hydrol, 610, 127791. doi: 10.1016/j.jhy­ drol.2022.127791 LOBOREC, J., KAPELJ, S. & NOVAK, H. (2015): Analysis of groundwater pollution hazard in karst: a case study of Jadro and Žrnovnica catchment area.– Građevinar, 67/11, 1093–1103. doi: 10.14256/JCE.1250.2015 LUKAČ REBERSKI, J., RUBINIĆ, J., TERZIĆ, J. & RADIŠIĆ, M. (2020): Climate Change Impacts on Groundwater Resources in the Coastal Kar­ stic Adriatic Area: A Case Study from the Dinaric Karst.– Nat Resour Res, 29/3, 1975–1988. doi: 10.1007/s11053­019­09558­6 LUKAČ REBERSKI, J., SELAK, A., LAPWORTH, D. J., MAURICE, L. D., TERZIĆ, J., CIVIL, W. & STROJ, A. (2023): Emerging organic conta­ minants in springs of the highly karstified Dinaric region.– J Hydrol, 621, 129583. doi: 10.1016/j.jhydrol.2023.129583 MALDINI, K., CUKROV, N., PIKELJ, K., MATIĆ, N. & MLAKAR, M. (2023): Geochemistry of Metals and Organic Matter in Water and Sedi­ ments of the Karst River Cetina, Croatia.– Water, 15/7. doi: 10.3390/ w15071429 MAAS, B., PETERSON, E.W., HONINGS, J., OBERHELMAN, A., OWA­ RE, P., RUSTHOVEN, I. & WATSON, A. (2019): Differentiation of Sur­ face Water and Groundwater in a Karst System Using Anthropogenic Signatures.– Geosciences, 9/4. doi: 10.3390/geosciences9040148 MATIATOS, I., ALEXOPOULOS, A. & GODELITSAS, A. (2014): Multiva­ riate statistical analysis of the hydrogeochemical and isotopic composi­ tion of the groundwater resources in northeastern Peloponnesus ( Greece).– Sci Total Environ, 476/477, 577–590. doi: 10.1016/j.scito­ tenv.2014.01.042 MATIĆ, N., MALDINI, K., CUCULIĆ, V. & FRANČIŠKOVIĆ-BILINSKI, S. (2012): Investigations of karstic springs of the Biokovo Mt from the Dinaric karst of Croatia.– Geochemistry, 72/2, 179–190. doi: 10.1016/j. chemer.2011.08.001 MUDARRA, M., ANDREO, B. & BAKER, A. (2011): Characterisation of dissolved organic matter in karst spring waters using intrinsic fluores­ cence: Relationship with infiltration processes.– Sci Total Environ, 409/18, 3448–3462. doi: 10.1016/j.scitotenv.2011.05.026 OGLE, D.H., DOLL, J.C., WHEELER, A.P., DINNO, A. (2023): FSA: Simple Fisheries Stock Assessment Methods. R package version 0.9.5, https:// CRAN.R­project.org/package=FSA. OKSANEN, J., SIMPSON, G., BLANCHET, F. G., KINDT, R., LEGENDRE, P., MINCHIN, P., HARA, R., SOLYMOS, P., STEVENS, H., SZÖCS, E., WAGNER, H., BARBOUR, M., BEDWARD, M., BOLKER, B., BORCARD, D., CARVALHO, G., CHIRICO, M., DE CÁCERES, M., DURAND, S. & WEEDON, J. (2022): Vegan community ecology pac­ kage version 2.6­2 April 2022. ONDRASEK, G., BAKIĆ BEGIĆ, H., ROMIĆ, D., BRKIĆ, Ž., HUSNJAK, S. & BUBALO KOVAČIĆ, M. (2021): A novel LUMNAqSoP approach for prioritising groundwater monitoring stations for implementation of the Nitrates Directive.– Environ Sci Eur, 33/1, 23. doi: 10.1186/s12302­ 021­00467­1 PADILLA, I.Y. & VESPER, D.J. (2018): Fate, Transport, and Exposure of Emerging and Legacy Contaminants in Karst Systems: State of Knowled­ ge and Uncertainty BT.– In: WHITE, W.B., HERMAN, J.S., HERMAN, E.K. & RUTIGLIANO, M. (eds.): Karst Groundwater Contamination and Public Health. Springer International Publishing, 33–49. PARKHURST, D.L. & APPELO, C.A.J. (1999): User’s Guide to PHREEQC (Version 2)—A Computer Program for Speciation, Batch­Reaction, One­Dimensional Transport, and Inverse Geochemical Calculations. – U.S. Geological Survey, Water Resources Investigations Report, 99­4259, Washington DC. PATEKAR, M., BAŠIĆ, M., POLA, M., KOSOVIĆ, I., TERZIĆ, J., LUCCA, A., MITTEMPERGHER, S., BERIO, L. R. & BOROVIĆ, S. (2022): Mul­ tidisciplinary investigations of a karst reservoir for managed aquifer re­ charge applications on the island of Vis (Croatia).– Acque Sotter – Ital J Groundw, 11/1 37–48. doi: 10.7343/as­2022­557 PERRIN, J. (2003): A conceptual model of flow and transport in a karst aqu­ ifer based on spatial and temporal variations of natural tracers. PIPER, M.A. (1944): A graphic procedure in the geochemical interpretation of water­analyses.– Eos, Trans Am Geophys Union, 25/6, 914–928. doi: 10.1029/TR025i006p00914 G eo lo gi a C ro at ic a 156 Geologia Croatica 77/2 PLANTAK, L., BIONDIĆ, R., MEAŠKI, H. & TEŽAK, D. (2021): Hydroc­ hemical Indicators Analysis of Seawater Intrusion in Coastal Karstic Aquifers on the Example of the Bokanjac-Poličnik Catchment Area in Zadar, Croatia.– Applied Sciences, 11/24. doi: 10.3390/app112411681 POLJAK, (2014): Određivanje srednjeg vremena zadržavanja podzemne vode u kršu primjenom prirodnih izotopa - primjer sliva Jadro i Žrnovnica. [Determination of mean residence time of groundwater in karst using natural isotopes – example of Jadro and Žrnovnica catchment – in Croatian].– PhD Thesis, University of Zagreb, Faculty of Geotechnical Engineering. Varaždin. PULIDO-BOSCH, A. (2021): Principles of Karst Hydrogeology.– Springer International Publishing. doi: 10.1007/978­3­030­55370­8 R CORE TEAM (2023): R: A Language and Environment for Statistical Com­ puting.– Foundation for Statistical Computing, Vienna, Austria, https:// www.R­project.org RICHARDS, L.A., GUO, S., LAPWORTH, D.J., WHITE, D., CIVIL, W., WILSON, G.J.L., LU, C., KUMAR, A., GHOSH, A., KHAMIS, K., KRAUSE, S., POLYA, D.A. & GOODDY, D.C. (2023): Emerging organic contaminants in the River Ganga and key tributaries in the middle Gangetic Plain, India: Characterization, distribution & controls. ­Environ Pollut, 327, 121626. doi: 10.1016/j.envpol.2023.121626 SCHAIDER, L.A., RUDEL, R.A., ACKERMAN, J.M., DUNAGAN, S.C., & BRODY, J.G. (2014): Pharmaceuticals, perfluorosurfactants, and other organic wastewater compounds in public drinking water wells in a shallow sand and gravel aquifer.– Sci Total Environ, 468–469, 384–393. doi: 10.1016/j.scitotenv.2013.08.067 SCHMIDT, S., GEYER, T., MAREI, A., GUTTMAN, J. & SAUTER, M. (2013): Quantification of long-term wastewater impacts on karst groun­ dwater resources in a semi­arid environment by chloride mass balance methods.– J Hydrol, 502, 177–190. doi: 10.1016/j.jhydrol.2013.08.009 SCHRIKS, M., HERINGA, M.B., VAN DER KOOI, M.M.E., DE VOOGT, P. & VAN WEZEL, A.P. (2010): Toxicological relevance of emerging contaminants for drinking water quality.– Water Res, 44/2, 461–476. doi: 10.1016/j.watres.2009.08.023 SELAK, A., LUKAČ REBERSKI, J., BOLJAT, I. & TERZIĆ, J. (2024): Cha­ racterizing occurrence of emerging organic contaminants in Dinaric karst catchment of Jadro and Žrnovnica springs, Croatia.– Emerg Contam, 10/3, 100327. doi: 10.1016/j.emcon.2024.100327 SELAK, A., LUKAČ REBERSKI, J., KLOBUČAR, G. & GRČIĆ, I. (2022a): Ecotoxicological aspects related to the occurrence of emerging conta­ minants in the Dinaric karst aquifer of Jadro and Žrnovnica springs.– Sci. Total Environ., 825, 153827. doi:10.1016/j.scitotenv.2022.153827 SELAK, A., LUKAČ REBERSKI, J., KLOBUČAR, G. & GRČIĆ, I. (2022b): Data on occurrence and ecotoxicological risk of emerging contaminants in Dinaric karst catchment of Jadro and Žrnovnica springs.– Data Br., 42, 108157. doi:10.1016/j.dib.2022.108157 SHEIKHY NARANY, T., BITTNER, D., DISSE, M. & CHIOGNA, G. (2019): Spatial and temporal variability in hydrochemistry of a small­scale do­ lomite karst environment.– Environ Earth Sci, 78/9, 273. doi: 10.1007/ s12665­019­8276­2 SPEARMAN, C. (1904): The Proof and Measurement of Association between Two Things.– Am. J. Psychol. 15, 72–101. ŠEGOTA, T., & FILIPČIĆ, A. (2003): Köppen’s classification of climates and the problem of corresponding Croatian terminology.– Geoadria, 8/1, 17–37. http://hrcak.srce.hr/index.php?show=clanak&id_clanak_ je­ zik=14843 TERZIĆ, J., MARKOVIĆ, T. & LUKAČ REBERSKI, J. (2014): Hydrogeolo­ gical properties of a complex Dinaric karst catchment: Miljacka Spring case study.– Environ Earth Sci, 72/4, 1129–1142. doi: 10.1007/s12665­ 013­3031­6 TORAN, L. & REISCH, C. E. (2013): Using Stormwater Hysteresis to Cha­ racterize Karst Spring Discharge.– Groundwater, 51/4, 575–587. doi: 10.1111/j.1745­6584.2012.00984.x TORRES-MARTÍNEZ, J. A., MORA, A., KNAPPETT, P. S. K., ORNELAS- SOTO, N. & MAHLKNECHT, J. (2020): Tracking nitrate and sulfate sources in groundwater of an urbanized valley using a multi­tracer ap­ proach combined with a Bayesian isotope mixing model.– Water Res, 182, 115962. doi: 10.1016/j.watres.2020.115962 USGS (2021): PHREEQC v. 3.7.3.15968. Software release 5 December 2021. URL: https://www.usgs.gov/software/phreeqc­version­3 VASIĆ, L., ŽIVOJINOVIĆ, D. & RAJAKOVIĆ-OGNJANOVIĆ, V. (2020): Hydrochemical changes and groundwater grouping data by multivariate statistical methods within one karst system: recharge–discharge zone (Eastern Serbia case study).– Carbonates and Evaporites, 35/1, 15. doi: 10.1007/s13146­019­00548­6 WEI, T. & SIMKO, V. (2021): R package ‘corrplot’: Visualization of a Corre­ lation Matrix. (Version 0.92), https://github.com/taiyun/corrplot. WHITE, W. (1997): Thermodynamic equilibrium, kinetics, activation barri­ ers, and reaction mechanisms for chemical reactions in Karst Terrains.– Environ Geol, 30, 46–58. doi: 10.1007/s002540050131 YURTSEVER, Y. & GAT, J. (1981): Atmospheric waters. Stable Isotope Hydrology: Deuterium and Oxygen­18 in the Water Cycle.– In: GAT, J.R. & GONFIANTINI, R. (eds.): Technical Reports Series 210, IAEA, Vienna, 103–142. ZEMANN, M., WOLF, L., GRIMMEISEN, F., TIEHM, A., KLINGER, J., HÖTZL, H. & GOLDSCHEIDER, N. (2015): Tracking changing X­ray contrast media application to an urban-influenced karst aquifer in the Wadi Shueib, Jordan.– Environ Pollut, 198, 133–143. doi: 10.1016/j.en­ vpol.2014.11.033 ŽIVANOVIĆ, V., JEMCOV, I. & DRAGIŠIĆ, V. (2022): Vulnerability meth­ ods in hard rock formation as a basis for groundwater risk assessment – from resource to source.– Geol Croat, 75, 381–392. doi: 10.4154/ gc.2022.23 https://www.usgs.gov/software/phreeqc-version-3 G eologia C roatica 157Selak et al.: Hydrochemical characterization of a Dinaric karst catchment in relation to emerging organic contaminants Figure S1. Correlation matrix of physicochemical parameters for Jadro (A), Žrnovnica (B), Gizdavac (C), and Cetina (D) (correlation coefficients order with AOE algorithm – the angular order of the eigenvectors). Values marked with X are not significant. Supplement 1. – Figure S1 G eo lo gi a C ro at ic a 158 Geologia Croatica 77/2 Table S1. Results of Shapiro-Wilk test (values in bold indicate normally distributed data i.e. p-value>0.05) (data October 2019 – October 2022). Jadro Žrnovnica Gizdavac Cetina Parameter W p-value W p-value W p-value W p-value Temperature 0.96 0.35 0.98 0.76 0.90 0.40 0.96 0.42 pH 0.89 0.01 0.92 0.04 0.89 0.34 0.96 0.49 Electrical conductivity 0.98 0.83 0.97 0.52 0.93 0.61 0.89 0.01 Na+ 0.86 2.36E-03 0.85 1.41E-03 0.89 0.30 0.84 8.02E-04 K+ 0.97 0.73 0.84 1.03E-03 0.80 0.06 0.87 2.93E-03 Mg2+ 0.92 4.75E-02 0.94 0.14 0.71 0.01 0.99 1.00 Ca2+ 0.97 0.62 0.90 0.02 0.90 0.36 0.87 3.01E-03 Cl– 0.84 7.73E-04 0.82 4.68E-04 0.78 0.04 0.85 1.25E-03 HCO3 – 0.95 0.28 0.97 0.50 0.95 0.77 0.81 2.28E-04 NO3 – 0.76 3.74E-05 0.79 1.02E-04 0.79 0.05 0.78 7.93E-05 SO4 2– 0.96 0.40 0.97 0.54 0.87 0.23 0.92 0.04 δ18O 0.96 0.48 0.96 0.43 0.83 0.11 0.86 2.40E-03 δ2H 0.97 0.71 0.94 0.16 0.86 0.18 0.92 0.05 O2 (%) 0.76 4.48E-05 0.93 0.07 0.95 0.77 0.91 0.02 Jadro Žrnovnica Gizdavac Cetina Table S2. Results of Kruskal-Wallis post hoc Dunn test (values in bold indicate significant difference i.e. p.adjusted<0.05) (data October 2019 – October 2022). Comparison P. adjusted Ca2+ Mg2+ Na+ Cl– HCO3 – SO4 2– NO3 – Cetina – Gizdavac 1.72E-08 1 1.28E-04 1.94E-03 8.49E-08 5.06E-06 5.11E-06 Cetina – Jadro 5.14E-12 1 0.34 0.34 2.17E-08 0.13 1.44E-07 Cetina – Žrnovnica 5.40E-03 0.73 0.3 0.37 6.85E-03 4.60E-05 0.8 Gizdavac – Jadro 0.71 1 0.01 0.09 0.24 2.64E-03 0.13 Gizdavac – Žrnovnica 5.63E-04 0.42 0.01 0.08 1.40E-03 0.17 2.64E-03 Jadro – Žrnovnica 7.63E-04 0.04 1 1 0.05 0.18 6.21E-03 Table S3. Main statistical descriptors of saturation indices (SI) calculated with PHREEQC for Jadro, Žrnovnica, Gizdavac, and Cetina samples (data October 2019 – October 2022). Site Statistics Calcite (CaCO3) Dolomite (CaMg(CO3)2) Anhydrite (CaSO4) Gypsum (CaSO₄·2H₂O) Ja dr o sp rin g Min –0.06 –1.22 –3.41 –2.97 Max 0.71 0.65 –2.60 –2.16 Mean 0.17 –0.51 –2.89 –2.45 Median 0.16 –0.49 –2.87 –2.43 Žr no vn ic a sp rin g Min –0.02 –1.19 –3.49 –3.05 Max 1.16 1.49 –2.76 –2.32 Mean 0.43 –0.07 –3.04 –2.60 Median 0.39 –0.10 –3.03 –2.59 Ce tin a Ri ve r Min 0.37 –0.32 –3.04 –2.55 Max 1.11 1.50 –2.56 –2.13 Mean 0.79 0.77 –2.82 –2.36 Median 0.82 0.77 –2.84 –2.37 G iz da va c bo re ho le Min 0.17 –0.93 –3.47 –3.03 Max 0.45 0.04 –3.13 –2.68 Mean 0.27 –0.43 –3.31 –2.86 Median 0.26 –0.42 –3.31 –2.86 Supplement 2. – Tables S1, S2 and S3.