Come calcolare la quantità di carbonio in un albero? Metodi, incertezza e contabilità

Circola un dato molto popolare: un albero sequestrerebbe circa 22 chilogrammi di CO₂ all'anno. È rassicurante, semplice, e sbagliato abbastanza spesso da essere pericoloso come base per la contabilità del carbonio. La stessa specie arborea che cresce in due climi diversi può differire di un fattore tre. Un albero giovane e un albero maturo della stessa specie possono differire di un fattore dieci. La domanda non è «quanto», ma «come lo calcoliamo per questo albero, in questo luogo, a questa età, con un'incertezza difendibile attorno al numero?».

Il calcolo non ha nulla di misterioso. È una catena di quattro passaggi, più un quinto: la propagazione dell'incertezza, diventata non negoziabile con l'impianto normativo del 2026. Ogni passaggio ha un metodo difendibile e un punto ben noto in cui i professionisti prendono scorciatoie. Nella letteratura sottoposta a revisione paritaria, l'incertezza sulla biomassa epigea a scala di plot varia ancora dal 10 al 40%, e il solo errore allometrico a livello di albero può contribuire per il 30–75% all'incertezza totale delle stime di carbonio dal popolamento al paesaggio1. Nel 2026 la domanda dell'audit non è soltanto «qual è la Sua stima?», ma «qual è il Suo intervallo di confidenza, e come lo ha costruito?».

Passaggio 1: dalle dimensioni dell'albero alla biomassa epigea

La biomassa epigea (AGB) si stima a partire da dimensioni misurabili dell'albero, più comunemente il diametro a petto d'uomo (DBH, misurato a 1,3 m), l'altezza dell'albero (H) e la densità del legno (ρ). La forma funzionale dominante nella letteratura forestale tropicale è l'equazione allometrica pantropicale di Chave2:

Si tratta di un'equazione pantropicale, cioè adattata su molte specie e regioni. È un punto di partenza utilizzabile ma un punto d'arrivo pericoloso per le applicazioni ad alto rischio, perché il suo errore è sistematico anziché casuale. L'altezza è la colpevole abituale. Poiché la relazione altezza–diametro di una regione può discostarsi nettamente dalla media pantropicale, l'altezza degli alberi ha dovuto essere integrata esplicitamente nei modelli pantropicali di biomassa, anziché restare implicita nel diametro3. Nelle foreste dell'Amazzonia nord-orientale, la relazione pantropicale sovrastima la biomassa della foresta inondata (várzea) di circa il 17% e sottostima quella di terra firme di un margine analogo, mentre l'impiego di modelli altezza–diametro adattati localmente riduce la differenza a circa l'1%4. Il problema si aggrava sul campo: la misurazione di routine sottostima sistematicamente l'altezza degli alberi alti, e quell'errore si propaga direttamente nella biomassa5. La conseguenza è un paradosso familiare: una mappa pantropicale regionale può essere priva di distorsioni nel complesso, mentre ogni singolo progetto locale che essa copre è gravemente sbagliato.

La scansione laser aviotrasportata ha reso evidente la posta economica in gioco. In Guyana francese, sostituire l'allometria pantropicale dell'altezza con relazioni locali derivate dall'ALS ha quasi dimezzato l'errore medio sull'altezza e spostato le stime di AGB di 40–54 t/ha, uno scarto dell'11–13% superiore al margine tipico di un progetto di carbonio6. Le equazioni allometriche specifiche per specie o per sito superano costantemente la forma pantropicale, in particolare per le specie con un'architettura insolita (palme, felci arboree, specie con forti contrafforti) e per le specie agroforestali che non corrispondono alla struttura dominata dal fusto degli alberi di foresta naturale7.

Quando l'allometria specie-specifica vale il suo costo

Un'equazione allometrica su misura richiede un campionamento distruttivo: abbattere un insieme di alberi di riferimento che copra l'intervallo dimensionale del popolamento, pesare la biomassa fresca e misurare la frazione secca. Per una singola specie, sull'arco di vita di un progetto, si tratta in genere di 15–25 alberi. Il costo è reale, ma la riduzione dell'incertezza è di solito maggiore di quella di qualsiasi altro intervento in un progetto di carbonio arboreo, e spesso restringe l'intervallo di previsione attorno alle stime di biomassa a livello di popolamento del 40–60 per cento8.

