Kratki znanstveni prispevek

Napovedi o zdravju gozdov, 2025

DOI: 10.20315/NZG.74

Preverjanje verjetnostne napovedi sanitarnega poseka smreke zaradi podlubnikov v Sloveniji v 2025

Nikica OGRIS¹*, Maarten de GROOT¹

¹Gozdarski inštitut Slovenije, Večna pot 2, 1000 Ljubljana

*nikica.ogris@gozdis.si

Datum izdaje: 02.12.2025
Veljavnost: 2025

Ključne besede: navadna smreka, Picea abies, sanitarni posek, napoved, ogroženost, model, validacija, zmogljivost, zanesljivost, točnost, natančnost, AUC, občutljivost, specifičnost

Izvleček

Preverili smo zanesljivost verjetnostne napovedi sanitarnega poseka smreke zaradi podlubnikov v Sloveniji v 2025. Verjetnostni model za napoved sanitarnega poseka smreke zaradi podlubnikov je potrdil visoko zanesljivost (AUC modela = 0,83, AUC napovedi = 0,82). Napoved za leto 2025 smo naredili s pragom 0,50, ki pa se je izkazal za prenizkega, saj je bilo 23,4 % modelskih celic lažno pozitivnih in 7,1 % lažno negativnih. Optimalen prag za verjetnostni model v letu 2025 je bil 0,55. Optimalni prag med leti niha in ga ni mogoče določiti v naprej. Optimalni prag v obdobju 2020–2025 je imel vrednosti 0,30–0,55 s povprečjem 0,47. Že pri validaciji verjetnostne napovedi sanitarnega poseka smreke zaradi podlubnikov za leto 2024 smo ugotovili, da optimalni prag gravitira proti 0,5 in smo ga že takrat predlagali za izdelavo vseh nadaljnjih verjetnostnih napovedih sanitarnega poseka smreke zaradi podlubnikov v Sloveniji.

Validation of the probabilistic forecast of sanitation felling of spruce due to bark beetles in Slovenia in 2025

Abstract

We tested the reliability of the probabilistic forecast of the sanitation felling of spruce due to bark beetles in Slovenia in 2025. The probabilistic model for the prediction of the sanitation felling of spruce due to bark beetles showed high reliability (model AUC = 0.83, forecast AUC = 0.82). A threshold of 0.50 was used for the forecast of 2025, but this proved to be too low, as 23.4% of the model cells were false positives and 7.1 % false negatives. The optimal threshold for the probabilistic model in 2025 was 0.55. The optimal threshold in the period 2020–2025 ranged from 0.30 to 0.55, with an average of 0.47. Already during the validation of the probabilistic forecast of sanitary felling of spruce due to bark beetles for the year 2024, we found that the optimal threshold gravitates toward 0.50, and we proposed it at that time for the preparation of all subsequent probabilistic forecasts of sanitary felling of spruce due to bark beetles in Slovenia.

Uvod

Gozdarski inštitut Slovenije na podlagi modela (de Groot in Ogris, 2019) vsako leto izda verjetnostno napoved sanitarnega poseka smreke zaradi podlubnikov za tekoče leto. Validacija modelov je pokazala, da je zanesljivost napovedi relativno visoka, saj je vrednost AUC (angl. Area under the curve) za model znašala 0,83. Ko izdelamo napoved za novo obdobje, uporabimo obstoječi model za ekstrapolacijo. Zanesljivost ekstrapolirane napovedi pa se lahko razlikuje od zanesljivosti, ki je bila izračunana v procesu razvoja modela. Posledično so lahko napovedi iz leta v leto različno zanesljive. Cilj raziskave je bil preveriti zanesljivost verjetnostne napovedi sanitarnega poseka smreke zaradi podlubnikov v Sloveniji v 2025 (Ogris in de Groot, 2025).

Metode dela

