1

16MODELACIÓN DE CALIDAD DE AGUAS

PROYECTO DE CLASE
Karen Lorena Bernal Rincón 5400316
1. INFORMACIÓN:

1.1 DESCRIPCIÓN GENERAL:

El tramo de estudio seleccionado corresponde a la cuenca del río alto Ubaté, perteneciente a la
cuenca hidrográfica de los ríos Ubaté y Suárez, con una longitud de 7,74 Km; Dentro de esta
cuenca se encuentran localizados los municipios de Carmen de Carupa y Ubaté.
Cabe mencionar, que la cuenca del río alto Ubaté pertenece a la jurisdicción de la Corporación
Autónoma Regional de Cundinamarca, CAR, Regional Ubaté y sus principales afluentes son los
ríos Hato y Suta (CAR,2010)

2550 m.s.n.s

2588 m.s.n.s

Imagen 1: Tramo de estudio seleccionado

Fuente: Google Maps

1.2 CAUDALES Y CALIDAD DEL AGUA:

Las condiciones de calidad de agua en el tramo de estudio se ven afectadas por los vertimientos de
industrias lecheras , minería, agricultura, actividades pecuarias y domésticas, provenientes de los
municipios de Carmen de Carupa y Ubaté; Sin embargo, para conocer el diagnóstico de la calidad
y de las características físicas de la cuenca del río alto Ubaté, la Corporación Autónoma Regional
de Cundinamarca- CAR estableció estaciones de monitoreo presentes dentro de la zona de estudio
(CAR,2010), las cuales se mencionan a continuación:
No. Estación
1 LM La Malilla
2 Río Hato
3 Estación La Boyera
4 LG Pte. Barcelona
5 Río Suta.
Tabla 1: Estaciones de monitoreo
 2

Estaciones de monitoreo
Estación LG Pte. Río Suta
Parámetros Unidades LM La Malilla Río Hato La Barcelona
Boyera
Caudal 𝑚3 /𝑠 0,209 0,7025 0,7435 0,7365 0,6615
Temperatura °C 16 15,4 14 15 17
OD mg/L 6,3 6 5,9 5,7 2,0
DBO mg/L 4,5 4 6,4 4,7 25,7
SST mg/L 21,4 11,6 12,3 15,1 98,8
Coliformes NMP/100 4160 13898,5 7841,7 29146,2 592214,5
totales ml
Nitratos mg/L 0,4 0,5 1 0,5 0,5
Nitritos mg/L 0,004 0 0 0 0,1
P total mg/L 0,1 0,1 0,1 0,1 5,1
N-Amoniacal mg/L 0,7 0,6 0,6 0,6 25,3
pH 6,7 6,7 6,7 6,7 6,7
Tabla 2: Caracterización fisicoquímica del tramo de estudio.
Fuentes: (CAR, 2010); (Montañez,2019))

Los datos anteriores fueron recolectados de la campaña de monitoreo efectuada por la Corporación
Autónoma Regional de Cundinamarca- CAR en el año 2009 y de un trabajo que se realizó en la
Pontificia Universidad Javeriana, acerca del análisis fisicoquímico de los ríos que desembocan la
laguna de Fúquene. Se considera en términos generales que la calidad de la cuenca del río alto
Ubaté es aceptable, pero presenta un deterioro, después de la desembocadura del río Suta, donde
la calidad se evidencia como mala (Montañez,2019).

1.3 BALANCE HÍDRICO Y CALIDAD DEL AGUA PARA EL TRAMO

SELECCIONADO:
El balance hídrico para la cuenca alta del río Ubaté considera que a partir de las entradas menos
las salidas de la cuenca se obtienen, la variación en el almacenamiento. Es decir, compara los
aportes de agua que entran al sistema mediante la precipitación y vertimientos, con respecto a las
salidas dadas por la evapotranspiración y captación de agua, teniendo en cuenta la ley de
conservación de masa, como se muestra a continuación:

𝑨𝒍𝒎𝒂𝒄𝒆𝒏𝒂𝒎𝒊𝒆𝒏𝒕𝒐 = 𝑬𝒏𝒕𝒓𝒂𝒅𝒂𝒔 − 𝑺𝒂𝒍𝒊𝒅𝒂𝒔

Precipitación Evapotranspiración

𝑸𝐸𝑛𝑡𝑟𝑎𝑑𝑎
Río Alto Ubaté 𝑸𝑆𝑎𝑙𝑖𝑑𝑎

Aportes de Caudal Captación

(Vertimientos, Ríos)
 3

