Introduzione 

La conoscenza puntuale dello stato delle risorse naturali è requisito decisivo per un’efficace gestione forestale. All’uopo, è essenziale disporre di banche dati ampie ed aggiornate e di strumenti in grado di rendere le informazioni disponibili facilmente accessibili ed interrogabili. In questa prospettiva, in molte discipline, tra cui la pianificazione e la gestione del patrimonio forestale, si sono affermati i Sistemi di Supporto alle Decisioni (DSS), che rappresentano uno degli strumenti fondamentali della Precision Forestry ([5]).

Una definizione di DSS ampiamente accettata nell’ambito delle scienze informatiche è proposta da Burstein & Holsapple ([2]), che ne individuano le seguenti componenti fondamentali: il linguaggio di sviluppo, il sistema di presentazione dei risultati, il sistema di gestione della conoscenza, il sistema di analisi dei problemi. La caratteristica principale dei DSS è di poter rendere esecutivi strumenti di modellistica previsionale con i quali supportare le attività di gestione e pianificazione. Una rassegna analitica in merito ai principali DSS in uso nel settore forestale è reperibile in Gordon ([7]) e Borges et al. ([1]).

L’integrazione fra DSS e moderne tecnologie in ambito geografico, come ad esempio i Sistemi di Posizionamento Globale Satellitare (GNSS), i dati telerilevati (acquisiti da satellite, da aereo o da drone) e i Sistemi Informativi Geografici (GIS), ha portato negli ultimi anni allo sviluppo dei cosiddetti Sistemi Spaziali di Supporto alle Decisioni (Spatial Decision Support System - SDSS). Alcuni esempi di SDSS sviluppati in Italia per il settore forestale, e in particolare nello specifico ambito della determinazione delle risorse legnose e della loro ritraibilità in termini di legname da opera e di biomasse per la produzione di energia, sono Biomasfor ([12]), PoodyBLEP ([6]), FopP ([8]).

Un quadro dei DSS forestali sviluppati in Italia è fornito da Portoghesi et al. ([11]), che evidenziano la problematica della “manutenzione” nel tempo dei DSS, spesso confinati entro i limiti temporali del progetto di ricerca nel cui ambito sono stati sviluppati, e la difficoltà per il loro trasferimento nella pratica gestionale. Con l’intento di ovviare a questi limiti, nel presente lavoro vengono illustrati i metodi e le caratteristiche architetturali di un nuovo SDSS chiamato CFOR (Calabrian FORests), messo a punto nell’ambito del progetto PON “Ambi.Tec.Fil.Legno” (ALForLab). Attraverso una serie di strumenti operanti sul GIS open source QGIS (Fig. 1), l’obiettivo principale è fornire al gestore forestale un utile e pratico supporto per migliorare la fornitura di materiale legnoso nella Regione Calabria.

Fig. 1 - L’interfaccia grafica di CFOR su QGIS. In rosso l’extent della Area di interesse selezionata dall’utente.

 Enlarge/Shrink  Open in Viewer

Architettura del sistema 

Il plugin di QGIS CFOR si presenta come una suite costituita da otto moduli, ciascuno dei quali deputato all’esecuzione di una specifica fase dell’elaborazione (Fig. 2):

Fig. 2 - Schema dell’architettura del modello CFOR. All’interno di ogni box, ciascuno relativo a uno dei diversi moduli implementati, i rettangoli rappresentano gli algoritmi principali; gli ovali rappresentano le variabili di input e di output degli algoritmi; gli esagoni individuano gruppi di algoritmi già definiti in altri moduli. In basso a destra il logo di CFOR.

 Enlarge/Shrink  Open in Viewer

• area di interesse (delineazione da parte dell’utente dell’area geografica di lavoro all’interno della quale verranno eseguite le successive elaborazioni);
• esposizione (calcolo dell’esposizione, in gradi sessagesimali);
• pendenza (calcolo della pendenza percentuale);
• accessibilità (delineazione delle aree accessibili al personale forestale);
• VBP (stima di Volume legnoso e Biomassa legnosa epigea Presenti);
• classi dimensionali (stima della consistenza percentuale del volume legnoso e della biomassa legnosa ripartiti in tre classi di dimensione dei fusti arborei);
• VBR (stima di Volume legnoso e Biomassa legnosa epigea Ritraibili, tenuto conto dei vincoli di accessibilità, normativi e gestionali);
• intervento colturale (quantificazione del volume legnoso e della biomassa legnosa ritraibili secondo le esigenze dell’utente);
• esbosco (definizione delle aree tecnicamente esboscabili sulla base del sistema di utilizzazione forestale).

