GEOLOGIJA 59/2, 179-192, Ljubljana 2016 http://dx.doi.org/10.5474/geologija.2016.010
© Author(s) 2016. CC Atribution 4.0 License
Hydrogeochemical and isotopic characterization of Pesnica
River, Slovenia
Hidrogeokemične in izotopske značilnosti reke Pesnice
Tjaša KANDUČ1*, Zoran SAMARDŽIJA1b, Nataša MORI3, Andreja JEREBIC3, Ines LEVAČIČ4, Miha KRAČUN4, Johanna A. ROBINSON1- 2, Stojan ŽIGON1, Željko BLAŽEKA4 & David KOCMAN1
1Department of Environmental Sciences, 1bDepartment for Nanostructured Materials, Jožef Stefan Institute, Jamova cesta 39, SI-1000 Ljubljana, Slovenia; e-mail: tjasa.kanduc@ijs.si*, zoran.samardzija@ijs.si, johanna.a.robinson@ijs.si, stojan.zigon@ijs.si, david.kocman@ijs.si Jožef Stefan International Postgraduate School, Jamova cesta 39, SI-1000 Ljubljana, Slovenia
3Department of Organisms and Ecosystem Research, National Institute of Biology, Večna pot 111, SI-1000 Ljubljana, Slovenia; e-mail: natasa.mori@nib.si, andreja.jerebic@nib.si 4Institut za ekološki inženiring d.o.o., Ljubljanska ulica 9, SI-2000 Maribor, Slovenia; e-mail: ines.levacic@iei.si,
miha.kracun@iei.si, zeljko.blazeka@iei.si
Prejeto / Received 11. 10. 2016; Sprejeto / Accepted 9. 11. 2016; Objavljeno na spletu / Published online 23. 12. 2016
Key words: River Pesnica, hydrogeochemical composition, stable isotopes, carbon, biogeochemical processes
Ključne besede: reka Pesnica, hidrogeokemijska sestava, stabilni izotopi, ogljik, biogeokemijski procesi
Abstract
The surface water geochemistry and carbon cycling studied here depend mainly on geological composition of the river catchment. The investigated surface waters in the river Pesnica catchment that are heavily hydromorphologicaly altered (reservoirs, channelization, land melioration systems) represent waters influenced by chemical weathering of carbonates and of mostly clastic rocks (claystone, sandstones, marlstones and siltstones). The objectives of our study were to analyze hydrogeochemical and isotopic composition of dissolved (¿13Cdic) and particulate carbon (á13Cpoc) with characterization of suspended matter and evaluate biogeochemical processes in Pesnica River in Slovenia before and after the retention. Surface waters are generally close to saturation regarding calcite and dolomite, dissolved CO2 is 49 to 1000 times oversaturated relative to atmosphere. ¿13CdiC was in the range from -14.8 to -4.2 %o and shows following biogeochemical processes in river system: degradation of organic matter, dissolution of carbonates and biological activity, which was confirmed with SEM microscopy and EDXS microanalysis. Results of SEM/EDXS showed that suspended matter is composed of K, Mg and Ca from locations Pesnica 1 and Pesnica 2 (above retention Perniško lake). Sample from Pesnica 3 (below retention) shows higher biological activity, while sample from Pesnica 4 (tributary of Pesnica, drainage channel Biš) contains a lot of particles of microporous structure and fiber structure, which are probably of anthropogenic origin. ¿13CpoC is changing from -29.5 to-27.6%o and showed different stages of degraded terrestrial material. More negative ¿13CpoC are shown at location at drainage channel Biš and show higher terrestrial input in river system as locations 1, 2 and 3 Pesnica, which have higher c>13CpoC values. This first results about suspended matter in Pesnica watershed serve for evaluation of anthropogenic influences specially in relation with further influences and longer time planning of sustainable river water management.
Izvleček
Geokemija površinskih voda in kroženje ogljika v naši študiji v glavnem zavisita od geološke podlage v porečju. Raziskane površinske vode v porečju Pesnice, ki je hidromorfološko močno preoblikovano (vodni zadrževalniki, urejanje struge, melioracije), predstavljajo vode, na katere vpliva kemijsko preperevanje karbonatov in v glavnem klastičnih kamnin (glinavcev, peščenjakov, laporovcev in meljevcev). Namen naše raziskave je bil analizirati hidrogeokemično in izotopsko sestavo raztopljenega in partikulatnega ogljika skupaj s karakterizacijo suspendirane snovi in oceno biogeokemijskih procesov v porečju Pesnice v Sloveniji pred zadrževalnikom in za njim. Površinske vode so večinoma blizu nasičenja glede na kalcit in dolomit, raztopljen CO2 je od 49 do 1000 krat prenasičen relativno na atmosfero. ¿13CdiC je bil v razponu od -14.8 do -4.2 %o in kaže na sledeče biogeokemijske procese v rečnem sistemu: razgradnja organske snovi, raztapljanje karbonatov in biološka aktivnost, kar je bilo potrjeno s SEM mikroskopijo in EDXS mikroanalizo. Rezultati SEM/EDXS so pokazali, da je suspendirana snov
180 T. KANDUČ, Z. SAMARDŽIJA, N. MORI, A. JEREBIC, I. LEVAČIC, M. KRAČUN, J. A. ROBINSON, S. ŽIGON, Ž. BLAŽEKA & D. KOCMAN
iz Pesnice 1 in Pesnice 2 (nad zadrževalnikom Perniško jezero) sestavljena iz alumosilikatov s K, Mg in Ca. Vzorec iz Pesnice 3 (pod zadrževalnikom) kaže na veliko biološko aktivnost (diatomeje), medtem ko vzorec iz Pesnice 4 (pritok Pesnice, drenažni kanal Biš) vsebuje veliko delcev mikroporozne strukture in vlaknatih struktur, ki so po vsej verjetnosti antropogenega izvora. ¿13CPOC se spreminja od -29.5 do -27.6 %o in kaže na različne stopnje razgrajenega terestričnega materiala. Bolj negativne vrednosti ¿13CPOC opazimo na lokaciji v drenažnem kanalu Biš in kažejo na večji terestrični vnos v rečni sistem kot lokacije Pesnice 1, 2 in 3, ki imajo višje vrednosti ¿13CPOC. Ti prvi rezultati o suspendirani snovi v porečju Pesnice služijo za evalvacijo antropogenih vplivov, posebno v povezavi s prihodnjimi vplivi in dolgoročnim načtrovanjem trajnostnega upravljanja s porečjem.
Introduction
It is well known that rivers reflect the biogeo-chemical processes occurring in their catchment areas (Meybeck, 1981; Meybeck, 1987; Popescu et al., 2015; Li et al., 2016). Understanding the organic and inorganic carbon cycles in river systems is important because it helps to understand the biogeochemical processes in a river and its catchment basin (Karim & Veizer, 2000; Barth et al., 2003; Sameer et al., 2016). Several studies evaluating biogeochemical processes in Slovenian rivers have been carried out in the past ten years, contributing substantially to a general understanding of their interactions with their catchments (Kanduc & Ogrinc, 2007; Kanduc et al., 2008; Zavadlav et al., 2013; Kanduc et al., 2013).
Suspended organic matter is known to be complex mixture of molecules such as carbohydrates, amino acids, fatty acids and phenols, particles from wastewaters and industrial water, soil organic matter and biological material like phytoplankton and plant parts (Hope et al., 1994). Suspended matter in rivers is derived mostly from soil and plant material, so suspended organic matter (c>13CPOC) in rivers can be used to determine the contribution of terrestrial vegetation and soil matter in the river ecosystem (Hedges, 1992).
Pesnica is a river that was extensively regulated in the 60's when it was channeled and when multipurpose retention basins were built to mitigate floods (Leitinger, 2012). Nowadays, droughts in combination with severe loads of nutrients and other pollutants from households and agriculture adversely affect the flora and fauna in the area (Markoli, 2012).
The main objective of our study was to understand the hydrogeochemical characteristics and prevailing biogeochemical processes in the Pesnica catchment waters using geochemical and stable carbon isotopes analysis, respectively.
Geological and geomorphological setting
River Pesnica is located in northeastern Slovenia (Fig. 1). It is a left tributary of River Drava, and is the main watercourse in Slovenske gorice. The river originates above Pesnica (in Austria) at an altitude of 300 m; at town Ormož (altitude 190 m) it flows out to the Drava River (Fig. 1). Its length is 69 km, of which 65 km being in Slovenia (Statistical Office of Republic of Slovenia, 2002). Its catchment area comprises 539 km2. Its mean average discharge between 1926 and 1965 was 6.72 m3/s. It has a mixed snow/rain river regime with high water levels in March and November and low in August and September. In Pesnica catchment area, 3 gauging stations are monitored by the Slovenian Environmental Agency (SEA): Ranca, Gočova, Zamušani I Kobold, 2012).
Land cover in Pesnica catchment is dominated by agriculture areas (~75 %), followed by forests and semi-natural areas (~23 %) and with less than 2% of artificial surfaces (CLC, 2012).
Since Pesnica flows through a lower bushy environment, it has a low specific runoff (10.6 l/s/ km2; gauging station Zamušani) and a runoff coefficient of 32.2 % (gauging station Zamušani).
In past years it was found, that mainly in summer, the discharge of Pesnica is mostly below 0.5 m3/s (Urbanek, 1994). Since it flows through an intensively agricultural area, it contains large amount of atrazine and steroid organic substances, together with a high concentration of nitrates in past decade. Pollution with sewage sludge was moderate, and in recent years the river's ecological status has improved (Urbanek, 1994).
Pesnica River is heavily hydromorphologically altered (from Zgornja Kungota to outflow). For the purposes of flood prevention, several reservoirs have been constructed to retain the water: Pernica and tributaries were artificial lakes were constructed for water retention.
Hydrogeochemical and isotopic characterization of Pesnica River, Slovenia
181
Drainage channel at settlement Bis-the municipality Trnovska vas, is a right hand tributary of Pesnica River (P4 in Fig. 1).
Geologically, Pesnica catchment area belongs to the Panonian basin and is of Tertiary and Quaternary ages. The catchment is composed of Miocene and Pliocene sediments (siltstones, marlstones, claystones, dolomite/calcite) and quaternary sediments. The area of the Panonian basin was split into more tectonic units with the Slovenske gorice as the largest tectonic unit. In the Pesnica catchment more faults were found in different directions, belonging to the category of tectonic faults. The most important are: Globov-nica, Pesnica, Dravinja faults (named after the tributaries that flow through them). Surface water creates channels, in which the rocks are the most fractured (Mioc & Znidarcic, 1989; Perko et al., 2001). Following Miocene and Quaternary sediments are present in the River Pesnica watershed:
Miocene sediments (M21) that outcrop up in the northeastern part of the catchment are composed of conglomerate, sandstone and marl.
M22 - Miocene sediments are composed of gravel, sand, marl, limestone and sandstone.
M,,1 - In the lower part of the catchment sandy marl sediments prevail and within them are mica, sandstone, claystones and marlstones. Gravels are composed of siliciclastic and other metamorphic and magmatic rocks.
M„
Miocene sediments are located in the
eastern part of Pesnica catchment and are composed of sandy clays, sandy clay marls, sands and sandy gravels.
Pliocene sediments are located in the southern part of the River Pesnica catchment (Pliocene) and are composed of clay, sand, conglomerate and clayey marl (Mioc & Znidarcic, 1989).
Quaternary sediments are located along surface water flow in the valley and are composed of clastic material (gravels, sands, marls, clays).
2
Fig. 1. Land use map of River Pesnica catchment area (Corine Land Cover 2012)
182 T. KANDUČ, Z. SAMARDŽIJA, N. MORI, A. JEREBIC, I. LEVAČIC, M. KRAČUN, J. A. ROBINSON, S. ŽIGON, Ž. BLAŽEKA & D. KOCMAN
Materials and Methods
In order to determine how artificial reservoirs influence biogeochemical processes in the river system, water was sampled at the inflow and outflow of one of the accumulation lakes (Pernisko jezero) during the crop growing season. In addition, a small drainage channel nearby Bis village was chosen to investigate how an area dominated by intensive agriculture contributes to the biogeochemical processes of Pesnica River.
Discharge data (Miha Supek, personal communication) were obtained from the Environmental Agency of the Republic of Slovenia (Tab. 1). Samples were collected four times in each month from February to September (except July) 2016 at four gauging stations: 1 (before dam), 2 (before dam), 3 (after dam), 4 (drainage channel Bis). Temperature and dissolved oxygen (DO) of surface waters were measured with a WTW Multi 340, CellOx 325. pH was measured with a Seven Compact pH/ ion Metler Toledo and electroconductivity with an Ultrameter II MYRON L COMPANY 6 PFCE, serial 6218180, both in the laboratory after sampling.
Samples for characterization of suspended matter from each location (Fig. 1) were collected (September 2016 sampling) using standard representative sampling techniques (Schuster & Reddy, 2001) in HDPE bottles. From each location, 3 l of surface water were sampled for microscopic (SEM/EDXS-scanning electron microscopy/energy - dispersive X-ray spectroscopy) analyses (suspended matter characterization) and for isotopic composition of carbon in suspended matter (¿13CpOG). Surface water samples were filtered through a Whatman GF/F (pore size 0.7 ^m) glass fiber filter. The isotopic composition of carbon in suspended matter was determined using a Europa 20-20 continuous flow IRMS-SL preparation module. The filter was treated with 1M HCl to dissolve any carbonate present. Approximately 1 mg of sample was scrapped from the GF/F filter and dried at 60 °C. The isotopic composition of carbon was determined after combustion of the capsules in a hot furnace at 1000 °C. IAEA CH3 and IAEA CH6 standards were used to relate the analytical results to the VPDB. Sample reproducibility for carbon was ±0.2
Sample aliquots collected for cation, anion and alkalinity analyses were passed through a 0.20 ^m nylon filter into HDPE bottles (30 ml) and kept refrigerated until analyzed. Samples for
¿13C
analysis were passed through a 0.45 ^m
nylon filter and stored in glass serum bottles filled with no headspace and sealed with septa caps.
Alkalinity was measured using Gran titrations (Clesceri et al., 1998). Cations (Ca2+, Mg2+, Na+, K+) and anions (SO42-, Cl-, and NO3-) were measured by ion chromatography (Metrohm, 761 Compact IC).
The stable isotope composition of dissolved inorganic carbon (¿13CDIC) was determined with a Europa Scientific 20-20 continuous flow IRMS (isotope ratio mass spectrometer) with an AN-CA-TG preparation module. Phosphoric acid (100 %) was added (100-200 ^l) to a septum-sealed vial which was then purged with pure He. The water sample (5 ml) was injected into the septum tube and headspace CO2 was measured (modified after Kanduc, 2006). In order to determine the optimal extraction procedure for surface water samples, a standard solution of Na2CO3 (Carlo Erba) with a known ¿13CDIC of -10.8±0.2 %o was prepared with a concentration of 0.66 g/l.
All stable isotope results for carbon are expressed in the conventional delta (¿) notation, defined as per mil (%o) deviation from the reference standard VPDB. Precision of working standards was ±0.2 %o for c>13Cdic and c>13Cpoc.
Samples of suspended matter from Pesnica catchment area from September sampling were examined to study the morphology of particles and their element composition, using a field-emission-gun scanning electron microscope (FEGSEM) JEOL JSM 7600-F equipped with Oxford Instruments INCA energy -dispersive X-ray spectroscopy (EDXS) with 20-mm2 silicon-drift detector (SDD). SEM/EDXS analyses were performed at 20kV accelerating voltage, 15 mm working distance and 100 s acquisition time. The internal INCA EDS standard database was used for the analyses.
Thermodynamic modeling was used to evaluate pCO2 (partial pressure) and the saturation state of calcite (SI , ) and dolomite (SI, , )
v calcite	v dolomite'
using pH, alkalinity and temperature as inputs to the PHREEQC speciation program (Parkhurst & Appelo, 1999). Version 2.16.03 Phre-eqc for Windows was used to calculate SI , .4 ,
calcite
SIdo,om,e and pcO2.
Hydrogeochemical and isotopic characterization of Pesnica River, Slovenia
183
Results and discussion	„	, J,
for example the congruent dissolution of cal-
Theoretical background of dissolution of	cite (CaCO3) in limestone: CaCO3+ CO2(g) + H2O
minerals	= Ca2++ 2HCO„-
Dissolved Ca2+ and Mg2+ are largely supplied by the weathering of carbonates, with smaller contributions from silicate weathering, as indicated by the relatively high HCO3- and low Si concentrations. Most of the water samples in our study deviate from a 2:1 molar ratio of HCO3-(Fig. 2A) and are the consequence of weathering of clastic rocks composed of clay minerals as follows (1-5):
CaAl2Si2O8 + 3H2O + 2CO2 (Anorthite) ^ Al2Si2O5(OH)4 + Ca2+ + 2HCO3- (Kaolinite)
(1)
In the case of weathering anorthite, only Ca2+ is present in carbonate reactions (precipitation of calcite).
NaAlSi3O8 + CO2 + 11/2H2O (Albite) ^ Na+ + y2 Al2Si2O5(OH)4 + 2H4SiO4 + HCO3-(Kaolinite)	(2)
In case of weathering of albite, Na+ and HCO3-are in 1:1 stoichiometry ratio. The origin of K+ can be mica shist.
Weathering of calcite and dolomite is presented below:
Calcite :
CaCO, + CO,., + HO » Ca2+ + 2HCO,-	(3)
3 2(g) 2 3
Dolomite :
Ca05Mg05(CO3) + CO2 + H2O » Ca2+ + Mg2+ + 2HCO3-
(4)
The carbon dioxide content of rain and soil water is particularly important in weathering processes (Gaillardet et al., 1999). The chemical composition of rainwater changes markedly after entering the soils (Gaillardet et al., 1999). The upper part of the soil is a zone of intense biochemical activity. One of the major biochemical processes of the bacteria is the oxidation of organic material, which leads to an increase of carbon dioxide in the soil gas (Gaillardet et al., 1999). Above the zone of water saturation, the soil gases may contain 10 to 40 times as much as carbon dioxide as the free atmosphere (reference value for partial pressure in atmosphere is 316 ppm (Clark & Fritz, 1997)). This CO2 gives rise to a variety of weathering reactions,
and the incongruent reaction with K-spar:
(KAlSi3O8): 2KAlSi3O8+ 2CO2+ 11H2O = Al2Si2O5(OH)4+ 2K++ 2HCO3 + 4H4SiO2	(5)
The amount of CO2 dissolved according to reaction (3) depends mainly on the temperature and its partial pressure. For example, for an atmospheric carbon dioxide pressure of 10-2 atmosphere and for a soil atmosphere of nearly pure carbon dioxide, the amount of calcium that can be dissolved (at 25 °C) until saturation is 65 and 300 ppm, respectively. The calcium and hy-drogencarbonate ions released into soil water and groundwater eventually reach the river system. The water resulting from reaction (5) contains hydrogencarbonate, potassium, and dissolved silica in the ratios 1:1:2, and the new mineral, kaolinite, is the solid weathering product. The dissolved constituents of reactions (calcite) and (K-feldspar) eventually reach the river systems. The global silicate weathering fluxes and associated CO2 consumption fluxes have been freshly estimated on the basis of data on the 60 largest rivers of the world (Gaillardet et al., 1999). Only active physical denudation of continental rocks was found to be able to maintain high chemical weathering rates and significant CO2 consumption rates.
General geochemical characteristics of river water
Mean daily discharge data (Q), parameters measured in the field (T-temperature, D.O.-dissolved oxygen (mg/l)), geochemical and isotopic data of water (Pesnica catchment area) are presented in Table 1.
Most of the studied rivers in Slovenia (Kanduc et al., 2008) have a Ca2++Mg2+: HCO3-ratio 2:1 indicating that geological composition (carbonates prevail) controls the composition of river waters. In Pesnica River catchment (studied locations Pesnica 1, Pesnica 2, Pesnica 3 and Pesnica 4) clastic rocks prevail (Fig. 1) meaning that most of the cations (Ca2+, Mg2+, Na+, K+) are originated in clastic rocks (e.g. sandstones, claystones, marls) and also limestone, but to a minor extent.
184 T. KANDUČ, Z. SAMARDŽIJA, N. MORI, A. JEREBIC, I. LEVAČIC, M. KRAČUN, J. A. ROBINSON, S. ŽIGON, Ž. BLAŽEKA & D. KOCMAN
Table 1. Field and laboratory measurement data (T, pH, conductivity, DO), discharge data (Q), major geochemistry data (Ca2+, Mg2+, Na+, total alkalinity after Gran, SO42-, NO3-), stable isotopic composition of dissolved inorganic carbon (513Cdic), stable isotopic composition of particulate organic carbon (513CpOC).
Ca2++Mg2+ (mM)	r^. <N	m co <N	(N CO <N	& UD rH	r^. c^	r^. rH et	C^ rd	o un <N	r^. c^ rd	I.Û C^ rd	m m <N	co r-^	& r^ r^	Ln <N rd	O UD rH	Ln UD rH	Ln r-^ rd	O UD rd	et et rH	r-^ <N	C^ rH	(N rH rH	co co CD	(N et rH	et U^ C^	(N <N	et <N rH	(N m rH
<¿1	oo ct C) O	m r^. et o	r^. (N O m (N	m co (N co	Ln	m co Ln (N	& co et co Ln	(N et Ln r^. (N	I.D et Ln o	1.0 (N Ln	et (N co I.Û	r^. et	(N m (N O	o (N	m m et r^.	co co et Ln	I.D r^. m co	m co I.Û	co m	Ln co co co Ln	et r^. et O (N	& & et r^.	et m (N 1.0 Ln	1.0 (N (N O (N	et Ln (N et	I.Û (N (N O (N	co (N (N I.Û m	1.0 (N (N O (N
O	rH	co en rH	en rH	o <N	^ rH	o rH	O co rH	(.o <N	(N r-^ rH	O r-^ rH	r^. rH	m rH	<J) Ci rH	et C^ CM	o rH <N	C^ CM	co O <N	I.Û r^ CM	o CM	m <N CM	& UD rH	m c^ rH	Ln <N CM	(N C^ rH	Ln c^ rH	(N C^ rH	o u^ rH	(N C^ rH
"o -c uo	Ln CD	CO et CD	<N rH	m CD		r^. rH	I.Û O <N	et CD	m CD	C^ CD	et r^	r^. r-^ CD	I.Û rH	1.0 rH	et rH	UD CM	et <N CD	m un CD	m <N CD	rH	U^ CD	U^ CD	& r-^ CD	co r-^ CD	Ln en CD	UD CD	et UD CD	& 00 rH
JJ re u uo	<N O	<N CD	& <N CD	et en CD	m CD	et rH CD	Ln CD	co et O	Ln o CD	I.Û c^ CD	r^ CD	co rH O	r^ CD	(N C^ CD	C^ CD	I.Û CD	r^. CD	et U^ CD	et C^ CD	Ln m CD	CD	I.Û CD	co CD	& CD	& m CD	(N m CD	Ln CD	(N O
5 1uC POC (%o)	to d	to d	to d	d d	d	ré d	ré d	d d	!Z	ré d	d	d d	d d	ré d	ré d	!Z	ré d	ré d	d d	d d	d d	d d	TO !Z	TO !Z	UD r^		U^ co r^	
1 i«	rd tH	<N et	t-- rd	m rd	rH	oo et	m <N	et	rH	rH	m r-^	rd	t-- rd	rd	r^	C^ rd	m r^	r^		rd	oo rH	r^	co	rd	C^ rd	r^	en	rd
f££	oo 00 rd	r-^ rd	UD rd	r^. r-^ rH	Ln CO <t	Ln un	o <N un	d) rH	O un	O r^ un	o et	O c^ rd	O rd	O u^ rd	o r^ rd	et <N rH	Ln C^ un	C^ un	C^	et m et	O	o r-^ r^	o	O UD rH		et ud	r^. e^ r^	& et
Mg2+ (mM)	O Ln CD	o un CD	m CD	co m CD	& r-^ CD	o oo CD	r^. r-^ CD	Ln UD CD	m co CD	m co CD	& CD	r^. r-^ CD	m un CD	et UD CD	m CD	co m CD	m r-^ CD	O r-^ CD	O CD	co r-^ CD	Ln m CD	co r^ CD	co m CD	CD	m et CD	m et CD	co <N CD	1.0 CD
Ca2+ (mM)	r^. <n <n	m m <N	m <n	m rH	co rd	r^. m rd	Ln rd	rH	et C^ rd	m rd	et rH	c^ <N	r-^ rH		& r^ rH	r^. r^ rH	(N c^ rd	o r^	co m rH	c^	et u^ rH	et c^ CD	o un CD	C^ rH	r^	co C^ r^	Ci CD	Ln C^ rH
5 + 1	Ln O o	o CD	O CD	(N O CD	Ln C^ CD	Ln o CD	Ln o CD	C^ CD	1.0 C^ CD	1.0 C^ CD	I.Û C^ CD	C^ CD	Ln c^ CD	I.D C^ CD	m c^ CD	o CD	o CD	& o CD	co o CD	C^ CD	1.0 C^ CD	Ln O CD	r^. c^ CD	O O CD	Ln O CD	et O CD	O CD	O C^ CD
+ ? +1	(N et CD	et et CD	un CD	CD	^ CD	r^. un CD	m u^ CD	O <N CD	I.Û co CD	c^ CD	Ln UD CD	& rH CD	CD	(N u^ CD	Ln <N CD	(N CD	m c^ CD	& Ci CD	CD	r^. CD	co CD	et m CD	r^. r^ CD	co c^ CD	r^. CD	r^. CD	& CD	& O CD
o E UO -b-	Ln <n CD	Ln <n CD	& <n CD	^ CD	^ m CD	et m CD	m m CD	m r^ CD	m m CD	(N m CD	(N m CD	et C^ CD	co r^ CD	& CD	I.Û r^ CD	et r^ CD	r^. r^ CD	r^. r^ CD	r^ CD	m o CD	r^. r^ CD	r^ CD	r^. rH CD	(N O CD	CD	co m CD	r^. CD	& <D CD
irl	& tH CD	et CD	m CD	cfr CD	Ln c^ CD	Ln o CD	Ln c^ CD	<D CD	et O CD	Ln o CD	o CD	O C^ CD	O CD	r^. o CD	m o CD	Ln <N CD	r^. o CD	r^. o CD	O C^ CD	O C^ CD	r^. o CD	c^ CD	O C^ CD	O C^ CD	r^. o CD	O CD	Ln O CD	O O CD
ul	m m CD	(N cfr CD	en CD	r^. CD	Ln m CD	I.Û m CD	Ln m CD	co r^ CD	I.Û CD	r^. CD	o et CD	O <N CD	<J) <N CD	m m CD	& CD	& CD	r^. CD	co CD	<J) CD	(N CD	r^. CD	O CD	(N m CD	O CD	co et r^. et CD	r^. co (.o et CD	o I.Û m CD	& O Ln CD
m) u 3	to d	to d	TO d	TO !Z	ré d	ré d	ré d	d d	ré d	ré d	ré d	d d	!Z	ré d	ré d	TO !Z	ré d	ré d	d d	d d	ré d	ré d	ré d	d d	O O Ln	(N Ln	et m	co m
X O.	<N	(N <N	et	c^ r^	m r^	o	o m	O	(N	r^	m	m	m r^ r^	<N	m	Ln m r^	o un	(N un	co un	& un	o et	o m	o <N	m	r^	Ln et	(N <N	Ln <N
	TO C	to d	to d	TO !Z	ré d	ré d	ré d	d d	ré d	ré d	!=	d d	d d	ré d	ré d	TO !Z	o UD pC	o co	O (N rH (N	O c^ K	O Ln rH (N	o ro ro (N	o <JD (N	o en	o r^ et	O c^ et	o r^ CD I.Û	o m Ln
C i—	m <N un	m <N un	o <N un	O r^ un	m m <t	m m et	m m et	O c^ UD	o rd	o UD rd	O UD rd	O en un	o o et	o o et	o c^ et	O c^	o r^	O r^	O r^	o r^	m r^	u^ r^		r-^	o c^ co	O	o r^ (N	o m S
ah	O <N UD	O <N IJD	O <N IJD	(N <N	c^ CD	c^ CD	en CD	d) Ci	Ln r^ CD	Ln r^ CD	Ln r^ CD	r^. rH	rH	rH	& rH	I.Û	co CD	co CD	co CD	et CD	et m CD	et m CD	co Ln c^ CD	TO !Z	I.Û co c^ rH	I.Û co c^ rH	& r^. CD	d d
QJ TO Q	rH <N UD	rH <N	rH <N lO	<N	I.Û r^ co	I.Û rH od OO	(.o r^ co	rH od co	1.0	I.D rH et	I.D r^ et	I.Û r^	r^ U^ (JD	I.Û rH un	1.0 r^ U^ (JD	I.Û rH un	I.Û r^ rd	I.Û r^ rd	<■6 rd	I.Û rH 1.0 rd	I.Û r^ C^ L^	1.0 r^ C^ L^	(.o rH oo L^	C^ Ln	1.0 a^ o	1.0 d O r^	1.0 Gï O r^	(.o rH en <D
Location	TO U 'c fcO (U CL	(N TO U 'c w (U CL	m TO u C w (U CL	TITO U 'c w M	TO U iZ w (U CL	(N TO U 'c fcO (U CL	m TO U 'c w (U CL	et TO U 'c	TO iZ W (U CL	(N TO U 'c w (U CL	m TO *!= fcO (U CL	et TO U iZ	TO !Z W (U CL	(N TO U 'c w (U CL	m TO u 'c w (U CL	et TO U 'c cï	TO !Z W (U CL	(N TO !Z W (U CL	m TO u 'c w (U CL	et TO U 'c CÏ	TO U !Z W (U CL	(N TO U 'c w (U CL	m TO U 'c w (U CL	et TO U 'c cï ë	TO u !Z W (U CL	(N TO u 'c fcO (U CL	m TO u 'c w (U CL	et TO U 'c w M cï
o z		(N	m			(N	m	et		(N	m	et		(N	m	et		(N	m	et		(N	m	et		(N	m	
Hydrogeochemical and isotopic characterization of Pesnica River, Slovenia
185
6.00
• February 2016 ■ March 2016
Total alkalinity (mM)
Fig. 2A. Ca2++Mg2+ ratio versus alkalinity with lines 1:2 indicating weathering of carbonates; deviation from the 1:2 line indicates weathering of clay minerals (albite, anortite) and mica schist (origin of K+), which mostly contribute to alkalinity rather than to weathering of calcite/dolomite.
7.00
Dolomite only
0.75,,
s s" ✓ .y
/V
0.5 _____0:33"
___________________________
• February, 2016 ■ March, 2016 A April, 2016 XMay, 2016 O June, 2016 A August, 2016 — September, 2016
0.00	1.00	2.00	3.00	4.00	5.00	6.00	7.00
Ca2+ (mM)
1
Fig. 2B. Mg2+ and Ca2+ concentrations indicating the dominance of dolomite or calcite in the Pesnica catchment
The relation of Mg2+ versus Ca2+ indicates the relative contributions of calcite and carbonate to weathering intensity in the Pesnica catchment. Most of the samples indicate that weathering of calcite is dominant over the entire Pesnica River catchment (Fig. 2B). A Mg2+/Ca2+ ratio of around 0.33, typical for weathering of calcite with magnesium, is characteristic at all investigated locations in River Pesnica, only some samples having Mg2+/Ca2+ of around 0.5. A Mg2+/Ca2+ ratio of 0.33 is also characteristic of rivers comprising the Danube catchment, e.g., Sava River and its tributaries, Tisa River, Ilz River and Inn River (Kanduc et al. 2007; Szramek et al. 2007; Kanduc et al, 2013).
Thermodynamic modelling (pCO2, SI , SI )
From Fig. 3 it can be seen that the highest pCO2 was observed in September 2016 with a value of 316227.8 ppm (-1.50 bar) at location Pesnica 3 and the lowest value of 15488.2 ppm (-2.81 bar) at location Pesnica 4 (Fig. 3 A, Table 1). pCO2 concentrations in the river system depend mainly on conditions, such as temperature, pH and alkalinity. High concentrations of CO2 can arise from higher groundwater discharge and/or respiration in streams because of the input of organic matter. All the samples were supersaturated with CO2 with respect to the atmosphere (10-35 bar = 316 ppmv) according to thermodynamic conditions in the investigated period in 2016 (Fig. 3A).
186 T. KANDUČ, Z. SAMARDŽIJA, N. MORI, A. JEREBIC, I. LEVAČIC, M. KRAČUN, J. A. ROBINSON, S. ŽIGON, Ž. BLAŽEKA & D. KOCMAN
PH
0.00
-0.50
-1.00
-1.50
ft-2.00
M
O
-2.50
-3.00
A
Oversaturation with CO,
•february 2016 ■ March 2016 A April 2016 XMay 2016 O June 2016 A August 2016 -September 2016
Fig. 3A. pCO2 versus pH in river water system. All samples have pCO2 values greater than atmospheric (10-35 bar = 316 ppmv) causing increase in calcite saturation B. Saturation index of
calcite (SI ) versus dolov calcite7
mite (SI ) in river water
v dolomite7
system. Saturation for cal-cite/dolomite is defined as log IAP/K=0.
• 0.5
0
-1.2
-1 -0.8 -0.6
A
-0.4 -0.2
A -0.5 A
■ -1
• X
-1.5
-2.5
oversaturation with calcite
SW
0.2
►J?
rjl
0.4 0.6 O
0.8

-5°
a
oi
• February 2016
■ March 2016
A April 2016
xMay 2016
O June 2016
A August 2016
— September 2016
Fig. 3B. Saturation indices of calcite (SI ) versus
calcite
dolomite (SI, , .. ) in a river
dolomite
water system. Saturation for calcite/dolomite is defined as log IAP/K=0.
All samples in River Pesnica catchment were undersaturated with calcite and dolomite, only some samples in September 2016 season were oversaturated with calcite (Fig. 3B). SIcalcite ranged from -0.98 to 0.54 and SIdolomite ranged from -2.61 to 0.48. It has been clearly shown that most of the investigated rivers in the world constitute sources of CO2 to the atmosphere (May-orga et al., 2005). Oversaturation with SI , is
'	'	calcite
traced only in April at locations 1 and 2 before dam and oversaturation with dolomite only in June.
Isotopic composition of dissolved inorganic carbon (¿13CdiC) and particulate organic carbon
(¿13CPOC)
Values of the isotopic composition of dissolved inorganic carbon (¿13CdiC) can help decipher the contributions of organic matter decomposition, carbonate mineral dissolution, equilibration with atmospheric CO2 to DIC and in-stream biological processes (Kanduc et al., 2007). Figure 4 shows ¿13CdiC versus alkalinity for Pesnica river water from February 2016. ¿13C in Pesnica catchment ranges from -14.8
7.20
7.50
7.60
-3.50
-4.00 L
1
Hydrogeochemical and isotopic characterization of Pesnica River, Slovenia
187
0.00
1.00
2.00
Alkalinity (mM)
3.00	4.00	5.00
.00
7.00
.00
	4.0
	2.0
	0.0
	-2.0
	-4.0
	
IS	-6.0
y	
3	
U	-8.0
to	
	-10.0
	-12.0
	-14.0
	-16.0
	-18.0
4. Open system equilibration with the atmosphere with an average value of 1.6%o
3. Dissolution of carbonates with S13Ccarb of +2.4%o	o February 2016
■ March 2016 A April 2016
Pesnica 3	x May 2016
O
O June 2016
_ . _	A August 2016
Pesnica 2 Pesnica 3 Pesnica 3 A x	a	- September 2016
A Pesnica 3
2. Nonequilibrium carbonate dissolution by carbonic a^id produced from soil zone
with a S13Cco2 of -26.6%o	Pesnica 3
----------------------------------
Am	xx.A.
x ■ •
■
1. Open system equilibration of DIC originating from degradation of organic matter with S13CCO2 of -26.6%o
Fig. 4. Seasonal variation of i513Cdic values versus alkalinity in Pesnica river water system, with lines indicating biogeochemical processes in river water. These include (see equilibration lines): 4. Open system equilibration with the atmosphere with an average valueof 1.6 %o 3. Dissolution of carbonates according to the average <513CCaCO3 (+2.4 %o, Kanduc et al., 20133 of all measured samples composing the recharge area of river water, causing +1%o±0.2 %o enrichment in 12C in DIC (Romanek et al., 1992), 2. Nonequilib-rium carbonate dissolution by carbonic acid produced from soil zone CO2, 1. Open-system equilibration of DIC with soil CO2 originating from degradation of organic matter with <513Csoil = -26.6 %o.
to -4.2 %o (Table 1). An average ¿13Cpoc value of -26.6 %o (Kanduc et al., 2013) was assumed to represent the isotopic composition of DOC that was transferred to DIC (Dissolved Inorganic Carbon) by in-stream respiration. Open-system equilibration of DIC with CO2 enriches DIC in 13C by about 9 %o (Mook et al., 19 74), which corresponds to the value of -17.6 %o shown in Figure 4. Non-equilibrium dissolution of carbonates, with one part of DIC originating from soil CO2 (- 26.6 %o), and the other from carbonates with an average ¿13CCa of +2.4 %o, leads to an intermediate ¿13CDIC value of -12.1 %o (Kanduc et al., 2012). Given the isotopic composition of atmospheric CO2 (-7.8 %o, Levin et al., 1987) and the equilibrium fractionation with ¿13C of +9 %o, DIC in (¿13C ) was measured only in September and
2016) indicating that respiration can alter ¿13C-
to more positive ¿13CdiC values in-stream, especially in summer months (Fig. 4). Controls of longitudinal variation in ¿13CdiC in rivers have been studied by Roach et al., 2016. At location Pesnica 2 (before dam) and Pesnica 3 (after dam) ¿13CdiC values are less negative indicating CO2 reduction in river system and enrichment with 12C isotope (Fig. 4).
According to the model of Devol & Hedges (2001), suspended matter (particulate in continuation) may be divided into three fraction (Kanduc, 2010 & references therein). Isotopic composition of carbon in particulate organic matter
equilibrium with the atmosphere should have a ¿>13CdiC of about +1 %o (Fig. 4). Biogeochemical processes such as dissolution of carbonates and open-system equilibration with the atmosphere can be neglected in the Pesnica catchment since river water is composed mostly of clastic rocks (sandstones, marls, claystones).
Higher values of ¿13CdiC are observed at Pesnica 3 location (March, 2016) and slightly more negative values ¿13CdiC at Pesnica 2 (above the Pernica lake). From Figure 4 it can be seen that most of the samples from Pesnica River catchment fall around the line of nonequilib-rium carbonate dissolution by carbonic acid produced from the soil zone with ¿13CCO2 of -26.6 %o. Samples from Pesnica 3 (retention) have ¿13C values less negative up to -4.2 %o (June
ranged from -29.5 to -27.6 %o. At Biš location more negative ¿13Cpoc (-29.5 %o) was observed in comparison to characteristic degraded organic matter (-26.6 %o). For comparison, the c>13C value of particulate organic carbon (c>13CpoC) in the main channel of the Kamniška Bistrica River, a pre-alpine gravel-bed river with prevailing intensive agriculture in the catchment, varies from -28.6 to -25.4 %o and from -29.9 to -25.2 %o in the tributaries (Kanduc et al., 2013). Phyto-plankton binds CO2 from DIC in water leading to a 13C fractionation of around -20 %o, producing biomass with c>13CpoC ranging from -32.0 to -45.0 %o. At the River Pesnica catchment ¿13CpoC ranges from -27.6 to -29.5 %o, meaning that pho-tosynthetic activity was not so significant at our investigated locations.
188 T. KANDUČ, Z. SAMARDŽIJA, N. MORI, A. JEREBIC, I. LEVAČIC, M. KRAČUN, J. A. ROBINSON, S. ŽIGON, Ž. BLAŽEKA & D. KOCMAN
Composition of suspended matter differentiates by SEM/EDXS
The composition of suspended matter differentiates from Pesnica catchment area at investigated locations: Pesnica 1, Pesnica 2 (above Pernica lake), Pesnica 3 (below the Pernica lake) and Pesnica 4- drainage channel Biš (Fig. 1). At Pesnica 1 and Pesnica 2 sampling sites (Figs. 5A, B) we can observe that aluminosilicates (K, Al, Si) prevail in suspended matter, while at sampling location Pesnica 3 we observe high biological activity (diatoms) (Fig. 5C-b) and at Pesnica 4 (drainage channel Biš) inorganic fibers prevail (Fig. 5D), which could be of anthropogenic origin. All SEM Figures have corresponding spectra, meaning that elements: Al, Si, O, K, Na, Mg prevail with traces of Ti and Fe representing detrital geological material composed of alumi-
nosilicates (e.g. anorthite, albite, mica schist). From Fig. 5B with a corresponding spectrum, a pattern similar to that in Figure 5A is seen. Figure 5C-a with a corresponding graph shows that suspended matter is composed of diatoms; spherical structures are also observed. The spectrum shows that material is composed of Si, Al, O, Fe, K, Ca, Na, Mg (Fig. 5C-c). Diatoms are composed of Si-oxide, while spherical structures with corresponding spectrum have a mi-croporous structure and are composed of Fe, Si, O, Ca. Fig. 5D with corresponding spectrum represents suspended matter composed of Si, O, Fe, Mg, Na, K, Ca with traces of Zn and Fe (Fig. 5D). Fibrous and microporous particles that were observed and analysed under higher magnification are composed of Al, Si, Fe and P, which could be attributed to anthropogenic input in the river system.
Fig. 5A. SEM micrograph and corresponding EDXS analysis/spectrum of samples of suspended matter at Pesnica 1.
.	* V
*—. i '	<*.*- *
g' ' rt» " •	^ . ' C ,
; '	' ,r
'	, -
t	• » ^
•*■ > ^ *	; j
?■ --ft/	. i^H
' * -v v \ ^
••> - •• *	- ••••
X 250	20 . OkV LEI
Fig. 5B. SEM micrograph and corresponding EDXS analysis/spectrum of samples of suspended matter at Pesnica 2 (above reservoir).
Hydrogeochemical and isotopic characterization of Pesnica River, Slovenia	189
Fig. 5C-a. SEM micrograph and corresponding EDXS analysis/spectrum of samples of suspended matter at Pesnica catchment area from location Pesnica 3 (below the reservoir, retention).
Fig. 5C-b. SEM micrograph and corresponding EDXS analysis/spectrum of samples of suspended matter at Pesnica catchment area from location Pesnica 3 (below the reservoir, retention).
Fig. 5C-c. SEM micrograph and corresponding EDXS analysis/spectrum of samples of suspended matter at Pesnica catchment area from location Pesnica 3 (below the reservoir, retention).
190 T. KANDUČ, Z. SAMARDŽIJA, N. MORI, A. JEREBIC, I. LEVAČIC, M. KRAČUN, J. A. ROBINSON, S. ŽIGON, Ž. BLAŽEKA & D. KOCMAN
Fig. 5D. SEM micrograph and corresponding EDXS analysis/spectrum of samples of suspended matter at Pesnica catchment area from location Pesnica 4 (drainage channel Biš).
Conclusion
Pesnica River water solute chemistry is dominated by HCO3-, Ca2+ and Mg2+. Seasonal (spring, summer and autumn) concentrations of total alkalinity ranged from 1.24 to 6.69 mM. From a geochemical point of view, we found out that dissolution of clastic material influences the geochemical composition of river water. Most locations of surface water were near equilibrium with respect to calcite (SI , .. from -0.98 to
j.	v calcite
0.54) and dolomite (SI, , .t from -2.61 to 0.48)
'	v dolomite	'
in all sampling seasons or, otherwise, generally undersaturated. According to the calculated pCO2 values, surface water is a source of CO2 to the atmosphere. pCO2 in Pesnica catchment ranged from 15488 to 316228 ppm. Alkalinity and isotopic composition of ^13Cdic values were used to trace carbon cycling in surface waters. Seasonal ^13CDIC values (from February 2016 to September 2016) ranged from -14.8 to -4.2 %0 indicating different biogeochemical processes in Pesnica River catchment. ^13CDIC values around and below -12 %0 indicate that the most important process in Pesnica River catchment in our study is degradation of organic matter. At a location below the reservoir Pernica lake (Pesnica 3 location), higher ^13CDIC values have been observed since April 2016. The highest value of ^13CDIC of -4.2 %0 was observed in June 2016. This was due to the presence of large amounts of algae from phytoplankton and phytobenthos in reservoir upstream of the sampling location, as confirmed using SEM microscopy. £13Cpoc ranged from -29.5 to -27.6 %o and indicated different degree of
suspended material at investigated locations.
Microscopic SEM/EDXS analyses revealed that at Pesnica 1 and 2 detrital clastic material prevails in the samples, while at location Pesnica 3 phytoplantonic and/or phytobenthic algae prevail, presumably originating from Pernica Lake, while at location near drainage channel Biš, particles originating from wetland plants are observed in addition to degraded organic material.
These first results on suspended matter in the River Pesnica catchment provide first data, that will help in evaluation of anthropogenic impacts, especially in relation to future impacts and for long term planning for sustainable management (e.g. regulation of river channel, deforestation and agriculture activity). The application of stable isotopes in river management studies is useful tracer for determining the origin of dissolved solutes and evaluation of biogeochemical processes, and has been widely used in other catchments in Slovenia (e.g. Sava River, Kamniška Bistrica River, Krka River, Idrijca River).
Acknowledgements
This research was conducted in the framework of the project L2-6778 and program research founding P2-0084 and P1-0143 funded by the Slovenian Research Agency (ARRS). Sincere thanks to Roger Pain for improving the English of the manuscript.
Hydrogeochemical and isotopic characterization of Pesnica River, Slovenia
191
References
Barth, J.A.C., Cronin, A.A., Dunlop, J. & Kalin, R.M. 2003: Influence of carbonates on the riverine carbon cycle in an anthropogenically dominated catchment basin: evidence from major elements and stable carbon isotopes in the Lagan River (N. Ireland). Chemical Geology, 200: 203-216.
Clesceri, L.S., Greenberg, A.E. & Eaton, A.D. 1998. Standard methods for the examination of water and wastewater. 20 th edition. APHA, AWWA, WEF, Baltimore.
Corine Land Cover 2012: http://gis.arso.gov.si/ geoportal/catalog/search/resource/details
Devol, A.H. & Hedges, J.I. 2001: Organic matter nutrients in the mainstem Amazon River. In: McClain, M.E., Victoria, R.L. & Richey, J. E. (eds.): The biogeochemistry of the Amazon Basin. Oxford University Press, Oxford: 275-306.
Gaillardet, J., Dupre, B., Louvat, P. & Allegre, C.J. 1999: Global silicate weathering and CO2 consumption rates deduced from the chemistry of large rivers. Chemical geology, 159: 3-30.
Hedges, J.I. 1992: Global biogeochemical cycle: progress and problem. Marine chemistry, 39: 67-93.
Hope, D., Billet, M.F. & Cresser, M.S. 1994: A review of the export of carbon in river water: fluxes and processes. Environmental pollution, 84: 301-304.
Kanduč, T. 2006: Hydrogeochemical characteristics and carbon cycling in the River Sava catchment in Slovenia. Ph.D. Thesis, University of Ljubljana (in Slovene), Ljubljana: 141 p.
Kanduč, T. 2010: Characterization of suspended matter in river systems: river Sava in Slovenia case study = Značilnosti suspendirane snovi v rečnih sistemih: primer reke Save v Sloveniji. Geologija, 54/1: 55-66, doi:10.5474/ geologija.2011.004.
Kanduč, T., Szramek, K., ogrinc, N. & Walter, L.M. 2007: Origin and cycling of riverine inorganic carbon in the Sava River watershed (Slovenia) inferred from major solutes and stable carbon isotopes. Biogeochemistry, 86: 137-154.
Kanduč T. & ogrinc N. 2007: Hydrogeochemical characteristics of the River Sava watershed in Slovenia = Hidrogeokemične značilnosti porečja reke Save v Sloveniji. Geologija, 50/1: 157-177, doi:10.5474/geologija.2 007.013.
Kanduč, T., Kocman D. & ogrinc N. 2008: Hydrogeochemical and stable isotope characteristics of the River Idrijca (Slovenia), the
boundary watershed between the Adriatic and Black Seas. Aquatic Geochemistry, 14: 239-262.
Kanduč, T., Mori, N. Kocman D., Stibilj V. & Grassa, F. 2012: Hydrogeochemistry of alpine springs from North Slovenia: insights from stable isotopes. Chemical geology, 300-301: 40-45.
Kanduč, T., Sturm Burnik, M. & McIntosh, J. 2013: Chemical dynamics and evaluation of biogeochemical processes in Alpine River Kamniška Bistrica, North Slovenia, Aquatic Geochemistry, 19: 323-346, doi:10.1007/ s10498-013-9197-4.
Karim, A. & Veizer, I. 2000: Weathering process in the Indus River basin: implications from river in carbon, sulphur, oxygen and strontium isotopes. Chemical geology, 120: 133-177.
Kobold, M. 2012: Hidrološki letopis Slovenije 2009, Ljubljana, Ministrstvo za kmetijstvo in okolje, Agencija Republike Slovenije za okolje.
Levin, I., Kromer, B, Wagenback, D. & Minnich, K.O. 1987: Carbon isotope measurements of atmospheric CO2 at a coastal station in antartica. Tellus 39, B B: 89-95.
Leitinger, V. 2012: Hidrogeografske značilnosti porečja reke Pesnice s pudarkom na vodnogospodarskih ureditvah. Diplomsklo delo, Univerza na Primorskem, Koper: 92 p.
Li, Y., Zhang, H., Tu, C., Fu, C., Xue, Y. & Luo, Y. 2016: Sources and fate of organic carbon and nitrogen from land to ocean: Identified by coupling stable isotopes with C/N ratio. Estuarine, Coastal and Shelf Science 181: 114-122.
Markoli, T. 2012. Kvaliteta vode v akumulacijskih jezerih na reki Pesnici, Diplomsko delo, Univerza v Mariboru, Maribor: 87 p.
Mayorga, E. Aufdemkampe, A.K., Masiello, C.A., Krusche, A.V, Hedges, J.I, Quay, P.D, Richey, J.E. & Brown T.A. 2005: Young organic matter as a source of carbon dioxide outgassing from Amazon rivers. Nature, 436: 538-541. doi:10.1038/nature03880.
Mioč, P. & Žnidarčič, M. 1989: Osnovna geološka karta 1:100.000, Tolmač za lista Maribor in Leibniz-L33-56, L33-44, Zvezni geološki zavod, Beograd: 60 p.
Meybeck, M. 1981: River transport of organic carbon to the ocean. V Flux of organic carbon to the Oceans, Edited by Likens G.E., Mackenzie, F.T., Richey J.E., Seddell, J.R. Turekian K.K.U.S.D.O.E. C0NF-8009140: str. 219-269.
192 T. KANDUČ, Z. SAMARDŽIJA, N. MORI, A. JEREBIC, I. LEVAČIC, M. KRAČUN, J. A. ROBINSON, S. ŽIGON, Ž. BLAŽEKA & D. KOCMAN
Meybeck, M. 1987: Global chemical weathering of surficial rocks estimated from river dissolved loads. American Journal of Science, 287: 401-428.
Mook, W.G., Bommerson, J.C. & Staverman, W.H. 1974: Carbon isotope fractionation between dissolved bicarbonate and gaseous carbon dioxide. Earth Planet Sci Lett, 22: 169-176.
Parkhurst, D.L. & Appello 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. Water-Resources Investigation Report 99-4259.
Perko, D. & Orazen Adamic, M. (eds.) 2001: Slovenija: pokrajina in ljudje. Ljubljana, Mladinska knjiga.
Popescu, R., Mimmo, T., Dinca, O.R., Carici, C., Costinel, D., Sandru, C., Ionete, R.E., Stefanescu, I. & Axente, D. 2015: Using stable isotopes in tracing contaminant sources in an industrial area: A case study on the hydrolog-ical basin of the Olt River, Romania.
Roach, K.A., Rodriguez, M.A., Paradis, Y. & Cabana, G. 2016: Controls of longitudinal variation in c>13Cdic in rivers: A global meta-analy-sis. Aquatic ecology in press.
Romanek, C.S., Grossman, E.L. & Morse, J.W. 1992: Carbon isotopic fractionation in synthetic ar-agonite and calcite: effects temperature and precipitation rate. Geochim. Cosmochim. Acta, 46: 419-430.
Sameer, K.T., Santosh, K.R. Sukesh, K.B., Aril, K.G. & Negi, H. 2016: Stable isotopes (¿13C
¿D, c>18O) and geochemical characteristics of geothermal springs of Landakh and Himahal (India): Evidence for CO2 discharge in northwest Himalaya. Geothermics 64: 314-330.
Schuster, P.F. & Reddy, M.M. 2001: Particulate Carbon (PC) and Particulate Nitrogen (PN). In: Water and Sediment Quality in the Yukon River Basin, Alaska, During Water Year 2001, Open-File Report 03-427, National Research Program, USGS, available online: http:// pubs.usgs.gov/of /2003/ofr03427/.
Statistical office of Republic of Slovenia (Statistični urad RS) 2002: Reke dolge nad 25 km in njihova padavinska območja.
Szramek, K., McIntosh, J.C., Williams, E.L., Kanduč, T., Ogrinc, N. & Walter, L.M. 2007: Relative weathering intensity of cal-cite versus dolomite in carbonate-bearing temperature zone watersheds: carbonate geochemistry and fluxes from catchment within the St. Lawrence and Danube river basin. Geochem. Geophys. Geosys. 8/4:1-26, doi:10.1029/2006GC001337.
Urbanek, J. 1994: Pesnica, In Encyclopedia Slovenia No. 8, pp. 314, Mladinska knjiga, Ljubljana,
Zavadlav, S., Kanduč, T., McIntosh, J. & Lojen, S. 2013: Isotopic and chemical constraints on the biogeochemistry of dissolved inorganic carbon and chemical weathering in the karst watershed of Krka River (Slovenia). Aquatic geochemistry, 19: 3:209-230, doi:10.1007/ s10498-013-9188-5.