Introduzione 

Le foreste hanno un ruolo importante nel ciclo globale del carbonio (C - [19]). Inoltre la biomassa presente negli alberi dei nostri boschi rappresenta, assieme a quella nel suolo, la frazione più importante di anidride carbonica (CO₂) fissata negli ecosistemi terrestri ([3]). Il C circola fra gli oceani, la biosfera e l’atmosfera. Le attività umane come lo sfruttamento di combustibili fossili e la deforestazione influenzano la concentrazione di CO₂ nell’atmosfera ([17]). Se ne deduce che una stima precisa della biomassa arborea (fitomassa) è essenziale per poter valutare il contributo delle foreste al bilancio planetario del C ([39]). La quantificazione del C in foresta deve essere fatta tenendo conto dei 5 pools indicati dall’United Nation Framework Convention on Climate Change ([55]): biomassa epigea, biomassa ipogea, legno morto, lettiera e C organico nel suolo (SOC, Soil Organic Carbon). La concentrazione di C nel suolo è più elevata che nell’atmosfera e ciò indica l’importanza di questo pool di C. Inoltre, vi sono studi, come quello di Tonolli & Salvagni ([54]) in Italia, che confermano, se pur con un certo margine di incertezza, il contributo quasi paritetico fra il SOC e la biomassa epigea. Nonostante ciò, a causa delle ancora scarse informazioni riguardo al SOC, è possibile considerare nullo il contributo allo scambio di C di questo pool ai fini della contabilità del Protocollo di Kyoto, purché si dimostri che esso non rappresenti una fonte di emissione di CO₂ in atmosfera. La capacità delle foreste di fissare la CO₂ atmosferica è fondamentale sia ai fini della conoscenza e della gestione degli ecosistemi, sia per gli aspetti economici che da essi derivano (e.g., mercato dei crediti di carbonio disciplinato dal PK). La contabilizzazione delle foreste come sink di C è una sfida, ma le stime comportano alti gradi di incertezza ([12]) e di discrepanza fra metodi ([16]). Houghton ([16]) suggerisce che le differenze fra le varie stime del sink di C ottenute con metodi diversi non sono dovute ad alti livelli di incertezza, ma semplicemente a conteggi incompleti. Secondo Bravo ([3]) gli strumenti attualmente riconosciuti dalla comunità scientifica per quantificare la fissazione del C nella fitomassa sono riconducibili a: 1) inventari forestali a piccola o a grande scala); 2) stime del cambio di superficie; 3) modelli empirici: tavole dendrometriche (regionali o locali), fattori di espansione della biomassa (Biomass Expansion Factors, BEFs), relazioni allometriche; 4) modelli empirici di crescita e/o di produzione (a scala di ecosistema, di popolamento o del singolo individuo); 5) telerilevamento. A queste vanno aggiunte le tecniche di micrometeorologia che fanno uso di torri (eddy covariance towers) per la valutazione dei flussi di C fra ecosistema ed atmosfera.

La fissazione o capacità di sink degli ecosistemi forestali è collegata con le attività addizionali indicate dall’articolo 3.4 del Protocollo di Kyoto (PK) che prevedono quale strumento attuativo la gestione forestale (Forest Management, FM). Secondo l’IPCC ([18]) le terre classificabili come soggette a FM sono quelle che corrispondono ai requisiti richiesti dalla definizione di bosco adottata in sede nazionale e su cui sono svolte le attività di gestione forestale indicate dal nostro Paese. Il metodo per la contabilizzazione di C per le attività di FM è detto gross-net accounting (Schlamadinger et al. 2007) e considera esclusivamente variazioni di stock di C dovute alle differenze tra emissioni (source) e prelievi (sink) durante un periodo di conteggio, senza che queste siano confrontate con le variazioni dello stock avvenuto in un periodo (o anno) di riferimento. Il sink totale generato da questa tipologia di attività addizionale può essere contabilizzato ai fini del mercato dei crediti solo sino ad esaurimento del valore assegnato all’Italia dalla Decisione 8/CMP.2 ([56]) e che è stato calcolato applicando un fattore di sconto dell’85% alle stime presentate dal Paese al Segretariato dell’UNFCCC. Sebbene un tale fattore di sconto non sia previsto dalle regole successive di calcolo dei sink forestali tuttavia è opportuno considerare il 15% dell’assorbimento netto come quantità utilizzabile per crediti di carbonio dalle attività di FM, anche in considerazione di possibili procedure di ricalcolo e aggiustamento.

