Sinal versus ruído: quando uma mudança de carbono no solo é real

Um comprador pergunta ao consultor de um produtor de lavoura se o programa de culturas de cobertura está sequestrando carbono. O consultor relata +0.4 tonelada de C por hectare por ano ao longo de cinco anos, com 60 testemunhos pareados por ponto no tempo. O verificador pede o intervalo de confiança. Ele cruza o zero. A afirmação é abandonada em silêncio, o orçamento de MRV da safra é dado como perdido, e o comprador volta à equipe de cadeia de suprimentos sem nada a mostrar.

O padrão já está bem documentado. Bradford e colegas (2023) testaram a viabilidade da detecção empírica de mudança no COS em 45 talhões comerciais de lavoura e constataram que a detecção, no talhão individual, de um acúmulo relevante (cerca de 3 Mg C/ha ao longo de dez anos) é pouco confiável mesmo em densidades de amostragem de 1.2 hectare por amostra1. O sinal existe. Ele apenas não é maior do que o ruído na escala que a maioria dos desenhos de campo pressupõe. O colapso, em 2024, da emissão australiana de créditos de pastagem que rendeu manchetes, na qual cerca de 250.000 ACCUs foram emitidos em grande parte por causa de um ano úmido de Decil 10, e não de mudança de manejo, fez o mesmo apontamento em escala de política pública2.

A detecção de mudança no COS é um problema de relação sinal-ruído. De um lado, três fontes de ruído. Do outro, um pequeno sinal de manejo. O tamanho amostral, o horizonte temporal e as escolhas de desenho decorrem dessa desigualdade; não são o ponto de partida.

Anatomia do piso de ruído

O ruído de medição do COS é a variância da mudança de estoque estimada em torno do seu valor verdadeiro. Ele se acumula a partir de três fontes, e a maior raramente é aquela com que os profissionais se preocupam.

Variância espacial dentro da unidade de inferência

O carbono do solo é heterogêneo em todas as escalas. Dentro de um único talhão de lavoura de aparência homogênea, Poeplau e colegas (2022) documentam desvios-padrão à escala da parcela de 5 a 8 Mg C por hectare nos estoques de 0–30 cm, mesmo ao reamostrar o mesmo perfil com um deslocamento posicional mínimo, um coeficiente de variação em torno de 10 por cento que é, em grande parte, irredutível3. Em paisagens heterogêneas, onde a posição na encosta, o material de origem e o histórico de manejo variam, CVs de 30 a 60 por cento são rotineiros14. Em perenes tropicais, os desvios-padrão publicados para os estoques de 0–30 cm em agrofloresta de cacau concentram-se em torno de 10 a 14 Mg C por hectare, determinados pela posição na encosta, pela arquitetura de sombreamento e pela idade na cronossequência5.

A variância espacial é a fonte de ruído para a qual a maioria dos projetos reserva orçamento insuficiente, porque a passagem-piloto normalmente amostra alguns poucos hectares de aparência uniforme. A variância se revela quando o desenho extrapola para a unidade de manejo completa. Se uma nova área está sequer dentro do domínio do modelo é uma questão relacionada que abordamos em Representativo de quê?

Variância temporal e sazonal

Mesmo no mesmo ponto da mesma fazenda, a reamostragem na mesma profundidade em datas diferentes devolve valores diferentes. Dois efeitos dominam. Primeiro, a umidade do solo e o reservatório de matéria orgânica recém-mineralizada flutuam com a chuva, e é por isso que a emissão australiana de créditos desapareceu no ano seguinte ao fim das condições de tríplice La Niña2. Segundo, a profundidade de amostragem no campo é uma operação estocástica, não um plano fixo: uns poucos centímetros de diferença entre a linha de base e a reamostragem podem alterar os estoques de 0–30 cm em vários Mg C por hectare quando a densidade aparente muda sob preparo reduzido ou aplicação de condicionador orgânico. Wendt e Hauser (2013) introduziram a correção por massa de solo equivalente (MSE) justamente para absorver esse artefato, e Raffeld et al. (2024) mostram que a escolha entre profundidade fixa e MSE pode alterar a mudança inferida em mais do que o próprio sinal de manejo, a partir dos mesmos dados brutos67.

