Como Modelar a Mudança do Carbono do Solo ao Longo do Tempo?

Um modelo de carbono do solo parte de uma descrição de um sistema (clima, solo, manejo) e retorna uma trajetória: como se espera que o estoque de carbono evolua sob aquela descrição. É por meio dessa trajetória que as organizações justificam a alocação de capital em programas regenerativos, projetam emissões e remoções de Escopo 3 sob a GHG Protocol Land Sector and Removals Guidance e, quando a estrutura do projeto o justifica, apresentam reivindicações nos mercados voluntários de carbono.

A pergunta “qual modelo devemos usar?” não tem resposta universal. A pergunta certa é mais estrita: qual modelo é adequado à decisão, para o sistema em questão, dados os dados que o projeto consegue obter de forma realista? Respondê-la é a decisão de julgamento mais importante em um trabalho de modelagem de carbono do solo.

Por baixo das siglas, os modelos de uso corrente compartilham um mesmo mecanismo: o carbono do solo repartido em compartimentos que decaem, cada um a seu próprio ritmo, alimentados por um fluxo contínuo de aportes frescos de plantas e de esterco. Compreendido esse mecanismo, a maioria das questões práticas se resolve por si, inclusive as duas que silenciosamente decidem uma previsão de décadas: onde se supõe que os compartimentos começam e quão ampla é a resposta honesta. As três ferramentas deste artigo executam uma versão compacta desse mecanismo, de modo que o comportamento passa a ser algo que se pode manipular, e não apenas algo sobre o que se lê.

O que a modelagem de fato entrega, e o que não entrega

Um modelo de carbono do solo é uma hipótese estruturada sobre como o carbono se move pelo sistema. Seus resultados são condicionais: dados estes dados de entrada, eis a trajetória esperada e sua incerteza. Um modelo não mede o solo; ele o prevê.

Isso importa porque boa parte da divergência entre desenvolvedores de projetos e verificadores nasce de confundir os dois1. Uma trajetória modelada usada no lugar de dados de campo no relato anual é um risco de conformidade. Uma trajetória modelada usada ao lado de dados de campo, com as faixas de incerteza explícitas, é prática consagrada sob a GHG Protocol Land Sector and Removals Guidance e sob a maioria das metodologias dos mercados voluntários.

Como esses modelos de fato funcionam

O carbono do solo não é uma única substância a decair a uma única taxa. Uma raiz fresca e um complexo húmico de séculos são ambos “carbono orgânico do solo”, mas um se renova em meses e o outro sobrevive a quem o mede. Todo modelo de uso corrente trata isso da mesma maneira. Reparte o carbono num punhado de compartimentos, cada qual com seu próprio tempo de renovação, e os acompanha separadamente2.

Dentro de cada compartimento, a decomposição segue uma cinética de primeira ordem: a quantidade perdida em um dado mês é proporcional à quantidade presente, de modo que um compartimento deixado a si mesmo decai ao longo de uma curva exponencial fixada por uma constante de taxa k.

A constante de taxa k é uma propriedade do compartimento: elevada para os resíduos frescos, baixa para a matéria humificada, praticamente nula para a fração inerte. A função f ajusta k para cima ou para baixo conforme as condições em que o solo de fato se encontra. O calor acelera a decomposição, o estresse hídrico a retarda e um dossel vivo a freia em relação ao solo descoberto. É por isso que o mesmo aporte de resíduos acumula carbono sob uma planta de cobertura numa estação temperada fria e o perde numa estação quente e seca.

Nem todo o carbono decomposto deixa o solo. Uma fração é respirada para a atmosfera como CO₂; o restante é estabilizado em biomassa microbiana e em compostos humificados que realimentam os compartimentos lentos. O teor de argila governa essa partição, porque solos de textura mais fina protegem uma parcela maior do carbono processado. É por isso que a fração de argila é um dado de entrada exigido em todos os modelos apresentados a seguir3.