Preverjali smo zanesljivost verjetnostne napovedi sanitarnega poseka smreke zaradi podlubnikov v Sloveniji v letu 2025 (Ogris in de Groot, 2025). Vrednosti verjetnostne napovedi se gibljejo od 0 do 1, kjer 1 pomeni največjo verjetnost, 0 pa najmanjšo. Napoved je bila narejena samo za območja, kjer se pojavlja smreka. Prostorska ločljivost napovedi je 1 km × 1 km. Za vsako modelsko celico smo preverili, ali je bila v njej zabeležena označba dreves za posek smreke zaradi podlubnikov, kar je vključevalo sanitarni posek lubadark (vrsta poseka 301), posek oslabelih dreves (vrsta poseka 901) in dodatno posekana drevesa (vrsta poseka 991). Upoštevali smo vse podatke iz zbirke o evidenci označbe dreves za posek iz 2025 (zbirka Timber), ki so bili zbrani do 1. 12. 2025 (ZGS, 2025).

Zanesljivost napovedi smo izračunali s pomočjo več kazalcev (Witten in Frank, 2005):

Točnost = (RP + RN)/(RP + RN + LP + LN),
Natančnost = RP/(RP + LP),
Občutljivost = RP/(RP + LN),
Specifičnost = RN/(RN + LP),
Koeficient podobnosti = 2×RP/(2×RP + LP + LN),
AUC = površina pod krivuljo Občutljivost (1 – Specifičnost),

kjer RP = resnično pozitiven, RN = resnično negativen, LP = lažno pozitiven, LN = lažno negativen.

Ker je napoved verjetnostna [0, 1], smo kazalce zanesljivosti napovedi izračunali večkrat (21- krat) in sicer na intervalu [0, 1] s korakom 0,05. S pomočjo kazalcev zanesljivosti napovedi smo določili kakšen je bil optimalni prag za verjetnostno napoved sanitarnega poseka smreke zaradi podlubnikov v 2025. Optimalni prag za napoved smo izbrali, tako da smo najprej za vsak kazalec zanesljivosti napovedi določili optimalni prag, nato pa smo izbrali optimalni prag za napoved, kjer se je večina krivulj kazalcev med seboj sekala.

Verjetnostno napoved za leto 2025 smo naredili pri pragu 0,50 (Ogris in de Groot, 2025). Zato smo preverili, kakšna je bila zanesljivost napovedi pri pragu 0,50.

Rezultati

Vrednost kazalca AUC je znašala 0,82. Točnost je naraščala do praga 0,50, kjer je dosegla vrednost 0,79, potem pa je upadla (slika 1). Natančnost je naraščala do verjetnosti 1, kjer je znašala 0,99. Občutljivost je bila največja pri pragu 0, kjer je imela vrednost 0,96, potem pa je do verjetnosti 1 postopoma padla na 0. Krivulja za specifičnost je imela obliko črke S. Specifičnost je naraščala od verjetnost 0 do 1. Koeficient podobnosti je imel pri pragu 0 vrednost 0,80, potem je naraščal do praga 0,45, pri katerem je dosegel vrednost 0,85. Nato je do verjetnosti 1 padel na vrednost 0. Optimalni prag za napoved v letu 2025 je bila verjetnost 0,55, kjer so se sekale tri krivulje, tj. natančnost, občutljivost in koeficient podobnosti.

Pri pragu 0,50 je bilo lažno pozitivnih modelskih celic 23,4 % in so se večinoma nahajale v neposredni bližini resnično pozitivnih modelskih celic (slika 2). Lažno negativnih modelskih celic pri pragu verjetnosti 0,50 je bilo 7,1 %, resnično negativnih celic 36,7 % in resnično pozitivnih 32,8 %.

Slika 1: Točnost, natančnost, občutljivost, specifičnost in koeficient podobnosti verjetnostne napovedi sanitarne sečnje smreke zaradi podlubnikov v 2025 glede na različen prag napovedi

Slika 2: Preverjanje verjetnostne napovedi sanitarne sečnje smreke zaradi podlubnikov v 2025 z upoštevanjem praga 0,50 (verjetnost napovedi)

