JULIO BELTRÁN-V.
J. ORLANDO RANGEL-CH.

Universidad Distrital Francisco José de Caldas. Bogotá D.C., Colombia.
Programa doctorado en Ciencias-Biología, Facultad de Ciencias, Instituto de Ciencias
Naturales, Universidad Nacional de Colombia, Apartado 7495, Bogotá D.C., Colombia. jebeltranva@unal.edu.co; jebelatran@udistrital.edu.co

Instituto de Ciencias Naturales, Universidad Nacional de Colombia, Apartado 7495, Bogotá D.C. Colombia. jorangelc@unal.edu.co; jorangelc@gmail.com

RESUMEN

Se propone un modelo de simulación dinámico para explicar los procesos hidrológicos del humedal Jaboque – Bogotá D.C., Colombia, en tres secciones con características físicas y químicas particulares. El modelo incluye los caudales de entrada y salida  de agua para cada una de las secciones del humedal. La modelación dinámica se basó en ecuaciones diferenciales. Se utilizó el método Euler de integración con el programa computacional Stella 9.1^®. El comportamiento hidrológico del humedal tiene una alta dependencia de los aportes de la cuenca aferente que a su vez obedecen fundamentalmente a las condiciones climáticas de la cuenca. En el año 1992 (año Niño) la precipitación fue 432.7 mm/año con una descarga de agua estimada  de la cuenca aferente de 905 225 m³/año, siendo el valor más bajo en los últimos 39 años. En el 2010 (año Niña) la precipitación anual fue de 1251.3. mm, con una descarga de la cuenca aferente de 2 471 634 m³/año, que corresponde al dato más alto para el mismo periodo. La comparación entre los caudales de entrada y salida del humedal Jaboque muestran una diferencia bruta de 82 165 m³/mes para un total anual de 985 980 m³/año. El tiempo hidráulico de retención modelado (THR) se estimó en quince para la primera sección, 36 y 42 días para la segunda y tercera sección. El promedio en el humedal fue de 93 días/año. El análisis de sensibilidad mostró que entre los limites de 800 mm/año y 900 mm/año, valores entre los cuales oscila el promedio multianual de la precipitación, el comportamiento del THR es de 13 a 15 días/año en la primera sección, en la segunda de 32 a 36 días/año y la tercera de 38 a 43 días. Se corrobora así la relación que existe entre la precipitación y los factores morfométricos del humedal.

Palabras clave: Humedales, hidroperiodo, modelación, tiempo hidráulico de retención, Jaboque.

ABSTRAC

We propose a dynamic simulation model to explain the hydrological processes in the wetland Jaboque - Bogota, Colombia, in three sections with particular physical and chemical characteristics. The model includes the inflow and outflow for each of the sections of the wetland. The dynamic modeling was based on differential equations. Euler method was used for integration with Stella 9.1 ® computer program. The hydrological behavior of the wetland has a high dependence on afferent input from the basin, which in turn is mainly due to the climatic conditions of the basin. In 1992 (El Niño year), rainfall was 432.7 mm / year with an estimated water discharge of the afferent basin of 905,225 m³/year, which is the lowest value in the last 39 years. In 2010 (La Niña year) the annual rainfall was 1251.3 mm, with a discharge of the afferent basin of 2,471,634 m³/year, that is the highest value ​​for the same period. The comparison between input and output streams of wetland Jaboque show a gross difference of 82 165 m³ per month, that reach for the whole year a total of 985,980 m³. The hydraulic retention time (THR) was estimated in 15 days for the first section, 36 days for the second and 42 days for the third section. The average retention time in the wetland was 93 days. Sensitivity analysis showed that between rain values of 800 mm/year and 900 mm/year, range that include the multiyear average precipitation, the behavior of the THR is 13 to 15 days in the first section, 32 to 36 days in the second section and the 38 to 43 days in the third section, These results corroborate the relationship between precipitation and wetland morphometric factors.

Key Words: Hydroperiod, hydraulic retention time, modeling, wetlands, Jaboque.

Recibido: 15/04/2012
Aceptado: 19/04/2013

INTRODUCCIÓN

Las condiciones hidrológicas pueden modificar o cambiar directamente las propiedades físicas y químicas del agua como la temperatura, el pH, la disponibilidad de nutrientes y los tiempos hidráulicos de retención (Wetzel 2001, Hakanson 2005, Jørgensen 2008) . El estudio de los procesos hidrológicos en humedales naturales y artificiales ha permitido desarrollar propuestas exitosas para la restauración y para la creación de este tipo de sistemas (Mitsch & Day 2004, Zhang & Mitsch 2005) . La modelación y la simulación de los procesos hidrológicos ayuda en la comprensión de los procesos bióticos, los ciclos geoquímicos, los ciclos de nutrientes, la dinámica de la vegetación y el comportamiento trófico (Bradley & Gilvear 2000, Hupp 2000, Quinn & Hanna 2003, Spieles & Mitsch 2003) . El hidroperiodo afecta el nivel de agua en un humedal y es considerado el factor dominante en el control de la composición  florística de la comunidades de plantas acuáticas (Duever et al. 1988) . Cuando las condiciones hidrológicas en un humedal cambian ligeramente, la biota puede responder con cambios  en la riqueza de especies y en la productividad del ecosistema (Mitsch et al. 2009) . Fundamentalmente el hidroperiodo está conformado por la intensidad de la inundación (altura del nivel del agua en el humedal), el tiempo de duración (tiempo que permanece cubierto de agua), la frecuencia (número de veces que el sistema se inunda) y la estacionalidad (época del año cuando se produce la inundación). El hidroperiodo, entonces, define el patrón estacional del nivel del agua en un humedal. 