Variância analítica e de processamento

O laboratório também acrescenta variância. Mesmo em um único laboratório credenciado, medições em réplica sobre a mesma amostra carregam um CV de 1 a 3 por cento para o carbono total por combustão a seco. A variância entre laboratórios e entre métodos (combustão a seco vs. Walkley-Black, deriva de calibração no infravermelho médio (MIR), diferenças no preparo das amostras) é maior e raramente é rastreada até o número final. Heuvelink e colegas (2021) defendem que a variância de predição de modelos de aprendizado de máquina pode dominar uma estimativa final de mudança de estoque em escala espacial fina, e que, sem a propagação explícita dos três termos de variância, um número reportado de mudança no COS é estruturalmente ininterpretável8.

Em uma linha de base típica de lavoura, as três contribuem com cerca de 80, 15 e 5 por cento do total. Em um desenho pareado de linha de base e reamostragem, no qual as mesmas parcelas são revisitadas, σ²_espacial colapsa (a componente entre parcelas se cancela), e σ²_temporal mais σ²_analítica passam a ser o ruído dominante. A razão entre as duas é a alavanca central do desenho, e a razão pela qual desenhos pareados são hoje o padrão em regimes sérios de crédito.

O sinal que tentamos detectar

O outro lado da desigualdade é o sinal impulsionado pelo manejo. A literatura enfim é rica o bastante para embasar faixas realistas, e elas são menores do que a maioria das apresentações de marketing admite.

• Culturas de cobertura. A metanálise de Poeplau e Don com 37 sítios encontra 0.32 ± 0.08 Mg C por hectare por ano no solo superficial9. Honesto, real, pequeno.
• Plantio direto e preparo reduzido. Ganhos aparentes de 0.3 a 0.7 Mg C por hectare por ano em 0–15 cm em grande parte desaparecem quando integrados até 0–40 cm. Powlson e colegas (2014) mostraram que isso é, sobretudo, um artefato de redistribuição em profundidade, e não sequestro líquido10.
• Agrofloresta temperada. Mayer et al. (2022) relatam uma média no solo superficial de 0.21 Mg C por hectare por ano ao longo de 51 estudos, com um desvio-padrão entre estudos de 0.79; as cercas-vivas alcançam 0.32, o cultivo em aleias 0.26, e o sistema silvipastoril fica ligeiramente negativo11. O DP entre estudos é maior do que a média, o que constitui o problema central.
• Agrofloresta de cacau. Adiyah et al. (2022, 2023) documentam grande variação por classe de idade e posição na encosta, com taxas de acúmulo no solo superficial de 0.2 a 0.5 Mg C por hectare por ano e ganhos substanciais no subsolo em povoamentos com mais de ~18 anos5.
• Biochar e condicionadores. Eventos pontuais e isolados de 5 a 15 Mg C por hectare seguidos de atenuação lenta; na prática, um degrau em vez de uma rampa.
• Pastejo adaptativo com múltiplos piquetes. Dependente do sítio e do equilíbrio; diferenças médias da ordem de 0.5 a 1.5 Mg C por hectare por ano onde existem, mas com grande variância entre sítios.

Duas implicações decorrem daí. O sinal realista de manejo está na faixa de 0.1 a 0.5 Mg C por hectare por ano para práticas do tipo rampa, menor, em termos absolutos, do que o erro de laboratório em uma única amostra. E o DP entre sítios do efeito, nas metanálises, é comparável ou maior do que a média, o que limita quanta compressão de variância um único projeto consegue extrair da literatura prévia.

A desigualdade de detectabilidade

Para um desenho pareado de linha de base e reamostragem entre dois pontos no tempo, a menor mudança média que uma campanha consegue distinguir estatisticamente de zero é a diferença mínima detectável (MDD):