Razprava

Z raziskavo smo ugotovili visoko zanesljivost verjetnostne napovedi sanitarne sečnje smreke zaradi podlubnikov v 2025, saj je vrednost AUC znašala 0,82, kar je le za malenkost manj kot v izvirni validaciji modela (de Groot in Ogris, 2019). Ugotovili smo, da je bil optimalni prag pri verjetnosti 0,55, tj. v bližini mesta na grafikonu (slika 1), ko se je sekala večina krivulj kazalcev zanesljivosti napovedi. Optimalni prag v obdobju 2020–2025 je imel vrednosti 0,30–0,55 s povprečjem 0,47. Že pri validaciji verjetnostne napovedi sanitarnega poseka smreke zaradi podlubnikov za leto 2024 smo ugotovili, da optimalni prag gravitira proti 0,50, kar je statistično pričakovano, in smo ga že takrat predlagali za izdelavo vseh nadaljnjih verjetnostnih napovedih sanitarnega poseka smreke zaradi podlubnikov v Sloveniji (Ogris in de Groot, 2024).

Obravnavan model ima še veliko možnosti za izboljšave, s katerimi bi lahko dodatno povečali zanesljivost napovedi sanitarne sečnje smreke zaradi podlubnikov. Upoštevati je treba dejstvo, da v času preverjanja napovedi v podatkovni zbirki Timber še niso bili vneseni vsi podatki za leto 2025 in zato morda nekaterih lokacij, kjer je prišlo do pojava žarišč smrekovih podlubnikov, nismo upoštevali. V letih 2020–2024 smo v validacijo vključili 2,0–12,8 % (povprečno 7,1 %) manj evidentiranega volumna poseka smreke zaradi podlubnikov (ZGS, 2025).

Znano je, da napovedna zanesljivost modelov s časom upada. Zato se morajo vsakih nekaj let modeli ponovno kalibrirati z novimi podatki. Predlagamo, da se model za kratkoročno napoved sanitarnega poseka smreke zaradi podlubnikov v Sloveniji ponovno kalibrira.

Zahvala

Članek je nastal v okviru Javne gozdarske službe, naloge 2 (Poročevalska, prognostično-diagnostična služba za gozdove) na Gozdarskem inštitutu Slovenije. Recenzentu se zahvaljujemo za koristne pripombe in predloge izboljšav.

Viri

e Groot M., Ogris N. 2019. Short-term forecasting of bark beetle outbreaks on two economically important conifer tree species. Forest Ecology and Management, 450: 117495. https://doi.org/10.1016/j.foreco.2019.117495

Ogris N., de Groot M. 2024. Preverjanje verjetnostne napovedi sanitarnega poseka smreke zaradi podlubnikov v Sloveniji v 2024. Napovedi o zdravju gozdov, 2024: 4. http://dx.doi.org/10.20315/NZG.70

Ogris N., de Groot M. 2025. Verjetnostna napoved sanitarnega poseka smreke zaradi podlubnikov v Sloveniji v 2025. Napovedi o zdravju gozdov, 2025. http://dx.doi.org/10.20315/NZG.73

Witten I.H., Frank E. 2005. Data mining: practical machine learning tools and techniques, 2^nd edition. San Francisco, Morgan Kaufmann: 524 str.

ZGS. 2025. Timber. Podatkovna zbirka o poseku gozdnega drevja. Zavod za gozdove Slovenije

Citiranje: Nikica OGRIS, Maarten de GROOT. 2025. Preverjanje verjetnostne napovedi sanitarnega poseka smreke zaradi podlubnikov v Sloveniji v 2025. Napovedi o zdravju gozdov, 2025. URL: https://www.zdravgozd.si/prognoze_zapis.aspx?idpor=74. DOI: 10.20315/NZG.74

Prispelo: 01. 12. 2025. Sprejeto: 02. 12. 2025. Objavljeno: 02. 12. 2025.

Št. ogledov: 30

na vrh strani