METODOLOGÍA

Generalidades del área de estudio

El humedal de Jaboque se localiza en la localidad de Engativá al occidente de la ciudad de Bogotá D.C, Colombia (Figura 1). Pertenece a la cuenca del río Juan Amarillo, entre el Aeropuerto Internacional El Dorado y la Autopista Medellín. Limita por el occidente con el río Bogotá, por el sur con los barrios Engativá, Las Mercedes, Puerto Amor, Bolivia, Villa del Mar y la carretera que une a Engativá con el parque La Florida. Por el Oriente limita con los barrios Álamos Norte, Álamos Sur y Bosques de Mariana entre otros. Por el Norte con los barrios Villas de Granada, Los Ángeles y áreas destinadas al pastoreo y cultivo. Abarca aproximadamente 57 hectáreas y presenta una forma alargada en dirección Sur-Oriente-Noroccidente. En las últimas cuatro décadas el humedal del Jaboque ha visto reducida su extensión por la variación del drenaje, la colmatación por cargas orgánicas de las aguas y vertimientos, ocupación por construcción y fraccionamiento del mismo. Esta situación ha incidido para que el humedal haya  disminuido su oferta de servicios ambientales y se hayan alterado atributos relacionados con su biota, su composición y su estructura (Rangel, 2005)

Las simulaciones se realizaron usando Stella 9.1®, software que se ha utilizado ampliamente para realizar modelos hidrológicos y ecológicos (Duever et al. 1988, Mancera et al. 2003, Zhang & Mitsch 2005, Hakanson 2005 entre otros) . Los datos morfométricos fueron tomados por investigadores de la Universidad Nacional de Colombia - la Empresa de Acueducto y Alcantarillado de Bogotá, EAAB, y  del Instituto de Hidrología, Meteorología y Estudios Ambientales, IDEAM, en el marco del proyecto  Investigación aplicada en restauración ecológica en el humedal de Jaboque , Rangel (2005). Se utilizó el sistema de interfaz de Stella para simular cambios mensuales y condiciones extremas de mayor inundación y menor precipitación. De igual manera, para realizar los análisis de sensibilidad del modelo se relacionó la precipitación con las variables dependientes como volumen, tiempo hidráulico de retención (THR) y profundidad. Para verificar su carácter determinístico y realizar la validación del modelo se realizaron regresiones simples lineares (Legendre & Legendre 1998) que permitieron comparar los datos medidos con los resultados modelados. Para evaluar la fortaleza del modelo, se determinó el error relativo (Bryhn & Hakanson, 2007) y para comprobar su carácter predictivo se realizó el análisis de sensibilidad (Jørgensen & Bendoricchio, 2001; Håkanson, 2004) .

A partir del análisis de balance de masas se estimaron los caudales  de entrada y salida del humedal (Tabla 1 y Tabla 2) y se determinó la estructura Anexo 1 de un modelo dinámico de simulación para estimar el comportamiento hidrológico mensual del humedal (Mitsch & Gosselink 2000, Zhang & Mitsch 2005, Acreman et al. 2007) Las ecuaciones se ajustaron  de la siguiente manera:

dv /dt = ∫ent (t) - ∫sal (t)                     (1)

∫ent (t) = P (t) + Si (t)                          (2)

∫sal = So (t) + Et (t)                             (3)

Donde:

dv/dt      = Tasa de cambio del volumen agua, si dv/dt =0

∫ent (t)    = Entradas

∫sal (t)    = Salidas

Si (t)       = Caudal mensual promedio

P (t)        = Precipitación media mensual

So (t)      = Caudal promedio salida

Et (t)       = Evapotranspiración mensual

El caudal de entrada al humedal se calculó a partir de los aportes de la cuenca aferente (m³/mes) que se recogen por medio de los canales de los Ángeles y Carmelo, y que transportan el caudal de agua hacia el interior del humedal. Para los caudales de entrada de cada una de las secciones se tuvo en cuenta el promedio mensual de precipitación efectiva (p.e.-mm) sobre cada sección, más el caudal proveniente de cada  sección.

Precipitación (P) mensual: Los datos de precipitación media mensual se tomaron de la estación hidrometeorológica Aeropuerto El Dorado (IDEAM) entre el período de enero 1972 a diciembre 2010. Se utilizaron para elaborar la caracterización climática de la zona del humedal y para alimentar el modelo de simulación.

Para la predicción de las pérdidas de agua por evapotranspiración para cada sección del humedal se utilizo la relación planteada por Kadlec (1990, 1994) y Werner & Kadlec (1996), la cual se basa en la determinación de las tendencias de evapotranspiración para sistemas de tratamiento en humedales. La evapotranspiración real (ET) es el volumen de agua que realmente sale del sistema en una unidad de tiempo y depende de que exista suficiente agua disponible para evaporar.

Et =        Ket^*ETF^* SA                                     (4)

K_e                 = Constante de evapotranspiración mm/mes

ETF        = Factor de corrección

SA          = Superficie de  la sección m²