Alcune funzionalità di base, come quella di rasterizzazione (trasparente all’utente) e quelle svolte dai moduli esposizione e pendenza, sono comuni a funzionalità analoghe presenti in tutti i software GIS. La scelta di implementare nella suite anche questi moduli basilari assicura all’utente il mantenimento di un unico ambiente di lavoro, garantendo al tempo stesso uniformità e caratteristiche degli elaborati intermedi per tutto il processo di analisi. Tutti i moduli lavorano su formato raster, con risoluzione spaziale a 30 metri tranne che per il modulo accessibilità, per il quale la risoluzione può essere invece definita dall’utente.

Nella prima finestra dell’interfaccia grafica (modulo area di interesse - Fig. 1) è richiesto all’utente di definire l’area di interesse secondo tre diverse modalità: (1) dalla extent di uno shapefile che l’utente fornisce; (2) dalla extent della view della finestra di QGIS; (3) sulla base di un input manuale da parte dell’utente (Fig. 3).

Fig. 3 - Esempio di output del modulo accessibilità in un’area di interesse definita all’interno della area test “Pollino”. La linea nera tratteggiata indica la viabilità disponibile, in scala di grigi i livelli previsti di accessibilità (1: “accessibile”; 2: “scarsamente accessibile”; 3: “non accessibile”).

 Enlarge/Shrink  Open in Viewer

Il modulo accessibilità permette di stabilire, sulla base del DTM e della viabilità, la potenzialità di accesso del personale forestale. Ad ogni pixel, sulla base della morfologia del terreno e dei tempi di percorrenza dalla rete stradale forestale, viene attribuita una delle seguenti tre classi di accesso: “accessibile”, “scarsamente accessibile” e “non accessibile” ([3]). Il modulo richiede come input un file vettoriale della viabilità presente nell’area di studio (file già presente in CFOR oppure fornibile dall’utente qualora ritenuto più preciso e/o più aggiornato) e il raster delle pendenze, che può essere generato anche utilizzando l’apposito tool messo a disposizione in CFOR.

Il modulo VBP fornisce una mappa raster della produzione legnosa dei soprassuoli forestali presenti nell’area in esame, sia in termini di provvigione legnosa (volume ad ettaro) che di biomassa epigea (peso secco della componente legnosa e totale e peso fresco della componente legnosa, ad ettaro). La provvigione e la biomassa sono stimate sulla base di modelli di regressione sviluppati a partire da metriche estratte dal Canopy Height Model (CHM) di derivazione LiDAR, calibrati con appositi rilievi al suolo. Per individuare il corretto modello di stima da applicare, CFOR utilizza una cartografia delle formazioni forestali opportunamente riclassificate secondo il sistema messo a punto da Scrinzi et al. ([13]), al quale si rinvia anche per i dettagli dei modelli di stima. Il modulo VBP permette, inoltre, di caricare anche dati LiDAR propri ed indipendenti qualora disponibili da parte dell’utente.

Con l’intento di permettere una valutazione orientantiva della possibile destinazione d’uso del materiale legnoso presente nell’area di interesse, è stato sviluppato uno specifico modulo chiamato classi dimensionali. Come per il precedente, anche questo modello è stato sviluppato calibrando relazioni empiriche tra il dato LiDAR e i dati rilevati a terra. In questo caso, seguendo un approccio basato su regressione lineare multipla, sono state messe in relazione alcune metriche estratte dal CHM con la distribuzione del diametro a petto d’uomo (dbh) dei fusti inclusi nei plot di campionamento a terra. In questo modo, è stato possibile sviluppare un modello di stima della percentuale di volume legnoso o di biomassa legnosa ripartita su tre grandi classi diametriche: fusti piccoli (dbh<12 cm), fusti medi (12 cm<dbh<28 cm), fusti grossi (dbh>28 cm).

Il modulo VBR sfrutta i raster prodotti dai due moduli precedenti (accessibilità e VBP) e li integra con strati informativi relativi a limitazioni d’uso della risorsa forestale potenzialmente disponibile; successivamente viene prodotta la mappa del volume legnoso e della biomassa legnosa effettivamente ritraibili per i boschi oggetto di interesse. Le limitazioni sono principalmente di carattere normativo, come ad esempio la presenza di vincoli protezionistici, o di tipo gestionale connessi alla forma di governo (ceduo o fustaia). Ad esempio, se per il ceduo la funzione offre la possibilità del “taglio a raso” con rilascio delle matricine secondo la normativa vigente, per la fustaia si è scelto di adottare come riferimento il principio della “provvigione minimale” ([4], [9]), criterio cautelativo che consente un uso sostenibile della risorsa forestale.