Lo scopo di questo lavoro è quello di quantificare i crediti di CO₂ equivalenti ad ettaro (CO₂ eq ha^-1) generabili dalle attività di FM di un bosco di faggio nel Primo Periodo di Impegno (1st Commitment Period) 2008-2012 con il modello CO2FIX v. 3.1. Questo modello rappresenta una fra le più moderne applicazioni per lo studio dei flussi di carbonio fra ecosistemi forestali ed atmosfera, inoltre, esso consente di quantificare dinamicamente il sequestro potenziale di C, nei 5 pools stabiliti dall’UNFCCC.

Si stima inoltre la potenzialità dell’ecosistema oggetto dello studio, espressa come superficie di bosco necessaria per la produzione di crediti di carbonio da FM in quantità sufficiente a neutralizzare le emissioni CO2 di 2 tipologie di automobili (Box 1).

Materiali e metodi 

Area di studio

Lo studio è stato condotto in un bosco di faggio delle Serre calabre (Appennino Meridionale - Fig. 1) recentemente assestato ([37]). I boschi di faggio di quest’area possono essere riferiti, secondo la classificazione tipologica di Mercurio & Spampinato ([28]), alla tipologia “Faggeta ad agrifoglio” che corrisponde alle seguenti classificazioni Europee:

Fig. 1 - Localizzazione dell’area studio.

 Enlarge/Shrink  Open in Viewer

• 41.18 Southern Italian beech forests, CORINE ([7]);
• 9210* Apennine beech forests with Taxus and Ilex, Natura 2000 ([8]);
• 6.7.3 Apennine-Corsican montane beech forests, European Forest Types ([6]).

Si tratta di faggete meso-termofile che presentano numerose specie nemorali in comune con quelle europee (Sanicula europea, ecc.) ed altre invece esclusive e tipiche delle faggete dell’Italia meridionale (Anemome apennina ecc.). Lo strato arboreo è dominato dal faggio. Lo strato arbustivo e scarsamente sviluppato e rappresentato esclusivamente da Ilex aquifolium. Lo strato erbaceo è caratterizzato da specie indicatrici di buona fertilità quali Galium odoratum, Sanicula europea, ecc. In questa tipologia si riscontrano popolamenti con alberi di buona forma, che in alcuni casi superano i 30 m. I suoli sono in genere evoluti, profondi, a tessitura franco-sabbiosa, acidi, ricchi di sostanza organica con alta capacità di ritenuta idrica, riferibili ai Typic Dystrudept (Inceptisols, [57]). Secondo la classificazione di Rivas Martinez ([45]) queste faggete rientrano nel macrolima temperato, ombrotipo iperumido, termotipo mesomediterraneo con caratteri di transizione. I principali parametri dendroa-auxometrici vengono riportati nella Tab. 1.

Il modello CO2FIX v. 3.1

CO2FIX è un modello a scala di ecosistema (ecosystem-level model) che consente di calcolare il C in foresta della biomassa (fitomassa), del suolo e dei prodotti legnosi con cadenza annuale o per un periodo di riferimento ([27]).

Fig. 2 - I moduli di CO2FIX ([49] - ridisegnato). Le frecce rappresentano i flussi di C: in verde le perdite dell’atmosfera (Atm loss); in rosso gli incrementi nell’atmosfera (Atm inc).

 Enlarge/Shrink  Open in Viewer

Il modello è diviso in 6 moduli: biomassa (biomass), suolo (soil), prodotti (products), bioenergia (bioenergy), modulo finanziario (finance module) e di accreditamento (carbon accounting module - Fig. 2). In questo studio sono stati utilizzati i primi due moduli. CO2FIX fornisce risultati in forma grafica e di tabelle ed è stato scelto perché consente di: 1) calcolare e stimare lo stock di C; 2) calcolare e stimare i flussi di C fra ecosistema ed atmosfera e quindi la capacità di sink e di source quale conseguenze di un protocollo di gestione (utilizzazioni). Lo stock totale di C calcolato dal modello è dato dalla somma del C nella biomassa (C_bio cheinclude la necromassa in piedi), del C nel suolo (C_suolo che include la necromassa giacente a terra) e dal C nei prodotti legnosi di cui non si è tenuto conto in questo lavoro (C_prodotti - eqn. 1):

