LA GESTIÓN DEL AGUA EN LAS CUENCAS HÍDRICAS
Jordi Casals Fernandez
Barcelona, 10 de junio de 2020
Máster en Sistemas de Información Geográfica
Datos SIG e Indicadores Medioambientales
Introducción
Nos han contratado para elaborar un informe analizando la evolución del embalse de Darnius-
Boadella y la correlación de ésta con la pluviometría. Para ello deberemos seguir el flujo de
trabajo del análisis espacial.
• Definir problema y formular pregunta
• Dividir el estudio en submodelos
• Explorar y preparar los datos
• Elegir métodos y herramientas
• Analizar los resultados
Para la realización del estudio se utilizarán datos que tengan un recubrimiento de la totalidad
del planeta con el objetivo de obtener un flujo de trabajo replicable en cualquier zona del
mundo. Además, se compararán los resultados con los datos reales obtenidos de la medida
directa por la estación meteorológica situada en la presa.
Definir problema i formular pregunta
En Catalunya, la Agencia Catalana del Agua (ACA) es la encargada de la gestión, mantenimiento
y análisis de los recursos hídricos de las cuencas internas. Para ello elaboran informes anuales
sobre la situación:
En Cataluña, el régimen pluviométrico se caracteriza por una gran irregularidad
general y por una variabilidad interanual elevada, fruto del clima típico
mediterráneo, lo que nos hace especialmente vulnerables a sufrir episodios de
sequía, especialmente en las cuencas internas que es donde se concentra el grueso
de la población (más del 80%) pero sólo el 40% de la disponibilidad de agua,
esencialmente la que se encuentra en los embalses, los acuíferos y dentro de los
pozos.
La Agencia Catalana del Agua (ACA) hace un control exhaustivo de las reservas
hídricas disponibles, de sus oscilaciones y de posibles tendencias con el fin de
garantizar una gestión flexible de todas las fuentes de abastecimiento y no agotar
una en detrimento de otra (Figura 1).
Cataluña vive desde primeros de 2016 un periodo hidrológicamente muy seco,
donde la pluviometría ha sido muy escasa y muy por debajo de la media de los
últimos 5 años. Actualmente se mantiene la sequía meteorológica y casi el 75% de
los municipios de cuencas internas se encuentran en sequía extrema por falta de
lluvia.
Figura 1. Evolución pluviometría 2016-2017. Fuente: Agencia Catalana del Agua (ACA) 1
En estas condiciones críticas se plantea la necesidad del análisis de la pluviometría en la
cuenca hidrográfica del embalse de Darnius-Boadella y estudiar la relación de estas con la
superficie de la lámina de agua y el estado de la vegetación que en épocas de sequía podrían
favorecer la propagación de incendios.
Dividir el estudio en submodelos
Para afrontar cualquier análisis espacial de complejidad debe subdividirse en submodelos que
permitan afrontar el estudio por partes. En este caso se dividirá en cuatro modelos:
• Establecer cuenca hidrográfica
• Evolución pluviométrica
• Evolución embalse de Darnius-Boadella
• Evolución estado de la vegetación
Explorar y preparar los datos
En el proceso del análisis territorial uno de los pasos de mayor importancia es la
documentación y la obtención de los datos base para el desarrollo del estudio. Para el
desarrollo del estudio se deberá obtener un Modelo Digital del Terreno (MDT), datos
pluviométricos e imágenes de satélite multibanda.
Además de descargar los datos también es importante obtener la información de cómo han
sido tomados los datos (formato, resolución, unidades…) para poder trabajar con ellos.
Paso 1. Documentación embalse Darnius-Boadella
Obtener información relacionada con el embalse (Situación, datos generales, historia…) antes
de empezar cualquier estudio es necesario conocer el ámbito de estudio.
1
http://aca-web.gencat.cat/aca/sequera/es/situacio-actual.jsp
Paso 2. Obtención MDT
El Servicio Geológico de los Estados Unidos genera un MDT de todo el planeta a una resolución
de 90 metros. Se puede obtener en:
https://www.sciencebase.gov/catalog/item/5920dd2de4b0ac16dbdf3a39
Paso 3. Obtención datos pluviométricos
El mismo servicio geológica de los EEUU mediante el grupo CHIRPS (Climate Hazards Group
InfraRed Precipitation with Station data) genera datos pluviométricos entre las latitudes 50ºS y
50ºN con una resolución de 0.05º.
Descargaremos las precipitaciones mensuales de todo el 2017. Se puede obtener en:
ftp://ftp.chg.ucsb.edu/pub/org/chg/products/CHIRPS-2.0/global_monthly_EWX/data/
Paso 4. Obtención imágenes por satélite
Las misiones Landsat permiten obtener imágenes con una resolución de 30 metros (15 metros
para les EEUU) de todo el planeta con una resolución temporal de 16 días.
Descargaremos imágenes Landsat8 de todo el 2017. Antes de solicitar la descarga es
importante analizar qué imágenes no contienen nubes que interfieran en el análisis. Como
mínimo se utilizará una imagen por estación, distribuidas lo más homogéneamente posible. Se
puede obtener en:
https://earthexplorer.usgs.gov/
Paso 5. Preparación de los datos. Recortar (Opcional)
Al tratarse de datos de gran extensión, las herramientas de geoprocesamiento pueden tardar
mucho tiempo en realizarse. Para ello se recomienda recortar los datos sólo para la zona de
estudio mediante la herramienta clip.
Data Management Tools / Raster / Raster Processing / Clip
Paso 6. Preparación de los datos. Reproyectar
Como se ha observado, todos los datos están generados en distintos sistemas de referencia y
proyección. Para ello se transformarán los datos, mediante la herramienta Project Raster, a
ETRS89 UTM 31N, el oficial en la zona de estudio.
Data Management Tools / Projections and Transformations / Raster / Project Raster
Paso 7. Preparación de los datos Landsat8
Para poder realizar el cálculo de los diferentes índices que nos permitirán analizar el territorio,
debemos transformar los valores digitales (Digital Numbers) a reflectancia mediante la
calculadora gráfica.
Data Analyst Tools / Map Algebra / Raster Calculator
Las ecuaciones para realizar la transformación de los DN a reflectancia se pueden obtener
en la documentación oficial de Landsat8 (página 60).
https://landsat.usgs.gov/sites/default/files/documents/Landsat8DataUsersHandbook.pdf
Elegir métodos y herramientas. Cuenca hidrográfica
El primer paso será delimitar la cuenca hidrográfica que delimitará todo aquel territorio que
desagüe en el embalse de Darnius-Boadella. Para ello se utilizarán herramientas de hidrología,
dentro de las herramientas de análisis espacial.
ArcToolbox / Spatial Analyst Tools / Hydrology
Paso 8. Eliminar imperfecciones
Para evitar que puedan producirse errores en los siguientes cálculos, debe corregirse el MDT
para eliminar huecos y sumideros que pueda contener el raster con la herramienta Fill.
Paso 9. Direccion del flujo
Establecer la dirección del flujo hidrológico de la pendiente con la herramienta Flow Direction.
Paso 10. Acumulación del flujo
Determinar la acumulación del flujo de las celdas que fluyen hacia cada celda descendiendo
sobre la pendiente por medio de la herramienta Flow Accumulation.
Paso 11. Punto de desfogue
El siguiente paso determinar el punto de desfogue o drenaje de la microcuenca, esto se puede
hacer creando un shapefile tipo punto y añadiendo manualmente un punto cerca de la presa
donde se acumula el flujo de agua. A continuación, interpolar ese punto con el modelo de
elevación digital (TIN o DEM) para obtener sus coordenadas en tres dimensiones.
ArcToolbox / 3D Analyst Tools / Functional Surface / Interpolate Shape
Figura 2. Situación aproximada del punto. Fuente: Propia
Paso 12. Determinación cuenca
Para culminar con la herramienta Watershed usar el ráster creado con Flow Direction y el
punto de desfogue interpolado.
Paso 13. Vectorización cuenca
Finalmente, para obtener el polígono de la microcuenca tan solo basta en convertir el ráster a
shapefile tipo polígono.
ArcToolbox / Conversion Tools / From Raster / Raster to Polygon
Figura 3. Resultado cuenca hidrográfica. Fuente: Propia
Elegir métodos y herramientas. Pluviometría
Para el análisis de la pluviometría en la cuenca hidrográfica se obtendrá la precipitación
acumulada en la región de la cuenca mediante los datos CHIRPS y se comparará con la estación
meteorológica más cercana.
Paso 14. Agregar campos
A la cuenca creada anteriormente, añadir dos campos al polígono de la cuenca: Mes (Tipo
texto) y CHIRPS (Tipo float). Iniciar una sesión de edición para introducir el nombre de los
meses (debe coincidir con el archivo Estacion.xlsx adjuntado a la documentación del ejercicio).
Paso 15. Duplicar entidades
Con la sesión de edición iniciada, selecciona el polígono de la cuenca y seguidamente copia y
pega, para agregar un nuevo elemento dentro de la misma clase de entidad de la cuenca.
Realiza este paso hasta obtener 12 polígonos, uno para cada mes del año.
Paso 16. Media CHIRPS
La Herramienta Zonal Statistics as Table permite crear una tabla con diferentes cálculos
estadísticos (Media, mínimo, máximo…) de los pixeles contenidos en un polígono.
Calcularemos la media para cada mes y la introduciremos al polígono de la cuenca.
Spatial Analyst Tools / Zonal / Zonal Statistics as Table
Paso 17. Valor Estación
Hacer un Join entre la cuenca y la tabla Estación.xlsx, utilizando el campo mes para realizarla.
La tabla asociada a la cuenca debería parecerse a la tabla 1.
OBJECTID Shape_Leng Shape_Area Mes CHIRPS Estación
1 70265,59 183619147,9 Enero XX,XX 100,8
2 70265,59 183619147,9 Febrero XX,XX 121,4
3 70265,59 183619147,9 Marzo XX,XX 70,8
4 70265,59 183619147,9 Abril XX,XX 32,4
5 70265,59 183619147,9 Mayo XX,XX 34,5
6 70265,59 183619147,9 Junio XX,XX 75
7 70265,59 183619147,9 Julio XX,XX 19,9
8 70265,59 183619147,9 Agosto XX,XX 22,2
9 70265,59 183619147,9 Septiembre XX,XX 48,5
10 70265,59 183619147,9 Octubre XX,XX 28
11 70265,59 183619147,9 Noviembre XX,XX 3,5
12 70265,59 183619147,9 Diciembre XX,XX 22,5
Tabla 1. Resultado aproximado de la tabla con las precipitaciones. Fuente: Propia
Elegir métodos y herramientas. Lámina de agua
Para obtener la superficie del embalse por lo menos se utilizará el índice normalizado de
diferencias de agua modificado (MNDWI. Modified normalized difference water index).
Paso 18. Calculo MNDWI
Se valorará probar distintos índices y comparar los distintos resultados, para ello se puede
utilizar el siguiente link o cualquier otro artículo donde se proponga un índice de agua.
https://www.int-arch-photogramm-remote-sens-spatial-inf-sci.net/XLII-4-
W4/185/2017/isprs-archives-XLII-4-W4-185-2017.pdf
Utilizar la calculadora gráfica para calcular el MNDWI siguiendo la siguiente fórmula que
combina las bandas 3 y 6.
𝑏𝑏𝑔𝑔𝑔𝑔𝑔𝑔𝑔𝑔𝑔𝑔 − 𝑏𝑏𝑠𝑠𝑠𝑠𝑠𝑠𝑠𝑠1
𝑀𝑀𝑀𝑀𝑀𝑀𝑀𝑀𝑀𝑀 =
𝑏𝑏𝑔𝑔𝑔𝑔𝑔𝑔𝑔𝑔𝑔𝑔 + 𝑏𝑏𝑠𝑠𝑠𝑠𝑠𝑠𝑠𝑠1
Paso 19. Reclasificar MNDWI
En un caso teórico el MNDWI genera un raster donde los valores positivos se corresponden
con agua así, con la herramienta Reclassify, clasificaremos el raster en dos clases: de -1 hasta 0
y de 0 hasta 1. En alguno de los casos puede que el resultado no sea el óptimo y tenga que
modificarse para obtener un resultado correcto.
Spatial Analyst Tools / Reclass / Reclassify
Paso 20. Vectorización embalse
Para obtener el polígono del embalse debemos convertir el raster reclasificado a una entidad
tipo polígono. Asegurarse que se exporta en una Geodatabase para que se generen los campos
Perímetro y Superficie.
ArcToolbox / Conversion Tools / From Raster / Raster to Polygon
Paso 21. Añadir superficie
Finalmente creamos un nuevo campo tipo float a la entidad de cuenca y añadimos la superficie
del polígono que se corresponde al embalse en una sesión de edición.
Elegir métodos y herramientas. Estado vegetación
Para obtener el estado de la vegetación se utilizará el índice normalizado de diferencias de
vegetación (NDVI. Normalized difference vegetation index).
Paso 22. Calculo NDVI
Utilizar la calculadora gráfica para calcular el NDVI siguiendo la siguiente fórmula que combina
las bandas 5 y 4.
𝑏𝑏𝑁𝑁𝑁𝑁𝑁𝑁 − 𝑏𝑏𝑅𝑅𝑅𝑅𝑅𝑅
𝑀𝑀𝑀𝑀𝑀𝑀𝑀𝑀𝑀𝑀 =
𝑏𝑏𝑁𝑁𝑁𝑁𝑁𝑁 + 𝑏𝑏𝑅𝑅𝑅𝑅𝑅𝑅
Paso 23. Reclasificar NDVI
Con la herramienta Reclassify, clasificaremos los raster en 5 clases con los siguientes rangos:
• NDVI < 0, correspondientes a agua o cubiertas artificiales.
• 0 < NDVI < 0.2, suelo desnudo o vegetación muerta
• 0.2 < NDVI < 0.4, vegetación dispersa o poco vigorosa
• 0.4 < NDVI < 0.6, vegetación abundante y vigorosa
• NDVI > 0.6, vegetación muy densa y vigorosa
Se valorará probar distintos índices y comparar los distintos resultados, para ello se puede
utilizar el siguiente link o cualquier otro artículo donde se proponga un índice de agua.
https://www.researchgate.net/profile/Hanqiu_XU/publication/232724035_A_new_index_f
or_delineating_built-
up_land_features_in_satellite_imagery/links/09e4150904c665d747000000/A-new-index-
for-delineating-built-up-land-features-in-satellite-imagery.pdf?origin=publication_detail
Paso 24. Recortar
Utilizando la entidad cuenca como geometría de corte, recortar el raster reclasificado con el
NDVI, para eliminar las zonas fuera de la cuenca hidrográfica.
Data Management Tools / Raster / Raster Processing / Clip
Paso 25. Vectorización
Convertir el raster recortado con el NDVI reclasificado a una clase de entidad tipo polígono con
la herramienta Raster to Polygon. Recordar que debe exportarse en una Geodatabase para
que se genere el campo Superficie.
ArcToolbox / Conversion Tools / From Raster / Raster to Polygon
Paso 26 Resumen tabla
Desde la tabla de atributos de la entidad generada, resumir la tabla mediante el campo
gridcode que se ha generado con las distintas clases reclasificadas añadiendo la suma de las
superficies.
Paso 27. Añadir superficies
Añadir los valores obtenidos a la tabla de la entidad cuenca a su mes correspondiente, así se
obtendrán los datos necesarios para realizar el análisis.
OBJECTID Mes CHIRPS Estación Sup_Embalse Sup_NDVI1 Sup_NDVI2 Sup_NDVI3 Sup_NDVI4 Sup_NDVI5
1 Enero XX,XX 100,8 XX,XX XX,XX XX,XX XX,XX XX,XX XX,XX
2 Febrero XX,XX 121,4
3 Marzo XX,XX 70,8
4 Abril XX,XX 32,4 XX,XX XX,XX XX,XX XX,XX XX,XX XX,XX
5 Mayo XX,XX 34,5
6 Junio XX,XX 75
7 Julio XX,XX 19,9 XX,XX XX,XX XX,XX XX,XX XX,XX XX,XX
8 Agosto XX,XX 22,2
9 Septiembre XX,XX 48,5
10 Octubre XX,XX 28
11 Noviembre XX,XX 3,5 XX,XX XX,XX XX,XX XX,XX XX,XX XX,XX
12 Diciembre XX,XX 22,5
Tabla 2. Estructura aproximada de la tabla (Pueden variar los meses rellenados). Fuente: Propia
Analizar los resultados
Finalmente, se debe analizar los resultados obtenidos. En una primera fase debe comprobarse
que los resultados son coherentes con los resultados esperados, de no ser así, deberá
canalizarse la razón y corregirlo si es posible. En una segunda fase se tratará de responder las
hipótesis realizadas al inicio del estudio.
Paso 28. Comparación precipitaciones
Para comparar las precipitaciones obtenidas por los datos CHIRPS y por las estaciones, desde la
tabla de atributos se puede realizar un gráfico para visualizar los datos con la opción Create
Graph… También se puede exportar la tabla i trabajar con programas office.
Paso 29. Analizar evolución lámina de agua
Analizar si la evolución de la superficie de la lámina de agua del embalse se corresponde
directamente con las precipitaciones obtenidas o si hay más factores que influyan en este
aspecto.
Paso 30. Analizar evolución cobertura vegetal
Analizar si la evolución de la cobertura vegetal se corresponde directamente con las
precipitaciones obtenidas o si hay más factores que influyan en este aspecto.
Se valorará cualquier análisis adicional que se elabore.
