Segnale contro rumore: quando una variazione del carbonio del suolo è reale

Un acquirente chiede al consulente di un agricoltore di seminativi se il programma di colture di copertura stia sequestrando carbonio. Il consulente riporta +0,4 tonnellate di C per ettaro all'anno su cinque anni, con 60 carote appaiate per data di misura. Il verificatore chiede l'intervallo di confidenza. Comprende lo zero. L'affermazione viene silenziosamente abbandonata, il budget MRV della stagione viene registrato come perdita, e l'acquirente torna dal team della filiera senza nulla in mano.

Lo schema è ormai ben documentato. Bradford e colleghi (2023) hanno verificato la fattibilità del rilevamento empirico delle variazioni di SOC su 45 appezzamenti commerciali a seminativo e hanno riscontrato che il rilevamento, a livello di singolo appezzamento, di un accumulo significativo (circa 3 Mg C/ha in dieci anni) non è affidabile nemmeno con densità di campionamento di 1,2 ettari per campione1. Il segnale esiste. Semplicemente non è più grande del rumore alla scala che la maggior parte dei disegni di campo presuppone. Il crollo nel 2024 della clamorosa emissione di crediti australiani legati al pascolo, in cui circa 250.000 ACCU sono stati emessi in larga parte a causa di un'annata umida da decimo decile piuttosto che di un cambiamento di gestione, ha dimostrato la stessa cosa alla scala delle politiche pubbliche2.

Il rilevamento delle variazioni di SOC è un problema di rapporto segnale/rumore. Da un lato, tre fonti di rumore. Dall'altro, un piccolo segnale di gestione. La dimensione del campione, l'orizzonte temporale e le scelte di disegno discendono da questa disuguaglianza; non ne sono il punto di partenza.

Anatomia della soglia di rumore

Il rumore di misura del SOC è la varianza della variazione di stock stimata attorno al suo valore vero. Si somma a partire da tre fonti, e la più grande è raramente quella che preoccupa i professionisti.

La varianza spaziale all'interno dell'unità di inferenza

Il carbonio del suolo è eterogeneo a ogni scala. All'interno di un singolo appezzamento a seminativo dall'aspetto omogeneo, Poeplau e colleghi (2022) documentano, a scala di plot, deviazioni standard da 5 a 8 Mg C per ettaro sugli stock 0–30 cm anche quando si ricampiona lo stesso profilo con uno scostamento di posizione minimo, un coefficiente di variazione intorno al 10 per cento in gran parte irriducibile3. In paesaggi eterogenei, dove variano la posizione sul versante, il materiale parentale e la storia gestionale, CV dal 30 al 60 per cento sono la norma14. Nelle colture perenni tropicali, le deviazioni standard pubblicate per gli stock 0–30 cm nell'agroforestazione del cacao si concentrano intorno a 10–14 Mg C per ettaro, per effetto della posizione sul versante, dell'architettura dell'ombreggiamento e dell'età lungo la cronosequenza5.

La varianza spaziale è la fonte di rumore che la maggior parte dei progetti sottodimensiona nel budget, perché il passaggio pilota campiona tipicamente pochi ettari apparentemente uniformi. La varianza si dispiega quando il disegno estrapola all'intera unità di gestione. Se una nuova area ricada davvero nel dominio del modello è una questione correlata che affrontiamo in Rappresentativo di che cosa?

La varianza temporale e stagionale

Anche nello stesso punto della stessa azienda agricola, un campionamento ripetuto alla stessa profondità in date diverse restituisce valori diversi. Dominano due effetti. Primo, l'umidità del suolo e il pool di sostanza organica appena mineralizzata fluttuano con le precipitazioni, ed è per questo che l'emissione di crediti australiani è scomparsa nell'anno successivo alla fine delle condizioni di tripla La Niña2. Secondo, la profondità di campionamento sul campo è un'operazione stocastica, non un piano fisso: pochi centimetri di differenza tra baseline e ricampionamento possono far oscillare gli stock 0–30 cm di diversi Mg C per ettaro quando la densità apparente cambia con la lavorazione ridotta o con l'apporto di ammendanti organici. Wendt e Hauser (2013) hanno introdotto la correzione per massa di suolo equivalente (ESM) proprio per assorbire questo artefatto, e Raffeld et al. (2024) mostrano che la scelta tra profondità fissa ed ESM può far oscillare la variazione inferita più del segnale di gestione, a parità di dati grezzi67.

