Seasonal rainfall predictability governed by interacting oceanic and terrestrial memory reservoirs

Resumen Técnico: La predictibilidad de las precipitaciones estacionales está gobernada por la interacción de reservorios de memoria oceánica y terrestre

Planteamiento del problema
La predictibilidad de las precipitaciones estacionales se ha atribuido tradicionalmente a anomalías oceánicas persistentes, particularmente al El Niño-Oscilación del Sur (ENOS) y al Dipolo del Océano Índico (DOI). El paradigma prevaleciente interpreta las anomalías de precipitación principalmente como respuestas a un forzamiento oceánico remoto. Sin embargo, esta visión suele confundir la sensibilidad al forzamiento con la dominancia predictiva. Los autores argumentan que la predictibilidad depende fundamentalmente de dónde se almacenan y se reintroducen las anomalías hidroclimáticas dentro del sistema climático. Si bien los océanos proporcionan una memoria a gran escala, las superficies terrestres (a través de la humedad del suelo y la evapotranspiración) también almacenan anomalías en escalas de tiempo semanales a estacionales. La pregunta central abordada es si la capacidad de predicción reside en el forzamiento oceánico remoto o en los reservorios terrestres que sostienen y reintroducen las anomalías, y cómo estas contribuciones se distribuyen entre diferentes regímenes hidroclimáticos.

Metodología
El estudio utiliza un marco de diagnóstico aplicado a la región del monzón del sur de Asia (60–110°E, 5–40°N) utilizando datos de reanálisis ERA5 de 1982 a 2025.