Se asumen los siguientes supuestos con el fin de delimitar el alcance del modelo y su posterior generalización o aplicación:

En  cada una de las secciones del humedal se asume su profundidad media.

El caudal de agua al humedal proviene básicamente de la descarga de la cuenca aferente y de la precipitación efectiva en cada una de las secciones.

La descarga de la cuenca aferente se calculó con base en la relación existente entre precipitación y área en km² de la cuenca aferente (Hakanson & Peters,1995).

No se tienen en cuenta los planos de desborde del humedal, debido a la construcción de canales perimetrales en los dos primeros tercios y las obras de infraestructura o terraplenes en el tercio bajo, calle sesenta y cuatro, costado sur y terraplén con ciclo vía, jarillón río Bogotá, barreras artificiales que impiden la inundación y generan un sistema hidrológico de volumen controlado.

Todas las simulaciones se realizaron con los promedios mensuales y anuales de los datos climatológicos entre los años 1972 – 2010.

Dadas las condiciones arcillosas en el sustrato de la zona (Van der Hammen, 1995)  Muñoz & Jaramillo 2005) , se asume que la infiltración es mínima y que no existen caudales de aguas subterráneas hacia el cuerpo de agua del humedal.

Se considera el tiempo de residencia hidráulica (THR) como uno de los factores críticos para explicar la retención de nutrientes en los humedales (Knight 2000, Werner & Kadlec 2000) . El THR se define como el tiempo de retención de un volumen de agua en un cuerpo de agua antes de salir del sistema (Werner & Kadlec 2000) . En la construcción de humedales para tratamiento de aguas residuales el THR es una de las variables principales utilizadas para determinar el diseño de un humedal artificial. El THR puede calcularse de acuerdo con los lineamientos de Almendinger (1999), Håkanson (1999) y Werner & Kadlec (1996).

THR = V = (L^*W^*D)/(Q^*T^-1)

Donde

V     = Volumen de la sección  

T     = Duración del caudal

L     = Longitud del humedal

W    = Ancho promedio

D     = Profundidad promedio

Q     =  Caudal que ingresa al sección

A partir de esta ecuación se simplifica el procedimiento y de acuerdo con Håkanson & Boulion (2002), Håkanson (2004) y Zhang & Mitsch (2005) , el THR se puede calcular  a partir de la siguiente ecuación:

THR = V/Q                                   (6)

Donde THR es el tiempo hidráulico de retención (días o años), V es el volumen de agua (m³), Q es el caudal promedio anual (m³/año), en términos de entradas y salidas del sistema.

Para la predicción de los aportes de agua provenientes de la cuenca aferente del humedal se ajustó el submodelo de descarga (Hakanson & Peters 1995). Se calcularon las descargas medias mensuales (Q) que llegan al humedal Jaboque a través de los canales Carmelo y Los Ángeles. Para esto se tuvieron en cuenta las siguientes variables: área de la cuenca aferente (5.8 km²) y escorrentía (0.01 m³/km²),calculados por los estudios técnicos para la restauración hidrodinámica y ecológica del humedal Jaboque (Universidad Javeriana 2007) , y la precipitación media anual (P) establecida a partir de los datos de precipitación de la Estación Aeropuerto El Dorado 1979 – 2010, del IDEAM.

Q = ADA^* (P/817.7) ^*0.01^*60^*60^*24^*365                                            (7)

Q = 5.8 km²  (P/817.7) ^*0.01^*60^*60^*24^*365                            (8)

Donde:

Q             = Caudal cuenca aferente m³/mes (Tabla 3)

ADA       = Área de la cuenca aferente

P             = Precipitación media anual

Descarga de la cuenca aferente

El comportamiento de la descarga de la cuenca aferente del humedal (Q en m³/mes), se asimila al patrón pluviométrico de la zona (Figura 2 y Tabla 3), donde se muestra la relación existente entre la magnitud de las descargas en m³ y la precipitación en mm, entre los años 1990 y 2010. Se evidencia que las variaciones de la descarga de la cuenca aferente responden al comportamiento de la precipitación; valores extremos de precipitación corresponden también a montos extremos en los caudales. En los últimos 39 años se han presentado variaciones significativas en el comportamiento de la precipitación que obviamente afectan la descarga de la cuenca aferente al humedal. Los años con mayores precipitaciones fueron 1979, 1984, 1990, 1995, 1999, 2004, 2006, 2008 y 2010. En estos años la precipitación promedio anual varió entre 920 mm y 1250 mm, con descargas superiores a 1 933 408 m³/año, que llegaron hasta 2 471 654 m³/año. Los años con menores  precipitaciones fueron 1977, 1978, 1980, 1985, 1989, 1996, 2001 durante los cuales la precipitación promedio anual osciló entre 598.1  mm y 683.3 mm, con descargas que fluctuaron entre 1206 273 m³/año y 1 432 008  m³/año. El promedio anual de la precipitación  para este periodo fue de 828.86 mm y el promedio de la descarga estimada de la cuenca aferente fue 1 718 646 m³/año, lo cual evidencia un alta variabilidad entre los valores extremos de la precipitación y la descarga de agua al humedal de Jaboque.