Lo stock di C e i flussi nella fitomassa vengono calcolati come bilancio fra l’incremento corrente (CAI, Current Annual Increment) da una parte, e turnover, mortalità e utilizzazioni dall’altra parte, tramite un approccio a coorti (cohort model approach - [43]). Una coorte può essere definita come un gruppo di alberi o di altre specie vegetali che presentano le stesse caratteristiche di accrescimento e per cui è possibile trattare come singola entità all’interno del modello ([58], [1]). Grazie a questo tipo di approccio è possibile applicare il modello a popolamenti plurispecifici e/o pluristratificati o tenere conto, nel conteggio dello stock di C, di altri tipi di vegetazione (ad es., vegetazione arbustiva, muschi ecc., che in casi particolari contribuiscono per buona parte allo stock totale di C di un ecosistema forestale). Lo stock totale (Cb) è dato dalla sommatoria dello stock di ogni singola coorte (Cb _it - eqn. 2):

dove Cb _it è il C stoccato nella fitomassa della coorte i al tempo t (Mg C ha^-1).

Per ogni nuovo periodo temporale (Cb _it +1) la fitomassa è calcolata come bilancio fra la fitomassa originaria (Cb _it), più l’incremento di fitomassa (Ib _it), meno il turnover dei rami, foglie e radici (T _it), meno la mortalità naturale (Mn _it), meno le utilizzazioni (U _it), meno la mortalità causata dalle operazioni di utilizzazione (Mu _it - eqn. 3):

dove Kc è una costante per trasformare la fitomassa in contenuto di C.

Per simulare Ib _it il modello usa come input il CAI che può essere calcolato attraverso le tavole alsometriche. Dal CAI viene calcolato il tasso di crescita per il compartimento foglie, rami e radici utilizzando coefficienti di allocazione. In questo studio i coefficienti di allocazione sono considerati costanti per il periodo di simulazione. Il modello offre due possibilità per parametrizzare l’accrescimento di ogni singola coorte: 1) in funzione dell’età del popolamento (CAI, tavole alsometriche); 2) in funzione della biomassa totale. Nel secondo caso l’accrescimento viene espresso come rapporto fra biomassa esistente e quella massima raggiungibile dall’ecosistema (utile per le foreste tropicali). In questo studio il CAI è considerato costante per il periodo di simulazione (Tab. 2).

La mortalità naturale (Mn _it), conteggiata dal modello come necromassa in piedi, può essere parametrizzata in funzione dell’età del popolamento o in funzione della biomassa totale per ogni singola coorte; la mortalità causata dalle operazioni forestali (Mu _it) viene parametrizzata come percentuale della biomassa che muore dopo le operazioni di utilizzazione. In questo studio la mortalità naturale è considerata nulla nel periodo di simulazione.

Il turnover (Tb_it) di ogni coorte è stimato come somma dei singoli turnover dei compartimenti foglie, rami e radici.

La biomassa asportata dalle utilizzazioni (U_it) viene parametrizzata come percentuale sottratta dalla biomassa totale nell’anno di riferimento nel periodo di simulazione. Questa può essere allocata nel modulo prodotti, inoltre, è possibile parametrizzare la quantità di residui di lavorazione che rimangono al suolo come percentuale sottratta dalla biomassa asportata.

Il SOC è stato valutato attraverso il modello YASSO incluso come modulo in CO2FIX ([27], [24]). Il modello, costituito da 5 compartimenti di decomposizione e 2 compartimenti per il legno morto giacente a terra (coarse woody debris - Fig. 3), è stato parametrizzato per stimare il SOC fino alla profondità di 1 metro di suolo minerale. I dati di base del modello, che vengono considerati costanti per il periodo di simulazione, sono: (1) i tassi di input di C nel suolo (produzione di lettiera, legno morto, radici) e (2) dati termo-pluviometrici (da cui dipendono i tassi di decomposizione della lettiera e del legno morto stimati da YASSO). YASSO nasce per fornire dati di riferimento a carattere generale, in quanto le conoscenze su questo indicatore sono ancora scarse: nonostante ciò è stato testato in regioni tropicali, temperate e boreali con esito positivo ([27], [10]).

