Parametrización de la cuenca hidrográfica del río Luján, Argentina

Resumen: El manejo de cuencas esta orientado a formular un conjunto integrado de acciones sobre el medio natural y social con fines de promover la conservación y la utilización sustentable de los recursos hídricos. La geomatica, como disciplina científica, aporta numerosas herramientas metodológicas que contribuyen al análisis de las cuencas hidrológicas como objeto de estudio. En la actualidad la posibilidad de contar con datos provenientes de sensores satelitales, disponibles en la web, junto con el creciente desarrollo de software libre permiten el libre acceso a la información y su tratamiento a través de software específico, contribuyendo a optimizar la toma de decisiones en ámbitos gubernamentales. El siguiente trabajo, enmarcado en el proyecto de doctorado denominado Manejo integral de cuencas hidrográficas basada en Geomática, tiene como principal objetivo generar insumos temáticos que serán utilizados en la elaboración de indicadores ambientales del área bajo análisis. La zona seleccionada corresponde a una cuenca de llanura con pendientes muy suaves ubicada en el noreste de la provincia de Buenos Aires cuyo principal curso de agua es el río Luján, que se encuentra amenazado por las presiones del crecimiento urbano del área metropolitana. Los resultados obtenidos corresponden a diferentes capas temáticas relacionadas con parámetros como: pendiente, área de captación, limitación de subcuencas, red hídrica secundaria y escurrimiento superficial, los que permitirán contribuir en la elaboración de una propuesta de manejo integral del recurso.

Palabras clave: Geomática, manejo de cuencas, software libre, modelo de elevación digital

1. Introducción

Actualmente el desarrollo y la utilización de herramientas informáticas basadas en la filosofía del software libre es cada vez más extendida. Los Sistemas de Información Geográfica (SIG) no son una excepción, planteándose como verdaderas alternativas frente al desarrollo comercial de estas tecnologías. La gran difusión que presentan hoy el uso de los SIG en pequeñas instituciones, municipios, oficinas públicas y organizaciones no gubernamentales, donde no es necesario un tratamiento comercial de los datos, ha favorecido la aplicación de este tipo de software ya que el costo de adquisición en estas entidades representa un factor determinante, presentando así a los SIG libres como una verdadera opción (Jiménez Berni, J., et al, 2005). El desarrollo de las tecnologías espaciales ha evolucionado de manera ininterrumpida desde la puesta en órbita de los primeros satélites; las imágenes o productos elaborados de resultantes de la observación de los recursos terrestres y la dinámica de ciertos fenómenos se encuentran disponibles en diferentes servidores y suponen una fuente de datos para los SIG (Buzai, G., 2000).

Su aplicación a los recursos hídricos y en particular al manejo de cuencas, resulta una herramienta indispensable para conocer el ambiente y de esta manera hacer un uso sustentable del mismo. Un ejemplo de ello lo configuran los parámetros geomorfológicos derivados de imágenes topográficas que permiten establecer el comportamiento de determinadas variables y ser incluidos en modelos más complejos como la ecuación universal de pérdida de suelos (USLE), entre otros. Las imágenes de la Shuttle Radar Topography Mission (SRTM) disponibles gratuitamente en la web favorecen la extracción de estos parámetros brindando información topográfica del área (NASA, 2005). El objetivo de este trabajo fue generar insumos temáticos a partir de parámetros geomorfométricos derivados de datos topográficos en la cuenca del río Luján. Los mismos serán utilizados para la elaboración de una propuesta de manejo integral del recurso.

2. Materiales y Métodos

Área de Estudio

El área bajo estudio se encuentra ubicada en el este de la provincia de Buenos Aires, como se puede observar en la Figura 1. Desde el punto de vista geomorfológico, pertenece a la superficie superior de acumulación de depósito del Pampeano, caracterizada por loess y limo de diversas fuentes, y presenta una leve elevación con respecto al nivel del mar de 30 metros, lo que le confiere un aspecto llano. Se ubica en la región denominada Pampa Ondulada.

Su red de drenaje superficial presenta como principal componente al río Luján, típico río de llanura, que nace en la confluencia de los arroyos El Durazno y Los Leones, en el partido bonaerense de Suipacha y luego de recorrer los partidos de Mercedes, Luján, Pilar, Exaltación de la Cruz, Campana, Escobar y San Fernando desemboca en el Río de la Plata. Este río recibe los aportes de numerosos arroyos, presentando un régimen regular a través de todo el curso, alimentándose principalmente de precipitaciones. Las especies forestales que se encuentran en la actualidad son en su mayoría exóticas introducidas al país en diferentes etapas o bien nativas de otras regiones ya que la región no se caracteriza por especies forestales nativas, con excepciones como el tala (celtis tala).