Por lo tanto, se tiene el siguiente balance:

𝑄𝑆𝑎𝑙𝑖𝑑𝑎 = 𝑄𝐸𝑛𝑡𝑟𝑎𝑑𝑎 + 𝑄𝐴𝑝𝑜𝑟𝑡𝑒𝑠 − 𝑄𝐶𝑎𝑝𝑡𝑎𝑐𝑖ó𝑛 + 𝑃𝑟𝑒𝑐𝑖𝑝𝑖𝑡𝑎𝑐𝑖ó𝑛 − 𝐸𝑣𝑎𝑝𝑜𝑟𝑎𝑐𝑖ó𝑛

Entonces, los principales aportes de caudal son el Río Suta, el Río Hato y las aguas tratadas de las
industrias lecheras, los cuales corresponden a la oferta hídrica; Por otra parte, la captación de agua
se da por los municipios de Ubaté y Carmen de Carupa para consumo humano, agrícola, pecuario,
industrial e institucional (CAR, 2010), lo cual corresponde a la demanda hídrica.
En el siguiente mapa topográfico del río alto Ubaté, se puede observar claramente la oferta y la
demanda hídrica.

Imagen 2: Topología del río alto Ubaté.

Elaboración propia

1.4DIMENSIONES DEL CANAL Y PENDIENTE PARA DIFERENTES SECCIONES

DEL RÍO:
Según el Plan de ordenamiento de la cuenca río alto Ubaté, a nivel fisiográfico y morfológico,
describe características generales con respecto a su forma según la topografía del terreno lo que
determina físicamente la estructura de la cuenca (POMCA, 2009). A continuación, se presenta un
consolidado aproximado de las características físicas de la cuenca río alto Ubaté en diferentes
 4

tramos:

Parámetros Unidades Tramo 1 Tramo 2 Tramo 3 Tramo 4 Tramo 5

Ancho Medio1 m 9,77

Pendiente 2 % 9,7 8,2 8,1 5,6 3
Profundidad m 0,18 0,18 0,32 0,36 0,34
Km2 0 0,645 2,885 6,335 7,74
Tabla 3: Dimensiones aproximadas del canal en diferentes puntos.
1
POMCA
2
Google earth pro.

1.5 CARACTERIZACIÓN CLIMÁTICA DE LA ZONA DE ESTUDIO:

Según el POMCA, para el caso de la cuenca río Alto Ubaté se tiene:
-Por clasificación de Caldas es una cuenca de piso térmico Frío.
-Por clasificación de Lang es una cuenca con una relación de precipitación/ temperatura de 48 se
deduce que el clima es semiárido
-Por Caldas –Lang se tienen una clasificación climática de Frío-semiárido Fsa para la cuenca.

Por otra parte, las condiciones de temperatura, humedad relativa, brillo solar y radiación solar en
la cuenca del río Ubaté son homogéneas:

Temperatura (°C) 14
Humedad relativa % 70%
Brillo Solar 157 h/d
Radiación Solar 471,5 𝑐𝑎𝑙/𝑐𝑚2
Velocidad del viento 1,3 m/s
Dirección del viento SE
Evaporación 22 mm-30 mm
Tabla 4: Caracterización climática del tramo.

2. CONFIGURACIÓN DEL MODELO:

2.2 METODOLOGIA:
- Headwater: La cabecera del tramo de estudio corresponde a los datos de la estación la Malilla
(Montañez,2019), donde se encuentra el nacimiento del río Hato, a la altura del Páramo de Salinas;
Cabe resaltar que el valor de la DBO en esta estación es de 4,5 mg/l, por lo tanto, se asume que el
50% de este valor corresponde a la DBOfast y el otro 50% a la DBOslow.
 5

Por otra parte, se asumieron los datos de conductividad y alcalinidad de la cuenca río alto Suarez
[5]. Se toma como referencia esta cuenca, ya que, se localiza al norte del tramo de estudio, por lo
tanto, se considera la más cercana.

- Reach: El criterio para la división de los tramos fue el cambio de pendiente, utilizando la
aplicación de google earth pro [7]; Se observó que, una parte del tramo de estudio, presenta una
superficie plana que varía entre 1% a 3 % y otra parte que presenta una pendiente fuertemente
ondulada entre 10% a 12%. Además, se definió los coeficientes y los exponentes para la velocidad
y la profundidad, mediante el método de curvas de potencia, asumiendo datos de aforos realizados
en diferentes ríos de la sabana, en base a algunos datos ya encontrados. ANEXO [1]