Se puede inferir una alta influencia de los valores extremos. El menor valor de precipitación anual se presentó en 1992 y fue de 432.7 mm/año. La descarga calculada de la cuenca aferente fue de 905 225 m³, que es el valor más bajo en los últimos 39 años. La precipitación más alta (1250.5 mm,) ocurrió en 2010, con una descarga de la cuenca aferente de 2 471 634 m³. Estos valores afectan de manera significativa la dinámica hidrológica del humedal, ya que se modifica el tiempo hidráulico de retención (Carpenter et al. 2008) , lo que incide en los procesos geoquímicos y biológicos del humedal. Una situación parecida se presenta actualmente en algunos lagos europeos, que se asocia con fenómenos de cambio climático (Håkanson et al. 2003) . En la Figura 3 y en la Tabla 3 se muestran los resultados de la modelación de la precipitación y la descarga de la cuenca aferente para el periodo de 1992 correspondiente al año de máxima sequía, típico representante del fenómeno de El Niño. Se puede apreciar que la curva de la precipitación coincide con la descarga a la cuenca aferente. Predomina un largo periodo seco desde enero hasta la mitad de octubre, cuando se incrementa súbitamente la precipitación. El periodo lluvioso permanece hasta noviembre, lo cual obviamente se refleja en el comportamiento de la descarga de la cuenca aferente.

En la Figura 4 se muestran los resultados de la modelación de precipitación en mm y la descarga de la cuenca aferente para el año 2010, correspondiente al año de máxima precipitación La Niña. Se puede apreciar que la marcha mensual de la precipitación coincide con la de la descarga a la cuenca aferente en diciembre – enero, febrero – marzo y agosto – septiembre. Las mayores descargas ocurren en abril, mayo, junio, octubre y noviembre. Al comparar las curvas de descarga al humedal en las dos situaciones extremas se puede apreciar que durante el fenómeno de La Niña fuerte, la descarga se aumentó aproximadamente en un 25% con relación al promedio histórico y de 150% con relación al comportamiento durante un fenómeno de El Niño. Es decir, que el hidroperiodo se caracteriza por una alta variabilidad que depende de las condiciones climatológicas globales. Esta situación confirma que los fenómenos climáticos globales afectan fenómenos físicos locales de suma importancia, como es el comportamiento del hidroperiodo en el humedal, lo que afecta los procesos biogeoquímicos y biológicos (Hakanson et al. 2003).

La simulación tuvo por objeto conocer y predecir el comportamiento de las fluctuaciones del volumen (m³) de agua del humedal en sus tres secciones, con relación a los caudales de entrada, y en especial al proveniente de los canales Ángeles y Carmelo. Estos recogen las aguas provenientes de sus respectivas cuencas aferentes en la sección uno (Figura 5). También se incluye la entrada de la precipitación efectiva (mm) para cada una de sus secciones. Todo esto se estableció para el periodo comprendido entre los años 2004 y 2005 (24 meses), época en que se realizaron los estudios de investigación para la restauración ecológica del humedal. Se determinó el THR para cada una de las secciones, así como para la totalidad del humedal Jaboque.

El volumen del tercio alto depende fundamentalmente de la descarga de la cuenca aferente, debido en parte a que las estructuras que recogen el caudal proveniente de los canales Carmelo y Los Ángeles se encuentran "encausadas" por estructuras perimetrales que dirigen los excesos de caudal hacia la parte media del humedal. En la Figura 5 se muestran los resultados de la simulación de la descarga de la cuenca aferente y del comportamiento del volumen de la sección para los años 2004 y 2005. En el año 2004 los meses con mayor almacenamiento de agua fueron abril (66 137 m³), mayo (70 643 m³), octubre (72 447 m³) y noviembre (67 643 m³). Estos corresponden a los meses con mayores precipitaciones. Los meses con menor almacenamiento de agua fueron enero (50 300 m³), febrero (53 976 m³) y diciembre (61 500 m³), meses con menores valores en precipitación. En el año 2005 los mayores almacenamientos de agua ocurrieron en abril (66 400 m³), mayo ), (74 643 m³), octubre (71 888 m³) y noviembre (70 500 m³), que fueron los meses con mayores lluvias. Los menores almacenamientos de agua se dieron en enero (58 300 m³), febrero (62 551 m³) y agosto (59 200 m³), meses con bajos valores en precipitación. Se compararon los datos medidos y simulados mediante un análisis de regresión entre los dos grupos. Se encontró una buena correspondencia, r² = 0.704 (y = 0.8692x+23033). La diferencia entre los valores estimados y observados podría deberse en parte a las alteraciones morfométricas que ha sufrido esta sección y en especial a la construcción de los canales perimetrales y de las materas centrales, que originaron una morfometría compleja en esta sección. Las diferencias entre el comportamiento del caudal y el volumen de esta sección se deben en parte al THR, que es de quince días en promedio.