com z_α/2 = 1.96 em α = 0.05 (bilateral) e z_β = 0.84 a 80 por cento de poder, de modo que (z_α/2 + z_β)² ≈ 7.85. σ_pareado é o desvio-padrão da mudança por parcela entre os pontos no tempo (muito menor do que o σ absoluto do estoque se as mesmas parcelas forem revisitadas e a autocorrelação dentro da parcela for alta). n é o número de parcelas pareadas.

Três propriedades dessa desigualdade importam para o desenho.

A MDD é simétrica. Reduzi-la à metade quadruplica o n necessário. Dobrar o σ também quadruplica o n. Não há caminho suave de uma campanha que detecta 1 Mg C por hectare por ano para outra que detecta 0.25; isso custa dezesseis vezes a amostra.

σ_pareado colapsa onde σ não colapsa. É por isso que os desenhos pareados são hoje a estrutura crível para verificação. Um desenho pareado de 60 parcelas com σ_pareado = 4 Mg C por hectare consegue detectar uma mudança cumulativa em 5 anos de cerca de 2 Mg C por hectare com 80 por cento de poder. As mesmas 60 parcelas, amostradas de forma independente na linha de base e na reamostragem, precisariam de um σ mais próximo de 2 Mg C por hectare para alcançar a mesma MDD, o que não é atingível em lavoura real.

O sinal se acumula com o tempo, não apenas com a contagem de amostras. Um acúmulo de 0.3 Mg C por hectare por ano acumula-se até 1.5 Mg C por hectare em cinco anos e 3 Mg C por hectare em dez. O mesmo n que não detecta mudança no ano cinco pode detectá-la com folga no ano dez sem acrescentar um único testemunho. Esse é o paradoxo do horizonte temporal: a forma mais barata de comprar poder de detecção é esperar.

Simulador do envelope do sinal

O primeiro diagrama torna a desigualdade concreta. Escolha uma taxa anual de acúmulo, um desvio-padrão pareado e um número de parcelas; o gráfico desenha a trajetória média verdadeira e o corredor de ruído em ±ε, e marca o ano em que a trajetória ultrapassa o corredor pela primeira vez. Reduza σ com um desenho pareado e o ano de detecção cai pela metade; aumente a taxa de acúmulo e a trajetória ultrapassa o corredor mais cedo; corte o n e o corredor se alarga até que a trajetória não consiga escapar dele dentro de 15 anos.

O simulador mostra por que a maioria das afirmações publicadas do tipo “sequestramos X” é inverificável nos próprios dados: os controles que resolvem um sinal crível ficam fora daquilo que o desenho de campo de fato custeou.

O paradoxo do horizonte temporal

A implicação da desigualdade da MDD é incômoda. Para a maioria das práticas de estilo regenerativo, a janela honesta de detecção é de 8 a 12 anos, não os 3 a 5 anos sobre os quais os planos corporativos e os ciclos de crédito costumam ser construídos112. Smith (2004) modelou isso há duas décadas: sob um aumento de 20 a 25 por cento nos aportes de carbono, a detecção com 90 por cento de confiança requer cerca de 6 a 10 anos se o regime tolerar uma variação de 3 por cento no COS de fundo, e é essencialmente indetectável em ciclos mais curtos se o regime tolerar apenas 15 por cento da variação de fundo12.

Três consequências decorrem para o desenho de projeto. Primeiro, aceitar uma contabilidade híbrida de medição e modelagem para os primeiros anos, em vez de exagerar na afirmação. A consulta do VM0042 v3.0 da Verra (o Soil Sampling and Analysis Handbook, publicado em março de 2026) e o argumento econômico de medição e remedição de Potash et al. (2025) ambos apontam nessa direção1314. Segundo, estruturar o contrato em torno da janela de verificação, não da de marketing: um relatório de 3 anos que diz “abaixo do limite de detecção, trajetória do sinal coerente com o modelo” é mais durável do que um relatório de 3 anos que diz “+1.4 Mg C por hectare” com um IC que cruza o zero. Terceiro, levar a sério a questão da profundidade. Os ganhos de COS em perenes de raízes profundas com mais de ~18 anos muitas vezes superam os do solo superficial; um desenho de 0–30 cm simplesmente não os enxerga57.

