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

Resumen Técnico: Respuestas de Escurrimiento Divergentes ante el Calentamiento de 1.5°C, 2.0°C y 3.0°C en las Regiones de Cabecera y Aguas Abajo de la Cuenca del Río Indo

Planteamiento del Problema
Las cuencas fluviales de montaña, particularmente en el Alto Asia Montañosa, son altamente sensibles al calentamiento climático, el cual altera la magnitud, el tiempo y la fase del escurrimiento. Si bien estudios previos han establecido la sensibilidad de la Cuenca Alta del Indo (UIB) —dominada por la nieve y los glaciares— al calentamiento, existe una cuantificación insuficiente de cómo estos signos se propagan hacia la Cuenca Baja del Indo (LIB). La LIB es un sistema fuertemente gestionado y dependiente de la irrigación, donde la dependencia de las aguas subterráneas y la regulación humana modifican significamente el flujo natural. Existe una brecha de conocimiento crítica sobre si el calentamiento global progresivo produce una divergencia hidrológica donde las ganancias de escurrimiento aguas arriba coinciden con déficits aguas abajo. Además, muchas evaluaciones dependen de periodos de tiempo fijos en lugar de Niveles de Calentamiento Global (GWL) equivalentes, lo que no logra contabilizar las variaciones en las tasas de calentamiento entre modelos y trayectorias de emisiones. Este estudio aborda estas brechas evaluando las respuestas del escurrimiento a los GWL de 1.5°C, 2.0°C y 3.0°C en todo el sistema completo de la cuenca superior-inferior.

Metodología
El estudio emplea un marco a escala de estación, integrando datos de 35 estaciones de medición hidrológica en la Cuenca del Río Indo (IRB).

• Forzamiento de Datos: El análisis utiliza datos de forzamiento hidroclimático con sesgo ajustado provenientes del Proyecto de Intercomparación de Modelos de Impacto Intersectorial 3b (ISIMIP3b) y del Proyecto de Intercomparación de Modelos Acoplados Fase 6 (CMIP6). Se utilizan cuatro Modelos de Circulación General (GCM: CNRM-CM6-1, CNRM-ESM2-1, EC-EARTH3, MIROC6) bajo tres Trayectorias Socioeconómicas Compartidas (SSP1-2.6, SSP3-7.0, SSP5-8.5).
• Marco de Modelado: Un Transformador de Fusión Temporal (TFT), una arquitectura de aprendizaje profundo diseñada para el pronóstico de horizontes múltiples, es entrenado con registros observados de descarga mensual y los predictores climáticos con sesgo ajustado. El modelo incorpora 11 predictores dinámicos (p. ej., precipitación, temperatura, radiación, humedad) y covariables estáticas (latitud, longitud) para capturar controles hidroclimáticos no lineales.
• Extracción de GWL: En lugar de años fijos, las respuestas hidrológicas se extraen basándose en el alcance de GWL específicos (1.5°C, 2.0°C, 3.0°C) en relación con la línea base 1850–1900. El año de alcance se define como el punto medio de la primera ventana de 20 años que excede el umbral.
• Diagnósticos: El estudio cuantifica anomalías anuales y estacionales del escurrimiento, la fracción de nieve, los Meses Positivos de Grados (PDM), la sensibilidad al calentamiento (m³/s por °C) y la incertidumbre (rango intercuartílico, acuerdo entre modelos).

Resultados Clave
El análisis revela una divergencia hidroclimática robusta y longitudinalmente diferenciada entre la UIB y la LIB bajo un calentamiento progresivo:

• Divergencia de Escurrimiento:
  • Cuenca Alta del Indo (UIB): El escurrimiento anual promedio aumenta de forma monótona con el calentamiento, pasando de +172.80 m³/s a 1.5°C a +725.49 m³/s a 3.0°C. La sensibilidad al calentamiento es de +289.40 m³/s por °C. Este aumento viene acompañado de un rango de incertidumbre más amplio (el IQR aumenta de 46.35 a 137.02 m³/s).
  • Cuenca Baja del Indo (LIB): En contraste, los déficits de escurrimiento anual promedio se profundizan, pasando de −78.40 m³/s a 1.5°C a −438.86 m³/s a 3.0°C, con una sensibilidad de −187.20 m³/s por °C.