La varianza analitica e di processamento

Anche il laboratorio aggiunge varianza. Persino in un singolo laboratorio accreditato, le misure replicate sullo stesso campione presentano un CV dall'1 al 3 per cento per il carbonio totale mediante combustione a secco. La varianza tra laboratori e tra metodi (combustione a secco contro Walkley-Black, deriva di calibrazione del medio infrarosso (MIR), differenze nella preparazione dei campioni) è maggiore e raramente viene tracciata fino al dato di punta. Heuvelink e colleghi (2021) sostengono che la varianza di previsione dei modelli di machine learning può dominare una stima finale di variazione dello stock a scala spaziale fine, e che senza una propagazione esplicita di tutti e tre i termini di varianza un dato di variazione del SOC riportato è strutturalmente non interpretabile8.

In una tipica baseline su seminativo, le tre componenti contribuiscono rispettivamente per circa l'80, il 15 e il 5 per cento del totale. In un disegno appaiato baseline-ricampionamento, in cui si ritorna sugli stessi plot, σ²_spaziale crolla (la componente tra plot si annulla), e σ²_temporale più σ²_analitica diventano il rumore dominante. Il rapporto tra le due è la leva di disegno centrale, ed è il motivo per cui i disegni appaiati sono ormai lo standard nei regimi di crediti seri.

Il segnale che cerchiamo di rilevare

L'altro lato della disuguaglianza è il segnale indotto dalla gestione. La letteratura è finalmente abbastanza ricca da fondare intervalli realistici, e sono più piccoli di quanto ammetta la maggior parte delle presentazioni commerciali.

• Colture di copertura. La meta-analisi di Poeplau e Don su 37 siti trova 0,32 ± 0,08 Mg C per ettaro all'anno nello strato superficiale9. Onesto, reale, piccolo.
• Semina su sodo e lavorazione ridotta. Guadagni apparenti da 0,3 a 0,7 Mg C per ettaro all'anno in 0–15 cm scompaiono in gran parte quando si integra fino a 0–40 cm. Powlson e colleghi (2014) hanno mostrato che si tratta principalmente di un artefatto di ridistribuzione in profondità, non di sequestro netto10.
• Agroforestazione temperata. Mayer et al. (2022) riportano una media nello strato superficiale di 0,21 Mg C per ettaro all'anno su 51 studi, con una deviazione standard tra studi di 0,79; le siepi arrivano a 0,32, le colture in filari (alley cropping) a 0,26, il silvopascolo è leggermente negativo11. La DS tra studi è maggiore della media, e questo è il problema centrale.
• Agroforestazione del cacao. Adiyah et al. (2022, 2023) documentano un'ampia variazione per classe di età e posizione sul versante, con tassi di accumulo nello strato superficiale da 0,2 a 0,5 Mg C per ettaro all'anno e guadagni sostanziali nel sottosuolo nei popolamenti con più di ~18 anni5.
• Biochar e ammendanti. Eventi discreti e una tantum da 5 a 15 Mg C per ettaro seguiti da una lenta attenuazione; di fatto un gradino anziché una pendenza.
• Pascolamento adattivo a rotazione su più parcelle. Dipendente dal sito e dall'equilibrio; differenze medie dell'ordine di 0,5–1,5 Mg C per ettaro all'anno dove esistono, ma con un'ampia varianza tra siti.

Ne derivano due implicazioni. Il segnale di gestione realistico si colloca nell'intervallo da 0,1 a 0,5 Mg C per ettaro all'anno per le pratiche di tipo pendenza, più piccolo, in termini assoluti, dell'errore di laboratorio su un singolo campione. E la DS dell'effetto tra siti, nelle meta-analisi, è paragonabile alla media o maggiore, il che limita quanta compressione della varianza un singolo progetto possa ricavare dalla letteratura precedente.

La disuguaglianza di rilevabilità

Per un disegno appaiato baseline-ricampionamento tra due date, la più piccola variazione media che una campagna può distinguere statisticamente da zero è la differenza minima rilevabile (MDD, minimum detectable difference):