Imagen 3. Segmentación del tramo de estudio (cambio de pendiente)

- Air Temperature: La temperatura del área de estudio es casi constante durante todo el año, por
lo tanto, la temperatura promedio es de 14 °C. (CAR,2010)

- Dew Point Temperature: En esta cuenca, se encuentra ubicada una estación pluviométrica
denominada “los novilleros” con el código 2401519; Teniendo en cuenta esto, se descargaron los
datos de punto de rocío (°C) de la página del ideam en el periodo comprendido entre 1978 al 2015
(IDEAM); Por lo tanto, se hizo un promedio con todos los datos, dando como resultado un punto
de rocío de 8,89 °C.

-Wind Speed: En el informe diagnóstico de la cuenca río alto Ubaté elaborado por la CAR, indica
que la velocidad del viento promedio en la estación de novilleros es de 1,3 m/s.

- Cloud Cover: Para representar el porcentaje razonable de nubes que cubre el cielo para el
tramo de estudio; Se tomaron fotos de cielo desde Zipaquirá- Cundinamarca, durante tres días en
horas diferentes y se les realizó una cuadrícula 4x4; Esto permitió identificar la fracción
aproximada del cielo cubierto por nubes sobre la cuenca río alto Ubaté.
Nota: Las fotos se encuentran en el ANEXO 2 -hoja “Fotos”.
 6

- Shade: Se consideró el porcentaje de radiación solar que es bloqueada por la vegetación en cada
segmento del tramo de estudio, mediante las imágenes satelitales de google earth pro.

-Diffuse Sources: En este tramo el cultivo principal es el de la papa; Por lo tanto, con ayuda del
google earth pro, se identificó un área de 0,1 𝐾𝑚2 de cultivo de papa.
Para hallar el caudal, se asume que hay una precipitación constante de 9 mm/h que cae
homogéneamente sobre la cuenca y que el 30% de esa precipitación se infiltra por escorrentía
con un coeficiente de 0,4 [13]; Aplicando la ecuación correspondiente, el caudal de esta fuente
difusa es de0,03𝑚3 /𝑠.
Por otra parte, se tomaron valores de referencia de diferentes fertilizantes utilizados para este
cultivo [9][10][11], con el fin de tener una caracterización fisicoquímica de estas aguas por
escorrentía.

-Point Sources: En este tramo hay dos fuentes puntuales que corresponden a corrientes naturales
tributarias que son el Río Hato y el Río Suta ; Adicionalmente las industrias lecheras (Doña leche,
Parmalat y lácteos Casaleta) vierten sus aguas residuales en este tramo [4]; Sin embargo, para
estos vertimientos solo se tiene el dato de caudal, por lo tanto, se buscaron valores de referencia
para la caracterización fisicoquímica de las aguas residuales de industrias lecheras que
tuvieran como tratamiento trampa de grasas [8], como es el caso de las industrias lecheras de
Ubaté.
Finalmente, en algunas partes del tramo hay extracción de caudal, que se asume como abstracción
de fuentes puntuales.

- Temperature Data: Se tenían los siguientes valores de temperatura media del agua: 16 °C,
15,6°C, 14 °C, 15 °C y 17 °C los cuales se distribuyeron en los diferentes segmentos del tramo,
teniendo en cuenta la ubicación de las estaciones de monitoreo.

- WQ Data: Sobre este tramo hay 5 estaciones de monitoreo (Tabla 1) que se basan en el análisis
de la calidad hídrica del río; Los datos de cada estación, fueron encontrados en un estudio que se
realizó en la laguna de Fúquene, donde se evaluó la calidad fisicoquímica del agua de los ríos que
desembocan en esta laguna (Tabla 2), como es el caso del río Ubaté. (Montañez,2019).

2.2 RESULTADOS Y ANÁLISIS DEL ESCENARIO ACTUAL: (Gráficas- ANEXO 2)

2.2.1 Hidráulica del río:

-Tiempo de Viaje y velocidad:

 7

El tiempo de viaje es un parámetro importante en la aplicación de modelos de calidad del agua o

