Global warming increases lightning activity over Brazil regardless of the trend in cloud ice content simulated by a GCM

Resumo Técnico: O Aquecimento Global Aumenta a Atividade de Raios sobre o Brasil, Independentemente da Tendência no Conteúdo de Gelo nas Nuvens Simulada por um GCM

Problema
A evolução futura da atividade de raios sob o aquecimento global permanece uma fonte significativa de incerteza na ciência climática. Enquanto muitos estudos projetam um aumento nos raios com base em proxies termodinâmicos como a Energia Potencial Convectiva Disponível (CAPE) ou a altura do topo das nuvens, outros sugerem uma potencial diminuição nas regiões tropicais quando os raios são parametrizados usando o fluxo de gelo nas nuvens. Essa discrepância surge de inconsistências fundamentais na forma como os Modelos de Circulação Global (GCMs) simulam a microfísica das nuvens, particularmente a redução no conteúdo de gelo simulado nas nuvens sob cenários de aquecimento. Além disso, modelos de regressão empírica enfrentaram limitações devido a vieses em conjuntos de dados observacionais (por exemplo, eficiência de detecção de satélite sobre a América do Sul) e suposições estatísticas quanto à independência de dados temporais. Existe uma escassez específica de projeções robustas de raios para a América do Sul, com estudos anteriores limitados a estados brasileiros específicos. Este estudo aborda a necessidade de quantificar como as mudanças projetadas na energia convectiva e no conteúdo de gelo nas nuvens influenciam a densidade de raios em todo o Brasil, investigando especificamente se as mudanças termodinâmicas podem compensar as reduções no conteúdo de gelo das nuvens.

Metodologia
O estudo utiliza uma abordagem de conjunto de dados de múltiplas fontes para projetar a densidade de raios sobre o Brasil para o período 2071–2100 sob cenários intermediários (SSP245) e de altas emissões (SSP585).

• Fontes de Dados:
  • Raios: Observações da World Wide Lightning Location Network (WWLLN) para o período 2013–2025, corrigidas pela Eficiência de Detecção Relativa (RDE).
  • Reanálise Atmosférica: Dados do ERA5 (1979–presente) para treinamento e validação do modelo.
  • Projeções Climáticas: Saídas do Modelo de Circulação Global (GCM) EC-Earth3-Veg do CMIP6, selecionado por sua consistência física com o ERA5.
• Seleção de Preditores: Um conjunto de variáveis atmosféricas foi triado com base em correlações espaciais, habilidade em reproduzir o ciclo sazonal e significância estatística. Variáveis termodinâmicas absolutas (ex: umidade específica) foram excluídas para evitar vieses associados a mudanças termodinâmicas uniformes. Os preditores finais incluíram CAPE, conteúdo total de gelo na coluna (ICE), taxa de precipitação convectiva (CPR), pressão ao nível do mar (MSLP) e profundidade de fase mista (MPD).
• Abordagem de Modelagem: Dois modelos de Regressão Linear Múltipla (MLR) foram desenvolvidos para projetar a densidade de raios:
  1. Modelo 1: Inclui explicitamente o conteúdo de gelo nas nuvens (ICE) junto com CAPE e MSLP.
  2. Modelo 2: Exclui o ICE, utilizando CPR e MPD como proxies para a convecção realizada e o potencial para formação de gelo, respectivamente, junto com CAPE e MSLP.
• Correção de Viés: Um método de escalonamento baseado em distribuição foi aplicado para alinhar as simulações do EC-Earth3-Veg com a climatologia do ERA5 antes de projetar tendências futuras.
• Análise: O estudo emprega testes não paramétricos de Mann-Kendall e estimadores de inclinação de Sen para análise de tendência. A sensibilidade dos raios é quantificada como uma mudança fracionária por Kelvin de aquecimento global (% K⁻¹).

Principais Contribuições