con z_α/2 = 1,96 ad α = 0,05 (bilaterale) e z_β = 0,84 all'80 per cento di potenza, per cui (z_α/2 + z_β)² ≈ 7,85. σ_appaiato è la deviazione standard della variazione per plot tra le due date (molto più piccola del σ assoluto dello stock se si ritorna sugli stessi plot e l'autocorrelazione intra-plot è alta). n è il numero di plot appaiati.

Tre proprietà di questa disuguaglianza contano per il disegno.

La MDD è simmetrica. Dimezzarla quadruplica l'n richiesto. Raddoppiare σ quadruplica anch'esso n. Non esiste un percorso graduale da una campagna che rileva 1 Mg C per ettaro all'anno a una che ne rileva 0,25; costa sedici volte il campione.

σ_appaiato crolla dove σ non crolla. È per questo che i disegni appaiati sono oggi la struttura credibile ai fini della verifica. Un disegno appaiato con 60 plot e σ_appaiato = 4 Mg C per ettaro può rilevare una variazione cumulata a 5 anni di circa 2 Mg C per ettaro con l'80 per cento di potenza. Gli stessi 60 plot, campionati in modo indipendente alla baseline e al ricampionamento, richiederebbero un σ più vicino a 2 Mg C per ettaro per raggiungere la stessa MDD, il che non è realizzabile su seminativi reali.

Il segnale si accumula con il tempo, non con il solo numero di campioni. Un accumulo di 0,3 Mg C per ettaro all'anno arriva a 1,5 Mg C per ettaro in cinque anni e a 3 Mg C per ettaro in dieci. Lo stesso n che non riesce a rilevare la variazione al quinto anno può rilevarla agevolmente al decimo senza aggiungere una sola carota. È il paradosso dell'orizzonte temporale: il modo più economico per acquistare potenza di rilevamento è aspettare.

Simulatore dell'inviluppo del segnale

Il primo diagramma rende concreta la disuguaglianza. Scelga un tasso di accumulo annuo, una deviazione standard appaiata e un numero di plot; il grafico traccia la traiettoria media vera e il corridoio di rumore a ±ε, e segna l'anno in cui la traiettoria esce per la prima volta dal corridoio. Abbassi σ con un disegno appaiato e l'anno di rilevamento si dimezza; alzi il tasso di accumulo e la traiettoria ne esce prima; riduca n e il corridoio si allarga finché la traiettoria non riesce più a uscirne entro 15 anni.

Il simulatore rende visibile perché la maggior parte delle affermazioni pubblicate del tipo «abbiamo sequestrato X» sia non verificabile sulla base dei propri stessi dati: i cursori che permettono di risolvere un segnale credibile si collocano al di fuori di ciò che il disegno di campo ha effettivamente finanziato.

Il paradosso dell'orizzonte temporale

L'implicazione della disuguaglianza della MDD è scomoda. Per la maggior parte delle pratiche di tipo rigenerativo, la finestra di rilevamento onesta è di 8–12 anni, non i 3–5 anni su cui sono tipicamente costruite le roadmap aziendali e i cicli di crediti112. Smith (2004) lo ha modellato due decenni fa: con un aumento del 20–25 per cento degli apporti di carbonio, il rilevamento con il 90 per cento di confidenza richiede all'incirca 6–10 anni se il regime tollera uno scostamento del 3 per cento del SOC di fondo, e non è in pratica mai rilevabile su cicli più brevi se il regime tollera solo il 15 per cento della variazione di fondo12.

Ne derivano tre conseguenze per la progettazione. Primo, accettare per i primi anni una contabilità ibrida di misura e modellazione anziché rivendicare troppo. La consultazione di Verra sulla VM0042 v3.0 (il Soil Sampling and Analysis Handbook pubblicato a marzo 2026) e l'argomentazione economica a favore di misura-e-rimisura di Potash et al. (2025) vanno entrambe in questa direzione1314. Secondo, strutturare il contratto attorno alla finestra di verifica, non a quella del marketing: un rapporto a 3 anni che dice «sotto il limite di rilevamento, traiettoria del segnale coerente con il modello» è più duraturo di un rapporto a 3 anni che dice «+1,4 Mg C per ettaro» con un IC che comprende lo zero. Terzo, prendere sul serio la questione della profondità. I guadagni di SOC nelle colture perenni a radicazione profonda con più di ~18 anni superano spesso quelli dello strato superficiale; un disegno 0–30 cm semplicemente non li vede57.

