How Uncertain Are Estimates of Marine Carbon Dioxide Removal? Insights from a Simplified Model

Resumo Técnico: Quão incertas são as estimativas de Remoção de Dióxido de Carbono Marinho? Insights de um modelo simplificado

Definição do Problema
As técnicas de Remoção de Dióxido de Carbono Marinho (mCDR), especificamente o Aumento da Alcalinidade Oceânica (OAE) e a Remoção Direta do Oceano (DOR), dependem da redução da pressão parcial de CO2 na superfície do oceano para absorver o CO2 atmosférico. Embora a eficiência máxima teórica dessas intervenções seja conhecida (ex: ~0,8 mol de CO2 por mol de alcalinidade para OAE), a evolução temporal da remoção de carbono é altamente sensível aos processos físicos e biogeoquímicos. As estimativas atuais da eficácia da mCDR dependem fortemente de modelos oceânicos, porém as incertezas inerentes a esses modelos — particularmente no que diz respeito à mistura vertical, troca gasosa ar-mar e propriedades das massas de água — são mal compreendidas. Os modelos globais existentes frequentemente empregam parametrizações fixas que podem subestimar a variabilidade natural, e sua complexidade torna difícil isolar os mecanismos específicos que impulsionam a incerteza nos cálculos de CDR. Essa falta de quantificação da incerteza impõe um desafio fundamental para as estruturas de Monitoramento, Relato e Verificação (MRV) e para o tempo de neutralidade líquida de carbono para projetos de mCDR.

Metodologia
Para examinar sistematicamente essas incertezas, o estudo utiliza um modelo oceânico abiótico unidimensional (1D). Esta abordagem de complexidade reduzida isola os processos verticais (difusão turbulenta e troca gasosa ar-mar) enquanto ignora a advecção horizontal, subducção e interações biológicas. O modelo evolui a Alcalinidade Total (TA) e o Carbono Inorgânico Dissolvido (DIC) ao longo do tempo, calculando o fluxo de CO2 superficial com base na troca gasosa ar-mar.

O estudo utiliza duas estruturas analíticas principais:

1. Simulações de Sensibilidade: Parâmetros individuais do modelo (ex: profundidade da camada de mistura, velocidade do vento, coeficiente de transferência de gás, difusividade da picnoclina, TA/DIC de fundo, profundidade de implantação) são perturbados por valores consistentes com a variabilidade natural observada. Estes são comparados contra uma simulação de referência para determinar a "razão de sensibilidade de CDR" ($r_{CDR}$), definida como a razão entre o CDR em uma execução de sensibilidade e a execução de referência.
2. Ensembles de Monte Carlo: Mais de 1.000 pares de simulações (referência e intervenção) são executados simultaneamente, com parâmetros extraídos de distribuições estatísticas que representam a variabilidade oceânica observada. Esta abordagem avalia como as incertezas combinadas nos processos físicos e nas propriedades das massas de água se manifestam nas estimativas de CDR.

As simulações são conduzidas para profundidades oceânicas de 200 m e 2.000 m ao longo de um período de 50 anos para avaliar o impacto da profundidade na persistência da incerteza.

Principais Resultados
O estudo identifica três regimes temporais distintos que regem a incerteza do CDR:

1. Equilíbrio Rápido Inicial ($t \ll \tau_{mix}$): Nos primeiros anos, as estimens de CDR são altamente sensíveis aos parâmetros que controlam a escala de tempo inicial da troca ar-mar ($\tau_{air-sea}$), incluindo a profundidade da camada de mistura, velocidade do vento, coeficientes de transferência de gás e a química de carbonatos de fundo (TA, DIC, temperatura). Durante esta fase, as estimativas de CDR podem variar por um fator de dois dependendo dos vieses dos parâmetros.
2. Perda de Alcalinidade da Camada de Mistura ($t \approx \tau_{mix}$): À medida que a mistura vertical remove a alcalinidade da superfície, as incertezas passam a ser dominadas pela interação entre a escala de tempo de troca ar-mar e a escala de tempo de mistura ($\tau_{mix}$). As incertezas na difusividade da picnoclina ($\kappa_{pyc}$) tornam-se críticas, fazendo com que as estimativas de CDR alcancem o pico de variabilidade (até 30%) 2 a 3 anos após a implantação.
3. Mistura Profunda ($t \gg \tau_{mix}$): A longo prazo, as incertezas persistem por décadas se a escala de tempo para a mistura no oceano profundo for comparável ou mais rápida que o equilíbrio ar-mar. O estudo descobre que a incerteza na difusividade da picnoclina é o fator dominante para a incerteza de CDR a longo prazo, pois controla tanto a fração de alcalinidade equilibrada antes da subducção profunda quanto a taxa de retorno das águas profundas não equilibradas à superfície.

Achados Quantitativos:

• Magnitude da Incerteza: Nos primeiros anos, as estimativas de CDR podem diferir em até 100% (um fator de dois) da simulação de referência. Após 20 a 50 anos, essa diferença diminui para 10–20%, embora uma incerteza significativa permaneça para implantações em águas profundas.
• Profundidade de Implantação: Implantar alcalinidade abaixo da camada de mistura (mesmo que por 20–40 m) resulta em eficiência significativamente menor e maior incerteza, com algumas estimativas permanecendo muito abaixo da simulação de referência mesmo após 50 anos.
• Tempo de Neutralidade Líquida: O tempo necessário para uma intervenção tornar-se líquida de carbono negativa (considerando os custos de Avaliação de Ciclo de Vida) é altamente variável. Enquanto um caso de referência sugere ~1,4 anos, o ensemble de Monte Carlo mostra uma ampla distribuição: ~30% das simulações atingem a neutralidade líquida em menos de um ano, 50% levam mais de 2 anos e ~4% levam mais de 50 anos. Esta "cauda longa" está primariamente associada a implantações profundas ou mistura rápida que impede o equilíbrio precoce.
• Dependência da Profundidade: As incertezas persistem por mais tempo em oceanos de 2.000 m de profundidade em comparação com oceanos de 200 m, onde toda a coluna de água pode se equilibrar dentro do intervalo de tempo da simulação.

Significância e Alegações
O artigo argumenta que as incertezas nas estimativas de mCDR não são meramente transitórias, mas podem persistir por décadas devido às escalas de tempo da mistura do interior em relação à troca ar-sea. Os autores afirmam que estes resultados identificam um desafio fundamental para as simulações de mCDR regionais e globais, onde os processos físicos governantes são frequentemente mal restringidos.

Especificamente, o estudo destaca que:

• Implicações para MRV: A incerteza na evolução temporal da remoção de carbono impacta diretamente o tempo de neutralidade líquida, o que é crítico para as estruturas de crédito de carbono. Sem incertezas bem quantificadas, a concessão de créditos de carbono pode ser prematura ou imprecisa.
• Limitações de Modelos: Os modelos oceânicos atuais frequentemente dependem de parametrizações para turbulência e troca de gases que estão sujeitas a subtanciais vieses espaciais e temporais. O estudo sugere que validar estes processos físicos específicos (mistura e transferência de gás) é mais crítico para a mCDR do que para projeções climáticas gerais.
• Estratégia de Implantação: Os resultados enfatizam a necessidade de garantir que a alcalinidade seja implantada estritamente dentro da camada de mistura superficial para maximizar a eficiência e minimizar a incerteza.

Os autores mantêm um escopo modesto, observando que seu modelo 1D exclui a advecção horizontal e os feedbacks biológicos, que poderiam introduzir outras incertezas. No entanto, eles concluem que, mesmo em um sistema abiótico, as incertezas decorrentes de processos físicos sozinhas são significativas o suficiente para justificar uma investigação séria adicional, a fim de garantir que as intervenções de mCDR sejam genuinamente negativas em carbono em escalas de tempo decenais.