Evaluación espacial de la presión agrícola del cultivo de café y anomalías climáticas en subzonas hidrográficas (szh) Cafeteras
1. Introducción
La gestión del riesgo agroclimático en sistemas de ladera y cuencas hidrográficas requiere herramientas analíticas que integren la presión antrópica (demanda productiva) con las variaciones climáticas coyunturales (oferta hídrica).
El presente documento detalla la metodología y el flujo de trabajo técnico desarrollado para evaluar la vulnerabilidad frente a temporadas de sequía en subzonas hidrográficas con presencia de cultivos de café. El proceso combina análisis vectorial en sistemas de información geográfica (QGIS) con procesamiento masivo de datos de teledetección en la nube mediante Google Earth Engine (GEE).
2. Metodología y Fases del Proceso
El desarrollo técnico se estructuró en tres fases principales: cálculo de la presión agrícola estructural, procesamiento climático multitemporal y cruce de variables para la generación de indicadores de riesgo.
Fase 1: Cálculo del Índice de Presión Agrícola del Café (I_pa)
Para determinar la intensidad del uso del suelo agropecuario en cada unidad territorial, se relacionó el área sembrada de café (datos de las evaluaciones agropecuarias EVA 2025) con la superficie total de cada sub zona hidrográfica (szh).
Procedimiento: Mediante operaciones geométricas y tabulares en QGIS, se intersectaron las capas vectoriales con las areas de café por municipio y las de las sub zonas hidrográficas. Se obtienen áreas de siembra del cultivo por cada sub zona hidrográfica y la relación area de siembra / área de la sub zona hidrográfica.
Interpretación: Este cociente adimensional (expresado de 0 a 1 o en porcentaje) cuantifica qué proporción del territorio de la cuenca está comprometida con la caficultura, reflejando indirectamente la magnitud de la demanda hídrica potencial para procesos productivos y de beneficio.
Fase 2: Análisis de Anomalías de Precipitación en Google Earth Engine (GEE)
El comportamiento de la oferta hídrica superficial se evaluó utilizando series de tiempo de precipitación diaria del sensor CHIRPS (Climate Hazards Group InfraRed Precipitation with Station data).
Extracción y procesamiento en GEE:
Acumulado actual: Se calculó la precipitación acumulada para el periodo de análisis del año en curso (hasta sep. 2026), recortada a la geometría de las subzonas de estudio.
Climatología de referencia: Se procesó una serie histórica de 20 años (2006-2025) para los mismos rangos de fechas y la misma geometría, obteniendo el promedio histórico esperado.
Cálculo de anomalías: Se derivó la anomalía relativa porcentual (
Rel_Anom_pct) para medir el grado de desviación frente a la normal histórica:
Fase 3: Integración de Variables y Estandarización
Los datos tabulares climáticos extraídos de GEE se importaron e integraron en QGIS con el geopackage original de las subzonas hidrográficas.
Normalización del Déficit Climático (I_dc): Los valores de anomalía relativa se clasificaron de manera que los excedentes de lluvia (valores positivos) se fijaran en cero, mientras que los déficits (valores negativos) se transformaron en un índice de severidad de sequía estandarizado entre 0 y 1.
Modelo de Riesgo Cruzado: Mediante geoprocesamiento en QGIS, se integró el I_pa (presión) y el I_dc (déficit climático) para identificar aquellas subzonas que experimentan una doble vulnerabilidad: alta densidad cafetera coincidiendo con las mayores anomalías negativas de precipitación.
3. Resultados y Aplicaciones Prácticas
El flujo de trabajo implementado permite obtener los siguientes productos cartográficos y analíticos:
Matriz Tabular Integrada: Un conjunto de datos a nivel de subzona que incluye área total, área cafetera, presión agrícola, precipitación acumulada, valores históricos y porcentajes de anomalía.
Mapas de Amenaza Climática: Visualización espacial discreta de las cuencas bajo estrés hídrico severo.
Focalización de Medidas de Mitigación: Permite a las entidades de planificación orientar recursos técnicos, asistencia en gestión eficiente del agua y planes de contingencia exclusivamente hacia las subzonas hidrográficas donde la sequía representa un riesgo real para la productividad cafetera.
Comentarios