Il disegno di campionamento come leva di riduzione della varianza

Se σ è fissato dal suolo e l'attesa è vincolata dal contratto, la leva che resta è il disegno. Tre scelte di disegno si cumulano, e un progetto serio le usa tutte e tre.

Stratificazione. Allocare lo sforzo di campionamento a strati che catturano la varianza (topografia, anzianità della gestione, tipo di suolo) riduce il σ effettivo e quindi n. de Gruijter et al. (2016) hanno introdotto l'algoritmo Ospats per un'allocazione stratificata a costo ottimale a partire da una mappa a priori del SOC, e hanno mostrato che il guadagno è sostanziale quando l'informazione a priori è buona15. Empiricamente, lavori successivi sui dati di Bradford et al. mostrano che la stratificazione su covariate facilmente disponibili consente una riduzione media del 17 per cento della dimensione campionaria (IC 95% dall'11 al 23 per cento) rispetto al campionamento casuale semplice. Informazioni a priori solide unite a un campionamento doppiamente bilanciato possono portare questo valore al 30–50 per cento1. Attenzione alla soglia minima: la VM0042 v2.2 fissa un minimo di 15 campioni per strato, il che può annullare il guadagno di Neyman su piccole unità di gestione.

Disegno appaiato baseline-ricampionamento. Tornare sugli stessi plot al ricampionamento fa crollare la componente di varianza tra plot. Le riduzioni di rumore tipiche di un disegno appaiato vanno dal 30 al 70 per cento, a seconda dell'autocorrelazione temporale del SOC a scala di microsito. Il Soil Carbon Monitoring Network dell'USDA NRCS lo rende operativo con sei campioni per appezzamento di 4 ettari secondo un calendario a disegno appaiato, all'incirca una carota ogni 0,7 ettari16.

Compositi all'interno dello strato, indipendenza tra strati. Accorpi le carote all'interno di uno strato per assorbire la varianza a microscala in una stima della media di strato a costo inferiore; mantenga gli strati come repliche indipendenti, così che il test conservi gradi di libertà. L'opposto (accorpare tra strati in un unico composito) è rapido ed economico e inutile per l'inferenza.

Selettore del disegno di campionamento

Il secondo diagramma combina le tre leve di disegno. Un disegno di riferimento casuale semplice, con σ impostato dall'utente, produce una stima degli anni necessari al rilevamento. La stratificazione riduce il σ effettivo di un fattore realistico secondo la letteratura. Il disegno appaiato baseline-ricampionamento fa crollare ulteriormente la varianza tra plot. Il grafico aggiorna gli anni fino al rilevamento, il numero per strato dopo l'applicazione della soglia minima di 15 in stile VM0042 e il numero totale di carote che ne risulta.

Lo scopo del selettore non è ottimizzare verso un unico numero. È far emergere dove ciascuna scelta di disegno ripaga, e dove la soglia minima di verifica (il vincolo di 15 per strato, il requisito di profondità, la correzione ESM) costa più di quanto Neyman faccia risparmiare.

Che cosa cambia la svolta della fusione modello-dati

Il cambiamento più rilevante del 2024–2026 è stato il passaggio da crediti basati sulla sola misura o sul solo modello verso disegni ibridi. L'architettura è ormai ben descritta: modelli di Tier 3 basati sui processi o empirici, calibrati localmente con reti di misura appaiate, con la varianza di previsione e la varianza dei parametri propagate congiuntamente nell'intervallo di incertezza riportato1714. Il documento di consultazione di Verra sulla VM0042 v3.0 e il GHG Protocol Land Sector and Removals Standard (in vigore dal 1º gennaio 2027) richiedono entrambi dati empirici, specifici per ciascun serbatoio, per qualsiasi rimozione di carbonio del suolo riportata, chiudendo la scappatoia del «solo modello» sfruttata da diversi programmi di prima generazione1318. La revisione indipendente di CarbonPlan, rivolta agli acquirenti, ha concluso nel 2023 che nessun protocollo disponibile sul mercato soddisfaceva contemporaneamente requisiti adeguati di rigore, durabilità, addizionalità e profondità; il ciclo di standard del 2026 è, in parte, una risposta esplicita19.