• Primeira Projeção em Escala Nacional: Este estudo fornece a primeira projeção abrangente de raios para todo o território brasileiro, cobrindo diversos regimes climáticos (monçônico tropical, equatorial e subtropical).
• Estrutura de Modelo Duplo: Ao desenvolver dois modelos distintos — um explicitamente dependente de gelo nas nuvens e outro dependente de proxies de convecção realizada — o estudo isola o impacto dos vieses do GCM na simulação de gelo nas nuvens nas projeções de raios.
• Resolução do Debate Gelo vs. Instabilidade: A pesquisa aborda diretamente o conflito entre a redução projetada de gelo nas nuvens e o aumento da instabilidade convectiva, oferecendo um arcabouço estatístico para avaliar seus efeitos combinados.
• Caracterização de Regimes Regionais: O estudo identifica três regimes atmosféricos distintos que controlam os raios no Brasil: um regime controlado por termodinâmica no Brasil central, um regime de saturação de umidade no extremo norte da Amazônia e um regime impulsionado por sinótica no sul do Brasil.

Resultos

• Tendências Atmosféricas: O modelo EC-Earth3-Veg projeta um aumento generalizado de CAPE (3–9% por década) em todo o Brasil, indicando maior potencial de energia convectiva. Inversamente, o conteúdo total de gelo nas nuvens (ICE) mostra uma redução generalizada no Brasil tropical, consistente com as tendências globais de GCM, embora aumente no sul do Brasil. A precipitação convectiva (CPR) geralmente diminui nos trópicos, sugerindo uma redução na frequência de eventos convectivos, apesar do potencial de maior intensidade.
• Projeções de Raios: Ambos os modelos projetam um aumento espacialmente disseminado na densidade de raios sobre o Brasil até 2071–2100, apesar da redução simulada no conteúdo médio de gelo nas nuvens.
  • Sensibilidade Térmica: O aumento projetado varia de 8% a 11,5% por Kelvin de aquecimento.
  • Mudanças Absolutas: Os aumentos absolutos na densidade de raios variam de 14% a 51%, dependendo do cenário e do conjunto de preditores.
  • Comparação de Modelos: O Modelo 2 (excluindo ICE) projeta sensibilidades mais altas (10,6–11,5% K⁻¹) em comparação ao Modelo 1 (7,6–8,1% K⁻¹), mas ambos concordam com a direção da tendência.
  • Distribuição Espacial: Os maiores aumentos são projetados para a Amazônia brasileira (atingindo 9–24% K⁻¹), o Nordeste e o sul/sudeste do Brasil. Reduções localizadas estão confinadas principalmente ao estado do Acre, no extremo oeste da Amazônia.
• Mecanismo: O estudo sugere que o aumento na intensidade dos raios compensa a redução na frequência de eventos convectivos. Valores mais altos de CAPE podem intensificar as correntes ascendentes e as taxas de colisão nas tempestades que ocorrem, mitigando parcialmente os efeitos da redução da disponibilidade média de gelo. Alternativamente, a redução do gelo médio nas nuvens pode refletir uma diminuição no número de tempestades, em vez do conteúdo de gelo dentro de tempestades individuais, que podem se tornar mais severas.

Significância
O artigo afirma que suas descobertas são significativas para entender os riscos associados às mudanças climáticas na Amazônia e em outros biomas brasileiros. O aumento projetado na atividade de raios, mesmo diante de uma redução simulada de gelo nas nuvens, sugere um risco elevado de incêndios florestais iniciados por raios em florestas que já estão se tornando mais suscetíveis ao fogo devido ao desmatamento e à seca. Além disso, os resultados têm implicações para o ciclo de carbono global; uma vez que os raios são um importante motor natural de mortalidade de árvores de grande porte (que armazenam carbono significativo), o aumento da frequência de raios poderia acelerar o declínio do sumidouro de carbono da Amazônia. O estudo conclui que o aumento da atividade de raios sobre a Amazônia brasileira é robusto à representação específica de gelo nas nuvens nos GCMs, destacando um ciclo de feedback crítico entre o aquecimento global, a eletrificação atmosférica e a resiliência florestal.