• Reorganización Estacional:
  • UIB: Las ganancias se concentran en las estaciones cálidas (JJA y SON), impulsadas por un deshielo incrementado, aunque la magnitud de las anomalías varía según el GWL.
  • LIB: Las pérdidas persisten durante MAM–SON, con solo cambios positivos modestos en invierno (DJF). El régimen estacional en la LIB se debilita estructuralmente y es cada vez más negativo.
• Impulsores Criosféricos:
  • En la UIB, el aumento del escurrimiento coincide con una disminución sistemática de la fracción de nieve (de −0.03 a 1.5°C a −0.06 a 3.0°C) y un aumento brusco de los PDM (de +9.40 a +31.60). Esto indica un cambio hacia regímenes dominados por la lluvia y un aumento de la energía de deshielo térmica.
  • La LIB muestra aumentos de temperatura similares pero carece de la amplificación criosférica vista aguas arriba, reflejando en su lugar los efectos compuestos de la evaporación y la regulación humana.
• Heterogeneidad Espacial:
  • El análisis a escala de estación confirma que las señales a escala de cuenca emergen de respuestas locales heterogéneas. Las estaciones de la UIB (p. ej., Gilgit, Shigar) muestran una fuerte sensibilidad positiva, mientras que las estaciones de la LIB (p. ej., Guddu, Kotri) muestran una fuerte sensibilidad negativa. Algunas estaciones intermedias muestran señales débiles o no significativas, resaltando la incertidumbre a escala local.

Significancia y Afirmaciones
El artículo afirma que la Cuenca del Río Indo no responde al cambio climático como un sistema espacialmente uniforme. En su lugar, el calentamiento impulsa una reorganización longitudinal pronunciada donde la mejora del escurrimiento en las cabeceras no se traduce en disponibilidad de agua aguas abajo.

• División Hidroclimática: El estudio demuestra que la IRB funciona como un sistema acoplado pero funcionalmente heterogéneo. La UIB actúa como un sistema condicionado por la criosfera y dirigido por la energía, donde el calentamiento amplifica temporalmente el escurrimiento mediante el deshielo. Por el contrario, la LIB actúa como un sistema altamente regulado donde la descarga aguas abajo se desacopla de las ganancias de aguas arriba debido a las extracciones para irrigación, las derivaciones de canales y la alta demanda evaporativa.
• Naturaleza Transitoria de las Ganancias: Los autores enfatizan que el aumento proyectado del escurrimiento en la UIB debe interpretarse como una fase transitoria de "pico de agua" asociada con un deshielo incrementado bajo un almacenamiento criosférico decreciente, en lugar de una ganancia permanente de recursos hídricos.
• Relevancia Política: Los hallazgos subrayan que la planificación de la adaptación no puede basarse únicamente en las proyecciones de escurrimiento de las cabeceras. La divergencia plantea riesgos específicos para la confiabilidad de la irrigación, la dependencia de las aguas subterráneas y la seguridad alimentaria en la LIB. El estudio proporciona un marco físicamente fundamentado para evaluar riesgos en umbrales de calentamiento específicos (1.5°C, 2.0°C, 3.0°C), ofreciendo información relevante para la toma de decisiones en uno de los sistemas fluviales más estresados del mundo.

Limitaciones
Los autores señalan que el marco TFT captura la respuesta neta del escurrimiento sin separar explícitamente los procesos hidrológicos (p. ej., partición lluvia-escurrimiento frente a deshielo glaciar). Además, aunque los GWL mejoran la comparabilidad entre modelos, persisten las incertidumbres respecto a la selección de los GCM y la corrección de sesgos. El estudio no modela explícitamente la dinámica del balance de masa de los glaciares ni los esquemas de operación de embalses, lo que sugiere la necesidad de una futura integración con modelos basados en la física.