Dai lavori recenti è emerso un punto più sottile: le equazioni allometriche dovrebbero essere scelte in base alla capacità predittiva a livello di plot, non in base all'RMSE o all'AIC a livello di albero. Gli errori residui si compensano in parte tra gli alberi di uno stesso plot, per cui il modello che meglio prevede il singolo albero spesso non è quello che meglio prevede il totale del plot, cioè la grandezza che conta per la contabilità del carbonio9. I progetti che scelgono i modelli in base alla statistica sbagliata stanno misurando la cosa sbagliata.

Passaggio 2: la densità del legno, il moltiplicatore trascurato

La densità del legno (ρ, in g/cm³) entra linearmente nelle equazioni di tipo Chave, per cui un errore del 10% su ρ si traduce direttamente in un errore del 10% sull'AGB. La densità varia tra le specie più di quasi ogni altro termine dell'equazione10. Una pioniera a rapido accrescimento come Ochroma pyramidale (balsa) ha una ρ di circa 0,15 g/cm³; un legno duro come Milicia excelsa ha una ρ di circa 0,65; alcune acacie dense superano 0,85. Un tasso generico di carbonio per albero che non si adegua alla ρ della specie presuppone una media, e in un sistema agroforestale plurispecifico questa ipotesi è spesso la principale fonte singola di errore sistematico nella stima dell'AGB.

Il Global Wood Density Database è il riferimento standard per i valori di ρ sottoposti a revisione paritaria, per circa 16.000 specie11, con aggiornamenti continui orientati all'agroforestazione. Quando non lo si usa bene, la posta in gioco è reale: appoggiarsi a una media di genere al posto della densità specifica della specie gonfia l'incertezza a livello di plot, e nei sistemi africani di miombo e di cacao le densità predefinite possono stimare in modo errato la biomassa di particolari gruppi di specie. Per le applicazioni prive di un valore a livello di specie, una ρ a livello di genere o di famiglia è un ripiego difendibile; una media globale non lo è.

Il cacao ne è un esempio chiaro. Le allometrie generiche e pantropicali stimano in modo errato l'AGB per albero di Theobroma cacao, e le equazioni specifiche per il cacao possono spostare in modo sostanziale le stime di sequestro a livello di popolamento: le differenze pubblicate vanno da pochi punti percentuali a decine di punti percentuali, a seconda del sito e dell'equazione di riferimento12. Questa dispersione ha implicazioni dirette per il mercato dei crediti dell'agroforestazione del cacao in Africa occidentale, dove la chioma d'ombra detiene la maggior parte del carbonio del sistema.

Passaggio 3: la biomassa ipogea e la frazione di carbonio

La biomassa epigea rappresenta all'incirca il 70–80 per cento della biomassa totale dell'albero. La biomassa ipogea (BGB) viene quasi universalmente stimata tramite un rapporto radici/fusto (root:shoot) anziché scavata, perché lo scavo su larga scala è impraticabile13. Il valore predefinito dell'IPCC per la foresta tropicale umida si è storicamente collocato intorno a 0,24–0,3714, ma una sintesi globale di dati a livello di singolo albero ha riancorato la disciplina: la vera media globale di R:S è più vicina a 0,25 ± 0,10, e diminuisce con le dimensioni dell'albero, per cui le stime costruite su piccoli plot e alberi giovani la sovrastimano sistematicamente15. Le foreste tropicali umide si collocano all'estremo inferiore (R:S ~0,20–0,24); i sistemi boreali e aridi sono più alti. Un unico valore predefinito su un portafoglio di progetti eterogeneo può introdurre un errore del 20–40% sul carbonio ipogeo. Nell'agroforestazione e nel ripristino delle terre aride, dove le radici possono avvicinarsi al 30–40% del carbonio totale, quell'errore è rilevante.

La conversione dalla biomassa al carbonio usa una frazione di carbonio (CF). Il valore predefinito dell'IPCC è 0,47, un numero che la chimica moderna del legno ha in gran parte mandato in pensione. Le sintesi globali rilevano che il carbonio volatile (terpeni, composti fenolici e altri composti persi durante l'essiccazione in stufa) è in media intorno all'1,4% della massa secca e sfugge sistematicamente ai metodi di laboratorio standard16. Su base reale, corretta per i volatili, le specie tropicali contengono in media circa il 48–52% di carbonio, con un intervallo da circa il 45% a oltre il 55% tra i taxa; le latifoglie temperate si collocano più vicino al 48%, le conifere al 50–52%17. Usare 0,47 dove il valore reale è 0,51 introduce una distorsione sistematica al ribasso di ~8% che si propaga direttamente nelle tonnellate accreditate. Le linee guida del GHG Protocol e il 2019 Refinement dell'IPCC indicano ora entrambi frazioni specifiche per gruppo di specie, e i database sottoposti a revisione paritaria (Doraisami et al., Thomas & Martin) dovrebbero essere il riferimento predefinito.

