GEOLOGIJA 50/1, 157–177, Ljubljana 2007 Hydrogeochemical characteristics of the River Sava watershed in Slovenia Hidrogeokemi~ne zna~ilnosti pore~ja reke Save v Sloveniji Tja{a KANDU^ & Nives OGRINC Department of Environmental Sciences, Jožef Stefan Institute, Jamova 39, 1000 Ljubljana, Slovenia Key words: River Sava, chemical composition, stable isotopes, sulphur, carbon, oxygen, deuterium Klju~ne besede: reka Sava, kemijska sestava, stabilni izotopi, žveplo, ogljik, kisik, dev-terij Abstract The River Sava is a typical HCO3- – Ca2+ – Mg2+ River. Total alkalinity increases in the part of the watershed composed of carbonate and clastic rocks, which are less resistant to weathering processes. Ca2+/Mg2+ ratios are around 2 in the carbonate part of the watershed and increase in the watershed composed of carbonate and clastic rocks, indicating dissolution of calcite with magnesium. According to PHREEQC for Windows calculations, the River Sava and its tributaries are oversaturated with respect to calcite and dolomite. ?18OH2O and ?DH2O are related to the meteorological patterns in the drainage basin. River water temperatures fluctuate annually following air temperatures. The relationship between the temperature and ?18OH2O and ?DH2O values primarily reflects the strong dependence of ?18O and ?D on precipitation and evaporative enrichment in heavy oxygen and hydrogen isotopes of infiltrating water recharging the River Sava from its slopes. The ?13CDIC values are controlled by processes in the terrestrial ecosystem and stream processes such as: (1) dissolution of carbonates, (2) soil derived CO2, and (3) equilibration with atmospheric CO2. Lower ?13CDIC values are observed in the spring sampling season due to abundant precipitation related to soil leaching of CO2 in the river system. From discharge and concentration measurements of sulphate and according to the drainage area of the River Sava basin, the annual sulphur flux at the border with Croatia was estimated to be 1.4 × 107 g SO4/km2. Assuming that the sources of SO42- to the Sava are its tributaries, precipitation and other sources, the contributions of these inputs were calculated according to steady state equations and estimated to be 52 : 8 : 40 %, respectively. Other sources are attributed to human influences such as industrial pollution and oxidation of sulphides. Izvle~ek Reka Sava je tipi~na HCO3- – Ca2+ – Mg2+ reka. Totalna alkalnost nara{~a po pore~ju navzdol, ki ga sestavljajo poleg karbonatnih tudi klasti~ne kamnine, ki so manj odporne na procese preperevanja. Razmerje Ca2+/Mg2+ je okrog 2 v karbonatnem delu pore~ja in naraste v pore~ju, ki ga sestavljajo karbonatne in klasti~ne kamnine, kar kaže na raztapljanje kalcita z magnezijem. Glede na izra~une PHREEQC for Windows so reka Sava in njeni pritoki prenasi~eni s kalcitom in dolomitom. ?18OH2O in ?DH2O vrednosti sta povezani z meteorolo{kimi spremenljivkami v drenažnem bazenu. Temperatura re~ne vode se letno spreminja glede na temperature v zraku. Povezava med temperaturo, ?18OH2O odraža mo~no odvisnost ?18O v padavinah in evaporativno obogatitev na težjih izotopih kisika infiltrirajo~e vode s pobo~ij, ki napaja reko Savo. ?13CDIC vrednosti kontrolirajo procesi v terestri~nih ekosistemih in v reki, kot so: (1) raztapljanje karbonatov, (2) preperinski CO2 ter (3) uravnoteženje z atmosferskim CO2. Bolj negativne vrednosti ?13CDIC so posledica ve~je koli~ine padavin in spiranja izotopsko lažjega preperinskega CO2 s terestri~nega ekosistema v spomladanskih mesecih. Iz podatkov o pretokih, koncentracijah sulfata v vodi in odvodnjevalni povr{ini smo izra~unali letni snovni tok žvepla na meji s Hrva{ko, ki zna{a 1,4 × 107 g SO4/km2. Glede nato, da smo predpostavili slede~e vire sulfata v reki Savi: pritoke, padavine in ostale vire, smo izra~unali deleže posameznih prispevkov iz ravnotežnih ena~b koncentracijske in izotopske masne bilance, ki zna{ajo 52 : 8 : 40 %. Ostale vire lahko pripi{emo industrijski onesnaženosti in oksidaciji sulfidnih mineralov. 158 Introduction The geochemical study of river water allows important information to be obtained on chemical weathering of rocks/soil and the chemical and isotopic compositions of the drainage basin (Gibbs, 1972; Reeder et al., 1972; Hu et al., 1982; Stallard & Edmond, 1983; Goldstein & Jacobsen, 1987; Elderfield et al., 1990; Zhang et al., 1995; Huh et al., 1998). Since carbonate weathering largely dominates the water chemistry of river waters, characterization of the water chemistry of rivers draining carbonate-dominated terrain is crucial to precisely identify the various contributions of the different sources to water solutes, and to estimate weathering rates of the continental crust and associated CO2 consumption (Fairchild et al., 1994, 1999, 2000; Gail-lardet et al., 1999; Liu & Zhao, 2000). Rivers also reflect biogeochemical processes occurring in their catchment areas and help to quantify material transport from land to ocean (Palmer et al., 2001). Within this context, understanding of the carbon cycle is particulary important because it helps to evaluate the health of the river and its catchment basin (Telmer & Veizer, 1999). Investigations of major elements and stable carbon isotopes of dissolved inorganic carbon (?13CDIC) are useful for such studies and to evaluate environmental influences (Ka-rim & Veizer, 2000; Barth et al., 2003). Rivers represent a linkage between precipitation and groundwater, which all together form the hydrogeological cycle. Stable isotopes of oxygen and hydrogen are used to trace processes of evaporation, condensation, snow melting, mixing of waters of different origin and recharge conditions in the studied watersheds (Clark & Fritz, 1997). The natural sources of SO42- in river water are rainfall, groundwater, and weathering of S-rich minerals. Human activities (e.g. air pollution, mining, smelting of sulphide ore, refining of petroleum, and chemical industries) also contribute to the SO42- content of river water. These multiple sources of sulphate can frequently be distinguished by their specific isotopic signatures. For example, river sulphate derived from dissolution of evaporates has positive ?34S values whereas sulphate from oxidation of sulphides or from biogenic emissions may have strongly negative ?34S. The latter is due to the fact Tja{a Kandu~ & Nives Ogrinc that sulphide minerals tend to be isotopi-cally lighter than sulphate minerals, and hence the oxidation of sulphides (and that of organic sulphur) transposes the depleted ?34S signal to sulphate products (Pearson & Rightmire, 1980). The main objectives of this study were to understand the hydrogeochemical characteristics of the River Sava watershed through time (different seasons and flow regimes), and to quantify sources, fluxes and sinks of carbon and sulphur from the isotopic characteristics of their dissolved species. Characterization of the River Sava drainage basin The River Sava, the largest river in Slovenia and a tributary of the Danube, originates in the Triassic carbonate hinterland at Zelenci (Figure 1, location 1) as the Sava Dolinka and the karst spring Savica (Figure 1, location 6) as the Sava Bohinjka. The confluence of these two sources is at Radovljica. From there the river is named Sava and finishes its course at Belgrade, merging with the Danube. From the source of the Sava Dolinka to the national border with Croatia its length is 219 km. Along its studied course of 219 km, as well as the Sava Bohinjka source, the river flows through agricultural areas and the urban centres of Jesenice, Kranj, Ljubljana, the mining area of Zasavje and Kr{ko (Figure 1). Discharge regimes along its flow are controlled by precipitation and the configuration of the landscape. In the upper part of the Sava a snow – rain regime prevails and turns over in the central and lower part to a rain – snow regime (Hrvatin, 1998). Annual maxima are characteristics in spring and late summer, while minima occur in summer and winter months. In the years 1961–1990 the mean annual amount of precipitation in the River Sava watershed was 1500 mm (Zupan~i~, 1998). Long term (from the years 1960–1991) the mean annual discharge varies at the gauging stations of Radovljica (in the upper part of the watershed), Hrastnik (in the central part of the watershed) and ~atež (in the lower part of the watershed) from 18.6, 181.7 and 290 m3/s, respectively (ARSO, 2004). Discharge ratios between high and low waters are 1 : 100 (extreme 1 : 250) and are also controlled by hydroelectric power plant outflows. Hydrogeochemical characteristics of the River Sava watershed in Slovenia 159 Figure 1. Detailed location map of the numbered sampling sites in the River Sava watershed. Sample sites are described in Table 1. The drainage area of the River Sava in Slovenia is 10881 km2 (ARSO, 2004). The valley of the Sava extends in a NW– SE direction and comprises almost half the Slovenian territory with variegated geological composition. At the confluence of the Sava Bohinjka and Sava Dolinka it accumulates Pleistocene fluvioglacial sediments and formed terraces, from Ljubljana the watershed is mainly composed of Permo-Carbo-nian clastic sediments, which alternate with Triassic carbonates in the Zasavje area and pass over to Miocene sandstones, clays and gravels on the left bank of the river. In the Kr{ko–Brežice area the watershed mainly consists of terrace Pleistocene sediments and from there the flow continues to Croatia. The watershed of the Sava’s tributaries is composed of Triassic and Jurassic carbonates, Permo-Carbonian, Oligocene, Miocene clastic rocks and Pleistocene sediments (Buser & Draksler, 1989). Methodology Surface water sampling locations were selected based on their relationship to the confluences of major and minor streams, at points before and after the confluence. Sampling locations of the Sava watershed are presented in Figure 1. Sampling was performed at 41 locations in different seasons (April – spring 2004, September – late summer 2004 and January – winter 2005), according to the discharge regimes of the Sava and its tributaries. Temperature, conductivity, dissolved oxygen (DO), and pH measurements were performed in the field. Sample aliquots collected for chemical analysis were immediately passed through a 0.45 µm nylon filter into bottles and kept refrigerated until analysed. Samples for cation (pre-treated with HNO3), anion and alkalinity analyses were preserved in HDPE bottles. Samples for DIC (dissolved inorganic carbon) and ?13CDIC analyses were stored in glass bottles, filled to the top, with no headspace. Samples for DOC analyses were acidified and sampled in glass bottles of volume 30 ml. Samples for isotopic composition of sulphur were collected in 4.5 l plastic bottles and were acidified after filtration to pH 2. Then BaCl2 was added to precipitate BaSO4, which was then collected by filtering through a 0.45 µm fil- 160 ter. Sampling of water samples for 518O and 5D analyses was performed in plastic bottles of volume 50 ml filled to the top. Total alkalinity was measured by Gran titration with a precision of ± 1 % within 24 h of sample collection. Major elements (Ca2+, Mg2+, Na+, K+, Sr2+, Si, NO3-, SO42-, Cl-) in water samples were determined by Inductively Coupled Plasma-Optical Emission Spectroscopy (ICP-OES) and ion chromatography, with a precision of ± 2 %. Dissolved organic carbon (DOC) concentrations were measured using a high-temperature platinum-catalyzed carbon analyzer via combustion, followed by infrared detection of CO2 (Shmidizu TOC-5000A) with an uncertainty of ± 2 %. The concentrations of DIC were determined on a UIC Coulometrics CO2 Coulometer with a precision of ± 2 %. Isotopic composition of light elements (H, C, N, O, S) in sample is expressed as 5 value as relative difference of stable isotope ratios with regard to default reference material as follows: g^=vz------EM.looo y‰\ (1) RRM Value R is ratio between heavier stable isotope and lighter stable isotope: D/H, 13C/12C, 15N/14N, 18O/16O, 34S/32S. Positive 5 values (expressed in ‰) mean that sample contains more heavier isotope than reference material, and negative that it contains less (O’Neil, 1979). Isotopic reference materials are determined by International Atomic Energy Agency (IAEA) and National Institute of Standards and Technology (NIST). For each of the element the reference materials are selected so that isotopic ratio in reference material is similar to average ratio of the same isotopes in nature. At measurments of isotopic composition of oxygen (hydrogen), carbon, nitrogen and sulphur the following reference materials are used: • VSMOW - Vienna Standard Mean Ocean Water • VPDB - Vienna Pee Dee Belemnite • AIR - Air • VCDT - Vienna Canon Diablo Troilite The stable isotope composition of dissolved inorganic carbon (513CDIC) was determined with a Europa Scientific 20-20 continuous flow IRMS ANCA - TG preparation module. Phosphoric acid (100 %) was added Tja{a Kandu~ & Nives Ogrinc (100-200 fj.l) to a septum-sealed vial which was then purged with pure He. The water sample (6 ml) was injected into the septum tube and headspace CO2 was measured. In order to determine the optimal extraction procedure for surface water samples, a standard solution of Na2CO3 (Carlo Erba) with a known 513CDIC of - 10.8 ± 0.2 ‰ was prepared with a concentration of either 4.8 mmol/l (for samples with an alkalinity above 2 mmol/l) or of 2.4 mmol/l (for samples with alkalinity below 2 mmol/l). Since no reference material for 513CDIC measurements exists, CO2 was used with 513CCO2 - 5.2 ± 0.2 ‰ as a laboratory standard, which is calibrated to international reference materials: NIST 8562, NIST 8563 and NIST 8564 with default 5 value relative to VPDB. The isotopic composition of oxygen and hydrogen in water (518O and 5D) were measured with a Varian Mat 250. The same reference materials (NIST) were used for 518O and for 513CDIC measurements as described above. For measurement control of 5D value IAEA OH4 and IAEA OH1 reference materials were used with - 109 ‰ ± 1.5 and - 3.9‰ ± 1.5 to relate analytical results to VSMOW. The stable isotope composition of sulphur (534SSO4) was determined with a Europa Scientific 20-20 continuous flow IRMS ANCA - SL preparation module. About 10 mg of BaSO4 was scraped from the filters and transferred to a tin capsule. 534SSO4 was determined after combustion (1000 °C) of the capsule and reduction in a Cu tube (600 °C). NBS 127 and NBS 22 reference materials for sulphur were used to relate analytical results to the VCDT standard. Partial pressure of CO2 (pCO2), as well as saturation indexes of calcite and dolomite (SIcalcite, SIdolomite), were calculated by the PHREEQC for Windows program (P ar k -hurst & Appelo, 1999). Results and discussion Major chemistry The concentration of Ca2+ in Sava water seasonally changes from 39.9 to 60.8 mg/l in spring 2004, from 36.5 to 62.8 mg/l in late summer and from 39.2 to 72.8 mg/l in winter 2005. The concentration of Ca2+ in water of Sava tributaries seasonally ranges from 5.95 to 93.5 mg/l in spring 2004, from 11.9 to 95.4 mg/l in late summer 2004 and from 6.76 to 96.2 mg/l in winter 2005, respectively (Tables 1, 2 and 3). Hydrogeochemical characteristics of the River Sava watershed in Slovenia 161 The concentration of Mg2+ in Sava water changes from 8.0 to 12.8 mg/l in spring 2004, from 10.1 to 13.7 mg/l in late summer and from 11.8 to 16.7 mg/l in winter. The concentration of Mg2+ in water of Sava tributaries seasonally ranges from 1.84 to 22 mg/l in spring 2004, from 3.4 to 28 mg/l in late summer 2004 and from 1.9 to 29.9 mg/l in winter (Tables 1, 2 and 3). A Ca2+/Mg2+ molar ratio around 2 is characteristic of carbonate weathering (Mey-beck, 1996) and therefore for the carbonate part of the Sava watershed (sampling locations 1–5, upper part of the watershed). Deviation from this ratio arises from additional mineralization of water due to erosion of mechanically less resistant clastics rocks in the central and lower part of the River Sava flow, where the ratios are higher and also concentrations of Ca2+ increase in comparison to the Sava Bohinjka and Sava Dolinka, where weathering of carbonates prevails (Figure 2). The Ca2+/Mg2+ ratio also increases at the confluence of the Sava Dolinka and Bohinjka (9) in the spring sampling season due to snow melting and higher amounts of precipitation. The Ca2+/Mg2+ ratio also indicates the relative contributions of calcite and dolomite which contribute to the chemical composition of water. Dolomite weathering give a Ca2+/Mg2+ molar ratio around 1, while carbonate weathering gives a ratio around 2 (Pawellek et al., 2002). In the case of weathering of calcite as the dominant mineral, the Ca2+/Mg2+ ratio rises above 5. From the Ca2+/Mg2+ ratio in the Sava watershed it can be concluded that its chemical composition is mainly controlled by weathering of calcite with magnesium. Precipitation of tufa is a common process (Ramov{, 1983) at location 28 in the late summer season (low discharge conditions), which is reflected in the Ca2+/Mg2+ ratio. In the late summer sampling season the Ca2+/Mg2+ ratio is 2.7, while in spring it increases to 4.5. Saturation indexes for the Sava watershed were calculated by the PHREEQC program for Windows (Parkhurst & Ap-pelo, 1999). Thermodynamically calculated of saturation indexes of carbonates (SIcalcite, SIdolomite) were in the range of 0.5 from equilibrium (Appelo & Postma, 1994) and are presented for the River Sava system in Tables 1, 2 and 3. These saturation indexes indicate precipitation of carbonate in the central and lower Sava River flow, while in the upper flow the index rarely reaches satura- tion. The stable isotope composition of par-ticulate inorganic carbon (?3CPIC) measured from the River Sava watershed varied from 0.5 to 2.1‰ with an average of 1.3 ± 0.5 ‰ (n = 9), virtually identical in value that of carbonate rocks (Kandu~, 2006). This suggests that suspended carbonate material in these streams is detrital in origin and does not represent authigenic carbonate formed in response to the high degrees of supersaturation with respect to calcite. The partial pressure of CO2 in the Sava watershed changed from 158.5 ppm to 2951 ppm in spring 2004, from 446 to 9120 ppm in late summer and from 549 to 2290 ppm in winter (Tables 1, 2 and 3). The highest partial pressure was observed in the late summer season probably due to degradation of organic matter in soils leaching into the river and/or in the river itself (Dever et al., 1983). The higher partial pressure in the winter sampling season could also be attributed to higher CO2 dissolution in water in the winter months (Atkins, 1994). The River Sava in Slovenia shows seasonal patterns of CO2 flux between the water – atmosphere interface similar to other large rivers, like the Rhine (3300 ppm), Amazon (4000 ppm) and Yangtze (from 3160 to 4700 ppm) (Gao & Kempe, 1987). Alkalinity in the River Sava ranged from 2.61 to 3.75 mmol/l in spring 2004, from 2.63 to 4.76 mmol/l in late summer 2004 and from 2.67 to 4.22 mmol/l in winter 2005, respectively. Alkalinity ranged from 0.38 to 5.22 mmol/l in spring 2004, from 0.84 to 6.02 mmol/l in late summer 2004 and from 0.44 to 6.30 mmol/l in winter 2005 (Tables 1, 2 and 3). This seasonal variability of alkalinity is due to discharge conditions as well as processes in the River Sava watershed (Figures 3 A and 3 B). Lower alkalinities are observed in the spring sampling season due to dilution of river water by precipitation with an alkalinity around 0 mmol/l, while in late summer the higher alkalinity could be attributed to oxic degradation of organic matter. In winter 2005 discharges, degradation processes and leaching from the terrestrial to the river ecosystem are limited and consequently alkalinities are higher in comparison with the spring sampling season. In late summer higher alkalinities are observed due to lower discharge conditions and higher degradation of organic matter. Alkalinity is also related to the geological composition of the River Sava watershed. The upper part 162 Tja{a Kandu~ & Nives Ogrinc Table 1. Chemical and isotopic data for the River Sava watershed, spring 2004. Locations are plotted on Figure 1. Table 1 Sampling point Q (m3/s) T (°C) D. O. (%) Conductivity (nS/cm) p H Alkalinity (mmol/l) DIC (mmol/l) DOC Ca Mg Na K (mg/l) (mg/l) (mg/l) (mg/l) (mg/l) 1 8.6 88.0 276 7.83 2.98 2.98 3.83 40.75 12.60 1.23 2 1.55 8.9 96.0 244 8.14 2.62 2.54 2.64 3 7.6 94.5 216 8.36 2.37 2.29 3.11 32.94 10.64 1.59 4 11.4 9.7 99.3 251 8.28 2.81 2.69 4.07 39.96 12.81 1.60 5 28.6 9.7 98.6 283 8.48 3.19 2.81 2.76 45.09 10.40 2.18 6 n.a. 6.2 102.5 172 8.32 2.00 1.65 3.03 30.32 4.39 0.16 7 3.75 12.3 112.7 201 8.26 2.35 2.07 2.66 36.01 4.88 1.11 8 10.5 9.5 96.0 203 8.43 2.30 2.24 1.48 40.50 4.73 0.66 9 39.9 10.3 102.6 254 8.57 2.91 2.47 2.23 44.04 8.00 1.60 10 2.37 11.0 111.4 277 8.95 2.60 2.27 48.73 10.74 1.99 11 42.1 12.1 100.5 264 8.10 2.64 2.45 1.51 43.20 10.52 1.71 12 31.9 10.5 108.8 246 8.59 2.89 2.60 3.45 45.93 8.53 1.77 13 67.2 10.6 86.5 305 8.42 2.79 2.73 3.11 47.38 9.19 2.93 14 4.32 10.2 94.7 260 8.37 2.48 2.43 1.18 37.48 10.45 2.75 15 10.4 102.8 302 8.43 3.16 2.77 1.67 46.21 9.49 2.23 16 0.731 11.2 107.0 381 8.17 3.90 2.09 56.88 11.26 8.09 17 7.22 10.9 97.0 391 8.10 4.01 3.65 2.36 58.89 12.21 3.90 18 10.9 110.2 338 8.23 3.27 3.15 1.45 54.42 10.94 3.97 19 13.0 99.4 62 7.79 0.39 0.44 3.88 5.95 1.84 2.99 20 11.1 105.5 324 8.20 3.60 3.15 1.58 53.16 10.92 3.85 21 105 12.1 114.0 323 8.27 3.28 2.41 53.58 11.04 3.74 22 13.4 99.0 272 8.18 2.59 2.51 32.12 15.09 4.02 23 11.2 90.0 344 8.18 3.00 2.94 1.63 50.44 10.96 3.52 24 8.4 93.1 208 8.17 1.87 22.82 10.85 3.43 25 1 8.7 91.0 491 8.16 4.15 3.98 65.88 18.52 6.69 26 10.4 87.4 380 8.10 3.32 3.27 53.78 11.59 3.96 27 8.7 98.0 424 8.34 4.26 4.15 59.17 22.01 3.35 28 8.4 95.0 318 8.19 3.34 3.32 63.14 8.26 0.74 29 10.1 89.2 389 8.31 3.56 2.69 36.39 18.36 8.89 30 10.7 84.2 395 8.08 3.57 3.27 54.98 11.95 4.49 31 92.2 10.8 84.6 349 8.08 3.38 2.95 53.03 11.31 4.04 32 11.0 93.7 573 8.10 4.48 4.10 33 25 14.4 127.8 375 8.98 3.04 2.74 6.53 52.07 10.56 7.40 34 11.6 66.7 375 8.10 3.66 3.27 53.84 11.38 4.66 35 1.94 12.9 95.6 443 8.48 4.57 4.27 66.37 20.03 2.48 36 13.4 94.0 383 8.42 3.17 3.17 2.44 55.46 11.73 4.37 37 14.2 90.2 396 8.04 3.39 3.25 56.82 12.63 4.86 38 26.9 13.5 105.5 423 8.32 4.33 3.75 66.81 12.24 2.55 39 120 13.4 98.0 412 8.28 3.52 3.53 60.81 12.24 3.60 40 1.8 14.0 91.5 570 8.15 5.22 5.09 3.85 93.46 16.61 8.56 41 144 13.8 106.3 398 8.23 3.75 3.47 59.39 12.26 4.02 of the watershed is composed of carbonate rocks, where low alkalinities are observed, while in the central and lower flow higher alkalinities are observed since the watershed is composed of carbonate and clastic rocks. The annual flow of dissolved inorganic carbon (FDIC), calculated from carbonate alkalinity and discharges (EIONET, 2005) was estimated to be 4.1 × 1011 g C/year or 4.0 × 107 g C/(year km2) at the Bregana (41) location, according to the drainage area of the River Sava in Slovenia. In Tables 1, 2 and 3 results of other cation (Na+, K+, Sr2+ and Si) and anion (NO3-, SO42-and Cl-) analyses are also presented. Hydrogeochemical characteristics of the River Sava watershed in Slovenia 163 Na K (mg/l) (mg/l) Sr (mg/l) Si (mg/l) NO3 (mg/l) SO4 (mg/l) Cl (mg/l) 0.00 0.02 0.72 ^calcite ^dolomite pCO2 (ppm) 813Cdic (%o) s18o (%o) 5D (%o) d (%o) 834Ss04 (%o) 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.24 1.12 0.35 0.05 1.30 0.13 0.15 1.14 0.47 0.14 1.23 0.24 1.26 2.11 0.19 0.48 0.09 1.73 0.53 0.32 1.13 0.74 0.14 0.44 0.72 0.11 0.16 2.00 0.49 0.10 0.15 0.12 0.01 0.01 0.02 0.07 0.32 0.12 0.09 0.09 0.07 0.09 0.06 0.06 0.08 0.03 0.07 0.08 0.03 0.08 0.05 0.14 0.08 0.06 0.02 0.07 0.09 0.08 0.15 0.10 0.11 0.09 0.10 0.05 0.07 0.39 0.08 0.35 0.72 0.70 0.14 0.30 0.37 0.55 0.81 0.45 0.62 0.68 1.31 0.72 1.03 0.70 0.78 4.45 0.78 0.74 3.03 0.82 3.99 2.76 0.72 1.49 0.88 3.34 0.76 0.66 0.89 0.74 1.87 0.74 0.90 1.02 0.92 2.43 0.90 2.11 1.39 2.24 2.05 2.53 2.55 2.55 2.64 2.76 2.46 2.62 2.94 3.95 5.32 4.62 13.00 6.13 1.84 1.84 1.84 1.84 1.84 1.84 1.84 1.84 1.84 1.84 1.84 1.84 1.84 1.84 1.84 1.84 1.84 1.84 1.84 1.84 1.84 1.84 1.84 6.01 9.24 6.39 12.29 14.57 2.86 3.76 4.02 8.78 34.95 18.02 11.97 12.98 11.15 12.97 16.46 14.15 6.35 13.73 14.22 13.36 13.69 14.51 36.24 15.43 24.26 16.13 26.85 19.79 15.38 108.37 28.35 18.61 17.97 17.05 21.43 9.81 15.12 39.56 18.33 0.04 5.27 0.04 0.05 0.16 0.02 0.02 1.17 0.09 0.13 0.05 2.97 0.08 0.05 0.07 0.07 0.06 2.65 0.06 0.06 0.06 5.72 3.33 0.13 5.85 5.77 1.74 0.18 6.22 6.19 7.50 0.08 0.08 0.08 5.77 0.08 0.05 5.62 10.40 0.07 0.04 0.36 0.45 0.75 0.21 0.38 0.54 0.8 1.12 0.33 0.83 0.67 0.48 0.72 0.64 0.6 0.61 -1.56 0.61 0.66 0.29 0.5 -0.06 0.66 0.48 0.8 0.6 0.53 0.5 0.46 1.27 0.53 1.07 0.82 0.5 0.92 0.76 0.95 0.74 -0.34 0.31 0.54 0.98 -0.37 0.06 0.26 0.99 1.73 0.2 1.07 0.76 0.52 0.88 0.72 0.67 0.66 -3.46 0.68 0.79 0.43 0.48 -0.36 0.87 0.41 1.28 0.41 0.89 0.46 0.39 2.06 0.54 1.8 1.14 0.54 1.3 1.02 1.35 0.98 2138.0 158.5 489.8 707.9 501.2 457.1 645.7 407.4 371.5 128.8 1047.1 354.8 512.9 512.9 562.3 1288.2 1548.8 933.3 338.8 1096.5 871.0 891.3 977.2 616.6 1349.0 1288.2 912.0 1023.3 831.8 1445.4 1380.4 1819.7 141.3 1445.4 724.4 588.8 2951.2 1023.3 912.0 1819.7 1096.5 -10.7 -10.3 -7.9 -8.6 -9.6 -5.8 -8.3 -9.1 -9.0 -6.8 -7.0 -9.5 -10.2 -11.2 -10.7 -12.4 -13.5 -12.7 -8.1 -12.5 -12.3 -13.2 -12.0 -11.6 -13.0 -11.0 -12.4 -12.6 -10.5 -11.5 -11.5 -9.5 -8.5 -11.1 -11.5 -10.0 -11.9 -12.9 -11.1 -11.8 -11.1 -10.4 -10.2 -10.5 -10.5 -11.0 -12.4 -10.2 -10.6 -10.8 -10.5 -10.3 -10.6 -10.2 -9.3 -10.1 -9.2 -10.0 -9.9 -9.1 -9.8 -9.8 -9.4 -9.7 -9.6 -9.8 -10.4 -10.2 -10.2 -10.3 -10.0 -10.7 -10.1 -10.2 -10.2 -10.0 -10.6 -9.7 -10.1 -10.5 -10.3 -10.2 -65.5 -64.1 -67.1 -67.9 -65.5 -80.0 -62.8 -71.5 -66.0 -66.0 -66.7 -66.4 -69.7 -61.5 -67.0 -60.9 -66.3 -68.0 -57.6 -63.3 -61.2 -61.0 -68.0 -63.4 16.66 16.48 15.85 15.05 21.4 17.96 17.78 12.24 19.32 16.95 14.67 17.34 10.88 11.97 12.79 11.78 12.7 10.21 14.29 13.33 14.64 16.42 14.16 17.18 9.5 14.5 9.4 10.4 9.6 3.9 1.8 7.8 -72.3 9.07 -70.5 -62.6 -67.0 -64.7 -65.0 -65.4 -55.7 -69.9 -64.8 -66.5 -68.8 14.03 17.19 13.58 15.88 14 18.34 20.93 9.89 18.15 14.87 11.78 7.1 6.1 8.5 5.5 4.2 8.5 4.3 3.7 6.3 3.9 3.8 -1.3 5.4 Carbon Isotopes Values of ?13CDIC in the River Sava ranged from – 12.7 to – 8.6 ‰ in spring 2004, from – 11.8 to – 7.3 ‰ in late summer 2004 and from – 10.6 to – 6.3 ‰ in winter 2006 (Tables 1, 2 and 3, Figure 3 C). In the Sava tributaries ?13CDIC values ranged from – 13.5 to – 5.8 ‰ in spring 2004, from – 12.8 to – 3.3 ‰ in late summer 2004 and from – 11.9 to – 4.2 ‰ in winter 2005 (Tables 1, 2 and 3, Figure 3 D). Lower ?13CDIC values in the Sava and its tributaries were observed in the spring season due to more abundant precipitation. More positive ?13CDIC values were observed in all sampling seasons in the carbonate 164 Tja{a Kandu~ & Nives Ogrinc Table 2 Table 2. Chemical and isotopic data for the River Sava watershed, late summer 2004. Locations are plotted on Figure 1. Sampling Q T D.O. Conductivity Alkalinity DOC Ca Mg Na K point (m3/s) (°C) (%) (jiS/cm) pH (mmo/l) (mg/l) (mg/l) (mg/l) (mg/l) (mg/l) 1 5.6 110.0 279 7.54 2.63 9.17 36.57 11.40 1.55 2 1.55 7.1 103.0 380 7.56 2.67 5.94 37.45 11.51 2.95 3 7.6 114.0 239 7.82 3.36 5.50 34.36 10.75 1.04 4 11.4 8.5 110.0 286 7.86 2.65 7.13 40.04 12.32 1.62 5 28.6 10.7 118.0 311 8.31 3.06 10.75 44.88 12.28 2.35 6 n.a. 6.0 116.0 190 8.08 2.19 13.72 30.14 5.49 0.21 7 3.75 9.5 104.0 369 7.85 4.11 7.17 93.62 12.41 6.39 8 10.5 12.0 118.0 247 8.11 2.69 5.23 43.15 5.80 0.83 9 39.9 12.2 112.0 305 8.26 3.24 6.99 47.22 11.33 2.46 10 2.37 11.3 107.0 378 8.44 2.81 4.89 55.22 14.50 2.43 11 42.1 10.7 109.0 333 8.20 3.22 9.71 49.24 13.05 2.15 12 31.9 14.7 107.0 632 7.57 3.22 6.57 55.21 12.55 3.88 13 67.2 12.8 109.0 307 8.17 4.76 9.66 48.82 10.14 2.61 14 4.32 12.5 104.0 379 7.48 4.51 5.48 55.20 15.39 4.45 15 13.2 102.0 379 7.32 2.99 6.26 16 0.731 13.8 99.0 554 7.74 4.6 6.35 68.81 14.19 16.95 17 7.22 15.7 96.0 500 7.93 4.79 13.85 68.95 18.03 11.87 18 14.0 106.0 366 8.08 3.44 6.92 53.78 12.04 5.20 19 11.9 100.0 118.5 7.24 0.84 11.98 3.39 3.72 20 11.5 100.0 387 7.29 3.58 5.99 58.14 12.01 4.03 21 105 12.1 100.0 393 7.81 3.48 6.42 59.23 12.03 5.22 22 13.6 96.0 416 7.89 3.93 6.53 23 13.6 106.0 378 8.32 3.81 6.53 57.56 11.84 3.99 24 12.7 108.0 316 8.34 2.59 12.38 25 1 18.5 99.0 574 8.64 4.89 3.58 67.73 24.30 9.35 26 13.7 95.0 403 8.35 3.63 59.77 12.53 4.10 27 14.7 85.0 487 8.74 5.15 5.11 60.90 2.72 28 15.0 95.0 414 8.72 4.03 3.01 61.80 13.44 1.44 29 12.5 94.0 512 8.04 4.22 18.66 56.42 24.22 16.40 30 17.2 96.0 403 8.48 3.54 4.09 59.39 12.73 4.95 31 92.2 13.5 100.0 376 8.52 3.45 3.20 57.34 12.34 4.48 32 17.3 64.0 575 8.34 3.42 11.10 66.85 24.65 17.68 33 25 14.3 473 8.82 3.42 64.97 13.70 15.71 34 14.4 397 8.54 3.36 3.61 58.39 12.94 7.36 35 1.94 14.4 511 8.99 5.49 2.95 69.46 28.03 5.31 36 14.6 453 8.63 3.6 5.00 37 14.4 428 7.86 3.29 4.86 62.83 13.67 8.46 38 26.9 14.5 459 8.77 4.65 3.47 39 120 14.7 412 8.20 3.74 4.50 59.63 13.07 7.91 40 1.8 12.4 630 8.61 6.02 7.67 95.43 20.29 15.36 41 144 14.0 88.0 424 7.93 3.51 6.18 61.15 13.46 6.89 part of the watershed, while in the central and lower Sava flow, which is composed of carbonate and clastic rocks, more negative ?13CDIC values were observed. In the carbonate part of the Sava watershed lower al-kalinities and lower saturation indexes were determined in comparison with the central and lower flow (Tables 1, 2 and 3). Also a negative trend of ?13CDIC values versus alkalinity was observed in all sampling seasons (Kandu~, 2006). The isotopic composition of the source material should be studied before making an interpretation of the origin of DIC in a river ecosystem. Plants represent the source of organic material in rivers, while carbo- Hydrogeochemical characteristics of the River Sava watershed in Slovenia 165 Na K Sr Si N03" S04 Cl PC02 813Cdic 8180 8D d mg/l) (mg/l) (mg/l) (mg/l) (mg/l) (mg/l) (mg/l) ^'calcite ^'dolomite (ppm) (%o) (%o) (%o) (%o) 0.47 0.03 0.62 2.10 4.58 2.52 -0.39 -1.25 3630.8 -7.8 -9.9 -66.5 11.7 0.06 0.07 1.27 1.82 9.91 10.54 -0.34 -1.12 3548.1 -8.6 -9.8 -67.5 9.9 0.78 0.09 0.36 1.99 5.24 1.55 -0.01 -0.44 2454.7 -5.4 -10.4 -72.0 10.2 0.81 0.12 0.71 3.54 38.87 1.95 -0.04 -0.5 1737.8 -7.3 -10.1 -66.6 13.2 0.29 0.15 0.74 3.70 18.45 4.39 0.58 0.73 724.4 -7.3 -9.8 -67.1 10.3 0.00 0.01 0.20 1.82 2.79 0.47 0.01 -0.68 871.0 -3.3 -8.8 -58.0 11.5 0.81 0.05 1.09 5.39 5.86 8.42 0.52 0.28 2754.2 -12.8 -9.2 -60.6 12.1 0.20 0.02 0.48 2.81 7.28 2.91 0.35 -0.01 1047.1 -7.2 -9.1 -58.9 13.0 1.29 0.11 0.80 2.88 14.23 3.55 0.6 0.74 871.0 -7.3 -9.6 -67.5 8.3 0.35 0.51 1.45 2.49 55.49 2.06 0.73 1.04 489.8 -6.4 -9.7 -66.5 10.1 0.47 0.16 1.28 3.00 20.14 2.25 0.53 0.61 977.2 -7.5 -9.3 -63.5 10.0 0.56 0.14 1.22 6.01 17.17 6.77 0.02 -0.41 4466.8 -8.8 -9.3 -68.6 4.9 0.62 0.11 0.86 3.66 12.50 3.57 0.69 0.87 1584.9 -8.1 -9.4 -68.9 5.4 1.88 0.11 2.03 9.01 14.84 6.75 0.03 -0.33 7413.1 -10.5 -8.8 -57.8 11.7 4.54 12.07 3.73 7762.5 -8.6 -9.1 -60.4 11.5 6.46 0.07 1.88 42.60 17.23 18.11 0.39 0.29 4168.7 -9.3 -9.0 -58.2 12.9 2.15 0.11 1.02 6.49 13.70 16.48 0.62 0.88 2884.0 -11.8 -8.7 -58.3 10.4 1.51 0.10 1.10 7.58 14.81 6.68 0.52 0.59 1445.4 -9.9 -8.9 -61.1 9.2 0.98 0.06 5.46 1.78 7.35 3.08 -1.51 -3.42 2511.9 -9.0 -61.0 1.70 0.09 1.48 6.25 14.00 4.97 -0.25 -1.03 9120.1 -11.8 -9.0 -65.3 5.8 1.78 0.09 1.53 7.45 6.52 14.74 14.89 6.43 10.64 0.27 0.01 2691.5 2691.5 -10.9 -11.6 -8.7 -59.2 -59.5 9.5 1.08 0.09 1.48 5.94 2.03 13.70 16.84 5.07 3.35 0.81 1.13 891.3 631.0 -10.2 -10.4 -9.0 -9.3 -66.4 4.7 3.67 0.13 2.97 7.31 29.78 9.21 1.31 2.44 537.0 -11.0 -9.4 -63.4 10.9 0.84 0.09 1.55 5.97 14.47 5.52 0.92 1.35 977.2 -9.8 -64.7 0.00 0.07 1.57 3.04 22.91 4.01 1.33 2.53 426.6 -10.1 -9.8 -65.2 12.2 1.02 0.03 1.07 4.46 21.70 2.37 1.24 2.04 363.1 -10.9 -10.1 4.04 0.10 4.16 3.66 12.50 3.57 0.55 0.9 1905.5 -11.4 -9.0 -58.0 13.1 1.57 0.10 1.51 6.68 16.85 5.68 0.99 1.56 588.8 -10.9 -9.1 1.44 0.09 1.39 10.32 14.59 5.29 0.96 1.44 501.2 -10.2 -9.1 -64.2 7.7 4.72 0.15 3.78 6.32 42.22 14.63 0.87 1.55 776.2 -11.6 63.6 3.06 0.20 2.83 7.58 48.86 11.90 1.25 2.04 234.4 -10.2 -8.8 -63.7 5.8 1.53 0.12 1.62 6.16 22.84 7.50 0.98 1.51 467.7 -9.9 -8.8 -57.0 12.5 3.48 0.12 2.60 3.27 17.55 4.92 1.59 3.01 234.4 -11.2 -9.1 -62.1 9.8 7.20 23.97 8.11 446.7 -10.4 -8.7 -58.0 10.7 2.29 0.11 1.84 6.79 27.25 8.66 0.35 0.23 2290.9 -10.8 -8.8 -55.7 13.8 7.27 9.78 5.66 407.4 -11.1 -9.3 -64.0 9.5 1.91 0.11 1.79 6.31 25.83 8.27 0.71 0.97 1174.9 -11.5 -8.8 -59.4 10.1 6.77 0.40 4.52 8.34 37.39 13.90 1.41 2.33 660.7 -12.2 -8.5 -60.2 7.0 1.82 0.10 1.67 2.49 21.40 2.06 0.42 0.37 2089.3 -10.8 -8.9 -59.4 10.9 nates represent the source of inorganic material. Plants in the Sava watershed have an average value of – 31.6 ‰ ± 1.5 (n = 22), while Mesozoic carbonates have an average value of 1.4 ‰ ± 1.3 (n = 13). The major inputs to the DIC flux (DICout) and ?13CDIC originate from tributaries (DICtri), degradation of organic matter (DICorg), exchange with the atmosphere (DICex), and dissolution of carbonates (DICca): DICout = DICtri – DICex + DICorg + DICca (2) = DIC · ?13C DIC0Ut • 513C0Ut • 513C + DIC DIC ex • S13C (3) 166 Tja{a Kandu~ & Nives Ogrinc Table 3. Chemical and isotopic data for the River Sava watershed, winter 2005. Table 3 Locations are plotted on Figure 1. Sampling Q T Conductivity Alkalinity DOC Ca Mg Na K point (m3/s) (°C) (p.S/cm) pH (mmo/l) (mg/l) (mg/l) (mg/l) (mg/l) (mg/l) Sr (mg 1 5.3 257 8.06 2.73 39.21 11.78 1.21 0.70 2 0.66 2.3 286 8.25 2.67 0.84 41.78 12.39 2.33 0.58 3 2.3 285 8.29 2.35 1.26 35.54 10.92 0.85 1.20 4 6.66 5.1 288 8.34 2.88 1.04 43.80 12.94 1.24 0.96 5 14.90 3.3 327 8.43 3.14 0.77 51.54 14.13 2.69 0.29 6 3 187 8.51 1.81 7 1.40 3.3 187 8.25 2.12 0.93 8 2.95 3.2 260 8.75 2.73 48.20 6.91 0.61 0.43 9 20.70 3.1 391 8.21 3.42 1.03 56.68 13.53 2.57 3.88 10 2.84 3 385 8.42 2.48 1.02 69.20 26.57 3.13 0.53 11 1.69 2.4 337 8.51 3.18 3.07 57.42 14.13 1.62 3.98 12 9.69 4.2 440 8.02 3.64 0.24 62.17 14.38 8.06 1.83 13 92.30 4.4 373 8.1 3.51 0.54 64.08 15.91 5.72 1.96 14 7.03 3.6 334 8.22 3.11 0.77 51.99 14.45 3.00 0.92 15 4.5 373 8.26 3.44 58.19 14.45 3.41 0.88 16 0.31 7.7 565 7.76 4.89 1.66 79.33 16.64 14.17 4.31 17 18.10 5.5 401 8.14 3.95 0.79 70.43 15.66 2.21 0.48 18 5.2 415 8.12 3.82 0.89 62.78 15.01 4.30 0.52 19 1.6 7.85 0.44 1.43 6.76 1.89 2.56 0.80 20 3.8 414 8.25 21 4.3 441 8.31 4.00 1.64 66.01 14.92 4.32 1.13 22 2.4 8 3.71 0.11 46.76 23.54 2.92 0.52 23 1.9 423 7.98 3.93 1.44 64.11 15.48 5.03 2.64 24 2.1 8.1 2.24 0.48 28.12 14.37 2.46 2.33 25 1.4 8.17 4.62 0.24 69.74 23.65 6.07 1.47 26 2.8 7.88 3.58 0.13 63.57 15.82 5.33 3.64 27 1.6 8.38 4.71 0.86 62.66 27.89 1.94 2.87 28 5.5 370 8.28 3.58 2.15 63.04 10.96 0.41 0.78 29 4.4 928 7.87 6.30 30 3.2 432 7.94 4.22 0.60 64.68 15.73 5.35 1.44 31 62.20 3.9 421 7.97 4.17 0.40 69.53 9.35 0.47 2.31 32 3.2 544 8.1 4.89 1.55 63.67 24.39 12.96 3.88 33 12.40 1.3 481 8.54 3.80 72.08 14.58 13.82 3.44 34 3.4 444 8.14 3.88 0.42 60.01 14.06 3.71 1.10 35 1.29 5.10 0.00 36 3.2 424 8.2 3.97 0.53 65.21 15.63 5.96 1.30 37 6.2 467 8.08 3.96 0.35 72.76 16.71 6.67 1.49 38 12.40 2.8 469 8.44 4.33 39 96.30 5.4 445 8.27 3.54 0.55 67.21 16.11 6.71 1.64 40 1.67 0.4 606 8.22 5.91 0.84 96.95 22.18 12.49 4.14 41 128.00 4.9 445 8.29 3.49 0.53 69.17 16.57 5.79 1.83 The process of photosynthesis was considered insignificant and therefore excluded from the mass balance calculations. Dissolved oxygen saturation reached a maximum of only 118 % in the late summer sampling and the ?13CPOC, with an average value of – 26.7‰, indicates that organic matter derived from aquatic photosynthesis is re- latively insignificant in the River Sava system. According to the concentration and isotope mass balance calculations (equations 2 and 3) of dissolved inorganic carbon, the proportions of the four processes (1) effect of tributaries (DICtrb), (2) degradation of organic matter (DICorg), (3) dissolution of Hydrogeochemical characteristics of the River Sava watershed in Slovenia 167 K Si N03 S04 Cl PC02 813Cdic 8D d (mg/l) Sr (mg/l) (mg/l) (mg/l) (mg/l) (mg/l) ^'calcite ^'dolomite (ppm) (%o) 8180 (%o) (%o) (%o) 0.00 0.25 3.44 4.34 2.78 0.17 -0.16 1174.9 -9.3 -9.1 -70 2.3 0.06 0.65 3.40 11.59 5.52 0.32 0.07 707.9 -8.2 -10.0 0.08 0.17 1.77 7.42 2.21 0.24 -0.07 575.4 -7.3 0.11 0.34 2.22 14.56 3.06 0.50 0.49 631.0 -7.5 -10.0 -69 9.7 0.16 0.48 3.00 20.81 5.75 0.65 0.73 549.5 -6.3 -9.8 -68 9.3 1.70 2.53 0.45 263.0 -4.2 -9.6 -59 17.1 2.18 2.91 0.96 575.4 -7.5 -9.6 -64 12.3 0.01 0.25 2.86 4.56 1.92 0.87 0.90 223.9 -7.3 -8.8 -59 10.3 0.09 0.50 6.14 16.30 9.68 0.51 0.38 1000.0 -8.9 -9.1 -63 9.1 0.09 0.90 3.38 52.98 5.01 446.7 -5.8 -9.6 -67 9.3 0.45 0.60 1.77 12.90 1.31 457.1 -6.4 -9.1 -61 10.8 0.14 0.61 8.31 23.06 12.58 0.40 0.17 1659.6 -9.4 -9.3 -63 10.5 0.10 0.53 7.77 20.96 6.80 1349.0 -8.7 -9.1 -61 10.8 0.07 1.03 8.21 14.04 4.77 0.45 0.35 891.3 -9.7 -8.9 -59 11.4 0.12 0.67 8.96 21.57 7.25 0.59 0.59 912.0 -8.5 -9.2 -63 10.0 0.06 0.93 39.21 19.23 19.94 0.41 0.22 4265.8 -10.9 0.06 0.58 5.22 13.15 4.68 1380.4 -11.9 -8.6 -60 7.5 0.09 0.68 8.21 16.94 7.00 0.54 0.48 1412.5 -10.0 -8.9 -59 10.9 0.01 2.32 2.52 6.92 2.61 -1.57 -3.74 302.0 n.a. -8.8 -61 8.3 0.09 0.65 933.3 -9.7 -8.8 -61 8.7 0.01 1.23 4.29 13.28 4.59 0.24 0.14 1737.8 -10.7 -8.7 -63 5.7 0.10 0.50 9.81 19.39 10.30 0.37 0.07 1949.8 -10.0 -8.3 -60 5.3 0.03 1.73 2.17 16.47 3.71 -0.07 -0.47 851.1 -10.4 -8.9 -57 13.7 0.11 1.01 5.06 30.37 7.86 0.63 0.74 1445.4 -10.2 -9.1 -63 9.1 0.10 0.54 8.29 19.76 10.54 0.24 -0.14 2238.7 -9.1 -8.9 -61 9.5 0.04 0.57 3.73 25.96 6.24 0.80 1.20 912.0 -9.0 -9.3 -66 8.0 0.00 0.45 4.45 18.66 3.16 0.68 0.63 912.0 4073.8 -10.6 -11.0 -9.8 -9.1 -62 -63 15.4 8.9 0.11 0.54 9.13 21.59 8.98 0.38 0.13 2290.9 -9.5 -8.8 -56 13.1 0.00 0.44 8.11 19.97 8.44 2138.0 -9.5 -8.9 -63 7.2 0.10 1.34 9.05 40.60 15.85 0.58 0.73 1819.7 -8.5 -9.4 -62 12.3 0.18 1.02 10.70 44.57 14.08 0.93 1.10 501.2 -8.8 -9.1 -62 9.6 0.13 0.61 9.08 23.00 9.65 0.52 0.40 1349.0 -10.5 -8.8 -61 9.0 0.00 0.00 3.83 18.00 5.63 -10.0 -9.2 -62 11.0 0.11 0.65 8.54 21.66 8.99 0.61 0.59 1174.9 -9.2 -8.9 -62 8.1 0.12 0.76 9.10 25.63 10.18 0.59 0.58 1621.8 -10.6 -8.9 -60 10.1 8.09 10.07 6.33 741.3 -10.3 -9.8 -65 12.5 0.11 0.63 10.40 23.60 9.97 912.0 -10.4 -8.8 -61 8.2 0.37 1.28 5.55 38.90 16.04 0.75 0.93 1621.8 -10.9 -9.4 -67 7.1 0.10 0.62 5.47 15.15 6.68 1096.5 -9.8 -9.2 -60 12.8 carbonates (DICcarb) and (4) equilibration of atmospheric CO2 (DICex) that influence ?13CDIC values in the river were calculated. The calculated proportions of each process, i. e. DICtrb : DICorg : DICex : DICcarb vary seasonally and at Bregana (41) are: 70 : 11 : – 1:19 % in spring 2004, 62 : 17 : – 5:26 % in late summer 2004 and 70 : 11 : – 4 : 23 % in winter 2005. Among the biogeochemical processes the dissolution of carbonates contributes a higher proportion than degradation of organic matter and is reflected in all sampling seasons (Kandu~, 2006, Kandu~ et al. in press). 168 Tja{a Kandu~ & Nives Ogrinc RIVER SAVA A 5 .5 5 .0 4 - 4 2 .5 2 .0 1 .5 1.0 0.5 6.0 5.0 4.5 4.0 spring 2004 late summer 2004 winter 2005 RIVER SAVA A CARBONATE DISSOLUTION 41 DOLOMITE DISSOLUTION Distance from the source (km) 20 40 60 80 100 120 140 Distance from the source (km) 160 180 200 220 0.5 DDDDDDDDDD T--------------------------------------1-------------------------------------- 20 40 60 80 100 120 140 Distance from the source (km) 160 180 200 220 uFirgeur2e 2. Ca2+/Mg2+ ratio versus distance to the source (A – River Sava, B – River Sava tributaries). Distance from the source (km) Oxygen and Hydrogen Isotopes FigDuirsech2arges at the sampling locations (Tables 1, 2 and 3) in the River Sava are dependent on precipitation, evaporation, eva-potranspiration, infiltration, equilibration with run-off and anthropogenic engineering structures like dams (Yee et al., 1990). The contributions of tributary discharges to the River Sava were estimated at locations where discharge data exists (Tables 1, 2 and 3). At the confluence of the Sava Bohinjka and Sava Dolinka (sampling point 9, Sava Oto~e) the Sava Bohinjka contributes 64 %, while the Sava Dolinka represents 55 % of the total discharge at the confluence in the spring sampling season. In the late summer and winter sampling season the Sava Bohinjka contributes 26 and 14 % of total discharge, while the Sava Dolinka contributes 72 %. The tributaries have higher contributions in the spring sampling season due to snow melting, the higher amount of precipitation and consequent surface and subsurface leaching from the slopes to the tributaries and finally to rivers. Those Sava 0 0 era MC < rt: 5> w 3 er i" SL wB' wf re a W^ EmLL f= e 3.< ro et t/3 ^ 7q -;: spring 2004 late summer 2004 winter 2005 18 iBON RIVER SAVA A 5-1 4- |2 1 4 RBO a 12, 13 CAF 21 23 A^134 ™ 3g 20 *-A26 / ^A._______A-----A™ V 37 " 4 ATE AND CLASTIC ROCKS I3" < 2-1 -0- NATE ROCKS 60 80 100 120 140 160 18 Distance from the source (km) Wit. ib ^ < !! IE ^ Ö 3d ^ !! tc e ¦? C/3 IB spring 2004 late summer 2004 winter 2005 RIVER SAVA C -4- -6- 4 -8- '*/ \9\12 » 36 10- ¦ 2 13^ 1 5 \ . /AA 34 /*\ 41 39A 12- 18 20 _¦ ,. 31 r^23 30 21 37 14- 100 120 140 160 180 200 220 Distance from the source (km) .-spring2oo4 R|VER SAVA TRIBUTARIES B •—late summer 2004 A-^/inter 2005 CARBONATE CARBONATE AND CLASTIC ROCKS CARBONATE AND CLASTIC ROCKS -4--5 -10--11 - -13--14- 60 80 100 120 140 160 Distance from the source (km) RIVER SAVA TRIBUTARIES 60 80 100 120 140 160 Distance from the source (km) rt SB rt m CB < CB "S CB w < re 170 Tja{a Kandu~ & Nives Ogrinc tributaries with higher drainage areas, namely the Ljubljanica (17), Savinja (33) and Krka (38), have contributions of 42 %, 14 % and 18 % in the spring sampling season, 7 %, 27 % and 22 % in late summer and 33 %, 20 % and 13 % in winter, respectively. ?18Owater values in the River Sava seasonally changed from – 11.0 to – 9.7 ‰ in spring 2004, from – 10.1 to – 8.9 ‰ in late summer 2004 and from – 10.0 to – 8.8 ‰ in winter 2006, respectively. Values of ?18Owater of Sava tributaries varyed from – 12.4 to – 9.1 ‰ in spring 2004, from – 10.4 to – 8.7 ‰ in late summer 2004 and from – 9.8 to – 8.6 ‰ in winter month (Tables 1, 2 and 3, Figure 4 A and 4 B). Lower ?18Owater was observed in the Sava watershed in spring 2004 while in the late summer and winter seasons lower ?18O values were observed due to the high amount of precipitation and water leached from slopes forming the watershed (Vre~a et al., 2004). In the winter sampling season more negative ?18O values would be expected due to the isotopically lighter precipitation, but in fact the values were similar to the late summer sampling season. Seasonal variation of ?18O values in the Sava watershed results from differences in water temperature and the height of the recharge area between sampling points in the watershed. The River water temperature is higher in all three seasons in the central and lower part of the River Sava flow in comparison with the upper part of the flow (sampling point 9). The height of the recharge area in the upper flow (sampling point 5) is above 1500 m a. s. l. in comparison to the central and lower part of the flow, where the height above sea level is up to 1000 m a. s. l. The longer residence time in dams could affect the evaporation process and enrichment with heavier oxygen isotope in the summer months in the lower part of the Sava flow (location 36, see Figure 4 A and 4 B). Kinetic effects during evaporation may be attributed to surface water temperature, wind velocity (shear at the surface water) and, most importantly the relative humidity of the air (molecular diffusion between water-atmosphere) (Gonfiantini, 1986). Evaporation between the River Sava and the atmosphere in the watershed is also dependent on the air temperature, which varied from 8.6 to 14.2 °C in spring 2004, from 10 to 19 °C in late summer 2004 and from – 4.9 to 2.2 °C in winter (EIONET, 2005). In the Sava watershed fog phenomena (condensation), which are more pronounced in late summer, could also daily affect ?18Owater and ?Dwater values (Petkov{ek, 1969). In the spring sampling season the ?18O value was lower above the confluence of the Sava Dolinka and Sava Bohinjka (sampling location 5) than below the confluence (sampling location 9) since the Sava Bohinjka (sampling location 8) contributes 64 % of the discharge, with higher ?18Owater than the Sava Dolinka. A higher deviation in ?18Owater values between the seasons was observed at sampling location 6 (Sava Bohinjka source). In spring the lower ?18Owater value is due to snow melting, isotopically lighter rain and the height of the recharge area, which is around 1800 m a. s. l. (Urbanc & Brancelj, 2002). The higher ?18Owater value at location 6 in comparison to other locations could be due to isotopically heavier precipitation in the late summer months since generation of water vapour and precipitation is a complicated process, which includes kinetic evaporation and mixing (Clark & Fritz, 1997). The lower ?18O at sampling location (7), which reflects mixing of water and bio-geochemical processes in lake Bohinj in the late summer season in comparison with sampling location (8), could be attributed to mixing of waters in the lake after stratification (mixing of cold and warm water) and/or anthropogenic influence. In the late summer season at sampling location (7) higher concentrations of other parameters were also observed (Ca2+, Mg2+, HCO3-, Na+, K+, SO42-and NO3-). ?Dwater values in the Sava River watershed ranged from – 80.0 to – 55.7 ‰ in spring 2004, from – 72.0 to – 55.7 ‰ in late summer 2004 and from – 69.6 ‰ to – 56.4 ‰ in winter 2005, respectively. The correlation between ?18Owater and ?Dwater for most meteorological and surface waters which do not reflect evaporation is expressed by the GMWL (Global Meteoric Water Line), which is determined by the equation: ?D = 8 · ?18O + 10 and is valid for a humidity of more than 85 % (Clark and Fritz, 1997). In precipitation from Ljubljana the local meteorological line (LMWL) for the period from 2001 to 2003 was expressed by ?D = 7.9 ?18O + 8.3 (Vre~a et al., 2004). Sampling points deviating from the LMWL in the Sava watershed are a result of kinetic evaporation/condensation processes and are dependent on conditions above the surface water (temperature, Hydrogeochemical characteristics of the River Sava watershed in Slovenia 171 - -11 -12 -13 20 40 60 Distance from the source (km) 80 100 120 140 160 180 200 220 -7 -8 - -9 -, -7 -8 -9 -10 -11 -12 - -13 20 40 60 Distance from the source (km) 80 100 120 140 160 180 200 220 Figure 4. ?18O versus distance to the source of the River Sava (A – River Sava, B – River Sava tributaries). relative humidity). Previous investigations of the isotopic composition of precipitation collected in the Mediterranean basin indicate a deuterium excess (d) for the western part of the Mediterranean basin of 14 ‰, while for eastern part it was 22 ‰ (IAEA, 2005). It was found that at locations where the value of d in precipitation is below 10 ‰ (Dansgaard, 1964) this is due to secondary processes (evaporation in periods of low relative humidity in the atmosphere), while values above 10 ‰ belong to precipitation from the Mediterranean Sea (Cruz -San et al., 1992). In the Sava watershed d values changed from 9.1 to 21.4 ‰ in spring 2004, from 4.7 to 13.8 ‰ in late summer 2004 and from 2.3 to 17.1 ‰ in winter 2005 (Tables 1, 2 and 3). 0 0 172 Tja{a Kandu~ & Nives Ogrinc Sulphur isotopes The concentration of sulphate in the Sava watershed changed from 6 to 39.6 mg/l in spring 2004, from 4.6 to 55.5 mg/l in late summer and from 4.3 to 25.6 mg/l in winter (Tables 1, 2 and 3, Figure 6). The highest concentrations were observed in the periods of low discharges (late summer 2004) while the opposite trend is observed in spring (Figure 6). The lowest sulphate concentration was observed at the Sava Dolinka spring (sampling point 1, Figure 1) during different sampling seasons. The highest sulphate concentrations were observed in all sampling seasons at the locations Trbovelj{~ica (29), Savinja (33) and Sotla (40) due to anthropogenic influences, arising from coal mining waste water at location 29, industry with sulphuric acid production at location 33 and spraying vinyards with Bordeaux mixture at location 40. High sulphate concentrations were also observed at Trži{ka Bistrica (10) due to dissolution of Palaeozoic carbonates with evaporates which form the watershed. Furthermore, higher concentrations of sulphate were observed in the River Sava at sampling locations 4 and 26. At all sampling locations in the Sava River watershed sulphate concentrations are below 250 mg/l, which represents the limiting concentration for tap water. Values of ?34S in the Sava watershed (Figure 7) which were sampled in spring 2004 varied from – 1.3 ‰ (sampling point 40) to 14.9 ‰ (sampling point 10). More positive values (14.9 ‰, sampling point 10) indicate dissolution of evaporates (Krouse, 1980). The ?34S depleted signature of the River Sotla suggests that its SO42- was derived from the oxidation of sulphide minerals. Oxidation of sulphide minerals composing sedimentary clastic rocks could contribute higher sulphate concentrations than atmospheric sources (precipitation) and result in ?34SSO4 values up to 10 ‰ (Clark & Fritz, 1997). Figure 7 also shows that most of the samples did not originate from precipitation (except sampling point 1) and from soil (except locations 11, 16, 36 and 39) due to their higher concentrations of sulphate. Precipitation and soil water are not the sole sources contributing to sulphate concentrations in Sava water samples. Assuming that the sources of SO42- to the Sava River are from tributaries (Ftri), precipitation (Fp) and other sources (Fother), the contributions of these inputs can be quantified by the following steady state equations: F = F + F + F SR tri p other 534S Fsr = 534Stri Ftr + 834Sother Fother + 534Sp Fp + (4) (5) -30 -40 -50 -60 -70 -80 -90 -100 s /• ' ^t ' s<' GMWL s ' /* 5D =8-818O + 10 *><_. *z* y "~~^L s MML x' sT 5D =8×818O + 22 X-"A *%&S^ ^ * 8D =7.9-818O + 8.3 River Sava spring 2004 ¦ River Sava late summer 2004 River Sava tributaries spring 2004 ^ ^*-^ ^ ^^ 0 River Sava tributaries late summer 2004 S A River Sava winter 2005 s A River Sava tributaries winter 2005 -14 -13 -12 -11 -10 -9 8180 (%») -8 -7 -6 -5 -4 Figure 5. 8D vs. 818O for the River Sava system in Slovenia. GMWL - Global meteorological water line, LMWL - Local meteorological water line, MML - Mediterranean meteorological water line. Hydrogeochemical characteristics of the River Sava watershed in Slovenia 173 70 60 50 40 30 20 10 RIVER SAVA A spring 2004 -¦-late summer 2004 4 -A-winter 2005 26 - / ¦ \ 5 12 15 ^^9L-^~^^\ 13 18 23/ ' 20 21 / * 34 . 30 ,^fc-= •so 3^ 39 36 J-7» ______*--^>^™5\ 41 31 ~" ^* 72 11 0 20 40 60 80 100 120 140 160 180 200 220 Distance from the source (km) 60 50 - 40 30 - 20 10 H 3 A RIVER SAVA TRIBUTARIES 33 20 40 60 80 100 120 140 160 180 200 220 Distance from the source (km) rFeig 6ur e 6. Sulphate concentrations versus dtriisbtauntacreietos) t.he source (A – River Sava, B – River Sava Where: Fi = Qi • [S042"]j (mol/s) (6) i = RS – River Sava, tri – tributaries, p – precipitation, other – other sources FRS (mol/s) – is the annual dissolved sulphate flow at the Bregana sampling location and contributions affecting flow of dissolved sulphate in water: from tributaries (Ftri), precipitation (Fp) and other sources (Fother) ?34SRS (‰) – is the measured isotopic composition of sulphur at the Bregana sampling location and contributions affecting isotopic composition of sulphate at this point: from tributaries (?34Stri), precipitation (?34Sp) and other sources (?34Sother) 0 0 0 174 Tja{a Kandu~ & Nives Ogrinc 20 18 16 -14 12 10 % 8 6 i 4 2 0 -2 dissolution of marine sulphates \ SULPHATE REDUCTION ¦ River Sava spring 2004 D River Sava tributaries spring 2004 ATMOSPHERIC SULPHATE SOIL WATER OXIDATION OF PYRITE oxidation of igneouus sulphate 10 15 20 25 30 35 40 45 SO42- (mg/l) gFuigruere 7 7. Scatter diagram of sulphur isotopic composition versus sulphate concentration of sulphate for the River Sava system (after Yang et al., 1996). Qi (m3/s) - is discharge in River Sava at the Bregana sampling location (RS) and contributions affecting discharge at this location from tributaries (tri), average annual precipitation amount in the watershed (p) and discharge of other sources (other) [SO42-]i - are the measured sulphate concentrations (mmol/l) at the Bregana sampling location (RS) and contributions affecting sulphate concentration at this location from tributaries (tri), precipitation (p) and other sources (other) Fothers (mol/s), 534Sothers (‰) are in equations (4) and (5) unknowns. Assuming that the sources of SO42- to the River Sava arise from the tributaries (Ftri), precipitation (Fp) and other sources (Fother), the contributions of these inputs were calculated to be 52 : 8 : 40 % for spring 2004 at the Bregana location (41). The three main tributaries contributing to the sulphate flux are the Trži{ka Bistrica (10.1 %), the Savinja (34 %) and the Krka (29.3 %). In precipitation a 534SSO4 value of 5.2 ‰ was measured in spring 2004. The annual sulphur flow in precipitation according to the average amount of precipitation of 1500 mm/a and concentration of sulphate in precipitation 1 mg/l in spring (Kandu~, 2006) is estimated to be 1.7 × 108 mol S/a. The flux from other sources, that accounts for 40 % of SO42-, may be derived from sources such as dissolution of evaporate minerals, the oxidation of sulphides, or from anthropogenic sources, such as air pollution, smelting of sulphide ores and the chemical industry. Since evaporates and sulphide minerals have a limited distribution in the bed rocks of the Sava channel, and the calculated ?34Sother value is 5.2 ‰, this suggests that evaporates are not the missing contributor. Industrial pollution probably accounts for the bulk of the “other” flux, amounting 8.2 × 108 mol/a of sulphur. Unfortunately, the isotopic composition of sulphur from these industrial operations is not available and this precludes a more precise assessment of industrial impact on the River Sava system. It is known that the Cinkarna Celje factory/industrial complex (sampling location 33) produced 100000 t/a of sulphuric acid from the year 1975 and represents a huge ecological problem for the Celje region (Podpe~an, 1988). From discharge and concentration measurements (EIONET, 2005), the annual sulphate flux was at Bregana (41) was calculated to be 1.5 × 109 mol/a. According to Ivanov (1983), the annual SO42- flux from continents to oceans is 6.8 × 1012 mol/a, with contributions from rivers, anthropogenic emissions, and groundwater of 48, 48, and 4 %, respectively. 0 5 Hydrogeochemical characteristics of the River Sava watershed in Slovenia 175 Conclusion Acknowledgements In the carbonate part of the Sava watershed lower Ca2+ and Mg2+ concentrations and lower alkalinity were found in comparison to the central and lower part of the river, which is composed of carbonate and clastic rocks. The partial pressure of CO2 ranges from 129 ppm to 9120 ppm and on average is 5 times oversaturated in comparison to atmospheric CO2 (316 ppm). The River Sava and its tributaries represent a source of CO2 to the atmosphere in all sampling seasons, especially in late summer months. Ther-modynamic calculations indicate precipitation of calcite and dolomite in the river system, but ?13CPIC values reveal detritical carbonate. ?13CDIC values change seasonally from – 13.5 to – 3.3 ‰ and reflect relations between the terrestrial, riverine and atmospheric reservoirs and the following processes: (1) influence of tributaries, (2) dissolution of carbonates, (3) degradation of organic matter and (4) equilibration with atmospheric CO2. The annual long term flux of dissolved inorganic carbon is estimated to be 4.1 × 1011 g C/year, and according to drainage area, 4 × 107 g C/year km2 at the border with Croatia. ?18OH2O values in the River Sava change seasonally from – 11.0 to – 8.3 ‰, while in the tributaries they change from – 12.4 to – 8.6 ‰. More positive values are attributed to the evaporation process. More negative ?18OH2O values are detected in spring months due to snow melting, while in the winter months ?18OH2O values are comparable to the late summer months. In the Sava watershed sulphate concentrations change seasonally from 4.34 to 55.5 mg/l. Concentrations of sulphate and ?34SSO4 values indicate that the source of the Sava Dolinka spring is mainly recharged from precipitation. From discharge data, concentrations of sulphate in water and the drainage area, the sulphate flux is estimated to be 1.4 × 107 g SO4/km2 at the border with Croatia. Assuming that the sources of SO42-to the River Sava are its tributaries, precipitation and other sources, the ratio of the contributions of these inputs are calculated to be 52 : 8 : 40 %. Other sources could be attributed to industrial pollution and oxidation of sulphide minerals. The calculated ?34SSO4 value for other sources is 5.1 ‰, which indicates that the greatest contribution of these other source is from industrial activity. The authors are thankful to Mr. Stojan Žigon for technical support. This research was conducted in the framework of the project L2-6458-792 funded by the Slovenian Research Agency (ARRS) and within the EU 6th Framework Specific Targeted Research Project – SARIB (Sava River Basin: Sustainable Use, Management and Protection of Resources), Contract No. INCO-CT-2004-509160. The project was also financially supported by the National Science Foundation, USA (NSF-EAR#0208182). Sincere thanks to Anthony Byrne for improving the English of the manuscript. References Agency of Republic of Slovenia for the Environment (ARSO), 2004: monthly bulletin. – No. 1–12. (in Slovene). Appelo, C. A. J. & Postma, D. 1994: Geochemistry, groundwater and pollution. – A. A. balkema, 536 pp., Rotterdam/Brookfield. Atkins, P. W. 1994: Physical chemistry. – Oxford University press, 1031 pp., Oxford. Barth, J. A. C., Cronin, A. A., Dunlop, J., & Kalin, R. M. 2003: Influence of carbonates on the riverine carbon cycle in an antropogenically dominated catchment basin: evidence from major elements and stable carbon isotopes in the Lagan River (N. Ireland). – Chemical Geology, 200, 203–216. Buser, S. & Draksler, V. 1989: Geological composition of Slovenia. In: Javornik, M., Voglar, D., Dermastia, A. (Eds.). Encyclopedia of Slovenia, 3, Mladinska knjiga Ljubljana, 200–201 (in Slovene). Clark, I. & Fritz, P. 1997: Environmental Isotopes in Hydrogeology. – Lewis Publishers, 328 pp., New York. Cruz – San, J., Araguas, L., Rozanski, K., Benavente, J., Cardenal, J., Hidalgo, M. C., Garcia – Lopez, S., Martinez – Garrido, J. C., Moral, F. & Olias, M. 1992: Sources of precipitation over South–Eastern Spain and ground-water recharge. An isotopic study. – Tellus, 44 B, 226–236. Dansgaard, W. 1964: Stable isotopes in precipitation. – Tellus, 16, 436–468. Dever, L., Durand, R., Fontes, J. Ch. & Vai-cher, P. 1983: Etude pédogénétique et isotopique des néoformations de calcite dans un sol sur craie. – Caractéristiques et origines. – Geochim. Cosmo-chim. Acta, 47, 2079–2090. Elderfield, H., Upstill-Goddard, R. & Sholkovitz, E. R. 1990: The rare earth elements in rivers, estuaries and coastal seas and their significance to the composition of ocean waters. – Geochimica Cosmochimica Acta, 54, 971–997. EIONET (2005) European Environment Information and Observation Network http://eio-net-eu.arso.gov.si Cited 11 Nov 2005. 176 Tja{a Kandu~ & Nives Ogrinc Fairchild, I. J., Bradby, L., Sharp, M. & Ti-son, J.–L. 1994: Hydrogeochemistry of carbonate terrains in alpine glacial setting. Earth Surfaces Processes Landforms, 19, 33–54. Fairchild, I. J., Killawee, J. A., Hubbard, B. & Dreybrodt, W. 1999. Ineractions of calcareous suspended sediment with glacial meltwater, a field test of dissolution behaviour. – Chemical geology, 155, 243–263. Fairchild, I. J., Borsato, A., Tooth, A. F. , Frisia, S., Hawkesworth, C. J. Huang, Y., McDermott, F. & Spiro, B. 2000: Controls on trace element (Sr–Mg) compositions of carbonate cave waters: implications for spelothem climatic records. – Chemical geology, 166, 255–269. Gaillardet, J., Dupré, B. & Allegre, C. J. 1999: Global silicate weathering and CO2 consumption rates deduced from the chemistry of large rivers. – Chemical geology, 159, 3–30. Gao, W. & Kempe, S. 1987: The Changjiang: its long – term change in pCO2 and pCO2 carbonate mineral saturation. In: Degens E. T. (Ed.), Transport of Carbon and Minerals in Major World Rivers, Mitt. Geol. – Paläont. Inst. Univ. Hamburg, 64, 207–216. Gibbs, R. J. 1972. Water chemistry of the Amazon River. – Geochimica Cosmochimica Acta, 36, 1061–1066. Goldstein, S. J. & Jacobsen, S. B. 1988: Rare earth element in river waters. – Earth Planetary Science Letters, 89, 35–47. Gonfiantini, R. 1986: Environmental isotopes in lake studies. In: Fritz, P. , Fontes J. Ch. (Eds.) Handbook of Environmental Isotope Geochemistry., 2, 113–168. Hrvatin, M. 998: Discharge regimes in Slovenia. – Geografski zbornik, XXXVIII, 60–87. Hu, M.–H., Stallard, R. F. & Edmond, J. M.. 1982: Major ion chemistry of some large Chinese rivers. – Nature, 298, 550 – 553. Huh, Y., Tsoi, M. Y., Zaitsev, A. & Edmond, J. M. 1998: The fluvial geochemistry of the rivers of Eastern Siberia: I. Tributaries of the Lena River draining the sedimentary platform of the Siberian Craton. – Geochimica Cosmochimica Acta, 62, 1657–1676. IAEA 2005: Isotopic composition of precipitation in the Mediterranean Basin in relation to air circulation patterns and climate. Final report of a coordinated research project 2000–2004. IAEA-TECDOC-1453, Vienna : IAEA, 223 pp. Ivanov, M. V. 1983: The sulfur cycle in continental reservoirs. In: Ivanov M. V. , Freney R. J. (Eds.) – The Global Biogeochemical Sulphur Cycle. 19, Wiley, SCOPE 19, 297–356. Javornik, M., Voglar, D. & Dermastia A., (Eds.) 1988: Encyclopedia of Slovenia, 2, pp. 68, Mladinska knjiga, Ljubljana. Kandu~, T. 2006: Hydrogeochemical characteristics and carbon cycling in the River Sava watershed in Slovenia. – Ph. D. Thesis, pp. 141, University of Ljubljana (in Slovene). Kandu~, T., Szramek, K., Ogrinc, N. & Walter M. L. 2007: Origin and cycling of riverine inorganic carbon in the Sava River watershed (Slovenia) inferred from major solutes and stable isotopes, Biogeochemistry in press. Karim, A. & Veizer, J. 2000: Weathering processes in the Indus River Basin: implications from riverine carbon, sulphur, oxygen, and strontium isotopes. – Chemical geology, 170, 153–177. Kempe, S. 1982: Long-term record of CO2 pressure fluctuations in fresh waters: In Degens E.T. (Ed.) Transport of Carbon and Minerals in Major World Rivers, Part 1, 52, pp. 91–332. Mitt. Geol.-Paläont. Inst. Univ. Hamburg. Krouse, H. R. 1980: Sulphur isotopes in our environment. In: Fritz P. , Fontes J. Ch. (Eds.) Handbook of Environmental Isotope Geochemistry I The Terrestrial Environment, Elsevier, Amsterdam, 435–472. Liu, Z. & Zhao, J. 2000: Contribution of carbonate rock weathering to the atmospheric CO2 sink. – Environmental geology, 39, 1053– 1058. Meybeck, M. 1996: River water quality: global ranges, time and space variabilities, proposals for some redefinitions. – Verh. Internat. Verein Limnol., 26, 81–96. O’Neil, J. R. 1979: Stable Isotope Geochemistry of Rocks and Minerals. – V: Lectures in Isotope Geology, Jager, E., Hunzinger, J. C., (Eds). – Springer Verlag, 235–263, Berlin. Palmer, S. M., Hope, D., Billett, M. F. , Dawson, J. J. & Bryant, C.L. 2001: Sources of organic and inorganic carbon in a headwater stream: evidence from carbon isotope studies. – Biogeochemistry, 52, 321–338. Parkhurst, D. L. & Appelo, C. A. J. 1999: User’s guide to PHREEQC for Windows (version 2) – a computer program for speciation, batch – reaction, one – dimensional transport, and inverse geochemical calculations. – Water – Resources Investigations Report, 99–4259. Pawellek, F. , Frauenstein, F. , & Veizer, J. 2002: Hydrochemistry and isotope geochemistry of the upper Danube River. – Geochimica et Co-smochimica Acta, 66/21, 3839–2854. Pearson, F. J. & Rightmire, C. T. 1980: Sulphur and oxygen isotopes in aqueous sulphur compounds. In: Fritz P. and Fontes J. Ch. (Eds.) Handbook of Environmental Isotope Geochemistry, 1, 227–259, Elsevier. Petkov{ek, Z. 1969: Fog occurrence in lowlands of Slovenia, Thesis. – Dru{tvo meteorologov Slovenije, 11, 57–89 (in Slovene). Podpe~an, D. 1988: Cinkarna Celje. In: Javornik M., Voglar D., Dermastia A., (Eds.) Encyclopedia of Slovenia, 2, Mladinska knjiga, Ljubljana, pp. 68 (in Slovene). Ramov{, A. 1983: Waterfalls in Slovenia. – Slovenska matica, 1983, 188 pp., Ljubljana. Reeder, S. W., Hitchon, B. & Levinson, A. A. 1972: Hydrogeochemistry of the surface waters of the Mackenzie River drainage basin, Canada: 1. Factors controlling inorganic composition. – Geo-chimica Cosmochimica Acta, 36, 181–192. Stallard, R. F. & Edmond, J. M. 1983: Geochemistry of the Amazon: 2. The influence of geology and weathering environment on the dissolved load. – Journal of Geophysical Research, 88, 9671–9688. Telmer, K. & Veizer, J., 1999: Carbon fluxes, pCO2 and substrate weathering in a large northern river basin, Canada: carbon isotope perspectives. – Chemical geology, 159, 61–86. Urbanc, J. & Brancelj, A. 2002: Hydrologi-cal Connections between some lakes in the Triglav Lakes Valley. In: Brancelj A. (Ed.) High – mountain lakes in the eastern part of the Julian Alps., Ljubljana : ZRC SAZU, National Institute for biology, pp. 77–90. Hydrogeochemical characteristics of the River Sava watershed in Slovenia 177 Vre~a, P. , Kandu~, T., Žigon, S. & Trkov, Z. 2004: Isotopic composition of precipitation in Slovenia. In: Isotopic composition of precipitation in the Mediterranean basin in relation to air circulation patterns and climate: final report of a coordinated research project 2000–2004, (IAEA-TECDOC, 1453). – Vienna: IAEA, 2004, pp. 157–172. Yang, C., Telmer, K., & Veizer, J. 1996: Chemical dynamics of the ’St. Lawrence’ riverine system: ?DH2O, ?18OH2O, ?13CDIC, ?34Ssulfate and dissolved 87Sr/86Sr. – Geochimica Cosmochimica Acta, 60, 851–866. Yee, P. , Edgett, R., & Eberhardt, A. 1990: Great lakes – St. Lawrence River regulation; what it means and how it works. Joint publication of Environment Canada Ontario Region, and the U. S. Army Corps of Engineers. Zhang, J., Quay, P. D. & Wilbur, D. O. 1995: Carbon isotope fractionation during gas – water exchange and dissolution of CO2. – Geo-chimica Cosmochimica Acta, 59/1, 107–1146. Zupan~i~, B. 1998: Precipitations. – Geografski atlas Slovenije, DZS, 99 pp., Ljubljana.