1. Marco de diagnóstico: Las anomalías de precipitación estacional ($P'$) se modelan como una combinación lineal de forzamiento oceánico (representado por el índice Niño3), anomalías de evapotranspiración ($ET'$) y anomalías de humedad del suelo en la zona radicular ($SM'$):
   $$P' = \alpha \cdot \text{océano}' + \beta \cdot ET' + \gamma \cdot SM' + \varepsilon$$
   Esta formulación es estadística y diagnóstica, destinada a particionar la influencia predictiva en lugar de establecer una causalidad mecánica.

2. Descomposición de la varianza: Utilizando la regresión de mínimos cuadrados ordinarios, la varianza total explicada ($R^2_{full}$) se particiona en:

   • Contribuciones únicas del forzamiento oceánico, la evapotranspiración y la humedad del suelo.
   • Covarianza compartida entre estos predictores.
     Las contribuciones únicas se aíslan mediante regresión parcial (por ejemplo, $R^2_{océano, única} = R^2_{full} - R^2(P' \sim ET', SM')$).

3. Índice de dominancia terrestre: Se define un índice interno ($D_{land}$) para distinguir entre la variabilidad impulsada por flujos (evapotranspiración) y la variabilidad mediada por el almacenamiento (humedad del suelo) dentro del bloque terrestre. Los valores negativos indican dominancia de la humedad del suelo.

4. Validación y pronóstico retrospectivo (Hindcasting):

   • Validación cruzada: Se utiliza una validación cruzada de dejar un año fuera (1982–2025) para evaluar la capacidad fuera de la muestra, evitando el sobreajuste dentro de la muestra.
   • Experimentos de pronóstico retrospectivo: El estudio compara una línea base de solo océano (utilizando el índice Niño3 de marzo–mayo antecedente) contra un modelo que incorpora la inicialización de la humedad del suelo en el pre-monzón.
   • Pruebas de robustez: Se prueban la sensibilidad a los esquemas de ponderación de la humedad del suelo vertical y la estabilidad temporal (mediante ventanas deslizantes de 25 años).

Resultados clave

• Redistribución espacial de la predictibilidad: Aunque la predictibilidad total de la precipitación estacional es sustancial en regiones semiáridas y de transición, la contribución única del forzamiento oceánico es espacialmente limitada y generalmente modesta ($<0.1$). En contraste, la humedad del suelo explica la mayor parte de la varianza independiente en regiones interiores.
• Jerarquía regional: Se analizaron tres regímenes hidroclimáticos distintos:
  • Meseta de Deccan (DP): Exhibe el control más fuerte mediado por el almacenamiento. La contribución única de la humedad del suelo es la más alta ($R^2 \approx 0.391$), mientras que la contribución única oceánica es mínima ($R^2 \approx 0.010$). El índice de dominancia terrestre es fuertemente negativo ($D_{land} = -0.652$).
  • Llanura Indo-Gangética (IGP): Muestra una estructura mixta pero ponderada por el almacenamiento. La humedad del suelo sigue siendo el predictor único dominante ($R^2 \approx 0.315$), con la influencia oceánica mediada en gran medida a través de la covarianza compartida con los estados terrestres.
  • Noroeste de la India (NWI): Representa un régimen fuertemente regulado por la tierra donde la variabilidad mediada por el almacenamiento supera al forzamiento oceánico, aunque la contribución oceánica es ligeramente más evidente que en la DP.
• Mediación de las señales oceánicas: El análisis de mediación revela que las asociaciones entre ENSO y la precipitación se debilitan sustancialmente tras condicionar por el almacenamiento terrestre. Esto sugiere que la variabilidad oceánica frecuentemente se proyecta sobre la precipitación a través de la covarianza compartida con la memoria del estado terrestre, en lugar de actuar como un predictor independiente.
• Ganancias en la capacidad de pronóstico: Los experimentos de pronóstico retrospectivo fuera de la muestra demuestran que la inicialización de modelos con la humedad del suelo del pre-monzón mejora significamente la capacidad de pronóstico estacional en regímenes dominados por el almacenamiento (específicamente en la Meseta de Deccan, con una ganancia de habilidad de correlación $\Delta r \approx 0.39$). En contraste, los modelos de solo océano exhiben una habilidad débil o negativa fuera de la muestra en estas regiones.
• Estabilidad temporal: La jerarquía de dominancia (DP > IGP > NWI) y la partición espacial del control predictivo se mantienen estables durante el período 1982–2025, lo que indica que son rasgos estructurales del sistema del monzón y no anomalías transitorias.

Significancia y afirmaciones
El artículo afirma establecer un principio unificador: la predictibilidad de la precipitación estacional no está gobernada únicamente por la magnitud del forzamiento remoto, sino por la persistencia y la organización espacial de los reservorios de memoria climática que interactúan.

• Distinción entre sensibilidad y dominancia: El estudio aclara que las fuertes correlaciones entre la precipificación y la variabilidad oceánica no implican que la capacidad de predicción resida en el forzamiento oceánico una vez que se consideran los componentes interactuantes. La ventaja predictiva se alinea con el reservorio que almacena y reintroduce las anomalías (la memoria terrestre en regiones semiáridas), no necesariamente con el forzamiento que inicia la variabilidad.
• Implicaciones para el pronóstico: Los resultados sugieren que los sistemas de pronóstico operativo, que a menudo priorizan la inicialización del estado oceánico, pueden beneficiarse de estrategias regionalmente adaptativas. En regímenes dominados por el almacenamiento, la representación e inicialización precisa de la memoria del estado de la tierra (humedad del suelo) puede ser tan crítica como la inicialización del estado oceánico.
• Contexto antropogénico: Debido a que la memoria terrestre es sensible al cambio de uso de suelo, el riego y la gestión hidrológica, las actividades humanas pueden redistribuir sistemáticamente la predictibilidad regional bajo un clima cambiante, incluso sin cambios en la amplitud del forzamiento oceánico.

Los autores concluyen que identificar dónde reside la ventaja predictiva requiere separar la sensibilidad de la precipitación al forzamiento de la persistencia de los reservorios que sostienen las anomalías, una distinción que no es capturada por los análisis convencionales de acoplamiento tierra-atmósfera.