O desenho amostral como alavanca de redução de variância

Se σ é fixado pelo solo e esperar mais tempo é restringido pelo contrato, a alavanca que resta é o desenho. Três movimentos de desenho se compõem, e um projeto sério usa os três.

Estratificação. Alocar o esforço amostral a estratos que capturam variância (topografia, idade de manejo, tipo de solo) reduz o σ efetivo e, portanto, o n. de Gruijter et al. (2016) introduziram o algoritmo Ospats para alocação estratificada de custo ótimo usando um mapa prévio de COS, e mostraram que o ganho é substancial quando a informação prévia é boa15. Empiricamente, trabalhos posteriores sobre os dados de Bradford et al. mostram que estratificar por covariáveis prontamente disponíveis rende uma redução média de 17 por cento no tamanho amostral (IC 95% de 11 a 23 por cento) em relação à amostragem aleatória simples. Priores fortes somados à amostragem duplamente balanceada podem elevar isso a 30 a 50 por cento1. Observe o piso: o VM0042 v2.2 estabelece um mínimo de 15 amostras por estrato, o que pode apagar o ganho de Neyman em unidades de manejo pequenas.

Linha de base e reamostragem pareadas. Retornar às mesmas parcelas na reamostragem colapsa a componente entre parcelas da variância. As reduções de ruído típicas de desenhos pareados são de 30 a 70 por cento, a depender da autocorrelação temporal do COS de microssítio. A Soil Carbon Monitoring Network do USDA NRCS operacionaliza isso com seis amostras por talhão de 4 hectares em um cronograma de desenho pareado, cerca de um testemunho por 0.7 hectare16.

Composição dentro do estrato, independência entre estratos. Combine testemunhos dentro de um estrato para absorver a variância de microescala em uma estimativa de média do estrato a menor custo; mantenha os estratos como réplicas independentes para que o teste ainda tenha graus de liberdade. O oposto (combinar entre estratos em uma única amostra composta) é rápido, barato e inútil para inferência.

Painel do desenho amostral

O segundo diagrama compõe as três alavancas de desenho. Um desenho de base aleatório simples, com σ definido pelo usuário, produz uma estimativa de anos até a detecção. A estratificação reduz o σ efetivo por um fator realista segundo a literatura. A linha de base e reamostragem pareadas colapsam ainda mais a variância entre parcelas. O gráfico atualiza os anos até a detecção, a contagem por estrato após aplicado o piso de 15 no estilo VM0042, e a contagem total de testemunhos resultante.

O objetivo do painel não é otimizar até um único número. É revelar onde cada escolha de desenho se paga, e onde o piso de verificação (a restrição de 15 por estrato, a exigência de profundidade, a correção MSE) custa mais do que Neyman economiza.

O que a virada da fusão modelo-dados muda

A mudança mais consequente em 2024–2026 foi a passagem do crédito apenas por medição ou apenas por modelo rumo a desenhos híbridos. A arquitetura já está bem descrita: modelos de Nível 3 baseados em processos ou empíricos, calibrados localmente com redes de medição pareadas, com a variância de predição e a variância de parâmetros propagadas conjuntamente para o intervalo de incerteza reportado1714. O documento de consulta do VM0042 v3.0 da Verra e o GHG Protocol Land Sector and Removals Standard (em vigor a partir de 1º de janeiro de 2027) ambos exigem dados empíricos, específicos por reservatório, para qualquer remoção de carbono do solo reportada, encerrando a brecha do “apenas modelo” que vários programas de primeira geração usaram1318. A revisão independente do guia do comprador da CarbonPlan concluiu, em 2023, que nenhum protocolo comercialmente disponível atendia simultaneamente aos requisitos adequados de rigor, durabilidade, adicionalidade e profundidade; o ciclo de normas de 2026 é, em parte, uma resposta explícita19.