Questo conta per l'impostazione segnale contro rumore. Un modello da solo fornisce una previsione deterministica con un limite di incertezza strutturalmente ottimistico. Una campagna di misura da sola è soggetta alla soglia di rumore descritta sopra. Un modello calibrato sulle misure può comprimere l'intervallo di incertezza al di sotto di entrambe le componenti, ma solo se le misure sono abbastanza dense da vincolare i parametri che contano e se l'errore strutturale del modello è incluso nella varianza propagata. Gli ibridi «modello leggero, misura leggera» non sono credibili; «la misura come calibrazione del modello» è lo standard operativo verso cui convergono i regimi del 2026.

Modelli di rendicontazione onesti

La disciplina che rende operativa un'impostazione segnale contro rumore sta nella rendicontazione. Quattro numeri dovrebbero accompagnare qualsiasi affermazione di variazione del SOC, e un verificatore dell'era 2026 li chiederà.

1. La stima puntuale della variazione, con unità e profondità (ad esempio «variazione del SOC 0–30 cm di +1,4 Mg C per ettaro in 5 anni, corretta ESM»).
2. L'intervallo di confidenza al 95 per cento calcolato dalla varianza effettiva del disegno, non da un CV generico tratto dalla letteratura.
3. La MDD che il disegno può risolvere con l'80 per cento di potenza, calcolata prima della campagna, e il risultato espresso come «affermazione entro / oltre il limite di rilevamento».
4. La tassonomia dell'incertezza propagata secondo le convenzioni del Volume 4 AFOLU dell'IPCC20: varianza dei dati di attività, varianza dello stock di riferimento e varianza della variazione di stock o dell'errore del modello, combinate in modo additivo.

Un'affermazione di +1,4 Mg C per ettaro con un IC da −0,2 a +3,0 e una MDD di 1,6 è onesta, utile e non verificabile; entra nel rapporto in quanto tale e il progetto pianifica la rimisura successiva. Un'affermazione di +1,4 Mg C per ettaro con un IC da +0,6 a +2,2 e una MDD di 0,8 è onesta, utile e difendibile; entra nell'inventario.

Che cosa ne consegue per la progettazione

L'impostazione segnale contro rumore semplifica diverse decisioni di acquisto che, negli ultimi cinque anni, sono state ripetutamente dibattute come se fossero una questione di opinione.

L'unità di inferenza è il progetto, non l'appezzamento. Il rilevamento per singolo appezzamento delle pratiche di tipo pendenza non è statisticamente affidabile a una densità sostenibile nei costi, mentre il rilevamento della media di progetto su 20–60 appezzamenti appaiati è operativamente gestibile. Progetti di conseguenza, e riconduca l'affermazione dell'acquirente al livello che i dati possono difendere.

La finestra di rilevamento onesta è l'orizzonte della gestione, non la roadmap aziendale. Per colture di copertura, semina su sodo e la maggior parte delle forme di agroforestazione, sono 8–12 anni dalla baseline. I cicli di rendicontazione triennali restano utili, purché riportino i risultati rispetto alla traiettoria della MDD anziché rispetto a un segnale presunto.

Il profilo profondo non è più facoltativo. I guadagni nel sottosuolo nei sistemi perenni con più di ~18 anni superano spesso quelli dello strato superficiale e sono invisibili a un disegno 0–30 cm5721. La correzione ESM è parimenti non negoziabile nei confronti tra lavorazione e semina su sodo67.

Il protocollo rispetto al quale l'acquirente viene sottoposto ad audit fissa la soglia minima. VM0042 v2.2/v3.0, il CRCF dell'UE (Regolamento 2024/3012), il GHG Protocol Land Sector and Removals Standard e SBTi FLAG convergono su requisiti di dati empirici più stringenti nel corso del 2026–202713182223. Progettare secondo lo standard più permissivo che il team conosce, quando l'acquirente sarà sottoposto ad audit rispetto al più rigoroso, è una modalità di fallimento al rallentatore.

Punti chiave