Junte essas peças e um solo sob manejo estável deriva rumo a um estado estacionário, o estoque no qual os aportes anuais equilibram exatamente as perdas anuais. Aumente os aportes e o estado estacionário sobe; aqueça o clima ou reduza a argila e ele cai. A aproximação é lenta, cadenciada pelo tempo de renovação do compartimento ativo mais lento, razão pela qual ganhos críveis de carbono do solo são expressos ao longo de décadas, e não de anos. A ferramenta abaixo executa esse mecanismo para a frente. Ajuste o aporte, a argila e a temperatura e observe onde o estoque se estabiliza e quanto tempo leva para chegar lá.

Duas características dessa curva importam para um projeto. O destino é fixado pelos aportes e pelo local, não pelo ponto de onde o solo por acaso parte, de modo que dois talhões conduzidos sob o mesmo manejo convergem para o mesmo estado estacionário a partir de direções diferentes. E a jornada é longa: mesmo uma grande mudança precisa de décadas para se manifestar, e os primeiros anos avançam devagar. Uma trajetória modelada que alcança um novo equilíbrio em cinco anos está descrevendo um artefato de sua configuração, não um solo.

Os quatro modelos candidatos

Quatro arcabouços dominam o trabalho prático. Cada um foi concebido para uma classe específica de pergunta e cada um carrega a marca do sistema sobre o qual foi construído. O RothC surgiu dos experimentos de longo prazo de Rothamsted em camada superficial arável de clima temperado2. O AMG foi calibrado inicialmente em ensaios de cultivo franceses e, em sua forma atual, reconstruído em mais de quarenta locais europeus de longo prazo45. O Century foi concebido para explicar os níveis de matéria orgânica em pastagens das Grandes Planícies, com nitrogênio, fósforo e enxofre acoplados ao ciclo do carbono6, e o DayCent é seu irmão de passo de tempo diário, construído para resolver os fluxos de gases-traço que um modelo mensal não consegue7. Usar um deles fora de seu domínio de origem é possível, mas caro em credibilidade. Sobre como testar se um novo local está sequer dentro desse domínio, veja Representativo de quê?

Os quatro foram aferidos contra dados experimentais de longo prazo em dezenas de locais. A comparação clássica confrontou nove modelos de matéria orgânica do solo com sete conjuntos de dados de longo prazo e não encontrou um vencedor único: o desempenho dependia do local e da pergunta8. Trabalhos de conjunto (ensemble) mais recentes chegam ao mesmo veredito e acrescentam um alerta: a dispersão entre modelos é, ela própria, uma fonte de incerteza que uma projeção de modelo único esconde11. RothC e Century têm o mais longo histórico internacional; o AMG é a opção mais forte para sistemas de cultivo franceses e demais europeus; o DayCent é a ferramenta especializada quando as emissões de N₂O estão no escopo.

Escolhendo o modelo certo para a pergunta

Quatro perguntas práticas determinam a escolha.

O nitrogênio está no escopo?

Se a decisão depende das emissões de N₂O, da eficiência de uso do fertilizante ou da interação entre o manejo do nitrogênio e a dinâmica do carbono, o RothC é estruturalmente a escolha errada; ele não tem compartimento de N. Century ou DayCent são as respostas honestas. Se a decisão é puramente sobre o carbono do solo sob manejo estável de nitrogênio, o RothC é defensável e mais simples.

O sistema é limitado por água?

O modificador hídrico padrão do RothC tem desempenho fraco em sistemas de terras secas e mediterrâneos, onde a atividade biológica fica limitada pela umidade durante boa parte do ano9. Existe uma calibração modificada do RothC (Farina et al.) que é a escolha prática nesses contextos. Century e DayCent tratam a limitação hídrica de forma nativa por meio de seus submodelos.

Qual resolução temporal é necessária?

O passo mensal é suficiente para uma trajetória de décadas ao longo de 20 anos. O diário é necessário quando o resultado de interesse é um fluxo rápido: emissões de gases-traço após uma aplicação de fertilizante, por exemplo. RothC, AMG e Century operam em passo mensal; o DayCent opera em passo diário.

Que dados o projeto consegue produzir de forma realista?