Esta sección se diferencia el tercio alto porque se encuentra medianamente intervenida por obras civiles y recibe agua por escurrimiento y a través de pequeñas aberturas localizadas en los camellones que delimitan la sección. En el año 2004 los meses con mayor almacenamiento de agua fueron abril con 204 071 m³, mayo con 197 066 m³, octubre con 210 747 m³ y noviembre con 201 478 m³; estos fueron los meses con mayores precipitaciones.  El menor almacenamiento de agua se presentó en enero (161 556 m³), febrero (198 801 m³), agosto (176 705 m³) y diciembre (187 598 m³), meses con menores lluvias. En el año 2005 los mayores almacenamientos de agua se dieron en los meses más lluviosos, es decir en abril 212 779 m³), mayo (223 328 m³), septiembre (191 882 m³), octubre (204 788 m³) y noviembre (203 130 m³). Los meses con menor precipitación fueron también los de menos almacenamiento de agua (enero con 182 411 m³, febrero con 189 295 m³, agosto con 188 581 m³ y diciembre con 190 496 m³). En la descarga proveniente del tercio alto, las variaciones del volumen fuero), n perceptibles y se presentaron algunas inundaciones en los barrios de marginales (ADESA 2004) . Sin embargo, su intensidad se vio reducida en el flanco sur debido a la construcción de un canal interceptor que recoge las aguas de escorrentía provenientes del costado Sur del humedal. Al comparar los volúmenes estimados y los simulados se encontró una correlación aceptable, con un nivel de explicación de r² = 0.7181 (y=0.3857x+84833). Las variaciones entre el volumen calculado y el volumen simulado en esta sección pueden obedecer en parte a las discontinuidades batimétricas debidas a la construcción de jarillones  posiblemente de origen precolombino (Van der Hammen 2003, Muñoz 2004). Este patrón se replica en el tercio bajo, donde hay múltiples puntos de acumulación y transporte de sedimentos.

Esta sección, a diferencia de las dos primeras, es la de mayor  extensión del sistema con 315 308.4 m² y un volumen de agua de 230 175 m³ (Izquierdo 2005). Muestra buena capacidad para almacenar agua y amortiguar crecientes en épocas de altas descargas. En el año 2004 los meses con mayor almacenamiento de agua fueron los más lluviosos (abril con 245 155 m³, mayo con 236 697 m³, junio con 222 476 m³, octubre con 253 216 m³ y noviembre con 242 024 m³). Los de menor almacenamiento de agua fueron a su vez los más secos (enero con 193 820 m³, febrero con 202 537 m³ y agosto con 212 111 m³). En el año 2005 los mayores almacenamientos de agua ocurrieron en abril (255 669 m³), mayo (265 992 m³), octubre (246 020 m³) y noviembre (239 621 m³), que corresponden a los meses con mayores precipitaciones. El menor almacenamiento se dio en los meses menos lluviosos (enero con 219 001 m³, febrero con 227 314 m³, agosto con 226 451 m³ y diciembre con 223 067 m³). Al igual que en las secciones anteriores, se encontró una buena correspondencia entre los valores calculados y simulados, con una explicación de r² = 0.7125 (y = 0.4627x+84609). Persiste una diferencia del 29% entre los valores calculados y simulados, que obedece a las complejas condiciones morfométricas (batimetría irregular coadyuvada por la construcción de diques y jarillones especialmente).

La fortaleza de predicción del modelo se verificó mediante la determinación del error relativo cercano a cero de 0.23 para la primera sección y de -0.17 para la segunda y tercera. El resultado permite considerar de manera positiva la utilización del modelo en estudios de caso similares, como son los demás humedales localizados al menos en la jurisdicción de Bogotá D.C.

El volumen total de agua del humedal se descarga hacia el río Bogotá a través de tres salidas o desfogues, dos localizadas en la parte Occidental del tercio alto, en el  jarillón que separa al humedal del río Bogotá y otra, construida en el margen Norte de esta misma sección y que hace parte del jarillón – ciclo vía, que aísla el humedal de su antigua área natural de desborde. Una comparación entre los caudales de entrada y salida  del humedal Jaboque (Figura 6) muestran una diferencia bruta de 82 165 m³/mes para un total anual de 985 980 m³/año, con lo cual el promedio del THR es de 93 días, es decir que el humedal, en las condiciones actuales recambia su volumen cuatro veces al año. Las diferencias entre el patrón de comportamiento del volumen con relación a los caudales de entrada se explican por el THR de la sección.  En la Figura 6 se muestran los resultados de la modelación del comportamiento de la descarga de la cuenca aferente, la descarga al río Bogotá y del volumen de la tercera sección. Nótese que existe coincidencia entre el patrón de picos y valles de cada curva, tanto para los meses de mayores precipitaciones como para los más secos.

Cada sección tiene un comportamiento del THR diferente, debido principalmente a las características morfométricas: volumen, área, profundidad de la cubeta, patrón irregular del fondo y la presencia de diques artificiales de origen precolombino y colonial. El THR mensual promedio del primer tercio corresponde a quince días,  para la sección intermedia el promedio es de 36 días y en el tercer tercio es de 42 días (Figura 7). Estas cifras superan los valores mínimos recomendados para épocas de verano según Zhang & Mitsch (2005) , establecidos para el tratamiento de aguas residuales en lagunas facultativas.

