UDK/UDC: 550.47:556.51(497.4) Izvirni znanstveni članek - Original scientific paper Prejeto/Received: 07.04.2015 Sprejeto/Accepted: 01.10.2015 SPREMLJANJE DINAMIKE VNOSA IN SPIRANJA NITRATNEGA DUŠIKA NA POREČJU POTOKA PADEŽ TRACING NITRATE NITROGEN INPUT-OUTPUT DYNAMICS AT THE 'Fakulteta za gradbeništvo in geodezijo, Univerza v Ljubljani, Hajdrihova 28, 1000 Ljubljana Izvleček Občutljivost ekosistemov na spremembe količin dušika v okolju je velika. V želji, da jih zaščitimo pred preobremenjenostjo s hranili in jim omogočimo zdrav razvoj, je poznavanje masne bilance dušika ključnega pomena. Spremembe v količini dušika se, poleg naraščajočih antropogenih vnosov dušikovih spojin v okolje, pojavijo zaradi delovanja meteoroloških, hidroloških, geoloških in bioloških dejavnikov oziroma biogeokemijskega kroženja med biotopom in biocenozo. S pomočjo podatkov o koncentracijah nitratnega dušika v vodotoku in padavinah ter podatkov o količini padavin in odtoku s prispevnega območja potoka Padež smo analizirali mesečno dinamiko vnosa in spiranja nitratnega dušika, ki je podana kot razlika med količino vnesenega in spranega nitratnega dušika na hektar prispevne površine. V časovni dinamiki je opazen sezonski vpliv. Porečje nekaj časa deluje kot izvor in nekaj časa kot ponor nitratnega dušika, odvisno od hidroloških razmer in potencialne količine mobilnega nitratnega dušika v sistemu tal, ki se nadalje spira v vodotok, ter njegove koncentracije v padavinah. Ključne besede: nitratni dušik, biogeokemija, hidrološke razmere, porečje Padeža. Abstract The sensitivity of ecosystems to changes in the amounts of nitrogen in the environment is large. In order to protect them against overloads of nutrients and to allow their healthy development, the knowledge of the nitrogen mass balance is crucial. Differences in the amount of nitrogen occur because of increasing inputs from anthropogenic sources, but they are also due to meteorological, hydrological, geological and biological factors or, in other words, to biogeochemical circulation between the biotope and organisms. The data on stream water nitrate concentrations and nitrate concentrations in precipitation together with measurements of precipitations and Padež stream discharge were used to analyze nitrate-nitrogen input-output dynamics. The analysis is presented as a difference between inputs and outputs of nitrate-nitrogen per hectare of catchment area. The impact of seasonality could be observed. Watershed behaves as a sink or as a source of nitrate-nitrogen depending on the hydrological conditions and the amount of potentially available nitrate-nitrogen in soil which is flushed towards the stream and its concentrations in precipitation. Keywords: nitrate nitrogen, biogeochemistry, hydrological conditions, the Padež stream watershed. *Stik / Correspondence: katarina.cesnik@gmail.com © Česnik K., Rusjan, S.;Vsebina tega članka se sme uporabljati v skladu s pogoji licence CreativeCommons Priznanje avtorstva -Nekomercialno - Deljenje pod enakimi pogoji 4.0. © Česnik K., Rusjan, S.; This is an open-access article distributed under the terms of the Creative Commons Attribution - Non Commercial - Share Alike 4.0 Licence. PADEŽ STREAM WATERSHED v 1 * Katarina Česnik ' , Simon Rusjan 1 1. Uvod Čista in neoporečna voda je ključna za zdravje ljudi in delovanje naravnih ekosistemov, zato je varstvo kakovosti voda eden od temeljev evropske okoljske politike. To potrjuje tudi dejstvo, da je bila že leta 1991 sprejeta direktiva o nitratih, ki je eden izmed prvih zakonodajnih aktov v EU, katerih namen je nadzorovati onesnaževanje voda in izboljševati njihovo kakovost. Če želimo sonaravno in trajnostno gospodariti z vodotoki, je potrebno poznati procese, ki se odvijajo na celotnem prispevne območju (porečju) vodotoka (Fogg in Wells, 1998). Procesi biogeokemijskih pretvorb in transporta hranil na nivoju posameznih prostorskih in ekosistemskih enot so že vrsto let predmet raziskav znanstvenikov z različnih znanstvenih področij. Ob tem je dušik zaradi splošne okoljske razširjenosti in kritične vloge kot dejavnik biološke produkcije bil, in prav gotovo bo tudi v prihodnje, deležen posebne pozornosti (Vitousek et al., 1997). V prispevku smo pozornost posvetili kroženju nitratnega dušika znotraj eksperimentalnega porečja potoka Padež. Pridobili smo vpogled v bilanco »naravnega ozadja« količin nitratnega dušika in ugotovili, v katerih primerih in kdaj porečje deluje kot potencialni izvor in kdaj kot ponor nitratnega dušika. Za obdobje dveh let (2006 in 2007), smo imeli na voljo podatke o hidrometeoroloških razmerah in kemizmu vode na gozdnatem prispevnem območju potoka Padež. S pomočjo pridobljenih podatkov smo izračunali masno bilanco vnesenega in spranega nitratnega dušika v tem času in opredelili letno, sezonsko in mesečno masno bilanco nitratnega dušika v gozdnatem porečju Padeža. 2. Metode 2.1 Kroženje dušika Kroženje dušika v naravi (slika 1) opišemo z dušikovim ciklom, ki predstavlja enega izmed pomembnejših biogeokemijskih ciklov v naravi. Glavni biološki procesi dušika v biosferi so (Stopar, 2013): • biološka fiksacija dušika, • mineralizacija dušika (amonifikacija), • nitrifikacija (heterotrofna in autotrofna), • denitrifikacija, • asimilacija dušika (imobilizacija). Slika 1: Shematični prikaz kroženja dušika (prilagojeno po Ward, 2005). Figure 1: Schematic diagram of the nitrogen cycle (adapted after Ward, 2005). Biološka fiksacija dušika omogoča vezavo atmosferskega dušika v anorgansko obliko, ki jo rastline lahko asimilirajo. Običajno je glavni produkt mineralizacije amonij, zato temu procesu rečemo tudi amonifikacija. Gonilno silo mineralizacije dušika predstavlja mineralizacija ogljika. Pri razgradnji ogljika mikroorganizmi razgradijo organsko snov, pri tem pa se sprostijo tudi drugi elementi, del teh se porabi za sintezo nove biomase. Mineralizacije dušika so sposobni skoraj vsi mikroorganizmi, zato je potencialna mineralizacija vedno prisotna, njena hitrost v tleh pa je odvisna od (Stopar, 2013): vlage, prezračenosti, pH vrednosti, temperature, anorganskih hranil, talnih koloidov. Razpadajoča organska snov se v procesu mineralizacije oziroma amonifikacije spremeni v amonijev ion NH4+, ki predstavlja izhodiščno spojino za nadaljnje pretvorbe. V vodnem okolju ob prisotnosti kisika poteče mikrobni proces nitrifikacija. Pri tem procesu pride do pretvorbe amonijevega dušika v oksidirane oblike. Mineralizacija organskih dušikovih spojin ter nadaljnja nitrifikacija se najintenzivneje odvijata v zgornjih talnih horizontih, posledica tega je, da z globino vsebnost dušikovih spojin upada (Rusjan, 2008). Velja, da se poraba dušika s strani heterotrofnih mikrobov ob višjih vrednostih razmerja C/N poveča. V kolikor je limitirajoči faktor dušik, pride do imobilizacije dušika v tleh in zmanjšanja razpoložljivih količin za nadaljnjo nitrifikacijo in odtok s pronicajočo vodo. Če je limitirajoči faktor ogljik, pride do neto mineralizacije. Nadaljnja nitrifikacija poteka v dveh stopnjah, v prvi stopnji kot produkt nastane nitrit NO2-, v drugi nitrat NO3-. Nitrat je hranilo za vse primarne proizvajalce, ki ga kot vir dušika asimilirajo in vežejo v organske snovi. Povečana asimilacija v vodnem sistemu pospeši rast alg in makrofitov ter potencialno evtrofnost. Pride lahko tudi do obratnega anaerobnega mikrobnega procesa denitrifikacije (redukcije), v katerem bakterije pretvarjajo organske dušikove spojine nazaj v elementarni dušik (N2) (Ward, 2005). Poti kroženja dušika v vodotoku so prikazane na sliki 2. Kroženje dušika je sezonsko pogojeno glede na stopnjo aktivnosti vegetacije. V obdobju neaktivnosti vegetacije pozno jeseni, pozimi in zgodaj spomladi so bile v številnih gozdnih ekosistemih v vodotokih opažene višje koncentracije raztopljenega neorganskega dušika (večinoma nitrata). Te dodatne količine predstavljajo presežek neorganskega dušika za gozdni ekosistem in se akumulirajo v sistemu gozdnih tal, dokler hidrološke razmere ne omogočijo njihove mobilizacije (Rusjan, 2008). Največ dušika se v porečje vnese s padavinami (mokra depozicija - odlaganje), manjši del pa s suho depozicijo nekaterih dušikovih spojin. Povprečne ocene količine vnesenega dušika iz atmosfere se gibljejo od 1 do 5 kg-N/ha/leto na območjih brez industrijskega onesnaženja do več kot 25 kg-N/ha/leto na intenzivno kmetijskih in industrijskih območjih v Evropi (Binkley et al., 2000). Potencialno zelo velik vir vnosa dušika v porečje je lahko biološka fiksacija dušika (Binkley et al., 1994). Na biološko fiksacijo dušika močno vplivajo okoljske razmere (biološka fiksacijska sposobnost nekaterih rastlinskih vrst) in zaloge ogljikovega substrata v tleh. Stopnja biološko fiksiranega dušika je letno večinoma manjša od nekaj kilogramov na hektar površine (Mitchell et al., 1996). V prisotnosti simbiotskega razmerja med bakterijami, sposobnimi fiksacije dušika, in rastlinskimi vrstami (od drevesnih vrst npr. jelša) pa lahko količina fiksiranega atmosferskega dušika preseže 100 kg/ha/leto (Binkley et al., 1994). Posredna nitrifikacija Denitrifikacija Slika 2: Kroženje dušika v vodotoku (povzeto po Heggie et al., 1999). Figure 2: Nitrogen cycle in surface waters (adapted from Haggie et al., 1999). Slika 3: Merska oprema na prispevnem območju potoka Padež (Rusjan, 2008). Figure 3: Monitoring system at Padež stream watershed. (Rusjan, 2008). Količine dušika, ki se z odtokom sperejo z gozdnatih prispevnih območij in dejansko predstavljajo izgubo ali presežek hranila za gozdni ekosistem, so zelo različne v odvisnosti od strukture gozda in klimatskih razmer, večinoma pa se gibljejo med 5 in 10 kg-N/ha/leto (Rusjan, 2008). V Sloveniji se višje vrednosti spranega dušika pojavljajo zelo prostorsko razdrobljeno v zelo majhnem obsegu. Podatki kažejo, da prihaja do večjih količin spranega dušika na manj kot 1% gozdnih območij Slovenije (Eler et al., 2010). 2.2 Opis porečja Padež Porečje Padež, ki obsega 42,1 km2, je del porečja reke Reke. Nahaja se na hribovitem območju Brkinov, ki predstavljajo klimatsko prehodno območje med primorsko in celinsko Slovenijo. Iz notranjosti na območju piha burja z morske strani pa toplejši veter jugo. Povprečna gostota poselitve znaša na območju Brkinov manj kot 5 ljudi/km2. Posledično več kot 80% prispevnega območja Padež prekrivajo gozdovi, preostali del pripada različnim kmetijskim površinam. Znotraj porečja se nahajata glavna vodotoka Padež in Suhorka s pritoki (slika 3). 2.3 Metode zajema podatkov Količine padavin na prispevnem območju Padež so bile izmerjene s pomočjo avtomatskih dežemerov ONSET RG2-M. Merjenje pretokov je potekalo s pomočjo dveh različnih inštrumentov, srednji in višji pretoki so bili merjeni z merilcem pretočnih hitrosti Flowtracker HandheldADV, ki deluje na principu Doppleijevega pojava, medtem ko je merjenje nizkih pretokov potekalo s pomočjo prenosnega merilca Flo-Tracker, ki pretok izračuna na podlagi sprememb v koncentraciji sledila v vodotoku. Merjenje fizikalnih in kemijskih parametrov vode je potekalo z minisondo Hydrolab DataSonder 4. Sonda je opremljena s senzorji za merjenje temperature in pH vrednosti vode, elektroprevodnosti, količine raztopljenega kisika, redoks potenciala in koncentracije nitratnega dušika. Za natančnost podatkov je potrebno sondo redno umerjati in vzdrževati. Kontinuirane meritve so potekale na vodomerni postaji Padež (slika 3), zajem podatkov se je izvajal s 15-, 30-minutnim in urnim časovnim korakom, izračuni količin spranega nitratnega dušika so bili opravljeni v urnem časovnem koraku. Za kontrolo delovanja multiparametrske sonde so služile laboratorijske analize, ki jih je v letih 2005 in 2006 izvajal Zavod za zdravstveno varstvo Maribor (Rusjan, 2008). S programskim orodjem za rudarjenje po podatkih WEKA (Waikato Environment for Knowledge Analysis) so bile na podlagi izmerjenih koncentracij modelirane koncentracije nitratnega dušika v [mg/l-N] v času, ko le-te niso bile izmerjene. Meritve kemijske sestave padavin so v okviru osnovne merilne mreže (meritve Agencije Republike Slovenije za okolje - ARSO) potekale na merilnem mestu Škocjan, ki se nahaja na nadmorski višini 420 m, v relativno čistem in podeželskem okolju. Zbiranje padavin je potekalo s pomočjo le-mokrih (wet-only) padavinskih vzorčevalnikov. Padavine se zberejo v posebni posodi, ki ne spremeni kemijske sestave padavin in daje zanesljive podatke o njihovi dnevni količini. Vzorči se torej samo mokro usedlino (kapljice padavin) zato podatki o vnosu suhih usedlin v analizi niso bili upoštevani. Posamezni dnevni vzorci padavin so bili združeni v tedenske vzorce (dnevni vzorci od ponedeljka do nedelje iste serije) in so bili nato kemijsko analizirani (Bolte, 2007). 2.4 Računanje masne bilance nitratnega dušika Na podlagi časovnega spremljanja količine vnosa in iznosa dušika v določenem okoljskem (prostorskem) okvirju lahko izračunamo masno bilanco dušika. Izražamo jo v kilogramih dušika na hektar prispevne površine (kg N/ha). Pri računanju sezonskih zalog dušika v ekosistemu upoštevamo dva kriterija. Kot prvo določimo usmerjenost izmenjave snovi, torej ali se dušik v ekosistem vnaša ali se iz njega izgublja, kot drugo določimo magnitudo izmenjave oziroma spremembo v količini med neto vnosom in izgubo dušika v ekosistemu (Likens in Bormann, 1995). Masna bilanca dušika v posameznih komponentah hidrološkega kroga letno, predvsem pa sezonsko, močno niha, zato so za analizo spremenljivosti potrebna dolgoletna merjenja in spremljanja hidrološke bilance in vsebnosti dušika (Likens in Bormann, 1995). Z vidika vpliva hidroloških mehanizmov na kroženje dušika so najbolj zanimive tiste kemijske oblike, ki so z vodo najbolj mobilne, obenem pa je njihova količina prevladujoča. Če zaradi odsotnosti antropogenih vplivov in rabe tal na porečju Padeža zanemarimo vnos nitratnega dušika zaradi delovanja človeka, se v obravnavan ekosistem, v našem primeru je to porečje, nitratni dušik v prispevno območje vnese iz ozračja preko padavin, spiranje iz sistema pa poteka pod vplivom tvorbe padavinskega odtoka. Dodatno na spiranje nitratnega dušika iz sistema močno vplivajo intenzivni biogeokemijski procesi pretvorb različnih kemijskih oblik dušika (procesi fiksacije, mineralizacije, denitrifikacije) v odvisnosti od fizikalnih okoljskih razmer, sestave rastlinstva, geološke in pedološke sestave tal (Aber et al., 1991; Band et al., 2001). Teoretično nam razlike med sezonskim vnosom nitratnega dušika (I -Input) in spiranjem (O - Output) povedo, ali se kvantitativno dušik v obliki nitrata v ekosistemu akumulira (I>O), se iz njega izgublja (IO). Obilnejše padavine v tem času lahko povzročijo, da se nitratni dušik kot produkt mineralizacije v gozdnih tleh prične iz porečja intenzivneje spirati, kvantitativno se nitratni dušik iz porečja izgublja. Kot primer lahko navedemo dogajanje avgusta 2006. Kot posledica podnebne spremenljivosti se napoveduje višanje povprečnih temperatur zraka, s tem pa podaljševanje biogeokemijsko najaktivnejših obdobij (Glavan et al., 2015). Ob tem lahko v prisotnosti hidroloških mehanizmov (padavinski dogodki z visoko intenziteto) pričakujemo potencialno povečanje količin nitratnega dušika, ki se z gozdnatih prispevnih območij sprostijo v vodotoke. Dodatno lahko zaradi naraščajočega antropogenega vnosa neorganskih dušikovih spojin v ozračje (povečevanje porabe umetnih gnojil, fosilnih goriv) pričakujemo povečan atmosferski vnos nitratnega dušika. Povečane količine spranega nitratnega dušika po vodotokih v vegetacijsko aktivnih obdobjih pomenijo dodatne količine hranil za mikroorganizme in vegetacijo, ki se nahaja dolvodno. Poveča se možnost evtrofikacije vodnih teles. V odsotnosti podrobne analize časovne spremenljivosti vsebnosti nitratnega dušika v tleh, kjer potekajo intenzivni biogeokemijski procesi tvorbe in odvzema nitratnega dušika, ponujajo predstavljeni rezultati raziskave zgolj grob vpogled v časovno spremenljivost dinamike vnosa in spiranja nitratnega dušika iz porečja potoka Padež. Ob tem se nakazujejo vzorci sezonske spremenljivosti dinamike vnosa in spiranja nitratnega dušika, ki so močno pogojeni s spremenljivostjo hidroloških razmer. Rezultati meritev nadalje kažejo, da lahko časovno zvezno spremljanje koncentracije nitratnega dušika v vodotoku v času padavinskih dogodkov veliko pripomore k boljši kvantifikaciji količin spranega nitratnega dušika z manjših prispevnih območij, kakršno je porečje potoka Padež. Za natančnejši vpogled v spremenljivost masne bilance dušika na območju eksperimentalnega porečja Padež bi bila potrebna dolgotrajnejša in večletna spremljanja koncentracij različnih oblik dušikovih spojin v vodotoku, atmosferi in tleh. Viri Aber, J.D., Melillo, J.M., Nadelhoffer, K.J., Pastor, J., Boone, R.D. (1991). Factors controlling nitrogen cycling and nitrogen saturation in northern temperate forest ecosystems, Ecological Application 1, 303-315. Band, L.E., Tague, C.L., Groffman, P., Belt, K. (2001). Forest ecosystem processes at the watershed scale: hydrological and ecological controls of nitrogen export, Hydrological Processes 15(10), 2013-2028. Binkley, D., Cromack, K., Baker, D. (1994). Nitrogen Fixation by Red Alder; Biology, Rates and Controls. V: Hibbs, D.E., DeBell, D.S., Tarrant, R., The Biology & Management of Red Alder, Oregon State University Press, Corvallis: 57-72. Binkley, D., Son, Y., Valentine, W.D. (2000). Do Forests Receive Occult Inputs of Nitrogen?, Ecosystems 3(4), 321-331. Bolte, T. (2007). Kakovost zraka v Sloveniji v letu 2006. ARSO: 106 str. Eler, K., Kobal, M., Simončič, P. (2010). Kritični vnos in prekomerne obremenitve ekosistemov z dušikom in žveplom. http://kazalci.arso.gov.si/print?ind id=576&lang id=30 2 (Pridobljeno 15. 03. 2014). Fogg, J., Wells, G. (1998). Stream Corridor Restoration, Principles, Processes, and Practices, U.S. Department of Agriculture, U.S. Environmental Protection Agency, Tennessee Valley Authority, Federal Emergency Management Agency, U.S. Department of Commerce, U.S. Department of Housing and Urban Development, U.S. Department on the Interior: 536 str. Glavan, M., Ceglar, A., Pintar, M. (2015). Assessing the impacts of climate change on water quantity and quality modelling in small Slovenian Mediterranean catchment - lesson for policy and decision makers, Hydrological Processes 29(14), 3124- 3144. Heggie, D. T., Skyring, G. W., Orchardo, J., Longmore, A. R., Nicholson, G. J., and Berelson, W. M. (1999). Denitrification and denitrifying efficiencies in sediments of Port Phillip Bay: direct determinations of biogenic N2 and N-metabolite fluxes with implications for water quality, Marine Freshwater Research 50, 589-596. Likens, G.E., Bormann, F.H. (1995). Biogeochemistry of a Forested Ecosystem. 2nd edition. New York, Springer-Verlag: 159 str. Mitchell, M.J., Driscoll, C.T., Kahl, J.S., Likens, G.E., Murdoch, P.S., Pardo, L.H. (1996). Climatic Control of Nitrate Loss from Forested Watersheds in the Northeast United States, New York, Enviromental science and technology 8, 2609-2612 Rusjan, S. (2008): Hidrološke kontrole sproščanja hranil v porečjih (Hydrological controls of nutrient mobilization in watersheds). Unpublished Doctoral Thesis, Univerza v Ljubljani, FGG: 206 p. (in Slovenian). Simončič, S., Kutner, L., Kovač, M., Skudnik, M., Žlindra, D., Ferlan, M., Levanič, T., Verlič, A. (2011). Raziskovalna ploskev Brdo. Gozdarski inštitut Slovenije, Ljubljana: 2 str. Stopar, D. (2013). Kroženje dušika. http://web.bf.uni-li.si/zt/mikro/homepage /krozenieDUSIK.pdf/ (Pridobljeno 13. 09. 2013). Vitousek, P.M., Aber, J.D., Howarth, R.W., Likens, G.E., Matson, P.A., Schindler, D.W., Schlesinger, W.H., Tilman, D.G. (1997). Human alteration of the global nitrogen cycle: Sources and consequences. Ecological Application 7, 737-750. Ward, B. (2005). The Global Nitrogen Cycle. V: Knoll, A.H., Canfield, D.E., Konhauser, K.O., Fundamentals of Geobiology, Blackwell publishing; 36-48.