GDK 56:228__015.5 Prispelo/ÄeceiVerf: 20. 07. 1999 , . . Sprejeto//icceprerf: 03 11 1999 Izvimi znanstven, članek Original scientific paper ANALIZA SPREMENJENOSTI KOLIČINSKE IN VRSTNE SESTAVE GOZDOV PO DVEH RAZLIČNIH METODAH TER NJUNA PRIMERJAVA Janez KRČ* Izvleček leZrLT^'^rt spremenjenost gozdnih sestojev po dveh metodah. S prvo ocenjujemo spremembe vrstne sestave rastlmskih skupnosti, druga metoda pa ugotavlja spremembe sestojev po višini in vrsti lesnih zalog Za .zračun ocene spremenjenosti na ravm osnovnih enot popisa gozdov smo uporabili modelno optimalno in dejansko strukturo lesnih zalog po drevesnih vrstah v gozdnih združbah, popisanih v gozdnih odsekih Zaradi podatkov potrebnih za izračun, smo v analizo vključili le odseke, v katerih so ugotovili samo eno gozdno zdruzbo. Prikazujemo rezultate anahze na ravni gozdnogospodarskih območij in gozdnih združb. Rezultate na višj.h ravneh velikosti gozdnih površin smo dobili s tehtano aritmetično sredino rezultatov osnovnih izračunov oLh metod Rezultate ocen spremenjenosti smo prikazali tudi prostorsko in naredili primerjavo Ključne besede: vrstna sestava, spremenjenost gozdnih sestojev, računalniški model, kazalec vrstne spremenjenosti ANALYSIS OF FOREST QUANTITIES AND SPECIES COMPOSITION ALTERATION USING TWO DIFFERENT METHODS WITH COMPARISONS Abstract The article analyses composition alterations in forest stands using two methods. In the first composition alterations m types of plant community are estimated and in the second, composition alterations in stands according to amounts and types of wood stock are established. To assess composition alterations at the level of basic units of forest inventory, the optimal and actual structure of wood stock was used in models according to tree types m forest communities within forest tracts. The data needed for calculations required analysis only of tracts containing one forest community Results are shown at the level of forest management region and forest community. Results at higher levels of forest area were obtained using a weighted arithmetic mean of results from basic calculations at the level of forest tracts. The results of composition alteration assessment are also shown by area and comparisons were made of both methods used. Key words: composition type, forest stand composition alteration, computer model, type alteration indicator dr. asist., Biotehniška fakulteta. Oddelek za gozdarstvo in obnovljive gozdne vire. Večna pot 83 1000 Ljubljana, SVN VSEBINA CONTENTS 1 UVOD INTRODUCTION............................................................................................................213 2 DELOVNA METODA WORKING METHOD............................................................................................214 3 REZULTATI RESULTS...................................................................................218 4 POVZETEK.................................................................231 5 SUMMARY....................................................................................232 6 VIRI REFERENCES..................................................................................................................................234 7 PRILOGA APPENDIX................................................................................ 235 1 UVOD INTRODUCTION Spremenjenost rastlinskih skupnosti v gozdu je kazalec človekovega neposrednega ir posrednega poseganja oz. vplivanja na gozdni ekosistem v prostoru in času Vzroki spremenjenosti najprej izvirajo iz različnih sistemov gospodarjenja z gozdom (teorija največje zemljiške rente in gospodarjenja z uporabo metode starostnih razredov), načma obnavljanja sestojev (naravna in umetna obnova) ter najrazličnejših biotskih in abiotskih vplivov na gozdni prostor (imisije, antropogene spremembe klime, .„). Spremenjena vrstna sestava nam da informacijo, ki je rezultat primerjave dejanskega m potencialno naravnega stanja določene rastlinske skupnosti. Služi kot osnova za presojo in vrednotenje preteklih ter določanje novih ukrepov. Hkrati je indikator dogajanja v preteklosti. V veliki meri vpliva na izkoriščenost proizvodne sposobnosti rastišč, ekološko primernost in ekonomiko gospodarjenja z gozdom S tehnološkega vidika m z vidika načrtovanja gozdne proizvodnje so spremenjeni gozdovi delovišča z višjo stopnjo verjetnosti opravljanja del, saj je potrebno v njih bolj pogosto opravljati sanacije škod (snegolomi, vetrolomi, gradacije insektov, bolezm, m redne sečnje zaradi posredne m neposredne premene sestojev. Spremenjenost vrstne sestave gozdnih sestojev ima zato pomembno vlogo pri gozdnogospodarskem načrtovanju, ker bodo zahtevam po trajnem m večfunkcionalnem gospodarjenju z gozdom ustrezali le gozdovi z ustrezno, rastišču primemo vrstno drevesno sestavo in lesno zalogo. Ugotavljanje spremenjenoiti vrstne sestave je podatkovno obsežno ter raziskovalno zahtevno delo, saj moramo poznati dejanske in naravne strukture rastlinskih skupnosti. Če imamo na voljo metodo določanja spremenjenosti in zahtevane podatke na računalniškem mediju, lahko izračun za obsežne podatkovne zbirke opravi računalnik. Znane so številne metode, ki primerjajo sestave rastlinskih skupnosti. Sodobne metode vse bolj uporabljajo računalniško tehnologijo m geografske informacijske sisteme. V Ameriki (ZDA) razvijajo geografski pristop k načrtovanju biološke pestrosti (SCOTT 1997), da bi zagotoviti informacije o stanju neogroženih rastlinskih in živalskih vrst. Te informacije bodo namenjene strokovnjakom, ki upravljajo naravne vire, znanstvenikom m politikom za skladnejše odločitve oz. odločitve, ki bodo sprejete na višji ravni informacijske podpore. Računalniško obdelavo podatkov popisa gozdov smo opravili tudi v našem primeru, da bi ugotovili stopnjo spremenjenosti gozdov na državni ravni. 2 DELOVNA METODA WORKING METHOD 2.1 METODA OCENJEVANJA SPREMENJENOSTI GOZDOV PO METODI SPREMENJENOSTI VRSTNE SESTAVE RASTLINSKIH SKUPNOSTI METHOD OF ASSESSING FOREST COMPOSITION ALTERATIONS USING COMPOSITION ALTERATIONS ESf PLANT COMMUNITY TYPE Raziskava obravnava spremenjenost slovenskih gozdov po dveh razhčnih metodah. Prva, metoda I, je metoda ocenjevanja spremenjenosti vrstne sestave rastlinskih skupnosti (BONČINA / ROBIČ 1998): D h=m*- Is =indeks spremenjenosti i =število različnih drevesnih vrst Mj =modelna vrednost za deleže lesne zaloge posamezne drevesne vrste S; —dejanska vrednost za deleže lesne zaloge posamezne drevesne vrste Izdelali smo računalniški program, ki je obdelal podatke popisa gozdov iz 1. 1990. Prostorski okvir za izračun je bil odsek kot osnovna ureditvena enota baze popisa gozdov. Analizo smo naredili z uporabo razmerij med deleži drevesnih vrst v lesni zalogi odsekov. Deleže smo primerjali z modelnimi vrednostmi deležev lesne zaloge po gozdnih združbah (KRČ 1999). V analizo smo vključili le odseke, ki imajo samo eno gozdno združbo, so gozdna zemljišča in imamo zanje podatke o lesni zalogi po drevesnih vrstah. Takih odsekov je v bazi popisa gozdov 44878 s skupno površino 448926 ha oz. 43 % v popisu zajetih površin. Če so odseki razdeljeni na več gozdnih združb (največ tri), ne dajo potrebnih podatkov za uporabo metode ocenjevanja spremenjenosti vrstne sestave rastlinskih skupnosti, saj nimamo podatka o lesnih zalogah po drevesnih vrstah in gozdnih združbah. pač pa le o lesnih zalogah na ravni odsekov. Pomanjkanje podatkov bi lahko premostili z modelno porazdelitvijo lesnih zalog po gozdnih združbah v odsekih, vendar bi približek zagotovo močno vplival na rezultat, kar bi zameglilo dejansko stanje na različnih prostorskih ravneh. Natančnost prikazanih rezultatov je torej posledica metode izračuna in uporabe (točnosti) razpoložljivih podatkov. Rezultate analize smo predstavili tudi z izdelavo slike spremenjenosti gozdov na državni ravni (sliki 2 in 3). 2.2 METODA UGOTAVLJANJA SPREMENJENOSTI SESTOJEV PO VIŠINI IN VRSTI LESNIH ZALOG METHOD OF ESTABLISHING STAND COMPOSITION ALTERATIONS ACCORDING TO AMOUNT AND TYPE OF WOOD STOCK Imamo analizo spremenjenosti slovenskih gozdov po metodi primerjave višin lesnih zalog med modelno optimalnim ter dejanskim stanjem po podatkih popisa gozdov (metoda II) (KRČ 1999). Izračun spremenjenosti je bolj kompleksen v primerjavi z metodo I. Modelno višino lesne zaloge za drevesno vrsto izračunamo z uporabo vrednosti za prvi bonitetni razred donosnih tablic (ČOKL 1961). Vrednosti v donosnih tablicah so določene za čiste sestoje in enotne dobrote (bonitete) rastišč, kar pa v naravi ni prav pogosto. Zato tablične vrednosti za lesne zaloge sestojev prve bonitete prilagodimo modelni rastni sposobnosti drevesne vrste v gozdni združbi. Za prilagoditev višin modelnih lesnih zalog smo uporabili korekcijski faktor, ki je kvocient med rastiščnim koeficientom (RK) gozdne združbe in modelno privzetim najvišjim RK za čiste sestoje prvega bonitetnega razreda. Za prilagoditev višin lesnih zalog s pomočjo rastiščnih koeficientov smo se odločili iz dveh razlogov: (a) rastiščni koeficienti temeljijo na podmeni, da proizvodnja organske mase narašča od ekstremnih proti optimalnim rastiščem (boniteta rastišč) po krivuljah, ki jih lahko izrazimo v povprečju z normirano normalno krivuljo (KOŠIR 1975) m (b) rastiščni koeficienti so določeni po gozdnih združbah - torej razpolagamo z vrednostjo RK na državni ravni oz. na vseh ravneh popisnih enot. RK je relativna ocena proizvodne sposobnosti gozdnih rastišč in je ne moremo istovetiti s potencialno prirastno sposobnostjo gozdov v kubičnih metrih, so si pa te vrednosti zelo blizu (KOŠIR 1975). Rastiščni indeks je tesno povezan s proizvodno sposobnostjo rastišč. Predstavlja zgornjo višino sestoja (po Pardeyu) pri starosti sto let. Raven proizvodnosti nato določimo s pomočjo korelacijskih povezav med zgornjo višino, gostoto in prsnim premerom dreves v sestoju. V praksi služi zgornja višina sestoja kot vhod v donosne tablice, kjer odčitamo vrednost povprečnega volumenskega prirastka v času kulminacije (m'/ha/leto). Rastiščni indeks so žal ugotavljali le v posameznih primerih ter v celoti na eni gozdarski območni enoti (Nazarje). Gre za zahtevno in obsežno delo, ki ga v sedanji organiziranosti gozdarstva opravi lahko le javna gozdarska služba (KOTAR 1996). Deleže lesnih zalog drevesnih vrst za modelno optimalno stanje smo določili na ravni gozdnih združb. Primerjavo lesnih zalog naredimo po razvojnih fazah gozda in sicer samo v primerih, kjer lahko lesno zalogo porazdelimo po razvojnih fazah gozda. Osnovo za razdelitev lesnih zalog po razvojnih fazah naredi model na podlagi Halajevih donosnih tablic (HALAJ 1987). Koeficient spremenjenosti je izračunan le v tistih odsekih, kjer ni prebiralnih in panjevskih gozdov ter listnikov in steljnikov. Skico modela primerjave lesnih zalog po drevesnih vrstah prikazuje slika 1. D =faza, - {K * * delež) * površina. LZ. 100% razlika (D) (model - stanje) 100% — LZ(trenutna) D K delež faza. vsota razlik lesnih zalog med obstoječim in modelno optimalnim stanjem / sum of wood stock dijjerence between existing and optimal situation model končna lesna zaloga po tablicah pri starosti Tmax let (mV ha) /final wood stock according to age tables Tmax years (m'/hn) ® ® rastiščno določen odstotek drevesne vrste v skupni lesni masi sestoja / natural site determined percentage of tree type in total wood mass of stand modelne lesne zaloge v razvojnih fazah / model wood stock in development stage povrsma; povximt^ nzvoin\h {3.7.! development stage areas LZ lesna zaloga (m'/ha) za posamezno drevesno vrsto v odseku / LZ wood stock (m'/ha) by tree type m the tract T starost dreves v razvojni fazi po podatkih popisa gozdov / tree age in development stage according to forest inventory data a s Slika 1: Prikaz izračuna razlik med obstoječim in rastišču primernim stanjem lesnih zalog po drevesnih vrstah Figure 1: Calculation of the difference between existing natural site and its appropriate wood stock situation by tree type Koeficient spremenjenosti {Ks) določimo na ravni odseka tako, da izračunamo kvocient evklidske razdalje razlik modela m dejanskega stanja ter maksimalne razdalje med modelom in dejanskim stanjem (enako, kot po metodi I). + Sf) Ks=100*E/E^ E., i M Ks - koeficient spremenjenosti E = evklidska razdalj a med modelom in dej anskim stanj em = maksimalna evklidska razdalja med modelom in dejanskim stanjem = število vseh različnih drevesnih vrst v odseku (model + dejansko stanje) = modelna višina lesne zaloge drevesne vrste v odseku S = dejanska višina lesne zaloge drevesne vrste v odseku Za obdelavo podatkov na državni ravni smo privzeli za vse gozdove modelno enako proizvodno dobo (120 let), ki pa je dejansko v naravi zelo različna in odvisna od številnih vplivnih dejavnikov (drevesne vrste, proizvodne sposobnosti rastišča, nadmorske višine, gozdne združbe in drugih). 3 REZULTATI RESULTS 3.1 REZULTATI OCENJEVANJA SPREMENJENOSTI GOZDOV PO METODI OCENJEVANJA SPREMENJENOSTI VRSTNE SESTAVE RASTLINSKIH SKUPNOSTI RESULTS OF FOREST COMPSITION ALTERATION ASSESSMENT USING THE METHOD OF ASSESSING COMPSITION ALTERATIONS IN PLANT COMMUNITY TYPE Rezultate obdelave smo predstavili z izračunom povprečnih vrednosti indeksa spremenjenosti (metoda I) po gozdnih združbah in gozdnogospodarskih območjih. Povprečja za gozdne združbe in gozdnogospodarska območja smo izračunali kot tehtane aritmetične sredine, pri čemer smo za utež uporabili površine odsekov. Hkrati smo podatke razvrstili po vrstnem redu od največjega proti najmanjšemu indeksu spremenjenosti vrstne sestave rastlinskih skupnosti. Preglednica 1: Indeks spremenjenosti drevesne vrstne sestave po gozdnogospodarskih območjih m gozdnih združbah na državni ravni Table 1: Index of tree type composition alterations by forest management region and Jorest communities at national level. gozdna zdr. FOREST COM. gozdnogospodarsko območje / FOREST MANAGEMENT REGION Preglednica 1: Nadaljevanje Table 1: Continuation gozdna zdr. gozdnogospodarsko območje / FOREST MANAGEMENT REGION FORESl ^COM. II 10 02 03 12 13 09 04 05 07 14 06 08 Ol ARF 56 85 64 44 31 39 22 19 53 24 13 HQC 55 99 95 72 73 59 76 78 54 66 44 ANF 54 95 55 62 55 38 33 QF 54 40 56 54 49 65 49 FSP 53 77 48 65 RC 53 57 58 50 78 69 74 49 QFL 52 53 78 66 71 80 50 CLA 52 52 OF 51 93 62 61 14 50 27 71 33 53 38 43 HF 50 96 83 69 46 65 60 61 46 45 30 60 SO 47 52 48 31 AF 46 67 60 43 57 45 68 NA 43 69 41 OFO 43 77 49 79 69 16 95 41 33 ADF 43 86 56 52 30 22 20 SAP 42 41 22 18 56 99 76 92 9 4 ASP 41 22 82 23 LÖF 41 59 48 16 PSI 40 29 54 65 13 22 SEF 40 64 69 18 58 ACF 38 29 53 EF 37 86 83 63 76 60 44 32 29 24 31 30 GP 37 32 52 5 38 39 38 2 Al 37 97 23 97 20 EP 35 35 VP 33 54 19 77 45 22 60 41 IF 33 33 MP 32 25 79 13 29 47 38 47 60 ORF 30 29 23 49 CO 28 28 FDF 27 27 LSP 23 27 7 AGP 23 24 23 14 25 CAP 21 21 44 8 10 LNF 21 21 BP 20 21 21 17 20 21 18 10 OP 18 18 FS 16 98 54 7 10 S 15 99 13 36 ApP 11 11 VPI 10 10 13 FSSA 9 9 DP 8 25 7 88 SP 3 3 |hp 0 0 Postopek izdelave slike na ravni Slovenije smo priredili prostorski predstavitvi podatkov Naredili smo rastrsko datoteko Slovenije s 500 m mrežo rastrskih celic, kjer predstavlja vsaka rastrska celica 25 ha površine. Posamezna celica rastra je obarvana v primeru da je centroid (središčna točka) odseka znotraj površine, ki jo predstavlja rastrska celica Na shki 2 smo rastrsko celico obarvah, če je znotraj površine, ki jo predstavlja odsek ki ima določeno samo eno gozdno združbo oz. izračunano spremenjenost drevesne sestave Odseki so v povprečju manjši od 25 ha (povprečna velikost odsekov, ki so vključeni v raziskavo je 10,35 ha). Zato slika prikazuje v povprečju za 38 % večji površinski delež posamezmh vrst spremenjenosti gozdnih sestojev, kakor je v resnici (preglednica 2) Uporabimo JO lahko le kot informacijo o lokaciji v prostoru za posamezne razrede vrstne spremenjenosti sestojev. Preglednica 2: Analiza prostorske predstavitve površin gozdov, za katere smo izračunali spremenjenost vrstne sestave Table 2: Space analysis of representative forest areas for which type composition alterations were calculated INDEKS SPREMENJENOST! [%] ALTERATION INDEX [%] ŠTEVILO RASTR. CELIC NUMBER OF GRID CELLS KUMULATIVA ŠTEVILA RASTRSKIH CELIC CUMULATIVE NUMBER OF GRID CELLS POVRŠINA [ha] AREA [ha] KUMULATIVA POVRŠINE [ha] CUMULATIVE AREA [ha] POVRŠINA PO POPISU GOZDOV [ha] AREA ACCORDING TO FOREST INVENTORY [hal % RAZLIKE V PRIKAZU POVRŠIN % DIFFERENCE IN AREA REVIEW 0-20 3.756 3.756 93.900 93.900 79.907 118% 21-40 4.004 7.760 100.100 194.000 77.605 129% 41-60 3.524 11.284 88.100 282.100 66.454 133% 61-80 4.406 15.690 110.150 392.250 80.086 138% 80-100 9.072 24.762 226.800 61^050 144.874 157% Skupaj / Total 448.926 138% Slika 2: Prostorski prikaz lokacije gozdov po razredih spremenjenosti drevesne sestave (metoda ocenjevanja spremenjenosti vrstne sestave rastlinskih skupnosti) Figure 2: Space review of forest locations by class of tree composition alteration (method of plant community type composition alteration) Če analiziramo frekvenčno porazdelitev števila odsekov po razredih spremenjenosti za osnovne podatke (grafikon 1), ugotovimo sorazmerno večji delež odsekov v najvišji stopnji spremenjenosti (preko 40 %). Rezultati kažejo tudi rahlo povezavo med velikostjo odseka in višino indeksa spremenjenosti, saj imajo majhni odseki sorazmerno večjo spremenjenost drevesne sestave v primerjavi s površinsko večjimi odseki, kar smo testirali tudi s korekcijsko odvisnostjo (Pearsonov koeficient korelacije je -0,132** oz. značilen na 1 % tveganju). 45% 0-20% 21-40% 41-60% 61-80% 81-100% Indeks spremenjenosti //(terarto« index Grafikon 1: Deleži frekvenčne porazdelitve števila odsekov po razredih indeksa spremenjenosti drevesne sestave Graph 1: Share of frequency distribution of tract number by class of tree composition alteration index Prostorska predstavitev na sliki 2 kaže tudi delno spremenjeno podobo razmerij posameznih stopenj indeksa spremenjenosti (preglednica 3). Vzrok je v tehniki izdelave rastrskih datotek z uporabo vektorskih podatkov za točkovne vrednosti. Ista rastrska celica lahko prekriva več odsekov, za katere smo izračunali indeks spremenjenosti drevesne sestave. V takih primerih bo vrednost rastrske celice enaka zadnjemu v zaporedju podatkov za odseke v računalniški zbirki, ki imajo središčno točko znotraj površine, predstavljene z rastrsko celico. Frekvenčna porazdelitev števila odsekov po velikosti površine kaže, da so pogostejši manjši odseki, saj je srednja vrednost za površino analiziranih odsekov 10,35 ha, mediana pa že pri velikosti površine odseka 5,36 ha. Večja je zato verjetnost, daje rastrsko celico določil manjši odsek, ki ima v povprečju višjo stopnjo spremenjenosti. Preglednica 3: Primerjava rastrske predstavitve z izračunanimi podatki Table 3: Comparison of grid representation and calculated data INDEKS SPREMENJENOSTI [%] ALTERA TION INDEX [%] RASTRSKE TOČKE [ŠTEVILO] GRID POINTS fNUMBERl DELEŽ SHARE POVRŠINE PO OSNOVNIH PODATKIH|| IZRAČUNA [%] AREA B Y BASIC CALCULA TION DA TA r%i 0-20 3.756 15,2% 17,8% 21-40 4.004 16,2% 17,3 % 41-60 3.524 14,2 % 14,8% 61-80 4.406 18,8% 17,8% j 80-100 9.072 37,6 % 32,3 % 3.2 REZULTATI ANALIZE PO METODI UGOTAVLJANJA SPREMENJENOSTI KOLIČINSKE IN VRSTNE STRUKTURE LESNIH ZALOG SESTOJEV ANALYSIS RESULTS ACCORDING TO THE METHOD OF ESTABLISHING ALTERATIONS IN QUANTITY AND TYPE OF STAND WOOD STOCK STRUCTURE Tudi pri izračunu koeficienta spremenjenosti (po metodi II) smo v analizo vključili le odseke s samo eno gozdno združbo. Manjši delež analiziranih gozdov je posledica metode izračuna, uporabimo pa jo lahko le v odsekih z razvojnimi fazami mladovja, drogovnjakov (mlajšega in starejšega), debeljaka ter pomlajenca. Koeficient spremenjenosti lahko zavzema vrednosti na intervalu od nič do sto. Tudi za metodo ugotavljanja spremenjenosti sestojev po količinski in vrstni strukturi lesnih zalog smo izdelali prostorsko predstavitev (slika 3), po enaki metodi kot pri metodi ocenjevanja spremenjenosti vrstne sestave rastlinskih skupnosti (slika 3), zato ima tudi enake lastnosti (25 ha površina rastrske celice in postopek rastrske prestavitve vektorskih točk za centroide odsekov). Slika 3: Prostorski prikaz lokacije gozdov po razredih spremenjenosti drevesne sestave (metoda spremenjenosti količinske in vrstne strukture lesnih zalog sestojev) Figure 3: Space review of forest locations by class of tree composition alteration (method of plant community type composition alteration) Kot osnovo za razvrstitev gozdnih združb in gozdnogospodarskih območij po vrstnem redu spremenjenosti sestojev smo uporabili tehtane aritmetične sredine vrednosti koeficienta spremenjenosti v odsekih, kjer je bila utež površina odseka. Preglednica 4: Povprečne vrednosti koeficienta spremenjenosti sestojev po gozdnogospodarskih območjih in gozdnih združbah Table 4: Mean values of stand alteration coefficient by forest management region and forest community FOREST ( ZOMMV' 14 02 11 10 03 04 Ol 09 12 05 08 07 06 13 Povprečje Mean 37,7 62,0 53,1 52,3 44,8 41,1 38,0 ' 35,2 : 33,9 ' 33,8 ! 31,5 29,2 28,S 1 28,3 23,4 RR 93,5 93 95 SO 79,4 79 88 70 PM 77,9 61 97 76 AG2 73,3 73 AGl 64,8 65 OC 60,8 63 57 LUQ 59,9 36 37 68 LAQ 58,4 47 44 62 COF 58,1 58 Al 56,1 61 46 FSS 55,8 56 TA 55,5 44 58 LXF 54,6 55 QLF 54,3 33 66 53 47 SF 54,1 61 37 64 47 36 40 FSP 53,9 57 34 SQ 51,4 69 37 LFP 51,2 33 66 UA 50,7 47 47 54 AFP 49,6 55 52 50 51 40 32 77 AFR 48,7 50 23 30 27 20 48 37 58 LA 48,3 49 39 46 DA 48,2 47 42 56 57 61 49 49 35 47 VPI 46,8 42 47 CQ 46,7 47 LYA 46,4 46 CF 46,0 58 60 36 23 38 17 16 BF 45,8 56 54 36 41 42 44 40 29 28 33 21 LF 45,4 57 57 53 57, 29 34 33 42 27 22 26 CVF 44,7 51 27 BA 44,5 45 58 56 38 34 33 43 ANF 43,6 59 67 54 50 28 AGl 42,6 1 15 22 77 Okrajšave imen gozdnih združb povzete po šifrantu popisa gozdov iz 1. 1990 (priloga 1) / Codes of abbreviated forest community names according to the forest inventory of 1990 (apendix 1) Imena gozdnogospodarskih območij so v preglednici 6 / Codes for forest management regions are in table 6 Preglednica 4: Nadaljevanje Table 4: Continuation gozdna združba FOREST COMMU. Preglednica 4: Nadaljevanje Table 4: Continuation GOZDNA ZDRUŽBA GOZDNOGOSPODARSKO OB VlOČJE / FORESTM4N4nFMPA,r ppnrn^r FOREST C :OMMU. 14 02 11 10 03 04 01 09 12 05 08 07 06 -iTj LAP 21,5 22 HF 19,6 26 20 16 20 13 17 RC 19,3 64 7,4 23 -- 6 LOF 19,2 28 18 CAP 19,0 22 10 16 ADF 17,7 18 FDF 17,4 17 LQC 17,3 16 31 25 27 9 CRF 16,8 24 3 EF 16,5 11 20 38 13 HP 15,7 16 SAP 15,5 8 56 9 9 18 DP 14,5 14 71 SF 14,4 15 14 ANF 14,2 15 14 SAP 14,1 11 20 AG3 11,2 87 20 7 AGI 9,0 5 11 18 MQ 8,3 8 CAP 6,1 6 S 5,7 13 5 5 Po metodi II je ocena spremenjenosti slovenskih gozdov bolj ugodna, saj ima največji delež analiziranih odsekov tretji razred spremenjenosti višine in vrstne strukture lesnih zalog (preglednica 5), v najvišjem razredu pa je manj kot 10 % odsekov, Korelacija med velikostjo odseka in vrednostjo razreda spremenjenosti je tudi pri tej metodi negativna, vendar neizrazita. Preglednica 5: Frekvenčna porazdelitev odsekov po razredih koeficienta spremenjenosti rable 5: Frequency distribution of tracts by class alteration coefficient koeficient spremenjenosti ALTERA TION COEFFICIENT [%] 3.3 PRIMERJAVA REZULTATOV OBEH METOD OCENJE V ANJA SPREMENJENOSTI GOZDOV COMPARISON OF RESULTS FROM BOTH METHODS OF FOREST COMPOSITION ALTERATION ASSESSMENT Metodi se razlikujeta predvsem v številu vplivmh dejavnikov, ki določajo stopnjo spremenjenosti gozda, kar ima za posledico tudi različen obseg vhodnih podatkov Izračun indeksa spremenjenosti {Is, metoda I) je sorazmerno preprost m upošteva zgolj razlike med deleži lesnih zalog dejansko prisotnih in modelno optimalne strukture drevesmh vrst na določenem rastišču. Koeficient spremenjenosti {Ks, metoda II) je bolj kompleksen, saj poleg ustreznosti vrstne strukture upošteva tudi absolutne višine lesmh zalog v gozdnem sestoju. Degradirane gozdove z nizko lesno zalogo in ustrezno drevesno sestavo prišteva med spremenjene, medtem ko z v takih gozdovih ne opozorimo na neustreznost sestojnih struktur. Prednost je v široki uporabi zaradi majhnega obsega potrebnih podatkov. Imeti moramo le numenčno informacijo o dejanskem m modelno optimalnem deležu v lesni zalogi, temeljmci, zastiranju, pokrovnosti itd. za posamezne rastlinske skupnosti. Uporabimo ga lahko tudi za ocenjevanje spremenjenosti prebiralnih panjevskih gozdov ter listnikov m steljnikov, čeprav je določitev optimalnih struktur v nekaterih izmed njih zelo zahtevno delo. Ne moremo pa med seboj neposredno primerjati vrednosti ocen spremenjenosti po obeh metodah, ker upoštevata pri izračunu različna izhodišča oz. različne vplivne dejavnike, ki določajo stopnjo spremenjenosti gozdnih sestojev. Zato smo uporabili metodo rangiranja, ki za posamezno metodo kaže vrstni red stopenj spremenjenosti na ravni gozdnogospodarskih območij. Preglednica 6: Rangi gozdnogospodarskih območij po obeh metodah ocenjevanja spremenjenosti gozdnih sestojev Table 6: Ranking of forest management regions by both methods of forest stand composition alteration assessment ŠIFRA CODE GOZDNOGOSPODARSKO OBMOČJE FOREST MANA CEMENT RECION METODA I / METHOD I METODA II / METHOD II II Rang {Is) Rank (Is) Število analiziranih primerov odsekov Number of analysed tract samples Rang (Ks) Rank (Ks ) Število analiziranih primerov odsekov Number of analysed tract samples 01 TOLMIN 14 1.409 7 740 02 BLED 3 824 2 187 03 KRANJ 4 1.509 5 1.155 04 UUBLJANA 8 2.297 6 1.955 05 POSTOJNA 9 1.892 10 561 06 KOČEVJE 12 4.597 13 2.977 07 NOVO MESTO 10 3.342 12 2.203 08 BREŽICE 13 2.398 11 1.633 09 CEUE 7 5.108 8 4.892 10 NAZARJE 2 2.390 4 2.110 11 SLOVENJ GRADEC 1 7.710 3 6.622 12 MARIBOR 5 3.480 9 3.345 13 MURSKA SOBOTA 6 2.959 14 2.733 14 KRAŠKO OBMOČJE 11 3.555 1 821 Velika razlika med vrednostmi rangov je predvsem v primerih, kjer smo zaradi potrebnih vhodnih podatkov ter metode izračuna analizirali zelo različno število primerov spremenjenosti gozda po odsekih. Rang se najbolj razlikuje na kraškem območju, kjer je tudi največ gozdov, za katere metoda II ne predvideva postopka ocene spremenjenosti sestojev (panjevci, listniki, steljniki, grmišča). Ocenjujemo, da ima za določeno, predvsem raziskovalno rabo kompleksnejša metoda računanja Ks veliko uporabno vrednost. Izračun Ks je zahtevnejši z večjim obsegom vhodnih podatkov, saj moramo poleg struktur lesnih zalog po drevesnih vrstah poznati tudi dejanske in modelno optimalne višine lesnih zalog pri določeni starosti sestojev oz. po razvojnih fazah gozda ter proizvodno dobo. V analizi smo uporabili zelo grob približek za proizvodno dobo ter prirejene vrednosti donosnih tablic za optimalne velikosti lesnih zalog. Primerjava višin lesnih zalog po drevesnih vrstah kaže lahko tudi na izkoriščenost proizvodne sposobnosti rastišč. Izračun se lahko bolje prilagodi posameznemu primeru oz. določitvi rastišču primernih višin in vrstnih struktur lesnih zalog. Za obe metodi je primerna računalniška obdelava podatkov iz popisa gozdov na različno velikih ravneh velikosti ureditvenih površin. Stopnja spremenjenosti struktur gozdnih sestojev je kazalec nujnosti m prioritete ukrepanja oz. načrtovanja del v gozdu, zato bi jo lahko uporabili kot koristno informacijo pri strateškem, taktičnem in operativnem gozdnogospodarskem načrtovanju. 4 POVZETEK Raziskava primerja rezultate izračuna spremenjenosti višine in vrstne strukture lesnih zalog v delu slovenskih gozdov, za katerega smo pridobili za izračun potrebne podatke. Izbiro gozdnih sestojev (odsekov) so narekovali podatki, ki so omogočih objektivno m čim bolj natančno pridobivanje in vrednotenje rezultatov dveh metod ocenjevanja spremenjenosti gozda. Prva je metoda ocenjevanja spremenjenosti vrstne sestave rastlmskih skupnosti (BONČINA / ROBIČ 1998), druga pa metoda ugotavljanja spremenjenosti sestojev po višini in vrsti lesnih zalog (KRČ 1999). Obe metodi izračunata kvantitativno oceno spremenjenosti, pri čemer upoštevata sedanje in potencialno naravno oz. optimalno stanje strukture lesnih zalog po gozdnih združbah. Metodi zahtevata podatke o višini in strukturi lesnih zalog po drevesnih vrstah na ravni gozdnih združb. Prvi metodi zadoščajo samo deleži v lesni zalogi, temeljnici ali pokrovnosti po drevesnih vrstah, medtem ko druga metoda analizira absolutne višine lesnih zalog po drevesnih vrstah in razvojnih fazah gozda. Druga metoda je podatkovno bolj zahtevna in postopek izračuna spremenjenosti je obsežnejši. Določili smo nekaj približkov, poenostavitev in modelnih postopkov za zagotovitev potrebnih vhodnih podatkov (proizvodna doba, prirejene vrednosti donosnih tablic za višino optimalnih lesnih zalog, modelna razporeditev lesnih zalog po razvojnih fazah gozda ter optimalna struktura lesnih zalog po drevesnih vrstah v gozdnih združbah), ki so omogočili izračun na državni ravni. V analizo smo vključili le odseke, ki imajo samo eno gozdno zdmžbo. so gozdna zemljišča in imamo zanje podatke o lesni zalogi po drevesnih vrstah. Odseki, ki so razdeljeni na več gozdnih združb (največ tri), ne dajo potrebnih podatkov za uporabo metode ocenjevanja spremenjenosti vrstne sestave rastlinskih skupnosti, saj nimamo podatka o lesnih zalogah po drevesnih vrstah in gozdnih združbah, pač pa le o lesnih zalogah na ravni odsekov. Za izračun po obeh metodah smo izdelali tudi model prostorskega prikaza stopenj spremenjenosti gozda na rastrski karti s 25 ha površino rastra. Model kaže približno sliko prostorske razporejenosti gozdov glede na izračunano oceno stopnje spremenjenosti gozdov. Izdelali smo jo z rastrsko prestavitvijo središčnih točk analiziranih odsekov. Za obe metodi smo izračunali tudi spremenjenost gozdnih sestojev po gozdnih združbah in gozdnogospodarskih območjih. Razmerja med stopnjami spremenjenosti na ravni gozdnogospodarskih območij, ki smo jih predstavili v obliki rangov, se v obeh primerih bistveno ne razlikujejo. Razlike med metodama so največje v območjih, kjer je povprečje spremenjenosti gozdnih sestojev na ravni gozdnogospodarskega območja pridobljeno z različnim obsegom osnovnih izračunov po odsekih. Veliko večja razlika je v frekvenčni porazdelitvi stopenj spremenjenosti. Strožjo oceno smo dobili z uporabo prve metode, kjer smo večino analiziranih odsekov ocenili z najvišjim razredom stopnje spremenjenosti gozdnih sestojev. Po metodi II je ocena spremenjenosti slovenskih gozdov bolj ugodna, saj ima največji delež analiziranih odsekov fretji razred spremenjenosti višine in vrstne strukture lesnih zalog (preglednica 5), v najvišjem razredu pa je le 5 % odsekov. Uporaba metod ocenjevanja spremenjenosti je odvisna od potrebe in podatkov, ki jih imamo. Obe metodi sta primerni za uporabo računalniške obdelave podatkov o gozdu, kar močno olajša analizo na različnih ravneh velikosti gozdnih površin. 5 SUMMARY The study compares results of calculated composition alterations of amount and type of wood stock structure in a part of Slovenian forests for which we obtained necessary data for calculation. The selection of forest stands (tracts) dictated the data enabling an objective and as accurate as possible acquisition of result values by two methods of forest fluctuation assessment. The first is the method of estimation of the species composition alteration in plant communities (BONČINA / ROBIČ 1998), the second, a method of establishing stand composition alteration by amount and type of wood stock (KRČ 1999). Both methods calculate the cumulative value of composition alteration in which the present and potential nature or optimal situation of wood stock structure by forest community is taken into account Both methods require data on the amount and structure of wood stock by tree type at the level of forest communities. The first method satisfies only the shares in wood stock and basal area by tree type, while the second method analyses absolute amounts of wood stock by tree type and forest development stages. The second method is more demanding as regards data and the procedure of calculating composition alterations is more comprehensive. We took certain short-cuts, simplifications and model procedures for establishing the required input data (production period increase values from profit tables for the amount of optimal wood stock model distribution of wood stock by forest development stage and the optimal wood stock structure by tree type in forest communities) that enabled calculation at national level In the analysis we included only tracts having only one forest community, are forest land and have data on wood stock by tree type. Tracts divided into more than one forest community (at most three) do not give the required data for using the method of assessing plant community type composition alteration as they do not have data on wood stock by type and forest community but only on wood stock at the level of tracts. For calculating by both methods we made a model of space review of the levels of composition alteration of forest on a card with 25 hectare basic raster cells. The model shows an approximate picture of space distribution of forests with respect to the calculated value of level of forest composition alteration. We made a raster presentation of mean points of analysed tracts. For both methods, we also calculated forest stand composition alteration by forest community and forest management region. The relationship between level offiuctuation at the level of forest management region is presented in the form of ranks that do not significantly differ in either case. Differences between methods are most pronounced in regions where the mean forest stand composition alteration at the level of forest management region is obtained with different extents of basic calculations by tract. There is a greater difference in the frequency of distributed fluctuation levels. The first method provided a more precise assessment where we assessed most analysed tracts with the highest class of forest stand composition alteration level. According to method II the assessment of Slovenia forest composition alteration was more favourable as most shares of analysed tracts have third class composition alteration amounts and types of wood stock structure (table 5) and only 5% of tracts were in the highest class. The use of method for assessing composition alterations is dependant on needs and data available. Computer processed data on forests are suitable for both methods as analysis is made easier at different levels of forest area size. 6 VIRI REFERENCES BONČINA, A. / ROBIČ, D., 1998, Ocenjevanje vrstne sestave rastlinskih skupnosti,-Zbornik gozdarstva in lesarstva, 57, s. 113-130. ČOKL, M. 1961. Gozdarski in lesnoindustrijski priročnik.- Ljublajna, Inštitut za gozdno in lesno gospodarstvo, Strokovna in znanstvena dela. DZS, 366 s. HALAJ, J. / GREK, J. / PANEK,F. / PETRAS, R. / 0EHÄK, J., 1987. Rastove tabul'ky hlavnyh dre vin ES SR.- Priroda, Bratislava, 362 s. KOŠIR, Ž., 1975. Zasnova uporabe prostora. Gozdarstvo.- Ljubljana, Zavod SRS za družbeno planiranje, 137 s. KOTAR, M., 1996. Poznavanje lesnoproizvodne sposobnosti gozdnih rastišč kot pogoj za kakovostne odločitve pri ravnanju z gozdovi.- Zbornik gozdarstva in lesarstva, 50, s. 221-231. KRČ, J., 1999. Večkriterialno dinamično vrednotenje tehnoloških, ekonomskih, socialnih in ekoloških vplivov na gospodarjenje z gozdovi,- Doktorska disertacija. Ljubljana, BF, Odd. za gozdarstvo in obnovljive gozdne vire, 173 s. SCOTT, J. M. / JENNINGS M. D., 1997. A Description of the National Gap Analysis Program.URL:http://www.gap.uidaho.edu/GAP/AboutGAP/GapDescription/Index.ht 7 PRILOGA APPENDIX Priloga 1: Šifrant okrajšav imen gozdnih združb po popisu gozdov iz leta 1990 Appendix 1: Code of abbreviated forest community names according to the forest inventory of 1990 gozdna združba / FOl naziv TITLE lEST COMMUNITY okrajšava ABBREVIATION gozdna združba / FORES', naziv TITLE r COMMUNITY okrajšava ABBREVIATIOI querco robori -carpinetum rc festuco drymeiae -fagetum fdf querco robori -ulmetum qu polygonatum verticillati-luzulo- fagetum lfp carici elatae -alnetum glut. ag i querco-luzulo-fagetum qlf carici elongatae-alnetum glut. ag2 blechno-fagetum bf carici brizoidi -alnetum glut. ag3 deschampsio - fagetum df alnetum glutinoso-incanae agi abieti - fagetum dinaricum af alnetum incanae ai abieti - fagetum prealp. din. afp salici - populetum sap abieti-fagetum praealpinum paf salicetum gr. s neckero - abietetum na querco - carpinetum v.hacq. hqc asplenio- abietetum asa querco - carpinetum v.luzula lqc festuco - abietetum fa ornithogalo pyrenaici-carpinetum oc clematido-abietetum cla latyro - quercetum laq lycopodio - abietetum lya orno - quercetum petr. -pub onq luzulo - abietetum la carici umbrosae-quercetum petraeae cq dryopterido - abietetum da seslerio autumnalis-quercetum petraeae sq oxalido - abietetum oxa luzulo - quercetum luq bazzanio- abietetum ba melampyro vulgati -quercetum mq asplenio - piceetum asp fagetum submontanum var.geogr. sesleria 1 autumnalis fssa carici albae-piceetum cap Priloga 1: Nadaljevanje Apendix 1: Continuation GOZDNA ZDRUŽBA / FOREST COMMUNITY GOZDNA ZDRUŽBA / FOREST COMMUNITY NAZIV / TITLE OKRAJŠAVA ABBREVIATION NAZIV / TITLE OKRAJŠAVA ABBREVIATION SESLERIO - FAGETUM SEF APOSERI - PICEETUM ApP HACQUETIO - FAGETUM HF ADENOSTYLO GLABRAE-PIC. -till_ AGP FAGETUM SUBMONTANUM SUBMEDITERRANEUM FSS ADENOSTYLO ALLIARIE-PIC. AAP FAGETUM SUBMONTANUM PRAEALPINUM FSP CAL.VILLOSAE- PICEETUM VPI ENNEAPHYLLO-FAGETUM EF LUZULO ALBIDAE-PICEETUM LAP ORVALO-FAGETUM ORF PICEETUM SUBALPINUM DINARICUM PSD ANEMONE - FAGETUM ANF PICEETUM MONTANUM PIM LAMIO ORVALAE-FAGETUM PRAEALPINUM LOF SORBO - PICEETUM SOP SAVENSI - FAGETUM SF BAZZANIO- PICEETUM BP ADENOSTYLO-FAGETUM ADF SPHAGNO - PICEETUM SP LARICI - FAGETUM LXF HOMOGYNO- PICEETUM HP LUZULO NIVEAE-FAGETUM LNF LUZULO SYLVATICAE-PICEETUM LSP CORYDALO OCHROLEUCAE-FAGETUM COF DESCHAMSIO-PICEETUM DP FAGETUM SUBALPINUM FS GENISTO - PINETUM GP OSTRYO-FAGETUM OF PINETUM SUBILLYRICUM PSI CARICI ALBAE - FAGETUM CF ORNO - PINETUM OP CALAMAGROSTIDO VARIAE - FAGETUM CVF ERJCO- PINETUM EP ARUNCO-FAGETUM ARF VACCINIO-VITIS IDEAE-PINETUM VP ISOPRYO - FAGETUM IF MYRTILLO-PINETUM MP ACERI - FAGETUM ACF TILIO- ACERETUM TA QUERCO-FAGETUM QF ULMO - ACERETUM UA QUERCO-FAGETUM VAR.LUZULA QFL ACERI - FRAX1NETUM AFR LUZULO - FAGETUM LF CARICI REMOTAE - FRAX. CRF SESLERIO - OSTRYETUM SO QUERCO - OSTRYETUM QO RHODODEND,- RHODOTH. RR OSTRYO - FRAXINETUM ORNII OFO PINETUM MUGHI PM CYTISANTO - OSTRYETUM CYO TILIO - OSTRYETUM TO