técnicas predictivas que se estima mediante estudio de trazadores, con el fin de conocer las posibles
consecuencias de un contaminante sobre el río.
En el caso del río alto de Ubaté, el tiempo de viaje varía entre 0, 0 y 4,01 días para una distancia
de 7,74 Km; Sin embargo, del kilómetro 7,74 al 6,34, el tiempo de viaje es relativamente “bajo”,
ya que, para recorrer 1,4 Km se está demorando aproximadamente 3 horas, además en el kilómetro
4,87 al 2,19 se está demorando 18 días y en total el tiempo de viaje en el tramo es de 36,33 días
(Gráfica 1) ; Por lo tanto, se estima que es un tiempo de viaje demasiado alto para una distancia
tan corta y esto se puede deber a la pendiente casi uniforme en la mayor parte del tramo que varía
entre 1% a 3%.
Por otra parte, se evidencia una relación inversamente proporcional entre el tiempo de viaje y la
velocidad, ya que, en el kilómetro 7,74 al 6,3, los tiempos de viaje son bajos y las velocidades son
altas, y en el kilómetro 4,87 al 2,19 los tiempos de viaje son altos y las velocidades son bajas
(Gráfica 2); Es decir, que en la mayor parte del río alto de Ubaté se tiene tiempo de viajes altos,
por ende, velocidades bajas.

- Caudal:
El perfil del caudal en el río, muestra una baja correlación entre el modelo y los datos encontrados
en la bibliografía, esto indica que el modelo no describe de forma adecuada las características
hidráulicas del río alto Ubaté, lo cual puede estar dado por los coeficientes de velocidad y
profundidad que no funcionaron con el modelo y arrojo datos de las dimensiones del canal que no
son consistentes, como el ancho; Por esta razón, se cambiaron los coeficientes por los que estaban
inicialmente en el modelo de QUAL2K, pero tampoco funcionaron.
Sin embargo, en el kilómetro 7,30, hay una extracción de caudal de 0,8 𝑚3 /𝑠 (Gráfica 3), lo cual
se asocia con lo descrito en el mapa topográfico (imagen 2) en el abscisado K7+00; Por otra parte,
el caudal se mantiene constante en la mayoría del tramo, pero en el kilómetro 1, hay un aporte de
caudal importante que está asociado a la entrada del afluente del río Hato y de las aguas residuales
de las industrias lecheras, por lo tanto, la velocidad en este punto es más alta (grafica 2).

-Profundidad
Las profundidades más altas en el Río alto Ubaté, se encuentran en el kilómetro 7,30 al 4,87, que
varía entre 0,43 a 0,57 m y las más bajas en el kilómetro 4,65 al 0,00 con una profundidad que
varía entre 0,29 y 0,32 m (Gráfica 4) ; Estos cambios de profundidad en el río pueden estar dados
por la velocidad y el caudal, por ejemplo, entre kilómetro 7,74 al 4,87, hay un aporte de caudal
importante del río Suta, lo que hace que la velocidad aumente (grafica 2) , por ende, también
aumenta la profundidad; También puede ocurrir que, en un punto se tenga un caudal alto, una
velocidad alta y una profundidad baja, como es el caso del kilómetro 1, puede ser que en este
punto, el ancho del canal haya cambiado; Cabe resaltar, que el río alto Ubaté es poco profundo,
por lo tanto, dominara el crecimiento de plantas fijas y va una haber mayor tasa de sedimentación
(Chin,2013).
 8

Estas relaciones hidráulicas son consecuencia de los datos de caudal, profundidad y velocidad que
se asumieron para definir las características hidráulicas mediante las curvas de potencia.

2.2.2 Calidad del Agua:

-Temperatura: La temperatura muestra una alta correlación entre el modelo y los datos obtenidos
de la bibliografía; En el principio y final del tramo de estudio la temperatura es un poco más alta
(Gráfica 5), debido a los vertimientos, por lo tanto, en estos dos puntos va a haber menos oxígeno
disuelto, lo cual va a afectar la actividad biológica y el crecimiento de las especies de plantas y
animales que se encuentran dentro del rio.

-Re-areación:
La re-areación, se da en los puntos donde hay vertimientos, como ya se mencionó anteriormente;
Estos vertimientos tienen concentraciones de DBO >4 mg/L que van a agotar el oxígeno disuelto
en el río para que las bacterias descompongan esa materia orgánica, cuando esto ocurra el río va a
llegar a una zona crítica donde la descomposición va a ser igual a la re-areación y posteriormente
va a ocurrir una recuperación del oxígeno por el intercambio de gases con la atmosfera y por el
proceso de fotosíntesis; Sin embargo, se debe tener en cuenta la capacidad que tiene el río para re
oxigenarse, la cual está dada por una tasa de re-areación.
En el modelo se aplicó el método de Owens y Gibbs, ya que relaciona un río que tiene
profundidades y velocidades bajas (Bowie .et al,1985) como es el caso del río alto de Ubaté; Por
lo tanto, el tramo inicial tiene una mayor capacidad de re-oxigenación que los demás puntos
(Gráfica 6), por su poca profundidad y por tener una mayor superficie expuesta a la atmosfera.