Isso importa para o enquadramento de sinal versus ruído. Um modelo sozinho tem uma predição determinística com um limite de incerteza estruturalmente otimista. Uma campanha de medição sozinha tem o piso de ruído descrito acima. Um modelo calibrado por medições consegue comprimir o intervalo de incerteza abaixo de qualquer dos componentes, mas apenas se as medições forem densas o bastante para restringir os parâmetros que importam e se o erro estrutural do modelo estiver incluído na variância propagada. Híbridos do tipo “modelo leve, medição leve” não são críveis; a “medição como calibração do modelo” é o padrão operacional para o qual os regimes de 2026 estão convergindo.

Modelos de relato honesto

A disciplina que torna operacional o enquadramento de sinal versus ruído está no relato. Quatro números devem viajar juntos com qualquer afirmação de mudança no COS, e um verificador da era de 2026 vai pedi-los.

1. A estimativa pontual da mudança, com unidades e profundidade (por exemplo, “mudança de COS em 0–30 cm de +1.4 Mg C por hectare ao longo de 5 anos, corrigida por MSE”).
2. O intervalo de confiança de 95 por cento calculado a partir da variância real do desenho, não de um CV genérico da literatura.
3. A MDD que o desenho consegue resolver a 80 por cento de poder, calculada antes da campanha, com o resultado expresso como “afirmação dentro / fora do limite de detecção”.
4. A taxonomia de incerteza propagada seguindo as convenções do Volume 4 (AFOLU) do IPCC20: variância dos dados de atividade, variância do estoque de referência, e variância de mudança de estoque ou de erro do modelo, combinadas de forma aditiva.

Uma afirmação de +1.4 Mg C por hectare com um IC de −0.2 a +3.0 e uma MDD de 1.6 é honesta, útil e inverificável; entra no relatório como tal e o projeto planeja a próxima remedição. Uma afirmação de +1.4 Mg C por hectare com um IC de +0.6 a +2.2 e uma MDD de 0.8 é honesta, útil e defensável; entra no inventário.

O que decorre para o desenho de projeto

O enquadramento de sinal versus ruído simplifica várias decisões de aquisição que, ao longo dos últimos cinco anos, foram repetidamente disputadas como se fossem questão de opinião.

A unidade de inferência é o projeto, não o talhão. A detecção por talhão de práticas do tipo rampa é estatisticamente pouco confiável em densidade acessível, ao passo que a detecção da média do projeto ao longo de 20 a 60 talhões pareados é operacionalmente tratável. Faça o desenho de acordo, e direcione a afirmação do comprador ao nível que os dados conseguem defender.

A janela honesta de detecção é o horizonte de manejo, não o plano corporativo. Para culturas de cobertura, plantio direto e a maioria dos tipos de agrofloresta, isso é de 8 a 12 anos a partir da linha de base. Ciclos de relato de três anos ainda são úteis, desde que relatem em relação à trajetória da MDD, e não a um sinal presumido.

O perfil profundo já não é opcional. Os ganhos no subsolo em sistemas perenes com mais de ~18 anos muitas vezes superam os ganhos do solo superficial e são invisíveis para um desenho de 0–30 cm5721. A correção MSE é igualmente inegociável em comparações entre preparo convencional e plantio direto67.

O protocolo contra o qual o comprador é auditado define o piso. O VM0042 v2.2/v3.0, o CRCF da UE (Regulamento 2024/3012), o GHG Protocol Land Sector and Removals Standard e o SBTi FLAG convergem para exigências mais rígidas de dados empíricos ao longo de 2026–202713182223. Desenhar para o padrão mais permissivo com o qual a equipe está familiarizada, quando o comprador será auditado contra o mais rígido, é um modo de falha em câmera lenta.

Pontos-chave