Il calcolatore qui sotto esegue l'intera catena su un singolo albero. Modifichi il diametro, l'altezza e la densità del legno e osservi come la stima epigea passa attraverso il rapporto radici/fusto e la frazione di carbonio fino a uno stock in chilogrammi di carbonio, e poi al suo equivalente in CO₂. Noti quanto direttamente la densità del legno scali il risultato: moltiplica la biomassa in modo esattamente proporzionale, per cui due alberi di dimensioni identiche ma di specie diverse possono contenere quantità di carbonio molto diverse.

Passaggio 4: dallo stock alla cattura annua

La cattura annua di carbonio è la derivata della curva di crescita dell'albero, non la sua media. Gli alberi non sequestrano carbonio a un ritmo costante, e i tassi di crescita iniziali variano enormemente a seconda della specie e del contesto. Un legno duro a crescita lenta può accumulare poca biomassa nei primi anni, mentre investe in radici e struttura, mentre le specie tropicali a rapido accrescimento (cacao, Terminalia, Albizia, Gliricidia) possono raggiungere due metri in due o tre anni con una cattura iniziale significativa. Ciò che resta costante tra le specie è la forma della curva, non il tasso di partenza: il sequestro raggiunge di solito il picco da qualche parte in una fase esponenziale prima di attenuarsi man mano che l'albero si avvicina alla chioma e alle dimensioni adulte. Usare un unico «tasso medio annuo» su un progetto di 30 anni colloca i crediti nel momento sbagliato, sovrastimando alcuni anni e sottostimandone altri, e questo è un campanello d'allarme per qualsiasi standard di carbonio serio.

I forestali dispongono di un vocabolario preciso per questo. L'incremento corrente annuo (CAI) è il carbonio che un albero aggiunge in un dato anno, la pendenza locale della sua curva di crescita. L'incremento medio annuo (MAI) è lo stock accumulato finora diviso per l'età, la media sull'intera vita fino a quel momento. Lo strumento qui sotto traccia entrambi in funzione dell'età per una curva di crescita di Chapman-Richards. Il CAI sale fino a un picco e poi ridiscende; il MAI raggiunge il picco più tardi, esattamente dove le due curve si incrociano, l'età che un forestale legge come il punto di crescita media più rapida. Un progetto che accredita un «tasso medio» costante sta prezzando il plateau del MAI e ignorando la gobba del CAI che sta sotto.

L'approccio difendibile consiste nello stimare lo stock di carbonio a due età utilizzando la catena allometrica descritta sopra, e nel rendicontare la variazione nell'intervallo come cattura incrementale. Per i progetti con orizzonti più lunghi, adattare ai dati di coorte una curva di crescita specifica per la specie (spesso nella forma di Chapman-Richards o di Von Bertalanffy) fornisce una funzione di cattura continua, anno per anno. È la differenza tra un profilo di assorbimento che corrisponde alla biologia e uno che è una linea retta per comodità.

Passaggio 5: il budget di incertezza

Nessuno dei quattro passaggi precedenti è sufficiente da solo. I verificatori non accettano più una stima puntuale dell'AGB; richiedono un intervallo di confidenza difendibile attorno a essa. La disciplina è convergita sulla propagazione dell'errore Monte Carlo come migliore pratica, implementata nel pacchetto R open source BIOMASS18 e nei suoi successori.

Un budget di incertezza rigoroso tiene oggi conto di quattro fonti: (1) l'errore di misura sull'albero, per DBH, altezza e identificazione della specie; (2) la variabilità residua del modello allometrico, la dispersione dell'adattamento originale di Chave; (3) l'incertezza sui parametri allometrici, la covarianza dei coefficienti adattati stessi; e (4) l'incertezza di campionamento, ossia quanto i plot misurati siano rappresentativi della popolazione. Due lezioni dovrebbero essere integrate in ogni piano di progetto. Primo, un'analisi dell'incertezza dall'albero al paesaggio ha rilevato che le sole equazioni allometriche contribuivano per il 30–75% al totale, spesso più dell'errore di telerilevamento, il che significa che calibrare l'allometria vale di solito più che perfezionare il satellite1. Secondo, ignorare la correlazione spaziale tra le osservazioni di alberi e plot gonfia la dimensione campionaria effettiva e sottostima l'incertezza reale, un'omissione comune e gravida di conseguenze.