-Oxígeno Disuelto:
Este parámetro es muy importante, debido a que determina la cantidad de oxígeno disponible para
que se dé la actividad química y biológica en el agua del río. En la gráfica 7, se puede deducir que,
el río alto de Ubaté, contiene la cantidad suficiente de oxígeno disuelto > 5 mg/L para la
supervivencia de las plantas y los animales acuáticos, que puede provenir de las corrientes de aire
y de las tasas de re-aireación; Además, se observa que los vertimientos al inicio y al final
disminuyen considerablemente los valores de oxígeno disuelto por el exceso de materia orgánica;
Sin embargo, no se considera preocupante, ya que, estos valores no están por debajo de 3mg/L,
por ende no se generan condiciones anóxicas (EPA, 2020).
Por otra parte, la saturación de oxigeno depende principalmente, de la presión atmosférica de río
Ubaté que es de 765 mmHg y de la temperatura promedio del agua que esta alrededor de 16 °C,
dando como resultado una saturación de oxigeno alrededor de 9,8 mg/L. Este valor genera cierta
 9

incertidumbre, porque, no es común tener este valor de saturación para la zona en la que se está
trabajando, por lo tanto, este valor se verifico en el software DOTABLE (USGS) como se muestra
a continuación:

Imagen 4: Verificación de la saturación de oxigeno.

Fuente: Software DOTABLES (USGS)
En la anterior imagen, se puede observar que los valores de saturación de oxígeno arrojados por el
modelo, si coinciden con las condiciones del río Alto Ubaté.

-Patógenos:
Los datos de patógenos ingresados, muestran una baja correlación con los datos arrojados por el
modelo; En la gráfica 8, se observa que en todo el tramo del río hay presencia de patógenos, esto
se puede deber, al fácil acceso que tiene el río para el ingreso de personas y animales o por las
corrientes naturales del río Hato y Suta, que pueden traer patógenos por lo vertimientos generados
allí.
-Sólidos suspendidos totales (SST):
Los Sólidos suspendidos totales pueden afectar la penetración de la luz en la superficie del agua y
por tanto interfiere en la vida acuática, en este caso los valores de SST para el río alto de Ubaté
son aceptables, ya que se encuentran por debajo de los 500 mg/l. (Gráfica 9)

3. ESCENARIOS DE TMDL: (Graficas y cálculos -ANEXO 3)

3.1 BREVE DESCRIPCIÓN DE LA METODOLOGIA:
Los escenarios de cargas diarias máximas totales (TMDL), se evaluaron principalmente para las
aguas residuales de la heladería san Jerónimo S.A, las cuales van a ser descargas a 50 metros
después de la cabecera del río alto Ubaté; Para esto, fue necesario obtener información de la
caracterización fisicoquímica del agua residual de esta empresa, que fue encontrada en un informe
realizado por la universidad de la sabana (Sastre, 2016). Los parámetros de interés fueron la
demanda bioquímica de oxígeno, fósforo total, oxígeno disuelto, nitritos y nitratos.
Por lo tanto, se procedió a realizar un balance de masas simple para cada uno de estos parámetros,
punto por punto, teniendo en cuenta la información ya obtenida de la calidad del agua del río alto
Ubaté (Tabla 2), para conocer como iba a ser el comportamiento del río en respuesta a este
vertimiento y se describirá como “el escenario actual”.
Dicho lo anterior, este escenario debe cumplir con los límites máximos permisibles establecidos
en la norma que regula el tema ambiental concerniente a recursos hídricos, establecida por la CAR,
que es el acuerdo 43 del 15 de octubre de 2007; En caso de que esto no fuera así, se realizaban las
reducciones necesarias en el vertimiento de las aguas residuales de la heladería San jerónimo para
 10

determinar la carga diaria maxima (TMDL) permitida por el rio para que se cumpliera con
estándares de calidad y finalmente se aplicaba la siguiente formula:
𝑘𝑔 𝑚3 𝑚𝑔
𝐶𝑎𝑟𝑔𝑎 ( ) = 𝐶𝑜𝑛𝑐𝑒𝑛𝑡𝑟𝑎𝑐𝑖ó𝑛 ( ) ∙ 𝐶𝑎𝑢𝑑𝑎𝑙( )
𝑑 𝑠 𝐿
Ecuación 1: Determinación de la carga contaminante.

3.2 ANALISIS DE RESULTADOS:

Nota: Cs: Concentración de salida del vertimiento de la heladería de San jerónimo S.A.

Según la resolución 3462 de 28 de diciembre del 2009, el río Ubate desde la estación limnimétrica
La Malinilla hasta la confluencia con el río Suta y sus afluentes corresponden a Clase II del acuerdo
43, que hace referencia a valores de los usos del agua para consumo humano y doméstico con
tratamiento convencional, uso agrícola con restricciones y uso pecuario, e industrial. Los
vertimientos de este río, incluido el de la heladería San Jerónimo S.A, son principalmente de
actividades productivas de agroindustrias y cumplen con los parámetros fisicoquímicos y sus
valores máximos permisibles definidos en el articulo 9° de la resolución 631 de 2015.

En base a lo anterior, el efluente de la heladería San Jerónimo S.A, ingresa al río y debe pasar al
otro tramo cumpliendo el objetivo de calidad de agua de los parámetros de interés (DBO, OD, PT,
Nitritos y nitratos) según la clase II del acuerdo 43, teniendo en cuenta las concentraciones y los
caudales que se tiene antes de este vertimiento; Con esto presente, se realizó el balance de masas
simple. En primer lugar, la concentración de DBO que lleva el río antes de la descarga de la
heladería san jerónimo es de 4,58 mg/L con un caudal de 0,209 𝑚3 /𝑠 y luego ingresa el
vertimiento de las aguas residuales de San jerónimo con una concentración de DBO de 425 mg/L
con un caudal de 0,00428 𝑚3 /𝑠 (Sastre, 2016), dando como resultado una concentración final
de DBO de 12,95 mg/L con un caudal de 0,214 𝑚3 /𝑠 , por ende, este vertimiento incumple con
el acuerdo 43, dado que: Cs(12,95 mg/L) > (10,0 mg/L) Acuerdo 43 como se observa en la gráfica
10.

Por lo tanto, se realizó las reducciones necesarias en la concentración de DBO y el caudal del
vertimiento, en este caso fue del 20% y el 30%, respectivamente, en base al escenario actual; De
tal manera, que el vertimiento de la heladería de San jerónimo S.A cumplió con el acuerdo 43, ya
que, se obtuvo que: Cs (9,25 mg/L) < (10,0 mg/L) Acuerdo 43 como se observa en la gráfica 11;
Lo que da lugar a que la carga contamínate diaria máxima soportada de DBO en este tramo (0,486
al 0,5 km) es de 0,088 kg/día, lo que representa una reducción del 56 % en comparación con
escenario actual. Cabe resaltar, que se evidencio que, a mayor caudal, mayor concentración de
DBO, lo cual está estrechamente relacionado con la disminución de oxígeno disuelto; Sin embargo,
la EPA (2009), menciona que las condiciones criticas de OD se dan es cuando hay bajo caudal, ya
que, en teoría habría una velocidad de flujo promedio más lenta, por ende, va a haber menos
aireación y mayor impacto proporcional por la DBO.
 11

No obstante, esto depende principalmente de la ubicación, las condiciones hidráulicas y climáticas

del rio. En este caso, el río alto de Ubaté tiene poca profundidad, una velocidad baja (gráfica 2) y
una mayor superficie expuesta a la atmosfera, lo que hace que allá una mayor aireación (gráfica
6) , por ende, la concentración de salida de OD en el vertimiento de la heladería de san jerónimo
S.A no se encuentra en condiciones críticas, a pesar de que la concentración de entrada de OD del
vertimiento fue de 2,1 mg/L (Sastre, 2016); Dicho esto, este vertimiento está ajustado al
cumplimiento del acuerdo 43, dado que: Cs(6,21 mg/L) > (4 mg/L) Acuerdo 43, como se observa
en la gráfica 12; Cabe adicionar, que incluso con una concentración de entrada del OD más baja
en el vertimiento, la concentración de salida de OD sigue siendo favorable.

En segundo lugar, se evaluó el fosforo total que se encuentra en el agua residual en forma de
fosfatos (Metcalf & Eddy, 2014), su principal fuente son los detergentes debido al uso inadecuado
de los productos de limpieza y desinfección utilizados dentro de las instalaciones de la heladería
de san Jerónimo S.A, contribuyendo con el crecimiento desmesurado de las plantas causando una
mala calidad del agua, ya que, la concentración de entrada de P Total es de 3,5 mg/L dando como
resultado: Cs (0,17 mg/L) > (0,1 mg/L) Acuerdo 43, incumpliendo con el estándar de calidad del
agua como se puede observar en la grafica 13; Por lo tanto, se redujo la concentración de P Total
en un 92% y el caudal en un 30% para que la Cs (0,1 mg/L) = (1,0 mg/L) Acuerdo 43.(gráfica
14) Lo que da lugar a que la carga contamínate diaria máxima soportada de P total es de
2,59 × 10−5 kg/día, lo que representa una reducción del 41,17 % en comparación con escenario
actual.

