Divergent Runoff Responses to 1.5°C, 2.0°C, and 3.0°C Warming Across Headwater and Downstream Regions of the Indus River Basin

Resumo Técnico: Respostas Divergentes de Escoamento para Aquecimento de 1,5°C, 2,0°C e 3,0°C entre Regiões de Cabeceira e Jusante da Bacia do Rio Indo

Problema
As bacias hidrográficas montanhosas, particularmente na Ásia Central de Alta Montanha, são altamente sensíveis ao aquecimento climático, que altera a magnitude, o tempo e a fase do escoamento. Embora estudos anteriores tenham estabelecido a sensibilidade da Bacia do Alto Indo (UIB) — dominada por neve e geleiras — ao aquecimento, há uma quantificação insuficiente de como esses sinais se propagam para a Bacia do Baixo Indo (LIB), situada jusante. A LIB é um sistema fortemente gerido, dependente de irrigação, onde a dependência de águas subterrâneas e a regulação humana modificam significativamente o fluxo natural. Existe uma lacuna de conhecimento crítica sobre se o aquecimento global progressivo produz uma divergência hidrológica onde os ganhos de escoamento a montante coincidem com déficits a jusante. Além disso, muitas avaliações baseiam-se em recortes temporais fixos em vez de Níveis de Aquecimento Global (GWLs) equivalentes, o que falha em considerar as diferentes taxas de aquecimento entre os modelos e as trajetórias de emissão. Este estudo aborda essas lacunas avaliando as respostas de escoamento em GWLs de 1,5°C, 2,0°C e 3,0°C em todo o sistema de bacia superior-inferior.

Metodologia
O estudo emprega uma estrutura de escala de estação, integrando dados de 35 estações de medição hidrológica na Bacia do Rio Indo (IRB).

• Forçante de Dados: A análise utiliza dados de forçante hidroclimática com correção de viés provenientes do Inter-Sectoral Impact Model Intercomparison Project 3b (ISIMIP3b) e do Coupled Model Intercomparison Project Phase 6 (CMIP6). Quatro Modelos de Circulação Geral (GCMs: CNRM-CM6-1, CNRM-ESM2-1, EC-EARTH3, MIROC6) são utilizados sob três Trajetórias Socioeconômicas Compartilhadas (SSP1-2.6, SSP3-7.0, SSP5-8.5).
• Estrutura de Modelagem: Um Temporal Fusion Transformer (TFT), uma arquitetura de aprendizagem profunda projetada para previsão de múltiplos horizontes, é treinado em registros mensais de vazão observados e nos preditores climáticos com correção de viés. O modelo incorpora 11 preditores dinâmicos (ex: precipitação, temperatura, radiação, umidade) e covariáveis estáticas (latitude, longitude) para capturar controles hidroclimáticos não lineares.
• Extração de GWL: Em vez de anos fixos, as respostas hidrológicas são extraídas com base no alcance de GWLs específicos (1,5°C, 2,0°C, 3,0°C) em relação à linha de base de 1850–1900. O ano de alcance é definido como o ponto médio da primeira janela de 20 anos que excede o limiar.
• Diagnósticos: O estudo quantifica anomalias anuais e sazonais de escoamento, fração de neve, Meses de Temperatura Positiva (PDM), sensibilidade ao aquecimento (m³/s por °C) e incerteza (intervalo interquartil, concordância entre modelos).

Principais Resultados
A análise revela uma divergência hidroclimática robusta e longitudinalmente diferenciada entre a UIB e a LIB sob aquecimento progressivo:

• Divergência de Escoamento:
  • Bacia do Alto Indo (UIB): O escoamento anual mediana aumenta monotonicamente com o aquecimento, subindo de +172,80 m³/s a 1,5°C para +725,49 m³/s a 3,0°C. A sensibilidade ao aquecimento é de +289,40 m³/s por °C. Este aumento é acompanhado por um alargamento da faixa de incerteza (o IQR aumenta de 46,35 para 137,02 m³/s).
  • Bacia do Baixo Indo (LIB): Em contraste, os déficits de escoamento anual mediana aprofundam-se, passando de −78,40 m³/s a 1,5°C para −438,86 m³/s a 3,0°C, com uma sensibilidade de −187,20 m³/s por °C.