Il budget qui sotto combina queste fonti in un unico numero. Imposti un coefficiente di variazione per la misura, per l'allometria (errore residuo più errore sui parametri) e per il campionamento, e legga l'intervallo combinato al 95% su uno stock nominale per ettaro. Poiché gli errori indipendenti si sommano in quadratura, la fonte singola più grande domina il totale: la barra della quota di varianza mostra perché rifinire un termine piccolo sposta appena l'intervallo, mentre ridurre il più grande, di solito l'allometria, è l'unica leva che lo restringe.

La Land Sector and Removals Guidance del GHG Protocol richiede esplicitamente che l'incertezza sia quantificata e resa nota19. I Core Carbon Principles dell'ICVCM penalizzano ora le metodologie che usano stime puntuali senza intervalli di confidenza. L'era del «±30% implicito» è finita.

Da un albero a un ettaro

Un progetto di carbonio non viene pagato per singoli alberi; viene pagato per un popolamento. Il calcolo per albero è l'operazione elementare, ma il dato rendicontato è uno stock per ettaro, costruito sommando gli stock degli alberi all'interno dei plot misurati ed estrapolando all'area del progetto. In questo passaggio cambiano due cose. Gli errori che appaiono allarmanti su un singolo albero si compensano in parte sull'insieme di un plot, perché un modello che sovrastima un albero grosso tende a sottostimarne uno sottile: è questa la vera ragione per cui le equazioni allometriche dovrebbero essere scelte per la previsione a livello di plot anziché per l'adattamento a livello di albero. E compare una fonte di incertezza che non ha nulla a che vedere con l'allometria: se la manciata di plot misurati sia davvero rappresentativa dell'intera area.

Questa questione di campionamento è una disciplina a sé. Quanti plot, di quali dimensioni e come disporli su un mosaico di età, mescolanze di specie e suoli: è questo a decidere se una stima a livello di popolamento è affidabile o è un valore puntuale con barre d'errore disegnate a posteriori. Le stesse statistiche che dimensionano una campagna sul suolo dimensionano una campagna sulla biomassa; affrontiamo direttamente il problema del quanti-e-dove in Quanto deve essere grande la Sua campagna di campionamento del carbonio del suolo? e la questione preliminare, cioè se un nuovo sito ricada davvero nel dominio su cui l'allometria è stata calibrata, in Rappresentativo di che cosa?

Dal suolo al cielo: il telerilevamento

Gli anni 2020 hanno trasformato il carbonio forestale da una disciplina basata sui plot a una disciplina basata su plot più satellite, con avvertenze importanti per gli sviluppatori di progetti:

• GEDI (il lidar spaziale della NASA) fornisce a scala globale l'altezza della chioma e la biomassa su un'impronta di 25 m, ma nell'agroforestazione dell'Africa occidentale le sue previsioni L4A presentavano errori circa 9× superiori all'AGB misurata sul campo20. Funziona bene nelle foreste a copertura chiusa sopra ~50 Mg/ha; per i sistemi di colture arboree dei piccoli produttori non è affidabile.
• La missione Biomass dell'ESA è stata lanciata il 29 aprile 2025. Il suo radar ad apertura sintetica in banda P (lunghezza d'onda di circa 70 cm) penetra le chiome per rilevare tronchi e grandi rami, dove si trova la maggior parte del carbonio forestale21. Test tomografici aviotrasportati in Guyana francese hanno ricostruito l'AGB con un errore inferiore a ~10% fino a 500 Mg/ha, un intervallo in cui i radar a lunghezza d'onda più corta saturano22. La missione esclude il Nord America e l'Europa a causa di restrizioni sullo spettro radar.
• La scansione laser terrestre con modellazione quantitativa della struttura (TLS + QSM) è l'attuale gold standard non distruttivo per l'AGB del singolo albero: una sintesi di studi di validazione distruttiva riporta una concordanza quasi unitaria (CCC ≈ 0,98) e un bias medio inferiore all'1% rispetto agli alberi abbattuti23. Sta rapidamente diventando il dataset di riferimento per la calibrazione satellitare.