En la Figura 8 y en la Tabla 4 se muestran los resultados de la modelación del comportamiento de la precipitación y del THR  durante el periodo El Niño de 1992. El THR correspondiente a una precipitación mínima de 8.8 mm fue de 108 días año para la primera sección y de 265 días año  y 314 días año para la segunda y tercera sección. Esto indica que a menores precipitaciones el THR aumenta, mientras que con la precipitación de 108.8 mm en diciembre, el THR fue de 9 días en la primera sección, de 22 en la segunda y de 27 en la tercera. En consecuencia, el THR disminuye con las mayores precipitaciones. En la Figura 9 y la Tabla 5 se muestran los resultados de la modelación del comportamiento de la precipitación y del THR en las tres secciones para el año 2010, en el que se presentó el fenómeno de La Niña. En enero y marzo las precipitaciones fluctuaron entre 6.3 mm y 23.2 mm. En estos meses los THR para la primera sección se estimaron entre 41 y 152 días, para la segunda sección entre 102 y 365 días y para la tercera entre 121 y 420 días. Cuando la precipitación se incrementó a valores entre 51.8 mm y 76.4 mm en agosto y septiembre, los THR fluctuaron entre 13 y 20 y días en la primera sección, 32 y 48 días en la segunda y 22 y 38 días en la tercera. En los meses de mayores  precipitaciones (abril, mayo, junio y julio) las precipitaciones fluctuaron entre 187.3 mm y 136.7 mm y los THR para la primera sección se estimaron entre 5 y 8 días, para la segunda sección entre 13 y 24 días/año y para la tercera entre 16 y 22 días. En octubre y noviembre las precipitaciones alcanzaron los 132.9 mm y 136.6 mm respectivamente; en estas condiciones retención hidráulica mostró los valores más bajos. La sección uno tuvo un THR que fluctúo entre 8 y 5 días, en la segunda varió entre 13 y 19 días y en la tercera osciló entre 16 y 22 días. Los resultados evidencian que en época húmeda disminuye notablemente el THR, lo que indica un mayor transporte de agua en las cubetas del humedal y esto posiblemente incrementa el lavado de nutrientes, como lo sugiere Jørgensen (2003).

Para el análisis de sensibilidad del modelo se tuvo en cuenta como variable independiente el promedio anual de precipitación en la zona (P) en la cuenca aferente, y como variables dependientes la descarga al río Bogotá (Tabla 5) y el THR para las tres secciones. Para ambos casos, las variaciones de la precipitación se establecieron 10 clases entre 400 mm/año y 1300 mm/año, a partir del comportamiento de los valores extremos de la precipitación en la cuenca aferente, entre los años de 1990 y 2010 (Tabla 3).  Se puede observar (Figura 10 y Tabla 7) que en la medida en que aumenta la precipitación aumenta la descarga hacia el rio Bogotá, debido a que el tercio bajo recibe el mayor impacto de la descarga de la cuenca aferente. Se corrobora así el supuesto de que en condiciones extremas de precipitación el humedal se convierte en un canal de paso de aguas lluvias, ya que el THR disminuye a medida que aumenta la velocidad de transporte de agua en la sección. Cuando aumenta la descarga de la cuenca aferente los tiempos de retención disminuyen en cada sección (Figuras 11, 12 y 13, Tabla 7), ya que se aumenta el volumen de agua que se transporta entre cada sección del humedal. La variabilidad entre valores mínimos de precipitación (por ejemplo los relacionados con el fenómeno de El Niño) muestra que a promedios de precipitación entre 400 mm/año y 600 mm/año, el THR oscila entre 20 y 30 días para la sección uno, 48 y 72 días  para la sección dos y 57 y 84 días para la sección tres (Tabla 7). En condiciones de alta precipitación, correspondientes al fenómenos de La Niña (entre 1.200 y 1.300 mm/año), el THR se reduce de manera significativa. Para la primera sección los valores son de 9 y 10 días, para la segunda de 22 y 24días y para la tercera sección de 27 y 29 días. Estos resultados concuerdan con las modelaciones del THR realizados para un año Niño 1992 y un año Niña 2010, con lo cual se comprueba el carácter predictivo del modelo.

DISCUSIÓN

Los resultados de la modelación hidrológica del humedal, permiten detectar una fuerte relación entre los montos y distribución de la precipitación de la cuenca aferente y los procesos hidrológicos del humedal. Esto se manifiesta en las tres secciones del humedal de manera distinta. Se definió un patrón estacional del nivel del agua para cada sección. El promedio anual de la precipitación para este periodo fue de 828.86 mm y el promedio de la descarga estimada de la cuenca aferente fue 1 718 646 m³, lo cual evidencia un alta variabilidad entre los valores extremos de la precipitación y la descarga de  agua al humedal de Jaboque. La dinámica hidrológica modifica el tiempo hidráulico de retención (Carpenter et al. 2008) , lo que incide en los procesos geoquímicos y biológicos del humedal. La comparación entre los caudales de entrada y salida  del humedal Jaboque muestran una diferencia bruta de 82 165 m³/mes para un total anual de 985980 m³/año. La modelación entre las entradas y salidas de agua al sistema reproducen de manera parcial el aporte de aguas de Jaboque al río Bogotá, pero también reconfirman el papel amortiguador que tienen todas las secciones del humedal, y en especial la segunda y tercera sección, al recibir y atenuar el impacto de las crecientes de la cuenca aferente del cuerpo de agua.

El volumen de agua en la primera sección depende fundamentalmente de la descarga de la cuenca aferente y la precipitación, en parte por las estructuras que recogen el caudal  proveniente de los canales Carmelo y Los Ángeles, la segunda sección recibe los aportes del tercio alto y la precipitación, y la tercera sección la descarga del tercio medio y amortigua las crecientes fuertes de las dos primeras secciones. Las condiciones climáticas de la zona se reflejan en el comportamiento del hidroperiodo. El mayor volumen de agua de cada sección se presentó en los meses de abril – mayo y octubre – noviembre, con descensos en los meses de diciembre – enero y julio – agosto.  La comparación de los volúmenes de agua estimados y calculados para todas las secciones indica que es necesario disponer de datos de los parámetros morfométricos con mayor detalle, en especial los que hacen referencia a la batimetría del humedal que es compleja y dinámica, debido a los procesos de colmatación que ha sufrido el humedal. La modelación del THR permitió establecer una aproximación al rol hidrológico y ecológico que juega cada sección en el contexto del humedal. Se infiere que la tercera sección cumple una función reguladora y quizá de estabilizador fisicoquímico y biológico. Se confirma así, que el comportamiento del THR permite amortiguar el impacto de la contaminación difusa ocasionada por múltiples vertimientos ilegales en la cuenca aferente.