El clima del partido, de acuerdo a la clasificación efectuada por Koeppen, pertenece al tipo CFa, templado húmedo con verano muy cálido. Los vientos predominantes son del NENO, correspondiendo a masas húmedas y SO a masas secas. Las precipitaciones presentan un promedio anual de 950 mm, el promedio de temperaturas estivales es de 25º C mientras que las invernales no superan, en promedio, los 9,5º C con baja amplitud térmica.

Materiales

Para la realización de este trabajo se utilizó una imagen correspondiente a la Shuttle Radar Topography Mission (SRTM) path/row 225/084 con una resolución espacial de 90 metros, de fecha febrero de 2000. Es generada a partir de imágenes radar en banda C cuya longitud de onda es de 5,6 cm. El tipo de SRTM seleccionada corresponde al Unfilled_finished_B, que de acuerdo a la bibliografía, presenta pocos valores sin completar y es el sugerido para este tipo de análisis (Sanders B., 2007; Verdin K., et al, 2007). Estos valores son representados por el número -32768 y deben ser interpolados en el procesamiento. El formato de distribución es Geotiff y su resolución radiométrica es de 16 bit (Farr, 2007).

Se contó con material cartográfico auxiliar correspondiente a coberturas vectoriales de diferentes organismos públicos, las mismas correspondieron a la red hídrica superficial, poblados, vías de comunicación y límites administrativos. La Figura 2 muestra la escena 225/084, procesada inicialmente en donde se puede observar la configuración topográfica y su red hídrica superficial.

Software Utilizado
Teniendo en cuenta la utilización de software de libre acceso del Open Geospatial Consortium (OGC) se utilizaron los siguientes programas:

– QUANTUM GIS
– GRASS GIS
– SAGA GIS

Métodos

La imagen inicial fue seleccionada y descargada del servidor online de la Universidad de Maryland (ESDI-GLCF-EE.UU). Inicialmente se utilizó el complemento de GRASS para QGIS que permitió generar los límites de la cuenca hídrica del río Luján y su red de drenaje principal y secundaria. Para ello fue necesario corregir las anomalías presentes en los datos del SRTM y crear subcapas para la delimitación, como la dirección de flujo de agua de cada celda y la cantidad acumulada de celdas drenando en cada una. Una vez definida y recortada el área, utilizando el software SAGA GIS, se convirtió a formato propio y se derivaron los diferentes parámetros geomorfométricos correspondientes al área de estudio utilizando la herramienta Basic Terrain Analysis. La Figura 3 muestra un esquema con la metodología empleada partiendo de la imagen SRTM inicial. A continuación se enumeran y describen brevemente los principales parámetros obtenidos (Verdin, K., et al, 2007; Lehner, B. 2008; Tzelepis, N., 2008).

Sombreado analítico: es una simulación de sombreado basada en valores de iluminación solar sobre la superficie terrestre con propósitos cartográficos. El efecto visual esta dado por la diferenciación en los tonos de la imagen resultado.

Pendiente: definida como el ángulo existente entre el vector normal a la superficie en ese punto y la vertical.

Aspecto: se puede definir como la dirección de la pendiente y se mide en grados (0 a 360) en tanto que las zonas planas poseen un valor de -1.

Plano de curvatura: parámetro utilizado para analizar las crestas y los valles. Su valor es positivo con celdas de contorno cóncavo mientras que negativo para las celdas con contornos convexos. Corresponde a la segunda derivada como el siguiente parámetro.

Perfil de curvatura: es la curvatura de la superficie en la dirección de la pendiente más pronunciada. Este parámetro esta relacionado con la deposición y la erosión de la red hídrica y sus márgenes.

Índice de convergencia: se calcula teniendo en cuenta las celdas adyacentes, su aspecto y su relación con la celda central. Determina como decrece o crece la celda circundante en relación a la central.

Zona de captación: es el área geográfica donde el agua superficial drena hacia un punto determinado.
Índice de humedad: se utiliza para describir patrones espaciales de humedad en una región. Se calcula a partir de la pendiente, siendo húmedas aquellas zonas donde las pendientes son bajas y sus áreas de contribución grandes, mientras que las secas corresponden a mayor pendiente y menor área de contribución.

Factor LS: factor de inclinación de la pendiente y longitud, muy utilizada en la ecuación revisada universal de perdida de suelo (USLE y RUSLE).

3. Resultados y Conclusiones