En tercer lugar, se tuvo en cuenta los nitritos y los nitratos que hacen parte natural del ciclo del
nitrógeno y es fundamental para el desarrollo de las plantas, por ende, el fosforo y el nitrógeno
contribuyen a la eutrofización. La oxidación del nitrógeno amoniacal a nitritos y nitrato, consume
oxígeno y, por lo tanto, puede agotar seriamente los niveles de oxigeno del cuerpo de agua (Chin,
2013); La concentraciones de nitritos y nitratos del vertimiento de interés corresponde a 0,14 mg/L
y 90 mg/L, respectivamente; Estos valores sobrepasan las concentraciones del acuerdo 43, sin
embargo, para los nitritos la Cs (0,0067 mg/L) < (1 mg/L) Acuerdo 43 .(gráfica 16) y para los
nitratos la Cs (2,25 mg/L) < (1 mg/L) Acuerdo 43. . (gráfica 15)

En base a todo lo anterior, se puede decir que el río alto Ubaté tiene una buena asimilación de los
parámetros de interés del vertimiento de la heladería de San jerónimo S.A, a excepción de la DBO
y el T total. Sin embargo, se puede mejorar el sistema de tratamiento de las aguas residuales para
cumplir con las reducciones propuestas y reevaluar el permiso de vertimientos; Por otra parte, los
valores establecidos en el acuerdo 43 de 2007, son acordes para representar los objetivos de calidad
del agua del río alto Ubaté, ya que, no tiene tanta intervención humana.
 12

CONCLUSIONES:
• El río alto Ubaté, en las condiciones actuales muestra valores de OD en casi todo el tramo
de estudio encima de 4 mg/l, lo que indica condiciones de un rio saludable; Sin embargo,
después de la descarga del río Suta, la concentración de OD es un poco preocupante, ya
que, este punto puede llegar a ser seriamente contaminado si no se toman las medidas
adecuadas.
• Con la evaluación de escenarios TMDL, puede decirse que hay una interacción aceptable
en materia de nutrientes en el rio alto de Ubaté, incluso con el vertimiento estudiado.
• Para la recuperación del río alto Ubaté se requiere el mejoramiento de sus afluentes, en
especial del río Suta, el cual es el receptor de los sistemas de alcantarillado de las
poblaciones existentes en dicha cuenca.

BIBLIOGRAFIA:
[1]Montañez Diana, (2019); ANÁLISIS DE LA CALIDAD FISICOQUÍMICA DE LA LAGUNA DE
FÚQUENE PROPUESTA DE LINEAMIENTOS PARA LA GESTIÓN AMBIENTAL DEL
RECURSO HÍDRICO; Pontificia Universidad Javeriana. Obtenido de:
https://repository.javeriana.edu.co/bitstream/handle/10554/42995/Definitivo%20BIBLIOTECA
%20%20DIANA%20PAOLA%20MONTA%C3%91EZ%20QUIROGA%20TESIS%20LAGUN
A%20DE%20FUQUENE%20MGA.pdf?sequence=1&isAllowed=y

[2] Corporación autónoma regional de Cundinamarca (CAR); DIAGNÓSTICO PROSPECTIVA Y

FORMULACIÓN DE LA CUENCA HIDROGRÁFICA DE LOS RÍOS UBATÉ Y SUÁREZ. Santa
Fe Bogotá D.C. Obtenido de:
https://www.car.gov.co/uploads/files/5ac692fb56934.pdf

[3] CORPORACIÓN AUTÓNOMA REGIONAL – CAR. Boletín de Calidad de Las Cuencas de

la Jurisdicción CAR 2007 – 2009. Bogotá: CAR. 2010. Obtenido de:
https://www.car.gov.co/uploads/files/5ada3cc8b8754.pdf

[4] ArcGIS Web Application. (2021). Ideam.gov.co. Obtenido de:

http://dhime.ideam.gov.co/atencionciudadano/

[5]Corporación autónoma regional de Cundinamarca (CAR); CUENCA RÍO ALTO UBATE 2401-
02; Informe de la fase de diagnóstico contrato de consultoría No.800 de 2005. Obtenido
de: https://www.car.gov.co/uploads/files/5ac693ff7f486.pdf
 13