• Reorganização Sazonal:
  • UIB: Os ganhos concentram-se nas estações quentes (JJA e SON), impulsionados pelo degelo intensificado, embora a magnitude das anomalias varie conforme o GWL.
  • LIB: As perdas persistem de MAM a SON, com apenas mudanças positivas modestas no inverno (DJF). O regime sazonal na LIB torna-se estruturalmente enfraquecido e cada vez mais negativo.
• Direcionadores Criosféricos:
  • Na UIB, o aumento do escoamento coincide com um declínio sistemático na fração de neve (de −0,03 a 1,5°C para −0,06 a 3,0°C) e um aumento acentuado no PDM (de +9,40 para +31,60). Isso indica uma transição para regimes dominados pela chuva e um aumento da energia térmica de degelo.
  • A LIB mostra aumentos de temperatura semelhantes, mas carece da amplificação criosférica vista a montante, refletindo, em vez disso, os efeitos combinados da evaporação e da regulação humana.
• Heterogeneidade Espacial:
  • A análise em escala de estação confirma que os sinais em escala de bacia emergem de respostas locais heterogêneas. Estações da UIB (ex: Gilgit, Shigar) mostram forte sensibilidade positiva, enquanto estações da LIB (ex: Guddu, Kotri) mostram forte sensibilidade negativa. Algumas estações intermediárias apresentam sinais fracos ou não significativos, destacando a incerteza em escala local.

Significância e Alegações
O artigo afirma que a Bacia do Rio Indo não responde às mudanças climáticas como um sistema espacialmente uniforme. Em vez disso, o aquecimento impulsiona uma reorganização longitudinal pronunciada onde o aumento do escoamento a montante não se traduz em disponibilidade de água a jusante.

• Divisão Hidroclimática: O estudo demonstra que a IRB funciona como um sistema acoplado, mas funcionalmente heterogêneo. A UIB atua como um sistema condicionado pela criosfera e movido por energia, onde o aquecimento amplifica temporariamente o escoamento via degelo. Por outro lado, a LIB atua como um sistema altamente regulado, onde a vazão a jusante é desacoplada dos ganhos a montante devido às retiradas para irrigação, desvios de canais e alta demanda evaporativa.
• Natureza Transitória dos Ganhos: Os autores enfatizam que o aumento projetado no escoamento na UIB deve ser interpretado como uma fase transitória de "pico de água" (peak water) associada ao degelo intensificado sob a redução do armazenamento criosférico, e não como um ganho permanente de recursos hídricos.
• Relevância para Políticas: Os resultados sublinham que o planejamento de adaptação não pode basear-se apenas nas projeções de escoamento a montante. A divergência apresenta riscos específicos para a confiabilidade da irrigação, dependência de águas subterrâneas e segurança alimentar na LIB. O estudo fornece um quadro fisicamente fundamentado para avaliar riscos em limiares de aquecimento específicos (1,5°C, 2,0°C, 3,0°C), oferecendo percepções relevantes para a tomada de decisão no gerenciamento de um dos sistemas fluviais mais estressados do mundo.

Limitações
Os autores observam que o framework TFT captura a resposta líquida do escoamento sem separar explicitamente os processos hidrológicos (ex: partição chuva-escoamento vs. degelo de geleiras). Além disso, embora os GWLs melhorem a comparabilidade entre modelos, permanecem incertezas relativas à seleção de GCMs e correção de viés. O estudo não modela explicitamente a dinâmica do balanço de massa das geleiras ou esquemas de operação de reservatórios, sugerindo a necessidade de integração futura com modelos baseados em processos físicos.