Los resultados del análisis de sensibilidad permitieron comparar situaciones de variabilidad climática extrema, un año seco con el fenómeno de El Niño y un año húmedo con el fenómeno de La Niña. Los tercios medio y bajo del humedal desempeñan un rol importante como depuradores en estas condiciones climáticas extremas. En condiciones cercanas al monto promedio de la precipitación las secciones dos y tres (que tienen la mayor superficie) estarían asumiendo la función depuradora del ecosistema; mantener esta condición sería la situación más apropiada para el humedal. Los resultados de la modelación indican que las variables de precipitación, caudal de la descarga de la cuenca aferente, profundidad, área de la cubeta y volumen almacenado, permiten comprender el comportamiento del THR, motivo por el cual estas variables se deberían tener en cuenta para futuras investigaciones relacionadas con medidas de restauración y conservación del humedal.

AGRADECIMIENTOS

Al Grupo de Investigación Biodiversidad y Conservación de la Universidad Nacional de Colombia – Instituto de Ciencias Naturales y a la Empresa de Acueducto y Alcantarillado de Bogotá – EAAB – ESP, por el aporte de la información relacionada con el estudio de caso. A la Universidad Distrital Francisco José de Caldas – Rectoría, a la Facultad del Medio Ambiente y Recursos Naturales por el apoyo económico otorgado para realizar el programa de doctorado. Al Centro de Investigaciones y Desarrollo Científico – CIIDC – y al Grupo de Investigación INDESOS, que apoyaron la financiación para la adquisición de los programas computacionales, el equipo de computo y la  bibliografía especializada que se utilizó  en la investigación.

LITERATURA CITADA

1. Acreman, M., J. Fisher, C. Stratford, D. Mould, D & J. Mountford, J. 2007. Hydrological science and wetland restoration: some case studies from Europe. Hydrology and Earth System Sciences 11(1): 158–169.         [ Links ]

2. Adesa. 2004. Plan de manejo integral y comunitario para la recuperación y protección del humedal El Jaboque ubicado en la localidad de Engativá. Bogotá D.C.         [ Links ]

3. Almendinger, J. 1999. A method to prioritize and monitor wetland restoration for water-quality improvement. Wetlands Ecology and Management 6(4): 241–252.         [ Links ]

4. IDEAM. 2010. Registros climatológicos del Sistema de Información Hidrometereológica. Instituto de Hidrología, Meteorología y Estudios Ambientales de Colombia IDEAM. Bogotá D.C.         [ Links ]

5. Bradley, C. & D. Gilvear, 2000. Saturated and unsaturated flow dynamics in a floodplain wetland. Hydrological Processes 14(16-17): 2945–2958.         [ Links ]

6. Bryhn, A. & L. Håkanson. 2007. A comparison of predictive phosphorus load-concentration models for lakes. Ecosystems 10(7): 1084–1099.         [ Links ]

7. Carpenter, S., N. Caraco, D. Correll, D. Howarth, R. Sharpley & V. Smith. 2008. Nonpoint pollution of surface waters with phosphorus and nitrogen. Ecological Applications 8(3): 559–568.         [ Links ]

8. Duever, M., W. Mitsch, M. Straskraba & S. Jørgensen, S. 1988. Wetland Modelling. Elsevier Scientific Publishing, Amsterdam. 238 pp.         [ Links ]

9. Hakanson, L. & H. Peters. 1995. Predictive limnology: methods for predictive modelling. SPB Academic publishing, Amsterdam. 464 pp.         [ Links ]

10. Håkanson, L. & V. Boulion. 2002. The lake foodweb: modelling predation and abiotic/biotic interactions. Backhuys Publishers, Amsterdam.         [ Links ]

11. Håkanson, L. 1999. On the principles and factors determining the predictive success of ecosystem models, with a focus on lake eutrophication models. Ecological Modelling 121(2-3): 139–160.         [ Links ]

12. Håkanson, L. 2004. Break-through in predictive modelling opens new possibilities for aquatic ecology and management a review. Hydrobiologia 518(1): 135–157.         [ Links ]

13. Hakanson, L. 2005. The importance of lake Morphometry for the structureand function of lakes. International Review of Hydrobiology 90(4): 433–461.         [ Links ]

14. Håkanson, L., A. Ostapenia, A. Parparov, K. Hambright & V. Boulion. 2003. Management criteria for lake ecosystems applied to case studies of changes in nutrient loading and climate change. Lakes & Reservoirs: Research & Management 8(2): 141–155.         [ Links ]

15. Hupp, C. 2000. Hydrology, geomorphology and vegetation of Coastal Plain rivers in the south-eastern USA. Hydrological Processes 14(16-17): 2991–3010.         [ Links ]