[6]] Corporación autónoma regional de Cundinamarca (CAR); CUENCA RÍO ALTO SUÁREZ
2401-10; Informe de la fase de diagnóstico contrato de consultoría No.800 de 2005. Obtenido
de:https://www.car.gov.co/uploads/files/5ac694dfa1c76.pdf

[7] (s.a) (s.f); Apéndice E: Calidad del agua y mecanismos de polución

https://openjicareport.jica.go.jp/pdf/11572898_11.pdf
[8] Google Earth Pro (2021)
[9] Ruiz Esparza, Ángela Carolina. et al ; Diagnóstico ambiental industrial y guía ambiental
para las industrias lecheras; Universidad libre
https://repository.unilibre.edu.co/bitstream/handle/10901/10917/Cuerpo%20del%20proyecto.pdf
?sequence=3&isAllowed=y

[10] Falak Naz. et al (2011); EFFECT OF DIFFERENT LEVELS OF NPK FERTILIZERS ON

THE PROXIMATE COMPOSITION OF POTATO CROP AT ABBOTTABAD; Department of Plant
Protection, Agricultural University Peshawar. pdf

[11] Spyridon A. Petropoulos. et al (2020); The Impact of Fertilization Regime on the Crop
Performance and Chemical Composition of Potato (Solanum tuberosum L.) Cultivated in Central
Greece. pdf

[12] Goméz Manuel, et al (2017); Accumulation of N, P, and K in the tubers of potato (Solanum
tuberosum L. ssp. andigena) under contrasting soils of the Andean region of Colombia. Obtenido
de: http://www.scielo.org.co/scielo.php?script=sci_arttext&pid=S0120-99652017000100059

[13] Hidrología Superficial (III): Relación Precipitación -Escorrentía. (n.d.).

https://hidrologia.usal.es/temas/Hid_Sup_3.pdf

[14] EPA. (2020). Indicators: Dissolved Oxygen.

https://www.epa.gov/national-aquatic-resource-surveys/indicators-dissolved-oxygen

[15] Chapra, Steven C. Surface Water-Quality Modeling. Reissued. Long Grove, Ill: Waveland
Press, 2008.
[16] David A. Chin (2013); Water-quality engineering in natural systems: fate and transport
processes in the water environment. Second Edition.
[17] Science for a changing world (USGS), Dissolved oxygen solubility tables DOTABLES
https://water.usgs.gov/cgi-bin/dotables
[18] Bowie, George L et al. (1985); RATES, CONSTANTS, AND KINETICS FORMULATIONS
IN SURFACE WATER QUALITY MODELING (SECOND EDITION); EPA/600/3-85/040. pdf
 14

[19] Sastre Torres, Diana Alejandra, (2016); MODELO DE GESTIÓN PARA EL TRATAMIENTO
DE AGUAS RESIDUALES DENTRO DEL ÁREA DE PRESTACIÓN DEL SERVICIO DE
ALCANTARILLADO EN EL MUNICIPIO DE CAJICÁ – CUNDINAMARCA. Universidad de la
Sabana.
https://intellectum.unisabana.edu.co/bitstream/handle/10818/30825/Diana%20Alejandra%20Sast
re%20Torres%20%28Tesis%29.pdf?sequence=1&isAllowed=y

[20] Corporación autónoma regional de Cundinamarca (CAR); Resolución 3462 de 28 de

diciembre de 2009. https://www.car.gov.co/uploads/files/5ac79c5b18c9e.pdf

[21] Ministerio de ambiente y desarrollo sostenible; Resolución 631 de 17 de marzo de 2015.

http://www.emserchia.gov.co/PDF/Resolucion631.pdf

[22] United States Environmental Protection Agency (EPA,2009); New River QUAL2K Water
Quality Model for the New River Dissolved Oxygen TMDL; San Francisco, California. PDF

[23] Metcalf & Eddy. (2014). Wastewater Engineering: Treament and Resource Recovery. Ed. Mc
Graw Hill. Education. PDF.

[24] Habes Suad, et.al (2015); Microbiological quality of processed water during processing of
milk and dairy products. Balkan Journal of Health Science.
https://www.researchgate.net/publication/303524323_Microbiological_quality_of_process
ed_water_during_processing_of_milk_and_dairy_products
[25] Lestari, et.al (2019); Determination of Pollution Load Capacity Using QUAL2Kw Program
on The Musi River Palembang; Environment and Water Resources Division. PDF