Fig. 3 - Diagramma di flusso del modello YASSO. I box rappresentano i compartimenti di C mentre le frecce rappresentano i flussi di C ([49]).

 Enlarge/Shrink  Open in Viewer

Il modello CO2FIX stima i flussi di C annuali di entrata (Atm loss) e di uscita (Atm inc) dall’ecosistema a cadenza annuale come bilancio fra la produzione (biomassa prodotta, CAI) da una parte, e tassi di decomposizione del suolo (lettiera e legno morto giacente a terra, stimati dal modello YASSO) e dei prodotti (utilizzazioni) dall’altra. Per il C stoccato nei prodotti il modello fornisce 3 opzioni (Fig. 2): 1) può essere conteggiato per intero come rilasciato nell’atmosfera (questo studio); 2) viene calcolata la progressiva decomposizione dei prodotti durante il ciclo di vita e le conseguenti emissioni di CO₂ attraverso il modulo prodotti; 3) viene utilizzato il modulo della bioenergia che prende in considerazione le emissioni nell’atmosfera generate dallo sfruttamento dei prodotti come biomasse. La capacità di sink o di source dell’ecosistema viene calcolata attraverso i dati dei flussi annuali stimati dal modello per il periodo di simulazione come differenza fra le perdite dell’atmosfera (Atm loss _t = carbonio sottratto all’atmosfera al tempo t) e gli incrementi nell’atmosfera (Atm inc _t = carbonio immesso nell’atmosfera al tempo t - eqn. 4):

I valori di Atm inc e Atm loss rappresentano il bilancio netto di tutti i cambiamenti che avvengono a livello della biomassa, del suolo e dei prodotti mentre i dati dei singoli compartimenti (fusto, foglie, ecc.) forniti dal modello non tengono conto di tali cambiamenti. Ciò vuol dire che i valori di riferimento per applicare il metodo della differenza di stock di inventari successivi ([50]) sono quelli di Atm inc e Atm loss forniti da CO2FIX e non quelli dello stock dei compartimenti della fitomassa (fusto, foglie, ecc.; per ulteriori chiarimenti si rimanda a [49]).

Il sink totale del periodo di simulazione è dato dalla sommatoria dei sink annuali (sink _t - eqn. 5):

Secondo l’IPCC ([18]) nella valutazione del bilancio del C è necessario giustificare l’esclusione di stocks di carbonio in foresta. In questo studio sono state escluse dal conteggio tutte le piante con diametro a petto d’uomo inferiore a 2.5 cm, soglia minima di cavallettamento ([37]).

L’attendibilità della stima di CO2FIX dipende soprattutto da due aspetti: 1) dalla natura stocastica dei coefficienti stimati dal modello; 2) dalla mancanza o dalla cattiva natura dei dati di input utilizzati per la stima del C ([49]). Secondo Nabuurs & Mohren ([32]) per ottenere risultati attendibili è necessario: a) utilizzare il modello a scala locale; b) utilizzare il modello per popolamenti omogenei sotto il profilo strutturale; c) utilizzare dati a valenza locale ed aggiornati. Secondo Nabuurs et al. ([35]) i parametri di input che dimostrano di avere maggiore influenza sul sequestro di C calcolato tramite CO2FIX sono la densità basale del legno, il contenuto di C del legno ed il CAI. Queste considerazioni valgono soprattutto per il calcolo delle variazioni di stock della fitomassa.

Per ulteriori dettagli su CO2FIX e YASSO è possibile consultare il lavori di Mohren et al. ([30]), Nabuurs & Schelhaas ([34]), Masera et al. ([27]), Schelhaas et al. ([49]) e Liski et al. ([24]).

In questo studio il modello è stato testato a livello di compresa forestale utilizzando i dati dendrometrici e delle utilizzazioni forniti nel piano dei tagli del Piano di Gestione Forestale ([37]). In Tab. 3 viene riportata un breve rassegna dei modelli sviluppati con lo scopo di valutare quantitativamente e qualitativamente il C in foresta.

Dati e parametri di input del modello

I parametri di input per la stima del C nella fitomassa sono stati ottenuti attraverso l’analisi bibliografica ed analisi condotte direttamente in bosco (Tab. 2).