In aggiunta, per permettere all’utente di interagire nella definizione dei livelli di prelievo per le diverse classi dimensionali, è stato sviluppato uno specifico modulo chiamato intervento colturale, attivabile dall’utente. Tale modulo trova la sua migliore applicazione nell’ambito della pianificazione comprensoriale, dove stima delle masse e valutazione del legname ritraibile, al netto da vincoli, sono informazioni che servono alla definizione di linee di indirizzo per la gestione forestale multifunzionale di un’area vasta.

Infine, per supportare il gestore forestale nella valutazione della fattibilità tecnica e dell’onere, anche economico, derivante dalla asportazione del materiale legnoso, è stato implementato un modulo chiamato esbosco, deputato alla produzione della mappa delle aree tecnicamente esboscabili. Stabilita dall’utente un’area da esboscare, CFOR, dopo aver eseguito un aggancio automatico (snapping) al più vicino punto su strada (imposto), definisce un’area di pertinenza dell’imposto all’interno della extent. Sulla base della pendenza del terreno e della viabilità, il modulo realizza dapprima una mappa della distanza reale di ciascun pixel dall’imposto, applicando l’algoritmo r.cost di GRASS. Questa viene successivamente classificata in funzione dei parametri intrinseci (distanza massima, pendenza minima e massima) e di operabilità del sistema di esbosco scelto dall’utente: strascico con animali, trattore a strascico, trattore con verricello, gru a cavo. Incrociando questo prodotto con l’output del modulo intervento colturale, l’applicativo genera, inoltre, una tabella che riassume indicazioni sulle quantità di legname asportabili verso l’imposto e sui costi derivanti dalle diverse modalità di esbosco stabilite ([10]).

Aree test 

Tutti i moduli sono stati sperimentati in due aree del Pollino e delle Serre. Poiché gli strati informativi necessari al funzionamento dell’applicativo sono attualmente disponibili solo per le quattro grandi aree di studio del progetto ALForLab (Pollino, Serre, Sila, Catena costiera), il tool viene fornito come funzionante in tutte le sue parti solamente in questi contesti. Al di fuori di questi scenari, CFOR può comunque essere utilizzato dall’utente che, però, deve disporre dell’insieme delle informazioni necessarie e caricarle autonomamente nei diversi moduli.

Conclusioni 

I Sistemi di Supporto alle Decisioni Spaziali sono strumenti che possono favorire lo sviluppo e l’ammodernamento del settore operativo forestale. La presente nota illustra le caratteristiche architetturali e le metodologie di CFOR, un SDSS sviluppato per la regione Calabria nell’ambito del progetto ALForLab. Il tool, attraverso una serie di moduli, guida i professionisti forestali a decisioni di carattere tecnico allo scopo di favorire l’uso sostenibile delle risorse legnose.

I punti di forza di CFOR sono: semplicità di utilizzo; integrazione in un ambiente di lavoro GIS open-source, quale plug-in di QGIS liberamente distribuito; utilizzo di modelli di stima della provvigione legnosa su base LiDAR, permettendo di effettuare valutazioni di precisione su vaste aree; estrazione in tempo reale e in modo versatile di informazioni di carattere gestionale. In particolare, disponendo del dato LiDAR, e a condizione di mantenere la stessa risoluzione geometrica dei raster di output, sarebbe possibile usare CFOR per l’inventariazione delle masse evitando la fase dei rilievi dendrometrici di un piano di assestamento forestale.

Ringraziamenti 

Si ringraziano i signori Carlo Di Marco, Franco Calabrese, Gino Scarpelli, Massimo Salerno ed Ettore Salerno per il supporto tecnico. Il lavoro è stato finanziato dal Progetto “ALForLab” (PON03PE_000 24_1) cofinanziato dal Programma Operativo Nazionale Ricerca e Competitività (PON R&C) 2007-2013, attraverso il Fondo Europeo di Sviluppo Regionale (FESR) e risorse nazionali (Fondo di Rotazione (FDR) - Piano di Azione e Coesione (PAC) MIUR).

References