Hidro-Estudfio Regional-Calamar-Boca Ceniza-U.n
Hidro-Estudfio Regional-Calamar-Boca Ceniza-U.n
Hidro-Estudfio Regional-Calamar-Boca Ceniza-U.n
componentes hidrológicos e
hidráulicos del complejo cenagoso
del Bajo Magdalena entre Calamar y
Bocas de Ceniza con apoyo de
herramientas de modelación
matemática y SIG
Tesis de investigación presentada como requisito parcial para optar al título de:
Magister en Ingeniería Civil – Recursos Hidráulicos
Director:
Ph.D. Erasmo Alfredo Rodríguez Sandoval
Línea de Investigación:
Modelación de Fenómenos y Amenazas Naturales
Grupo de Investigación:
GIREH - Grupo de Investigación en Ingeniería de los Recursos Hídricos
Al Instituto Geográfico Agustín Codazzi – IGAC, por haber facilitado los productos
provenientes de diferentes sensores remotos (imágenes satelitales), utilizadas en esta
investigación.
A los ingenieros Leonardo Alfonso y Fabio Amador de UNESCO – IHE en Holanda, por
amablemente haber facilitado los resultados de su investigación, adelantada en un sector
de la misma zona de estudio, dentro del marco del mismo proyecto COLCUENCAS.
Estudio a escala regional de los componentes hidrológicos e hidráulicos del complejo
cenagoso del Bajo Magdalena ente Calamar y Bocas de Ceniza con apoyo de
modelación matemática y SIG
A la Universidad del Norte, en Barranquilla, y allí al profesor Manuel Alvarado por sus
amables comentarios e información suministrada sobre la zona de estudio y a la
Universidad Nacional sede Medellín (UNALMED), por haber facilitado la herramienta
FLOODEM, utilizada en esta investigación y sugerida para adelantar estudios preliminares
de inundación a partir de parámetros geomorfológicos.
Resumen
El Bajo Magdalena, entre Calamar y Bocas de Ceniza, corresponde a una región con
importantes intercambios río – ciénagas, que ha sido reconocida por la convención
internacional Ramsar como zona estratégica de preservación mundial. Además, constituye
una fuente importante para el desarrollo regional, no sólo por la navegabilidad del río
Magdalena, sino por la sostenibilidad de las ciénagas circundantes y conectadas al río, las
cuales sirven de sustento importante para las comunidades asentadas en esta zona del
país. En la presente década, esta zona experimentó una gran afectación a causa del
fenómeno de La Niña (2010 – 2011), evidenciando el poco conocimiento sobre el
hidrosistema y la escasez de información para una adecuada gestión del mismo. En esta
investigación se estudian a nivel regional, y a partir de la información existente, los
componentes hidrológicos e hidráulicos de la zona entre Calamar y Bocas de Ceniza,
incluyendo su conexión con la Ciénaga Grande de Santa Marta (CGSM), complementando
la información a partir de productos derivados de sensores remotos. Esto ha permitido no
solo conocer preliminarmente cómo son los intercambios hídricos entre las ciénagas y el
río Magdalena, sino además identificar las zonas susceptibles a inundación, como un
aporte al entendimiento del sistema y a la adecuada gestión de este. Para ello, se ha
seguido una propuesta metodológica en aras de abordar la calibración de modelos
hidráulicos a partir de la información extraída de sensores remotos la cual se ha aplicado
en la zona de estudio. Los resultados obtenidos se consideran satisfactorios, calibrando el
parámetro sensible en la modelación hidráulica que casi siempre resulta ser el coeficiente
de rugosidad “n” de Manning. El mejor valor de este parámetro se encontró en un rango
entre 0.020 a 0.024 para el cauce principal, mientras que para la zona inundable se
obtuvieron valores del orden de 0.028 a 0.032. La evaluación del desempeño de la
modelación hidráulica se realizó a partir del cálculo del parámetro “F”, ampliamente
reportado en la bibliografía científica de inundaciones y con el cual es posible establecer
una tendencia adecuada de los valores del coeficiente de rugosidad en la zona analizada
a partir de sensores remotos. La zona con mayor susceptibilidad a la inundación resultó
ser la de Suan y Santo Tomás, al costado occidental, en jurisdicción de la Corporación
Estudio a escala regional de los componentes hidrológicos e hidráulicos del complejo
cenagoso del Bajo Magdalena ente Calamar y Bocas de Ceniza con apoyo de
modelación matemática y SIG
Autónoma Regional del Atlántico (CRA) y los municipios de El Piñón, Salamina y Remolino
en el costado oriental, en la jurisdicción de la Corporación Autónoma Regional del
Magdalena (CORPAMAG). Para la identificación de las zonas susceptibles a inundación,
además de la modelación matemática, se utilizó un criterio geomorfológico a partir del
análisis de las geoformas y dinámica fluvial existentes en la zona, con lo cual se logró
establecer el cauce histórico del río Magdalena y presentar una propuesta de delimitación
del complejo cenagoso del departamento del Atlántico entre el K45 y el K20, medido a
partir del cero en el Mar Caribe.
En cuanto a los promedios anuales de intercambios hídricos entre el río Magdalena y las
ciénagas circundantes en el período 2009 – 2012, estos muestran un comportamiento
alternado y/o cíclico, río – ciénaga y viceversa, con volúmenes máximos de (+168 Mm3)
hacia la CGSM, hacia la margen oriental. Hacia el costado occidental los mayores aportes
se dan hacia el complejo cenagoso del Atlántico, del orden de (+40 Mm3). Con respecto a
estos valores promedio de intercambio, las cuantificaciones aquí realizadas muestran que
para el primer evento La Niña (2010 – 2011), los intercambios hídricos fueron en dirección
río – ciénagas con valores mayores en comparación con el período seco (2009 – 2010)
que lo antecedió (+627 Mm3), donde la dirección de los aportes fue ciénaga – río (-12
Mm3). En comparación con el segundo período La Niña (2011-2012), cuando las ciénagas
ya estaban llenas, los volúmenes de intercambio disminuyeron su magnitud y
permanecieron en dirección río – ciénagas con valores menores al primer período La Niña
(+184 Mm3). También se demostró la capacidad de regulación de caudales que tiene todo
el complejo cenagoso del Bajo Magdalena, donde se registra una regulación cercana al
10% del caudal entrante a través del río Magdalena en épocas húmedas y de los cuales
el 6% lo regulan las ciénagas del costado oriental (departamento del Atlántico) y el 4% es
regulado por las ciénagas del costado occidental (departamento del Magdalena). Esta
regulación se presenta en las ciénagas conectadas o no directamente con el río
Magdalena. Es así, como esta investigación aporta al conocimiento de la dinámica hídrica
de este valioso hidrosistema con interacción río – ciénagas y contribuye a una mejor toma
de decisiones sobre el mismo.
Abstract
Palabras clave: ENOS, Factor de desempeño "F", FLOODEM, HEC RAS, HEC-GeoRAS,
HEC HMS, HEC-GeoHMS, Iber, inundación, MDT, modelación matemática en hidráulica,
modelación matemática en hidrología, Sensores remotos, SIG, Regionalización.
Estudio a escala regional de los componentes hidrológicos e hidráulicos del complejo
cenagoso del Bajo Magdalena ente Calamar y Bocas de Ceniza con apoyo de
modelación matemática y SIG
Abstract
The lower Magdalena, between Calamar and Bocas de Ceniza, corresponds to a region
with important river - swamp exchanges, which has been recognized by the international
Ramsar convention as a strategic zone of world preservation. In addition, it is an important
source for regional development, not only for the navigability of the Magdalena River, but
also for the sustainability of the surrounding marshes connected to the river, which serve
as important sustenance for the communities settled in this area of the country. In the
present decade, this area experienced a great affectation due to the phenomenon of La
Niña (2010 - 2011), evidencing the scarce knowledge about the hydrosystem and the
scarcity of information for its adequate management. In this research, the hydrological and
hydraulic components of the area between Calamar and Bocas de Ceniza, including its
connection with the Ciénaga Grande de Santa Marta (CGSM), are complemented with
existing information derived from remote sensors. This has allowed not only to know how
water exchanges occur between the marshes and the Magdalena River, but also to identify
areas susceptible to flooding, as a contribution to the understanding of the system and its
proper management. For this, a methodological proposal has been followed to approach
the calibration of hydraulic models using the information extracted from the remote sensors
that has been applied in the study area. The results obtained are considered satisfactory,
calibrating the sensitive parameter in the hydraulic modeling that almost always turns out
to be the roughness coefficient "n" of Manning. The best value of this parameter was found
in a range of 0.020 to 0.024 for the main channel, while for the flood zone the values ranged
between 0.028 and 0.032. The evaluation of the performance of the hydraulic modeling
was made from the calculation of the parameter "F", widely reported in the scientific
literature of floods and with which it is possible to establish an adequate trend of the values
Abstract
of the roughness coefficient in the area analyzed. The area with the highest susceptibility
to flooding was Suan and Santo Tomas, on the western side, in the jurisdiction of the
Corporación Autónoma del Atlántico (CRA) and the municipalities of El Piñón, Salamina
and Remolino on the east side, jurisdictión of the Corporación Autónoma del Magdalena
(CORPAMAG). For the identification of areas susceptible to flooding, in addition to
mathematical modeling, a geomorphological criterion was used based on the analysis of
geoforms and river dynamics in the area, which established the historic channel of the
Magdalena River and presented a proposal to delimit the Atlantic Swampy complex of the
coast between K45 and K20. Measured from zero in the Caribbean Sea.
As for the annual averages of water exchanges between the Magdalena River and the
surrounding marshes in the 2009-2012 period, these show an alternative and / or cyclical
behavior, river - swamp and vice versa, with maximum volumes of (+168 Mm3) to the
CGSM, on the eastern bank. Towards the western side, the largest contributions are to the
Atlantic marsh complex, on the order of (+40 Mm3). With respect to the mean values of
exchange, the quantifications performed here show that for the first La Niña event (2010-
2011), water exchanges went towards river swamps with higher values compared to the
dry period (2009-2010) (+627 Mm3), where the direction of the contributions was swamp –
river (- 12 Mm3). Compared with the second period of La Niña (2011-2012), when swamps
were already full, swamp volumes decreased in magnitude and remained in the direction
of marshes with values lower than the first period of La Niña (+184 Mm3). It was also
demonstrated the capacity of regulation of flows that has the entire marsh complex of the
low Magdalena, where about 10% of the incoming flow from the Magdalena River in wet
times and of which 6% is regulated by marshes. (Department of the Atlantic) and 4% is
regulated by the marshes on the western side (Department of Magdalena). This regulation
occurs in wetlands connected or not directly with the Magdalena River. Thus, this research
contributes to the knowledge of the water dynamics of this valuable hydrosystem with River
- swamp interactions and contributes to its better decision-making process.
Contenido
Pág.
1. Introducción ............................................................................................................. 1
1.1 Generalidades .....................................................................................................1
1.2 Objetivo general ..................................................................................................5
1.3 Objetivos específicos ..........................................................................................5
1.4 Pregunta de investigación ...................................................................................6
1.5 Organización del documento ...............................................................................6
2. El Bajo Magdalena.................................................................................................... 9
2.1 Resumen de la información primaria y secundaria disponible y recolectada en la
zona de estudio. ..........................................................................................................12
2.2 Climatología y geomorfología en la zona de estudio .........................................16
2.2.1 Climatología .................................................................................................... 16
2.2.2 Dinámica fluvial en la zona de estudio ............................................................. 27
2.3 Comentarios finales del Capítulo 2 ....................................................................34
Lista de figuras
Pág.
Figura 1-1: Localización general de la zona de estudio...............................................5
Figura 2-1: Áreas de drenaje aferentes al río Magdalena en los costados oriental y
occidental de la zona de estudio. Fuente: Elaboración propia. ........................................10
Figura 2-2: Zona de estudio y sus principales componentes (ríos y ciénagas) en el Bajo
Magdalena entre Calamar y Bocas de Ceniza. Se localizan aquí las principales
estaciones hidrológicas del hidrosistema. Fuente: Elaboración propia. ...........................12
Figura 2-3: Localización de 61 estaciones climatológicas en la zona de estudio ............13
Figura 2-4: Mosaico de Planchas IGAC en la zona de estudio georreferenciadas sobre
mapa de alturas sombreado. ...........................................................................................15
Figura 2-5: Precipitación total multianual estación Los Campanos departamento del
Atlántico. .........................................................................................................................17
Figura 2-6: Precipitación total multianual para nueve estaciones ubicadas en el costado
occidental departamento del Atlántico. ............................................................................17
Figura 2-7: Precipitación total multianual estación Monterrubio departamento del
Magdalena. .....................................................................................................................18
Figura 2-8: Precipitación total multianual 28 estaciones costado oriental departamento
del Magdalena. ................................................................................................................19
Figura 2-9: Precipitación promedio mensual multianual en estaciones representativas del
norte, centro y sur del área de estudio. a) Costado Oriental y b) Costado Occidental.
Fuente: Elaboración propia..............................................................................................20
Figura 2-10: Isoyetas de precipitación anual multianual en la zona de estudio. Fuente:
Elaboración propia. .........................................................................................................21
Figura 2-11: Precipitación total en la estación Ernesto Cortissoz (1986 – 2011) Vs.
anomalía ONI2. (fase neutra). Fuente: Elaboración propia. .............................................22
Figura 2-12: Precipitación mensual acumulada (trimestral) en la estación Ernesto
Cortissoz (1997 – 1998) Vs. anomalía ONI2. (fase neutra). Fuente: Elaboración propia. 23
Figura 2-13: Precipitación mensual acumulada (trimestral) en la estación Ernesto
Cortissoz (1997 – 1998) Vs. anomalía ONI2. (fenómeno El Niño). Fuente: Elaboración
propia. .............................................................................................................................24
Figura 2-14: Precipitación mensual acumulada (trimestral) en la estación Ernesto
Cortissoz (2010 – 2011) Vs. anomalía ONI2. (fenómeno La Niña). Fuente: Elaboración
propia. .............................................................................................................................24
Figura 2-15: Valores de temperatura promedio anual (1980 – 2011). Fuente: Elaboración
propia. .............................................................................................................................25
Estudio a escala regional de los componentes hidrológicos e hidráulicos del complejo
cenagoso del Bajo Magdalena ente Calamar y Bocas de Ceniza con apoyo de
modelación matemática y SIG
Figura 4-16: Hidrograma de caudal en la entrada y relación nivel contra tiempo aguas
abajo, para la modelación 1D. La línea con cuadros representa la fecha de la imagen del
22 de octubre de 2010 y la línea con triángulos la imagen del 28 de diciembre de 2011.
..................................................................................................................................... 111
Figura 4-17: Resultados: (a) Imagen 28/12/2011 vs resultado de modelación 1D
(nc y nll = 0.020) - medida de desempeño F (Aronica) = 56.08%., (b) Imagen 28/12/2011
vs resultado de modelación 2D (nc = 0.020, nll1 = 0.028, nll2 = 0.032) - medida de
desempeño F (Aronica) = 56.20%., (c) Imagen 28/12/2011 vs resultado de modelación
2D (nc = 0.036, nll1 = 0.036, nll2 = 0.036) - medida de desempeño F (Aronica) = 53.52%.
..................................................................................................................................... 115
Figura 4-18: Sección transversal típica para los caños Renegado, Aguas Negras y
Clarín. Fuente: Adaptado de (CORPAMAG-DNP-PROCIÉNAGA, 1993). ..................... 117
Figura 4-19: Curvas nivel vs almacenamiento en ciénagas con batimetría (costado
occidental). Fuente: Elaboración propia. ....................................................................... 117
Figura 4-20: a) MDT conformado con la información secundaría (representación global) y
topobatimetría del cauce principal y cuerpos de agua (ciénagas más representativas).
b) Representación en modelo hidrodinámico 1D (HEC RAS V4.1)................................ 118
Figura 4-21: Infograma de modelación 1D (HEC-RAS V4.1) en la zona general de
estudio. ......................................................................................................................... 120
Figura 4-22: Caudal de entrada Vs. niveles en la salida del área general de estudio. .. 121
Figura 4-23: Relación caudal de entrada (Calamar-Incora K7) Vs. niveles en la salida del
área general de estudio. ............................................................................................... 121
Figura 4-24: Análisis de frecuencia de valores extremos en la estación Calamar a partir
de análisis de la serie anual de excedencias (período 1940 – 2011)............................. 122
Figura 4-25: Análisis de frecuencia de valores extremos en la estación Incora K7 a partir
de valores medios máximos anuales (período 1972 – 2011)......................................... 123
Figura 4-26: Comparación resultado de inundación en el modelo 1D Vs imagen satelital
sensor RapidEye 28/12/2011. ....................................................................................... 126
Figura 4-27: Multitemporal de imágenes en la zona de estudio. ................................... 127
Figura 4-28: Dinámica fluvial del río Magdalena entre Calamar y Bocas de Ceniza, zona
inundable por geomorfología (elaboración propia). ....................................................... 128
Figura 5-1: Dinámica general de los flujos de agua entre el río Magdalena, las ciénagas
circundantes y los caños de interconexión en la zona de estudio.................................. 133
Figura 5-2: Perfil longitudinal y curva de remanso sobre el río Magdalena diciembre 21
de 2010 (fase La Niña).................................................................................................. 134
Figura 5-3: Visita a la zona de estudio. Diques y estructuras de control alrededor de las
ciénagas. Complejo cenagoso del Departamento del Atlántico. Sobre imagen Google
Earth. ............................................................................................................................ 135
Figura 5-4: Dinámica general del flujo entre el río Magdalena y los caños Renegado
(K54), Aguas negras (K39) y Clarín (K13) - Escenario fase La Niña del ENOS............. 136
Figura 5-5: Caudales netos (ΔQ=ƩQentrada – ƩQsalida) en la CGSM y el complejo cenagoso
Pajaral - Escenario fase La Niña del ENOS. ................................................................. 137
Contenido
Figura 5-6: Caudales netos de intercambio en las ciénagas Buenavista, Tiogollo, San
Antonio - Escenario fase La Niña del ENOS. .................................................................138
Figura 5-7: Caudales netos en las ciénagas costado occidental del río Magdalena -
Escenario fase La Niña del ENOS. ................................................................................139
Figura 5-8: Perfil longitudinal sobre el río Magdalena febrero 16 de 2010 (Fase El Niño).
......................................................................................................................................140
Figura 5-9: Dinámica general del flujo entre el río Magdalena y los caños Renegado
(K54), Aguas negras (K39) y Clarín (K13) - Escenario fase El Niño del ENOS. .............141
Figura 5-10: Caudales netos en la CGSM y el complejo cenagoso Pajaral - Escenario
fase El Niño del ENOS. .................................................................................................142
Figura 5-11: Caudales netos en las ciénagas Buenavista, Tiogollo, San Antonio -
Escenario fase El Niño del ENOS..................................................................................143
Figura 5-12: Caudales netos en las ciénagas costado occidental del río Magdalena -
Escenario fase El Niño del ENOS..................................................................................143
Figura 5-13: Volúmenes de intercambio hídrico en las ciénagas del Bajo Magdalena
período JJA/2009 – DEF/2012. ..................................................................................... 147
Figura 5-14: Volúmenes de intercambio hídrico en las ciénagas del Bajo Magdalena
período JJA/2009 – DEF/2012 y condiciones de frontera (entrada y salida). .................149
Figura 5-15: Volúmenes de intercambio hídrico en las ciénagas del Bajo Magdalena
período JJA/2009 – DEF/2012 y condiciones de frontera aguas abajo aumentada 1.0 m.
......................................................................................................................................152
Figura 5-16: Condiciones de frontera para escenario de inundación con períodos de
ocurrencia de 1 cada 50 años y 1 cada 100 años. Condiciones de frontera (entrada y
salida). .......................................................................................................................... 153
Figura 5-17: Inundación caudal Tr 100 años en el Bajo Magdalena. ............................ 155
Figura 5-18: Inundación caudal Tr 100 años en el Bajo Magdalena, sobre imagen Google
Earth. ............................................................................................................................ 156
Figura 5-19: Inundación caudal Tr 100 años en el Bajo Magdalena, Nivel de agua hasta
de 1.00 m en el Municipio de Suan (K100) en la margen occidental. Al costado oriental
desborde de la conexión entre ciénagas San Antonio (K100) y El Piñón (K85). b) Registro
fotográfico municipio de Suan (26/01/2011) en la ocurrencia de la Fase Niña del ENOS
en 2011. Parte baja de la zona de estudio. Fuente: IDEHA - Gobernación del Atlántico,
2011. ............................................................................................................................. 156
Figura 5-20: Inundación caudal Tr 100 años en el Bajo Magdalena, Nivel de agua hasta
de 1.30 m en el Municipio de El Piñón (K85) en la margen oriental. .............................. 157
Figura 5-21: Inundación caudal Tr 100 años en el Bajo Magdalena, Nivel de agua hasta
de 1.30 m en el Municipio de Salamina (K80) en la margen oriental (Departamento de
Magdalena). ..................................................................................................................157
Figura 5-22: Inundación caudal Tr 100 años en el Bajo Magdalena, riesgo de desborde
aguas arriba del Municipio de Remolino (K40) en la margen oriental (Departamento de
Magdalena). Ciénaga Paraíso en su máxima capacidad. ..............................................158
Figura 5-23: Inundación caudal Tr 100 años en el Bajo Magdalena, riesgo de desborde
aguas arriba del Municipio de Sitio Nuevo (K34) en la margen oriental (Departamento de
Magdalena). Ciénaga Paraíso en su máxima capacidad. ..............................................158
Estudio a escala regional de los componentes hidrológicos e hidráulicos del complejo
cenagoso del Bajo Magdalena ente Calamar y Bocas de Ceniza con apoyo de
modelación matemática y SIG
Figura 5-24: Inundación caudal Tr 100 años en el Bajo Magdalena, extensión del
complejo cenagoso del Atlántico entre el K40 - K20 las ciénagas: Paraíso, Santo Tomás,
Sabanagrande, El convento y Malambo. ....................................................................... 159
Figura 5-25: Inundación caudal Tr 50 años en el Bajo Magdalena. .............................. 160
Figura 5-26: Cota de lámina de agua en la zona de estudio. Tr 50 años. ..................... 161
Figura 6-1: Salinidad intersticial (0.5 m) promedio registrada en los suelos de manglar de
cinco estaciones de monitoreo durante el período 1995-2013. Las flechas rojas indican la
apertura de los caños Clarín (año 1996), Aguas Negras y Renegado (año 1998). Fuente:
INVEMAR (2013) .......................................................................................................... 168
Figura 6-2: Delimitación propuesta para la ciénaga Malambo (K20) a partir de criterio
hidráulico. Se identifican en color cian la extensión propuesta, 3.4 km a partir del río
Magdalena, en color rojo la afectación en la zona con una cota de agua 6.60 m.s.n.m.
Tr 50 y Tr 100 años. Predios afectados considerados dentro de la ciénaga.................. 170
Figura 6-3: Curva de duración de caudales (CDC) estación Calamar. Fuente:
Elaboración propia. ....................................................................................................... 177
Figura 6-4: Curva de duración de caudales (CDC) estación Incora K7. Fuente:
Elaboración propia. ....................................................................................................... 177
Figura 6-5: Resultado de los intercambios hídricos para el día 16 de febrero de 2010 en
la modelación matemática implementada. Período seco. .............................................. 180
Figura 6-6: Resultado de los intercambios hídricos para el día 21 de diciembre de 2010
en la modelación matemática implementada. Período húmedo (La Niña 1). ................. 181
Figura 6-7: Resultado de los intercambios hídricos para el día 28 de diciembre de 2011
en la modelación matemática implementada. Período húmedo (La Niña 2). ................. 182
Figura 6-8: Resultado mancha de inundación modelación matemática 1D (TR 100) +
inundación por geomorfología Vs Zona inundación periódica IDEAM (2011-2012). ...... 184
Figura 6-9: Resultado mancha de inundación modelación matemática 1D (TR 50) +
inundación por geomorfología Vs Zona inundación periódica IDEAM (2011-2012). ...... 185
Figura 6-10: Propuesta de delimitación de ciénagas del costado occidental a partir del
criterio conjunto geomorfología y modelación matemática TR 100 años. ...................... 186
Figura 6-11: Propuesta de reingeniería de la red hidrológica en el hidrosistema Bajo
Magdalena entre Calamar – Bocas de Ceniza. Inicialmente medición de niveles en
ciénagas principales y caudal en cauce principal del río Magdalena, caño Renegado,
caño Aguas Negras y caño Clarín. ................................................................................ 187
Contenido
Lista de tablas
Pág.
Tabla 2-1: Relación de planchas IGAC (Esc 1: 25.000) disponibles para la zona de
estudio ............................................................................................................................14
Tabla 2-2: Relación de principales imágenes satelitales utilizadas. .................................16
Tabla 3-1: Perspectivas de modelos hidrológicos considerados potencialmente
aplicables en esta investigación. .....................................................................................46
Tabla 3-2: Características de los modelos hidrológicos considerados potencialmente
aplicables en esta investigación. .....................................................................................47
Tabla 3-3: Panorama de herramientas computacionales (modelación hidráulica)
potencialmente aplicables en la zona de estudio. ............................................................58
Tabla 3-4: Resumen de productos globales para la representación de la tierra ..........66
Tabla 3-5: Sistemas de Información Geográfica populares, tabla adaptada de (Bedient,
Huber, et al., 2008). .........................................................................................................67
Tabla 3-6: Principales misiones y sensores para observaciones de la tierra. ..................73
Tabla 4-1: Principales parámetros morfométricos en las microcuencas del área de
estudio (Bajo Magdalena). ...............................................................................................93
Tabla 4-2: Reporte trimestral del comportamiento de fenómeno ENOS en el período 1981
– 2011, identificando fases Niño (color rojo), Niña (color azul) y Normal (color gris).
Fuente: NOAA. ..............................................................................................................100
Tabla 4-3: Comparación de resultados 1D contra imagen clasificada RapidEye
(28/12/2011). .................................................................................................................112
Tabla 4-4: Comparación de resultados 2D contra imagen clasificada RapidEye
(28/12/2011). .................................................................................................................113
Tabla 4-5: Caudales entrantes a la zona de estudio para diferentes períodos de retorno
en el río Magdalena. ......................................................................................................124
Tabla 4-6: Caudales en las cuencas del Bajo Magdalena para diferentes períodos de
retorno. .......................................................................................................................... 125
Tabla 4-7: Exigencia computacional de la modelación hidráulica implementada. ..........125
Tabla 5-1: Resumen volúmenes de intercambio en las 31 ciénagas del área de estudio.
En color rojo los volúmenes de salida de las ciénagas y viceversa en color azul. .........145
Tabla 5-2: Volúmenes de intercambio trimestral en las 31 ciénagas del área de estudio.
El color rojo indica la salida de volúmenes de agua de las ciénagas y el color azul la
entrada de volúmenes de agua a las ciénagas. ............................................................. 148
Estudio a escala regional de los componentes hidrológicos e hidráulicos del complejo
cenagoso del Bajo Magdalena ente Calamar y Bocas de Ceniza con apoyo de
modelación matemática y SIG
1.1 Generalidades
La PNGIRH reconoce que existe baja capacidad técnica en el ámbito municipal en gestión
del riesgo y manejo integral del recurso hídrico y de esta forma, parte de las acciones se
orientan al apoyo técnico a los municipios en el desarrollo de planes de mitigación. En las
directrices para el cumplimiento de este objetivo, la política en cuanto al riesgo prevé
“Generación y divulgación de información y conocimiento sobre riesgos que afecten la
oferta y disponibilidad hídrica”, siendo una línea de acción “Generar conocimiento sobre
riesgos asociados al recurso hídrico, mediante acciones de identificación y caracterización
de la vulnerabilidad de los ecosistemas clave para la regulación hídrica”, mientras que para
el fortalecimiento institucional la estrategia busca fomentar y desarrollar acciones de
investigación y de manejo para la información relacionada con el recurso hídrico (MAVDT,
2010). La Guía de POMCAs (MADS, 2013), incluye entre otros aspectos los análisis de
riesgo por inundación, y sin lugar a dudas, el conocimiento generado en esta tesis sobre
el hidrosistema objeto de estudio apoyará en este y otros aspectos la gestión del mismo..
Esta tesis, hace parte del componente de investigación del proyecto de acompañamiento
liderado por UNESCO – IHE y responde a la necesidad e importancia de mejorar el
entendimiento del comportamiento hidrológico e hidráulico regional del sistema cenagoso
ubicado en la parte baja del río Magdalena y que es de vital importancia, no solo para las
actividades de la región, sino como zona de amortiguamiento, con ecosistemas específicos
asociados. Este sistema cenagoso fue designado en 1998 dentro del convenio Ramsar, o
convención relativa a los humedales de importancia internacional y hace parte de la
jurisdicción de la Corporación Autónoma Regional del Atlántico (CRA) en el costado
occidental y de la Corporación Autónoma Regional del Magdalena (CORPAMAG) en el
costado oriental.
1
(MADS), Ministerio de Ambiente y Desarrollo Sostenible antes (MAVDT), Ministerio de Ambiente,
Vivienda y Desarrollo Territorial.
4 Introducción
El río Magdalena se ha dividido en los sectores geográficos, alto, medio y bajo, para su
caracterización. La parte baja, que es objeto del trabajo aquí reportado, va desde la
difluencia del Canal del Dique (municipio de Calamar) hasta su desembocadura en el Mar
Caribe (municipio de Barranquilla), en una longitud de aproximadamente 100 km y con un
área del orden de 8,000 km2 (incluyendo la Ciénaga Grande de Santa Marta y sus afluentes
desde la Sierra Nevada de Santa Marta). Esta área constituye una zona plana, poco
pendiente e inundable, con un complejo sistema de ríos y humedales y de interacciones
entre ellos. La planicie de inundación del Magdalena apoya diversos tipos de ecosistemas
asociados a la formación de humedales o ciénagas con espejos de agua o no 2; también
da sustento a las áreas agrícolas y pobladas, así como a importantes actividades
económicas en la región. El proyecto piloto de ajuste del POMCA, de acuerdo con la nueva
PNGIRH (realizado por la Universidad del Magdalena) se centró en las cuencas locales y
en los humedales ubicados sobre la orilla occidental del río Magdalena en el Departamento
del Atlántico, justo aguas arriba de la ciudad de Barranquilla (U. del Magdalena, 2011);
mientras que, a escala regional, en esta tesis, el costado oriental es también el objeto de
estudio, debido a las conexiones existentes entre el río Magdalena y la Ciénaga Grande
de Santa Marta (CGSM) (ver Figura 1-1).
2
Los ecosistemas con espejos de agua permanente están asociados a los ciclos biológicos de los
peces, aves migratorias, actividades de subsistencia y turismo de diferentes grupos sociales. Los
segundos, que no tienen espejo de agua permanente, son importantes por el tipo de vegetación en
donde prevalecen plantas herbáceas y en donde igualmente interrelaciona una fauna típica (aves,
insectos, reptiles y otros) (UNAL, 2004). Además de la pesca existe una gran cantidad de especies
vegetales acuáticas, de pantano y terrestres (CORMAGDALENA, 2007). Existe gran cantidad de
especies vegetales características del espejo de agua como la lechuga de agua, la taruya o batata
de agua y la hoja de raya. En las zonas próximas al litoral las ciénagas están acompañadas de
manglar.
Introducción 5
▪ Consolidar en la zona de estudio, y con ayuda de una herramienta SIG, una base
de datos de información primaria y secundaria, complementada con información
6 Introducción
¿Cuál es el método más adecuado que permite a nivel regional estimar intercambios –río-
ciénaga en una zona baja con información limitada y cómo esta metodología puede apoyar
la definición de rondas hídricas en un complejo cenagoso de escala regional?
ocurrencia del fenómeno macroclimático ENOS3 y sus dos fases: La Niña y el Niño. En el
Capítulo 3, se describe el estado del arte y las características de los modelos matemáticos
hidrológicos e hidráulicos existentes y que se consideran aplicables en la zona de estudio,
además, se especifica la importancia de las fuentes de información secundarias a partir de
sensores remotos y Sistemas de información Geográfica (SIG) y cómo estas apoyan
notablemente la investigación de los recursos hidráulicos; seguidamente se discute el
papel de la geomorfología y más específicamente de la dinámica fluvial en la delimitación
de zonas susceptibles a inundación y se particulariza la metodología utilizada para el
estudio, destacando el desarrollo de una propuesta metodológica para la calibración de
modelos matemáticos en hidráulica a partir de la información aportada por sensores
remotos, entre otros aspectos. En el Capítulo 4, se presenta la aplicación de la
metodología discutida en el Capítulo 3, enfatizando en la hidrología, la hidráulica y la
dinámica fluvial de la zona de estudio. Mientras que el resultado de los intercambios
hídricos río Magdalena – ciénagas se presenta y discute en el Capítulo 5 a través de
distintos escenarios de análisis, un año normal, año Niño y año Niña. En este Capítulo 5,
se presentan también los resultados de las zonas susceptibles a inundación, aplicando los
criterios geomorfológicos e hidráulicos descritos en el Capítulo 3. Finalmente, en el
Capítulo 6 de esta investigación se incluyen las conclusiones generales y específicas,
además de las recomendaciones a las diferentes entidades del orden local y nacional que
tienen injerencia directa o no, en la zona de estudio (el Bajo Magdalena entre Calamar y
Bocas de Ceniza) y cómo a través de nuevas investigaciones en la zona se pueden
complementar los resultados aquí reportados, para tener cada día mejores herramientas
de gestión de hidrosistemas complejos, como el analizado en esta investigación.
3
ENOS: El Niño Oscilación del Sur.
2. El Bajo Magdalena
En este Capítulo se describe la zona de estudio, que corresponde al Bajo Magdalena entre
Calamar y Bocas de Ceniza, sus características socioeconómicas principales y la
problemática existente con relación a las inundaciones. Además, se detalla la información
primaria y secundaria que fue posible conseguir de la zona; información principalmente de
tipo cartográfico, climatológico, hidrológico y geomorfológico.
El río Magdalena es quizás la fuente hídrica más importante del país y comúnmente suele
dividirse en tres grandes regiones con características geomorfológicas, hidrológicas e
hidráulica distintas; en la zona del Alto Magdalena el río nace a 3.865 msnm en el páramo
de las Papas y finaliza en el Municipio de Honda (Tolima), en donde comienza el
Magdalena Medio que se extiende hasta el Municipio del Banco. A partir de este municipio,
en el Departamento del Magdalena y hasta Bocas de Ceniza (Departamento del Atlántico)
el cauce se conoce como el Bajo Magdalena. En este último tramo el río desciende
aproximadamente 33 m en una longitud de 400 km, con un área de drenaje aproximada
de 105.250 km2 y un caudal promedio de 7100 m3/s a la altura del municipio de Calamar
(Bolívar) (CORMAGDALENA, 2007). Dentro del Bajo Magdalena se encuentra enmarcada
el área particular de estudio de esta investigación, la cual inicia en cercanías del municipio
de Calamar, inmediatamente aguas abajo de la derivación del Canal del Dique, en donde
se deriva cerca del 7.5% de los caudales que transporta el río Magdalena a esta altura
(CORMAGDALENA-UNAL, 2007) y finaliza al norte, en Bocas de Ceniza en Barranquilla,
en donde se encuentra la desembocadura del río Magdalena en el Mar Caribe (ver Figura
2-1). Al costado occidental la zona de estudio limita con el Departamento del Atlántico,
donde se encuentra la zona de humedales del río Magdalena, la cual tiene una superficie
aproximada de 124,676 Ha y se encuentra subdividida en trece subcuencas que drenan
sus aguas al río Magdalena (U. del Magdalena, 2011). Ver Figura 2-1.
10 Estudio a escala regional de los componentes hidrológicos e hidráulicos del complejo
cenagoso del Bajo Magdalena ente Calamar y Bocas de Ceniza con apoyo de
modelación matemática y SIG
Figura 2-1: Áreas de drenaje aferentes al río Magdalena en los costados oriental y
occidental de la zona de estudio. Fuente: Elaboración propia.
De otra parte, el límite oriental de la zona de estudio, como se aprecia en la Figura 2-1, lo
constituye una vasta zona de descanso del río Magdalena. Esta es una zona muy compleja
en donde se encuentran múltiples madreviejas o meandros abandonados que han formado
ciénagas, algunas de las cuales se encuentran interconectadas entre sí o con el río
Magdalena. Debido a esta característica, el costado occidental del área de estudio resultó
ser la franja más compleja para establecer los límites de la zona de investigación asociada
con el río Magdalena. Sin embargo, después de varias visitas técnicas y visitas de campo
a la región (Barranquilla y Santa Marta), revisión de distintos documentos e informes
técnicos de proyectos elaborados en esta área por entidades como INVEMAR,
CORPAMAG, distintas Universidades e investigadores tales como: (Díaz-Granados,
Camacho, et al., 1996; CORMAGDALENA-UNINORTE, 2007), etc., se pudo verificar que
El Bajo Magdalena 11
existen varios caños de derivación del río Magdalena, que lo conectan con la Ciénaga
Grande de Santa Marta (CGSM), ver Figura 2-2, y por tanto el área de estudio debe incluir
esta última. Por su tamaño, son tres los caños más representativos: a) Caño Renegado,
ubicado aproximadamente en el K54 del río Magdalena, en las inmediaciones del
corregimiento del Guaimaro (Magdalena), tiene una longitud aproximada de 6 km hasta la
ciénaga Buenavista. b) Caño Aguas Negras, ubicado aproximadamente en el K39 del río
Magdalena, entre los municipios de Remolino y Sitionuevo (Magdalena), tiene una longitud
aproximada de 18.6 km en el brazo más largo que llega al complejo cenagoso Pajaral.
c) Caño Clarín, es la derivación más aguas abajo de todos los caños, se encuentra ubicado
aproximadamente en el K13 del río Magdalena, aguas abajo del Puente Pumarejo que
cruza el río Magdalena para unir los departamentos del Magdalena y Atlántico; tiene una
longitud de 33 km hasta la CGSM. Cabe resaltar que el cero en el abscisado del río
Magdalena se encuentra en la desembocadura de éste en el Mar Caribe.
Estas conexiones antrópicas, fueron creadas entre los años 80 y 90 con varios objetivos,
entre ellos, llegar con agua dulce a la Ciénaga Grande de Santa Marta (CGSM) y así evitar
la salinización de este hidrosistema y la amenaza a los ecosistemas naturales de la CGSM,
además de servir al riego de fincas y otros. (Albarracín et.al,1995; Múnera et al., 2005). Lo
anterior, obligó lógicamente a considerar un área mayor como área de estudio en
comparación con la inicialmente planteada en la propuesta de investigación, para incluir la
CGSM y los aportes principales que a ésta llegan desde la Sierra Nevada de Santa Marta,
a través de los ríos Fundación y Aracataca principalmente, con el propósito de representar
lo mejor posible los intercambios hídricos de todo el hidrosistema.
MAR CARIBE
CIÉNAGAS DEL
COMPLEJO DE
CGSM
ATLÁNTICO
COMPLEJO
CENAGOSO
PAJARAL
COMPLEJO
CENAGOSO
UVERO C. BUENAVISTA
Figura 2-2: Zona de estudio y sus principales componentes (ríos y ciénagas) en el Bajo
Magdalena entre Calamar y Bocas de Ceniza. Se localizan aquí las principales estaciones
hidrológicas del hidrosistema. Fuente: Elaboración propia.
Pumarejo y el K100, en la derivación del Canal del Dique). Esta última información fue
levantada en campo por la Universidad Nacional de Colombia sede Bogotá para
CORMAGDALENA en el año 2012 y se caracteriza por tener secciones transversales del
río cada 200 metros. De otra parte, fue posible obtener la información de batimetrías
realizadas en las ciénagas: Sabanagrande, Santo Tomás, El Espinal y Palmar de Varela,
al costado occidental; esta información fue suministrada amablemente por la CRA 4 (U. del
Magdalena, 2011). Se adquirió también la cartografía oficial de la zona, comprendida por
planchas IGAC5, las cuales datan del año 1956 (ver Tabla 2- 2).
Las planchas IGAC son importantes como referencia inicial para los análisis de dinámica
fluvial y están a escala 1:25000, mientras que el levantamiento batimétrico en la zona se
utilizará para consolidar junto con información secundaria un modelo digital de terreno
(MDT) que permita representar lo mejor posible la zona de estudio. Información primaria
adicional no fue encontrada en la zona, a pesar de las visitas a entidades como
CORPAMAG, INVEMAR, etc.
Tabla 2-1: Relación de planchas IGAC (Esc 1: 25.000) disponibles para la zona de estudio
4
Corporación Regional del Atlántico.
5
Instituto Geográfico Agustín Codazzi
El Bajo Magdalena 15
2.2.1 Climatología
3500
3000
Precipitación (mm)
2500
2000
1500
1000
500
0
2000
2004
2008
1980
1982
1984
1986
1988
1990
1992
1994
1996
1998
2002
2006
2010
Año
29040290_CAMPANOS_LOS PROMEDIO
Figura 2-5: Precipitación total multianual estación Los Campanos departamento del
Atlántico.
2000
1500
1000
500
0
2002
1980
1982
1984
1986
1988
1990
1992
1994
1996
1998
2000
2004
2006
2008
2010
Año
29030140_SAN_JOSE 29035070_REPELON
29035080_NORMAL_MANATI 29035120_LIMON_EL
29040020_MONTEBELLO 29040070_PONEDERA
29040080_POLONUEVO' 29040190_SABANALARGA
29040200_LENA PROMEDIO
Figura 2-6: Precipitación total multianual para nueve estaciones ubicadas en el costado
occidental departamento del Atlántico.
Hacia el costado oriental, departamento del Magdalena, la estación con mayor registro de
precipitación total anual es la Estación Monterrubio, con un valor máximo cercano a los
2300 mm/año, registrado para el año 2010, mientras que la precipitación total media de
esta zona (costado oriental) se estima en 1238 mm/año, sin considerar los aportes
18 Estudio a escala regional de los componentes hidrológicos e hidráulicos del complejo
cenagoso del Bajo Magdalena ente Calamar y Bocas de Ceniza con apoyo de
modelación matemática y SIG
provenientes de la Sierra Nevada de Santa Marta (ver Figura 2-7 y Figura 2-8). Se aprecia
que la precipitación media anual en ambos costados del río Magdalena, sin contar la Sierra
Nevada (SN), es similar, con un promedio de 1150 mm/año y se destaca que se tiene una
máxima altura sobre el nivel del mar aproximada de 100 msnm, siendo la Sierra Nevada
(SN) un caso especial en el Bajo Magdalena entre Calamar y Bocas de Ceniza, ya que
corresponde a una altura máxima de 5678 msnm y cuya climatología se puede considerar
diferente a la del resto de la zona de estudio por sus distintos pisos térmicos, vegetación,
etc. Según (Mesa, 1997), la precipitación media en la SN es del orden de 1500 mm/año y
cae principalmente de mayo a octubre con dos máximos separados en mayo - junio y
septiembre-octubre; superior al promedio de 1150 mm/año estimado en esta investigación
para la zona baja sin contar la SN.
3000
2500
2000
1500
1000
500
0
1998
1980
1982
1984
1986
1988
1990
1992
1994
1996
2000
2002
2004
2006
2008
2010
Año
28040100_MONTERRUBIO PROMEDIO
3500
3000
2500
2000
1500
1000
500
0
1985
1986
1987
2009
2010
1980
1981
1982
1983
1984
1988
1989
1990
1991
1992
1993
1994
1995
1996
1997
1998
1999
2000
2001
2002
2003
2004
2005
2006
2007
2008
Año
29065080_MEDIA_LUNA PROMEDIO
Figura 2-8: Precipitación total multianual 28 estaciones costado oriental departamento del
Magdalena.
En la Figura 2-8 se aprecia que los valores máximos de precipitación son reportados en el
año 2010, coincidentes con la época de marcadas inundaciones. Se reportan valores
excepcionalmente muy por encima del promedio en la estación La Unión (cód. IDEAM
29060240) cerca de la CGSM, la cual, solo para el mes de octubre de 2010 reportó una
precipitación mensual de 1209 mm/mes, similar al promedio anual de la zona.
448 mm/trimestre, cerca del 37% del total de precipitación anual en la zona; mientras que
los menores valores de precipitación se presentan los primeros meses del año, en el
trimestre enero, febrero y marzo (EFM), con un valor medio de 36 mm/trimestre, solo el 3%
de la precipitación anual.
29045020_APTO_E_CORTISSOZ 15015020_LAYE
(1980 - 2010) (1980 - 2010)
200 200
PPT (mm)
PPT (mm)
150 150
100 100
50 50
0 0
FEB
SEP
NOV
ENE
MAR
MAY
DIC
OCT
JUL
JUN
ABR
AGO
FEB
SEP
NOV
ENE
MAY
MAR
OCT
DIC
JUL
JUN
AGO
ABR
Mes Mes
a. b.
29040290_CAMPANOS_LOS 29060350_BAYANO
(1980 - 2010) (1980 - 2010)
250 300
PPT acumulada (mm)
200 250
200
150
150
100 100
50 50
0 0
FEB
SEP
NOV
ENE
MAY
MAR
OCT
DIC
JUL
JUN
AGO
ABR
FEB
SEP
NOV
MAY
ENE
MAR
JUN
OCT
JUL
DIC
ABR
AGO
Mes Mes
a. b.
29040250_CAMPO_DE_LA_CRUZ 29040020_MONTEBELLO
(1980 - 2010) (1980 - 2010)
200 200
PPT (mm)
PPT (mm)
150 150
100 100
50 50
0 0
FEB
SEP
FEB
SEP
NOV
NOV
MAY
ENE
ENE
MAY
MAR
MAR
OCT
OCT
JUN
JUL
DIC
JUL
DIC
JUN
ABR
AGO
ABR
AGO
Mes Mes
a. b.
Figura 2-9: Precipitación promedio mensual multianual en estaciones representativas del
norte, centro y sur del área de estudio. a) Costado Oriental y b) Costado Occidental.
Fuente: Elaboración propia.
El Bajo Magdalena 21
Una comparación entre la precipitación y el índice ONI26 (por sus siglas en inglés)
particulariza cómo esta variable se ve notablemente afectada por la ocurrencia de
fenómenos macroclimáticos como el ENOS, en sus fases Niño y Niña. Se tomó como
referencia la estación Ernesto Cortissoz ubicada en el área metropolitana de la ciudad de
Barranquilla, que hace parte del área de estudio y se compararon los valores de
precipitación acumulada trimestralmente versus el índice ONI2. Se realiza la comparación
inicialmente para la totalidad de la serie de datos en el período (1986 – 2011), Figura 2-11,
luego se efectuó una comparación para un año normal (1989 -1990), Figura 2-12, año Niño
(1997 - 1998), Figura 2-13, y año Niña (2010 – 2011), Figura 2-14; los resultados para las
distintas fases se presentan a continuación y pueden ser considerados representativos, en
orden de magnitud, para toda la zona de estudio, debido a la enorme correlación existente
en el comportamiento de esta variable en ambos costados del río, como fue discutido
anteriormente.
Figura 2-11: Precipitación total en la estación Ernesto Cortissoz (1986 – 2011) Vs.
anomalía ONI2. (fase neutra). Fuente: Elaboración propia.
6
ONI2: Índice Oceánico del Niño – Corrección 2012. Valores menores a -0.5 oC indican fenómeno
de la Niña y valores mayores 0.5 oC indican fenómeno del Niño. Promedio de temperatura en el
Océano Pacífico.
El Bajo Magdalena 23
• Fase Neutra: Ver Figura 2-11, donde se detallan los valores de referencia o
normales en la estación, registrando un máximo de precipitación acumulada en el
trimestre, agosto, septiembre, octubre, (ASO), con un valor cercano a los
570 mm/trimestre, mientras que el trimestre con menores precipitaciones es enero,
febrero, marzo, (EFM), con valores cercanos a 2 mm/ trimestre, es decir,
precipitaciones nulas teniendo en cuenta la precisión de los aparatos de medida,
en su gran mayoría pluviómetros.
Anomalia ONI2
500 0.2
400
0
300
-0.2
200
100 -0.4
0 -0.6
1989 1989 1989 1989 1989 1989 1989 1989 1990 1990
Precipitación ONI2
2.5
Anomalia ONI2
400
2
300
1.5
200
1
100 0.5
0 0
1997 1997 1997 1997 1997 1997 1997 1997 1998 1998
Precipitación ONI2
700
Anomalia ONI2
600 -0.5
500
400 -1
300
200 -1.5
100
0 -2
2010 2010 2010 2010 2010 2010 2010 2010 2011 2011
Precipitación ONI2
Junto con la precipitación, la temperatura es una de las variables más importantes para
describir la climatología. La zona de estudio cuenta con 13 estaciones que reportan esta
variable y se caracteriza por registrar mínimos del orden de 24 oC, máximos que pueden
superar incluso los 38 oC y valores promedio de 28 oC, aproximadamente. Esto se resume
en la Figura 2-15 y Figura 2-16, para el período 1980 – 2011, donde se puede apreciar que
la temperatura media anual varía en +/- 10C. Según el método Caldas para clasificación de
clima (U. del Magdalena, 2011), las zonas con elevaciones menores a 1000 msnm y
temperaturas mayores o iguales a 24 oC son clasificadas como clima Cálido, como es el
caso de la zona de estudio; esto trae consigo que los tipos de cobertura vegetal, entre
otros aspectos morfológicos, tengan características especiales, razón por la cual en estos
análisis no ha sido incluida estrictamente la Sierra Nevada de Santa Marta, por tener pisos
térmicos diferentes.
28.50
Temperatura (0C)
28.00
27.50
27.00
1980
1994
2010
1982
1984
1986
1988
1990
1992
1996
1998
2000
2002
2004
2006
2008
Año
Figura 2-15: Valores de temperatura promedio anual (1980 – 2011). Fuente: Elaboración
propia.
26 Estudio a escala regional de los componentes hidrológicos e hidráulicos del complejo
cenagoso del Bajo Magdalena ente Calamar y Bocas de Ceniza con apoyo de
modelación matemática y SIG
7
Era Cenozoica, Período cuaternario, Época Holoceno/pleistoceno.
28 Estudio a escala regional de los componentes hidrológicos e hidráulicos del complejo
cenagoso del Bajo Magdalena ente Calamar y Bocas de Ceniza con apoyo de
modelación matemática y SIG
ser fluvial y lacustre, con la presencia de importantes barras, algunas ya estabilizadas que
forman islas y diques naturales, que mantienen tramos del río invariables en el tiempo,
terrazas intermedias y bajas que originan un valle de inundación pronunciado con
importante presencia de cuerpos cenagosos interconectados y meandros abandonados,
con evidencias de paleocanales del río Magdalena, principalmente hacia el costado
oriental. En esta geomorfología actual, resalta la presencia de importantes ejes viales y
diques (intervención antrópica), los cuales, en la mayoría de los casos, alteran y cambian
el funcionamiento natural del hidrosistema, limitando la conexión río – ciénagas a pasos
controlados, y la extensión de estas últimas.
CGSM
.
El Bajo Magdalena 29
otro elemento morfológico destacable es la isla fluvial (isla Rondón) ubicada en el K20,
mientras que, al costado oriental de la misma se encuadra un dique aluvial que controla el
Magdalena en el sector. Este dique aluvial se hace evidente al comparar el cauce activo al
año 2011 (color azul), con el cauce identificado para el año 1976 (color rojo), este último a
partir de la digitalización de la plancha IGAC (C_171vB). Hacia el costado occidental se
encuentran ciénagas conectadas al río Magdalena y entre sí, estas constituyen el complejo
cenagoso del Departamento del Atlántico; en la figura están las ciénagas Sabanagrande y
Malambo, más al norte. En la Figura 2-20 se presenta la franja de divagación del río
Magdalena a través de un análisis multitemporal de imágenes de distintos sensores
remotos, la cual está sobre las planchas IGAC de 1976. Se aprecia para los años 1976,
1990, 2002, 2005 y 2011, los cauces activos del río Magdalena junto a los espejos de agua
en las ciénagas. Es notable la dinámica del río Magdalena en la zona, donde los espejos
de agua en las ciénagas también han variado; además, se identifica un paleocauce del río
Magdalena más al oriente del cauce activo al año 2011, esta geoforma representa un
antiguo cauce del río Magdalena (Vargas Cuervo, 2015), con lo cual queda claro que toda
esa zona perteneció al río Magdalena. Otra zona del río se muestra en la Figura 2-21, esta
corresponde al (K50), entre los municipios de Ponedera (Atlántico) y Remolino
(Magdalena). Nótese la dinámica de las orillas en este sector al comparar la cartografía
análoga del IGAC, con fecha de 1956, con imágenes satelitales del sensor Landsat TM,
más recientes, años 1998 y 2002. Finalmente, como resultado de estos análisis se obtiene
el cauce histórico de río Magdalena entre Calamar y Bocas de Ceniza a partir de un criterio
geomorfológico, el cual corresponde a la vega del río Magdalena.
El Bajo Magdalena 31
Figura 2-19: Geomorfología fluvial en la zona de estudio entre Calamar - Bocas de Ceniza.
Fuente: Elaboración propia.
32 Estudio a escala regional de los componentes hidrológicos e hidráulicos del complejo
cenagoso del Bajo Magdalena ente Calamar y Bocas de Ceniza con apoyo de
modelación matemática y SIG
Figura 2-20: Dinámica fluvial río Magdalena entre Calamar – Bocas de Ceniza.
Municipios de Sabanagrande y Soledad. Fuente: Elaboración propia.
El Bajo Magdalena 33
1956
Figura 2-21: Dinámica fluvial río Magdalena río Magdalena entre Calamar – Bocas de
Ceniza. Municipio de Ponedera. Fuente: Elaboración propia.
34 Estudio a escala regional de los componentes hidrológicos e hidráulicos del complejo
cenagoso del Bajo Magdalena ente Calamar y Bocas de Ceniza con apoyo de
modelación matemática y SIG
río Magdalena y la Ciénaga Grande de Santa Marta (CGSM), ver Figura 2-22. Todos estos
aspectos, zona de regulación de caudales, interacción río – Ciénagas, etc., dan cuenta de
la importancia en la zona del sistema río Magdalena – ciénagas.
Hacia finales de 2010 con la denominada ola invernal que afectó todo el país y esta zona
en particular, quedó en evidencia la limitada planeación sobre el hidrosistema y la
vulnerabilidad de este ante escenarios extremos, con las inundaciones más críticas de la
historia reciente en esta zona del país, que ocasionaron incalculables pérdidas económicas
y efectos sociales. Ver Figura 2-23.
Distintos investigadores (Poveda, 2004b; Poveda, Gil, et al., 1998; Múnera, Vélez, et al.,
2003a), han demostrado que la hidroclimatología del país se ve afectada, entre otros, por
la presencia de fenómenos macroclimáticos como el ENOS en sus fases El Niño y La Niña,
lo cual, ha quedado en evidencia para la zona de estudio, a través de los análisis detallados
en este capítulo y basados en los datos de precipitación recibidos. Se concluye por tanto
que cualquier estudio de intercambio hídrico río Magdalena – ciénagas y viceversa, debe
contemplar por lo menos estos tres escenarios de funcionamiento del hidrosistema (seco,
húmedo, normal), como fue comentado anteriormente, en aras de poder disminuir la
vulnerabilidad de las comunidades asentadas en la zona, ya que la carencia de
conocimiento en cuanto a volúmenes de intercambio hídrico interanual e intranual y sus
efectos, limita la adecuada planeación de los entes territoriales y autoridades regionales
encargadas de la gestión del recurso hídrico en la zona. Teniendo en cuenta lo anterior, la
multiplicidad de soluciones individuales ante problemáticas como inundaciones (jarillones,
compuertas, diques, etc.) y la carencia de herramientas adecuadas que permitan una mejor
gestión y planeación de acciones sobre el hidrosistema regional del Bajo Magdalena, esta
investigación pretende contribuir al entendimiento hídrico de la macroregión, a través del
conocimiento de los intercambios hídricos río Magdalena – Ciénagas, y viceversa,
aplicando modelación matemática y otras herramientas.
Distintos estudios han sido realizados de forma aislada en la zona de estudio, lo que no ha
permitido una adecuada articulación de toda la información topográfica, y de
funcionamiento de todo el macro sistema. Haciendo un recorrido iniciando en el costado
occidental de la zona de estudio, en la jurisdicción de la CRA, allí se destacan los estudios
realizados en los POMCAS (Planes de Ordenamiento y Manejo de Cuencas) del año 2007
y 2011 (U. del Magdalena, 2011; CRA, 2007). En estos documentos se detallan
importantes aspectos hidrológicos, geomorfológicos e hidráulicos desde la ciénaga
Sabanagrande al norte, hasta el complejo cenagoso del Uvero, más al sur del
departamento del Atlántico y se incluyen además descripciones sociales y las interacciones
de la comunidad con los recursos hídricos de la zona. También en el departamento del
Atlántico, (Amador, 2013) analiza, a partir de la modelación matemática en hidráulica, la
generación de mapas probabilísticos de inundación en el sector comprendido entre las
ciénagas El Convento y Paraíso, entre los municipios de Ponedera y Palmar de Varela,
utilizando el modelo hidrodinámico SOBEK (1D-2D), el cual resuelve la hidrodinámica del
cauce principal en 1D y de la zona de inundación en 2D. La evaluación de la modelación
efectuada se realizó a partir de la comparación de la mancha de inundación simulada, con
el resultado de la mancha de inundación extraída de la imagen de sensor remoto Landsat
para el año 2007, mostrando en sus resultados el comportamiento hidráulico del sector
38 Estudio a escala regional de los componentes hidrológicos e hidráulicos del complejo
cenagoso del Bajo Magdalena entre Calamar y Bocas de Ceniza con apoyo de
modelación matemática y SIG
Una descripción detallada de la historia de los caños de conexión entre el río Magdalena y
la CGSM es presentada por (Múnera, Vélez, et al., 2003b), quienes describen el análisis
de la CGSM utilizando datos de sensores remotos, incluyendo los datos del sensor
TOPEX-POSEIDON, para calibración del modelo implementado aplicando un balance
hídrico de largo plazo para la estimación de los intercambios hídricos . También (Munera,
Vélez, et al., 2003a), detallan en este estudio la aplicación del balance hídrico en la zona
en la cual se considera cada ciénaga como un tanque donde se cumple el principio de
conservación de masa, siendo las variables precipitación y evapotranspiración
consideradas como principales, mientras que la infiltración, en suelos arcillosos que forman
diques y basines, como una variable secundaria con magnitudes dos a tres veces menores
que las de las variables principales. De esta manera los intercambios con el subsuelo
Estado del arte y propuesta metodológica 39
Una mención especial sobre estudios realizados al complejo lagunar CGSM, merece el
proyecto Colombo – Alemán pro-ciénaga 1993, el cual tuvo como objetivo la recuperación
ambiental de la ecorregión CGSM, la calidad de las aguas, la reforestación del mangle y
la reactivación de los flujos de agua entre el rio Magdalena y la CGSM; a este respecto se
destacan las dimensiones de las estructuras proyectadas, la sección trasversal y la
pendiente de los canales de conexión entre el río Magdalena y el complejo lagunas Pajaral
y la CGSM (CORPAMAG-DNP-PROCIÉNAGA, 1993).
8
Laboratorio de Ensayos Hidráulicos de Las Flores (Barranquilla).
40 Estudio a escala regional de los componentes hidrológicos e hidráulicos del complejo
cenagoso del Bajo Magdalena entre Calamar y Bocas de Ceniza con apoyo de
modelación matemática y SIG
Figura 3-1: Esquema de flujos de entrada y salida de una ciénaga. Fuente: Adaptado de
(CORMAGDALENA-UNAL, 2007).
Los modelos de caja negra hacen referencia a aquellos modelos en los que no se busca
detallar el comportamiento del fenómeno sino reproducir los valores de entrada y salida
del sistema: las redes neuronales son una representación de este tipo de modelos. Estos
modelos no son aplicables en aquellos casos en los que es necesario el estudio de la
variabilidad de clima, dado que se requiere información estructurada de los procesos. En
cambio, se requerirán en estos casos los modelos conceptuales que aplican conceptos
simplificados de los procesos físicos de generación de escorrentía (OMM, 2011)
44 Estudio a escala regional de los componentes hidrológicos e hidráulicos del complejo
cenagoso del Bajo Magdalena entre Calamar y Bocas de Ceniza con apoyo de
modelación matemática y SIG
Existen muchas herramientas computacionales que pueden ser utilizadas para realizar
modelación matemática en hidrología. (Beven, 2011), propone la siguiente metodología
para la selección del modelo, basada esencialmente en consideraciones de la posible
estructura del mismo:
1. Realizar una lista de los modelos en consideración. Esta lista puede incluir dos
partes: aquellos modelos que son fácilmente aplicables y aquellos que, aunque
complejos, deben considerarse en un proyecto si la inversión en tiempo y dinero no
son una restricción.
2. Preparar una lista de las variables previstas y requeridas por cada modelo. Decidir
si el modelo bajo consideración producirá las salidas requeridas para satisfacer el
objetivo particular de un proyecto.
3. Preparar una lista con las suposiciones hechas por el modelo. ¿Son las
suposiciones probablemente limitantes en términos del conocimiento que se tiene
acerca de la respuesta de la cuenca?
Estado del arte y propuesta metodológica 45
4. Hacer una lista de las entradas requeridas por el modelo para especificar el flujo
dominante, las fronteras y condiciones iniciales, y para la especificar los valores de
los parámetros. Decidir si toda la información requerida puede ser provista dentro
de las limitaciones de tiempo y costo del proyecto.
5. Determinar si se tiene el posible modelo en la lista. Si no, revisar los tres pasos
anteriores, relacionados con los criterios usados.
En la Tabla 3-1 se presenta, a manera de resumen, una breve descripción de los modelos
hidrológicos considerados aplicables en el presente estudio, destacando principalmente
los modelos de uso libre. En la Tabla 3-2 se presentan las características principales de
los modelos listados en la Tabla 3-1.Entre las características evaluadas se encuentran, por
ejemplo, el tipo de modelo hidrológico (determinístico o estocástico) y si tienen alguna
herramienta para la calibración de los parámetros del modelo, entre otros factores.
46
Tabla 3-1:
IHACRES significa Identificación de unidad Hidrográfica y componentes de flujo a partir de lluvias, Evaporación y datos de corrientes.
TOPMODEL es un modelo lluvia - escorrentía que basa sus predicciones distribuidas en un análisis de la topografía de la cuenca. El
modelo predice el exceso de saturación y el exceso de infiltración de la escorrentía superficial y subsuperficial. Desde el primer
TOPMODEL artículo que se publicó sobre el modelo en 1979 ha habido muchas versiones diferentes. La idea siempre ha sido que el modelo debe (K. Beven et al) , 1995; (K. Beven, 2011)
ser lo suficientemente simple para ser modificado por el usuario, de manera que las predicciones se ajusten, en la medida de lo
posible a las percepciones de los usuarios de cómo funciona una cuenca.
Es un modelo Australiano, el cual brinda la posibilidad de modelar escorrentía, tránsito de crecientes y calcular pérdidas en cuencas
RORB tanto urbanas como rurales a partir de eventos de tormenta. Este modelo es empleado comúnmente en gestión de inundaciones y (Laurenson, E. M., and Mein, 1995)
diseño de vertederos.
Estudio a escala regional de los componentes hidrológicos e hidráulicos del complejo
cenagoso del Bajo Magdalena entre Calamar y Bocas de Ceniza con apoyo de
modelación matemática y SIG
*ADAPTABLE
TODOS LOS MODELOS SON DETERMINISTICOS
Descritos los modelos hidrológicos considerados con mayor potencial de aplicación en esta
investigación, queda claro que la selección final del tipo de modelo aplicable en la zona,
dependerá en gran medida de la cantidad y calidad de información disponible, pensando
en la implementación rigurosa de un protocolo de modelación que incluya las fases de
calibración y validación. Es la carencia de datos una importante limitante para la aplicación
efectiva de los modelos lluvia – escorrentía y el uso adecuado de los resultados del mismo,
dado que se requiere de información coincidente de precipitación dentro de la cuenca y de
caudales en la salida de la misma, para la correcta implementación del protocolo de
modelación. Es precisamente la carencia de esta información la que obliga a adoptar
metodologías alternativas que permitan evaluar el componente hidrológico en zonas con
información limitada, siendo las ecuaciones regionalizadas una alternativa viable de
aplicación, teniendo en cuenta el uso directo en función de parámetros morfométricos
como el área o la pendiente, para describir el potencial de generación de escorrentía en
las cuencas aledañas a cuerpo hídricos.
A este respecto, en el año 1995, el antiguo HIMAT hoy en día Instituto de Hidrología
Meteorología y Estudios Ambientales (IDEAM), elaboró para todo el país un ejercicio de
regionalización con el cual se buscaba obtener una herramienta adecuada para la
estimación de caudales extremos en cauces no principales, a partir de algunas
características morfométricas, específicamente el área.
Qt
At Ecuación 1
Q2.33
Donde At es el valor establecido para cada región en función del período de retorno, Q2.33
es el caudal promedio de los caudales máximos en m3/s y Qt el caudal máximo anual con
períodos de retorno de T años.
Adicionalmente, el HIMAT obtuvo para cada región expresiones que ligan el Q2.33 (m3/s)
con el área de drenaje de cuencas dentro del rango en el cual fueron deducidas estas
ecuaciones (km2), como se presenta en la Ecuación 2.
Donde A representa el área en km2, los coeficientes C4 y C1, se encuentran tabulados para
cada una de las 21 zonas. En el caso particular de la zona de estudio, la cual
correspondería a la denominada región 14 (Occidente de la Sierra Nevada), según la
clasificación HIMAT (1995) se tendría la Ecuación 3 con los valores de C4=84.39 y C1=0.12.
Finalmente, es posible deducir los valores de Q t para Tr = 2, 2.33, 5, 10, 20, 50 y 100 años
mediante la Ecuación 4.
Qt Q2.33 At Ecuación 4
Estado del arte y propuesta metodológica 49
Cabe recordar que la aplicabilidad de la ecuación está limitada a los afluentes de las
fuentes principales del país, tales como los ríos Magdalena, Cauca, Atrato, Meta, Arauca,
etc. (IDEAM, 1995)
Por otra parte, los estudios de balance hídrico y principalmente el balance hídrico de largo
plazo (BHLP), constituyen otras herramientas importantes para el análisis hidrológico de
zonas con información limitada. La Figura 3-2, resume las principales variables de la
ecuación de balance hídrico aplicada a una parcela, se destaca: a) Lluvia total. b) Agua
almacenada en la cobertura vegetal por intercepción. c) Pérdidas por evaporación desde
la cobertura vegetal. d) Agua almacenada en las depresiones de la superficie e)
Escorrentía superficial o neta en el área de estudio. f) Infiltración. g) Agua almacenada en
la zona no saturada.
50 Estudio a escala regional de los componentes hidrológicos e hidráulicos del complejo
cenagoso del Bajo Magdalena entre Calamar y Bocas de Ceniza con apoyo de
modelación matemática y SIG
Figura 3-2: Algunos términos de la ecuación del balance hídrico de la interface suelo –
atmósfera – vegetación. Fuente: Adaptado de (Hingray, Picouet, et al., 2014)
Q K (P E) A Ecuación 5
Queda claro entonces que los estudios hidrológicos en zonas con poca instrumentación e
información (fundamentalmente carencia de información coincidente entre precipitación y
caudal a la salida de las microcuencas), como el caso del área de estudio, se ven limitados
en una rigurosa implementación de los modelos lluvia – escorrentía, así como en el uso de
otras técnicas hidrológicas alternativas, sobre todo, debido a la dificultad en las fases de
calibración y validación. Sin embargo, la aplicación de los modelos lluvia escorrentía, junto
con técnicas como el BHLP, resulta ser una muy buena alternativa para la estimación de
los caudales medios aportantes por escorrentía superficial. Esto aún sin una rigurosa
calibración, ya que, el empleo combinado de los dos tipos de métodos permite establecer
órdenes de magnitud de los caudales aportantes por las cuencas involucradas en todo el
hidrosistema y su variabilidad. De esta forma, en esta investigación se ha realizado para
cada microcuenca:
9
HidroSIG 4.0 potenciado por MapWindow (http://poseidon.unalmed.edu.co/~hidrosig)
es un Sistema de Información Geográfica que ofrece una serie de herramientas para el
procesamiento y análisis de información hidrológica y climatológica. (UNAL-MED, 2011).
Estado del arte y propuesta metodológica 53
caudal medio anual (L3/T). Este proceso se aplicó para cada uno de los años en el período
comprendido entre (1980 – 2011), eliminando de la interpolación las estaciones sin
registros.
La modelación extendida en HEC-HMS puede ser realizada de dos formas: a) A través del
modelo de contabilidad de humedad en el suelo (SMA), por sus siglas en inglés (Soil-
Moisture Accounting), el cual simula el movimiento vertical del agua utilizando
almacenamientos o tanques para la intercepción de agua por vegetación, depresión en el
suelo, infiltración y percolación (U.S. Army Corps of Engineers, 2000). Este proceso se
resume en la Figura 3-3.
54 Estudio a escala regional de los componentes hidrológicos e hidráulicos del complejo
cenagoso del Bajo Magdalena entre Calamar y Bocas de Ceniza con apoyo de
modelación matemática y SIG
Los modelos hidráulicos tienen como objetivo el cálculo de la lámina de agua en el cauce,
por lo cual reciben como insumo los caudales generados en las cuencas aferentes al cauce
principal y que son obtenidos de los análisis hidrológicos previamente elaborados. Cuando
se piensa en la delimitación de zonas inundables y zonas de riesgo, se deben tener en
cuenta las fuentes de incertidumbre que se desprenden de la modelación hidráulica para
identificar cual tiene mayor relevancia y poder así evaluar la confiabilidad y
representatividad de los resultados obtenidos. La Figura 3-5 ilustra potenciales fuentes de
incertidumbre de las modelaciones hidráulicas en canales abiertos, las cuales deben
tenerse en cuenta tanto en la etapa de selección del modelo, como en la etapa de
valoración de los resultados.
La física de los flujos de fluidos está bien entendida y es generalmente capturada por la
formulación matemática de las ecuaciones de conservación de masa, energía y
momentum. En la práctica general de la hidráulica de canales abiertos, el promedio de la
velocidad en la profundidad es una buena representación de la velocidad del flujo. El flujo
puede ser aproximado utilizando modelos en una dimensión (1D), modelos en dos
dimensiones (2D) y modelos en tres dimensiones (3D). Existen también los modelos Cuasi
2D y los modelos 1D/2D, estos últimos se caracterizan principalmente por realizar una
aproximación (1D) en el cauce principal y (2D) en la llanura inundable. En todos estos
modelos se puede considerar el análisis mediante flujo permanente cuando se establece
un caudal invariante en el tiempo o mediante flujo no permanente cuando se utiliza como
condición de entrada un hidrograma.
56 Estudio a escala regional de los componentes hidrológicos e hidráulicos del complejo
cenagoso del Bajo Magdalena entre Calamar y Bocas de Ceniza con apoyo de
modelación matemática y SIG
Los modelos bidimensionales (2D), conocidos también como modelos de aguas someras,
calculan las componentes de velocidad en el plano horizontal (x e y) y promedian la
velocidad en la profundidad (z). Estos modelos resuelven completamente las ecuaciones
de Saint-Venant utilizando algoritmos numéricos que avanzan iterativamente en la solución
sobre cuadriculas o triángulos. La solución de los flujos en dos dimensiones es
computacionalmente más demandante, sin embargo, ofrece importantes ventajas con
respecto a los modelos unidimensionales, como son la resolución más precisa de la
elevación de la superficie de agua y la determinación directa de la llanura de inundación,
a medida que avanza el cálculo. Un estudio comparativo entre modelos unidimensionales
y modelos bidimensionales encontró que los modelos (2D) tiene mayor precisión para
determinar la magnitud y dirección de la velocidad (TRB10, 2006), este último es un
elemento importante en la estimación del daño por inundación. El estudio de (TRB, 2006)
10
The Transportation Research Board, 2006
Estado del arte y propuesta metodológica 57
también encontró que las diferencias entre modelos (2D) y (1D) son pequeñas en el canal
principal, incrementan entre el canal y la llanura inundable y son mayores en las
derivaciones. Estas divergencias entre el canal y la zona inundable resultan de la
incapacidad de los modelos 1D para representar complejas características y fenómenos,
tales como corrientes trenzadas, múltiples aperturas, etc. Por consiguiente, la elección del
modelo puede afectar significativamente los niveles de inundación. FEMA (Federal
Emergency Management Agency) recomienda hacer un mayor esfuerzo hacia el uso de
los modelos de dos dimensiones (2D) cuando así lo justifiquen la geometría de la llanura
aluvial, las vías preferentes de flujo y las estructuras existentes y previstas (FEMA, 2009).
Por otra parte, la modelación en tres dimensiones implica naturalmente una mayor
demanda computacional que la de un modelo (2D) o (1D), previamente comentados,
siendo el uso de este tipo de modelación bastante restringido a zonas y casos puntuales y
zonas no muy extensas. La aplicación principal de este tipo de modelación está asociada
al estudio de estructuras hidráulicas y modelación de cauces con fondo móvil, además del
estudio del transporte de sedimentos (Di Baldassarre, 2012).
Cabe destacar que el uso de un modelo o software libre permite que tanto los modelos
hidráulicos como sus resultados, sean fácilmente utilizados y apropiados por las entidades
o en este caso Corporaciones Autónomas Regionales que tienen injerencia directa en la
zona de estudio (el Bajo Magdalena). El uso de estos modelos libres permite mejoras en
los mismos en función de una mejor calidad de información que puede ser generada a
posteriori de esta investigación, por lo tanto, han sido priorizados en esta investigación. En
la Tabla 3-3 se presenta el panorama de algunas herramientas computacionales para
modelación hidráulica, comerciales y libres, consideradas potencialmente aplicables en el
proyecto.
Tabla 3-3:
Modelo Descripción Licencia Refencias
ArcGis Un SIG permite la modelación en OD basada en algebra de mapas, no constituye un modelo físico Comercial (Di Baldassarre, 2012)
Sobek Entorno integral para la modelación hidrodinámica 1D y 2D, desarrollado por Delft Hydraulics Comercial https://www.deltares.nl/en/software/sobek/
Mike 11 Entorno integral para la modelación hidrodinámica 1D desarrollado por Danish Hydraulics Institute Comercial (Di Baldassarre, 2012)
Mike 21 Resuelve las ecuaciones de onda de superficie en 2D, desarrollado por Danish Hydraulics Institute Comercial (Di Baldassarre, 2012)
Combinación de hidrología e hidrodinámica, transporte de fondo y transporte de contaminantes, en ríos, estuarios,
RiverFLOW-2D Comercial (Llc, 2005)
costas y llanuras de inundación - desarrollado por Hydronia LLC
FLO-2D Modelo hidrodinámico para modelación 2D, es un modelo físico. Cuenta con un módulo hidrológico (lluvia - escorrentía). Comercial (Flo-2D, 2009)
Modelo hidrodinámico de aguas promediadas en la profundidad. Este modelo está adaptado al estudio de hábitat de
peces. En total el modelo se compone de 4 programas, los cuales interactúan entre si y deben ser aplicados de forma
potencialmente aplicables en la zona de estudio.
En esta investigación se han utilizado los modelos hidráulicos HEC-RAS (1D o 2D) y el
modelo bidimensional Iber (2D), principalmente por ser software libre, ampliamente
documentados, sin embargo, podría ser cualquier otro software; estos modelos en
particular cuentan con extenso soporte teórico y aplicado. Las características de ambos
modelos y los principios para el análisis de flujo en superficie libre se resumen a
continuación:
En donde Z es el nivel del fondo del canal aguas arriba (1) y abajo (2) del tramo,
denominado este término cabeza de posición, Y es la lámina de agua aguas arriba (1) y
abajo (2) del tramo, denominado este término cabeza de presión, V2/2g es la cabeza de
velocidad aguas arriba (1) y abajo del tramo (2), he son las pérdidas de energía en el tramo,
dividiéndose en pérdidas por fricción y pérdidas localizadas. Las pérdidas por fricción en
el tramo 𝐿𝑆̅𝑓 para flujo gradualmente variado en un tramo de longitud L del canal se pueden
60 Estudio a escala regional de los componentes hidrológicos e hidráulicos del complejo
cenagoso del Bajo Magdalena entre Calamar y Bocas de Ceniza con apoyo de
modelación matemática y SIG
𝑛2 𝑉 2 Ecuación 11
𝑆𝑓 = 𝑅 4/3
En el caso del flujo no permanente HEC-RAS utiliza un modelo implícito para resolver la
ecuación de continuidad para un volumen de control dado (ver Ecuación 12) y la ecuación
de cantidad de movimiento o momentum (ver Ecuación 13).
∂A ∂Q Ecuación 12
+ =0
𝜕𝑡 𝜕𝑥
∂Q ∂(𝑄𝑣) 𝜕𝑧 Ecuación 13
+ + 𝑔𝐴 ( + 𝑆𝑓 ) = 0
𝜕𝑡 𝜕𝑥 𝜕𝑥
La implementación del modelo hidráulico (1D) en la zona de estudio ofrece una importante
ventaja computacional, teniendo en cuenta que la longitud a modelar del río Magdalena es
de un poco más de 80 km (entre el K20 y K100), los caños principales de conexión río –
ciénaga suman una longitud aproximada de 90 km (caño Renegado, caño Aguas Negras
Estado del arte y propuesta metodológica 61
y caño Clarín), y las ciénagas principales son alrededor de 30, donde se destacan el
conjunto cenagoso del Atlántico, el complejo cenagoso Pajaral y la CGSM.
• Iber, por otra parte es un modelo matemático bidimensional para la simulación del
flujo en ríos y estuarios, desarrollado a partir de la colaboración del Grupo de
Ingeniería del Agua y del Medio Ambiente, GEAMA (Universidade da Coruña), del
Grupo de Ingeniería Matemática (Universidade de Santiago de Compostela), del
Instituto Flumen (Universitat Politècnica de Catalunya y Centre Internacional de
Métodos Numèrics en Enginyeria) y promovido por el Centro de Estudios
Hidrográficos del CEDEX (España). El modelo Iber consta de diferentes módulos
de cálculo acoplados entre sí. Incluye un módulo hidrodinámico, un módulo de
turbulencia y un módulo de transporte de sedimentos, de granulometría uniforme,
por carga de fondo y por carga en suspensión (Bladé, Cea, et al., 2012).
𝜕ℎ 𝜕ℎ𝑈𝑥 𝜕ℎ𝑈𝑦
𝜕𝑡
+ 𝜕𝑥
+ 𝜕𝑦
=0 Ecuación 14
𝜕 𝜕 2 ℎ2 𝜕
(ℎ𝑈𝑥 ) + (ℎ𝑈𝑥 + 𝑔 ) + (𝑈 𝑈 )
𝜕𝑡 𝜕𝑥 2 𝜕𝑦 𝑥 𝑦
Ecuación 15
𝜕𝑧𝑏 𝜏𝑠,𝑥 𝜏𝑏,𝑥 𝜕 𝜕𝑈𝑥 𝜕 𝜕𝑈𝑥
= −𝑔ℎ + − + (𝑣𝑡 ℎ )+ (𝑣𝑡 ℎ )
𝜕𝑥 𝜌 𝜌 𝜕𝑥 𝜕𝑥 𝜕𝑦 𝜕𝑦
𝜕 𝜕 ℎ2 𝜕
(ℎ𝑈𝑦 ) + (ℎ𝑈𝑦 2 + 𝑔 ) + (𝑈 𝑈 )
𝜕𝑡 𝜕𝑥 2 𝜕𝑦 𝑥 𝑦
Ecuación 16
𝜕𝑧𝑏 𝜏𝑠,𝑥 𝜏𝑏,𝑥 𝜕 𝜕𝑈𝑦 𝜕 𝜕𝑈𝑦
= −𝑔ℎ + − + (𝑣𝑡 ℎ )+ (𝑣𝑡 ℎ )
𝜕𝑦 𝜌 𝜌 𝜕𝑥 𝜕𝑥 𝜕𝑦 𝜕𝑦
62 Estudio a escala regional de los componentes hidrológicos e hidráulicos del complejo
cenagoso del Bajo Magdalena entre Calamar y Bocas de Ceniza con apoyo de
modelación matemática y SIG
𝑛2 𝑈𝑥 |𝑈|2 Ecuación 17
τ𝑏,𝑥 = 𝜌𝑔ℎ
ℎ 4/3
𝑛2 𝑈𝑦 |𝑈|2 Ecuación 18
τ𝑏,𝑦 = 𝜌𝑔ℎ
ℎ 4/3
Una gran discusión ha existido sobre el nombre a utilizar para referirse a la representación
digital del terreno, por ejemplo, se habla de modelos digitales de terreno (MDTs), Modelos
Digitales de Elevaciones (DEMs, por sus siglas en inglés Digital Elevation Models) y
Modelos Digitales de Superficie (MDS). (Bedient, Huber, et al., 2008), por ejemplo, definen
los MDT como cualquier representación digital del terreno tal como una serie de puntos,
elevaciones, líneas de contorno y composiciones triangulares o red triangular Irregular (TIN
por sus siglas en inglés), conocido también como un modelo digital de elevación (DEM);
en otras palabras, se considera que el MDT es igual al DEM. Sin embargo, otros autores
discrepan de esta definición y establecen una diferencia entre ambas terminologías de
modelos; para (Felicísimo, 1994), “Un MDT es una estructura numérica de datos que
representa la distribución espacial de una variable cuantitativa y continua”, con esta
definición indica además que puede ser representada cualquier otra variable del terreno
como la pendiente, etc., y no solo la elevación, señalando que cuando la variable a ser
representada a partir del MDT es la elevación, se constituye un modelo digital de elevación
(DEM). Por tanto, el DEM es un caso particular del MDT (Olaya, 2011).
a) b)
Según (Bedient, Huber, et al., 2008) uno de los datos geoespaciales más utilizados es el
terreno digital. Un modelo de terreno tradicionalmente se basa en contornos de igual
elevación que definen las características de la cuenca, pendiente del canal y otras
características. A través de un SIG, la representación del terreno puede ser contornos,
elevación de celdas (también conocidas como Raster DEM), y red irregular triangular (TIN).
Estado del arte y propuesta metodológica 65
Por otra parte, los sistemas de información geográfica (SIG), representan una gran
herramienta de gestión y manejo de información espacial. Existen diversas definiciones
para los SIG dependiendo del enfoque. Una definición completa puede ser la propuesta
por (Korte, 2001) y citada en (Olaya, 2011), en la cual “un SIG es un sistema que integra
tecnología informática, personas e información geográfica cuya principal función es
capturar, analizar, almacenar, editar y representar datos georreferenciados”.
La aplicación de los SIG en los recursos hídricos es variada, desde apoyar la modelación
hidrológica e hidráulica en las etapas de pre-proceso o captura de información requerida
por modelos hasta el post proceso, o representación de los resultados de modelaciones.
Además, ayudan por ejemplo en la identificación de zonas susceptibles a desastre natural
y zonas de riesgo. Algunos programas comunes de Sistemas de Información Geográfica
(SIG) se presentan en la Tabla 3-5 adaptada de (Bedient, Huber, et al., 2008) y
complementada con información de SIG más recientes y de uso libre como son
MapWindows11 y SagaGIS12. En la actualidad, la Universidad Nacional de Colombia sede
Bogotá, cuenta con licencias para la utilización del Software ArcGIS de ESRI y ERDAS
imagine. Sobre ArcGIS funcionan las extensiones GeoHMS y GeoRAS del US ARMY,
módulos ampliamente utilizados para el pre-proceso y post-proceso de las herramientas
11
http://live.osgeo.org/es/overview/mapwindow_overview.html
12
http://www.saga-gis.org/en/index.html
Tabla 3-5:
Software Descripción Módulos
ArcGIS incluye otros productos que proporcionan mapeo
Un software popular para enseñanza e investigación. Distribuido por Clark Labs, IDRISI Andes es un software integrado de SIG y
IDRISI
Clark University, Worcester, Massachusetts. procesamiento de imágenes.
Microestation Geographics provee soporte para análisis,
Microestation Geographics provee herramientas para análisis y visualización de presentación, transformación de coordenadas, creación
Microstation
información espacial en formato vector y raster. de datos, etc. Este software ofrece una interacción entre
SIG y sistemas CAD.
Interfaz gráfica de usuario amigable, fácil visualización de
muchos formatos vectoriales y raster, permite crear,
MapWindow GIS es un cliente GIS basado en Windows y Open Source (código editar y exportar datos espaciales. Realizar análisis
MapWindow GIS[1] abierto), con el que se puede visualizar, gestionar, editar y analizar datos, así espacial usando plugins (álgebra de mapas, análisis del
como componer mapas para imprimir. (software libre) terreno, Modelos hidrológicos, análisis de redes, etc.). La
herramienta HidroSIG se encuentra sobre esta plataforma
a partir de su versión 4.0
SAGA (Sistema para Análisis Automatizados Geocientíficos). Es un sistema de
arquitectura modular en el cual se pueden visualizar datos espaciales, editar, etc. Software libre. Recibe diferentes formatos, y permite
SAGA GIS[2] Creado por investigadores del Departamento de Geografía Física, Göttingen. En realizar múltiples proceso incluidos en software (s)
2007, el centro del desarrollo SAGA se movió hacia Hamburgo, donde esta ahora comerciales.
los datos requeridos en la implementación y análisis de la zona de estudio.
Para Ritchie et al., 1996, citado por (Schultz and Engman, 2001), “La Teledetección ha
sido definida como la ciencia y arte de obtener información acerca de objetos, áreas, o
fenómenos a través del análisis de datos adquiridos por un sensor que no está en contacto
directo con el objetivo de la investigación”. Un sensor remoto, según (Menéndez and
Núñez, 2009), se puede definir como un instrumento sensible, capaz de captar y evaluar
objetos ubicados a la distancia; siendo en este orden de ideas las cámaras convencionales,
los detectores infrarrojos, el radar y los sentidos de la vista, el oído y el olfato, ejemplo de
sensores remotos. Esta definición resume el potencial de los productos derivados de estas
tecnologías para ser aplicados en los recursos hídricos, siendo posible obtener información
en zonas de difícil acceso, o zonas con información escasa para eventos determinados,
etc.
Figura 3-7: Configuración gráfica planeada para la misión EOS. Fuente: Adaptado de
(Campbell and Wynne, 2011).
Un sensor remoto puede ser clasificado como pasivo o activo según la tecnología utilizada
para la captura de la información. Dentro de los sensores activos encontramos aquellos
como la tecnología radar, los cuales son capaces de emitir su propia energía, mientras que
los sensores pasivos aprovechan la energía emitida por el sol principalmente y reciben la
radiación de las coberturas.
La diferencia fundamental entre estos dos tipos de sensores radica en que los activos
pueden funcionar en cualquier instante y no dependen de las condiciones atmosféricas (no
Estado del arte y propuesta metodológica 71
se ven afectados por la presencia de nubes) o el momento del día. Los sensores activos
suelen trabajar en el rango de las microondas (frente a los sensores pasivos, que los hacen
en las regiones del visible y el infrarrojo principalmente), ya que estas son capaces de
atravesar la atmósfera en prácticamente todas las condiciones, presentando una ventaja
frente a los pasivos en este aspecto, (Olaya, 2011).
correspondientes a las fuentes rojo, verde y azul, conocido como el RGB en el rango del
visible, sin embargo, algunos sensores pueden tener bandas adicionales en el infrarrojo;
con lo cual sería posible mejorar la resolución ya que cada banda cubriría un rango de
frecuencias menor. Entre más bandas existan en el sensor mejor será la resolución
espectral, siendo las imágenes hiperespectrales con más de cien bandas las que ofrecen
la mejor resolución, permitiendo una caracterización espectral sumamente precisa (Olaya,
2011). Finalmente, la resolución radiométrica representa la capacidad del sensor para
detectar las variaciones en la radiación espectral que recibe. Tanto las imágenes en blanco
y negro como las imágenes en color trabajan con 256 niveles o 8 bits. Sin embargo, esta
clasificación depende del sensor; por ejemplo, los sensores ópticos como LANDSAT y
SPOT, trabajan con 256 niveles de codificación, mientras que las imágenes IKONOS son
de 2048 niveles o 11 bits. La mayor resolución radiométrica conlleva en teoría a una mayor
resolución. (Olaya, 2011).
La Tabla 3-6, resume algunas de las principales misiones científicas que producen la
información, sensores y las características de los productos derivados de los mismos,
destacando aspectos como resolución espacial, temporal, espectral y radiométrica. Esta
tabla se elaboró a partir de información obtenida en
http://gdsc.nlr.nl/FlexCatalog/catalog.html# (base de datos de misiones espaciales para
observación de la tierra).
Son múltiples las aplicaciones reportadas de los sensores remotos y sus productos
derivados en el estudio de los recursos hídricos, teniendo en cuenta que su potencial
aumenta por la facilidad de integración con otras tecnologías como los SIG. Es así, como
los sensores remotos pueden ser de gran utilidad para la delimitación de áreas afectadas
por inundaciones y permitir la evaluación de modelos hidráulicos (Di Baldassarre, 2012),
delimitar la evolución temporal de cuerpos de agua tales como ciénagas, lagos, lagunas,
etc., apoyar la geomorfología y dinámica fluvial a partir de estudios multitemporales e
identificar zonas de movimientos en masa, etc. (Vargas Cuervo, 2015; Zhang, Guo, et al.,
2015; Sarhadi, Soltani, et al., 2012). Sobre el recurso hídrico, (Gao, Birkett, et al., 2012),
presentan una aplicación de sensores remotos para reconstruir las curvas de
Sensor óptico Sensor radar Resolución
Longitud de onda (nm) Resolución (m) Frecuencia (GHz) temporal
Resolución espacial (m)
alcance.
Setelite Mín Máx Mín Máx Mín Máx Mín Máx (días) Operador Misión País Inicio Fin
LANDSAT 1 80 80 480 1100 28 USGS Landsat Landsat USA Jul 23, 1972 Ene 6, 1978
LANDSAT 2 80 80 480 1100 18 USGS Landsat Landsat USA Ene 22, 1975 Ene 1, 1983
LANDSAT 3 80 80 480 1100 28 USGS Landsat Landsat USA Mar 5, 1978 Ene 7, 1983
LANDSAT 4 30 120 450 12500 28 USGS Landsat Landsat USA Jul 16, 1982 Dic 1, 1993
LANDSAT 5 30 120 450 12500 16 USGS Landsat Landsat USA Mar 1, 1984 Ene 1, 2013
LANDSAT 6 15 120 450 12500 16 USGS Landsat Landsat USA
LANDSAT 7 15 120 450 12500 16 USGS Landsat Landsat USA Abril 15, 1999
LANDSAT 8 15 120 433 12500 16 USGS LDCM USA Feb 11, 2013
Ikonos - 2 1 4 450 900 3 GeoEye Ikonos - 2 USA Sep 24, 1999
SPOT 1 10 20 500 890 26 SPOT - image SPOT Francia Feb 22, 1986 Nov 1, 2003
SPOT 2 10 20 500 890 0.4 2.4 26 SPOT - image SPOT Francia Ene 22, 1990
SPOT 3 0 20 352 1104 0.4 2.4 26 SPOT - image SPOT Francia Sep 26, 1993 Nov 14, 1996
SPOT 4 0 1000 349 1750 0.4 2.4 26 SPOT - image SPOT Francia Mar 24, 1998
SPOT 5 5 1000 450 1750 0.4 2.4 26 SPOT - image SPOT Francia May 3, 2002
SPOT 6 1.5 6 455 890 1 CNES SPOT Francia Sep 9, 2012
SPOT 7 1.5 6 455 890 1 CNES SPOT Francia Jun 30, 2014
Alemania,
SRTM 40 40 450 750 30 30 5.3 9.6 NASA SRTM Italia, USA Feb 11, 2000 Feb 22, 2000
Quickbird 2 0.6 2.4 450 900 3.5 Digital globe Quickbird USA Oct 18, 2001
Estado del arte y propuesta metodológica
Aqua Mission
(EOS/PM-1) 250 20000 300 125000 5 40 6.93 183.31 2 NASA EOS USA May 4, 2002
Terra Mission
(EOS/AM-1) 15 22000 300 125000 16 NASA ESE USA Dic 18, 1999
Tabla 3-6: Principales misiones y sensores para observaciones de la tierra.
almacenamiento de varios embalses en los Estados Unidos a partir de un análisis
de campo arrojaron errores del 4 al 15%, considerados bastante aceptables, por citar un
multitemporal de las imágenes del sensor MODIS, los cuales al ser comparados con datos
73
74 Estudio a escala regional de los componentes hidrológicos e hidráulicos del complejo
cenagoso del Bajo Magdalena entre Calamar y Bocas de Ceniza con apoyo de
modelación matemática y SIG
a) b)
Figura 3-10: a) Imagen 22/11/2010, sensor Radarsat 2, (b) Imagen 22/11/2010, resultado
de la clasificación supervisada sobre una máscara o fondo para destacar la cobertura de
agua hacia la margen izquierda.
2013). Llevar a cabo la comparación de las manchas es comúnmente una actividad visual
que adolece de subjetividad. Por ello es deseable complementar este análisis visual con
algunos estadísticos que ayudan a cuantificar el grado de ajuste de la información espacial,
basado bien sea en la comparación de polígonos que definen los límites de la inundación
o en un mapa raster (MARM, 2011). De acuerdo con (Mason, Bates, et al., 2009; García-
Pintado, Mason, et al., 2015), un método tradicional para analizar el desempeño de un
modelo determinístico de inundación, cuando se encuentran disponibles datos de sensores
remotos de la extensión de la inundación, es comparar la extensión de la inundación
simulada mediante la aplicación del modelo, con la observada.
(Aronica, Bates, et al., 2002; Horritt and Bates, 2002), describen la medida de desempeño
“F” mediante la función objetivo formulada en la Ecuación 19, mientras que (Di Baldassarre,
Schumann, et al., 2009; Amador, 2013), utilizan una variación de la Ecuación 19, para el
cálculo de F, la cual se presenta en la Ecuación 20.
A Ecuación 19
F
( A B C)
( A B) Ecuación 20
F
( A B C)
(MARM, 2011) propone calcular el estadístico Kappa, el cual utiliza la valoración de mapas
lógicos [0,1], donde el valor 1 representa agua y 0 otras coberturas. Este estadístico mide
la relación entre el nivel de acuerdo observado [P(O)] sobre el mínimo nivel de acuerdo
posible [P(E)], respecto al máximo nivel de acuerdo posible (Ecuación 21 a Ecuación 23).
Estado del arte y propuesta metodológica 77
El valor máximo del estadístico Kappa es 1 y el valor mínimo -1. (MARM, 2011) recomienda
además el uso del valor del parámetro K límite (KLIM) cuando se utiliza este estadístico, el
cual está relacionado con el estadístico “F” mediante la Ecuación 24, considerando un valor
satisfactorio para valores de KLIM ≥ 0.90, lo cual correspondería a un valor “F” de 0.82.
P (O ) P ( E ) Ecuación 21
K
1 P( E )
nss n AB Ecuación 22
P(O)
nT
2F Ecuación 24
K LIM
1 F
En esta investigación, otra aplicación relevante de los sensores remotos incluyó el análisis
de la geomorfología y la dinámica fluvial específicamente aplicada al río Magdalena, con
la cual se realizó la delimitación de zonas susceptibles a inundación, zonas de flujo
preferente y llanuras de inundación. Con la delimitación histórica del cauce fue posible,
además de establecer las zonas de flujo preferente, obtener las zonas de orillas con mayor
actividad y la envolvente de orillas máximas identificadas, también la dinámica de los
cuerpos lenticos. Los resultados del cauce histórico posteriormente se han contrastado con
los obtenidos de la modelación hidráulica, analizando las discrepancias y similitudes
encontradas, en aras de establecer zonas inundables a partir de un criterio conjunto.
78 Estudio a escala regional de los componentes hidrológicos e hidráulicos del complejo
cenagoso del Bajo Magdalena entre Calamar y Bocas de Ceniza con apoyo de
modelación matemática y SIG
El proceso se inicia con la identificación de las entidades que para la zona de estudio
suministran información referente a: precipitación, caudales, topografía, batimetrías,
coberturas, imágenes satelitales, etc.; luego se realiza la integración de la información
cartográfica primaria y secundaria, con la cual fue posible generar el MDT integrado de la
zona de estudio. La integración de la información se debe realizar teniendo como criterio
una misma resolución en todos los productos, en caso de no ser posible, se integra la
información teniendo como referencia la información con menor resolución, dado que, esto
garantiza una mayor validez de la información generada. Posteriormente, es abordado el
análisis hidrológico de la zona, donde se ha implementado primero un balance hídrico de
largo plazo (BHLP) con el cual se establecen en todas las cuencas aferentes al
hidrosistema los caudales con agregación anual, seguidamente la modelación matemática
utilizando HEC-HMS (modelación extendida), determinando así los caudales de
escorrentía superficial con agregación diaria, para el año normal, año Niño y año Niña,
previamente establecidos. Por otra parte, se determina el cauce máximo a partir del
análisis de la dinámica fluvial y la geomorfología, utilizando información de sensores
Estado del arte y propuesta metodológica 79
mancha de inundación y así tener un criterio objetivo para la obtención de los mejores
valores de los parámetros sensibles, se propone evaluar el factor de desempeño “F”
(discutido en el ítem 3.5 de este Capítulo) que relaciona las coincidencias o no en la
representación de la mancha de agua entre las simulaciones del modelo hidráulico (1D) y
la clasificación supervisada para la fecha de calibración, destacando que la evaluación se
realiza a la resolución espacial de la imagen satelital y para la llanura inundable. Esto último
con el objetivo de ser más críticos en la evaluación de los resultados del modelo hidráulico,
ya que al incluir en la evaluación la coincidencia de la mancha en el cauce se considera
que no se estaría evaluando adecuadamente la respuesta hidráulica del modelo, pues
aumentaría significativamente el valor del factor de desempeño “F”, pudiéndose
enmascarar deficiencias en la representación de la mancha de inundación. Con esta
estrategia se definen los mejores valores del coeficiente de rugosidad “n” de Manning en
el cauce principal y “Cd” en la estructura de derivación; adicionalmente, se obtienen los
niveles en el cauce y en las ciénagas con los cuales se alimenta ahora el modelo hidráulico
(2D), en caso de contar con la información requerida para hacer una modelación de este
tipo.
entendimiento hidrodinámico en las ciénagas y/o la llanura inundable, ya que este permite
establecer la distribución en el plano horizontal (x e y) de la velocidad. La Figura 3-14
resume la metodología propuesta.
Delimitación del
Acondicionar MDT
cauce en la imagen
con HEC -GeoDozer
más antigua.
Clasificación
Aplicación del supervisada de la
modelo Floodem y imagen satelital
reportar inundación destacando
cobertura agua
No
No Nuevas zonas
Reacondicionar ¿De acuerdo con ¿De acuerdo con
el DTM los resultado ? de
los resultado ?
entrenamiento
Si Si
Almacenar
resultados
Reportar zonas
susceptible por
inundación
Si
Existen más
imágenes ?
No
Elaborar
multitemporal de
resultados
Reportar envolvente
máxima de terrenos
pertenecientes al río.
Cabe destacar que los aportes de caudales esperados para las áreas aferentes en la zona
de estudio, se estiman (a priori) considerablemente mucho menores a los que transporta
el río Magdalena desde Calamar, una vez se deriva el caudal para el Canal del Dique, por
lo que la hidrología de las cuencas aferentes no impactaría dramáticamente los resultados
de modelaciones hidráulicas para zonas de amenaza por inundación en el área de estudio
aferente al río Magdalena en su parte baja. Sin embargo, estos caudales aportantes si
permiten realizar una modelación más completa de los escenarios de intercambio hídrico
río-ciénaga o viceversa, bajo diferentes condiciones de variabilidad climática (caso ENOS).
Este Capítulo se inicia con la descripción de los resultados del estudio hidrológico para la
zona de estudio, delimitación de cuencas y estimación de la escorrentía superficial con
agregación anual, a partir de la aplicación del BHLP y la estimación de la escorrentía con
agregación diaria implementando la herramienta computacional HEC-HMS v4.1
(modelación continua o extendida). Seguidamente, se presenta la modelación hidráulica y
en ella se implementa la calibración del modelo hidráulico a partir de los resultados de la
mancha de inundación extraída de la interpretación de una imagen satelital. Esta
calibración se efectúa utilizando como referencia una zona reducida dentro del área de
estudio (20 km), la cual presenta desborde lateral e interacción río Magdalena – ciénagas.
Posteriormente se aplica la modelación a la zona amplia de estudio, utilizando el
coeficiente “n” de Manning que permitió la mejor aproximación a la mancha de inundación
de referencia. Continúa este Capítulo con el detalle de los resultados del análisis de la
dinámica fluvial del río Magdalena entre Calamar y Bocas de Ceniza (Barranquilla), y la
dinámica de los espejos de agua de las ciénagas, información que posteriormente es
utilizada en la estimación de zonas susceptibles a inundación.
92 Estudio a escala regional de los componentes hidrológicos e hidráulicos del complejo
cenagoso del Bajo Magdalena ente Calamar y Bocas de Ceniza con apoyo de
modelación matemática y SIG
Tc (min) -
Tc (h) Tiempo de Tiempo de
2
ID COUNT Area (km ) Elev. Min. (mnsnm) Elev. Max. (msnsm) Elevación promedio (msnsm) Longitud cauce ppal (Km) Pendiente (m/m) CN (SCS) Uso de tierra (SCS) concentración concentración Tlag (min) Estación PM Código
2 15737796 393.4 1 4252 1360.79 35.2 0.121 55 Bosques cubierta buena 13.37 802.33 481.40 El Enano 29060160
4 19068048 476.7 1 4453 1175.41 47.89 0.093 55 Bosques cubierta buena 19.50 1169.86 701.91 Prado Sevilla 29065030
estudio (Bajo Magdalena).
8 21291840 532.3 24 5580 1800.01 65.61 0.085 55 Bosques cubierta buena 26.28 1576.78 946.07 Carmen_El 29060140
12 4501548 112.5 4 413 47.33 9.53 0.043 55 Bosques cubierta buena 7.89 473.21 283.92 Padelma 29065020
13 30085632 752.1 32 5673 2858.15 79.98 0.071 55 Bosques cubierta buena 33.74 2024.37 1214.62 San Sebastian 29060090
16 7756192 193.9 1 54 5.48 18.03 0.003 58 Vega de ríos 46.53 2791.89 1675.13 San Rafael 29060540
17 3576161 89.39 1 60 6.26 3.3 0.018 58 Vega de ríos 4.85 290.99 174.60 San Rafael 29060540
18 93382357 2336 1 542 37.74 57.4 0.009 58 Vega de ríos 65.63 3937.57 2362.54 Doña María 29060560
19 42918490 1073 32 3917 1262.22 82.19 0.047 55 Bosques cubierta buena 42.12 2527.32 1516.39 San Sebastian 29060090
21 7405305 185.2 1 74 21.02 14.28 0.005 58 Vega de ríos 29.28 1756.82 1054.09 Tiogollo 29050010
26 5904584 147.7 1 124 33.80 9.12 0.013 58 Vega de ríos 12.59 755.58 453.35 Tiogollo 29050010
32 1667529 41.67 1 53 5.55 2.34 0.022 58 Vega de ríos 3.30 198.25 118.95 Apto_Ernesto_C 29045020
34 27220614 680.8 1 113 4.17 19.8 0.006 58 Vega de ríos 36.15 2169.19 1301.51 San Rafael 29060540
51 1397308 34.91 0 65 6.06 1.7 0.038 58 Vega de ríos 1.95 117.05 70.23 Sevillano 29060310
52 2543586 63.59 1 131 39.20 13.08 0.010 79 Pastizales pobres 11.09 665.20 399.12 Polonuevo 29040080
53 2129709 53.24 1 152 61.24 16.12 0.009 79 Pastizales pobres 13.50 809.87 485.92 Apto_Ernesto_C 29045020
54 4612982 115.3 1 139 50.86 17.36 0.008 79 Pastizales pobres 15.55 932.84 559.70 Polonuevo 29040080
Hidrología, hidráulica y dinámica fluvial en la zona de estudio.
55 2681298 67.03 1 61 19.44 9.26 0.006 79 Pastizales pobres 10.42 624.96 374.97 Montebello 29040020
56 3626846 90.67 1 201 61.31 18.8 0.011 79 Pastizales pobres 14.32 859.45 515.67 Lena 29040200
57 4146054 103.7 1 137 21.23 13.89 0.010 79 Pastizales pobres 11.72 703.19 421.92 Candelaria 29040260
58 3979857 99.5 1 124 33.22 19.54 0.006 79 Pastizales pobres 19.21 1152.34 691.40 Montebello 29040020
59 4287416 108.1 1 144 30.58 7.74 0.018 79 Pastizales pobres 5.34 320.65 192.39 Apto_Ernesto_C 29045020
60 8286894 207.2 1 311 83.81 36.67 0.008 79 Pastizales pobres 27.42 1645.33 987.20 Montebello 29040020
61 2184575 55.35 1 130 46.71 12.91 0.010 79 Pastizales pobres 10.94 656.51 393.90 Apto_Ernesto_C 29045020
62 1239610 31.36 1 27 5.06 6 0.004 58 Vega de ríos 15.89 953.57 572.14 Campo de la Cruz 29040250
63 3402812 85.07 1 205 50.26 23.78 0.009 79 Pastizales pobres 19.25 1155.02 693.01 Montebello 29040020
64 2894235 72.36 1 132 32.86 17.26 0.008 79 Pastizales pobres 15.84 950.27 570.16 Montebello 29040020
Principales parámetros morfométricos en las microcuencas del área de
93
94 Estudio a escala regional de los componentes hidrológicos e hidráulicos del complejo
cenagoso del Bajo Magdalena ente Calamar y Bocas de Ceniza con apoyo de
modelación matemática y SIG
Los resultados de caudal producto de la aplicación del BHLP se presentan en la Figura 4-5.
Nótese que los mayores valores de caudal medio anual corresponden a las microcuencas
del costado oriental, con un valor total acumulado cercano a los 22 m3/s (año 2000). Los
caudales medios de escorrentía en el costado occidental son mucho menores y suman
cerca de 3.0 m3/s (año 2000), aproximadamente un 14% del caudal medio total del costado
oriental. Cabe recordar que la resolución espacial del álgebra de mapas corresponde a un
tamaño de celda de 30 x 30m (ASTER G DEM).
Figura 4-5: Resultado de la implementación del Balance Hídrico de Largo Plazo (BHLP)
para el año 2000. Caudal en m3/s.
En total, el caudal medio por escorrentía superficial esperado en el tramo del área de
estudio es de aproximadamente 25 m3/s (año 2000). Al comparar estos valores con los
caudales de entrada al hidrosistema a través del río Magdalena después de la derivación
del Canal del Dique, se aprecia que el caudal promedio para el río Magdalena en el año
Hidrología, hidráulica y dinámica fluvial en la zona de estudio. 97
2000 fue de 9427 m3/s, siendo la escorrentía de la zona tan solo un 0.30% del caudal
medio del río Magdalena para este año, aproximadamente.
50 200
40 150
1993
2000
2007
1980
1981
1982
1983
1984
1985
1987
1988
1989
1990
1991
1992
1994
1995
1996
1997
1998
1999
2001
2002
2003
2004
2005
2006
2008
2009
2010
Año
0 1 2 3
4 5 6 7
8 9 10 11
12 13 14 15
16 17 18 19
20 21 22 23
Figura 4-6: Caudal medio anual en cada microcuenca con BHLP y caudal acumulado
período 1980 – 2010.
160 12000.00
cuencas bajo Magdalena m3/s
Caudal promedio acumulado
140
10000.00
Caudal promedio
8000.00
100
80 6000.00
60
4000.00
40
2000.00
20
0 0.00
1989
1990
2000
2010
1980
1981
1982
1983
1984
1985
1986
1987
1988
1991
1992
1993
1994
1995
1996
1997
1998
1999
2001
2002
2003
2004
2005
2006
2007
2008
2009
AÑO
Figura 4-7: Caudal promedio anual en el Bajo Magdalena estimado con BHLP y caudal
promedio anual río Magdalena obtenido de información hidrométrica período 1980 – 2010.
Una vez conocidas las magnitudes promedio de los caudales en la zona de estudio y
efectuada su comparación con los caudales en el río Magdalena, se ha implementado la
modelación matemática en hidrología con el propósito de estimar el posible aporte por
escorrentía superficial hacia las ciénagas, en año Normal, año Niño y año Niña, para una
resolución diaria. Lo anterior, considerando la influencia de los fenómenos macroclimáticos
en la zona, descritos en el Capítulo 2 de este documento.
Los resultados para cada una de las microcuencas se presentan en la Tabla 4-1. En esta
misma tabla se presenta también la estación de precipitación asociada a cada
microcuenca, la cual, se determinó a partir de polígonos de Thiessen, utilizando las
100 Estudio a escala regional de los componentes hidrológicos e hidráulicos del complejo
cenagoso del Bajo Magdalena ente Calamar y Bocas de Ceniza con apoyo de
modelación matemática y SIG
estaciones con mejor cantidad de información para el período comprendido entre 2008 y
2011, período durante el cual y de acuerdo con reportes del ONI (ver Tabla 4-2), el país
pasó de una condición normal en el trimestre AMJ (Abril, Mayo, Junio) del 2009, a una fase
Niño que se extendió hasta el trimestre MAM (Marzo, Abril, Mayo), del 2010, luego siguió
una fase Neutra hasta el trimestre MJJ, y una fase Niña desde el trimestre JJA (Junio,
Julio, Agosto), del 2010, que se extendió hasta (Marzo, Abril, Mayo), del 2011. Finalmente,
en el trimestre ASO (agosto, septiembre, octubre), de 2011, nuevamente se presentó una
fase Niña, que se extendió hasta comienzos del año 2012. Estas últimas fechas
corresponden a la fase la Niña que ocasionó una gran catástrofe en la zona, incluida la
rotura del Canal del Dique en 2011, popularmente conocidas como la “ola invernal” 2010-
2011. En la Tabla 4-2 se aprecia el promedio de temperatura en el Océano Pacífico.
Tabla 4-2: Reporte trimestral del comportamiento de fenómeno ENOS en el período 1981
– 2011, identificando fases Niño (color rojo), Niña (color azul) y Normal (color gris). Fuente:
NOAA.
Hidrología, hidráulica y dinámica fluvial en la zona de estudio. 101
La Figura 4-11 resume los parámetros de calibración, así como los resultados de los
parámetros calibrados (optimizados) que han sido utilizados en todas las microcuencas del
hidrosistema Bajo Magdalena. Cabe destacar que el parámetro más sensible resulta ser
el almacenamiento máximo de la cobertura (Simple Canopy – Max Storage). El lag Time
no es un parámetro de la cuenca.
Una vez han sido establecidos los mejores parámetros para la representación de
abstracciones y cálculo de la ETR en la cuenca piloto o calibración (ID 13), y teniendo en
cuenta que cuatro de las microcuencas provenientes de la SNSM (ID 4, 8, 13 y 19), ver
Figura 4-1, cuentan con estaciones de caudales en la salida de las mismas y por tanto
tienen datos de caudal, se procede a implementar la modelación matemática a las 23
cuencas restantes del área de estudio, en aras de conocer su respuesta hidrológica en el
período 2008 – 2011, utilizando los valores optimizados de las variables calibradas y
presentadas en la Figura 4-11. En la Figura 4-12, se presenta el modelo hidrológico general
implementado para la zona de estudio, mientras que los resultados de la modelación se
sintetizan en la Figura 4-13. Esta última figura detalla los resultados de las microcuencas
para el período 2008 – 2011, destacando que estos hidrogramas de respuesta serán
incluidos en la modelación hidrodinámica de la zona de estudio como aportes de
escorrentía superficial hacia las ciénagas, las cuales se conectan directa o indirectamente
(a través de caños) con el río Magdalena, entre Calamar y Bocas de Ceniza.
106 Estudio a escala regional de los componentes hidrológicos e hidráulicos del complejo
cenagoso del Bajo Magdalena ente Calamar y Bocas de Ceniza con apoyo de
modelación matemática y SIG
Figura 4-14: Zona para ejercicio de calibración del parámetro sensible “n” de Manning y
discretización en el modelo hidráulico 1D.
13
El K0 corresponde a la desembocadura del río Magdalena en el mar Caribe (Bocas de Ceniza).
Hidrología, hidráulica y dinámica fluvial en la zona de estudio. 109
Amador en el año 2013, (Amador, 2013), reporta para este período la Ecuación 27,
obtenida para las condiciones de salida aguas abajo de la zona de estudio, en función de
los caudales registrados aguas arriba (diferencia entre los caudales en la estación Calamar
y los caudales en la estación Incora K7). Esta ecuación se obtuvo a partir de datos de la
estación Puerto PIMSA del LEH-Flores, la cual se encuentra en el K38, inmediatamente
aguas abajo de la zona escogida para calibración y análisis de sensibilidad. Analizando la
110 Estudio a escala regional de los componentes hidrológicos e hidráulicos del complejo
cenagoso del Bajo Magdalena ente Calamar y Bocas de Ceniza con apoyo de
modelación matemática y SIG
3.5
7
NIVEL DE AGUA (m.s.n.m.)
50%
3
6
2.5
5
90%
4 2
3 MIN. 1.5
2010
2 1 2011
2012
1
0.5
ENE FEB MAR ABR JUN JUL AGO OCT DIC
MAY SEP NOV
0
0
2010 2011 2012 ENE FEB MAR ABR MAY JUN JUL AGO SEP OCT NOV DIC
Figura 4-15: Gráfica de niveles de agua en las estaciones Calamar y Puerto PIMSA (2010-
2012). Fuente: LEH – Flores, plano 17. Ref: E-004-013-12, mayo de 2012.
No obstante, cabe anotar que más aguas abajo esta condición cambia y por tanto esta
consideración no aplicaría estrictamente para la totalidad de la zona de estudio hasta el
Puente Pumarejo (objeto de la modelación total – zona de estudio) y por esta razón la
condición de frontera más adelante tiene un tratamiento diferente. A partir de esta ecuación
se establece la condición de frontera aguas abajo, para todas las modelaciones en el
escenario de calibración de parámetros. En la Figura 4-16 se presentan las condiciones de
entrada y salida utilizadas para la modelación unidimensional con flujo no permanente.
Hidrología, hidráulica y dinámica fluvial en la zona de estudio. 111
En donde h es el nivel del agua en el sector más aguas abajo de la zona de estudio (m) y
Q es el caudal de entrada en el río Magdalena (diferencia Calamar- Incora K7) en m3/s.
20000 5
Nivel (m.s.n.m)
4
Caudales m3/s
15000
3
10000
2
5000 1
0 0
0 100 200 300 400 500 0 100 200 300 400 500
Dias Dias
Figura 4-16: Hidrograma de caudal en la entrada y relación nivel contra tiempo aguas
abajo, para la modelación 1D. La línea con cuadros representa la fecha de la imagen del
22 de octubre de 2010 y la línea con triángulos la imagen del 28 de diciembre de 2011.
En síntesis, los mejores parámetros para la zona de estudio se presentan en la Tabla 4-3
y Table 4-4, mientras que el contraste de los mejores resultados sobre la imagen
clasificada para el evento ocurrido el 28 de diciembre de 2011, se aprecian en la Figura
4-17.
112 Estudio a escala regional de los componentes hidrológicos e hidráulicos del complejo
cenagoso del Bajo Magdalena ente Calamar y Bocas de Ceniza con apoyo de
modelación matemática y SIG
Los resultados presentados en la Tabla 4-3, sugieren, de acuerdo con las diversas
formulaciones del factor de desempeño utilizado para evaluar los reportes de la modelación
hidráulica 1D contra la imagen clasificada para el día 28 de diciembre de 2011, que
independientemente de la formulación, la mejor combinación de parámetros se mantiene
dentro de un mismo rango. Por ejemplo, utilizando el factor de desempeño propuesto por
(Aronica, Bates, et al., 2002), la mejor combinación del coeficiente de rugosidad “n” de
Manning en el cauce “nc”, “nll” Manning en la llanura y coeficiente de descarga para el
vertedero lateral de pared gruesa “Cd”, resulta estar en el rango (0.020 - 0.024), para el
coeficiente de Manning y entre 1.30 – 1.70, para el coeficiente “Cd”, resultando un valor de
56.08%, mientras que, cuando se usa la expresión propuesta por (Di Baldassarre,
Schumann, et al., 2009), el mayor valor obtenido es de 33.89%, con el mismo rango de
parámetros. Para el caso del factor Kappa, al relacionarse con el factor “F” propuesto por
(Aronica, Bates, et al., 2002), la mejor combinación se encuentra también en este mismo
rango. Las diferencias entre los valores absolutos de los distintos valores de factores de
desempeño “F” calculados, se deben al castigo de la sobre estimación en el modelo
hidráulico realizado en la expresión propuesta por (Di Baldassarre, Schumann, et al.,
Hidrología, hidráulica y dinámica fluvial en la zona de estudio. 113
Por su parte, la Figura 4-17 muestra la diferencia entre los mejores resultados de los
factores de bondad de ajuste “F” en los modelos hidráulicos 1D y 2D, figuras a y b
respectivamente, y el peor resultado de todas las modelaciones figura c. El mejor resultado
corresponde a la combinación nc=0.020, nll1=0.028, y nll2 = 0.032, m-1/3s, con la cual se
obtuvo una medida de desempeño “F” del 56.20% con la expresión propuesta por Aronica,
del 34.19% con la expresión propuesta por Di Baldassarre y del 71.95% según el
parámetro Kappa (Figura 4-17 (b)). A partir de esta combinación de rugosidades se
realizaron las modelaciones en una y dos dimensiones para toda la zona de estudio
considerada inundable, teniendo en cuenta que se presentan características de coberturas
vegetales muy homogéneas en casi toda la zona baja del hidrosistema regional.
Hidrología, hidráulica y dinámica fluvial en la zona de estudio. 115
a) b) c)
caños de derivación, este tamaño de celda permite tener mínimo 6 puntos para definir el
cauce en el modelo hidráulico 1D, permitiendo un adecuado balance entre la información
disponible y el tiempo computacional requerido.
Figura 4-18: Sección transversal típica para los caños Renegado, Aguas Negras y Clarín.
Fuente: Adaptado de (CORPAMAG-DNP-PROCIÉNAGA, 1993).
25
20
15
Elevación (m)
10
5
0
-20000 0 20000 40000 60000 80000 100000 120000 140000 160000
-5
-10
Volumen (m) - 1000 m3
sí fue posible contar con la ubicación y dimensión de los mismos (Departamento del
Atlántico).
R
A IN
CL salida_CGSM
PALMAR_BARR2
12 CLARIN
3257
Malambo 3254
3282_CGSM
GD LEN A
A2
3255 AC
AR
A
1
AN2 CATUCURINCA
ElConvento 9AN3
MA
PAJARAL A
AR
1
PALMAR_BARR
Sabanagrande 3256
AS
Sabanagrande2 ARACA1
R
EG
SantoTomas
AN1 3280
AN 3258
2 MA GD ALEN A FUNDACION_BAJA
LuisaManati 3283
Paraiso
Paraiso2 3261
3284
REN_33262
6REN_4
3259
3265 3264
PONEDERA_PALM REN_2 3260
Buenavista
1 REN_1 3263
ComplejoUvero
Pinon
A
L EN
DA
Tiogollo
G
MA
N
SATiogollo2
SA_1
CALAMAR_PONED
SanAntonio
a)
b)
En el caso del río Magdalena la información del Puente Pumarejo K20 hacia aguas abajo,
se representó utilizando como base las últimas secciones de la cual se tuvo información
batimétrica, dado que la información de batimetría levantada por la Universidad Nacional
sede Bogotá, para CORMAGDALENA, y amablemente suministrada para esta
investigación, inicia desde este puente hacia aguas arriba. La extensión de esta
información hacia aguas abajo se realiza por la necesidad de incluir las derivaciones o
Hidrología, hidráulica y dinámica fluvial en la zona de estudio. 119
entradas desde el caño Clarín que se encuentra en el K13 del río, aproximadamente. En
consecuencia, las secciones transversales generadas en la modelación ascienden a más
de 2,700, estas secciones se generaron cada 100 metros, y en algunos casos para
procurar la estabilidad numérica de la modelación (cumplir condición de Courant) fueron
cada 10m, principalmente en el caso de los caños. La razón estas secciones es buscar
representar la mayor cantidad de información posible en los cauces en balance con la
demanda computacional y se realizó directamente con ayuda de HEC - GeoRAS.
La condición de frontera aguas arriba se definió como la diferencia entre los registros de
caudal de la estación Calamar (29037020) ubicada sobre el río Magdalena y la estación
Incora K7 (29037360) ubicada sobre el Canal del Dique, como ha sido comentado
anteriormente. La condición de frontera aguas abajo resultó ser más compleja de definir,
por la falta de información; la importancia de esta radica en la necesidad de poder
representar en el hidrosistema el efecto de la marea proveniente del Mar Caribe. Para ello
se contó con los datos de la estación limnimétrica Dársena Barranquilla (29047080),
ubicada aguas abajo del Puente Pumarejo, aproximadamente en el K18 del río Magdalena,
la cual tiene información de niveles medios del río a partir del año 2008. Esta estación
resultó ser la mejor opción para la condición de frontera aguas abajo ya que representa el
cierre del sistema modelado y por tanto delimita el período de análisis de la modelación
hidráulica entre 2008 – 2011, pues en esta fecha se tienen caudales medios de entrada y
niveles medios de salida en el Bajo Magdalena. Al representar los caudales de entrada vs
los niveles en la estación Dársena Barranquilla, Figura 4-22, es posible apreciar un
comportamiento casi lineal entre el hidrograma de entrada y los niveles en la salida,
mostrando que cuando los caudales de entrada son altos los niveles de salida también lo
son, lo cual sugiere que aún en el K18, condición de frontera aguas abajo, no se evidencia
efecto de marea o es poco perceptible, lo cual es complementado con la Figura 4-23 donde
el factor de correlación es 0.94, muy cercano a 1.0.
De otra parte, a través del análisis de frecuencia realizado a las estaciones Calamar e
Incora K7, ubicadas sobre río Magdalena y el Canal del Dique respectivamente, se
establecen los caudales esperados para los diferentes períodos de retorno (Tr). Estos
caudales resultan de gran utilidad al momento de ser asociados con los presentados en el
río Magdalena, permitiendo establecer la cobertura o mancha de inundación a partir de la
120 Estudio a escala regional de los componentes hidrológicos e hidráulicos del complejo
cenagoso del Bajo Magdalena ente Calamar y Bocas de Ceniza con apoyo de
modelación matemática y SIG
Figura 4-22: Caudal de entrada Vs. niveles en la salida del área general de estudio.
CALAMAR vs DARSENA B.
3
Nivel estación Dársena Barranquilla
2.5
2
(m.s.n.m)
1.5 y = 0.0043x0.67
R² = 0.9378
1
0.5
0
0 2000 4000 6000 8000 10000 12000 14000 16000 18000
Caudal de entrada m3/s
Figura 4-23: Relación caudal de entrada (Calamar-Incora K7) Vs. niveles en la salida del
área general de estudio.
Los caudales entrantes a la zona de estudio para los períodos de retorno de 5, 10, 50, y
100 años, son resumidos en la Tabla 4-5, para el período de retorno más comúnmente
utilizado en el país (Tr 100 años) se puede esperar en la estación Calamar un caudal
aproximado de 18,100 m3/s (el máximo del intervalo de confianza), mientras que en la
estación Incora K7 este caudal sería aproximadamente de 1,400 m3/s; por consiguiente, el
caudal entrante en la zona de estudio sería del orden de 16,700 m3/s, para este período.
Acorde, en orden de magnitud, con el caudal entrante a la zona de estudio reportado por
(Amador, 2013), quien reporta un valor de 16,069 m3/s y con los caudales reportados por
(Consorcio Dique, 2005), donde se reporta un caudal de 18,429 m3/s en Calamar y 1,322
m3/s en Incora K7, es decir 17,107 m3/s hacia el área de estudio. De otra parte, producto
124 Estudio a escala regional de los componentes hidrológicos e hidráulicos del complejo
cenagoso del Bajo Magdalena ente Calamar y Bocas de Ceniza con apoyo de
modelación matemática y SIG
Así, la exigencia computacional para cada uno de los modelos hidráulicos implementado
en la zona y sus principales características se resume en la Tabla 4-7, el computador
utilizado fue Core i 7, RAM 12. Mientras que los productos: intercambios hídricos y
manchas de inundación, se presentan en el anexo 2 y se discuten los resultados en el
Capítulo 5. Los órdenes de Magnitud de caudales en el río Magdalena superan
notablemente la sumatoria de caudales en todas las cuencas aportantes de la margen
izquierda, estos últimos en su totalidad, constituyendo un máximo del 28% del caudal del
río para todos los períodos de retorno.
Tabla 4-5: Caudales entrantes a la zona de estudio para diferentes períodos de retorno en
el río Magdalena.
Tr 2.33 Tr 5 Tr 10 Tr 50 Tr 100 Tr 500
Estación
Tabla 4-6: Caudales en las cuencas del Bajo Magdalena para diferentes períodos de
retorno.
Modelación / Exigencia
herramienta Características principales HidroDEM Periodo de modelación computacional
• 30 x 30 (m) Cauce
principal
2D • Longitud de cauces = 80 Km
• 90 x 90 (m) área Evento (1día) 24 horas
IBER v 2.3 • Área aferente = Buffer 2Km
aferente
126 Estudio a escala regional de los componentes hidrológicos e hidráulicos del complejo
cenagoso del Bajo Magdalena ente Calamar y Bocas de Ceniza con apoyo de
modelación matemática y SIG
Figura 4-28: Dinámica fluvial del río Magdalena entre Calamar y Bocas de Ceniza, zona
inundable por geomorfología (elaboración propia).
resultados de modelación con una imagen reportada en la zona para el mes de diciembre
2011, arrojando resultados satisfactorios, ya que, se logra evidenciar la zona de desborde
en la isla Rondón, K18 aproximadamente; sin embargo, la extensión de la mancha sobre
la llanura inundable no se representa de forma adecuada por causas asociadas con la
limitada calidad de la información topográfica disponible. En este Capítulo también se
detallan los resultados del análisis de frecuencias elaborado para las estaciones
hidrológicas Calamar e Incora K7, y a partir de las cuales fue posible establecer los
caudales máximos entrantes en la zona de estudio para diferentes períodos de retorno y
se implementó también una ecuación de regionalización elaborada por el HIMAT con los
parámetros de la zona, en aras de conocer los valores máximos de estos caudales
esperados en las microcuencas asociadas a las ciénagas; se encontró de la comparación,
que los caudales instantáneos esperados en el río Magdalena, como era de esperarse,
son varios órdenes de magnitud más grandes a los generados en esta microcuencas por
lo que no fueron tenidos en cuenta para establecer el alcance de la franja de inundación
asociada a cada período de retorno a través del criterio hidráulico; de otra parte, en este
Capítulo, se presenta también el resultado de la franja de inundación con criterio
geomorfológico, el cual fue elaborado a partir del análisis multitemporal de imágenes sobre
el río Magdalena, permitiendo así establecer la vega de divagación del río entre Calamar
y Bocas de Ceniza, esto también se detalló en el Capítulo 2. Además, se implementó la
herramienta computacional FLOODEM (UNALMED) en la zona de estudio; sin embargo,
los resultados no fueron satisfactorios, nuevamente debido a la incertidumbre existente en
la información topográfica, que limitó la generación de una adecuada red de drenaje.
Insumo indispensable para la adecuada implementación de este modelo geomorfológico,
por esta razón, estos resultados no se presentan aquí.
La dinámica general de los flujos en el río Magdalena se resume en la Figura 5-1, donde
se aprecia que la respuesta de los intercambios es directamente proporcional al pulso o
caudal recibido a la entrada del hidrosistema, principalmente por los caudales del río
Magdalena en Calamar, con algunas ciénagas en donde ingresan flujos (valores positivos
132 Estudio a escala regional de los componentes hidrológicos e hidráulicos del complejo
cenagoso del Bajo Magdalena ente Calamar y Bocas de Ceniza con apoyo de
modelación matemática y SIG
en la Figura 5-1), y otras por el contrario en donde ocurren descargas al río Magdalena
(valores negativos). Por ejemplo, en la ciénaga Luisa-Manatí (K50), costado occidental, las
mayores ganancias de caudal se presentan entre los meses de septiembre y diciembre de
2010, mientras que la ciénaga San Antonio (K100) y Tiogollo (K90) costado oriental y
Malambo (K27) costado occidental, están descargando la mayor cantidad de caudal al río
en este mismo período. Este comportamiento cíclico se aprecia también en los períodos
mayo – julio de 2011 y diciembre 2011 - enero de 2012; períodos asociados a los caudales
máximos en el río Magdalena, coincidentes con el período húmedo La Niña. El
comportamiento contrastante de las ciénagas, se debe principalmente a la interconexión
entre ellas y en parte a los caudales aportantes de las microcuencas asociadas a cada una
de ellas; las ciénagas ubicadas más aguas arriba aumentan primero su nivel y reciben la
mayor cantidad de caudal desde el río Magdalena, luego derivan caudal a las ciénagas
circundantes a través de los diques naturales que las dividen y caños de interconexión y
estas últimas, luego descargan al río Magdalena más aguas abajo, donde se presenta un
gradiente hidráulico favorable hacia el río. De otra parte, los menores intercambios son
apreciables en el período diciembre 2009 – febrero 2010, período seco El Niño, donde se
evidencia una disminución del caudal en el río Magdalena y la desconexión hidráulica del
conjunto de ciénagas y caños ubicados alrededor de él, nótese en este período la
coincidencia entre el hidrograma de entrada en Calamar (RS 101071), línea punteada color
morado y el hidrograma de salida (RS 100), línea negra en la Figura 5-1.
El perfil de flujo de la Figura 5-2 (diciembre de 2010) también se relaciona con la Figura
5-4 en la fase La Niña del ENOS, donde se aprecia como los caudales derivados por los
caños Renegado (K54), Aguas Negras (K39) y Clarín (k13), costado oriental, también
aumentan significativamente, registrando valores pico de alrededor de 96 m 3/s, 305 m3/s,
y 245 m3/s, respectivamente; magnitudes que superan el caudal de 30 m3/s proyectado
para ser derivado por cada caño. (CORPAMAG-DNP-PROCIÉNAGA, 1993).
Volumen de
intercambio
hídrico
Figura 5-1: Dinámica general de los flujos de agua entre el río Magdalena, las ciénagas
circundantes y los caños de interconexión en la zona de estudio.
134 Estudio a escala regional de los componentes hidrológicos e hidráulicos del complejo
cenagoso del Bajo Magdalena ente Calamar y Bocas de Ceniza con apoyo de
modelación matemática y SIG
Figura 5-2: Perfil longitudinal y curva de remanso sobre el río Magdalena diciembre 21 de
2010 (fase La Niña).
Intercambios hídricos río Magdalena - Ciénagas e identificación e zonas susceptible a 135
inundación entre Calamar y Bocas de Ceniza
Vista a la zona de
estudio. Diques y
estructuras en las
ciénagas, El Paraiso,
Santo Tomás,
Luisa/Manatí.
Departamento del
Atlántico
Estructuras de
regulación de caudales
en la Ciénaga Santo
Tomás. Departamento
del Atlántico.
Estructuras de
regulación de caudales
entre la Ciénaga Paraiso
– Ciénaga Santo Tomás.
Departamento del
Atlántico.
Figura 5-3: Visita a la zona de estudio. Diques y estructuras de control alrededor de las
ciénagas. Complejo cenagoso del Departamento del Atlántico. Sobre imagen Google
Earth.
136 Estudio a escala regional de los componentes hidrológicos e hidráulicos del complejo
cenagoso del Bajo Magdalena ente Calamar y Bocas de Ceniza con apoyo de
modelación matemática y SIG
Figura 5-4: Dinámica general del flujo entre el río Magdalena y los caños Renegado (K54),
Aguas negras (K39) y Clarín (K13) - Escenario fase La Niña del ENOS.
Intercambios hídricos río Magdalena - Ciénagas e identificación e zonas susceptible a 137
inundación entre Calamar y Bocas de Ceniza
Con relación a las ciénagas, bajo este mismo escenario (fase La Niña), en la CGSM
predomina considerablemente la entrada de caudales, es decir, sumatoria de caudales
entrantes menos sumatoria de caudales de salida positiva (ΔQ=ƩQentrada – ƩQsalida),
destacándose el mes de diciembre de 2010, acorde con lo expresado anteriormente; algo
similar ocurre en el complejo cenagoso Pajaral como se aprecia en la Figura 5-5. Cabe
anotar que el complejo Pajaral es alimentado desde el río Magdalena por el caño Aguas
Negras al oriente y se conecta con la CGSM al occidente, mientras que el caño Clarín llega
a la CGSM por el costado nororiental. La CGSM a su vez se conecta con el Mar Caribe al
norte y también con el complejo Pajaral al occidente, además de contar con los aportes
reportados de la escorrentía superficial y los caudales de los ríos Aracataca y Fundación
procedentes de la SNSM. Estos balances son extraídos directamente de la herramienta
HEC-RAS para cada cuerpo de agua y/o estructura de derivación de caudal (vertedero
lateral).
Entradas (+)
Salidas (-)
Para la ciénaga Buenavista, que recibe los aportes del río Magdalena a través del caño
Renegado el comportamiento es más cíclico, como se aprecia en la Figura 5-6, donde
trimestralmente el caudal cambia de sentido río-ciénaga a ciénaga-río. Los caudales netos
máximos entrantes a la ciénaga se presentan en el mes de mayo de 2011 (49.98 m3/s),
mientras que las salidas máximas son cercanas a los 57 m3/s (febrero de 2011). La ciénaga
San Antonio (K100), ubicada también al costado oriental del área de estudio, pero más
aguas arriba, tiene un comportamiento similar al descrito para la ciénaga Buenavista,
siendo la que presenta mayores picos de caudales de intercambio, con salidas de
alrededor de 180 m3/s (octubre 2010) y máximos entrantes de 138 m3/s (junio 2011), ver
Figura 5-6, donde también se representan otras ciénagas del costado oriental del área de
estudio.
Entradas (+)
Salidas (-)
Figura 5-6: Caudales netos de intercambio en las ciénagas Buenavista, Tiogollo, San
Antonio - Escenario fase La Niña del ENOS.
para estos contrastes de caudal se asocia con las interconexiones existentes entre las
ciénagas, ver Figura 5-7. Nótese también en esta figura, que algunas ciénagas de este
costado generan muy poco aporte de caudal y/o recepción del mismo, reciben el caudal,
lo transitan y transportan a las ciénagas vecinas o nuevamente al río Magdalena a través
de las estructuras de conexión con el río Magdalena y entre ciénagas, lo cual, también
regula la entrada de caudales a las mismas. En el anexo 2, se presentan los intercambios
hídricos de caudales para los días 21 de diciembre de 2010 (máximo perfil de flujo en La
Niña 1) y 28 de diciembre de 2011 (máximo perfil de flujo en La Niña 2).
Entradas (+)
Salidas (-)
Figura 5-7: Caudales netos en las ciénagas costado occidental del río Magdalena -
Escenario fase La Niña del ENOS.
De otra parte, los resultados del modelo hidráulico para la fase El Niño del ENOS, se
resumen en el perfil de flujo o curva de remanso presentado en la Figura 5-8, el cual detalla
el comportamiento del río Magdalena en época de caudales bajos. En comparación con la
Figura 5-2, es claro que el nivel de la lámina de agua no supera los diques artificiales,
principalmente sobre la margen izquierda (Atlántico) y en algunos casos los albardones o
diques naturales de la margen derecha de la zona de estudio, con lo cual se configura un
140 Estudio a escala regional de los componentes hidrológicos e hidráulicos del complejo
cenagoso del Bajo Magdalena ente Calamar y Bocas de Ceniza con apoyo de
modelación matemática y SIG
sistema desconectado en algunos sitios. El período de análisis para esta fase se presenta
entre los trimestres de JJA de 2009 y MAM de 2010, mientras que en el anexo 2 se detalla
el intercambio hídrico en tiempo seco, tomando como referencia el día 16 de febrero de
2010, la misma fecha de la imagen del perfil o curva de remanso presentada en la Figura
5-8.
Figura 5-8: Perfil longitudinal sobre el río Magdalena febrero 16 de 2010 (Fase El Niño).
Intercambios hídricos río Magdalena - Ciénagas e identificación e zonas susceptible a 141
inundación entre Calamar y Bocas de Ceniza
La Figura 5-9 presenta los caudales derivados del río Magdalena a través de los caños
Renegado (K54), Aguas Negras (K39) y Clarín (K13) durante este período influenciado por
la Fase Niño del ENOS. El caño Renegado tiene un comportamiento cíclico del caudal,
primero en dirección hacia la ciénaga Buenavista, con un máximo alrededor de 50 m 3/s y
posteriormente, se invierte el flujo, esta vez en dirección ciénaga – río Magdalena con un
valor máximo de alrededor de 70 m3/s; el Caño Aguas Negras en este período muestra un
máximo de 92 m3/s en dirección hacia las ciénagas y zonas circundantes; mientras que,
por el caño Clarín se deriva hacia la CGSM un caudal máximo de 65 m3/s. Caudales
notablemente menores a los máximos reportados en la fase La Niña, previamente
analizada.
Figura 5-9: Dinámica general del flujo entre el río Magdalena y los caños Renegado (K54),
Aguas negras (K39) y Clarín (K13) - Escenario fase El Niño del ENOS.
En el caso de las ciénagas principales del costado oriental (departamento del Magdalena),
el complejo Pajaral y CGSM, el caudal neto (ΔQ=ƩQentrada – ƩQsalida), presenta un pico
máximo de caudal entrante de 8 m3/s y 32 m3/s, respectivamente; mientras que la mayor
salida de caudal neto es de 3 m3/s y 15 m3/s, respectivamente. (ver Figura 5-10).
142 Estudio a escala regional de los componentes hidrológicos e hidráulicos del complejo
cenagoso del Bajo Magdalena ente Calamar y Bocas de Ceniza con apoyo de
modelación matemática y SIG
Entradas (+)
Salidas (-)
Figura 5-10: Caudales netos en la CGSM y el complejo cenagoso Pajaral - Escenario fase
El Niño del ENOS.
En la Figura 5-12, se presentan los caudales netos de intercambio de las ciénagas del
costado occidental del río Magdalena, subrayando que son menores en términos generales
a las magnitudes de caudales netos que ocurren en las ciénagas del costado oriental. En
este sector se destacan las magnitudes de los caudales en las ciénagas Malambo y el
complejo cenagoso el Uvero, ubicadas respectivamente en los extremos aguas arriba y
aguas abajo de este costado; nótese que las magnitudes de caudal máximo entrante no
superan los 8 m3/s, mientras que los caudales máximos de salida no superan los 12 m3/s,
también, bastante menores a los caudales reportados para la fase La Niña (ENOS).
Intercambios hídricos río Magdalena - Ciénagas e identificación e zonas susceptible a 143
inundación entre Calamar y Bocas de Ceniza
Entradas (+)
Salidas (-)
Figura 5-11: Caudales netos en las ciénagas Buenavista, Tiogollo, San Antonio -
Escenario fase El Niño del ENOS.
Entradas (+)
Salidas (-)
Figura 5-12: Caudales netos en las ciénagas costado occidental del río Magdalena -
Escenario fase El Niño del ENOS.
144 Estudio a escala regional de los componentes hidrológicos e hidráulicos del complejo
cenagoso del Bajo Magdalena ente Calamar y Bocas de Ceniza con apoyo de
modelación matemática y SIG
En la Tabla 5-1, se resumen los volúmenes de intercambio para las ciénagas de la zona
de estudio, en las fases año neutro, El Niño (2009-2010) y La Niña (2010-2012) del ENOS
para el período modelado hidráulicamente. En color azul los volúmenes hacia las ciénagas
o ganancia de caudales y en color rojo los volúmenes salientes de la mismas o pérdida de
caudales. Nótese que para el período año neutro y El Niño, junio (2009) – mayo (2010), se
presenta una disminución de los caudales entrantes a las ciénagas en comparación con
los períodos La Niña. En total en este período sale de las ciénagas un volumen aproximado
de 12 Mm3 (millones de metros cúbicos). De otra parte, para los períodos La Niña, años
2010 – 2011 (Niña 1) y años 2011 – 2012 (Niña 2), es apreciable que los volúmenes de
entrada en la mayoría de ciénagas aumentan, se encuentran estas directamente
conectadas o no con el río Magdalena, mostrando mayores valores de intercambio en el
primer período (La Niña 1), junio (2010) – Mayo (2011), con 627 Mm3, este volumen llena
las ciénagas, ya que estas se encuentran precedidas del tiempo seco (El Niño); para el
segundo período (La Niña 2), junio (2011) – febrero (2012), la magnitud de los intercambios
disminuye en comparación con el período inmediatamente anterior, debido a que las
ciénagas se encontrarían llenas por el período La Niña 1. Sin embargo, siguen
prevaleciendo los flujos aportantes hacia estas, con un volumen neto total aproximado de
184 Mm3. Como es apreciable también en la Tabla 5-1, los mayores cuerpos receptores
de volúmenes de agua son las ciénagas San Antonio (K100), Buenavista (K54), complejo
Pajaral (K39) y la CGSM, en el costado oriental; también llama la atención en este costado
los volúmenes entrantes a las ciénagas El Piñón y Tiogollo, los cuales se deben a que
estas reciben los aportes de la ciénaga San Antonio. En la margen izquierda (departamento
del Arlántico) el mayor receptor de caudales es la ciénaga Malambo ubicada aguas abajo
(K20), la cual está conectada totalmente con río Magdalena, además, se conecta con la
ciénaga El Convento, que es la segunda mayor receptora de caudales en este costado.
En todo el sistema analizado, también llama la atención, que debido a la topografía la
ciénaga identificada como 3255, perteneciente al complejo cenagoso Pajaral
(departamento del Magdalena), queda completamente desconectada, sin recepción ni
aporte de agua.
Intercambios hídricos río Magdalena - Ciénagas e identificación e zonas susceptible a 145
inundación entre Calamar y Bocas de Ceniza
Tabla 5-1: Resumen volúmenes de intercambio en las 31 ciénagas del área de estudio. En
color rojo los volúmenes de salida de las ciénagas y viceversa en color azul.
entre el período DEF 2011 – JJA 2011, y posteriormente el período JJA 2011 – DEF 2012.
Las ciénagas del costado occidental presentan una dinámica de intercambios hídricos
similar a los descritos para el costado oriental, pero con menores magnitudes.
En todo el macrosistema se destaca al final del período de modelación que las ciénagas
presentan un superávit o recepción de agua fresca desde el río Magdalena. Lo anterior,
construye una evidencia de lo importante y significativo de las conexiones naturales entre
el río Magdalena y las ciénagas circundantes, con comportamiento multianual, las cuales
se deben preservar. Mientras que, los caudales de agua derivados a través de los caños
Renegado, Aguas Negras y Clarín, conducen desde el río Magdalena importantes
volúmenes de agua que tienen como receptor final la CGSM. Estos volúmenes de agua
garantizan el adecuado funcionamiento del sistema CGSM, como ha sido detallado
también por (Díaz-Granados, Camacho, et al., 1996) y (Múnera, Vélez, et al., 2003a),
requiriendo mantenimiento y adecuada información sobre las condiciones actuales de los
mismos. La Figura 5-14, es un infograma resumen de los resultados de intercambio hídrico
más representativos del Bajo Magdalena entre Calamar y Bocas de Ceniza.
Figura 5-13: Volúmenes de intercambio hídrico en las ciénagas del Bajo Magdalena
período JJA/2009 – DEF/2012.
148
VOLUMEN / TRIMESTRAL
NEUTRO / EL NIÑO (2009-2010) LA NIÑA 1 (2010 - 2011) LA NIÑA 2 (2011 - 2012)
ID CIÉNAGAS Margen / Departamento JJA-2009 SON-2009 DEF-2009/2010MAM-2010 JJA-2010 SON-2010 DEF-2010/2011 MAM-2011 JJA-2011 SON-2011 DE-2011/2012 TOTALES
Mm3 Mm3 Mm3 Mm3 Mm3 Mm3 Mm3 Mm3 Mm3 Mm3 Mm3 Mm3
3254 Derecha / Magdalena -0.19 0.36 -0.08 -0.05 0.61 0.45 -0.12 0.61 0.62 1.06 -0.17 3.11
3255 Derecha / Magdalena 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00
3256 Derecha / Magdalena -0.06 0.52 -0.05 0.02 1.40 0.74 0.05 0.68 0.58 1.09 -0.28 4.69
3257 Derecha / Magdalena -0.31 0.60 -0.12 -0.09 1.02 0.73 -0.18 1.00 1.03 1.72 -0.27 5.14
3258 Derecha / Magdalena -0.01 0.09 -0.02 0.00 0.21 0.09 -0.01 0.13 0.07 0.18 -0.05 0.69
3259 Derecha / Magdalena -1.11 1.00 -1.16 1.29 4.91 0.92 -6.72 5.59 -2.46 2.59 -0.95 3.89
3260 Derecha / Magdalena 4.96 1.55 -1.56 1.47 6.81 1.25 -9.37 7.74 -3.52 3.62 -1.29 11.65
3261 Derecha / Magdalena -6.26 7.16 -18.54 17.75 18.34 3.82 -30.29 24.55 -10.33 10.35 -3.42 13.13
3262 Derecha / Magdalena -8.86 11.54 -27.13 25.49 26.94 3.33 -43.86 36.62 -14.54 14.57 -3.69 20.40
de volúmenes de agua a las ciénagas.
3263 Derecha / Magdalena 21.63 4.90 -4.69 4.52 21.23 4.50 -29.21 24.04 -11.05 11.28 -4.40 42.75
3264 Derecha / Magdalena -2.11 1.91 -2.19 2.47 9.21 1.37 -12.63 10.51 -4.63 4.84 -1.62 7.14
3265 Derecha / Magdalena -1.06 0.93 -1.07 1.21 4.54 1.00 -6.18 5.21 -2.41 2.43 -0.92 3.68
3280 Derecha / Magdalena -0.27 0.82 -0.34 0.20 2.36 2.02 -1.70 1.54 -0.88 0.79 -0.94 3.59
3282-CGSM Derecha / Magdalena -9.03 19.29 -1.40 -5.02 28.12 18.35 4.49 28.28 39.34 53.49 -7.85 168.06
3283 Derecha / Magdalena -0.02 0.05 -0.03 0.02 0.19 0.17 -0.22 0.15 -0.12 0.03 -0.11 0.11
3284 Derecha / Magdalena -3.25 4.17 -9.14 8.62 9.86 1.10 -16.02 13.43 -5.28 5.33 -1.18 7.63
Buenavista Derecha / Magdalena -39.66 47.67 -117.63 112.95 116.27 11.27 -187.99 157.43 -62.16 62.37 -12.17 88.36
Complejo Uvero Izquierda / Atlántico -7.74 5.00 -19.40 18.14 14.83 1.02 -23.27 19.26 -8.01 7.86 -1.09 6.60
El convento Izquierda / Atlántico -0.51 -0.11 -0.07 0.00 10.78 2.54 -13.29 9.52 -5.96 6.39 -1.69 7.61
Luisa-Manatí Izquierda / Atlántico 0.00 0.00 0.00 0.00 6.24 1.33 -8.19 4.61 -4.60 5.80 -1.04 4.16
Malambo Izquierda / Atlántico -9.76 5.19 -16.61 11.75 38.59 9.25 -57.28 46.19 -20.74 20.90 -5.32 22.17
Complejo Pajaral Derecha / Magdalena -2.41 4.51 -0.97 -0.56 7.70 5.42 -1.58 7.67 7.65 13.09 -2.02 38.49
Pariso Izquierda / Atlántico 0.00 0.00 0.00 0.00 3.28 0.68 -4.29 2.42 -2.46 3.06 -0.55 2.14
Pariso 2 Izquierda / Atlántico 0.00 0.00 0.00 0.00 0.88 0.20 -1.16 0.65 -0.68 0.85 -0.16 0.58
El Piñon Derecha / Magdalena -13.87 10.38 -34.88 34.07 42.21 -9.48 -53.63 51.35 -13.04 29.96 -4.78 38.28
Sabanagrande Izquierda / Atlántico -1.29 0.70 -1.12 0.86 3.60 0.62 -5.07 4.24 -1.90 1.87 -0.30 2.23
Sabanagrande 2 Izquierda / Atlántico -0.84 0.44 -0.29 0.42 2.67 0.45 -3.56 2.97 -1.39 1.37 -0.28 1.97
San Antonio Derecha / Magdalena -105.38 75.80 -252.54 249.95 297.25 -53.19 -437.90 371.62 -174.25 164.62 -8.42 127.55
Santo Tomás Izquierda / Atlántico -3.24 1.95 -1.38 1.93 9.56 1.26 -12.96 10.83 -4.94 4.88 -1.13 6.76
Tiogollo Izquierda / Atlántico -3.97 2.71 -8.11 7.69 50.64 46.73 12.93 27.27 40.52 34.07 41.50 251.96
Tiogollo 2 Izquierda / Atlántico -2.43 1.57 -3.27 3.43 -9.92 -27.80 -41.84 -18.24 1.15 7.43 -4.64 -94.55
TOTALES -197.04 210.73 -523.78 498.53 730.34 30.14 -991.07 857.88 -264.41 477.89 -29.23 799.99
cenagoso del Bajo Magdalena ente Calamar y Bocas de Ceniza con apoyo de
modelación matemática y SIG
Estudio a escala regional de los componentes hidrológicos e hidráulicos del complejo
color rojo indica la salida de volúmenes de agua de las ciénagas y el color azul la entrada
Tabla 5-2: Volúmenes de intercambio trimestral en las 31 ciénagas del área de estudio. El
Intercambios hídricos río Magdalena - Ciénagas e identificación e zonas susceptible a 149
inundación entre Calamar y Bocas de Ceniza
-
Figura 5-14: Volúmenes de intercambio hídrico en las ciénagas del Bajo Magdalena
período JJA/2009 – DEF/2012 y condiciones de frontera (entrada y salida).
150 Estudio a escala regional de los componentes hidrológicos e hidráulicos del complejo
cenagoso del Bajo Magdalena ente Calamar y Bocas de Ceniza con apoyo de
modelación matemática y SIG
Figura 5-15: Volúmenes de intercambio hídrico en las ciénagas del Bajo Magdalena
período JJA/2009 – DEF/2012 y condiciones de frontera aguas abajo aumentada 1.0 m.
Intercambios hídricos río Magdalena - Ciénagas e identificación e zonas susceptible a 153
inundación entre Calamar y Bocas de Ceniza
Caudal (m3/s)
NIvel (m)
NIvel (m)
2
14000 2
14000
1
12000 0
Tiempo 12000 0
Tiempo (48 horas)
Nivel TR 50 Caudal TR 50
Caudal TR 100 Nivel TR 100
En las Figura 5-18 a Figura 5-24, se presentan como apoyo a los resultados descritos
anteriormente, una serie de mapas detallados con la mancha de inundación TR100 años
sobre imagen Google Earth de la zona de estudio, resaltando las zonas susceptibles a
inundación en este período de retorno. Los colores más oscuros representan mayor
profundidad en la lámina de agua y los colores más claros una menor profundidad de agua.
Cabe anotar, que se debe tener especial cuidado con la interpretación de estas imágenes
dado que la información topográfica limita la extensión de la mancha de inundación y puede
enmascarar el riesgo hacia las poblaciones allí asentadas.
SOLEDAD
GALAPA
CGSM
1,700,000
1,700,000
± TUBARA
MALAMBO
JUAN DE ACOSTA
BARANOA SABANAGRANDE
SITIONUEVO
POLONUEVO SANTO TOMAS
REMOLINO
EL RETEN
SABANALARGA PONEDERA
LURUACO
Rio Magdalena
ARACATACA
inundación entre Calamar y Bocas de Ceniza
FUNDACION
REPELON
SALAMINA
CANDELARIA
PIVIJAY
MANATI
1,650,000
1,650,000
SUAN
SANTA LUCIA CERRO DE SAN ANTONIO
Calamar PEDRAZA
lámina de agua_TR100
cabeceras
Zona inundable_geomorfología 0 5 10 20 30 40
Inundación Canal del Dique 2011
de este, claramente identificable hidráulica y geomorfológicamente. Ver Figura 5-17.
Km
Modelo de sombras total.img
occidental (departamento del Atlántico) y que constituyen por ende una zona de descanso
155
156 Estudio a escala regional de los componentes hidrológicos e hidráulicos del complejo
cenagoso del Bajo Magdalena ente Calamar y Bocas de Ceniza con apoyo de
modelación matemática y SIG
Profundidad de agua
Figura 5-18: Inundación caudal Tr 100 años en el Bajo Magdalena, sobre imagen Google
Earth.
a) b)
Figura 5-19: Inundación caudal Tr 100 años en el Bajo Magdalena, Nivel de agua hasta
de 1.00 m en el Municipio de Suan (K100) en la margen occidental. Al costado oriental
desborde de la conexión entre ciénagas San Antonio (K100) y El Piñón (K85). b) Registro
fotográfico municipio de Suan (26/01/2011) en la ocurrencia de la Fase Niña del ENOS en
2011. Parte baja de la zona de estudio. Fuente: IDEHA - Gobernación del Atlántico, 2011.
Intercambios hídricos río Magdalena - Ciénagas e identificación e zonas susceptible a 157
inundación entre Calamar y Bocas de Ceniza
Profundidad de agua
Figura 5-20: Inundación caudal Tr 100 años en el Bajo Magdalena, Nivel de agua hasta
de 1.30 m en el Municipio de El Piñón (K85) en la margen oriental.
Profundidad de agua
Figura 5-21: Inundación caudal Tr 100 años en el Bajo Magdalena, Nivel de agua hasta
de 1.30 m en el Municipio de Salamina (K80) en la margen oriental (Departamento de
Magdalena).
158 Estudio a escala regional de los componentes hidrológicos e hidráulicos del complejo
cenagoso del Bajo Magdalena ente Calamar y Bocas de Ceniza con apoyo de
modelación matemática y SIG
Profundidad de agua
Figura 5-22: Inundación caudal Tr 100 años en el Bajo Magdalena, riesgo de desborde
aguas arriba del Municipio de Remolino (K40) en la margen oriental (Departamento de
Magdalena). Ciénaga Paraíso en su máxima capacidad.
Profundidad de agua
Figura 5-23: Inundación caudal Tr 100 años en el Bajo Magdalena, riesgo de desborde
aguas arriba del Municipio de Sitio Nuevo (K34) en la margen oriental (Departamento de
Magdalena). Ciénaga Paraíso en su máxima capacidad.
Intercambios hídricos río Magdalena - Ciénagas e identificación e zonas susceptible a 159
inundación entre Calamar y Bocas de Ceniza
Profundidad de agua
Figura 5-24: Inundación caudal Tr 100 años en el Bajo Magdalena, extensión del complejo
cenagoso del Atlántico entre el K40 - K20 las ciénagas: Paraíso, Santo Tomás,
Sabanagrande, El convento y Malambo.
Para el caudal con período de retorno de 50 años, la afectación hacia el departamento del
Magdalena persiste en el canal Tiogollo y los municipios de El Peñón. Salamina, Remolino
y Sitio Nuevo. Los caños Clarín, Aguas Negras y Renegado presentan también
desbordamientos por falta de capacidad. Mientras que, hacia el departamento del
Atlántico, se sigue presentando el riesgo de inundación en las poblaciones de Suan,
Palmar de Varela, Santo Tomás, y Sabanagrande, esta vez con disminución aproximada
de 0.05 m. en la lámina de agua en comparación con el TR 100 años, lo cual se explica
por la variación en el nivel de la condición de frontera aguas abajo. En la Figura 5-25 la
mancha de inundación y en la Figura 5-26 las cotas de agua en el sitio de estudio, ambos
con período de retorno de 50 años. Ver también el anexo 3.
900,000 950,000 1,000,000
160
SOLEDAD
± TUBARA
GALAPA
CGSM
1,700,000
1,700,000
MALAMBO
JUAN DE ACOSTA
BARANOA SABANAGRANDE
SITIONUEVO
POLONUEVO SANTO TOMAS
REMOLINO
EL RETEN
SABANALARGA PONEDERA
LURUACO
Rio Magdalena
SANTA CATALINA ARACATACA
CLEMENCIA
FUNDACION
REPELON
SALAMINA
CANDELARIA
VILLANUEVA PIVIJAY
MANATI
1,650,000
1,650,000
SUAN
SANTA LUCIA CERRO DE SAN ANTONIO
Calamar PEDRAZA
Leyenda
lámina de agua_TR50 900,000 950,000 1,000,000
cabeceras
Zona inundable_geomorfología
Inundación Canal del Dique 2011
0 5 10 20 30 40
cenagoso del Bajo Magdalena ente Calamar y Bocas de Ceniza con apoyo de
modelación matemática y SIG
Estudio a escala regional de los componentes hidrológicos e hidráulicos del complejo
Low : 0
Intercambios hídricos río Magdalena - Ciénagas e identificación e zonas susceptible a 161
inundación entre Calamar y Bocas de Ceniza
Se debe así tener prudencia con los resultados reportados principalmente, por las
limitantes de la modelación, donde se destacan la falta de una mayor serie de datos
consistentes entre las condiciones de entrada aguas arriba y la condición de control aguas
abajo, caudales hacia los caños, medición de niveles en ciénagas, topobatimetría de todas
las ciénagas y batimetría actualizada de la franja de inundación del río, con lo cual se
podría realizar una mejor validación de los resultados. Sin desmedro de esto, los resultados
164 Estudio a escala regional de los componentes hidrológicos e hidráulicos del complejo
cenagoso del Bajo Magdalena ente Calamar y Bocas de Ceniza con apoyo de
modelación matemática y SIG
de manchas aquí reportados deben ser tomados como una importante referencia, dado
que integran, a escala regional, el análisis de todo el hidrosistema de forma integral.
6. Conclusiones y recomendaciones
Para la calibración y validación del modelo hidráulico se han utilizado datos de imágenes
satelitales clasificadas para la cobertura agua. Si bien es cierto por las limitaciones de
información y por el desafío que significa abordar el análisis regional de este enorme
hidrosistema, los resultados estrictamente no resultan ser totalmente concluyentes, sí son
166 Estudio a escala regional de los componentes hidrológicos e hidráulicos del complejo
cenagoso del Bajo Magdalena ente Calamar y Bocas de Ceniza con apoyo de
modelación matemática y SIG
Intercambio hídrico.
Mm3), (+7.61 Mm3), (6.76 Mm3) y (+6.60 Mm3), respectivamente; priman caudales de
entrada desde el río Magdalena, los cuales resultan relevantes de mantener, evitando que
los diques alrededor de estas ciénagas las aíslen del sistema río Magdalena – Ciénagas.
La importancia de estos volúmenes de intercambio ha sido ampliamente documentada en
la literatura. A modo de ejemplo los impactos medidos en la salinidad de agua en el
complejo cenagoso Pajaral y la CGSM, antes y después de la reactivación de los caños
Clarín y Renegado, indican que los valores de salinidad han venido disminuyendo como
resultado de los mayores intercambios (ver Figura 6-1) (INVEMAR, 2013). Esto confirma
aún más la importancia de la modelación integral aquí realizada, la cual, se puede seguir
alimentando en función de la mejor información generada por las entidades que tienen
jurisdicción en el macrosistema.
Figura 6-1: Salinidad intersticial (0.5 m) promedio registrada en los suelos de manglar de
cinco estaciones de monitoreo durante el período 1995-2013. Las flechas rojas indican la
apertura de los caños Clarín (año 1996), Aguas Negras y Renegado (año 1998). Fuente:
INVEMAR (2013)
última a partir de un modelo 1D. Los resultados a este respecto son satisfactorios pero
limitados por la calidad de la información topográfica de detalle existente. La dinámica
fluvial indica que las ciénagas circundantes al río Magdalena, principalmente entre el K60
y K13, hacen parte activa del hidrosistema natural y que estas son parte de la llanura aluvial
de inundación, la cual ha sido altamente intervenida por diques y otras estructuras de
regulación, por lo cual deben mantenerse conectadas al río Magdalena, se reitera. Estos
cuerpos de agua con funciones de regulación de caudales hacia aguas abajo son
representados también en los resultados de la modelación matemática implementada para
la zona, para los diferentes períodos de recurrencia, 1 en 100 años y 1 en 50 años.
Además, se han identificado claramente cuatro zonas de desborde para ambos períodos
de retorno, los cuales son el K19 a la altura de la isla Rondón, K34 municipio de Sitio
Nuevo, K40 municipio de Remolino, K80 municipio de Salamina, y el K100 entre la ciénaga
San Antonio y la ciénaga Tiogollo cuya afectación se extiende aguas arriba hacia el
municipio de El Piñón K85, estas en la margen oriental. Hacia la margen occidental el
riesgo más representativo es hacia el municipio de Suan (100), pero deben ser incluidos
en detalle los diques para ser concluyentes. Como resultado final y respondiendo a la
pregunta de investigación planteada, también con la modelación integral aquí presentada
fue posible elaborar una propuesta para la delimitación del complejo cenagoso del costado
oriental, evidenciando que este en algunos casos tiene presencia de infraestructuras, las
cuales están en riesgo de inundación para los períodos de retorno de 100 y 50 años. La
Figura 6-2 resume el caso de la Ciénaga Malambo, mientras que en el anexo 3 se
encuentra la totalidad del complejo cenagoso del departamento del Atlántico.
170 Estudio a escala regional de los componentes hidrológicos e hidráulicos del complejo
cenagoso del Bajo Magdalena ente Calamar y Bocas de Ceniza con apoyo de
modelación matemática y SIG
Figura 6-2: Delimitación propuesta para la ciénaga Malambo (K20) a partir de criterio
hidráulico. Se identifican en color cian la extensión propuesta, 3.4 km a partir del río
Magdalena, en color rojo la afectación en la zona con una cota de agua 6.60 m.s.n.m.
Tr 50 y Tr 100 años. Predios afectados considerados dentro de la ciénaga.
Cabe destacar además que no se ha incluido información sobre flujo sub-superficial (aguas
subterráneas) y aquí habría un amplio espacio para acoplar una modelación aguas
superficiales – aguas subterráneas.
Futuros desarrollos.
A continuación, una serie de propuestas de investigación que son tendientes a lograr
complementar el trabajo aquí expuesto y a lograr mayor entendimiento integral de la
macrorregión.
Figura 6-3: Curva de duración de caudales (CDC) estación Calamar. Fuente: Elaboración
propia.
Figura 6-4: Curva de duración de caudales (CDC) estación Incora K7. Fuente: Elaboración
propia.
Anexo 2. Infogramas de intercambios hídricos. 179
K 13 - CAÑO CLARÍN
11.00 m3/s
CIÉNAGA MALAMBO K 20 0.46%
H= 1.10 m
3
0 m /s
0.00%
CIÉNAGA EL CONVENTO K 29
H= 3.09 m
0 m3/s
0.00%
CIÉNAGA ID CIÉNAGA ID 3257
CIÉNAGA 3254 H= 0.50 m
K 32 0.26 m3/s 1.61 m3/s
SABANAGRANDE
H= 2.00 m 0 m3/s
RÍO MAGDALENA
0.00%
0.17 m3/s 1.34 m3/s
CIÉNAGA
K 34 13.47 m3/s
SABANAGRANDE 2
H= 2.54 m
0 m3/s
0.00% CIÉNAGA GRANDE
CIÉNAGA SANTO TOMÁS CIÉNAGA ID 3255 DE SANTA MARTA
K 37
H= 2.00 m 3 3 H= 0.93 m 3 H = 0.50 m
6.53 m /s 0 m /s 0 m /s
0 m3/s K 39 - CAÑO AGUAS NEGRAS
0.00%
CIÉNAGA ID 3256
3 COMPLEJO CENAGOSO
16.92 m /s H= 0.50 0m
0.70%
PAJARAL 0 m3/s 0 m3/s
H = 0.50 m
3
m3/s
m3/s
-0.19
CIÉNAGA PARAISO K 44
3
0.88
H= 0.51 m 0 m /s H= 0.51 m
H= 3.98 m
0 m3/s
0.00%
CIÉNAGA ID 3284
CAUDAL NETO
CAUDAL NETO
-0.18 m3/s
-19.12 m3/s
CIÉNAGA ID 3263
H= 3.07 m
CIÉNAGA SAN
ANTONIO
CALAMAR - INCORA K7 K 101
Q IN 2408 m3/s
H = 2.88 m
Figura 6-5: Resultado de los intercambios hídricos para el día 16 de febrero de 2010 en la
100%
K 13 - CAÑO CLARÍN
3
245.47 m /s
CIÉNAGA MALAMBO K 20 1.63%
H= 5.29 m
-413.18 m3/s
-2.74%
-406.65 m3/s
CIÉNAGA EL CONVENTO K 29
H= 5.36 m
254.44 m3/s
1.69%
3
-150.80 m /s
CIÉNAGA ID CIÉNAGA ID 3257
3254 H= 0.58 m
CIÉNAGA K 32 H= 0.58 m 2.90 m3/s 15.37 m3/s
SABANAGRANDE
H= 5.66 m -592.61 m3/s
RÍO MAGDALENA
-3.93%
-743.02 m3/s CIÉNAGA GRANDE
DE SANTA MARTA
2.82 m3/s 12.34 m3/s 274.49 m3/s H = 0.55 m
K 34
CIÉNAGA
SABANAGRANDE 2
H= 5.70 m -191.63 m3/s
|
-1.27%
67
33
-924.31 m3/s
m3/s
m3/s
CIÉNAGA PARAISO
H= 5.81 m
3
0 m /s
0.00%
-81.23 m3/s
CIÉNAGA ID 3284 CIÉNAGA ID 3261
CIÉNAGA PARAISO 2 H= 6.55 m H= 6.55 m
K 45 -1.95 m3/s
H= 5.83m
80.92 m3/s K 54 - CAÑO RENEGADO -10.2 m3/s -2.86 m3/s
CIÉNAGA
0.54% BUENAVISTA
H= 6.55 m
3 3
CAUDAL NETO
CAUDAL NETO
3
-1.03 m /s
995.16 m3/s
6.60% CIÉNAGA ID 3263
H= 6.55 m
CIÉNAGA SAN
ANTONIO
CALAMAR - INCORA K7 K 101
Q IN 15085 m3/s
H = 8.29 m
100%
Figura 6-6: Resultado de los intercambios hídricos para el día 21 de diciembre de 2010 en
181
CAUDAL SALIDA MODELO 13661.40 m3/s % CAUDAL REGULADO DEP. ATLÁNTICO 6.5 %
ERROR % 0.25 % % CAUDAL REGULADO DEP. MAGDALENA 3.2 %
INTERCAMBIO HÍDRICO DESDE EL RÍO MAGDALENA EN TIEMPO HÚMEDO (NIÑA 2)
FECHA 28/12/2011
3 3
Q OUT 13221.9 m /s Q OUT 593 m /s
K 13 - CAÑO CLARÍN
3
233.61 m /s
CIÉNAGA MALAMBO K 20 1.60%
H= 5.18 m
3
-336.42 m /s
-2.30%
-334.77 m3/s
CIÉNAGA EL CONVENTO K 29
H= 5.25 m
207.8 m3/s
1.42%
3
-126.75 m /s
CIÉNAGA ID CIÉNAGA ID 3257
3254 H= 0.56 m
CIÉNAGA 3 3
K 32 H= 0.56 m 3.03 m /s 15.33 m /s
SABANAGRANDE
RÍO MAGDALENA
H= 5.55 m 3
-524.36 m /s
-3.58%
-651.01 m3/s CIÉNAGA GRANDE
3 3 3
DE SANTA MARTA
3.52 m /s 13.66 m /s 275.97 m /s
CIÉNAGA K 34
H = 0.54 m
SABANAGRANDE 2
H= 5.58 m
-162.84 m3/s |
-1.11%
63
37
-813.77 m3/s
m3/s
m3/s
904
K 44 CIÉNAGA ID 3280 CIÉNAGA ID 3283
CIÉNAGA PARAISO
H= 0.56 m H= 0.56 m
DEPARTAMENTO DEL MAGDALENA
0.01 m3/s
527.74
CAUDAL NETO
CAUDAL NETO
-0.02 m3/s
962.06 m3/s
6.58% CIÉNAGA ID 3263
H= 6.41 m
CIÉNAGA SAN
ANTONIO
CALAMAR - INCORA K7 K 101
Q IN 14627 m3/s
H = 8.13 m
100%
Figura 6-7: Resultado de los intercambios hídricos para el día 28 de diciembre de 2011 en
3
CAUDAL A CIÉNAGAS 6928.43 m /s TODAS LAS CIÉNAGAS CONECTADAS E INTERCAMBIANDO CAUDAL DESDE Y HACIA EL RÍO MAGDALENA
CAUDAL A RÍO MAGDALENA -5496.69 m3/s CAUDAL TOTAL REGULADO 1432 m3/s
CAUDAL SALIDA 13195.26 m3/s % TOTAL CAUDAL REGULADO 9.8 %
CAUDAL SALIDA MODELO 13221.90 m3/s % CAUDAL REGULADO DEP. ATLÁNTICO6.2 %
cenagoso del Bajo Magdalena ente Calamar y Bocas de Ceniza con apoyo de
modelación matemática y SIG
Estudio a escala regional de los componentes hidrológicos e hidráulicos del complejo
Figura 6-10: Propuesta de delimitación de ciénagas del costado occidental a partir del
criterio conjunto geomorfología y modelación matemática TR 100 años.
880000 900000 920000 940000 960000 980000 1000000 1020000
LEYENDA
ESTACIIONES HIDROLÓGICAS
Cat
MAREOGRAFO $ LG_ EXISTENTE
! ! LG_ PROYECTADO
1720000
1720000
LM _ EXISTENTE
Anexo 3. Mapas
± "
)
CONERTUA DE AGUA
CANO CLARÍN MAREOGRAFO SANTA MARTHA TIPO
DÁRSENA - BARR
! ! Cuerpos de agua
"
) Inundacion
C. PAJARAL 1 Zonas inundables periodicamente
complejo Cenaoso MALAMBO ! CGSM spot_642_328_050215_m_b.img
1700000
1700000
del Atlántico !
!C. SABANAGRANDE RGB
so
C. LUISA MANATÍ C. PAJARAL 5 C. PAJARAL 4
! Red: Layer_1
1680000
1680000
! ! !
CAÑO AGUA NEGRAS
n ao
e
C. PAJARAL 10 Green: Layer_2
! Blue: Layer_3
jo C
C. PAJARAL!8 TREBOL EL
RÍO MAGDALENA K50 C. PAJARAL 7 Modelo de sombras total.img
C. PAJARAL !
9
ple
!C. PAJARAL 6 $
PTE FERROCARRIL Value
! Hi g h : 2 5 4
com
C.!UVERO ! "
)
! CAÑO RENEGADO Lo w : 0
er o
FUNDACION
1660000
1660000
Uv
"
)
1640000
1640000
INCORA K-7
SAN PEDRITO $ !C. SAN ANTONIO
CALAMAR
"
)
"
)
1620000
1620000
Bogota.
CORMAGDALENA (2007) Atlas de la cuenca Magdalena-Cauca, Bogotá D.C.
CORMAGDALENA-UNAL (2007) Estudios e investigaciones de las obras de restauración
ambiental y de la navegación del Canal del Dique, Bogotá D.C.
CORMAGDALENA-UNINORTE (2007) Actualización de la caracterización hidráulica y
sedimentológica del río Magdalena para los estudios e investigaciones de el canal
navegable de acceso al puerto de Barranquilla entre Bocas de Ceniza (K0) y Puente
Pumarejo (K22).
CORPAMAG-DNP-PROCIÉNAGA (1993) Plan de recuperación del complejo lagunar de
la Ciénaga Grande de Santa Marta - Análisis de alternativas, Santa Marta.
CRA (2007) Documentación del estado de las cuencas hidrográficas en el departamento
del Atlántico,
Díaz-Granados, M., Camacho, L. A., and Maestre, A. (1996) Modelación de Balances
Hídricos de ciénagas fluviales y costeras colombianas. Revista de Ingeniería
(versión electrónica) No1., 2–20.
DNP-UNAL (2012) Estudios, análisis y recomendaciones para el ordenamiento ambinetal
y el desarrollo territorial de la Mojana, Bogotá D.C.
Dutta, D., Teng, J., Vaze, J., Lerat, J., Hughes, J., and Marvanek, S. (2013) Storage-
based approaches to build floodplain inundation modelling capability in river system
models for water resources planning and accounting. Journal of Hydrology, 504, 12–
28. [online] http://www.sciencedirect.com/science/article/pii/S0022169413006860
(Accessed June 2, 2015).
FEMA (2009) Federal Emergency Management Agency - Flood maps., [online]
//msc.fema.gov.
Figueroa Ortiz, C. A. (2016) Caracterización del campo de precipitación sobre la
microcuenca urbana del campus Universidad Nacional de Colombia - sede Bogotá.
Gao, H., Birkett, C., and Lettenmaier, D. P. (2012) Global monitoring of large reservoir
storage from satellite remote sensing data products in comparison with gauge
observations for the five largest reservoirs in the. , 48(February), 1–12.
García-Pintado, J., Mason, D. C., Dance, S. L., Cloke, H. L., Neal, J. C., Freer, J., and
Bates, P. D. (2015) Satellite-supported flood forecasting in river networks: A real
case study. Journal of Hydrology, 523, 706–724. [online]
http://www.sciencedirect.com/science/article/pii/S0022169415001031 (Accessed
June 2, 2015).
Hingray, B., Picouet, C., and Musy, A. (2014) Hydrology:A Science for Engineers,
Horritt, M. S. and Bates, P. D. (2001) Effects of spatial resolution on a raster based model
of ¯ ood ¯ ow. , 253.
Horritt, M. S. and Bates, P. D. (2002) Evaluation of 1D and 2D numerical models for
predicting river flood inundation. , 268, 87–99.
Bibliografía 191
Hunter, N. M., Bates, P. D., Horritt, M. S., and Wilson, M. D. (2007) Simple spatially-
distributed models for predicting flood inundation: A review. Geomorphology, 90(3–
4), 208–225.
IDEAM (1995) Regionalización de crecientes máximas, Bogotá D.C.
IDEHA-Cormagdalena (2012) Evolución morfológica al comportamiento de las orillas y
canal navegable utilizando el modelo MIKE -21C, informe final. , Informe fi, 98.
INVEMAR (2013) Monitoreo de las condiiones ambientales y los cambios estructurales y
funcinales de las comunidades vegetales y de los recuros pesqueros durante la
rehabilitación de la Ciénaga Grande de Santa Marta. , 163.
Javaheri, A. and Babbar-Sebens, M. (2014) On comparison of peak flow reductions, flood
inundation maps, and velocity maps in evaluating effects of restored wetlands on
channel flooding. Ecological Engineering, 73, 132–145. [online]
http://www.sciencedirect.com/science/article/pii/S0925857414004194 (Accessed
June 2, 2015).
MADS (2013) Guía técnica para la formulación de los Planes de Ordenación y Manejo de
Cuencas Hidrográficas POMCAS, Bogotá D.C.
MARM (2011) Guía metodológica para el desarrollo del Sistema Nacional de Cartografía
de Zonas Inundables (Ministerio de Medio Ambiente y Medio Rural y Marino, ed.),
Madrid.
Mason, D. C., Bates, P. D., and Amico, J. T. D. (2009) Calibration of uncertain flood
inundation models using remotely sensed water levels. Journal of Hydrology, 368(1–
4), 224–236. [online] http://dx.doi.org/10.1016/j.jhydrol.2009.02.034.
MAVDT (2010) Política Nacional para la Gestión del Recurso Hídrico (G. de R. H.
Ministerio de Ambiente, Vivienda y Desarrollo Territorial, Viceministerio de Ambiente,
Dirección de Ecosistemas, ed.), Bogotá D.C. [online]
http://www.minambiente.gov.co/documentos/5774_240610_libro_pol_nal_rec_hidric
o.pdf.pdf.
Menéndez, M. and Núñez, V. (2009) Uso de los Sensores Remoros en los Recurcos
Naturales, Salta, Universidad Nacional de Salta - Argentina.
Mesa, O. J. et al. (1997) Introcucción al clima de Colombia, Medellín, Universidad
Nacional de Colombia sede Medellín.
Múnera, J. C., Vélez, J. I., Poveda, G., Posada, J. E., Montoya, J. D., and Cardona, Y. M.
(2003a) Dinámica hidrológica de la ciénaga grande de Santa Marta. Avances en
recursos hidráulicos, 10(1), 47–62. [online]
http://www.bdigital.unal.edu.co/6096/1/No._10-2003-4.pdf.
Múnera, J. C., Vélez, J. I., Poveda, G., Posada, J. E., Montoya, J. D., and Cardona, Y. M.
(2003b) Dinámica hidrológica de la Ciénaga Grande de Santa Marta., Bogotá D.C.
Neitsch, S. ., Arnold, J. ., Kiniry, J. ., and Williams, J. . (2011) Soil & Water Assessment
Tool Theoretical Documentation Version 2009. Texas Water Resources Institute, TR-
406, 1–647.
Olaya, V. (2011) Sistemas de Información Geográfica, [online]
http://wiki.osgeo.org/wiki/Libro_SIG.
192 Estudio a escala regional de los componentes hidrológicos e hidráulicos
del complejo cenagoso del Bajo Magdalena ente Calamar y Bocas de
Ceniza con apoyo de modelación matemática y SIG
U.S. Water Resources Council (1982) Guidelines for Determining Flood Flow Frequency.
Bulletin 17B: Reston, Virginia, Hydrology Subcommittee, Office of Water Data
Coordination, U.S. Geological Survey, 182.
UNAL-MED (2011) HidroSIG 4.0,
Vargas Cuervo, G. (2015) Guía y Catálogo de Unidades Geomorfológicas en Colombia
por Sensores Remotos, Bogotá D.C, Universidad Nacional de Colombia Sede
Bogotá.
Warren Viessman, J. and Gary L., L. (2008) Introduccion to Hydrology (Prince-Hall, ed.),
Werner, M. G. F. (2001) Impact of grid size in GIS based flood extent mapping using a 1D
flow model. Physics and Chemistry of the Earth, Part B: Hydrology, Oceans and
Atmosphere, 26(7–8), 517–522. [online]
http://www.sciencedirect.com/science/article/pii/S1464190901000430 (Accessed
June 2, 2015).
Zhang, S., Guo, Y., and Wang, Z. (2015) Correlation between flood frequency and
geomorphologic complexity of rivers network – A case study of Hangzhou China.
Journal of Hydrology, 527, 113–118. [online]
http://www.sciencedirect.com/science/article/pii/S0022169415003194 (Accessed
May 27, 2015).