I parametri per il calcolo del C organico tramite il modello YASSO sono stati stimati attraverso: (a) l’analisi delle serie storiche di 5 stazioni termo-pluviometriche prossime all’area studio (per il calcolo del tasso di decomposizione della lettiera); (b) analisi bibliografica per calcolare i tassi di input di C nel suolo (Tab. 4).

Il calcolo dei crediti di CO ₂ eq

Il 15 % del valore totale del sink (Mg C ha^-1) è stato moltiplicato per il fattore di conversione (44/12 per la CO₂) in modo da ottenere i crediti di CO₂ equivalenti (Mg CO₂ eq ha^-1 - eqn. 6):

Risultati 

Il sink generato dall’ecosistema ed i CO₂ eq

Il modello indica che il C fissato dall’ecosistema nel periodo 2008-2012 dalle attività di FM è pari a 9.77 Mg C ha^-1 (in media 1.95 ± 0.91 Mg C ha^-1 anno^-1 nel periodo di simulazione - Tab. 5).

La dinamica del sink viene riportata nella Fig. 4 per due scenari: con e senza utilizzazioni. Nello scenario senza utilizzazioni il valore medio del sink nel periodo di simulazione è di 2.37 (± 0.2) Mg C ha^-1 anno^-1.

Fig. 4 - Dinamica del sink di C dell’ecosistema nel Primo Periodo di Impegno (2008-2012).

 Enlarge/Shrink  Open in Viewer

Il totale crediti di carbonio generabili nel primo periodo d’impegno dalle attività di FM (art. 3.4 del PK) è di 5.36 Mg CO₂eq ha^-1 mentre la media annua per il primo periodo è di 1.07 (± 0.5) Mg CO₂eq ha^-1 anno^-1. I valori non scontati sono riportati in (Tab. 6).

Lo stock e il sink di C dell’ecosistema

La sommatoria dello stock dei 5 pools di C è detta anche C ecosistemico (Ceco). L’output del modello indica che il Ceco varia durante il periodo di riferimento ed in media risulta essere pari a 285.71 (± 2.3) Mg C ha^-1 (Tab. 5). Il modello indica che il pool che contribuisce maggiormente allo stock di C è la biomassa epigea (fusto, rami e foglie: 54.89 %). Il solo compartimento del fusto contribuisce per il 45.52 % seguito da quello dei rami (8.63 %) e delle foglie (0.74 % - Tab. 5). Il pool suolo contribuisce per il 29.22 % (Tab. 5). La ripartizione del C nei vari pools viene riassunta nella Fig. 3.

Il SOC ottenuto dal modello per il periodo di simulazione è in media di 83.49 (± 1.29) Mg C ha^-1. La dinamica del SOC è diversa da quella del C nella biomassa. Nel periodo simulato esso mostra un trend decrescente mentre il C nella biomassa epigea tende ad aumentare (Fig. 5).

Fig. 5 - Dinamica dello stock di C nel suolo (SOC) e nel fusto per il periodo di simulazione.

 Enlarge/Shrink  Open in Viewer

Confronto fra volume atteso e volume osservato nell’anno di inizio simulazione

La massa dendrometrica attesa calcolata dal modello per la compresa boschi di faggio per il primo anno di simulazione (2008) è di 495.4 m³ ha^-1 mentre quella osservata per lo stesso anno è di 491.1 m³ ha^-1 ([37]).

Discussione 