Un'avvertenza vale per qualsiasi metodologia di crediti basata su mappe: la maggior parte delle mappe di biomassa ottenute con il machine learning riporta un R² elevato con la validazione incrociata casuale e fallisce con la validazione a blocchi spaziali, il che significa che l'accuratezza apparente è in gran parte un artefatto dell'autocorrelazione spaziale24. Gli acquirenti che valutano tali metodologie dovrebbero pretendere una validazione a blocchi spaziali prima di fidarsi dei numeri. Se un nuovo plot ricada davvero nel dominio calibrato di un modello è la questione preliminare che affrontiamo in Rappresentativo di che cosa?

Un esempio svolto: due specie, stesso plot, risposte diverse

Consideriamo due alberi d'ombra in un sistema agroforestale a cacao dell'Africa occidentale, entrambi di dieci anni, entrambi misurati nella stessa parcella.

Stesso plot, stessa età, una differenza di un fattore due nella cattura annua. La Terminalia a crescita rapida accumula volume velocemente; la Khaya a crescita lenta sta formando un legno denso che dominerà lo stock del popolamento al trentesimo anno. Un programma che applica a entrambe un unico tasso per albero sbaglia della metà su una delle due, e sovrastima o sottostima gli assorbimenti dichiarati a ogni ciclo di verifica.

Allargando lo sguardo di un livello emerge una seconda lezione. Gli stock di carbonio riportati nelle agroforeste a cacao coprono un ordine di grandezza, da circa 10 Mg C/ha nei sistemi giovani in pieno sole a ben oltre 100 Mg C/ha in quelli ombreggiati maturi, determinati molto più dalla chioma d'ombra che dal cacao25. Le piante di cacao in sé detengono di solito una minoranza del carbonio del sistema; il fattore dominante è la composizione degli alberi d'ombra e, soprattutto, la presenza di alberi relitti rimasti dalla foresta precedente, che possono contenere diverse volte il carbonio degli alberi d'ombra piantati. Ciò ha implicazioni dirette per l'addizionalità rispetto alla data limite del 31 dicembre 2020 del Regolamento UE sulla deforestazione e per il modo in cui i progetti segmentano i loro inventari di baseline, un punto che affrontiamo in Pochi alberi detengono la maggior parte del carbonio.

Che cosa richiede ora l'impianto normativo del 2026

Quattro quadri definiscono le regole contabili che gli assorbimenti basati sugli alberi devono rispettare. Il GHG Protocol Land Sector and Removals Standard e la relativa Guidance di accompagnamento (finalizzati nel 2026, applicabili alla rendicontazione dal 2027)19 sono le regole fondamentali di contabilità aziendale: rendicontazione separata delle emissioni lorde e degli assorbimenti lordi, tracciamento dettagliato dei cambiamenti di uso del suolo e quantificazione difendibile con incertezza resa nota su tutti e cinque i serbatoi di carbonio. SBTi FLAG v1.2 (marzo 2026) richiede una riduzione lineare annua del 3,03%, il tracciamento separato degli assorbimenti biogenici e un impegno di deforestazione zero. I Core Carbon Principles dell'ICVCM hanno ora approvato otto programmi di crediti come idonei CCP e 38 metodologie come approvate CCP, respingendone 22, con la metodologia ARR di ACR condizionata alla sola costituzione di foresta naturale. Il Regolamento UE sugli assorbimenti di carbonio e la carboniocoltura (2024/3012), in vigore dalla fine del 2024, dovrebbe adottare il proprio atto delegato sulle metodologie di carboniocoltura (imboschimento, agroforestazione, riumidificazione delle torbiere) nell'estate 2026, con le prime certificazioni alla fine del 2026.

Un progetto che voglia superare l'audit in tutti e tre i regimi deve ora fare sei cose contemporaneamente: usare equazioni allometriche calibrate localmente o specifiche per la specie, con errore residuo documentato; applicare densità del legno e frazione di carbonio specifiche per la specie, tratte da database sottoposti a revisione paritaria anziché dai valori predefiniti 0,47/0,50; propagare l'incertezza tramite Monte Carlo attraverso la misura degli alberi, l'allometria e il campionamento; quantificare il carbonio ipogeo con rapporti R:S stratificati per clima anziché con un unico valore predefinito; validare qualsiasi prodotto di telerilevamento con una validazione incrociata a blocchi spaziali; e separare gli assorbimenti dalle riduzioni, sia nella rendicontazione sia nella definizione degli obiettivi.

Punti chiave