En el caso particular del sombreado analítico de la cuenca se observa, a través de las diferentes tonalidades, dos grandes áreas de nacimiento y desembocadura, no obstante, debido a las características de la zona fue necesaria la implementación de un elevado factor de exageración. En la desembocadura se observa un área plana que corresponde a una zona de humedales. La pendiente presenta sus mayores valores en la zona media del área correspondiendo a valores máximos que no superan el dos por ciento; en el aspecto, que expresa la dirección de la pendiente, se observan las zonas planas del nacimiento y la desembocadura, mientras que la zona media presenta la mayor variabilidad. Tanto el plano como el perfil de curvatura indican que la zona presenta valores muy bajos de erosión potencial de suelos. El índice de convergencia presenta una alta variabilidad, en particular en la zona media, observándose las mayores áreas en ambos extremos. Los valores correspondientes a la zona de captación indican la presencia de subcuencas atomizadas mientras que el patrón espacial de humedad, establecido como índice de humedad, no muestra grandes zonas diferenciables, salvo las dos áreas de los extremos, mencionadas anteriormente.
Tanto el factor LS como la red de canales no aportan datos relevantes. La aplicación de la metodología planteada permitió derivar del SRTM parámetros morfométricos que ayudaron a la caracterización del área de la cuenca. Por su escasa pendiente y la dinámica que presenta la zona de estudio convierte a los eventos hidrológicos, como las inundaciones, en amenazas tanto para la población como para los recursos productivos.

4. Bibliografía

Bates, P.D. and De Roo A.P.J., 2000. A simple raster-based model for flood inundation simulation. J. Hydrol. 236, 54-77. United Kingdom
Bruno, J., 2000. Manejo de Cuencas Hidrográficas, Análisis Conceptual. Magister Scientiae en Manejo de
Cuencas Hidrográficas. Facultad de Ciencias Agrarias y Forestales, Universidad Nacional de La Plata. Argentina
Buzai, G., 2000. La exploración geodigital. Ed. Lugar Editorial, Buenos Aires, Argentina.
Cabrera, Á.L., 1976. Regiones Fitogeográficas Argentinas; fascículo 1; Editorial ACME S.A.C.I. Argentina.
Castellon Pineda, N., 2004. Análisis socioambiental del uso y manejo del agua en la subcuenca del Río Aguas
Calientes, Somoto, Nicaragua Tesis de Maestría de la Escuela de Posgrado de CATIE. Turrialba, Costa Rica. p 168.
Chuvieco, E., 1996. Fundamentos de Teledetección Espacial. Tercera ed. rev., Ediciones Rialp S.A., Madrid, España. 568 p.
Dourojeanni, A., 2009. Los desafíos de la gestión integrada de cuencas y recursos hídricos en América Latina y el Caribe, en Boletín electrónico REDLACH, Santiago de Chile, Chile.
Farr, T. G., 2007. The Shuttle Radar Topography Mission, Rev. Geophys., 45, RG2004, doi:10.1029/2005RG000183. p 21- 22.
Goldberg, S., Cirera I., Parella M., Benitez A., Bulos L, y Troncoso A.,1995. Caracterización climática y agroclimática de la cuenca del Río Luján. Resúmenes Jornada sobre la cuenca del Río Luján; pp13-19. Argentina.
Jiménez Berni, J., Aguilera Ureña, M.J., Meroño de Labriva, J., 2005. Alternativas de software libre a los Sistemas de Información Geográfica comerciales. Congreso Internacional Conjunto XVII Ingegraf. Sevilla, España.
Koppen, W., 1936. Das geographisca System der Klimate, in: Handbuch der Klimatologie, edited by: Koppen
W. and Geiger G., 1. C. Gebr, Borntraeger, 1–44 p.
Lehner, B., Verdin, K., and Jarvis, A., 2008. New Global Hydrography Derived From Spaceborne Elevation
Data: Eos, Transactions, American Geophysical Union, v. 89, no. 10, p. 93 – 100.
Manual SAGA: http://volaya.es/pdf/SagaManual.pdf (consulta 2-10-2012)
Manual Quantum GIS: http://qgis.org/en/documentation/manuals.html (consulta 3-9-2011)
Manual Grass: http://grass.fbk.eu/gdp/manuals.php (consulta 17-1-2011)
NASA, 2005. Shuttle Radar Topography Mission: Instruments. http://www2.jpl.nasa.gov/srtm/instr.htm (Version 06 August 2008).
Sanders, B. F., 2007. Evaluation of on-line DEMs for flood inundation modeling: Advances in Water Resources,
v. 30, p.1831 – 1843.
Tzelepis, Ν., Loisios, D., Nakos Β., 2008. A new method of analytical hillshading by using different lighting directions. Proceedings of the 10th National Cartographic Conference, Hellenic Cartographic Society, Ioannina.
Verdin, K.L., Godt, J.W., Funk, C., Pedreros, D., Worstell, B., Verdin, J., 2007. Development of a global slope dataset for estimation of landslide occurrence resulting from earthquakes: Colorado: U.S. Geological Survey, Open-File Report 2007-1188, 25 p.

Por: Leonardo Di Franco
Alfredo Cuello
María Cristina Serafini
Fuente: www.selperbrasil.org.br

VER ARCHIVO DE ARTÍCULO