Um modelo que precisa de dados de entrada que o projeto não consegue fornecer é um modelo que rodará sobre suposições. Se dados meteorológicos diários não estiverem disponíveis, o DayCent operará sobre uma redução de escala (downscaling), o que é defensável se essa redução de escala for declarada. Se os aportes de resíduos e o histórico de manejo forem desconhecidos, todo modelo está na mesma posição: a incerteza dos dados de entrada faltantes domina a incerteza do modelo.

Calibração, validação e saber quando um modelo se ajusta

Um modelo que nunca foi testado contra carbono medido em um sistema comparável é uma hipótese, não evidência. A calibração ajusta os parâmetros de um modelo para reproduzir um conjunto de dados conhecido; a validação faz a pergunta mais difícil: se ele então prevê dados que nunca viu. As duas são rotineiramente confundidas, e um modelo que se ajusta primorosamente aos seus dados de calibração ainda pode falhar em um local independente.

Os instrumentos para isso são os experimentos de longo prazo, alguns em curso há mais de um século, nos quais o carbono do solo foi acompanhado sob tratamentos fixos. A habilidade de um modelo é reportada em relação a eles com estatísticas conhecidas: a raiz do erro quadrático médio para o tamanho do desvio e a eficiência do modelo para saber se ele supera a simples predição da média. A comparação de nove modelos que estabeleceu o padrão para esse trabalho encontrou concordância em alguns locais e divergência acentuada em outros, sem nenhum arcabouço dominante em toda parte8. Um conjunto (ensemble) de 26 modelos em solos de pousio nu chegou a uma conclusão ainda mais contundente: ao longo de projeções longas, a escolha do modelo pode importar tanto quanto a escolha do manejo11.

Duas regras práticas decorrem disso. Valide com dados de séries temporais da região e do sistema que se está modelando, não com a média global, porque um modelo calibrado em solos aráveis de clima temperado não oferece garantia alguma para a agrofloresta tropical. E trate o domínio de aplicabilidade como uma fronteira real: um modelo instado a extrapolar para além das condições em que foi testado devolve um número com falsa confiança, que é o tema de Representativo de quê? Quando um projeto precisa de contabilidade baseada em modelos no nível que os inventários nacionais usam, o IPCC chama isso de Tier 3 e exige exatamente essa evidência: um modelo validado contra medições representativas do sistema12. As revisões que acompanham o rumo da área fazem a mesma observação a partir do lado da pesquisa10.

Inicialização: a fonte silenciosa de erro

Todo modelo de carbono do solo reparte o estoque em compartimentos com diferentes taxas de renovação: rápido, lento e quase inerte. A divisão inicial entre os compartimentos raramente é mensurável e costuma ser estimada. Para o RothC, a equação de Falloon estima o compartimento inerte a partir do COS total; é uma aproximação estatística com efeito material sobre projeções de 30 anos13.

A equação acima é um ajuste estatístico, não uma medição física. Seus parâmetros carregam uma incerteza que se acumula ao longo de projeções de décadas. A alternativa mais defensável é o spin-up: rodar o modelo até o equilíbrio sob um regime histórico de manejo reconstruído e então usar a distribuição de compartimentos resultante como ponto de partida. A reconstrução é a parte difícil, e também é a etapa mais comumente negligenciada na documentação de projetos. Uma retroprojeção histórica de 50 anos sob uma rotação realista reduz o efeito da inicialização sobre as trajetórias previstas em uma ordem de grandeza, em comparação com a inicialização por compartimentos padrão14.

A ferramenta abaixo mostra por que isso não é um detalhe técnico. Fixe um estoque medido, fixe o manejo futuro e mude apenas como o compartimento inerte é suposto. A fração inerte nunca ganha nem perde carbono, de modo que colocar mais do estoque ali deixa menos carbono capaz de responder. A mesma medição, repartida de duas maneiras, pode produzir um ganho projetado sob uma suposição e uma perda projetada sob a outra. Essa é a decisão inteira, decidida por um número que ninguém mediu.

Ler os resultados como envelopes, não como linhas