I dati stimati dal modello CO2FIX indicano che nel periodo di simulazione la compresa forestale si comporta come sink in maniera differenziata nel tempo. L’effetto principale di tale variazione è da attribuire alle utilizzazioni programmate per l’anno 2010 (sink = 0.36 Mg C ha^-1 - Tab. 5, Fig. 4) in quanto il C asportato viene considerato come rilasciato nell’atmosfera. Queste osservazioni mettono in evidenza come la gestione forestale ed in particolare la selvicoltura influiscano sulla capacità di fissazione del C degli ecosistemi forestali. Il valore di assorbimento medio ottenuto dal modello, di 1.95 (± 0.91) Mg C ha^-1 anno^-1 nel periodo di simulazione, differisce di poco da quello riportato da Pregitzer & Euskirchen ([41]) di 1.7 (± 3.2) Mg C ha^-1 anno^-1 per le foreste temperate. Nabuurs & Schelhaas ([34]) hanno valutato il sink in 16 tipi forestali Europei dislocati in 10 Stati (Italia compresa, rappresentata dalle faggete dell’Appennino centrale, CAI = 7 m³ ha^-1 anno^-1) attraverso il software CO2FIX ottenendo un valore medio di culminazione della capacità di assorbimento di (1.15) 2.98 (4.1) Mg C ha^-1 anno^-1, che per la maggior parte delle tipologie si verifica dopo il 38esimo anno mentre i valori scendono a (0.13) 0.8 (1.4) Mg C ha^-1 anno^-1 dopo il 200esimo anno. A causa dell’alta variabilità dei dati, i valori forniti da Nabuurs e Schelhaas sono da considerare indicativi; inoltre gli autori non riportano i valori per singola tipologia. Mund ([31]) riporta per una faggeta trattata a tagli successivi in Germania valori di 4.6 Mg C ha^-1 anno^-1. Tonolli & Salvagni ([54]) riportano valori di 1.84 Mg C ha^-1 anno^-1 per una faggeta di 45 anni in provincia di Vicenza. Pilli et al. ([40]) riportano per una fustaia transitoria di faggio in provincia di Vicenza valori di 2.36 Mg C ha^-1 anno^-1 (questo valore è vicino a quello predetto dal modello in questo studio per il periodo di simulazione di 2.37 Mg C ha^-1 anno^-1 nello scenario senza utilizzazioni - Fig. 6, Tab. 7).

Fig. 6 - Ripartizione del C nei pools dell’ecosistema nel periodo di simulazione (media Mg C ha^-1).

 Enlarge/Shrink  Open in Viewer

Il Ceco ottenuto dal modello di 285.71 (± 2.3) Mg C ha^-1, nel periodo di simulazione, risulta maggiore rispetto ai valori riportati da Tonolli & Salvagni ([54]) di 176.2 (± 43.91) Mg C ha^-1 (per le faggete della Provincia di Trento), da Prentice ([42]) di 218 Mg C ha^-1, da Pregitzer & Euskirchen ([41]) di 239 (± 101) Mg C ha^-1 e da Mund ([31]) di 255.28 Mg C ha^-1. Pregitzer & Euskirchen ([41]) riportano che la variabilità di questo indicatore aumenta nei popolamenti di età avanzata (old growth forests).

Il valore medio del C della biomassa epigea ottenuto dal modello del periodo di simulazione di 156.8 (± 2.99) Mg C ha^-1 risulta: a) maggiore rispetto ai valori riportati da Pilli et al. ([40]), Visonà et al. ([59]), Joosten et al. ([20]) e Tonolli & Salvagni ([54]); b) minore rispetto ai valori riportati da Calamini et al. ([4]), Masci et al. ([26]) e Mund (2004 - Tab. 8). Dalla Tab. 8 si evince inoltre come i valori delle faggete di età comparabile ([26], [31]) sono molto simili. Il C nella biomassa epigea tende ad aumentare nel periodo di simulazione (Tab. 5) e ciò indica che le utilizzazioni previste non superano il saggio naturale di accrescimento. In un’ottica di lungo periodo questo tipo di gestione tende a far incrementare lo stock di C nella fitomassa.

Il valore medio del SOC nel periodo di simulazione ottenuto dal modello di 83.49 (± 1.29) Mg C ha^-1 è in linea con i valori riportati da Mund ([31]), Pregitzer & Euskirchen ([41]), Tonolli & Salvagni ([54]), mentre risulta minore ai valori riportati da Eswaran et al. ([9]) e Prentice (2001 - Tab. 9). La Fig. 5 e la Tab. 5 mettono in evidenza una diminuzione dello stock di C del suolo. Fermo restando che, in generale, i valori di questo indicatore, nonché di quelli forniti da CO2FIX, sono da considerare di larga massima e i dati di input del modello YASSO sono costanti nel periodo di simulazione, il trend negativo può essere attribuito a due aspetti che meritano ulteriori approfondimenti: 1) al regime di utilizzazioni; infatti, su questo aspetto vi sono studi ([5], [21], [23]) che dimostrano come il SOC diminuisce in conseguenza ad interventi selvicolturali; 2) alla risposta lenta del modello YASSO per raggiungere un equilibrium ; in questo caso avendo a disposizioni dati storici si sarebbe potuto inizializzare il modello partendo da un periodo precedente a quello di simulazione. Nonostante ciò, come già riportato, i valori medi forniti dal modello sono in linea con altri studi analoghi.