16. Izquierdo, D. 2005. Modelación hidráulica del humedal Jaboque. Informe final convenio de cooperación Empresa de Acueducto y Alcantarillado de Bogotá D.C. - Universidad Nacional de Colombia - Instituto de Ciencias Naturales. Bogotá D.C.         [ Links ]

17. Jørgensen, S. & G. Bendoricchio. 2001. Fundamentals of ecological modeling. Elsevier Science Ltda, Amsterdam.         [ Links ]

18. Jørgensen, S. 2008. Overview of the model types available for development of ecological models. Ecological Modelling 215(1-3): 3–9.         [ Links ]

19. Jørgensen, S. E. 2003. The aplication of models to find the relevance of residence time in lake and reservoir management. Papers from Bolsena Conference en Papers from Bolsena Conference: Residence time in lakes: science, Management, Education. Limnology 62: 16–20.         [ Links ]

20. Kadlec, R. 1990. Overland flow in wetlands: vegetation resistance. Journal of Hydraulic Engineering 116: 691–707.         [ Links ]

21. Kadlec, R. 1994. Detention and mixing in free water wetlands. Ecological Engineering 3(4): 345–380.         [ Links ]

22. Knight, R. 2000. Constructed wetlands for livestock wastewater management. Ecological Engineering 15(1-2): 41–55.         [ Links ]

23. Legendre, P. & L. Legendre. 1998. Numerical Ecology. Vol. 20 Elsevier Science, Ámsterdam. 870p.p         [ Links ]

24. Mancera, J., E. Peña, R. Giraldo & A. Santos. 2003. Introducción a la modelación ecológica. Principios y aplicaciones. Universidad Nacional de Colombia. San Andrés. 109 p.p        [ Links ]

25. Mitsch, W. & J. Day. 2004. Thinking big with whole-ecosystem studies and ecosystem restoration-a legacy of H.T. Odum. Ecological Modelling 178(1-2): 133–155.         [ Links ]

26. Mitsch, W. & J. Gosselink. 2000. The value of wetlands: importance of scale and landscape setting. Ecological Economics 35(1): 25–33.         [ Links ]

27. Mitsch, W., A. Nahlik, P. Wolski, B. Bernal, L. Zhang. & L. Ramberg. 2009. Tropical wetlands: seasonal hydrologic pulsing, carbon sequestration, and methane emissions. Wetlands ecology and management 18(5): 573–586.         [ Links ]

28. Muñoz, J. M. 2004. Humedal Jaboque, evolución geomorfológica y geológica, y su relación con las culturas prehispánicas. Universidad Nacional de Colombia, Departamento de Geología. Bogotá D.C.         [ Links ]

29. Muñoz, J. M., & J. Jaramillo. 2005. Evaluación geológica – geomorfológica del humedal Jaboque. Plan de restauración ecológica del humedal del Jaboque. Informe final. Instituto de Ciencias Naturales, Universidad Nacional de Colombia. Bogotá D.C.         [ Links ]

30. Quinn, N. & W. Hanna. 2003. A Decision Support System for Adaptive Real-Time Management of Seasonal Wetlands in California. Environmental Modelling and Software 18(6): 503–511.         [ Links ]

31. Rangel,Ch., J. O. (Compilador). 2005. Investigación aplicada en restauración ecológica en el humedal de Jaboque. Convenio Acueducto de Bogotá – Universidad Nacional de Colombia, Informe final. Bogotá D.C.         [ Links ]

32. Sierra, O. R. & C. A., Monsalve. 2005. Zonificación trófica del humedal Jaboque, con base en la familia Bacillariophyceae, Engativá. En: Rangel-Ch., J. O. (Compilador). Investigación aplicada en restauración ecológica en el humedal de Jaboque. Convenio Acueducto de Bogotá – Universidad Nacional de Colombia, Informe final. 1 – 23 p.         [ Links ]

33. Spieles, D. & W. Mitsch. 2003. A model of macroinvertebrate trophic structure and oxygen demand in freshwater wetlands. Ecological Modelling 161(3): 181–192.         [ Links ]

34. Van Der Hammen, T. 1995. Plioceno y Cuaternario del altiplano de Bogotá y alrededores. Análisis Geográficos. Instituto Geográfico Agustín Codazzi IGAC. Bogotá D.C.         [ Links ]

35. Van Der Hammen, T. 2003. Los humedales de la sabana. Origen, evolución, degradación y restauración. En: A. Guarnizo, & B. Calvachi. Los humedales de Bogotá y la Sabana. Empresa de Acueducto y Alcantarillado de Bogotá – Conservación Internacional, Bogotá D.C.         [ Links ]

36. Werner, T. & R. Kadlec. 1996. Application of residence time distributions to stormwater treatment systems. Ecological Engineering 7(3): 213–234.         [ Links ]

37. Werner, T. & R. Kadlec. 2000. Wetland residence time distribution modeling. Ecological Engineering 15(1-2): 77–90.         [ Links ]

38. Wetzel, R. 2001. Fundamental processes within natural and constructed wetland ecosystems: short-term versus long-term objectives. Wetland Systems for Water Pollution Control 2000 44(11-12): 1–8.         [ Links ]

39. Zhang, L. & W. Mitsch. 2005. Modelling hydrological processes in created freshwater wetlands: an integrated system approach. Environmental Modelling & Software 20(7): 935–946.         [ Links ]