Uma trajetória de carbono do solo com um intervalo de confiança de 95 por cento estreito reflete ou um sistema bem calibrado e rico em dados, ou subnotificação. A propagação de Monte Carlo da incerteza dos dados de entrada (variabilidade climática, teor de argila, aportes de resíduos, inicialização) produz um envelope de trajetória que costuma ser bem mais largo que a linha média. Esse envelope é a entrega honesta. Retire-o e o número que resta é, na melhor das hipóteses, uma estimativa central plausível.

Para o trabalho de comparação de cenários (linha de base versus intervenção, por exemplo), o resultado útil é a diferença entre cenários, com sua própria incerteza propagada. Boa parte do ruído dos dados de entrada é compartilhada entre os cenários e se cancela na diferença, que muitas vezes é mais estreita que qualquer das trajetórias individuais. Esse é o enquadramento certo para quase todas as questões de viabilidade de projetos de carbono.

A ferramenta abaixo torna isso concreto. Alterne entre mostrar os dois cenários separadamente, cada um com um envelope largo o suficiente para engolir o sinal, e mostrar a diferença entre eles, na qual o ruído compartilhado de clima e de solo se cancela e um número defensável emerge. Quanto mais a linha de base e a intervenção compartilham as mesmas condições, mais estreita essa diferença se torna. Essa é a razão estatística pela qual arcabouços de MRV sérios recompensam desenhos pareados e locais de controle em vez de reivindicações absolutas de antes e depois15.

A fronteira: o que os modelos de compartimentos deixam de fora

Os quatro modelos de uso corrente descendem todos da mesma ideia: compartimentos conceituais definidos pela taxa de renovação. É uma abstração poderosa, mas os compartimentos não são coisas que se possa ir a campo e medir. Não se peneira um solo em suas frações decomponível e humificada; são compartimentos contábeis inferidos a partir do comportamento. Há vinte anos essa lacuna impulsiona uma segunda geração de modelos construídos sobre grandezas que se podem medir.

A mudança na ciência subjacente vai da química para a biologia e a mineralogia. A persistência, na visão atual, tem menos a ver com uma molécula ser intrinsecamente difícil de decompor e mais com se os micróbios conseguem alcançá-la e se as superfícies minerais a retêm18. Modelos microbianamente explícitos como MIMICS e Millennial colocam a biomassa microbiana e o carbono associado a minerais no centro e substituem os compartimentos inferidos por frações que podem, em princípio, ser medidas em laboratório1617. Eles respondem de forma diferente dos modelos de primeira ordem sob aquecimento, e a diferença não é pequena.

Nada disso aposenta o RothC ou o Century para o trabalho de projetos. Os modelos mais novos são mais difíceis de parametrizar, menos aferidos contra dados de campo de longo prazo e ainda não incorporados às normas de MRV. Mas eles assinalam onde a incerteza de uma projeção de décadas de fato reside, e um profissional que trate a estrutura de compartimentos como verdade física, e não como uma ficção útil, será surpreendido pela próxima década da literatura.

Da trajetória à reivindicação defensável

Uma trajetória só se torna um ativo quando uma norma a aceita, e os arcabouços que importam aqui convergem para uma mesma postura: modelar para planejar e para interpolar entre medições, mas medir para reivindicar. A GHG Protocol Land Sector and Removals Guidance trata as remoções modeladas como passíveis de relato apenas quando acompanhadas de um programa de medição e monitoramento e com a incerteza declarada1. Os inventários nacionais que usam contabilidade Tier 3 baseada em modelos carregam a mesma exigência de validação contra dados representativos12. O artigo de referência sobre MRV de carbono do solo é explícito ao afirmar que os modelos reduzem o custo de amostragem, mas não eliminam a necessidade de medir15.

Assim, a entrega honesta de um trabalho de modelagem raramente é um número único. É uma trajetória com seu envelope, uma declaração clara de qual modelo foi usado e por quê, uma inicialização documentada e um plano para as medições que a confirmarão ou corrigirão. Estudos em escala continental modelam milhões de hectares exatamente com essa lógica19, e o mesmo vale na escala de talhão: o modelo conduz a projeção entre os pontos onde o solo é de fato amostrado.

Pontos-chave