Il rapporto fra il C contenuto nella biomassa epigea e SOC (C biomassa epigea/SOC) in questo studio è uguale a 2.16; secondo i dati ottenuti da Mund ([31]) è di 1.9, mentre secondo i dati di Tonolli & Salvagni ([54]) è di 0.65. Ciò indica una grande variabilità del parametro.

Conclusioni 

Il modello CO2FIX si è dimostrato uno strumento valido per la stima del sink di carbonio di un popolamento forestale e dei relativi crediti generabili dalle attività di FM. La flessibilità di questo strumento si traduce: (a) nella relativa facilità di reperire i dati di input sia per il modulo della biomassa che del suolo; (b) nella possibilità di parametrizzare con precisione le attività di utilizzazione; (c) nella possibilità di inizializzare il modello senza la necessità di disporre di dati storici (purché vi sia corrispondenza fra i dati osservati e quella forniti dal modello). Gli accorgimenti necessari per l’uso del modello sono quelli di utilizzare dati sito-specifici ed aggiornati. La stima deve essere considerata a valenza locale considerando che il contenuto di C negli ecosistemi forestali è molto variabile. Il valore di sink fornito dal modello CO2FIX è dello stesso ordine di grandezza di altri lavori presenti in letteratura e dimostra che la compresa contribuisce positivamente (ma in maniera differenziata nel periodo di simulazione) alle attività di FM nel primo periodo d’impegno. I dati mettono in evidenza che, per il popolamento oggetto di questo studio, per mantenere positiva la fissazione di C dell’ecosistema forestale non si deve utilizzare una massa superiore a quella generata del saggio naturale di accrescimento. Nel caso di FM sarebbe opportuno poter conteggiare il carbonio fissato nei prodotti legnosi, in questo modo la quantità di C fissato aumenterebbe. Per fare questo sono necessari studi sui cicli di vita dei prodotti legnosi e strumenti attuativi quali la certificazione forestale (gestione forestale sostenibile e di filiera). In questo caso per i proprietari forestali, oltre al costo del Piano di Gestione Forestale, si deve sommare quello della certificazione che potrebbe essere ammortizzato nel lungo periodo, visti i valori effettivamente quantificabili ai fini dell’FM, dalla vendita di crediti di C. In questo senso la possibilità di sfruttare la capacità di assorbimento di C degli ecosistemi forestali, da parte di privati e delle amministrazioni, diventa un investimento di lungo periodo che potrebbe avere dei risvolti positivi o negativi a seconda delle decisione politiche e delle tendenze del mercato. Dal 1 gennaio 2008 è partita la contabilizzazione degli assorbimenti di carbonio con l’istituzione del Registro dei Serbatoi di Carbonio Agro-Forestali. Viste le prospettive future, in questa fase è opportuno indirizzare le azioni di ricerca verso la sperimentazione di sistemi selvicolturali e forme di trattamento che consentono la massimizzazione della capacità di fissazione di C dei popolamenti forestali mentre, per i proprietari e le amministrazioni, è opportuno prepararsi ad un’eventuale entrata “in gioco”, munendosi degli strumenti di pianificazione e di certificazione disponibili. In questo senso l’accorpamento delle proprietà in strumenti di gestione quali il consorzio forestale consente di contenere i costi. Come per tutte le attività imprenditoriali, c’è un margine di rischio da tenere in considerazione.

Per quanto riguarda il pool suolo esso è stato scarsamente considerato nel passato, in quanto privo di interesse economico, e perciò la carenza di dati e conoscenze è maggiore rispetto alla fitomassa. Le attuali azioni di mitigazione dei cambiamenti climatici hanno rivalutato il suolo come sink e stock di C e per questo motivo è necessario approfondire le conoscenze a riguardo.

Per quanto riguarda l’applicazione di modelli per la simulazione delle quantità di carbonio residente e dei flussi tra i vari serbatoi di carbonio, si noti che nel caso di piantagioni come ad esempio quelle all’interno delle attività denominate AR-CDM (Afforestation/Reforestation-Clean Development Mechanism) le stime ottenute dai modelli sono meno affette da errore.

Ringraziamenti 

Grazie a due anonimi revisori per i loro commenti utili alla stesura di questo manoscritto ed al Dr. Palosuo Taru per le informazioni fornite sul modello YASSO.

References