IMPACTO DE LOS ESCENARIOS DE CAMBIO CLIMÁTICO EN LA RELACION

OFERTA-DEMANDA EN LA CUENCA DEL RÍO BALSILLAS - DEPARTAMENTO DE

CUNDINAMARCA.

AUTORES

Natalia Laverde Alarcón

Nicole Fernanda Moyano Martínez

UNIVERSIDAD SANTO TOMAS

FACULTAD DE INGENIERIAS
INGENIERIA AMBIENTAL
BOGOTA D.C.
2019

1
 IMPACTO DE LOS ESCENARIOS DE CAMBIO CLIMÁTICO EN LA RELACION
OFERTA-DEMANDA EN LA CUENCA DEL RÍO BALSILLAS - DEPARTAMENTO DE
CUNDINAMARCA.

NATALIA LAVERDE ALARCÓN

NICOLE FERNANDA MOYANO MARTÍNEZ

TRABAJO PARA OPTAR AL TITULO DE INGENIERIA AMBIENTAL

DIRECTOR
Darwin Mena Rentería
Ingeniero Ambiental y Sanitario
MSc en Evaluación de Recursos Hídricos

CO-DIRECTOR
Miguel Angel Cañon Ramos
Ingeniero Ambiental
MSc (c) Hidrosistemas

UNIVERSIDAD SANTO TOMAS

FACULTAD DE INGENIERIAS
INGENIERIA AMBIENTAL
BOGOTA D.C.
2019

2
 Nota de aceptación

_________________________________________
_________________________________________
_________________________________________
_________________________________________
_________________________________________
_________________________________________

_________________________________________
Firma director trabajo de grado

_________________________________________
Firma del jurado

_________________________________________
Firma del jurado

3
 AGRADECIMIENTOS

Agradecemos a Dios por ser nuestro guía a lo largo de nuestra carrera, ser el apoyo y
fortaleza en aquellos momentos de dificultad y de debilidad, por brindarnos paciencia y
sabiduría para culminar con éxito nuestras metas propuestas
Gracias a nuestros padres: Liliana Alarcón, Stanly Laverde, Gilma Martínez y Anatolio
Moyano, por ser los principales promotores de nuestros sueños, por confiar y creer en
nuestras expectativas, por los consejos, valores y principios que nos han inculcado y ante
todo por su apoyo incondicional.
Agradecemos de manera especial, al Docente Darwin Mena Rentería Director de nuestro
proyecto de investigación quien con su experiencia, conocimiento y motivación nos ha
orientado, así mismo por su paciencia y por la confianza que depositó en nosotras para
culminar satisfactoriamente este proyecto.
Agradecemos a nuestros docentes de la Facultad de Ingeniería Ambiental de la Universidad
Santo Tomas, por haber compartido sus conocimientos a lo largo de la preparación de
nuestra profesión, especialmente a Miguel Cañón, como Codirector y cuyo papel fue de
vital importancia en el desarrollo de este proyecto de investigación

4
 TABLA DE CONTENIDO

1. RESUMEN ............................................................................................................... 11
2. ABSTRACT ............................................................................................................. 12
3. INTRODUCCIÓN ..................................................................................................... 13
4. OBJETIVOS ............................................................................................................ 14
4.1. Objetivo general ................................................................................................. 14
4.2. Objetivos específicos ........................................................................................ 14
5. PLANTEAMIENTO DEL PROBLEMA ..................................................................... 15
6. JUSTIFICACIÓN ..................................................................................................... 16
7. MARCO TEÓRICO .................................................................................................. 17
7.1. Descripción general de la cuenca ..................................................................... 17
7.2. Topografía y morfometría .................................................................................. 19
7.3. Climatología ....................................................................................................... 21
7.4. Recursos hídricos e hidrología ......................................................................... 30
7.5. Cambio climático ............................................................................................... 30
7.6. Escenarios de cambio climático ....................................................................... 31
7.7. Modelos de recirculación global (GCM) ........................................................... 33
7.8. Reducción a escala ............................................................................................ 34
7.8.1. Reducción a escala basada en la teoría del caos .................................... 35
7.9. Experimento regional coordinado de reducción climática (CORDEX) ........... 37
7.10. Sistema de indicadores hídricos....................................................................... 38
7.10.1. Índice de retención y regulación hídrica (IRH) ......................................... 39
7.10.2. Índice de uso del agua superficial (IUA) ................................................... 39
7.10.3. Índice de vulnerabilidad por desabastecimiento (IVH) ............................ 39
7.11. Modelación hidrológica ..................................................................................... 39
7.11.1. Hydro-BID.................................................................................................... 40
7.11.1.1. Calibración del modelo con el software Hydro-BID ......................... 41
8. METODOLOGÍA ...................................................................................................... 43
9. RESULTADOS ........................................................................................................ 44
9.1. Estimación de la oferta ...................................................................................... 44
9.1.1. Definición y procesamiento de la cuenca ................................................. 44
9.2. Obtención y tratamiento de los datos de entrada ............................................ 45
9.3. Calibración y validación con el software Hydro-BID ....................................... 48
9.4. Análisis cambio climático.................................................................................. 52
9.4.1. Modelo de circulación global (GCM) ......................................................... 52
9.4.2. Reducción de escala basada en la teoría del Caos .................................. 57
10. OBTENCIÓN DE LA DEMANDA ............................................................................. 64
10.1. Obtención de la demanda actual y futura ......................................................... 64
10.1.1. Sector agrícola ........................................................................................... 66
10.1.2. Sector pecuario .......................................................................................... 67
10.1.3. Sector industrial ......................................................................................... 68
10.1.4. Sector domestico ....................................................................................... 69
11. ANÁLISIS DE LA VARIABILIDAD DE OFERTA HÍDRICA NATURAL ................... 70
11.1. Curva de duración de caudales......................................................................... 70
11.2. Índice de retención y regulación hídrica (IRH) ................................................. 72
11.2.1. Caudal ambiental ........................................................................................ 74

5
11.3. Índice de uso de agua (IUA) .............................................................................. 74
11.3.1. Demanda hídrica......................................................................................... 75
11.3.2. Oferta hídrica .............................................................................................. 76
11.4. Índice de vulnerabilidad por desabastecimiento (IVH) .................................... 77
12. IMPACTO SOCIAL .................................................................................................. 78
13. CONCLUSIONES .................................................................................................... 79
14. RECOMENDACIONES ............................................................................................ 81
15. BIBLIOGRAFÍA ....................................................................................................... 82
16. ANEXOS.................................................................................................................. 87
ANEXO A. Comparación de datos históricos y datos futuros para el modelo
CORDEX MPI-ESM-MR para la variable de precipitación ............................................ 87
ANEXO B. Comparación de datos históricos y datos futuros para el modelo GCM
MPI-ESM-MR para la variable de temperatura ............................................................. 90
ANEXO C. Puntos de demanda en la cuenca por municipios y sectores según
registros de la CAR ....................................................................................................... 93
ANEXO D. Calculo de la demanda por sectores económicos .................................. 133

LISTA DE TABLAS

Tabla 1 Distribución territorial dentro de la subcuenca .............................................. 18

Tabla 2 Porcentaje de población que se ubica en la Cuenca del Rio Balsillas ......... 18
Tabla 3 Coberturas y especies forestales en la Cuenca del Rio Balsillas ................. 19
Tabla 4 Características del suelo de la Cuenca del Rio Balsillas ............................... 19
Tabla 5 Distribución de los diferentes rangos de pendiente de la Cuenca del Rio
Balsillas. ......................................................................................................................... 20
Tabla 6 Características morfométricas de la Cuenca del Rio Balsillas...................... 20
Tabla 7 Comparación de ventajas y desventajas de cada tipo de técnica ................ 35
Tabla 8 Parámetros de calibración del software Hydro-BID ....................................... 41
Tabla 9 Métricas de desempeño para la calibración del software Hydro-BID ........... 42
Tabla 10 Datos de entrada del modelo en el software Hydro-BID. ............................. 45
Tabla 11 Parámetros de calibración para el modelo en Hydro-BID. ........................... 48
Tabla 12 Métricas de desempeño para la calibración del modelo en Hydro-BID. ..... 48
Tabla 13 Métricas de desempeño para la validación del modelo en Hydro-BID. ...... 51
Tabla 14. Resultados criterios de comparación MPI-ESMMR –CORDEX ................... 53
Tabla 15. Resultados criterios de comparación MPI-ESM-MR ................................... 55
Tabla 16 Variación de la precipitación bajo la reducción de escala. ......................... 58
Tabla 17.Variación de la temperatura bajo la reducción de escala ............................ 61
Tabla 18 Calculo de la población proyectada hasta el 2020 por Municipios. ............ 65
Tabla 19 Clasificación de demanda hasta el año 2020 por sectores en litro/año..... 65
Tabla 20 Demanda de agua en m3/s del sector agrícola para cada municipio ......... 66
Tabla 21 Demanda de agua en m3/s del sector pecuario para cada municipio ........ 67
Tabla 22 Demanda de agua en m3/s del sector industrial para cada municipio ....... 68
Tabla 23 Determinación de dotación neta................................................................... 69
Tabla 24 Demanda de agua en m3/s del sector doméstico para cada municipio ..... 69

6
Tabla 25 Rangos y categorías del Índice de Retención y Regulación Hídrica ......... 73
Tabla 26 Índice de retención y regulación hídrica para la cuenca del Rio Balsillas 73
Tabla 27 Caudal ambiental en m3/s para cada periodo de tiempo en la Cuenca del Rio
Balsillas .......................................................................................................................... 74
Tabla 28 Demanda total de agua en m3/s para cada periodo de tiempo ................... 75
Tabla 29 Oferta hídrica para cada periodo de tiempo en la Cuenca del Rio Balsillas
........................................................................................................................................ 76
Tabla 30 Rangos y categorías del Índice del Uso del Agua ....................................... 76
Tabla 31 Índice de uso de agua para la cuenca del Rio Balsillas .............................. 77
Tabla 32 Índice de vulnerabilidad al desabastecimiento hídrico IVH ........................ 77
Tabla 33 Índice de uso de agua extremo..................................................................... 78
Tabla 34 Índice de vulnerabilidad por desabastecimiento del Rio Balsillas............. 78

LISTA DE ILUSTRACIONES

Ilustración 1 Descripción de los escenarios de cambio climático ............................. 32

Ilustración 2 Evolución de la resolución espacial de los GCM a través del tiempo.. 33
Ilustración 3 Esquema del concepto downscalling espacial ...................................... 34
Ilustración 4 Comparación del espacio de fase de la serie de tiempo de precipitación
........................................................................................................................................ 36
Ilustración 5 Relación entre los puntos de los atractores .......................................... 37
Ilustración 6 Sistema Hydro-BID y sus componentes................................................. 41
Ilustración 7 Esquema metodológico........................................................................... 43
Ilustración 8 Datos observados de precipitación Vs caudales observados y simulados
........................................................................................................................................ 50

LISTA DE MAPAS

Mapa 1 Localización de la Cuenca del Rio Balsillas ................................................... 17

Mapa 2 Elevación medio de la Cuenca del Rio Balsillas............................................. 21
Mapa 3 Temperatura media en °C de la Cuenca del Rio Balsillas ............................. 23
Mapa 4 Precipitación media en mm de la Cuenca del Rio Balsillas .......................... 24
Mapa 5 Humedad Relativa en % de la Cuenca del Rio Balsillas ................................ 26
Mapa 6 Brillo solar en horas de la Cuenca del Rio Balsillas ..................................... 28
Mapa 7 Delimitación de la Cuenca del Rio Balsillas con el software QGIS ............... 45
Mapa 8 Ubicación de las estaciones de temperatura dentro de la Cuenca del Rio
Balsillas .......................................................................................................................... 46
Mapa 9 Ubicación de las estaciones de precipitación dentro de la Cuenca del Rio
Balsillas .......................................................................................................................... 47
Mapa 10 Ubicación de las estaciones de caudal dentro de la Cuenca del Rio Balsillas
........................................................................................................................................ 47

7
Mapa 11 Interpolación de precipitación bajo el escenario RCP 2.6 ........................... 59
Mapa 12 Interpolación de precipitación bajo el escenario RCP 4.5 ........................... 60
Mapa 13 Interpolación de precipitación bajo el escenario RCP 8.5 ........................... 60
Mapa 14 Interpolación de temperatura bajo el escenario RCP 2.6 ............................. 63
Mapa 15 Interpolación de temperatura bajo el escenario RCP 4.5 ............................. 63
Mapa 16 Interpolación de temperatura bajo el escenario RCP 8.5 ............................. 64

LISTA DE GRAFICAS

Grafica 1 Temperatura media en °C de la Cuenca del Rio Balsillas. .......................... 23

Grafica 2 Precipitación media en mm de la Cuenca del Rio Balsillas. ....................... 25
Grafica 3 Diagrama de cajas de precipitación de la Cuenca del Rio Balsillas .......... 25
Grafica 4 Promedio mensual de Humedad relativa de la Cuenca del Rio Balsillas .. 27
Grafica 5 Diagrama de cajas de Humedad relativa de la Cuenca del Rio Balsillas ... 27
Grafica 6 Promedio mensual de Brillo solar de la Cuenca del Rio Balsillas ............. 29
Grafica 7 Diagrama de cajas de Brillo Solar de la Cuenca del Rio Balsillas ............. 29
Grafica 8 Curva de duración de caudales de la calibración del modelo ................... 49
Grafica 9 Curva de duración de caudales de la validación del modelo .................... 51
Grafica 10 Comportamiento mensual multianual GCM vs datos históricos de
precipitación de la estación 301388800 ....................................................................... 53
Grafica 11 Comportamiento mensual multianual GCM vs datos históricos de
precipitación de la estación 301388900 ....................................................................... 54
Grafica 12 Comportamiento mensual multianual GCM vs datos históricos de
precipitación de la estación 301413300 ....................................................................... 54
Grafica 13 Comportamiento mensual multianual GCM vs datos históricos de
temperatura de la estación 301388800 ......................................................................... 55
Grafica 14 Comportamiento mensual multianual GCM vs datos históricos de
temperatura de la estación 301388900 ......................................................................... 56
Grafica 15 Comportamiento mensual multianual GCM vs datos históricos de
temperatura de la estación 301413300 ......................................................................... 56
Grafica 16 Precipitación promedio - Datos históricos VS series de escala reducida
........................................................................................................................................ 58
Grafica 17 Temperatura promedio - Datos históricos VS series de escala reducida62
Grafica 18 Curva de duración de caudales medios diarios del año 2015 a 2040 ..... 71
Grafica 19 Curva de duración de caudales medios diarios del año 2041 a 2070 ..... 71
Grafica 20 Curva de duración de caudales medios diarios del año 2071 a 2100 ..... 72

8
 LISTA DE ECUACIONES

Ecuación 1 Balance hídrico. ......................................................................................... 30

Ecuación 2 Cálculo de μ................................................................................................ 36
Ecuación 3 Cálculo de la demanda doméstica. ........................................................... 66
Ecuación 4 Cálculo de probabilidad de excedencia.................................................... 71
Ecuación 5 Cálculo del índice de retención hídrica. ................................................... 73
Ecuación 6 Cálculo de índice de uso de agua. ............................................................ 74
Ecuación 7 Cálculo de demanda hídrica sectorial. ..................................................... 75
Ecuación 8 Cálculo de oferta hídrica superficial disponible. ..................................... 76

9
 ABREVIATURAS

AHD: Base de Datos de Hidrología Analítica

ALC: América Latina y El Caribe
BID: Banco Interamericano de Desarrollo
CMIP: Proyecto de comparación de modelos acoplados
CMNUCC: Convención Marco de Las Naciones Unidas sobre el Cambio Climático
DANE: Departamento Administrativo Nacional de Estadística
ENA: Estudio Nacional del Agua
ENSO: Oscilación del Niño Sur
GCM: Modelo de Circulación Global
GEI: Gases de Efecto Invernadero
IA: Índice de aridez
IARC: Índice de agua no retornada a la cuenca
IDEAM: Instituto de Hidrología, Meteorología y Estudios Ambientales
IIASA: Instituto Internacional para el Análisis de Sistemas Aplicados
IVH: Índice de vulnerabilidad al desabastecimiento hídrico
IUA: Índice de uso del agua
IPCC: Panel Intergubernamental de Cambio Climático
IRH: Índice de retención hídrica
IPH: Índice de presión hídrica a los ecosistemas
JGCRI: Joint Global Change Research Institute
m.s.n.m: Metros sobre el nivel del mar
mm: Milímetros de precipitación
OMM: Organización meteorológica mundial
PIB: Producto Interno Bruto
PNUMA: Programa de las Naciones Unidas para el Medio Ambiente
POMCA: Plan de Ordenamiento y Manejo de Cuencas
RCM: Modelo de Circulación Regional
USGS: Registros del servicio geológico de los Estados Unidos
°C: Grados centígrados

10
 1. RESUMEN

En este trabajo de investigación se efectuó un análisis de la relación entre la oferta y la

demanda de la cuenca del río Balsillas ubicado en la Sabana de Bogotá, incluyendo las
variaciones de precipitación y temperatura asociados a los escenarios de cambio climático
planteados para Colombia; las cuencas andinas como la del río Balsillas son vulnerables a
los cambios del clima, por lo que la probabilidad a no satisfacer las demandas de los
diferentes sectores hidrodependientes (como, los sistemas de abastecimiento domésticos
e industriales y la agricultura), tiende a acentuarse, por efecto de la posible reducción en la
oferta y por una mayor demanda de agua debido al crecimiento de la población y de las
actividades económicas.

Para dichos efectos de no abastecimiento se estimó el régimen hidrológico futuro a partir

de la herramienta de modelación Hydro-BID, la cual permite estimar la disponibilidad de
agua, así como incorporar las variaciones del clima futuro de una zona, con el fin de
identificar los impactos de los escenarios de cambio climático sobre los recursos hídricos a
nivel de cuenca.

Para los análisis de la relación oferta-demanda de agua y de la demanda insatisfecha

obtenida de dicha modelación, se realizó el análisis de balances hídricos mensuales en los
nodos de demanda y oferta utilizando modelos de gestión del agua, los cuales permiten el
análisis de una serie de escenarios e incertidumbres definidas por el usuario, incluyendo
las relacionadas con el clima, las condiciones de las cuencas, las demandas proyectadas,
las necesidades de los ecosistemas, cambios en los usos del suelo, y las necesidades de
infraestructura.

Palabras clave: Hydro-BID, oferta-demanda hídrica, Rio Balsillas, cambio climático,

recurso hídrico

11
 2. ABSTRACT

In this research work an analysis was made of the relationship between supply and demand
of the Balsillas river basin located in the Sabana of Bogota, including precipitation and
temperature variations associated with climate change scenarios in Colombia; Andean
basins such as the Balsillas River are vulnerable to climate change, the probability of not
meeting the demands of the different hydro-dependent sectors (such as domestic and
industrial supply systems and agriculture) tends to increase, the effect of a possible
reduction in supply and increased demand for water due to population growth and economic
activities.

For these effects of no supply, the future hydrological regime was estimated from the Hydro-
BID modelling tool, which allows to estimate the availability of water, as well as to incorporate
the variations of the future climate of an area, in order to identify the impacts of climate
change scenarios on water resources at the basin level.

For the analysis of the supply-demand ratio of water and the unsatisfied demand obtained
from such modelling, the analysis of monthly water balances at demand and supply nodes
was carried out using water management models, which allow the analysis of a number of
user-defined scenarios and uncertainties, including those related to climate, watershed
conditions, projected demands, ecosystem needs, land use changes, and infrastructure
needs.

Keywords: Hydro-BID, water supply - demand, Balsillas river, climate change, hydric
resource.

12
 3. INTRODUCCIÓN

El agua contribuye a la estabilidad del funcionamiento del entorno y de los seres y

organismos que en él habitan, es, por tanto, un elemento indispensable para la subsistencia
de la vida animal y vegetal del planeta. "El agua es un bien de primera necesidad para los
seres vivos y un elemento natural imprescindible en la configuración de los sistemas
medioambientales" [¡Error! No se encuentra el origen de la referencia.], por ende la
contaminación del agua y su escasez plantean amenazas para la salud humana y la calidad
de vida; razón por la cual, se ha convertido en un tema de interés para el sostenimiento de
los ecosistemas tanto húmedo como terrestre, sometiendo a una presión todavía mayor a
la flora y la fauna, que padecen ya las repercusiones de la urbanización y el cambio
climático [¡Error! No se encuentra el origen de la referencia.].

Teniendo en cuenta que Colombia se ubica entre los países con mayor riqueza en recursos
hídricos en el mundo, se han realizado estimaciones por medio del IDEAM, mediante las
cuales se afirma que la demanda de uso de agua para el desarrollo de las actividades
socioeconómicas en el país se representa principalmente mediante los siguientes usos:
agrícola, doméstico, industrial, pecuario y servicios; el uso que presenta un mayor
porcentaje es el agrícola con el 54%, con el 29% el doméstico y con el 13% el industrial, en
menor escala el pecuario y el de servicios con porcentajes del 3%, y el 1% respectivamente,
en las cuales existen necesidades hídricas insatisfechas de los ecosistemas y que cada vez
es mayor el número de impactos de origen antrópico sobre el agua, por lo que la
disponibilidad del recurso es cada vez menor [¡Error! No se encuentra el origen de la
referencia.].

Con base a esta información y estudios realizados a nivel nacional para los diferentes
sectores económicos esta investigación emplea herramientas de modelación hidrológica y
de planificación del recurso hídrico, con el fin de evaluar el impacto del cambio climático
sobre la disponibilidad hídrica y su repercusión en la demanda de la cuenca del río Balsillas.
Se implementó la herramienta Hydro-BID para generar las series de tiempo de caudales
bajo el clima actual y bajo las proyecciones climáticas futuras, así como la información
planteada en el documento titulado Estudio Nacional del Agua (ENA).

13
 4. OBJETIVOS

4.1. Objetivo general

Analizar la relación oferta-demanda de agua en la cuenca del río Balsillas bajo escenarios
de cambio climático a partir de indicadores.

4.2. Objetivos específicos

 Evaluar la disponibilidad de agua de la cuenca del río Balsillas bajo escenarios de

cambio climático implementando la herramienta Hydro-BID.

 Estimar la demanda actual y futura de agua en la cuenca del río Balsillas para los
diferentes sectores económicos.

 Estimar los indicadores de sostenibilidad del recurso hídrico para la cuenca del Rio
Balsillas.

14
 5. PLANTEAMIENTO DEL PROBLEMA

Millones de personas sufren ya los efectos catastróficos de desastres naturales extremos

debido al cambio climático que no sólo conlleva un aumento de las temperaturas, sino
también fenómenos meteorológicos extremos, la elevación del nivel del mar y cambios en
las poblaciones y los hábitats de flora y fauna silvestres, entre otros efectos [¡Error! No se
encuentra el origen de la referencia.].

Actuar frente al cambio climático requiere la acción inmediata de todos los países del
mundo, Colombia ha sido uno de los 195 países que firmó el primer acuerdo mundial
vinculante sobre el clima “Acuerdo de París”, el cual plantea que de aquí a 2030 debemos
reducir las emisiones de gases con efecto invernadero a la mitad de su nivel en 2010 para
no alcanzar el aumento de 1,5 °C. Por tanto, nuestros gobiernos deben tomar medidas de
forma inmediata para dar un cambio de rumbo [¡Error! No se encuentra el origen de la
referencia.]. En cuanto a la gestión del recurso hídrico Colombia desarrollo medidas de
protección en respuesta a una creciente necesidad de lograr la conservación y el uso
sostenible del agua. Bajo esta premisa, es evidente como en los últimos años se han
intensificado los esfuerzos por el planteamiento de políticas encaminadas al manejo del
agua en el país. Sin embargo, existen limitantes como la falta de conocimiento y
caracterización del ciclo hidrológico en las diferentes cuencas y la variación del régimen
hidrológico como consecuencia de los inminentes cambios del clima [¡Error! No se
encuentra el origen de la referencia.].

Como respuesta a la falta de información identificada se han desarrollado herramientas,

como la modelación hidrológica mediante el software Hydro-BID, que contribuye a obtener
un mayor conocimiento sobre los cuerpos hídricos [¡Error! No se encuentra el origen de
la referencia.]. Este software “Hydro-BID” presenta varias ventajas como conocer las
condiciones actuales y futuras del recurso, evaluar el impacto del cambio climático sobre la
disponibilidad hídrica su grado de vulnerabilidad, prevención y reducción del riesgo
ambiental y social en las diferentes zonas, esto según la Ley 1523 de 2012 la cual rige
mediante el articulo 3 por el principio de protección donde especifica que los residentes de
Colombia deben ser protegidos por las autoridades competentes frente a posibles desastres
o fenómenos peligrosos que amenacen o infieran a la integridad física y mental de la
persona. En adición a esto, también facilita la caracterización en cuencas con escasez de
información o sin información como condiciones de difícil acceso o demasiado extensas de
una forma más eficiente y sencilla; con menor inversión de dinero y tiempo. De esta manera
es posible lograr una mejor gestión del recurso cumpliendo los objetivos planteados en la
Política Nacional para la gestión del Recurso Hídrico [¡Error! No se encuentra el origen
de la referencia.].

15
 6. JUSTIFICACIÓN

El análisis del cambio climático y sus repercusiones representan una preocupación

persistente a nivel global que ha llevado a los organismos gubernamentales a realizar un
esfuerzo por hacer frente a los efectos adversos; planteando acciones de mitigación y
adaptación y fortaleciendo el desarrollo de investigación [¡Error! No se encuentra el origen
de la referencia.]. En el caso particular de Colombia, a pesar de ser una nación con
economía en desarrollo, que no genera grandes aportaciones de gases de efecto
invernadero a la atmósfera, se ha determinado que hace parte de un conjunto de naciones
que se consideran más vulnerables a los efectos del cambio climático [¡Error! No se
encuentra el origen de la referencia.].

Según estudios realizados por parte del grupo de modelamiento de tiempo, clima y
escenarios de cambio climático perteneciente al Instituto de Hidrología, Meteorología y
Estudios Ambientales (IDEAM) la precipitación actual presenta un aumento de entre el 10%
y el 40% comparado con las décadas anteriores, además la temperatura presenta un
comportamiento acorde a la región, con anomalías entre 0.5 y 1.5 [¡Error! No se encuentra
el origen de la referencia.]. Adicionalmente, este instituto ha estimado que para los
escenarios de cambio climático planteados en el periodo 2011-2040 se prevé un aumento
de la precipitación que varía entre el 10% y el 30% [¡Error! No se encuentra el origen de
la referencia.].

Debido a la heterogeneidad del territorio colombiano, es necesario el desarrollo de

investigaciones que impliquen análisis a nivel regional y local, para poder determinar de
manera más concisa la vulnerabilidad, reducir los posibles riesgos por sequías o
inundaciones y asegurar la disponibilidad de agua para abastecimiento [¡Error! No se
encuentra el origen de la referencia.].

El efecto de estas modificaciones en el comportamiento de la precipitación y la temperatura

en el recurso hídrico se ha estudiado a escala nacional y de manera muy general, para el
caso de la cuenca del río Balsillas ubicada en la Sabana de Bogotá, el Plan de Ordenación
y Manejo Integral de Cuencas (POMCA) estableció los Índices de vulnerabilidad y
desabastecimiento hídrico (IVH) que mide el grado de fragilidad del sistema hídrico para
mantener una oferta en el abastecimiento de agua, índice de uso del agua (IUA) el cual
determina la cantidad de agua utilizada por los diferentes sectores usuarios y el Índice de
Retención y Regulación Hídrica (IRH) que mide la cantidad de humedad que pueden retener
las cuencas. Según el estudio nacional del agua el cual mide el grado de fragilidad del
sistema hídrico para mantener una oferta en el abastecimiento de agua, que ante amenazas
como periodos largos de estiaje o eventos como el fenómeno cálido del Pacífico podría
generar riesgos de desabastecimiento [¡Error! No se encuentra el origen de la
referencia.]. Sin embargo, no se determinó la vulnerabilidad con respecto a los escenarios
de cambio climático, por lo cual es muy difícil la aplicación de estrategias de gestión,
adaptación y mitigación, esto basado en la política nacional de cambio climático y el Panel
Intergubernamental de Cambio Climático, los cuales establecen que la mitigación es la
intervención que realiza el hombre para reducir los impactos o mejorar las condiciones de
calidad ambientales y condiciones de vida a las comunidades [¡Error! No se encuentra el
origen de la referencia.].

16
El planteamiento de medidas para la reducción de los impactos del cambio climático sobre
los recursos hídricos requieren tener el conocimiento del comportamiento futuro de los
sistemas de recursos hídricos, por lo tanto los resultados obtenidos de este proyecto servirán
de base para que las autoridades encargadas del manejo y administración del recurso
hídrico, así como los usuarios cuenten con las herramientas técnico-científicas para afrontar
los impactos que se puedan presentar en un futuro, debido al cambio en las condiciones
climáticas de la zona.

7. MARCO TEÓRICO

7.1. Descripción general de la cuenca

La cuenca del río Balsillas se encuentra ubicada al este de la Sabana y en la zona centro
occidental de la cuenca media del Rio Bogotá, en el departamento de Cundinamarca en los
municipios de Albán, Anolaima, Bojaca, El Rosal, Facatativá, Funza, La Vega, Madrid,
Mosquera, Pacho, San Francisco, Sasaima, Soacha, Subachoque, Tabio, Tenjo, Zipaquirá;
a una altura 2730 msnm con una temperatura de 13°C aproximadamente lo cual permite la
presencia de climas fríos, húmedos y secos [¡Error! No se encuentra el origen de la
referencia.]. Posee un área de 62441,6 Ha, con un cauce principal de 68,1 Km de longitud
abarcando una población cercana a los 197486 habitantes. Tiene como tributarios
importantes el río Subachoque, el río Bojacá, las quebradas San Joaquín, Yerbabuena,
Pozo Amarillo y Los Curubitos [¡Error! No se encuentra el origen de la referencia.].

Mapa 1 Localización de la Cuenca del Rio Balsillas

17
 Fuente Autores

La cuenca del Rio Balsillas hace parte de la Cuenca del Rio Bogotá, a continuación, se
presentan los municipios que hacen parte de esta cuenca y así mismo se especifican las
áreas en hectáreas (Ha) que están dentro esta:

Tabla 1 Distribución territorial dentro de la subcuenca

ÁREA (Ha)
MUNICIPIO
Urbana Rural Total
Albán 147.46 147.46
Anolaima 8.29 8.29
Bojaca 203.91 5241.53 5445.44
El Rosal 68.84 8422.96 8491.8
Facatativá 1156.01 13926.71 15082.72
Funza 633.08 633.08
La Vega 16.51 16.51

18
 Madrid 602.48 10724.81 11327.29
Mosquera 45.21 3244.94 3290.15
Pacho 65.4 65.4
San Francisco 147.96 147.96
Sasaima 75.13 75.13
Soacha 119.33 119.33
Subachoque 81.37 16840.91 16922.28
Tabio 129.56 129.56
Tenjo 31.28 31.28
Zipacon 446.07 446.07
Zipaquirá 61.87 61.87
Fuente Autores

En la Tabla 2 se presenta la información más precisa de la población que se encuentra

ubicada en la cuenca, en la Tabla 3 y Tabla 4 se presenta información característica de la
cuenca como lo es el paisaje, especies de mayor peso ecológico, y datos específicos del
suelo por el cual se encuentra rodeada la cuenca.

Tabla 2 Porcentaje de población que se ubica en la Cuenca del Rio Balsillas

POBLACION
SUBCUENCA Urbana % Rural % Total %
RIO BALSILLAS 125593 13.39 19491 5.92 145084 11.42
Fuente [¡Error! No se encuentra el origen de la referencia.]

Tabla 3 Coberturas y especies forestales en la Cuenca del Rio Balsillas

SUBCUENCA Coberturas vegetales Especies de mayor peso ecológico

Bosque Primario Hyeronima rufa

Oreopanax bogotensis
Diplostephium rosmarinifolium
Weinmannia tomentosa
Rastrojo - Matorral
RIO Myrsine guianensis
BALSILLAS Miconia squamulosa
Pernettya prostata
Hypericum thuyoides
Vegetación de Paramo Miconia salicifolia
Blechunum sp.
Paepalanthus columbiensis
Fuente [¡Error! No se encuentra el origen de la referencia.]

19
 Tabla 4 Características del suelo de la Cuenca del Rio Balsillas

SUBCUENCA RIO BALSILLAS

Plano en una condición climática fría
Paisaje mayormente seca.

pH ligeramente ácidos, moderada saturación

de aluminio, moderada a alta saturación de
bases, moderados contenidos de carbón
Características
orgánico, moderados contenidos de fosforo,
químicas
potasio y calcio. En general, los suelos son de
fertilidad moderada a alta.

Texturas moderadamente finas, estructuras

moderadas, densidades aparentes medias en
suelos originados a partir de rocas y bajas en
Características aquellos que son producto de la evolución de
físicas cenizas volcánicas y una alta susceptibilidad
a la erosión hídrica cuando son desprovistos
de vegetación.

Dominancia equilibrada de los suelos de los

Taxonomía órdenes andisol e inceptisol

Fuente [¡Error! No se encuentra el origen de la referencia.]

7.2. Topografía y morfometría

La Cuenca del Rio Balsillas cuenta con altitudes que oscilan entre los 2600 y 3500 msnm,
se caracteriza por relieves montañosos en su parte nort-oriental con predominio en relieves
ondulados fuertemente ondulados a fuertemente quebrados y escarpados constituyendo
42.5% del área total [¡Error! No se encuentra el origen de la referencia.].

En la zona centro-sur de la subcuenca se localizan la planicie lacustre aluvial a lo largo del

Rio Subachoque – Balsillas con predominio de relieves planos a ligeramente ondulados.
Las zonas de relieve plano a ligeramente inclinado constituyen el 43.2% del área total
[¡Error! No se encuentra el origen de la referencia.].
Sobre las zonas planas y ligeramente inclinadas se asientan las cabeceras municipales de
Facatativá, Madrid, El Rosal y Subachoque. Las zonas con predominio de relieve ondulado
a inclinado se asientan sobre el piedemonte limitado con la planicie aluvial y constituyen el
14.2% del área. En la Tabla 5 se presenta información de la distribución de los rangos de
pendiente para la Cuenca del Rio Balsillas [¡Error! No se encuentra el origen de la
referencia.].

Tabla 5 Distribución de los diferentes rangos de pendiente de la Cuenca del Rio

Balsillas.

GRADO Relieve Símbolo Área (ha) Porcentaje

20
 Plano, plano cóncavo y ligeramente
0-3 % a 12733.6 20.39
plano
Ligeramente inclinado, ligeramente
3-7 % b 14256.05 22.83
ondulado
7-12 % Ondulado, inclinado c 8856.67 14.18
Fuertemente ondulado, fuertemente
12-25 % d 15005.67 24.03
inclinado
25-50 % Fuertemente quebrado e 9846.25 15.77
50-75 % Escarpado f 1574.73 2.52
> 75 % Muy escarpado g 168.63 0.27
Fuente [¡Error! No se encuentra el origen de la referencia.]

En la Tabla 6 se presenta información de las características morfométricas de la cuenca, la

cual cuenta un área total de 644 km2 y una longitud de 68 km, es de tipo triangular, su
cabecera se ubica sobre la cota 3700 msnm y su cota más baja sobre los 2500 msnm,
presentando un drenaje superficial rápido de tipo sinuoso [¡Error! No se encuentra el
origen de la referencia.].

Tabla 6 Características morfométricas de la Cuenca del Rio Balsillas

CUENCA Rio Balsillas

Cota mayor (m) 3,700
Cota menor (m) 2,500
Pendiente media (m) 0.02
Altitud media 2730 msnm
Área (km2) 644
Perímetro (km) 173
Longitud cauce principal (km) 68
Ancho cuenca (m) 14,302
Fuente Autores

Mapa 2 Elevación medio de la Cuenca del Rio Balsillas

21
 Fuente Autores

7.3. Climatología

Como es bien sabido el clima es el conjunto de condiciones atmosféricas que caracterizan

el tiempo atmosférico y la evolución de una determinada región. Esto se determina bajo el
análisis espacio temporal de los elementos que lo definen como lo es la precipitación, la
temperatura, la humedad relativa, el brillo solar, entre otros, y factores que lo afectan como
las pendientes, la altitud, el relieve o cobertura vegetal.

Así mismo, el BoxPlot es un instrumento gráfico de la estadística descriptiva que permite

realizar un análisis más detallado y conciso respecto a la distribución de los datos en la
muestra. Esto se complementa, de manera cualitativa, con los resultados cuantitativos
obtenidos a través de los estadísticos de la muestra. Adicionalmente un Diagrama de Cajas
permite determinar si la muestra tiene elementos “outliers” y si presenta un sesgo a la
izquierda a la derecha o izquierda.

22
El gráfico de caja es una buena alternativa a la presentación tradicional de datos medidos
con escala cuantitativa: el histograma. En éste algunas características de la distribución
pueden quedar enmascaradas por la construcción de los intervalos de clase; como, por
ejemplo: los valores de la escala que corresponden al 50% de las observaciones centrales
y al 25% de las observaciones de cada lado, el valor de la mediana, los valores de
observaciones muy apartadas de las posiciones centrales, etc. Es decir, es una
presentación que agrega a las ventajas de los gráficos (rápida comprensión visual) algunas
referencias estadísticas Por otra parte, permite cotejar varias series de datos medidas con
la misma escala y ubicadas en posiciones parecidas de ésta; siendo, en tal sentido, más
claro y de mayor información que el polígono de frecuencias (utilizado para comparar dos
o más histogramas) [¡Error! No se encuentra el origen de la referencia.].

La Cuenca del Rio Balsillas se caracteriza por la presencia de climas fríos, húmedos y secos
debido a su topografía y sus variaciones altimétricas, estas características están
condicionadas como anteriormente se mencionaba por las variaciones de temperatura la
cual oscila generalmente entre los 10 °C y los 15 °C, teniendo como temperatura media
aproximada 13 °C, esta distribución de temperatura se debe a su vez a la elevación a la
que se encuentra la cuenca ya que su altitud varia de los 2672 m.s.n.m. a los 3860 m.s.n.m.

De acuerdo a lo anterior se puede decir que la temperatura varia inversamente con la altura,
teniendo una variación de 1 °C por cada 100 m de altura [¡Error! No se encuentra el origen
de la referencia.]. A continuación, se presentan los mapas de interpolación de las variables
climatológicas como lo son temperatura, precipitación, brillo solar y humedad relativa, las
cuales se realizaron bajo el método de interpolación de distancia inversa ponderada (IDW)
el cual fue seleccionado por su amplio uso en el campo climatológico y meteorológico
además de su precisión en los resultados arrojados con el fin de adaptarse a los
requerimientos de quien lo utilice.

En el Mapa 3 se observa que para la variable de temperatura dentro de la cuenca se

presentan los valores más bajos en la parte alta de esta, en las zonas con elevaciones bajas
se observan temperaturas de 13.60°C a 14 °C bajo las cotas 2600 m.s.n.m. a 2800 m.s.n.m.
Las partes altas de la cuenca, presentan temperaturas inferiores a los 13 °C entre las cotas
3000 m.s.n.m. y 3300 m.s.n.m aproximadamente. En la Grafica 1 Temperatura media en
°C de la Cuenca del Rio Balsillas.Grafica 1 de temperatura promedio de la cuenca, se
observa un comportamiento bimodal registrando los valores más bajos en los meses de
julio, agosto y septiembre, siendo julio el que presente la temperatura más baja con un valor
de 12.5 °C. los valores más altos se presentan en los meses de abril y mayo para el primer
semestre del año con un valor de 14 °C y 14.10 °C respectivamente y el mes de octubre
con un valor de 13.86 °C, en general se puede decir que la variación de temperatura a lo
largo del año no supera los 3°C.

23
Mapa 3 Temperatura media en °C de la Cuenca del Rio Balsillas

Fuente Autores

Grafica 1 Temperatura media en °C de la Cuenca del Rio Balsillas.

24
 TEMPERATURA PROMEDIO
14,50
14,00
13,50
PROMEDIO

°C
13,00
12,50
12,00
11,50

MESES

Fuente Autores

En el Mapa 4 se puede observar que la precipitación varia de 769.57 mm a 1039.3 mm, se

presentan los valores máximos en la parte norte de la cuenca con variaciones de 940.83
mm a 1039.3 mm, y los valores mínimos en la parte sur con variaciones de 769.57 mm a
838.07 mm. En la Grafica 2 de precipitación media se puede observar un comportamiento
bimodal registrando los valores más bajos en los meses de enero y febrero con 650 mm y
710 mm para el primer semestre del año y los meses de agosto y septiembre con 705.96
mm y 730 mm para el segundo semestre del año. El mes que presenta el valor máximo es
octubre con un registro de 1000 mm.

25
Mapa 4 Precipitación media en mm de la Cuenca del Rio Balsillas

Fuente Autores

Grafica 2 Precipitación media en mm de la Cuenca del Rio Balsillas.

26
 PRECIPITACIÓN PROMEDIO
1000
800
600

mm
400
PROMEDIO
200
0

MESES
Fuente Autores

En la Grafica 3 se observa el diagrama de cajas para la variable de precipitación que a

diferencia de la humedad relativa y el brillo solar, presenta una distribución asimétrica, es
decir, tiene un sesgo negativo donde la mediana está más cerca de la zona superior lo que
implica que el 25% de las observaciones comprendidas entre el segundo y el tercer cuartil
ocupa un menor tramo en la escala, así mismo, se observa gran variación y distancia entre
el LI (límite inferior) y el dato atípico que se presenta además de que este dato es negativo
con un valor de -133.45, este dato negativo se puede asociar a que posiblemente se
presentó un error en la toma de datos o se puede asociar a sucesos naturales presentados
como el fenómeno del niño dentro de la cuenca del Rio Balsillas.

Grafica 3 Diagrama de cajas de precipitación de la Cuenca del Rio Balsillas

DIAGRAMA DE CAJAS DE PRECIPITACIÓN

55000,00
Series3
45000,00 Series2
Series1
35000,00 33735,39
6732,69 MEDIA
25000,00 23636,35 MEDIANA
13404,82
16903,66 3° QUARTIL
15000,00
9625,90 11028,27
LI
5000,00 13537,32
6804,62 LS
-5000,00 1 DALi
DALs
-15000,00 -13393,45

Fuente Autores

27
En el Mapa 5 se presenta la distribución de humedad relativa dentro de la cuenca del Rio
Balsillas. El comportamiento de la humedad relativa es de tipo bimodal compartiendo
algunos patrones de comportamiento con la variable de la precipitación, en la Grafica 4 se
observa el registro de valores de humedad relativa más altos en los meses de abril y mayo
con 80.30% y 79.70% respectivamente, para el primer semestre del año y en los meses de
octubre y noviembre con valores de 80.10% y 81.90% respectivamente, para el segundo
semestre del año. La humedad más baja se registra en los meses de julio, agosto y
septiembre con valores de 76.70%, 75.80% y 76.70%. El porcentaje promedio anual de
humedad relativa es de 78.7% para la cuenca del Rio Balsillas.

Mapa 5 Humedad Relativa en % de la Cuenca del Rio Balsillas

Fuente Autores

28
 Grafica 4 Promedio mensual de Humedad relativa de la Cuenca del Rio Balsillas

HUMEDAD RELATIVA
83
82
81
80
79
78
%
77
76 PROMEDIO
75
74
73
72

MESES
Fuente Autores

En la Grafica 5 se observa el diagrama de cajas para la variable de humedad relativa el

cual presenta una distribución asimétrica, es decir, tiene un sesgo positivo donde la
mediana está más cerca de la zona inferior lo que implica que el 25% de las observaciones
comprendidas entre el segundo y primer cuartil ocupa un menor tramo en la escala, se
observa un que el LI (límite inferior) es un dato negativo con un valor de -23.98, al igual que
el atípico con un valor de -86.30 que se encuentra por fuera de la caja y los bigotes que la
componen, estos datos negativos se pueden asociar a que posiblemente se presentó un
error en la toma de datos o se puede asociar a sucesos naturales presentados como el
fenómeno del niño dentro de la cuenca del Rio Balsillas.

Grafica 5 Diagrama de cajas de Humedad relativa de la Cuenca del Rio Balsillas

DIAGRAMA DE CAJAS DE HUMEDAD RELATIVA

300,00
Series3
250,00
Series2
200,00 204,55
Series1
150,00 142,23 LI
100,00 LS
79,90
50,00 MEDIA
38,35
MEDIANA
0,00
-23,98 1 1° QUARTIL
-50,00
3° QUARTIL
-86,30
-100,00
Fuente Autores

29
En el Mapa 6 se presenta la distribución de brillo solar para la Cuenca del Rio Balsillas. En
la Grafica 6 se observa que esta variable presenta un comportamiento bimodal, para el
mes de enero y de diciembre se presentan valores máximos de 180.4 horas y 158.8 horas
respectivamente, para el mes de abril y octubre se presentan valores mínimos de 110.60
horas y 112.9 horas respectivamente, esta variable tiene una relación inversa con respecto
a la precipitación y a la humedad relativa.

Mapa 6 Brillo solar en horas de la Cuenca del Rio Balsillas

Fuente Autores

30
 Grafica 6 Promedio mensual de Brillo solar de la Cuenca del Rio Balsillas

BRILLO SOLAR
200
HORAS 150

100
PROMEDIO
50

MESES
Fuente Autores

En la Grafica 7 se observa el diagrama de cajas para la variable de brillo solar el cual

presenta una distribución asimétrica, es decir, tiene un sesgo positivo donde la mediana
está más cerca de la zona inferior lo que implica que el 25% de las observaciones
comprendidas entre el segundo y primer cuartil ocupa un menor tramo en la escala, se
observa un dato atípico ya que se encuentra por fuera de la caja y los bigotes que la
componen, este dato atípico podemos asociarlo a que posiblemente se presentó un error
en la toma de datos o se puede asociar a sucesos naturales presentados dentro de la
cuenca del Rio Balsillas.

Grafica 7 Diagrama de cajas de Brillo Solar de la Cuenca del Rio Balsillas

DIAGRAMA DE CAJAS DE BRILLO SOLAR

320,00

270,00 LI
LS
220,00 213,75 MEDIANA
1° QUARTIL
170,00 177,15
3° QUARTIL
140,55 133,33
120,00 MEDIA
116,15
DALi
79,55
70,00 DALs
42,95
20,00

Fuente Autores

31
 7.4. Recursos hídricos e hidrología

Para efectos de este trabajo, es necesario analizar los procesos hidrológicos que componen
el ciclo del agua. En primer lugar, debe determinarse el balance hídrico, que es una
herramienta que implementa el principio de conservación de masa en el cual el agua
almacenada será determinada por la diferencia entre los flujos de agua de entrada y salida
[¡Error! No se encuentra el origen de la referencia.].
𝛥𝑠
En la Ecuación 1, es el almacenamiento en el tiempo, las entradas consideradas son la
𝛥𝑡
precipitación, el aporte de aguas subterráneas y el aporte de aguas superficiales y las
salidas son la infiltración, la evapotranspiración y en algunos casos se consideran las
captaciones [¡Error! No se encuentra el origen de la referencia.].

Ecuación 1 Balance hídrico.

𝜟𝒔
= 𝛴 𝑒𝑛𝑡𝑟𝑎𝑑𝑎𝑠 ± 𝛴 𝑠𝑎𝑙𝑖𝑑𝑎𝑠
𝜟𝒕
Fuente [¡Error! No se encuentra el origen de la referencia.]

Por otro lado, al calcular el balance hídrico se obtiene la oferta hídrica, esta es la porción
de agua precipitada que no sufre el proceso de infiltración ni de evapotranspiración y que
escurre a los cauces de los ríos para posteriormente depositarse en un área mayor llamada
cuenca; otras interpretaciones de la oferta hídrica indican que es el agua disponible para
satisfacer la demanda [¡Error! No se encuentra el origen de la referencia.]. Sin embargo,
está afirmación no es del todo cierta pues a la oferta hídrica superficial debe restarse el
caudal ecológico, es decir la cantidad suficiente de agua para mantener el funcionamiento
y la composición de los ecosistemas hídricos y sus funciones ecológicas [¡Error! No se
encuentra el origen de la referencia.].

Por otro lado, después de asegurar el caudal ecológico, es posible acceder al recurso
hídrico para suplir los requerimientos de la demanda, es decir la cantidad de agua necesaria
para cubrir las necesidades producto de las actividades sociales y económicas del hombre
[¡Error! No se encuentra el origen de la referencia.]. Cabe denotar que el cálculo de esta
variable incluye la sumatoria de las demandas por sector. Finalmente, las relaciones entre
la oferta y la demanda pueden establecer el índice de escasez de agua superficial, como
será detallado posteriormente.

7.5. Cambio climático

Para poder abarcar el concepto de cambio climático, se debe comprender el clima como la
medida de la media estadística de las variables que componen al sistema climático durante
una serie temporal de mínimo 30 años y el estado del tiempo como “el conjunto de
fenómenos atmosféricos en un momento determinado”, por lo que se puede decir que el
clima es la estadística del estado del tiempo según la Organización meteorológica mundial
[¡Error! No se encuentra el origen de la referencia.]. Y la variabilidad climática como las
fluctuaciones de los componentes del sistema climático en periodos relativamente cortos
respecto a una media determinada por una serie climatológica. Teniendo en cuenta lo
anterior, el cambio climático se define como la alteración de la composición de la atmósfera

32
que puede ser atribuida a acciones antrópicas y que sumada a la variabilidad climática
puede conllevar a cambios en las condiciones predominantes del clima [¡Error! No se
encuentra el origen de la referencia.].

El IPCC es un organismo creado en el año 1988 en conjunto con la OMM y el PNUMA que
proporciona obras de referencia para los responsables de políticas los cuales son informes
de evaluación y documentos técnicos que proveen una herramienta científica para ser
usada en la elaboración de información y servicios climáticos que ayuden a afrontar las
consecuencias del cambio climático [¡Error! No se encuentra el origen de la referencia.].
El IPCC se constituye entonces como el referente científico y técnico en lo que a cambio
climático se refiere.

Según esa organización se ha demostrado con un 95% de seguridad que el ser humano y
sus actividades son la principal causa del calentamiento evidenciado desde hace más de
cinco décadas. Fenómeno que ha provocado efectos visibles como el calentamiento de la
atmosfera y el océano, la disminución de los volúmenes de nieve y hielo, la elevación del
nivel del mar y el aumento en las concentraciones de GEI [¡Error! No se encuentra el
origen de la referencia.].

Lo anterior, producirá cambios futuros proyectados en los fenómenos meteorológicos y

climáticos extremos, siendo muy probable para finales del siglo XXI que los días y noches
fríos sean más cálidos, de igual forma aumentará la temperatura en los días y noches
calurosos en relación con la actualidad. Además, se presentarán más episodios de olas de
calor, aumento de la intensidad y duración de la sequía, mayor intensidad de actividad de
los ciclones tropicales y aumento en la magnitud de los niveles del mar [¡Error! No se
encuentra el origen de la referencia.].

Con respecto a lo referente al recurso hídrico se pueden observar disminuciones en la

precipitación mundial y un aumento en sus intensidades, debido al cambio progresivo en la
temperatura reflejándose en un incremento año a año. Estos fenómenos representan un
gran impacto en el ciclo hidrológico y por lo tanto en los flujos de agua disponibles para
satisfacer las necesidades socio-económicas de los seres humanos [¡Error! No se
encuentra el origen de la referencia.].

7.6. Escenarios de cambio climático

En las bases físicas del quinto informe de evaluación del IPCC realizadas en el año 2013,
se proponen cuatro escenarios que tienen en cuenta las políticas encaminadas a reducir
las emisiones de GEI según los objetivos del protocolo de Kyoto. Estos escenarios son
función de las Trayectorias de Concentración Representativas más conocidas como RCP
2.6, 4.5, 6.0 y 8.5, respondiendo al forzamiento radiactivo en W/m2, de cada uno de dichos
valores. El RCP 2,6 es el escenario más optimista, con un incremento de temperatura menor
a 1,5°C para el año 2100 y el RCP 8,5 el caso pesimista, debido a un incremento en la
temperatura mayor a los 4°C para ese mismo año [¡Error! No se encuentra el origen de
la referencia.].

En la Ilustración 1 se encuentran la descripción de los 4 escenarios de cambio climático

según el instituto de hidrología, meteorología y estudios ambientales; programas de las
naciones unidas para el desarrollo; ministerio de ambiente y desarrollo sostenible;

33
departamento nacional de planeación y la cancillería [¡Error! No se encuentra el origen
de la referencia.].

Ilustración 1 Descripción de los escenarios de cambio climático

Fuente [¡Error! No se encuentra el origen de la referencia.]

Las cuatro trayectorias RCP comprenden un escenario en el que los esfuerzos de

mitigación conducen a un nivel de fortalecimiento muy bajo (RPC 2.6), dos escenarios de
estabilización (RCP 4.5 y RCP 6.0) y un escenario con un nivel muy alto de emisiones de
GEI (RCP 8.5), a continuación, se explicarán detalladamente cada uno de los escenarios
RCP [¡Error! No se encuentra el origen de la referencia.]:

 RCP 2.6: Es desarrollado por el equipo de modelado IMAGE de la Agencia de

Evaluación Ambiental de los Países Bajos. Es la vía de emisión representativa para
los escenarios en la literatura que conducen a niveles de concentración de gases
de efecto invernadero muy bajos. Es un escenario denominado "pico" debido a que
su nivel de forzamiento radiactivo alcanza primero un valor de alrededor de 3,1 W/m2
a mediados del siglo, volviendo a 2,6 W/m2 en 2100. Para alcanzar tales niveles de
forzamiento radiactivo, las emisiones de gases de efecto invernadero se reducen
sustancialmente con el tiempo [¡Error! No se encuentra el origen de la
referencia.].

 RCP 4.5: Es desarrollado por el equipo de modelado de MiniCAM en el Instituto de

Investigación de Cambio Global Conjunto del Pacífico Noroeste (JGCRI). Es un
escenario de estabilización en el que el forzamiento radiactivo total se estabiliza
antes de 2100 mediante el empleo de una gama de tecnologías y estrategias para

34
 reducir las emisiones de gases de efecto invernadero [¡Error! No se encuentra el
origen de la referencia.].

 RCP 6.0: Es desarrollado por el equipo de modelado de AIM en el Instituto Nacional

de Estudios Ambientales (NIES), Japón. Es un escenario de estabilización en el que
el forzamiento radiactivo total se estabiliza después de las 2100 sin excederse por
el empleo de una gama de tecnologías y estrategias para reducir las emisiones de
gases de efecto invernadero [¡Error! No se encuentra el origen de la referencia.].

 RCP 8.5: Es desarrollado por el equipo de modelado de MESSAGE y el Marco de

Evaluación Integrado de IIASA en el Instituto Internacional para el Análisis de
Sistemas Aplicados (IIASA), Austria. Este RCP se caracteriza por el aumento de las
emisiones de gases de efecto invernadero a lo largo del tiempo, es representativo
de los escenarios en la literatura que conducen a altos niveles de concentración de
gases de efecto invernadero [¡Error! No se encuentra el origen de la referencia.].

7.7. Modelos de recirculación global (GCM)

Los modelos GCM simulan la respuesta del clima a diferentes forzamientos para proyectar
su comportamiento futuro. En su aparición en la década de los 70s solo tenían en cuenta
las variables asociadas a la atmósfera, la superficie de la tierra y el océano. En la actualidad
tienen en cuenta esas variables y las asociadas con los aerosoles, el ciclo del carbón, la
vegetación, y el hielo terrestre. En los reportes del IPCC los modelos son sometidos a una
revisión de acuerdo con una serie de simulaciones realizada por el Proyecto de
Comparación de Modelos Acoplados (CMIP) que establece su confiabilidad en función de
que tan bien reproducen el clima observado. En la actualidad se encuentra vigente el CMIP5
que incluye más de 30 modelos producidos por investigadores de diferentes
nacionalidades, destacando Estados Unidos, Japón, China y el Reino Unido [¡Error! No se
encuentra el origen de la referencia.].

Los GCM representan el clima utilizando una cuadrícula tridimensional en todo el mundo
como se observa en Ilustración 2, generalmente tienen una resolución horizontal de entre
250 km y 600 km, de 10 a 20 capas verticales en la atmósfera y a veces hasta 30 capas en
los océanos. Su resolución es grande en relación con la escala necesaria en la mayoría de
las evaluaciones de impacto; muchos procesos físicos, como los relacionados con las
nubes ocurren a escalas más pequeñas y no pueden modelarse adecuadamente, esta es
una fuente de incertidumbre en las simulaciones del clima futuro basadas en GCM [¡Error!
No se encuentra el origen de la referencia.].

Los modelos numéricos no son siempre exactos, por lo cual hay algunos que son mejores
para determinadas zonas que para otras, lo que hace necesario que en la evaluación de
impactos usando GCM se valoré el modelo en la representación del clima actual local
[¡Error! No se encuentra el origen de la referencia.].

Ilustración 2 Evolución de la resolución espacial de los GCM a través del tiempo

35
 Fuente [¡Error! No se encuentra el origen de la referencia.]

7.8. Reducción a escala

Considerando las dimensiones de las celdas de los GCM, es necesario aplicar la

metodología de downscalling, con el fin de utilizar los datos generados por estos modelos,
en zonas donde los detalles como la topografía, las coberturas de suelo, los cuerpos de
agua, entre otros, deben ser tenidos en cuenta para la toma de decisiones. Esta
metodología es efectiva para utilizar los datos de los modelos climáticos globales de una
escala de 100 a 500 km, en una escala de 10 a 50 km, espacialmente hablando, como
observa en la Ilustración 3, de una escala anual o mensual, en una escala diaria o
subdiaria, en lo que a escala temporal se refiere [¡Error! No se encuentra el origen de la
referencia.].

Ilustración 3 Esquema del concepto downscalling espacial

36
 Fuente [¡Error! No se encuentra el origen de la referencia.]

Existen dos tipos de técnicas de reducción de escala; la dinámica y la estadística, cada uno
con sus ventajas y desventajas asociadas. Para la reducción de escala dinámica es
necesario contar con una gran cantidad de datos a escala más fina la cual se basa en un
modelo climático regional (RCM) para representar el proceso atmosférico físico, el RCM
toma la información atmosférica suministrada por el GCM en los límites laterales e incorpora
topografía más compleja, el contraste océano-superficie, heterogeneidades de superficie y
descripciones detalladas de procesos físicos para generar información climática realista a
una resolución espacial de aproximadamente 20 a 50 km [¡Error! No se encuentra el
origen de la referencia.].

Por otro lado, la reducción de escala estadística se enfoca en establecer relaciones

empíricas entre las características históricas a gran escala y las condiciones locales
observadas. Una vez que se ha determinado y validado una relación, se utilizan variables
atmosféricas futuras que son proyectadas por el GCM y se predicen las variables climáticas
locales futuras. El downscaling estadístico puede producir proyecciones climáticas
específicas al sitio de interés, mientras que RCM al estar limitadas a una resolución espacial
de 20 a 50 km no pueden generar información a este nivel de detalle [¡Error! No se
encuentra el origen de la referencia.]. En la Tabla 7 se realiza una comparación de las
ventajas y desventajas de cada tipo de técnica

Tabla 7 Comparación de ventajas y desventajas de cada tipo de técnica

REDUCCION DE ESCALA REDUCCION DE ESCALA DINAMICA

ESTADISTICA
VENTAJAS

37
 Puede proveer variables climáticas en Genera resultados basados en procesos
escala puntual a partir de los resultados de físicos consistentes.
GCM
Permite incorporar directamente las Produce información en resolución final de
observaciones locales los GCM que puede resolver procesos
atmosféricos en escalas más pequeñas
DESVENTAJAS
Requiere series históricas para calibración Calculo computacional intensivo
de larga duración y alta confiabilidad.
Necesita una alta calidad de datos para Número limitado de escenarios disponibles
calibrar
Tamaño del dominio, región climática y la Fuertemente dependiente de los límites del
temporada afectan el nivel de la reducción GCM
de escala.
Dependiente de los límites del GCM Altamente demandante de recursos
informáticos
Fuente [¡Error! No se encuentra el origen de la referencia.]

7.8.1. Reducción a escala basada en la teoría del caos

El método de reducción de escala basado en la teoría del caos evalúa la presencia de caos
determinístico para diferentes periodos de acumulación, estos sistemas se entienden como
aquellos que cuentan con una predictibilidad a corto plazo ya que son especialmente
sensibles a las variaciones en las condiciones iniciales. El sistema predictivo caótico se
construye con los resultados del tiempo de retardo, la dimensión de embebimiento y los
exponentes de Lyapunov encontrados para los sistemas en el intervalo de acumulación de
precipitación óptimo de cada sistema dinámico. El modelo predictivo se basa en la
sincronización entre las estaciones climáticas y el modelo climático global, dadas por el
parámetro μ del método de los vecinos falsos mutuos más próximos [¡Error! No se
encuentra el origen de la referencia.], [¡Error! No se encuentra el origen de la
referencia.].

En primer lugar se lleva a cabo evaluación de la presencia del caos determinístico para
obtener el menor periodo de acumulación en el cual se presenta caos; es necesario
reconstruir el espacio de fase por medio del método del Time Delay para calcular los valores
del tiempo de retardo (t) y la dimensión de embebimiento (m) necesarios para obtener los
exponentes de Lyapunov [¡Error! No se encuentra el origen de la referencia.], [¡Error!
No se encuentra el origen de la referencia.], [¡Error! No se encuentra el origen de la
referencia.], [¡Error! No se encuentra el origen de la referencia.], [¡Error! No se
encuentra el origen de la referencia.].

El espacio de fase es esencialmente un gráfico, cuyas coordenadas representan las variables

necesarias para describir completamente el estado del sistema en cualquier momento, la trayectoria
del diagrama de espacio de fase describe la evolución del sistema para un estado inicial. La "región
de atracción" de estas trayectorias en el espacio de fase proporciona información cualitativa
importante para determinar el grado de complejidad del sistema ", como se muestra en la figura 8
donde el espacio de fase se observa en tres dimensiones para diferentes intervalos de acumulación
de precipitación. El conjunto limitante que reúne trayectorias asintóticas cercanas al equilibrio se
conoce como "atractor". Los atractores de sistemas caóticos deterministas pueden exhibir un tipo

38
inusual de auto-similitud y estructuras presentes en todas sus escalas y, por lo tanto, es necesario
encontrar una dimensión apropiada del plano de fase, de modo que la estructura del atractor
permanezca invariante [¡Error! No se encuentra el origen de la referencia.], [¡Error! No se
encuentra el origen de la referencia.], [¡Error! No se encuentra el origen de la
referencia.].

Ilustración 4 Comparación del espacio de fase de la serie de tiempo de precipitación

Fuente [¡Error! No se encuentra el origen de la referencia.]

Posteriormente, se realiza la sincronización de los sistemas por medio del método del
vecino falso más cercano que es una técnica estadística basada en el cálculo del parámetro
μ como se observa en la Ecuación 2, que evalúa los vecindarios locales entre dos series
de tiempo, de modo que μ toma valores del orden de 1 si existe una sincronización general
completa , de lo contrario μ tiene que ser un número cuya magnitud sea comparable con el
producto del tamaño del atractor dividido por el producto de la distancia entre los vecinos
más cercanos en la serie de tiempo [¡Error! No se encuentra el origen de la referencia.].

Ecuación 2 Cálculo de μ
𝑁
1 ‖𝑦𝑛 − 𝑦𝑛𝑁𝑁𝐷 ‖ ‖𝑥𝑛 − 𝑥𝑛𝑁𝑁𝑅 ‖
𝜇= ∑
𝑁 ‖𝑥𝑛𝑁𝑁𝐷 ‖ ‖𝑦𝑛 − 𝑦𝑛𝑁𝑁𝑅 ‖
𝑛=1
Fuente [¡Error! No se encuentra el origen de la referencia.]

Donde:
 xn y y n , son puntos de la unidad y el sistema de respuesta en un instante n
 xn NND y yn NNR son los vecinos más cercanos de xn y y n en sus respectivos sistemas
 xnNNR y y nNND son puntos del sistema opuesto en el tiempo n para ynnR y xnnD

Ilustración 5 Relación entre los puntos de los atractores

39
 Fuente [¡Error! No se encuentra el origen de la referencia.]

7.9. Experimento regional coordinado de reducción climática (CORDEX)

CODEX es un proyecto mundial promovido por el WCRP (World Climate Research

Programme). El WCRP está conformado por miembros de la Organización Meteorológica
Mundial, del Consejo Internacional para la Ciencia y de la Comisión Oceanográfica
Intergubernamental de la UNESCO.

La misión de CORDEX, es producir un conjunto de escenarios regionales de cambio

climático a nivel global, para contribuir con el quinto reporte de evaluación del IPCC y a la
vez, crear un espacio para que la comunidad mundial interesada en la realización de
escenarios regionales de cambio climático comparta información, resultados, experiencias,
tecnología y encuentre soporte técnico. CORDEX busca mejorar la interacción y la
comunicación entre modeladores del clima global, los que están calculando escenarios
regionales y los que están interesados en realizar estudios de vulnerabilidad, adaptación y
mitigación debido al cambio climático [¡Error! No se encuentra el origen de la
referencia.].

La mayoría de la información que existe se ha basado en el uso de modelos de circulación

general océano- atmosfera (AOGCMs) que resulta ser una herramienta muy valiosa para
comprender los procesos que determinan la respuesta del clima frente a los procesos fuerza
antropogénica tales como los aumentos de gases efecto invernadero (GEI),
concentraciones atmosféricas y los cambios en el uso del suelo [¡Error! No se encuentra
el origen de la referencia.], [¡Error! No se encuentra el origen de la referencia.], [¡Error!
No se encuentra el origen de la referencia.]. En los últimos años se han tenido grandes
avances y mejoras en los modelos con respecto al comportamiento de la atmosfera y de
los procesos de la superficie de la tierra además de un avance tecnológico; el modelo
AOGCMs aun presenta un desfase de cientos de km en la precisión de los resultados ya
que tiene resoluciones más gruesas lo cual dificulta la obtención de la información de los
efectos que ocurren a nivel local, impidiendo la descripción precisa de acontecimientos
extremos que son fundamentales para entender la variabilidad y el cambio climático a nivel
regional y local por lo que existe una gran brecha entre el comportamiento del clima y la

40
información proporcionada por este modelo [¡Error! No se encuentra el origen de la
referencia.].

Con respecto a la problemática de “regionalización" o "reducción de escala" se han

desarrollado varias técnicas para refinar los resultados obtenidos en los modelos GCM en
especial en el AOGCMs para reducir o eliminar la brecha que existe entre la información
del clima y escala espacial. El downscaling dinámico (DD) es una de las técnicas
implementadas la cual hace uso de modelos basados en variables físicas como alta
resolución y variables resolución, así mismo, el downscaling estadístico (SD) es otra técnica
en donde las relaciones estadísticas son las que desarrollan y predicen a gran escala, a
escala regional y local. Por medio de estas herramientas se diseñó el Experimento Regional
Coordinado de Reducción del Clima (CORDEX), el cual tiene como objetivo proporcionar
información accesible de reducción de escala y del clima para ser aplicada por la ciencia ya
que hace uso máximo de sus dominios que cubren toda África, Australia, Sur América,
Norte América y Europa, por lo tanto abarcar la mayor parte de las zonas del mundo; de
esta forma logra solucionar los problemas de regionalización que se presentan con otros
modelos brindando información más precisa a nivel local y regional para entender mejor el
comportamiento del cambio climático y lograr formular mejores medidas [¡Error! No se
encuentra el origen de la referencia.].

7.10. Sistema de indicadores hídricos

El índice de escasez ha sido reemplazado por el índice de Uso de agua superficial como
un indicador de la intervención antrópica, consignado en el Estudio Nacional del Agua
(IDEAM, 2015). Estos índices, miden el estado del recurso hídrico ya sea cualitativa o
cuantitativamente, para el proceso de toma de decisiones. La evaluación nacional del agua
tiene soporte en conceptos y productos temáticos y en un análisis integrado de cambios
con respecto a la situación de referencia, que se sintetizan en gran medida en el conjunto
de indicadores hídricos. El sistema de indicadores está integrado por dos grandes grupos:
los que dan cuenta del régimen del sistema hídrico natural y los indicadores asociados con
la intervención antrópica [¡Error! No se encuentra el origen de la referencia.].

Los indicadores del sistema hídrico se componen de:

 El índice de retención y regulación hídrica (IRH)
 El índice de aridez (IA)

Los indicadores de intervención antrópica se componen de:

 Indicadores de presión por uso de agua que abarca los índices de uso del agua
superficial (IUA), índice de presión hídrica al ecosistema (IPHE), índice de agua no
retornada a la cuenca (IARC) y el índice de eficiencia en el uso del agua (IEUA).
 Indicadores de estado de calidad y presión por contaminación que abarca el índice
de calidad de agua (ICA) y el índice de alteración potencial a la calidad del agua
(IACAL)
 Indicador de vulnerabilidad que abarca el índice de vulnerabilidad por
desabastecimiento (IVH)

El grupo de indicadores de intervención antrópica a su vez se aborda desde tres temas

fundamentales en la evaluación: la presión por uso de agua, el estado de la calidad hídrica

41
y la presión por contaminación sobre las condiciones de calidad, y, por último, la
vulnerabilidad al desabastecimiento [¡Error! No se encuentra el origen de la referencia.].

7.10.1. Índice de retención y regulación hídrica (IRH)

El Índice de Retención y Regulación Hídrica (IRH), mide la cantidad de humedad que

pueden retener las cuencas. En las áreas del Pacífico, Amazonía y Bajo Cauca, se presenta
alta capacidad para retener humedad y mantener condiciones de regulación, mientras que
son evidentes las condiciones de baja regulación en sectores de las áreas del Magdalena-
Cauca, Orinoco y Caribe [¡Error! No se encuentra el origen de la referencia.].

Mide la capacidad de retención de humedad de las cuencas con base en la distribución de

las series de frecuencias acumuladas de los caudales diarios. Evalúa la capacidad de la
cuenca para mantener un régimen de caudales, producto de la interacción del sistema
suelo-vegetación con las condiciones climáticas y con las características físicas y morfo
métricas de la cuenca; este indicador permite evaluar la capacidad de regulación del
sistema en su conjunto. Este índice se calcula con base en la curva de duración de caudales
medios diarios, curva de frecuencias que permite reconocer las condiciones de regulación
de la cuenca y los valores característicos de caudales medios e interpretar, de manera
general, las características del régimen hidrológico de un río y el comportamiento de la
retención y la regulación de la humedad en la cuenca [¡Error! No se encuentra el origen
de la referencia.].

7.10.2. Índice de uso del agua superficial (IUA)

El índice del Uso del Agua (IUA), corresponde a la cantidad de agua utilizada por los
diferentes sectores usuarios, en un período determinado (anual, mensual) y por unidad
espacial de análisis ya sean áreas, subzonas hidrográficas y cuencas abastecedoras de
acueductos en relación con la oferta hídrica superficial disponible para las mismas unidades
de tiempo y espaciales [¡Error! No se encuentra el origen de la referencia.].

Para este índice se tiene en cuenta la categorización de recurso hídrico por presiones de
demanda y variabilidad de la oferta hídrica, se emplea como insumo la categorización de
variabilidad de la oferta hídrica natural y los usos de agua [¡Error! No se encuentra el
origen de la referencia.].

7.10.3. Índice de vulnerabilidad por desabastecimiento (IVH)

Este indicador, mide el grado de fragilidad del sistema hídrico para mantener una oferta en
el abastecimiento de agua, que ante amenazas –como periodos largos de estiaje o eventos
como el fenómeno cálido del Pacifico (El Nino) podría generar riesgos de
desabastecimiento [¡Error! No se encuentra el origen de la referencia.].

7.11. Modelación hidrológica

En primer lugar, hay que destacar que un modelo es una representación de la realidad de
forma simplificada que puede ser físico o matemático y que sirve para apoyar la toma de
decisiones en diversos contextos [¡Error! No se encuentra el origen de la referencia.].

42
Específicamente, en un modelo hidrológico se representa el sistema físico y real
considerado como la cuenca de estudio, que es caracterizada por los componentes del ciclo
hidrológico. Sus entradas y salidas dependen de la finalidad con la que se haga el modelo,
pero en general sus resultados deben ser representativos y coherentes con las condiciones
simuladas que dependen de la precipitación o de los caudales provenientes de la
escorrentía superficial [¡Error! No se encuentra el origen de la referencia.].

Para la modelación hidrológica, existen diversas herramientas que se usan a nivel global y
local según sea el caso y en condición de los datos que se tengan disponibles y que el
modelo a usar requiera como información de entrada, se implementan usualmente los
modelos WEAP, SWAT y TETIS que han sido ampliamente desarrollados y aplicados
[¡Error! No se encuentra el origen de la referencia.].

Unos de los software más utilizado en Colombia son QGIS, SAGA GIS, HEC HMS,
MODFLOW, PYTHON; según estudios realizados estos software tienen un margen de error
amplio debido a que las variables que los componen no son suficientes a la hora de realizar
la modelación por ende, PepsiCo y el Banco interamericano de Desarrollo han anunciado
que a partir del presente año se espera implementar el software Hydro-BID puesto que, es
una herramienta con mayor ventaja en relación a los anteriores ya que tiene más variables
como tipo de suelo, uso de la tierra, temperatura, lluvia, caudales, y arroja resultados
concisos para determinar un análisis pertinente de riesgo y diseño de infraestructura
adaptativa.

7.11.1. Hydro-BID

Hydro-BID es un software diseñado por la fundación de PepsiCo y el Banco Interamericano

de Desarrollo (BID), el cual fue creado en América Latina y el Caribe (LAC), con el fin de
facilitar la modelación de aguas en áreas donde esta sea escasa, este programa se
encuentra implementado en Perú desde el 2014 aplicándose al 65% de las cuencas de
agua de este país; luego llego en 2017 a Guatemala, Ecuador y Argentina. Antes del 2025
se cree que Hydro-BID será implementado en al menos tres países de LAC con el objetivo
de beneficiar 725.000 personas aproximadamente [¡Error! No se encuentra el origen de
la referencia.].

Es un programa que permite evaluar el potencial de cambio futuro en la disponibilidad de

agua. Este software, cuenta con un apoyo de sistemas de información geográfica que
permite la navegación aguas arriba y aguas abajo de las cuencas a estudiar por medio de
la herramienta QGIS acompañado de la posibilidad de seleccionar el tiempo a simular,
además de trabajar con la base de datos de clima e hidrología analítica diseñada para
América Latina y el Caribe, conocida como AHD-LAC. También, ofrece una interfaz en
donde se aplican las entradas de datos climatológicos para la obtención de datos de
precipitación y temperatura para la zona y periodo de interés con la posibilidad de ingresar
la variación en los escenarios de cambio climático, acompañado de los factores que
dependen del tipo de suelo como el porcentaje de impermeabilización, el número de curva
y el porcentaje de agua que se percola a los cuerpos de agua subterránea [¡Error! No se
encuentra el origen de la referencia.].

43
El procesamiento anterior arroja resultados con estructuras de serie de datos diarias con
los flujos estimados para la zona y tiempo de estudio, lo que facilita la interpretación de los
datos, dándole las herramientas necesarias a las agencias de planificación y manejo del
recurso hídrico para una toma de decisiones más acertadas que reduzcan la vulnerabilidad
de las poblaciones aledañas [¡Error! No se encuentra el origen de la referencia.]. Como
se observa en la Ilustración 6.

Ilustración 6 Sistema Hydro-BID y sus componentes

Fuente [¡Error! No se encuentra el origen de la referencia.]

7.11.1.1. Calibración del modelo con el software Hydro-BID

Hydro-BID incluye una serie de parámetros para realizar la calibración de cualquier modelo
los cuales, permiten ajustar la serie de datos observados, equivalentes a los caudales
obtenidos de las estaciones limnimétricas, las cuales son asociadas con las series
resultantes de la simulación de esta forma, se garantiza el éxito de un modelo confiable.
Este proceso tiene unos parámetros de calibración sugeridos, según el modelo o los datos
a ingresar al software, la calibración se realiza por el método de ensayo y error ingresando
los datos de forma manual. En la Tabla 8 se describen los parámetros requeridos y los
valores sugeridos por el software, en la Tabla 9 se expresan los valores óptimos
correspondientes a las métricas calculadas por el software.

Tabla 8 Parámetros de calibración del software Hydro-BID

Métrica de
Descripción Valor optimo
desempeño

44
 Caracteriza el tipo de suelo de la cuenca y su
Número de curva hidrología; oscila entre 0 y 1, en donde un valor de 1
-
(NC) representa que toda el agua que llueve se convierte en
escorrentía y 0 representa una completa infiltración.
Representa la cantidad de agua que puede ser
almacenada en el suelo para ser usada por las
Contenido de
plantas, afectando la infiltración hacia el agua
agua en el suelo 10
subterránea. En la interfaz puede ingresarse como un
(AWC)
multiplicador que altera el valor anidado en la AHD
puede ingresarse de forma directa.
Evidencia cómo el agua subterránea que se encuentra
cerca de la superficie contribuye a los caudales de los
ríos o corrientes superficiales después de un evento
Coeficiente de
de lluvia. En la interfaz puede ingresarse como un 0.01
Recesión (R)
multiplicador que altera el valor anidado en la AHD
puede ingresarse de forma directa como resultado de
una medición en campo.
Representa el intercambio entre las aguas
subterráneas cerca de la superficie con los cuerpos de
Seepage o
agua más profundos, en donde la cantidad de agua
coeficiente de 0.005
que se le atribuye a este parámetro se considera como
permeabilidad
una salida de la modelación. También puede ser
ingresado como un multiplicador o en un valor directo.
Factor ET de la Caracteriza el factor de evapotranspiración de las
estación de plantas durante la etapa de cultivo. Se puede ingresar 1
cultivo como multiplicador o como un valor directo.
Corresponde al factor de evapotranspiración en la
etapa de latencia o en la etapa en la que no se
Factor ET de la
encuentra cultivada la zona. En la interfaz puede
estación de 1
ingresarse como un multiplicador que altera el valor
latencia
anidado en la AHD puede ingresarse de forma directa
como resultado de una medición en campo.
Porcentaje de la Porcentaje estimado de la porción impermeable de la
cobertura cuenca y está disponible en algunas de las versiones 2%
impermeable del programa.
Fuente [¡Error! No se encuentra el origen de la referencia.]

Tabla 9 Métricas de desempeño para la calibración del software Hydro-BID

Parámetro Valor sugerido

Overall volumen error, ove (%) Más cercano a 0

Correlation, r (-) Cercano a 1

Modified correlation, Rmod Cercano a 1

45
 < 0.2 insuficiente
0.2 – 0.4 satisfecho
Nash-Sutcliffe Efficiency, R^2 0.4 – 0.6 bueno
0.6 – 0.8 muy bueno
> 0.8 excelente
Fuente [¡Error! No se encuentra el origen de la referencia.]

8. METODOLOGÍA

La metodología y desarrollo de este proyecto se divide en tres etapas las cuales tienen
como finalidad dar cumplimiento a los objetivos planteados inicialmente y la obtención de
información apta para recomendaciones pertinentes que sean útiles para la comunidad.

Ilustración 7 Esquema metodológico

Fuente Autores

A continuación, se describen detalladamente cada una de las etapas:

46
Primera etapa: Fase preliminar

a. Se inicia con la definición de la cuenca en la que se va a desarrollar el proyecto

teniendo en cuenta la disponibilidad de información, en este caso se trata de la
Cuenca del Rio Balsillas.
b. Se identifican las estaciones climatológicas activas ubicadas dentro de la cuenca o
en una zona de influencia de 15 metros.
c. Al tener definidas las estaciones climatológicas con las que se van a trabajar se
hace el llenado de datos faltantes con ayuda de un código generado en el software
Python el cual nos arroja las series diarias para cada una de ellas

Segunda etapa: Fase de ejecución

a. Se realiza el procesamiento de los datos de las estaciones para obtener los

centroides requeridos por el software Hydro-BID y así, asignar la información de
precipitación, caudal y temperatura correspondientes a la calibración del modelo.
b. Seguido a esto, se ingresan los datos al software para que este nos genere la serie
de caudales proyectada hasta el año 2100.
c. Posteriormente, se definen los escenarios de cambio climático que se van a aplicar,
estableciendo las variaciones de precipitación y temperatura con la que se realizará
la modelación.
d. Luego de realizar estas acciones, se procederá a validar el modelo con el fin de
corroborar que el procesamiento realizado es exitoso y coherente con las
condiciones de la cuenca.

Tercera etapa: Fase de resultados

a. Al obtener la demanda actual y futura de la cuenca se procede a realizar la

proyección para los diferentes sectores económicos en el mismo periodo de tiempo
por municipio.
b. Ya habiendo realizado este paso, se comparan los resultados de oferta y demanda
con ayuda del sistema de indicadores hídricos establecidos en el Estudio Nacional
del Agua (ENA), donde se contemplan los indicadores que hacen parte del régimen
hídrico natural y los asociados con la intervención antrópica.
c. Finalmente, se plantean las conclusiones y recomendaciones pertinentes a la
situación de demanda hídrica de la Cuenca del Rio Balsillas para un futuro.

9. RESULTADOS

9.1. Estimación de la oferta

9.1.1. Definición y procesamiento de la cuenca

La definición de la cuenca se realizó por medio de la base de datos AHD del software Hydro-
BID, adicional a esto la delimitación de la esta se determinó mediante el sistema de
información geográfica QGIS y el complemento navegador AHD el cual, arroja las
subcuencas presentes en la cuenca y su respectivo COMID (numero identificador) como se
observa en el Mapa 7, además de los ríos los ríos principales, los datos geométricos como

47
área, la longitud del cauce principal y el perímetro los cuales fueron de suma importancia
para hallar los centroides de cada subcuenca.

Mapa 7 Delimitación de la Cuenca del Rio Balsillas con el software QGIS

Fuente Autores

9.2. Obtención y tratamiento de los datos de entrada

Los datos de entrada necesarios para la modelación se encuentran en la Tabla 10 en la

cual se describe la temporalidad y la escala de estos, así como la fuente de la cual fueron
tomados.

48
 Tabla 10 Datos de entrada del modelo en el software Hydro-BID.

Datos de Número de
Fuente Descripción Temporalidad
entrada estaciones
Precipitación Estaciones
Datos obtenidos del
acumulada climáticas 26
IDEAM.
diaria en cm locales
Temperatura Estaciones
Datos obtenidos del 01/01/1983
media diaria climáticas 14
IDEAM. hasta
en °C locales
31/12/2015
Registros del servicio
Uso de la Base de datos geológico de los
NA
tierra AHD Estados Unidos
(USGS).
Fuente: Autores

A cada una de las estaciones de precipitación y temperatura se les realizó un análisis

estadístico, esto con el fin de determinar la calidad de la información y establecer el
porcentaje de datos faltantes. Mediante un código de interpolación del programa Phyton se
obtuvieron los datos faltantes con el fin de completar las series de datos desde el 01 de
enero de 1983 hasta el 31 de diciembre de 2015. A continuación, se presentan los mapas
con la ubicación de las estaciones climatológicas ubicadas dentro de la cuenca del Rio
Balsillas o en un área de influencia de 10 km, estas estaciones cuentan con un registro de
datos de más de 25 años y un porcentaje mínimo de datos faltantes.

Mapa 8 Ubicación de las estaciones de temperatura dentro de la Cuenca del Rio

Balsillas

49
 Fuente Autores

Mapa 9 Ubicación de las estaciones de precipitación dentro de la Cuenca del

Rio Balsillas

50
 Fuente Autores

Mapa 10 Ubicación de las estaciones de caudal dentro de la Cuenca del Rio

Balsillas

Fuente Autores

9.3. Calibración y validación con el software Hydro-BID

51
Para la calibración de la cuenca del Rio Balsillas se utilizó la estación 2120800
correspondiente a Puente Manrique la cual contiene datos de caudal desde 1983 hasta
2015. Para la calibración en el software Hydro-BID, se tiene influencia de ciertos parámetros
del modelo de caudal para ser más exactos en el resultado, como primera medida, se
evidenció una relación entre el Seepage y el volumen de caudal ya que estos son
inversamente proporcionales, es decir, al reducir el coeficiente de permeabilidad, el
volumen de los caudales reflejados en la curva de duración aumentan pues, se presenta
menor infiltración de agua superficial a los acuíferos; así mismo, se estableció una relación
directamente proporcional entre el coeficiente R y la pendiente de la curva y de igual, forma
para el valor del AWC y la pendiente de curva ya que estos parámetros tienen el mismo
comportamiento. Adicionalmente, cuando el valor del número de curva se reduce los
caudales disminuyen, siendo menos sensibles a los eventos de precipitación. En la Tabla
11, se encuentran los valores de los parámetros de calibración del software Hydro-BID para
la cuenca del Rio Balsillas.

Tabla 11 Parámetros de calibración para el modelo en Hydro-BID.

PARÁMETROS DE CALIBRACIÓN
COMID 301413300
ESTACIÓN 2120800
PARAMETROS VALORES
NC 0.15
AWC 1.5
R Coefficient 0.0047
Seepage 0.032
Grow Season ET Factor 1
Dormant Season ET Factor 1
Impervious Cover Percent 1
Fuente: Autores

Seguido a ingresar los parámetros de calibración al software y correr el modelo se obtienen

las métricas de desempeño a nivel diario y mensual para la cuenca, a continuación, en la
Tabla 12 se evidencian los valores arrojados por el software a nivel diario y mensual.

Tabla 12 Métricas de desempeño para la calibración del modelo en Hydro-BID.

MÉTRICAS DE DESEMPEÑO
COMID 301413300
ESTACIÓN 2120800
VALORES VALORES
PARAMETROS
DIARIOS MENSUALES
Volumen de Error
3.77 1.66
General (%)
Correlacion r(-) 0.24 0.3
Correlacion
0.13 0.18
modificada

52
 Eficiencia de Nash-
0.14[DIA1] 0.2
Sutcliffe R2
Fuente: Autores

Teniendo en cuenta los datos obtenidos, se observan mejores resultados a nivel mensual,
los cuales resultan ser más aceptables para visualizar el comportamiento de las variables
en la cuenca. Con la obtención de datos mensuales, menor será la probabilidad de error de
ajuste en las métricas debido a que existe un porcentaje menor de datos faltantes, adicional
a esto es importante mencionar que el llenado de datos diarios se realizó mediante el
programa Python el cual demuestra un margen de error más alto en el llenado debido a que
los datos faltantes son mayores que los mensuales, como se observa en la métrica de la
Eficiencia de Nash donde se expone un mejor ajuste para los valores mensuales que en los
valores diarios

Así mismo, se obtuvo una curva de duración de caudales seguido al proceso de calibración
para la cuenca, en la Grafica 8 se observa una línea de color roja la cual describe el caudal
modelado y una línea de color azul la cual describe el caudal observado; allí se evidencia
que a medida que se realiza la calibración, la curva de duración simulada se ajusta mejor a
las condiciones observadas, además de que los volúmenes de caudal se conservan,
validando el modelo.

Grafica 8 Curva de duración de caudales de la calibración del modelo

Fuente: Software Hydro-BID.

En la Ilustración 8, [DIA2]se observa con más detalle el comportamiento de los caudales

observados (línea morada) y los caudales modelados (línea azul oscura), en semejanza

53
con datos climáticos como lo es en este caso datos de precipitación (línea azul clara), esta
comparación se realiza con el fin de evidenciar la relación que tienen las variables ya que
a medida que la precipitación aumenta de forma directamente proporcional lo harán los
caudales. Para este tipo de gráfica, se debe tener en cuenta que la mayor precisión la
tienen los caudales bajos y medios debido a que los caudales observados y simulados se
relacionan más fácilmente en estos puntos y no en los picos altos ya que los eventos
extremos suceden de manera esporádica dejando claro que el comportamiento de los
caudales varia en épocas de sequía y de inundación.

Ilustración 8 Datos observados de precipitación Vs caudales observados y

simulados

Fuente Software Hydro-BID.

Seguido a esto se realiza el proceso de validación, al tener los datos de calibración se corre
el modelo nuevamente, esta vez el periodo de tiempo será el proyectado que corresponde
del 01 de enero de 2016 al 31 de diciembre de 2100, en la Grafica 9 se observa un
comportamiento de curva similar para los datos simulados y observados, cabe resaltar que
al igual que en el proceso de calibración se sobrestima en una porción los caudales
modelados.

54
 Grafica 9 Curva de duración de caudales de la validación del modelo

Fuente: Software Hydro-BID.

Con respecto a las métricas de desempeño del proceso de calibración del modelo, se
observa en la Tabla 13 que en la validación se obtiene un porcentaje de error más bajo y
una mejora en la eficiencia del número de Nash dando, así como resultado una variabilidad
de observación de simulación satisfactoria con respecto a los valores mensuales, lo que
quiere decir que este proceso es óptimo para fines que involucren valores medios
mensuales o anuales logrando así una mejor confiabilidad en sus resultados.

Tabla 13 Métricas de desempeño para la validación del modelo en Hydro-BID.

METRICAS DE DESEMPEÑO
COMID 301413300
ESTACIÓN 2120800
VALORES VALORES
PARAMETROS
DIARIOS MENSUALES
Volumen de Error
1.83 -0.24
General
Correlación r(-) 0.24 0.3
Correlación modificada 0.13 0.18
Eficiencia de Nash-
0.16[DIA3] 0.25
Sutcliffe
Fuente: Autores

55
Al comparar las métricas de desempeño para la validación y calibración del software Hidro-
BID se concluye que el comportamiento de los parámetros es similar, es decir, al obtener
datos mensuales el porcentaje de datos faltantes será menor por esta razón se evidencia
un mejor ajuste como en el caso de la eficiencia de Nash.

9.4. Análisis cambio climático

9.4.1. Modelo de circulación global (GCM)

Los escenarios de cambio climático tienen gran influencia sobre el comportamiento de la

oferta de la cuenca a futuro, para ello se selecciona un modelo de circulación global (GCM)
que se ajuste a los datos climáticos registrados en la zona, para posteriormente aplicar los
diferentes escenarios. Por medio un código en el software Spyder se hace la selección de
dicho modelo el cual, permite hacer un contraste entre diferentes criterios, para lograr una
semejanza entre los datos globales con las series de precipitación y de temperatura
generadas.

Con respecto al código de Spyder, se tienen en cuenta las series de datos de las variables
de precipitación y temperatura que se encuentran registradas en la cuenca, las
coordenadas de los centroides de las subcuencas y los datos históricos del modelo
climático global que se va a evaluar para la precipitación y la temperatura, los cuales se
obtienen de la página del Centro de Análisis de Datos Ambientales (CEDA). Inicialmente,
el código generado en este software lee los datos que se ingresan con el fin de identificar
la celda donde se encuentra ubicada la cuenca, básicamente son las coordenadas globales
donde se ubica la zona de estudio; seguido a esto se establecen los criterios a evaluar que
en este caso son: media, dato máximo, dato mínimo, desviación estándar y coeficiente de
correlación.

Con ayuda del software Spyder, se hace una evaluación de cada criterio calculando el
porcentaje de error que se obtiene al comparar los datos globales de los modelos y los
CORDEX con los datos de las series de las variables de precipitación y temperatura
históricos, así mismo, se aplican criterios estadísticos para realizar esta comparación. Se
estipularon dos criterios físicos que determinan los días secos/cálidos y los días
húmedos/fríos, debido a la relevancia que estos tienen frente al comportamiento del caudal
dentro de la cuenca; para la estructuración de dichos parámetros se realiza un algoritmo el
cual permite entender de forma más clara los procesos realizados y los resultados que
posteriormente se obtendrán.

Finalmente, para seleccionar el GCM se tienen en cuenta cuatro modelos globales y dos
experimentos CORDEX, estos últimos se caracterizan por ser de una menor escala, los
modelos que se evaluaran son: CanESM2, CNMR-MR5, MPI-ESM, CCSM4, CORDEX
CanESM2 y el CORDEX MPI-ESM. Para esto se genera un archivo Excel para realizar la
comparación de los porcentajes de error de las variables de media, dato máximo, dato
mínimo, desviación estándar, coeficiente de correlación, días húmedos/fríos y días
cálidos/secos para cada modelo y experimento CORDEX. Al tener la comparación de esos
criterios, se selecciona el modelo que mejor se ajusta y se evalúan las gráficas a nivel anual
y mensual multianual, en donde se observará la comparación entre los datos del modelo y
los datos de las series de las variables de temperatura y precipitación históricas.

56
Teniendo en cuenta lo anterior, el modelo que reflejó mejor ajuste para la variable de
precipitación es el MPI-ESM-MR CORDEX, en la Tabla 14 se pueden observar los valores
resultantes del código de comparación.

Tabla 14. Resultados criterios de comparación MPI-ESMMR –CORDEX

Modelo MPI (CORDEX)

Estacion 301388800 301388900 301143300
Media 263.99 271.05 392.15
Max 258.65 2.75.89 404.13
Min inf inf inf
Desviacion 349.31 363.5 423.76
Correlacion 0.03 0.03 0.04
Dias
0.16 0.48 0
Humedos/Frios
Dias
0.96 0.32 3.42
Secos/Calidos
Fuente Autores

Este modelo se eligió, debido a que en las gráficas de comportamiento mensual multianual
en donde se compara el GCM y los datos históricos de la variable de precipitación para la
estación 301388900 y 301413300, se presenta un comportamiento similar de la curva de
datos observados vs los datos del modelo; además, los valores que presenta el GCM son
más cercanos a los datos de la precipitación registrados, por otro lado, presenta la mejor
correlación de todos los modelos evaluados; en relación a los criterios de días secos y
húmedos el porcentaje de error es mínimo lo que indica que hay similitud entre los datos.

Grafica 10 Comportamiento mensual multianual GCM vs datos históricos de

precipitación de la estación 301388800

57
 Fuente Software Spyder

Grafica 11 Comportamiento mensual multianual GCM vs datos históricos de

precipitación de la estación 301388900

Fuente Software Spyder

Grafica 12 Comportamiento mensual multianual GCM vs datos históricos de

precipitación de la estación 301413300

58
 Fuente Software Spyder

Así mismo, para la variable de temperatura el modelo que mejor se ajusto es el GCM MPI-
ESM-MR, en la Tabla 15 se pueden observar los valores resultantes del código de
comparación.

Tabla 15. Resultados criterios de comparación MPI-ESM-MR

Modelo MPI
Estacion 301388800 301388900 301143300
Media 39.49 48 35.43
Max 42.93 49.67 40.22
Min 42.92 55.13 38.71
Desviacion 105.19 101.52 101.86
Correlacion 0.04 0.05 0.03
Dias
1.74 3.72 7.17
Humedos/Frios
Dias
0.64 0.16 0
Secos/Calidos
Fuente Autores

Este modelo se eligió, debido a que en las gráficas de comportamiento mensual multianual
en donde se compara el GCM y los datos históricos de la variable de temperatura para la
estación 301388900 y 301413300, se presenta un comportamiento similar de la curva de
datos observados vs los datos del modelo.

Grafica 13 Comportamiento mensual multianual GCM vs datos históricos de

temperatura de la estación 301388800

59
 Fuente Software Spyder

Grafica 14 Comportamiento mensual multianual GCM vs datos históricos de

temperatura de la estación 301388900

Fuente Software Spyder

Grafica 15 Comportamiento mensual multianual GCM vs datos históricos de

temperatura de la estación 301413300

60
 Fuente Software Spyder

A partir de la selección del GCM, se procede a establecer los escenarios de cambio

climático, con el fin de analizar el comportamiento del caudal en la zona de estudio para un
periodo de tiempo futuro, estos escenarios han sido estipulados en el quinto informe emitido
por el IPCC, los cuales se encuentran basados en las emisiones de GEI. El IPCC nombra
los escenarios RCP bajo ciertas condiciones de concentración dependiendo del
Forzamiento Radiativo (FR) y la concentración de CO2 para el año 2100. En el caso de
GCM MPI-ESM-MR y el CORDEX, se encuentran tres escenarios: RCP 2.6, RCP 4.5 y RCP
8.5; de los cuales se obtuvo información a través del portal del CEDA para los datos
correspondientes a los registros históricos para las variables de precipitación y temperatura
para cada uno de los RCP.

9.4.2. Reducción de escala basada en la teoría del Caos

Para establecer los escenarios de cambio climático en la zona de estudio, se estructura un

código con base al modelo estadístico de la teoría del Caos que tiene en cuenta la
climatología global contenida en el GCM elegido y los datos globales con los registros que
se tienen en las series históricas de la cuenca, es decir, compara y asemeja los datos
globales con los registros de la zona, teniendo en cuenta las condiciones establecidas en
cada uno de los RCP. A partir de ello, se generan datos climatológicos futuros para la
variable de precipitación y temperatura, para evaluar el comportamiento del caudal bajo
diferentes escenarios de cambio climático, que repercuten en la oferta de la cuenca hasta
el año 2100. Es importante realizar la reducción de escala, ya que por medio de esta se
extrae información a escala local o regional de variables climáticas a partir de modelos de
alta resolución que contienen variables atmosféricas, lo que permite disminuir la limitación
que existe en la calidad de la información a escalas reducidas y su precisión.

En las gráficas que se presentan en el Anexo A, se puede observar que el modelo

CORDEX MPI_ESM-MR para los datos de la variable de precipitación, tuvo mejor ajuste
con respecto a los datos históricos registrados para cada COMID; también, se presentan
los datos del modelo original y en el Anexo B, las gráficas en donde se observa que el
modelo GCM MPI-ESM-MR para los datos de la variable de temperatura tuvo un mejor
ajuste con respecto a los datos historias registrados para cada COMID. Con estos

61
resultados, se obtuvieron las series de datos correspondientes a las variables de
precipitación y temperatura para cada subcuenca, limitada por las características de cada
escenario RCP, de modo que se establezcan tres diferentes escenarios de cambio
climático, ya que estos presentan variaciones en diferentes grados debido al cambio en la
concentración de GEI.

Como se observa en la Tabla 16 para cada escenario y cada COMID se presenta un

incremento o disminución de las variables de precipitación y temperatura tanto a mediados
de siglo (MOC) como a finales de siglo (EOC). En el caso de la variable de precipitación
para cada escenario se presenta un incremento de los eventos tanto a mediados de siglo
como para finales, en el caso del RCP 2.6 se presenta un incremento promedio del 1.83%
para el MOC y 4% para el EOC, para el RCP 4.5 6.53% para MOC y 6.9% para EOC, por
ultimo para el escenario RCP 8.5 hay incrementos del 7.43% para MOC y 7.97% para EOC.

Tabla 16 Variación de la precipitación bajo la reducción de escala.

PRECIPITACION
ESCENARIO COMID
INC/DEC MOC (%) INC/DEC EOC (%)
301388900 -1.00 0.50
RCP 2.6 301388800 -1.50 0.70
301413300 8.00 10.80
PROMEDIO 1.83 4.00
301388900 2.00 6.00
RCP 4.5 301388800 7.20 3.90
301413300 10.40 10.80
PROMEDIO 6.53 6.90
301388900 10.40 7.60
RCP 8.5 301388800 0.10 2.70
301413300 11.80 13.60
PROMEDIO 7.43 7.97
Fuente Autores

En la Grafica 16 se puede observar, la variación de los datos históricos promedio de

precipitación en comparación con el promedio de las series de datos, obtenidas a través de
la reducción de escala para cada COMID, teniendo en cuenta ello, se evidencia que en
periodos futuros el promedio de los datos aumenta, lo que indica un incremento en la
precipitación dentro de la zona de estudio.

Grafica 16 Precipitación promedio - Datos históricos VS series de escala reducida

62
 4 PRECIPITACIÓN PROMEDIO
3

mm
2

0
301413300 301388900 301388800
COMID
Serie Historica RCP 2.6

Fuente Autores

A continuación, se presentan los mapas de interpolación del comportamiento de la

precipitación en cada escenario de cambio climático, en el Mapa 11 para el escenario RCP
2.6, los mayores valores de precipitación se concentran en el COMID 301388800, es decir,
que en esa porción de la cuenca se han de presentar eventos de mayor magnitud en
periodos futuros que corresponden entre los 3 mm/día a 3.1 mm/día, en comparación con
el escenario actual, en el COMID 301388900 la precipitación oscila entre los 2.9 mm/día a
3 mm/día y en el COMID 301413300, se presentan los valores menores de precipitación
donde los eventos registran eventos de 2.3 mm/día a 2.4 mm/día. Con respecto al escenario
de estabilización, el cual corresponde al RCP 4.5 se observa en el Mapa 12 que la
precipitación se concentra en la parte alta de la cuenca en los COMID 301388900 y
301388800, donde los valores de precipitación oscilan entre los 2.9 mm/día y 3 mm/día,
para el COMID 301413300, se presentan valores de precipitación entre los 2.3 mm/día a
2.5 mm/día. Para el escenario RCP 8.5 se observa en el Mapa 13 que el comportamiento
de la precipitación es similar al escenario que se plantea en el RCP 4.5, ya que los valores
se aproximan a los 3.1 mm/día en los COMID de la parte alta de la cuenca, sin embargo, el
área que se abarca es mayor en comparación con éste, en la subcuenca con COMID
301413300, se presentan las precipitaciones con menor valor oscilando entre los 2.3
mm/día a 2.5 mm/día.

Mapa 11 Interpolación de precipitación bajo el escenario RCP 2.6

63
 Fuente Autores

Mapa 12 Interpolación de precipitación bajo el escenario RCP 4.5

Fuente Autores

Mapa 13 Interpolación de precipitación bajo el escenario RCP 8.5

64
 Fuente Autores

En la Tabla 17 para la variable de temperatura, se observa un incremento para los

escenarios RCP 2.6 y RCP 4.5, tanto para mediados de siglo como para finales; en
promedio 0.17% y 0.23% para el escenario RCP 2.6 y, 0.57% y 0.47% para el escenario
RCP 4.5, respecto al escenario RCP 8.5 la temperatura presenta una disminución de -
0.43% para el MOC y -0.60% para el EOC.

Tabla 17.Variación de la temperatura bajo la reducción de escala

TEMPERATURA
ESCENARIO COMID
INC/DEC MOC (%) INC/DEC EOC (%)
301388900 0.4 0.4
RCP 2.6 301388800 -0.3 -0.2
301413300 0.4 0.5
PROMEDIO 0.17 0.23
301388900 -0.1 -0.2
RCP 4.5 301388800 0.7 0.6
301413300 1.1 1
PROMEDIO 0.57 0.47
301388900 0.5 0.2
RCP 8.5
301388800 -0.7 -0.9

65
 TEMPERATURA
ESCENARIO COMID
INC/DEC MOC (%) INC/DEC EOC (%)
301413300 -1.1 -1.1
PROMEDIO -0.43 -0.60
Fuente Autores

En la Grafica 17 se observa la variación de temperatura promedio histórica en comparación

con la temperatura promedio que se presentan en los escenarios futuros, estos datos
oscilan entre los 12.5°C y los 13.7°C; como se puede evidenciar para los tres COMID la
temperatura tiende a elevarse, sin embargo para el escenario RCP 8.5 los valores
disminuyen y se presentan promedios menores a los valores históricos.

Grafica 17 Temperatura promedio - Datos históricos VS series de escala reducida

14
TEMPERATURA PROMEDIO
13,5

13
°C

12,5

12

11,5
301413300 301388900 301388800
COMID
Serie Historica RCP 2.6 RCP 4.5 RCP 8.5

Fuente Autores

A continuación, se presentan los mapas de interpolación del comportamiento de

temperatura en cada escenario de cambio climático, en el Mapa 14 para el RCP 2.6 en el
COMID 301388900, se presentan los datos con menor valor, estos oscilan entre 12.56°C a
12.73°C principalmente, por otro lado, para el COMID 301388800 los datos se aproximan
a los 13°C, en el último COMID 301413300, los valores de temperatura se acercan a los
13.73°C. con respecto al escenario de estabilización correspondiente al RCP 4.5 se observa
en el Mapa 15 que la temperatura presenta un incremento significativo en comparación con
el escenario anterior ya que se evidencia que para el COMID 301388900 los valores
aumentan y éstos fluctúan entre los 13.26°C a los 13.53°C, para el COMID 301388800 la
temperatura oscila entre los 12.98°C y los 13.25°C, por último la subcuenca con COMID
301413300 presenta temperaturas de 13.54°C a 13.81°C. Para el escenario RCP 8.5, se
observa en el Mapa 16 que los valores de temperaturas se reducen, presentando un
comportamiento similar al del RCP 4.5, para el COMID 301388900 de la parte alta de la
cuenca, la temperatura oscila entre los 12.56°C y los 12.7°C, el COMID 301388800
presenta temperaturas de 13°C y el COMID 301413300 valores de hasta 13.56°C.

66
Mapa 14 Interpolación de temperatura bajo el escenario RCP 2.6

Fuente Autores

Mapa 15 Interpolación de temperatura bajo el escenario RCP 4.5

Fuente Autores

67
 Mapa 16 Interpolación de temperatura bajo el escenario RCP 8.5

Fuente Autores

10. OBTENCIÓN DE LA DEMANDA

10.1. Obtención de la demanda actual y futura

Para la obtención de la demanda actual de la cuenca, fue necesario recurrir a los registros
de la Corporación Autónoma Regional (CAR) los cuales se encuentran en el Anexo D;
esto, con el fin de obtener los valores correspondientes a las concesiones asignadas dentro
de la zona de estudio para cada sector económico, en este caso para el sector doméstico,
agrícola, pecuario e industrial.

Para la demanda futura se realizó la proyección de las demandas para cada uno de los
sectores económicos, teniendo en cuenta los registros históricos del Producto Interno Bruto
(PIB), datos del Departamento Administrativo Nacional de Estadística DANE y el Banco de
la Republica, en este caso tenemos un porcentaje para las actividades económicas que se
realizan en la cuenca del 2%.

PIB 0.02

68
 Tabla 18 Calculo de la población proyectada hasta el 2020 por Municipios.

Población Población Población Población

Proyectada de Proyectada de Proyectada de Proyectada
MUNICIPIO
1985 a 2020 2021 a 2040 2041 a 2070 2071 a 2100
(hab) (hab) (hab) (hab)
El Rosal 19253 24026 31186 38346
Facatativá 144149 194906 271042 347178
Mosquera 93461 137436 203398 269360
Subachoque 17712 17228 16502 15776[DIA4]
Madrid 85090 114401 158369 202336
Fuente Autores

Teniendo en cuenta los datos obtenidos en la Tabla 18 se observa que en el municipio de

Subachoque la población disminuye a lo largo del periodo de los años, esto se manifiesta
por factores externos. Adicional una parte de la población que se encuentra allí realizan
únicamente actividades laborales por lo cual se desplazan diariamente a este municipio,
cabe destacar que, aunque la población disminuye no se ve afectada la cuenca, es decir, a
menor población, menor consumo de agua y mayor oferta hídrica en este sector.

Como se mencionó anteriormente, para obtener la demanda actual se contó con

información más precisa y actualizada por parte de la Corporación Autónoma Regional
(CAR), para realizar la proyección se tuvo en cuenta el PIB el cual, es una magnitud
macroeconómica que expresa el valor monetario de la producción de bienes y servicios de
demanda final de un país o región durante un periodo determinado [¡Error! No se
encuentra el origen de la referencia.], así mismo, se contó con información de las
principales actividades económicas y se unifico por municipios como se observa en la Tabla
19.

Adicional se determino que la demanda Industrial se comporta de manera parecida para la

mayoría de los municipios con excepción de Mosquera, aunque este presenta presenta
mayor numero de industrias se refleja un menor consumo de agua, debido a que las
empresas que se encuentran ubicadas en este sector como Italcol S.A., Purina, Arinas
Lobo, entre otras tienen un proceso de producción en seco, es decir, no requieren de agua
para su producción. Además se evidencia que dichas empresas emplean métodos de
reutilización de agua o sistema de aguas lluvias utilizada para el área administrativa lo cual
hace que el consumo de agua por parte directa de la cuenca disminuya.

Tabla 19 Clasificación de demanda hasta el año 2020 por sectores en litro/año.

Demanda Demanda Demanda Demanda

MUNICIPIO Doméstica Agrícola Pecuario Industrial
(litro/año) (litro/año) (litro/año) (litro/año)
Bajaca 2343755520 13352610456 2557859731 110344464
El Rosal 2350551528 13488338246 2561707123 110344464
Funza 2192871528 13335130051 2554706131 110344464

69
 Madrid 2352128328 13508521286 2563283923 110344464
Mosquera 79921684.8 11530324771 2277126259 0[DIA5]
Subachoque 2352759048 13594929926 2567572819 110344464
Fuente Autores

El cálculo de la demanda de agua por uso doméstico se realizó con el método aritmético
según lo establecido en el artículo B “Sistemas de Acueducto” del Reglamento Técnico de
Agua Potable y Saneamiento Básico RAS 2017, en el cual se establece que el método de
cálculo para la proyección de la población depende del nivel de complejidad del sistema.
Se calcula la población utilizando el modelo matemático que mejor se ajuste al
comportamiento histórico de la población.

Se inició tomando la población correspondiente al censo inicial del año 1983 y la población
proyectada en el censo DANE para 2020.

Ecuación 3 Cálculo de la demanda doméstica.

Pf=(Puc+((Puc-Pci)/(Tuc-Tci)))*(Tf-Tuc)[DIA6]
Fuente Autores

Donde:
 Pf: población correspondiente al año de la proyección
 Pci: población correspondiente al censo inicial
 Puc: población correspondiente a la proyección del DANE
 Tuc: año correspondiente al último año proyectado por el DANE
 Tf: Tasa de crecimiento observado en el período
 Tci: año correspondiente al censo inicial con información

10.1.1. Sector agrícola

La demanda para este sector hace referencia a los requerimientos hídricos de los cultivos,
que deben ser suplidos de manera natural o artificial. Alrededor de la Cuenca del Rio
Balsillas, se encuentran cultivos de arroz, papa, hortalizas y frutales los cuales
corresponden al 80% del área de cultivo y los cultivos de flores e invernaderos
corresponden al 30% del área de cultivo dejando el 10% restante para el área de pastos
manejados.

A continuación, se presentan las demandas de agua en m3/s proyectadas del sector

agrícola para la Cuenca del Rio Balsillas.

Tabla 20 Demanda de agua en m3/s del sector agrícola para cada municipio

Demanda Agrícola (m3/s)

MUNICIPIO
2040 2070 2100
Bojaca 0.553 0.958 1.736
El Rosal 0.559 0.968 1.753
Funza 0.553 0.957 1.733

70
 Madrid 0.560 0.969 1.756
Mosquera 0.478 0.827 1.499
Subachoque 0.563 0.976 1.767
SUMATORIA 3.266 5.655 10.243
Fuente Autores

Como se puede observar en la Tabla 20, el municipio que requiere menor cantidad de agua
para suplir las necesidades de riego en sus cultivos es Mosquera, con una demanda de
agua de 2.804 m3/s ya que en esta zona se presentan mayores periodos de precipitación
lo cual ayuda a suplir esta demanda y su técnica de riego para cultivos de lechuga,
zanahoria, apio y arveja (predominan cultivos de verduras en la zona), ha mejorado con la
implementación de programas de ahorro y uso eficiente del agua lo que permite un mejor
uso del recurso hídrico. El municipio que requiere mayor cantidad de agua para suplir las
necesidades de riego en sus cultivos es Madrid, con una demanda de agua de 3.29 m3/s,
ya que este es el mayor municipio floricultor de Colombia por ende tiene una demanda de
agua más alta en comparación con los demás municipios que hacen parte de la cuenca del
Rio Balsillas. Para los municipios de Bojaca, El Rosal, Funza y Subachoque, la demanda
de agua es de 3.25 m3/s, 3.28 m3/s, 3.24 m3/s y 3.28 m3/s respectivamente, cuentan con
una demandan de agua similar ya que su sistema de riego comparte patrones debido a que
sus cultivos principales son de papa, zanahoria, maíz, lechuga, hortalizas y trigo.

10.1.2. Sector pecuario

La demanda para este sector hace referencia a los requerimientos hídricos de la actividad
ganadera, consistentes en la crianza de animales para diferentes fines como lo son
doméstico, comercialización y aprovechamiento económico. Alrededor de la Cuenca del Rio
Balsillas, las especies de mayor producción son los bovinos, porcinos, equinos, conejos y
aves correspondientes al 60%, 20%, 10%, 5% y 5% respectivamente.

A continuación, se presentan las demandas de agua en m3/s proyectadas del sector

pecuario para la Cuenca del Rio Balsillas.

Tabla 21 Demanda de agua en m3/s del sector pecuario para cada municipio
Demanda Pecuario (m3/s)
MUNICIPIO
2040 2070 2100
Bojaca 0.106 0.184 0.332
El Rosal 0.106 0.184 0.333
Funza 0.106 0.183 0.332
Madrid 0.106 0.184 0.333
Mosquera 0.094 0.163 0.296
Subachoque 0.106 0.184 0.334
SUMATORIA 0.625 1.082 1.960
Fuente Autores

Como se puede observar en la Tabla 21, la demanda del sector pecuario para el primer
periodo de tiempo del año 2015 al 2040 es similar para todos los municipios, exceptuando
Mosquera, para el periodo de tiempo del año 2041 a 2070 la demanda aumenta 0.078 m3/s

71
y para el periodo de tiempo del año 2071 a 2100 se evidencia un aumento de 0.227 m3/s.
Para este sector, los municipios que menor cantidad de agua requieren son Mosquera y
Funza, equivalentes al 15.1% y 16.9% respectivamente del consumo hídrico. Mosquera
tiene una demanda de agua de 0.553 m3/s, este municipio se destaca por contar con
grandes haciendas en donde aprovechan la cría y ceba de especies menores como cerdos
y pollos, en ganadería se explotan la raza Holstein y Normanda para lechería y carnes, así
mismo, se destacan por contar con buenas practicas del recurso hídrico, logrando una
disminución en sus consumos. Funza, tiene una demanda de agua de 0.621 m3/s, en este
municipio cuenta con un sistema de producción predominante referente a la ganadería de
leche explotando la raza Holstein (bovinos), seguidos por los porcinos, equinos y conejos.
Adicionalmente estos dos municipios cuentan con una serie de restricciones de consumo
de agua lo cual contribuye al bajo porcentaje de demanda hídrica que necesitan con
respecto a los demás municipios que hacen parte de la Cuenca del Rio Balsillas.

Para los municipios de Bojaca, El Rosal, Madrid y Subachoque se tiene una demanda de
agua de 0.622 m3/s, 0.623 m3/s, 0.623 m3/s y 0.624 m3/s respectivamente lo cual equivale
a un 17% de consumo hídrico para cada municipio.

10.1.3. Sector industrial

La demanda para este sector, hace referencia a los requerimientos hídricos de la actividad
industrial para cada municipio. Alrededor de la Cuenca del Rio Balsillas, las actividades
industriales representan un porcentaje mínimo de consumo de agua.

A continuación, se presentan las demandas de agua en m3/s, proyectadas del sector

industrial para la Cuenca del Rio Balsillas.

Tabla 22 Demanda de agua en m3/s del sector industrial para cada municipio

Demanda Industrial (m3/s)

MUNICIPIO
2040 2070 2100
Bojaca 0.005 0.008 0.014
El Rosal 0.005 0.008 0.014
Funza 0.005 0.008 0.014
Madrid 0.005 0.008 0.014
Mosquera 0.000 0.000 0.000
Subachoque 0.005 0.008 0.014[DIA7]
SUMATORIA 0.023 0.040 0.072
Fuente Autores

Como se puede observar en la Tabla 22, la demanda del sector industrial desde el primer
periodo de tiempo equivalente al año 2015 a 2040, hasta el último periodo de tiempo
equivalente al año 2071 a 2100, aumenta 0.003 m3/s y 0.009 m3/s, este sector representa
el 4% de la demanda hídrica total de la Cuenca del Rio Balsillas. Para las actividades
industriales (empresas), se debe tener en cuenta, que este consumo es bajo debido a que
la mayoría de las empresas que se encuentran ubicadas en estos municipios, tienen su

72
propio sistema de tratamiento de agua y de reusó, para actividades no domesticas
(sanitarios o riego), o para su mismo proceso productivo, así como el aprovechamiento de
las aguas lluvias. Las empresas que predominan son las dedicadas a la producción de
alimentos para humanos, alimentos para mascotas y materiales de construcción, aunque
algunas no requieren demanda hídrica debido a que sus procesos se realizan en seco, se
tiene un consumo mínimo de agua para su proceso productivo, como cuentan con su propio
sistema de tratamiento y reusó de agua, se les facilita suplir la necesidad de demanda
hídrica lo cual, baja los consumos referentes a fuentes aledañas. Dentro de este sector se
cuenta con empresas como Colceramiza S.A., Doria S.A.S, Talcos, Purina Nestle S.A.,
entre otras.

10.1.4. Sector domestico

La demanda para este sector, hace referencia a los requerimientos hídricos que necesita
un habitante para satisfacer sus necesidades básicas, para hallar la dotación máxima por
habitante, se revisó el Titulo B del RAS en donde se obtiene la dotación neta en L/hab*día,
adicional, se define el nivel de complejidad del sistema según lo establecido en el Titulo A
del RAS el cual, es determinado según el número de habitantes en una zona. Para
determinar el nivel de complejidad se tiene en cuenta la información de la Tabla 18, en
donde se observa la población proyectada al año 2100 por los municipios que hacen parte
de la Cuenca del Rio Balsillas, por lo tanto, los datos poblacionales corresponden de
complejidad Alto, asumiendo una dotación neta de 140 L/hab*día para climas templados y
fríos.

Tabla 23 Determinación de dotación neta

Nivel de complejidad Alto

Clima Templado y frio
L
Dotación neta 140 hab∗dia
Fuente Autores

A continuación, se presentan las demandas de agua en m3/s proyectadas del sector

doméstico para la Cuenca del Rio Balsillas.

Tabla 24 Demanda de agua en m3/s del sector doméstico para cada municipio
Demanda doméstica (m3/s)
MUNICIPIO
2040 2070 2100
Bojaca 0.051 0.075 0.100
El Rosal 0.399 0.587 0.775
Funza 0.291 0.474 0.656
Madrid 0.037 0.051 0.061
Mosquera 0.296 0.316 0.158
Subachoque 0.162 0.347 0.118
SUMATORIA 1.236 1.849 1.868
Fuente Autores

73
Como se puede observar en la Tabla 24, la demanda del sector doméstico se incrementa
sucesivamente al pasar de los años, desde el primer periodo de tiempo equivalente al año
2015 a 2040, hasta el último periodo de tiempo equivalente al año 2071 a 2100 aumenta
0.188 m3/s y 0.376 m3/s. Se debe tener en cuenta, que el consumo doméstico no solo hace
alusión al consumo realizado en las viviendas sino también, al consumo que se realiza en
los establecimientos prestadores de servicios para actividades de turismo (restaurantes,
hoteles, etc.). Se observa que el mayor requerimiento de demanda hídrica es para el
municipio de El Rosal y de Funza equivalente a un 35.5% y un 28.7% respectivamente, en
comparación con los demás municipios, esto se debe al consumo que se realiza por parte
de los habitantes de la zona y al consumo que realizan los turistas ya que particularmente,
estos municipios tienen un alto flujo de visitantes. Los municipios de Bojaca, Mosquera y
Subachoque tienen un requerimiento de agua del 4.6%, 15.5% y 12.7% respectivamente.
Con respecto al municipio que menor demanda hídrica requiere es Madrid, con un 3% ya
que, a pesar de tener una cantidad de habitantes considerable, cuentan con buenas
prácticas de uso y ahorro de agua y con ciertas restricciones para no exceder sus
consumos.

11. ANÁLISIS DE LA VARIABILIDAD DE OFERTA HÍDRICA NATURAL

Para el análisis de la variabilidad de oferta hídrica, se tuvo en cuenta la información

suministrada por el documento del IDEAM titulado Estudio nacional del agua- ENA 2014 en
el cual, se encuentran los indicadores que caracterizan el régimen de agua de una zona
determinada. El ENA, estima las categorías para IRH, IUA, IVH, proyectadas hasta el año
2020, en el Anexo 3 “Indicadores cabeceras municipales” se encuentra la información de
los rangos para la subzona hidrográfica para el departamento de Cundinamarca el rio
Bogotá. Según el análisis realizado y estipulado en el documento se tiene un Índice de
retención y regulación hídrica (IRH) bajo, lo cual indica que la cuenca cuenta con una baja
capacidad para retener y regular la humedad dentro de esta, un Índice de uso de agua (IUA)
muy bajo, lo que indica que la presión de demanda de la cuenca no es significativa con
respecto a la oferta disponible y un Índice de vulnerabilidad por desabastecimiento (IVH)
medio [¡Error! No se encuentra el origen de la referencia.]. Esta información se tuvo en
cuenta para la estimación de dichos indicadores aplicados a la cuenca del rio Balsillas para
tres periodos de tiempo diferentes, los cuales corresponden desde el año 2015 al 2040,
2041 al 2070 y 2071 al 2100.

11.1. Curva de duración de caudales

Las curvas de duración o distribución de frecuencias acumuladas de los caudales medios

diarios indican, el porcentaje (%) de tiempo durante el cual los caudales igualan o exceden
un valor dado; la forma y pendiente de estas curvas, refleja la capacidad de regulación de
la cuenca en la unidad de tiempo considerada, además, generan información sobre el
porcentaje del tiempo en el que el río lleva un caudal superior o inferior a un determinado
valor, pero no reflejan la distribución o secuencia de dicho periodo ni el momento del año
en que se produce [¡Error! No se encuentra el origen de la referencia.].

Al tener la serie de caudales proyectados hasta el año 2100, se dividen estos datos diarios
en tres periodos de tiempo, el primer periodo de tiempo es del año 2015 a 2040, el segundo

74
periodo de tiempo es del año 2041 al 2070 y finalmente el tercer periodo de tiempo es del
2071 al 2100. Posteriormente, para la construcción de la curva de duración de caudales
medios diarios se requiere:

a) Organizar los datos de caudal diario de mayor a menor valor

b) Enumerar los datos de caudal desde el número uno (1) para identificar su posición
en la serie, donde 1 corresponde al caudal de mayor valor. La letra m corresponde
al valor total de registros de la serie de caudal
c) Calcular la probabilidad de excedencia con las siguientes ecuaciones:

Ecuación 4 Cálculo de probabilidad de excedencia.

m m
𝐏 = ∗ 100 o 𝐏 = ∗ 100
n n+1
Fuente [¡Error! No se encuentra el origen de la referencia.]

Donde:
 P: probabilidad de excedencia
 m: posición de datos registrados
 n: total de datos
d) Graficar en el eje y los datos de caudal (m3/s) y en el eje x los datos obtenidos por
la función de la probabilidad de excedencia

A continuación, se presentan las curvas de duración de caudales para los tres periodos de
tiempo de 2015 a 2040, 2041 a 2070 y 2071 a 2100; en estas graficas es posible observar
que la curva se desplaza hacia arriba en el eje de las ordenadas del 0% al 30% de tiempo
de excedencia, indicando un posible aumento de los eventos extremos como se ha
mostrado también en las gráficas de precipitación; a su vez, se presenta una disminución
con respecto al valor observado del 45% al 100% del tiempo, dando como resultado una
reducción en el caudal.

Grafica 18 Curva de duración de caudales medios diarios del año 2015 a 2040

2015 al 2040
14
12
10
Q (m3/s)

8
6
4
2
0
0 16 32 48 64 80 96
% P=(m/n)*100
Fuente Autores

Grafica 19 Curva de duración de caudales medios diarios del año 2041 a 2070

75
 2041 al 2070
14
12
10

Q (m3/s)
8
6
4
2
0
0 16 32 48 64 80 96
% P=(m/n)*100

Fuente Autores

Grafica 20 Curva de duración de caudales medios diarios del año 2071 a 2100

2071 al 2100
14

12

10
Q (m3/s)

0
0 16 32 48 64 80 96
% P=(m/n)*100

Fuente Autores

11.2. Índice de retención y regulación hídrica (IRH)

Es un indicador que varía entre 0 y 1, siendo los valores más bajos los que se interpretan
como de menor regulación. Para el cálculo de este indicador, se tiene en cuenta la curva
de duración de caudales diarios, entre más cercano a 1 sea el valor, indica una mayor
regulación, representando una mayor homogeneidad en los caudales, y entre más cercano
a 0 sea el valor, indica un mayor contraste entre los caudales más bajos y más altos que
genera la cuenca, evidenciando rangos de tiempo con magnitudes características altas y/o
bajas y de mayor variabilidad.

76
 Ecuación 5 Cálculo del índice de retención hídrica.
𝑉𝑝
𝑰𝑹𝑯 =
𝑉𝑡
Fuente [¡Error! No se encuentra el origen de la referencia.]

Donde:

 IRH: índice de retención y regulación hídrica

 Vp: volumen presentado por el área que se encuentra por debajo de la línea de
caudal medio en la curva de duración de caudales diarios.
 Vt: volumen total representado por el área bajo la curva de duración de caudales
diarios

Tabla 25 Rangos y categorías del Índice de Retención y Regulación Hídrica

Rango de
Categoría Características
valores IRH
Capacidad de la cuenca para retener y regular muy
> 0.85 Muy Alto alta
0.75-0.85 Alto Capacidad de la cuenca para retener y regular alta
Capacidad de la cuenca para retener y regular
0.65-0.75 Medio media
0.50-0.65 Bajo Capacidad de la cuenca para retener y regular baja
Capacidad de la cuenca para retener y regular muy
< 0.5 Muy Bajo baja
Fuente [¡Error! No se encuentra el origen de la referencia.]

Con base en los rangos establecidos por el ENA y los datos obtenidos del proyecto se
determinó que la cuenca tiene una capacidad de retención de humedad media como se
observa en la Tabla 26, es decir, que para el periodo de tiempo del año 2015 al 2100 se
presentará un comportamiento continuo de retención, con un rango de 0.70 a 0.73
distribuida entre la parte alta, media y baja de la cuenca dando como resultado el equilibrio
del recurso hídrico para su disposición en los diferentes sectores económicos como lo son
el sector agrícola, pecuario, doméstico e industrial. Sin embargo, se debe tener en cuenta
la variación y posibles cambios que pueda sufrir la cuenca y estos cambios se inclinarían
hacia una retención hídrica alta, ya que por la ubicación geográfica de esta, se evidencia
un equilibrio entre los periodos secos y húmedos presentes.

Tabla 26 Índice de retención y regulación hídrica para la cuenca del Rio Balsillas

Índice de retención y regulación hídrica (IRH)

PERIODO DE TIEMPO VALOR CATEGORIA
2015 A 2040 0.71 Medio
2041 A 2070 0.72 Medio

77
 2071 A 2100 0.72 Medio
Fuente Autores

11.2.1. Caudal ambiental

El caudal ambiental, se define como el volumen de agua necesario en términos de calidad,

cantidad, duración y estacionalidad para el sostenimiento de los ecosistemas acuáticos y
para el desarrollo de las actividades socioeconómicas de los usuarios aguas abajo de la
fuente, de la cual dependen los ecosistemas. Los valores de la curva de duración de
caudales, fueron empleados adicionalmente para la estimación del caudal ambiental para
cada uno de los periodos de tiempo.

Para la determinación del caudal ambiental, debe ponderarse una variada gama de
resultados, desde la protección ambiental hasta las necesidades de industrias y de
personas, en este caso el caudal ambiental calculado para la cuenca en los tres periodos
de tiempo correspondientes al año 2015 a 2040, 2041 a 2070 y 2071 a 2100 es de 1.173
m3/s, 0.013 m3/s y 0.995 m3/s respectivamente. Como se observa en la Tabla 27, el periodo
de tiempo con menor uso y aprovechamiento del caudal ambiental máximo corresponde a
los años del 2041 a 2070 esto debido a factores naturales y antrópicos como crecimiento
poblacional, expansión de los sectores económicos entre otros factores; a diferencia de un
periodo actual, 2015 a 2040 donde se presenta una demanda de agua equilibrada y con
mayor volumen de caudal ambiental aprovechable ya que los municipios aledaños que se
abastecen cuenta con un porcentaje de población bajo y un mínimo porcentaje tomado para
sector industrial y usos agrícolas. Para el periodo de tiempo correspondiente al 2071 a 2100,
es posible que se vuelva aumentar el caudal ambiental, ya que se pueden presentar
periodos húmedos más largos debido a su ubicación geográfica, además, se espera hacer
uso y aprovechamiento de aguas lluvias por parte de las empresas aledañas que cada vez
más, invierten en proyectos para reutilización del agua y control de consumo energético.

Tabla 27 Caudal ambiental en m3/s para cada periodo de tiempo en la Cuenca del
Rio Balsillas

Caudal ambiental (m3/s)

2015 a 2040 1.173
2041 a 2070 0.013
2071 a 2100 0.995
Fuente Autores

11.3. Índice de uso de agua (IUA)

Se continúa considerando aun solamente la oferta superficial en el cálculo del indicador,

para lo cual se emplea la siguiente expresión:

Ecuación 6 Cálculo de índice de uso de agua.

Dh
𝐈𝐔𝐀 = ∗ 100
Oh
Fuente [¡Error! No se encuentra el origen de la referencia.]

78
Donde:

 Dh: demanda hídrica sectorial

 Oh: Oferta hídrica superficial disponible (esta última resulta de la cuantificación de
la oferta hídrica natural sustrayendo la oferta correspondiente al caudal ambiental).
 Dh: Σ (Volumen de agua extraída para usos sectoriales en un periodo determinado)

Ecuación 7 Cálculo de demanda hídrica sectorial.

Dh = Ch + Csp + Csm + Css +Cea + Ce +Ca +Aenc

Fuente [¡Error! No se encuentra el origen de la referencia.]

Donde:

 Dh: Demanda hídrica

 Ch: Consumo humano o domestico
 Csp: Consumo del sector agrícola
 Csm: Consumo del sector industrial
 Css: Consumo del sector servicios
 Ce: Consumo del sector energía
 Ca: Consumo del sector acuícola
 Aenc: agua extraída no consumida

11.3.1. Demanda hídrica

La demanda de agua estimada, corresponde a la cantidad o volumen de agua usado por

los sectores económicos y la población. Considera el volumen de agua extraído o que se
almacena de los sistemas hídricos y que limita otros usos; contempla el volumen utilizado
como materia prima, como insumo y el retornado a los sistemas hídricos.

El cálculo de la demanda hídrica para la cuenca del rio Balsillas se obtuvo mediante la suma
del consumo de cada sector, es decir, la demanda requerida para consumo generado por
el sector agrícola, pecuario, industrial y doméstico. En la Tabla 28 se observan las
demandas requeridas para los tres periodos de tiempo correspondientes al año 2015 a
2040, 2041 a 2070 y 2071 a 2100, en donde se obtuvo un valor 5.1496 m3/s, 8.6262 m3/s
14.1434 m3/s, para cada uno de ellos; se evidencia una un aumento con una variación de
3 m3/s a 4 m3/s entre los periodos de tiempo, esto se relaciona al aumento respectivo de
fenómenos tanto naturales como factores antrópicos.

Tabla 28 Demanda total de agua en m3/s para cada periodo de tiempo

Demanda total de la
cuenca (m3/s)
2015 a 2040 5.1496
2041 a 2070 8.6262
2071 a 2100 14.1434

79
 Fuente Autores

11.3.2. Oferta hídrica

La oferta hídrica disponible, es el volumen de agua promedio que resulta de sustraer a la

oferta hídrica total superficial, el volumen de agua que garantiza el uso para el
funcionamiento de los ecosistemas y de los sistemas fluviales, y en alguna medida un
caudal mínimo para usuarios que dependen de las fuentes hídricas asociadas a estos
sistemas. Esta se determina a partir de la oferta total menos el caudal ambiental.

Ecuación 8 Cálculo de oferta hídrica superficial disponible.

Oh = Ohtotal – OQamb
Fuente [¡Error! No se encuentra el origen de la referencia.]

Donde:

 Oh total: volumen total de agua superficial en una unidad de análisis espacial y

temporal determinada
 Oqamb: volumen de agua correspondiente al caudal ambiental en la misma unidad
de análisis espacial y de tiempo de la oferta total.

Como se observa en la Tabla 29, se obtuvieron las ofertas hídricas disponibles en m3/s
para cada periodo de tiempo, la menor oferta se da para el año 2015 a 2040;
correspondiente al 0.812 m3/s el cual equivale al 21% de la oferta total al año 2100. El
periodo de tiempo con mayor oferta es del año 2041 al 2070, correspondiente a 2.019 m3/s
equivalente al 52% de la oferta total y finalmente se tiene que para el periodo de tiempo del
año 2071 al 2100, la oferta es de 1.033 m3/s equivalente al 27% de la oferta total.

Tabla 29 Oferta hídrica para cada periodo de tiempo en la Cuenca del Rio Balsillas

Oferta hídrica
disponible (m3/s)
2015 a 2040 0.812
2041 a 2070 2.019
2071 a 2100 1.033
Fuente Autores

Tabla 30 Rangos y categorías del Índice del Uso del Agua

Rango (Dh/Oh)*100 Categoría

Significado
IUA IUA
> 100 Critico La presión supera las condiciones de la oferta
La presión de demanda es muy alta con
50.01-100 Muy Alto respecto a la oferta disponible
La presión de demanda es alta con respecto a
20.01-50 Alto la oferta disponible

80
 La presión de demanda es moderada con
10.01-20 Moderado respecto a la oferta disponible
La presión de demanda es baja con respecto a
1.0-20 Bajo la oferta disponible
La presión de demanda no es significativa con
<1 Muy Bajo respecto a la oferta disponible
Fuente [¡Error! No se encuentra el origen de la referencia.]

Mediante el cálculo del índice del uso del agua (IUA) que se obtiene de la demanda hídrica
sobre la oferta hídrica; se presentan los siguientes valores para diferentes periodos de
tiempo correspondientes, del año 2015 al 2040 con un valor de 633.881, del 2041 al 2070
con un valor de 427.286 y para el periodo de 2071 al 2100 con un valor de 1369.827; con
una clasificación de Critico, dentro de la categoría del IUA, ya que el rango se encuentra
por encima de 100, lo cual significa que la presión supera las condiciones de la oferta, es
decir, que la demanda hídrica por los sectores económicos es más alta que la oferta hídrica
y por ende, la cuenca no alcanza a suplir las necesidades requeridas.

Tabla 31 Índice de uso de agua para la cuenca del Rio Balsillas

Índice de uso de agua (IUA)

2015 a 2040 633.88 Critico
2041 a 2070 427.28 Critico
2071 a 2100 1369.82 Critico
Fuente Autores

11.4. Índice de vulnerabilidad por desabastecimiento (IVH)

El IVH se determina a través de una matriz de relación de rangos del Índice de regulación
hídrica (IRH) y el Índice de uso de agua (IUA). Las categorías de este Índice se presentan
en la siguiente tabla:

Tabla 32 Índice de vulnerabilidad al desabastecimiento hídrico IVH

Vulnerabilidad del recurso hidricorelacion IRH-IUA

IUA Extremo = porcentaje
Índice de regulación
(Oferta/Demanda)
Rango Categoría Alta Moderado Baja Muy Baja
<1 Muy Bajo Muy Baja Baja Media Media
1-10 Bajo Baja Baja Media Media
10-20 Moderado Media Media Alta Alta
20-50 Alto Media Alta Alta Muy Alta
50-100 Muy Alto Media Alta Alta Muy Alta
> 100 Critico Muy Alta Muy Alta Muy Alta Muy Alta
Fuente [¡Error! No se encuentra el origen de la referencia.]

81
 Tabla 33 Índice de uso de agua extremo
IUA Extremo= porcentaje
(Demanda/oferta)
2015 a 2040 6.33 Bajo
2041 a 2070 4.27 Bajo
2071 a 2100 13.69 Moderado
Fuente Autores

Tabla 34 Índice de vulnerabilidad por desabastecimiento del Rio Balsillas

Índice de vulnerabilidad por desabastecimiento (IVH)

Periodo de IUA Extremo= porcentaje Índice de
IVH
tiempo (Demanda/oferta) regulación
2015 a 2040 6.33 Bajo Moderado Baja
2041 a 2070 4.27 Bajo Moderado Media
2071 a 2100 13.69 Moderado Moderado Media
Fuente Autores

Con base en el estudio realizado por el ENA del año 2014, se determinó que la cuenca del
Rio Balsillas presento las siguientes características: una demanda hídrica que supera la
oferta, factor que incluye en la matriz del índice de vulnerabilidad, donde presenta una
categoría de nivel medio; índice del uso del agua y el índice de retención hídrica presentan
una categoría baja, al comparar estos resultados con los obtenidos mediante la proyección,
se evidencian aspectos similares y se obtiene que la cuenca presentará características
similares; es decir, un déficit para suplir las necesidades de abastecimiento requerido por
la población y en el transcurso del periodo de tiempo correspondiente al año 2015 a 2040
presentará un índice de vulnerabilidad bajo, debido a que en este periodo se presenta
mayor demanda en relación al periodo de tiempo de los años 2041 a 2070 y 2071 a 2100,
donde se prevé un índice de vulnerabilidad con un rango medio igual que el estimado en el
año 2014.
12. IMPACTO SOCIAL

El clima de la tierra ha tenido un cambio notorio a lo largo de su evolución, ya sea por

factores naturales o por factores antrópicos que producen cambios a mediano y largo plazo
sobre los ciclos y comportamientos naturales. Es por esto, que las progresivas
transformaciones ambientales requieren de estrategias colectivas para afrontar las
amenazas y riesgos, que conllevan a las poblaciones humanas y a los ecosistemas a
estados críticos de vulnerabilidad y cambio constante.

En Colombia, los impactos y transformaciones ambientales debido al cambio climático son

cada vez más evidentes, por lo cual es necesario generar estrategias globales, nacionales
y locales para reducirlos como vía hacia la adaptación, es por esto que se hace necesario
combatir los efectos del cambio climático a partir del análisis de estrategias ya existentes,
tanto a nivel de las comunidades en las que el conocimiento tradicional y ancestral ha sido
históricamente determinante en su adaptación al cambio climático, como en el ámbito
científico y técnico. En este orden de ideas, se ha podido evidenciar un crecimiento

82
constante a nivel poblacional el cual genera un alto consumo de recursos, en este caso, el
recurso hídrico, volviéndolo vulnerable hasta el punto de estimar un déficit en la oferta
producida por este y alteraciones en la calidad de vida de las personas.

En nuestro caso de estudio, la Cuenca del Rio Balsillas se tuvo en cuenta la demanda
requerida por los sectores económicos a nivel municipal, ya que estos, están directamente
asociados a la disponibilidad de agua. Por esta razón, se busca obtener la información
adecuada que permita dar un diagnostico a la comunidad con el fin de que se tomen las
medidas preventivas al estado y disponibilidad del recurso hídrico de la zona para un
adecuado desarrollo de actividades económicas y sociales sin mayores afectaciones. Es
por esto, que este proyecto se convierte en una herramienta de apoyo para las autoridades
ambientales competentes beneficiando a la comunidad con la toma de decisiones y
actividades que logren mitigar los impactos y logren un equilibrio en relación de oferta y
demanda hídrica, así, como la implementación de nuevos modelos y metodologías que
sean útiles a la hora evaluar los diferentes ecosistemas hídricos.

13. CONCLUSIONES

 Para la cuenca río Balsillas el modelo que más se ajustó en relación a la variable de
precipitación fue el MPI-ESM-MR CORDEX y para temperatura el GCM MPI-ESM-
MR, ya que presentaron los mejores resultados en cuanto a la comparación de las
medidas estadísticas y físicas.
 Bajo las condiciones de reducción de escala, para cada escenario la variable de
temperatura presento un cambio respecto a los datos históricos, para el escenario
RCP 2.6 la temperatura a mitad de siglo tendrá un incremento promedio de 0.17%
y a final de siglo del 0.23%, para el escenario RCP 4.5 el incremento será de 0.57%
para mitad de siglo y 0.47% para finales de siglo, finalmente se observara para el
escenario RCP 8.5, una disminución del -0.43% para mitad de siglo y -0.60% para
final de siglo; para estos escenarios la temperatura oscilara entre 12.5°C y 13.7° C.
 Bajo las condiciones de reducción de escala, para cada escenario la variable de
precipitación, se evidencio un cambio mínimo entre los escenarios y los datos
históricos. Para cada RCP se evidencio un incremento con comportamiento lineal,
para el escenario RCP 2.6 la precipitación a mitad de siglo presentará un aumento
promedio de 1.83% y a final de siglo del 4%, para el escenario RCP 4.5 el
incremento será de 6.53% para mitad de siglo y 6.9% para finales de siglo, por último
las condiciones del escenario RCP 8.5, presentan un incremento de 7.43% para
mitad de siglo y 7.97% para final de siglo; para estos escenarios la precipitación
oscilara entre 2.3 mm/día a 2.5 mm/día.
 El software Hydro-BID, presento en el proceso de modelación y calibración
condiciones medias y mínimas similares en la Cuenca del Rio Balsillas teniendo en
cuenta la información suministrada por la estación limnimetrica Puente Manrique,
esto se evidencio en los resultados de simulación de la escorrentía y el registro
histórico de caudales que presento esta estación. Cabe resaltar que las métricas
arrojadas por este software a nivel mensual presentaron de mejor manera las
condiciones de la zona en comparación con las métricas a nivel diario tanto en la
modelación como en la calibración.

83
 Con la determinación del caudal ambiental que se presenta en la cuenca para
diferentes periodos de tiempo, se concluyó que el caudal disponible para
aprovechamiento presentara un déficit a comparación de la demanda requerida, ya
que no cumple con dicha función y para el periodo comprendido en los años 2041
a 2070, el caudal ambiental disponible para aprovechamiento es mucho menor a los
demás periodos de tiempo, por lo que se estima que se presentaran problemas de
oferta hídrica y afectación al sector económico que requiere de este recurso para su
funcionamiento
 La demanda y la oferta disponible en la cuenca presentan una relación
inversamente proporcional; a medida que aumenta la demanda aproximadamente
4.000 m3/s la oferta disminuye 1.000 m3/s, afirmando el déficit de abastecimiento,
razón por la cual, al determinar el índice de uso de agua se encuentra clasificado
en una categoría critica.
 El resultado de la matriz de vulnerabilidad obtenido del cálculo del índice de
retención hídrica y de uso del agua, se determina que la cuenca del Rio Balsillas no
suple la necesidad de la demanda requerida por la comunidad aledaña que se
abastece, ya que la cuenca presentara un grado de fragilidad Muy Alto.
 Teniendo en cuenta los escenarios de cambio climático, para el RCP 2.6 el dato de
caudal se espera se presente un incremento del 5.2% que corresponde a un
promedio de 1.28 m3/s; para el RCP 4.5 se espera un incremento del 2.1% que
corresponde a un promedio de 1.24 m3/s y por último en el RCP 8.5, los datos de
caudal presentan una media de 1.29 m3/s, que indica un incremento de
aproximadamente del 6%, respecto a los caudales históricos registrados en la
cuenca.
 Con respecto a los valores obtenidos para los índices de variabilidad de oferta
hídrica, se determina que la demanda en la cuenca del río Balsillas no se satisface
completamente, debido a que en los escenarios aplicados se presenta una
probabilidad de fallo que genera un déficit en el suministro de agua.
 Según los resultados de los índices proyectados y contrastados a los valores
observados en la actualidad, se evidencia un aumento en el riesgo al
desabastecimiento dado que la cuenca presenta alta demanda hídrica en
comparación con la oferta de la misma; para los tres periodos de tiempo
correspondientes al año 2015 a 2040, 2041 a 2070 y 2071 a 2100 se presenta un
aumento progresivo con respecto al IRH con una categoría media, para IUA con
una categoría Critica, y para el IVH que es la relación de los índices anteriores una
categoría Baja a Media la cual aumenta a lo largo de los años.
 Al realizar la proyección a futuro para el índice de retención y regulación hídrica
(IRH) se hace evidente para los tres periodos de tiempo evaluados, una capacidad
media de regulación, lo que quiere decir que la cuenca presentara condiciones
estables de caudal, sin embargo, se hace evidente la necesidad de implementar
medidas de regulación y mitigación que generen una reducción de la demanda para
que esta capacidad se mantenga y no sufra un déficit o un cambio que altere en su
totalidad el ecosistema.
 Para la matriz de vulnerabilidad la cual se obtiene del cálculo del índice de retención
hídrica (IVH) y el índice de uso del agua (IUA) se determina que la cuenca del Rio
Balsillas no suple la necesidad de la demanda requerida por la comunidad aledaña
que se abastece, ya que la cuenca presentara un grado de fragilidad Muy Alto

84
 14. RECOMENDACIONES

 Para proyectos que se realicen a futuro, se recomienda contar con otros métodos
donde se haga una estimación de la demanda hídrica más precisa o con más
variables para asegurar un margen de error menor con respecto a las predicciones
no solo para los sectores económicos que están presentes en la zona sino también
para la estimación de la población.
 Se considera importante la participación activa de la comunidad implicada en los
procesos llevados a cabo dentro de la cuenca, con el fin de asegurar un buen uso y
manejo del recurso hídrico para diferentes actividades, creando conciencia de la
importancia del mismo.
 Se recomienda hacer uso de diferentes software que complementen la información
para dar más validez y veracidad de los datos de estimación de la oferta
 La educación como eje central para el fortalecimiento del conocimiento con respecto
a la vulnerabilidad y cambio climático del uso del recurso hídrico para evitar, mitigar
y reducir los impactos negativos que se puedan generar en la Cuenca del Rio
Balsillas.
 Control de las diferentes actividades económicas que se llevan a cabo alrededor de
la cuenca en usos dotacionales, definir áreas de conservación del potencial hídrico
y recuperación de las áreas de regulación natural.
 Implementar estrategias de reúso del agua generada en diferentes procesos en el
sector industrial, con el fin de disminuir el consumo y controlar los consumos que se
generan.
 Se recomienda hacer análisis periódicos para controlar el estado del agua y la
continuidad del caudal ambiental, con el fin de mantener el caudal ecológico para
preservar los ecosistemas dependientes de este

85
 15. BIBLIOGRAFÍA

[¡Error! No se encuentra el origen de la referencia.] Mg. Juana Paredes Díaz,

«Universidad de San Martin de Porres,» 04 2013. [En línea]. Available:
https://www.usmp.edu.pe/publicaciones/boletin/fia/info86/articulos/importanciaAgua.html.
[Último acceso: 15 06 2019].

[¡Error! No se encuentra el origen de la referencia.] PAWD, financiación de la Agencia

Canadiense para el Desarrollo Internacional (CIDA) en el marco del programa, «Planes de
Gestion Integrada del Recurso Hidrico - Manual de Capacitacion y Guia Operacional,» 03
2005. [En línea]. Available:
http://documentacion.ideam.gov.co/openbiblio/bvirtual/021381/PlanesdeGestion.pdf.
[Último acceso: 05 2019].

[¡Error! No se encuentra el origen de la referencia.] Corporación Ambiental Empresarial

(CAEM), filial de la Cámara de Comercio de Bogotá y la Corporación Autónoma Regional
de Cundinamarca (CAR), «CAEM - Corporacion ambiental empresarial,» Seminario
Internacional de Gestión del Riesgo Ambiental, 10 06 2019. [En línea]. Available:
http://www.caem.org.co/reducir-el-impacto-ambiental-y-aumentar-la-productividad-es-
posible-en-colombia-conozca-el-caso-de-la-caem-2/. [Último acceso: 07 2019].
[¡Error! No se encuentra el origen de la referencia.] Kumi Naidoo, Katharine Hayhoe,
Alicia Rivera, Greta Thunberg, Juliana vs United States Government, Chiara Liguori,
«Organizacion Amnistia Intenacional,» Amnesty, 02 2019. [En línea]. Available:
https://www.amnesty.org/es/what-we-do/climate-
change/?utm_source=google&utm_medium=cpc&gclid=EAIaIQobChMIuoGPpoTK4wIVDp
-fCh1zCwbLEAAYASAAEgJJgfD_BwE. [Último acceso: 06 2019].

[¡Error! No se encuentra el origen de la referencia.] Sociedad, Organizacion Ambiente y,

«Sociedad Ambiente y Sociedad,» Justicia Climatica, 08 11 2019. [En línea]. Available:
https://www.ambienteysociedad.org.co/justicia-climatica/. [Último acceso: 08 2019].

[¡Error! No se encuentra el origen de la referencia.] Calle, Efrain Antonio Dominguez,

«Protocolo Para La Modelacion Matematica De Procesos Hidrologicos,» IDEAM, Bogota,
2000.

[¡Error! No se encuentra el origen de la referencia.] Ministerio de Ambiente, vivenda y

desarrollo territorial, «Politica Nacional para la Gestion Nacional del Recurso Hidrico,» 2010.
[En línea]. Available: http://www.minambiente.gov.co/index.php/gestion-integral-del-
recurso-hidrico/direccion-integral-de-recurso-hidrico/politica-nacional-para-la-gestion-
integral-del-recurso-hidrico. [Último acceso: 01 03 2018].
[¡Error! No se encuentra el origen de la referencia.] IDEAM, Cambio climático más
probables para Colombia a lo largo del siglo XXI respecto al clima presente., «Cambio
Climatico mas probables para Colombia a lo largo del siglo XXI respecto al clima presente,»
Grupo de Modelamiento de Tiempo, Clima y Escenarios de Cambio Climático ISubdirección
de Meteorología – IDEAM, Bogota, 2012.

86
[¡Error! No se encuentra el origen de la referencia.] Torres, Blanca Elvira Oviedo,
«GENERACIÓN DE ESCENARIOS DE CAMBIO CLIMÁTICO REGIONALES Y LOCALES
A PARTIR DE MODELOS GLOBALES - GUÍA PARA TOMADORES DE DECISIONES,»
Subdirección de Meteorología, Instituto de Hidrología, Meteorología y Estudios Ambientales
– IDEAM, Bogota, 2010.

[¡Error! No se encuentra el origen de la referencia.] María Claudia García, Andrea

Piñeros Botero, Fabio Andrés Bernal Quiroga, Estefanía Ardila Robles, «Variabilidad
climática, cambio climático y el recurso hídrico en Colombia,» Revista de ingenieria
Universidad de los Andes, IDEAM, Bogota , 2012.

[¡Error! No se encuentra el origen de la referencia.] IDEAM, «Estudio Nacional del Agua

2014 - IDEAM (Capitulo 3),» 05 2015. [En línea]. Available:
http://documentacion.ideam.gov.co/openbiblio/bvirtual/023080/ENA_2014.pdf. [Último
acceso: 01 2019].

[¡Error! No se encuentra el origen de la referencia.] IDEAM, «Estudio Nacional del Agua

2014 - IDEAM (Capitulo 9),» 05 2015. [En línea]. Available:
http://documentacion.ideam.gov.co/openbiblio/bvirtual/023080/ENA_2014.pdf. [Último
acceso: 01 2019].

[¡Error! No se encuentra el origen de la referencia.] Corporacion Autonoma Regional -

CAR, «Elaboracion del Diagnostico, Proespectiva y Formulacion de la Cuenca Hidrografica
del rio Bogota. Subcuenca rio Alto Bogota,» 2006. [En línea]. Available:
https://www.car.gov.co/uploads/files/5ac26116d08ad.pdf. [Último acceso: 02 2019].

[¡Error! No se encuentra el origen de la referencia.] Corporacion Autonoma Regional -

CAR, «Elaboracion del Diagnostico, Prospectiva y Formulacion de la Cuenca Hidrografica
del rio Bogota. Subcuenca del rio Balsillas,» 2006. [En línea]. Available:
https://www.car.gov.co/uploads/files/5ac25b19243a8.pdf. [Último acceso: 02 2019].

[¡Error! No se encuentra el origen de la referencia.] Prof. Alberto C. Palladino , «Grafico

de Caja,» ATENCIÓN PRIMARIA DE SALUD, EPIDEMIOLOGÍA E INFORMÁTICA II,
Bogota, 2011.

[¡Error! No se encuentra el origen de la referencia.] Sokolov, Aleksei Aleksandrovich,

Chapman, Thomas Grandin, «Métodos de cálculo del balance hídrico: guia internacional de
investigación y métodos,» Version espanola : Rafael Heras, Madrid, 1981.

[¡Error! No se encuentra el origen de la referencia.] C. A. R. -. CAR, «Elaboracion del

diagnostico, prospectiva y formulacion de la Cuenca Hidrografica del Rio Bogotá,» 2006.
[En línea]. Available: https://www.car.gov.co/uploads/files/5ac24aeabc81c.pdf. [Último
acceso: 02 2019].
[¡Error! No se encuentra el origen de la referencia.] I. F. Draghici, G. V. Necco, R. W.
Riddaway, J. T. Snow, C. Billard, L. A. Ogallo, «Organizacion meteorologica mundial -
directivas de orientacion Volumen I: meteorologia,» Secretaría de la Organización
Meteorológica Mundial, Ginebra – Suiza, 2001.

87
[¡Error! No se encuentra el origen de la referencia.] IDEAM, «Instituto de Hidrología,
Meteorología y Estudios Ambientales - IDEAM,» Sistema Nacional Ambiental, 2010. [En
línea]. Available: http://www.ideam.gov.co/web/atencion-y-participacion-ciudadana/cambio-
climatico. [Último acceso: 16 02 2019]

[¡Error! No se encuentra el origen de la referencia.] Thomas F. Stocker, Dahe Qin,Gian-

Kasper Plattner, Melinda M.B. Tignor, Simon K. Allen, Judith Boschung, Alexander Nauels,
Yu Xia, Vincent Bex, Pauline M. Midgley , «Resumen para responsables de politicas,
resumen para responsables de politicas "Cambio Climatico 2013 - Bases fisicas",» IPCC-
OMM- PNUMA, Suiza, 2013.
[¡Error! No se encuentra el origen de la referencia.] IPCC, «Afirmaciones principales del
Resumen para responsables de políticas * Cambios observados en el sistema climático,»
2014. [En línea]. Available:
https://archive.ipcc.ch/news_and_events/docs/ar5/ar5_wg1_headlines_es.pdf. [Último
acceso: 2018].

[¡Error! No se encuentra el origen de la referencia.] IDEAM, PNUD, MADS, DNP,

CANCILLERIA, «Tercera comunicacion nacional sobre cambio climatico, herramientas
cientificas para la toma de desiciones,» IDEAM, Bogota, 2015.

[¡Error! No se encuentra el origen de la referencia.] Malte Meinshausen, Keywan Riahi,

Steve Smith, Detlef van Vuuren, «RCP Database,» IPCC, 2010. [En línea]. Available:
https://tntcat.iiasa.ac.at/RcpDb/dsd?Action=htmlpage&page=welcome. [Último acceso: 02
2019].

[¡Error! No se encuentra el origen de la referencia.] D. Q. G.-K. P. M. M. T. S. K. A. J. B.

A. N. Y. X. V. B. P. M. M. Thomas F. Stocker, «IPCC, 2013: Climate Change 2013: The
Physical Science Basis. Contribution of Working Group I to the Fifth Assessment Report of
the Intergovernmental Panel on Climate Change,» Cambridge University Press, New York,
2013.
[¡Error! No se encuentra el origen de la referencia.] IPCC, «Data distribution centre- What
is a GCM,» 2013. [En línea]. Available: https://www.ipcc-
data.org/guidelines/pages/gcm_guide.html. [Último acceso: 2019].

[¡Error! No se encuentra el origen de la referencia.] Rajendra K. Pachauri, Andy

Reisinger, Equipo de redacción principal , «Informe de sintesis. Contribucion de los Grupos
de trabajo I,II y III al Cuarto Informe de evaluacion del Grupo Intergubernamental de
Expertos sobre el Cambio Climatico,» IPCC, Suiza, 2007.

[¡Error! No se encuentra el origen de la referencia.] Georgios Zittis, Prof. Manfred!Lange,

Prof. Jos Lelieveld, Dr. Panos Hadjinicolaou, Prof. Zev Levin, Dr. Christos Giannakopoulos,
«REGIONAL CLIMATE MODELLING IN THE EASTERN MEDITERRANEAN AND THE
MIDDLE EAST: ASPECTS OF SUMMER WARMING AND RELATED FEEDBACKS
UNDER!A CHANGING CLIMATE,» The Cyprus Institute, Nicosia, 2015.

88
[¡Error! No se encuentra el origen de la referencia.] Wilby, Robert & Dawson, C. &
Barrow, E.M., «SDSM – A Decision Support Tool for the Assessment of Regional Climate
Change Impacts. ENVIRON MODELL SOFTW,» Elsevier, Canada, 2002.

[¡Error! No se encuentra el origen de la referencia.] United States Agency for

Intewrnational Development USAID, «The Center for International Earth Science
Information Network (CIESIN) - A review of downscaling methods for climate change
projections,» 03 2019. [En línea]. Available:
http://www.ciesin.org/documents/Downscaling_CLEARED_000.pdf. [Último acceso: 09
2017].

[¡Error! No se encuentra el origen de la referencia.] Jorge A. Amador, Eric J. Alfaro, «The

Remote Sensing and Geographic Information Science (RSGIS) - Métodos de reducción de
escala: aplicaciones al tiempo, clima,variabilidad climática y cambio climático,» 17 02 2009.
[En línea]. Available:
http://www.rsgis.envs.lsu.edu/crn/docs/AmadoryAlfaroRevibec2009.pdf. [Último acceso: 02
2019].

[¡Error! No se encuentra el origen de la referencia.] Freddy Duarte, Gerald Corzo,

Germán Santos, Oscar Hernández, «Chaotic Statistical Downscaling (CSD): Application
and Comparison in the Bogotá River Basin,» Goffredo La Loggia , Gabriele Freni , Valeria
Puleo y Mauro De Marchis , Bogota, 2018.

[¡Error! No se encuentra el origen de la referencia.] Duarte Prieto, Freddy Santiago ,

«Técnica de reducción de escala estadística basada en la teoría del caos: aplicación y
desempeño en la cuenca del río Bogotá,» Escuela Colombiana de Ingeniería Julio Garavito,
Bogota, 2017.

[¡Error! No se encuentra el origen de la referencia.] Bellie Sivakumar, Ronny Berndtsson,

Advances in Data-Based Approaches for Hydrologic Modeling and Forecasting, California,
Australia: Bellie Sivakumar (The University of New South Wales, Sydney, Australia &
University of California, Davis, USA) and Ronny Berndtsson (Lund University, Sweden),
2010.

[¡Error! No se encuentra el origen de la referencia.] Sharma Ashish and Mehrotre

Rajeshwar, «An information theoretic alternative to model a natural system using
observational information alone,» Water Resources Research, Estados Unidos, 2014.

[¡Error! No se encuentra el origen de la referencia.] M. Sano and Y. Sawada,

«Measurement of Lyapunov Spectrum from a Chaotic Time Series Vol.55,» Physical Review
Letters, Estados Unidos, 1985.

[¡Error! No se encuentra el origen de la referencia.] Blanca Elvira Oviedo Torres – Gloria

León Aristizábal, «Instituto de Hidrología, Meteorología y Estudios Ambientales (IDEAM) -
GUÍA DE PROCEDIMIENTO PARA LA GENERACIÓN DE ESCENARIOS DE CAMBIO
CLIMÁTICO REGIONAL Y LOCAL A PARTIR DE LOS MODELOS GLOBALES,» 2010. [En
línea]. Available:
http://www.ideam.gov.co/documents/21021/21138/Gu%C3%ADa+Escenarios+Cambio+Cli
matico.pdf/72eae24f-04ea-4ce2-9a4b-e551559c48fc. [Último acceso: 13 04 2019].

89
[¡Error! No se encuentra el origen de la referencia.] Blanca Elvira Oviedo Torres,
«Instituto de Hidrología, Meteorología y Estudios Ambientales (IDEAM)- GENERACIÓN DE
ESCENARIOS DE CAMBIO CLIMÁTICO REGIONALES Y LOCALES A PARTIR DE
MODELOS GLOBALES - GUÍA PARA TOMADORES DE DECISIONES,» 2010. [En línea].
Available:
http://www.ideam.gov.co/documents/21021/21138/Gu%C3%ADa+Escenarios+para+Toma
dores+de+Decisiones.pdf/fa7abe38-43cc-49c8-96a3-f2b5c24ecce3. [Último acceso: 06 03
2019].

[¡Error! No se encuentra el origen de la referencia.8] Clarke, J. Edmonds, H. Jacoby, H.

Pitcher, J. Reilly, R. Richels, «Scenarios of Greenhouse Gas Emissions and Atmospheric
Concentrations,» Editor Loel Kathmann, Technical Advisor David Dokken, Graphic
Production DesignConcept, Estados Unidos, 2007.

[¡Error! No se encuentra el origen de la referencia.] WCRP CORDEX IPOC, «SMHI -

CORDEX — a global collaborative initiative,» 2016. [En línea]. Available:
https://www.smhi.se/polopoly_fs/1.106822!/CORDEX_160707.pdf. [Último acceso: 15 04
2019].

[¡Error! No se encuentra el origen de la referencia.] IDEAM - Instituto de Hidrología,

Meteorología y Estudios Ambientales, «IDEAM - Modelacion Hidrologica,» Ministerio de
Ambiente, 2014. [En línea]. Available: http://www.ideam.gov.co/web/agua/modelacion-
hidrologica. [Último acceso: 05 04 2019].

[¡Error! No se encuentra el origen de la referencia.] Saul Montoya and filed under

Hidrología, «Gestion sostenible del Agua (gidahatari) - Caracteristicas y capacidades del
modelo hidrologico desitribuido TETIS,» gidahatari, 29 05 2016. [En línea]. Available:
http://gidahatari.com/ih-es/caracteristicas-y-capacidades-del-modelo-hidrologico-
distribuido-tetis. [Último acceso: 04 2019].

[¡Error! No se encuentra el origen de la referencia.] Fekadu Moreda, Ph.D. Fernando

Miralles-Wilhelm Raúl Muñoz Castillo, «Hydro-BID: Un sistema integrado para la simulacion
de impactos del Cambio Climatico sobre los recursos hidricos. Departamento de
infraestructura y medio ambiente, division de agua y saneamiento. Rio de Janeiro: Banco
Interamericano de Desarrollo,» Banco Interamericano de Desarrollo, Estados Unidos, 2014.

[¡Error! No se encuentra el origen de la referencia.] Fekadu Moreda, Ph.D., Fernando

Miralles-Wilhelm, Raúl Muñoz Castillo , «Hydro-BID: Un Sistema Integrado para la
Simulación de Impactos del Cambio Climático sobre los Recursos Hídricos. Parte 2,» Banco
Interamericano de Desarrollo, Estados Unidos, 2014.

[¡Error! No se encuentra el origen de la referencia.] Alan Wyatt Fekadu Moreda Eugene

Brantly Fernando Miralles-Wilhelm Raúl Muñoz Castillo , «Caso de estudio de Hydro-BID
N°1: Modelo de gestión del recurso hídrico en la cuenca del Río Grande en Argentina,»
Banco Interamericano de Desarrollo, Estados Unidos, 2014.

[¡Error! No se encuentra el origen de la referencia.] Banco Interamericano de Desarrollo

- Hydro-BID, «Hydro-BID,» Banco Interamericano de Desarrollo, 03 05 2016. [En línea].
Available: http://sp.hydrobidlac.org/sobre-hydrobid. [Último acceso: 06 05 2019].

90
[¡Error! No se encuentra el origen de la referencia.] Mauro Nalesso, Pedro Coli, «Guía
paso a paso: Manual de Hydro-BID,» Banco Interamericano de Desarrollo, Estados Unidos,
2017.

[¡Error! No se encuentra el origen de la referencia.] Banco Central de Colombia, «Banco

de la Republica - Producto Interno Bruto (PIB),» banco central de Colombia, 2015. [En
línea]. Available: https://www.banrep.gov.co/es/estadisticas/producto-interno-bruto-pib.
[Último acceso: 05 2019].

[¡Error! No se encuentra el origen de la referencia.] Fabio Andrés Bernal Quiroga, Jenny

Paola Santander Durán, «Instituto de Hidrología, Meteorología y Estudios Ambientales –
Formato Común de Hoja Metodológica de Indicadores Ambientales "Índice de retención y
regulación hídrica",» Ministerio de Ambiente y Desarrollo Sostenible, 2014. [En línea].
Available:
http://www.ideam.gov.co/documents/11769/646961/3.04+HM+Indice+Regulacion+hidrica.
pdf/d980a72f-9fa0-4a70-b626-8e99a743b465. [Último acceso: 05 2019].

16. ANEXOS

ANEXO A. Comparación de datos históricos y datos futuros para el modelo CORDEX

MPI-ESM-MR para la variable de precipitación

 RCP 2.6

91
 RCP 4.5

92
 RCP 8.5

93
94
ANEXO B. Comparación de datos históricos y datos futuros para el modelo GCM MPI-
ESM-MR para la variable de temperatura

 RCP 2.6

95
 RCP 4.5

96
 RCP 8.5

97
 ANEXO C. Puntos de demanda en la cuenca por municipios y sectores según registros de la CAR

 Municipio de Bojaca

Actividad Uso Uso Uso Uso

No. Expediente Vereda Predio Clase Total
Económica Doméstico Pecuario Agricola Industrial
Abastecimiento
San Cría de ganado
1 1312 Guarenta doméstico, 0.01 0.01 0 0 0.02
Antonio bovino y bufalino.
Abrevaderos
Abastecimiento Explotación mixta
San
2 17520 La Suiza doméstico, (agrícola y 0.01 0.01 0 0 0.02
Antonio
Abrevaderos pecuaria).
Captación,
San
3 26600 Las Luisitas Industrial depuración y 0.039 0 0 0.499 0.538
Antonio
distribución de agua
Abastecimiento
Barro Hacienda el Actividad mixta
4 32728 doméstico, 0.03 0.03 0 0 0.06
Blanco Triunfo (agricola y pecuaria)
Abrevaderos
Actividad mixta
5 35873 Chilcal Finca Pilitas Abrevaderos 0 0.01 0 0 0.01
(agricola y pecuaria)
Abastecimiento
San Finca la Actividad pecuaria
6 37391 doméstico, 0.01 0.1 0 0 0.11
Antonio Libertad no especializada
Abrevaderos
Actividades de
servicios, agricolas
Abastecimiento
Finca Santa y ganaderos
7 40479 Chilcal doméstico, 0.01 0.4 0 0 0.41
Cecilia exceptos las
Abrevaderos
actividades
veterinarias
Riego y
San Cría de ganado
8 40516 Buenavista silvicultura, 0 0.01 0.05 0 0.06
Antonio bovino y bufalino.
Abrevaderos
Abastecimiento
San Captación,
9 41727 la Libertad doméstico, 1.38 0.57 3.37 0 5.32
Antonio tratamiento y
Riego y

98
 silvicultura, distribución de
Abrevaderos agua.
Abastecimiento
Elaboracion de
San doméstico,
10 42700 El Mirador otros productos 0.01 0.01 0 0 0.02
Antonio Abrevaderos,
alimenticios
Industrial
Riego y Explotación mixta
Finca el
11 43843 Chilcal silvicultura, (agrícola y 0.01 0.01 0.02 0 0.04
Encenillo
Abrevaderos pecuaria).
San Lote San Abastecimiento
12 45009 0.01 0 0 0 0.01
Antonio Antonio doméstico
Finca la
Roble Abastecimiento Actividad mixta
13 46342 Atlantida 0.02 0 0 0 0.02
Hueco doméstico (agricola y pecuaria)
S.A.S.
Finca la Riego y
14 52360 Chilcal 0 0 0.25 0 0.25
Palma silvicultura
La Laguna, Riego y Cría de ganado
15 59655 Bobace 0 0 0 0 2.4
lote silvicultura bovino y bufalino.
Finca Riego y
Cría de ganado
16 64324 Chilcal Calichana y silvicultura, 0 2 37 0 39
bovino y bufalino.
Calucata Abrevaderos

 Municipio de El Rosal

Actividad Uso Uso Uso Uso

No. Expediente Vereda Predio Clase Total
Económica Doméstico Pecuario Agrícola Industrial

Actividades no
diferenciadas de los
Abastecimiento hogares
1 2131 Cruz Verde El Jardin doméstico, individuales como 0.01 0.01 0 0 0.02
Abrevaderos productores de
bienes para uso
propio.
Abastecimiento Cría de ganado
2 3511 Buenavista Santa Inés 0 0.01 0.13 0 0.14
doméstico bovino y bufalino.

99
 El
Recuerdo Riego y Explotación mixta
3 3572 La Piñuela lote a , el silvicultura, (agrícola y 0 0.02 1.12 0 1.14
recuerdo Abrevaderos pecuaria).
lote b
Abastecimiento
doméstico, Cultivo de frutas
La
4 12241 Buenavista Riego y tropicales y 0 0.01 0.01 0 0.02
Esperanza
silvicultura, subtropicales.
Abrevaderos
Cultivo de especias
Riego y
Santa y de plantas
5 19544 Santa Rita silvicultura, 0 0.01 0.52 0 0.53
Bárbara aromáticas y
Abrevaderos
medicinales.
Cultivo de especias
Riego y
Santa y de plantas
6 19544 Santa Rita silvicultura, 0 0.01 0.52 0 0.53
Bárbara aromáticas y
Abrevaderos
medicinales.
Cultivo de
Puerta de Puente el Riego y
7 24374 hortalizas, raíces y 0.01 0 0.27 0 0.28
Cuero Rosal silvicultura
tubérculos.
Abastecimiento
Cultivo de frutas
doméstico,
8 24983 Buenavista El Tibar tropicales y 0.01 0 0.09 0 0.1
Riego y
subtropicales.
silvicultura
Abastecimiento
Cultivo de frutas
doméstico,
9 24983 Buenavista El Tibar tropicales y 0.01 0 0.09 0 0.1
Riego y
subtropicales.
silvicultura
Santa Abastecimiento Actividad mixta
10 24985 Buenavista María, Villa doméstico, (agricola y 0 0 0.05 0 0.05
Real Abrevaderos pecuaria)
Abastecimiento Actividades de
11 24996 Buenavista Salitre doméstico, apoyo a la 0.01 0.01 0 0 0.02
Abrevaderos ganadería.

100
 Produccion
Riego y especializada de
12 25064 El Rodeo La Fortuna silvicultura, frutas, nueces, 0 0 0 0 0.24
Abrevaderos plantas bebestibles
y especias
Riego y Cultivo de frutas
13 25336 Buenavista Villa María silvicultura, tropicales y 0 0.01 0.13 0 0.14
Abrevaderos subtropicales.
Explotación mixta
Villa Riego y
14 25398 Buenavista (agrícola y 0 0 0.01 0 0.01
Tatiana silvicultura
pecuaria).
Abastecimiento Cria especializada
doméstico, de ovejas, cabras,
Finca la
15 25400 El Rodeo Riego y caballos, asnos, 0.009 0.025 0.26 0 0.294
Hamaca
silvicultura, mulas y
Abrevaderos burdeganos

LA
Riego y Actividad pecuaria
16 25858 PORQUER El Encanto 0.05 0.05 0.15 0 0.25
silvicultura no especializada
A

Abastecimiento
Cultivo de frutas
San doméstico,
17 27379 Buenavista tropicales y 0.009 0 0.18 0 0.189
Gregorio Riego y
subtropicales.
silvicultura
Riego y
Los Actividad pecuaria
18 29979 El Rodeo silvicultura, 0.0002 0 0 0 0.0002
Tamarindos no especializada
Abrevaderos
Abastecimiento
doméstico, Actividad mixta
19 30060 El Caucho La Maria Riego y (agricola y 0.01 0.01 0.003 0 0.023
silvicultura, pecuaria)
Abrevaderos
Manoa San Actividad mixta
Riego y
20 30661 Tibagota Rafael , (agricola y 0 0 4.41 0 4.41
silvicultura
Vicenza pecuaria)

101
 Pozo Azul
Riego y Cultivo de flor de
21 33063 Buenavista Cruz Verde 0 0 1.19 0 1.19
silvicultura corte.
y las Monas
Actividad mixta
Finca las Riego y
22 33208 Buenavista (agricola y 0 0 0.61 0 0.61
Monas silvicultura
pecuaria)
San Riego y Cria especializada
23 33592 El Tabor 0 0 0.12 0 0.12
Antonio silvicultura de ganado vacuno
Abastecimiento
doméstico, Actividad mixta
Finca la
24 33879 El Caucho Riego y (agricola y 0.01 0.01 0.12 0 0.14
Alcira
silvicultura, pecuaria)
Abrevaderos
Abastecimiento
doméstico,
Cría de aves de
25 34174 Cruz Verde El Triunfo Riego y 0.01 0.01 0.03 0 0.05
corral.
silvicultura,
Abrevaderos
El
Recuerdo ,
el Retiro , la Riego y Actividad mixta
26 35260 Tibagota Esmeralda, silvicultura, (agricola y 0 0.05 1.1 0 1.15
la Piñuela , Abrevaderos pecuaria)
las
Margaritas
Abastecimiento
Santa doméstico, Cultivo de flor de
27 36026 Indiana II 0.2 0 0.93 0 1.13
Bárbara Riego y corte.
silvicultura
Abastecimiento
Lote No.3
doméstico, Explotación mixta
Finca el
28 36568 Buenavista Riego y (agrícola y 0.02 0.01 0.16 0 0.19
Verde o la
silvicultura, pecuaria).
Cartuja
Abrevaderos
Abastecimiento Actividades de
Vereda
29 37264 Buenavista doméstico, apoyo a la 0.01 0.01 0.02 0 0.04
Buenavista
Riego y ganadería.

102
 silvicultura,
Abrevaderos
Abastecimiento
San Gabriel Cria especializada
30 37265 El Rodeo doméstico, 0.02 0.01 0 0 0.03
Vacolandia de ganado vacuno
Abrevaderos
Carrera 62
Santa Riego y Cría de ganado
31 40062 # 129 b 25 0 0 0 0 0.39
Bárbara silvicultura bovino y bufalino.
apto 201
Abastecimiento
Cultivo de especias
Mitad San doméstico,
y de plantas
32 40454 Buenavista Antonio y/o Riego y 0.01 0 0.03 0 0.04
aromáticas y
Zerta silvicultura,
medicinales.
Abrevaderos
Encenillo Produccion agricola
Riego y
33 41397 El Rodeo Camilo y NCP en unidades 0 0 0.04 0 0.04
silvicultura
Chital especializadas
Riego y Explotación mixta
Puerta de Hacienda
34 42994 silvicultura, (agrícola y 0 0.06 0.42 0 0.48
Cuero Tatambó
Abrevaderos pecuaria).
Riego y Actividad mixta
35 43187 El Caucho Tasue silvicultura, (agricola y 0 0.01 0.01 0 0.02
Abrevaderos pecuaria)
Santa Maria
de los
Angeles,
Santa Maria Riego y Explotación mixta
Hondura
36 44334 de los silvicultura, (agrícola y 0.01 0.01 0.01 0 0.03
Tibagota
Angeles, Abrevaderos pecuaria).
Santa Maria
de los
Angeles
Explotación mixta
Monroy I y Riego y
37 49348 Cruz Verde (agrícola y 0 0 0.06 0 0.06
II silvicultura
pecuaria).

103
 Lote 1 San
Carlos, lote
Abastecimiento
2 San
doméstico,
Alberto, lote
38 56270 Cruz Verde Riego y 0.001 0.001 0.0009 0 0.0029
3 el Portal
silvicultura,
del Rosal,
Abrevaderos
lote 4 Santa
Isabel
Abastecimiento
Explotación mixta
Lote 2 San doméstico,
39 60033 Cruz Verde (agrícola y 0.001 0.001 0.0009 0 0.0029
Alberto Riego y
pecuaria).
silvicultura
Riego y Explotación mixta
Lote 3 el
40 60034 El Rodeo silvicultura, (agrícola y
Portal
Abrevaderos pecuaria).
Actividades no
Abastecimiento diferenciadas de los
doméstico, hogares
Lote 4 Sta
41 60035 Cruz Verde Riego y individuales como 0.001 0.001 0.0009 0 0.0029
Isabel
silvicultura, productores de
Abrevaderos bienes para uso
propio.
Abastecimiento
Hondura Cría de ganado
42 60191 La Alhaja doméstico, 0.0017 0.0054 0 0 0.0071
Chingafrio bovino y bufalino.
Abrevaderos
Explotación mixta
Los Riego y
43 62517 El Rodeo (agrícola y 0 0 0.98 0 0.98
Angeles silvicultura
pecuaria).
Cultivo de frutas
Riego y
44 62564 Buenavista Rossana tropicales y 0 0 0.1 0 0.1
silvicultura
subtropicales.
Abastecimiento
Hondura doméstico,
45 63437 El Triángulo Ganadería. 0.008 0.243 0 0 0.251
Chingafrio Riego y
silvicultura

104
 Abastecimiento
Hondura Lote 1
46 63918 doméstico, Ganadería. 0 0.0253 0 0 0.0253
Chingafrio Predio rural
Abrevaderos
Cría de ganado
47 66598 El Rodeo Villa Julia Abrevaderos 0 0.01 0 0 0.01
bovino y bufalino.

 Municipio de Madrid

Actividad Uso Uso Uso Uso

No. Expediente Vereda Predio Clase Total
Económica Doméstico Pecuario Agrícola Industrial
1 880 Boyero El Imperio Riego y Cría de ganado 0 0.02 0.38 0 0.4
silvicultura, bovino y bufalino.
Abrevaderos
2 25415 Laguna El Colegio Riego y Cría de ganado 0 0.05 0.43 0 0.48
Larga silvicultura, bovino y bufalino.
Abrevaderos
3 27098 CORTIJO Planta Principal Abastecimiento Captación, 37.43 0 0 0 37.43
de Tratamiento doméstico depuración y
distribución de
agua
4 28624 La Finca Coralina Abastecimiento Explotación mixta 0.02 0 0 0 0.02
Cuesta doméstico (agrícola y
pecuaria).
5 30208 Laguna Barley del Rio Riego y Explotación mixta 0 0 1.83 0 1.83
Larga silvicultura (agrícola y
pecuaria).
6 30226 Los El Darien Riego y Actividad mixta 0 0.04 1.46 0 1.5
Arboles silvicultura, (agricola y
Abrevaderos pecuaria)
7 33870 Puente Segregado Industrial Elaboración de 0 0 0 2 2
Piedra Bosque II otros productos
alimenticios n.c.p.

105
8 33984 Los Finca la Riego y Otros cultivos 0 0 0.54 0 0.54
Arboles Yerbabuena silvicultura permanentes n.c.p.
lote N° 2
9 35125 Santa El Diamante Riego y Actividad mixta 0 0.14 0.6 0 0.74
Cruz silvicultura, (agricola y
Abrevaderos pecuaria)
10 37601 Los Finca Altagracia Riego y Cría de ganado 0 0 0.48 0 0.48
Arboles silvicultura bovino y bufalino.
11 43462 Los El Bosque Riego y Cultivos agrícolas 0 0 0.13 0 0.13
Arboles silvicultura transitorios.
12 46064 El Corzo Las Monjitas Riego y Explotación mixta 0 0 0.17 0 0.17
silvicultura (agrícola y
pecuaria).
13 49856 Santa Los Sauces Abastecimiento Explotación mixta 0.05 0.02 0.1 0 0.17
Cruz doméstico, (agrícola y
Riego y pecuaria).
silvicultura,
Abrevaderos
14 54221 Potrero Lote 8 Riego y Cultivo de 0 0 240 0 240
Grande Hacienda Las silvicultura hortalizas, raíces y
Mercedez tubérculos.
15 54797 El Corzo El Porvenir, Las Riego y Otros cultivos 0 0 0.15 0 0.15
Delicias, lote silvicultura transitorios n.c.p.
Nº8
16 59163 Santa Lote Riego y Cultivos agrícolas 0 0 0.13 0 0.13
Cruz silvicultura transitorios.
17 59165 Santa San Ignacio del Riego y Cultivos agrícolas 0 0 0.0013 0 0.0013
Cruz Rio silvicultura transitorios.
18 59166 Santa El Triunfo Riego y Explotación mixta 0 0.04 0.13 0 0.17
Cruz silvicultura, (agrícola y
Abrevaderos pecuaria).
19 62519 BARRO Los Napoles Riego y Explotación mixta 0 0.15 0.83 0 0.98
BLANCO silvicultura, (agrícola y
Abrevaderos pecuaria).
20 62907 Los Yerbabuena Riego y Actividades de 0 0 4.08 0 4.08
Arboles Lote 2 silvicultura apoyo a la
agricultura.

106
  Municipio de Mosquera

Actividad Uso Uso Uso Uso

No. Expediente Vereda Predio Clase Total
Económica Doméstico Pecuario Agrícola Industrial
Riego y Otros cultivos
1 49631 San Jorge Chamberi 0 0 0.0005 0 0.0005
silvicultura permanentes n.c.p.
Extracción de
Riego y piedra, arena,
2 52369 Balsillas Andes 1 0 0.03 0.31 0 0.34
silvicultura arcillas comunes,
yeso y anhidrita.
Riego y Explotación mixta
Lote San
3 53719 San José silvicultura, (agrícola y 0 0.0002 0.0027 0 0.0029
Pablo
Abrevaderos pecuaria).
Lote # 1
alfandoque,
lote b la Riego y Cría de caballos y
4 54212 San José 0 0.134 0 0 2
esperanza, silvicultura otros equinos.
lote de terreno
holanda alta
Riego y
5 54447 San José Las Goteras 0 0 1.21 0 1.21
silvicultura
Educación
Riego y
6 56908 San José Sena secundaria y de 0 0 2.06 0 2.06
silvicultura
formación laboral.
San Lote 1 (la Riego y Cultivos agrícolas
7 59367 0 0 0.44 0 0.44
Francisco valentina) silvicultura transitorios.
8 59368 San José El Trebol Abrevaderos Ganadería. 0 0.103 0.104 0 0.207
Riego y
San Los Locales
9 59369 silvicultura, Ganadería. 0 0.09 0.3 0 0.39
Francisco Yerbabuena
Abrevaderos
Riego y
San Los Locales
10 59370 silvicultura, Ganadería. 0 0.0009 0.003 0 0.0039
Francisco Yerbabuena
Abrevaderos

107
 Riego y
San
11 59376 El Carretonal silvicultura, 0 0.184 0.296 0 0.48
Francisco
Abrevaderos
Camino la
Alameda, la
San Riego y Cría de ganado
12 59726 Isabela I, la 0 0.252 3.648 0 3.9
Francisco silvicultura bovino y bufalino.
Isabela II, Lote
Casa
El Cigarrillo, el
San Riego y Cría de ganado
13 59734 Recreo, la 0 0.063 0 0 3.855
Francisco silvicultura bovino y bufalino.
Isabela
Riego y Explotación mixta
San
14 60219 La Herradurita silvicultura, (agrícola y 0 0.1 1.1 0 1.2
Francisco
Abrevaderos pecuaria).
Riego y Actividades de
San
15 60239 La Herradura silvicultura, apoyo a la 0 0.1 2.2 0 2.3
Francisco
Abrevaderos ganadería.
Riego y Cultivo de flor de
16 60640 San José Lote 2 0 0 0 0 3.06
silvicultura corte.
Captación,
Riego y tratamiento y
17 60645 San José La Playa 0 0 3.96 0 3.96
silvicultura distribución de
agua.
Cultivo de
Riego y
18 60878 San José La Pradera hortalizas, raíces y 0 0 0.5 0 0.5
silvicultura
tubérculos.
Establecimientos
San Abastecimient que combinan
19 61049 La Escuela 0.16 0 0 0 0.16
Francisco o doméstico diferentes niveles
de educación.
San Riego y Cría de ganado
20 61229 Venecia 0 0.08 2.7 0 2.78
Francisco silvicultura bovino y bufalino.
San Riego y Cría de ganado
21 61230 El Sociego 0 0.08 3.34 0 3.42
Francisco silvicultura bovino y bufalino.
Explotación mixta
San Abastecimient
22 61234 Normandia (agrícola y 0 0.18 2.4 0 2.58
Francisco o doméstico,
pecuaria).

108
 Riego y
silvicultura
Abastecimient
San o doméstico, Cría de ganado
23 61236 El Establo 0 0.09 3.21 0 3.3
Francisco Riego y bovino y bufalino.
silvicultura
Abastecimient
San o doméstico, Cría de ganado
24 61237 El Pensil 0 0.19 4.24 0 4.43
Francisco Riego y bovino y bufalino.
silvicultura
Explotación mixta
San Riego y
25 61238 Los Sauces (agrícola y 0 0.18 0 0 2.34
Francisco silvicultura
pecuaria).
Abastecimient
San o doméstico, Cría de ganado
26 61239 El Establito 0.05 0.1 3.86 0 4.01
Francisco Riego y bovino y bufalino.
silvicultura
Los Cristales,
Riego y Actividades de
San Maria
27 65335 silvicultura, apoyo a la 0 0.104 0.619 0 0.723
Francisco Fernanda,
Abrevaderos agricultura.
Pubenza

 Municipio de Subachoque

Actividad Uso Uso Uso Uso

No. Expediente Vereda Predio Clase Total
Económica Doméstico Pecuario Agrícola Industrial
1 271 Santuario Riego y Explotación mixta 0 0.01 0.57 0 0.58
Las Manas silvicultura (agrícola y
pecuaria).
2 430 Casa Riego y Produccion agricola 0.02 0.04 1.01 0 1.07
Blanca Iscata silvicultura, en unidades no
Abrevaderos especializadas

109
3 548 Rincón Riego y Cría de ganado 0 0.06 1.16 0 1.22
El Ttriunfo, La
Santo silvicultura, bovino y bufalino.
Esmeralda
Abrevaderos
4 1175 Páramo Abastecimient Actividad mixta 0.01 0.01 0 0 0.02
o doméstico, (agricola y
Abrevaderos pecuaria)
5 1246 Galdamez Abrevaderos Explotación mixta 0 0 0.11 0 0.11
Alejandria (agrícola y
pecuaria).
6 1538 Galdame Abastecimient Cria especializada 0.01 0.01 0 0 0.02
El Recuerdo
z o doméstico de ganado vacuno
7 1779 Canica Abastecimient Actividad mixta 0.02 0.01 0.15 0 0.18
o doméstico, (agricola y
Bocachica Riego y pecuaria)
silvicultura,
Abrevaderos
8 2943 Galdame Abrevaderos Explotación mixta 0.01 0.04 0.54 0 0.59
z Las Celicias 3 (agrícola y
pecuaria).
9 3219 Centro Abastecimient Explotación mixta 0.03 0.03 1.41 0 1.47
Suatia, o doméstico, (agrícola y
Tacarigua, Riego y pecuaria).
Zyruma silvicultura,
Abrevaderos
10 3496 Cascajal Riego y Actividad pecuaria 0 0.04 0.07 0 0.11
Palos Verdes silvicultura, no especializada
Abrevaderos
11 3623 Los Abastecimient Actividades no 0.06 0.03 0 0 0.09
Llanitos o doméstico, diferenciadas de
Abrevaderos los hogares
El Rodeo individuales como
productores de
bienes para uso
propio.
12 3874 Cascajal Abastecimient Actividad mixta 2.46 0 0 0 2.46
Had El Vergel o doméstico (agricola y
pecuaria)

110
13 3928 El Tobal Abastecimient Propagación de 0 0 0.0002 0 0.0002
o doméstico plantas (actividades
Arboreto de los viveros,
excepto viveros
forestales).
14 5329 Galdamez Riego y Actividad mixta 0 0.032 1.588 0 1.62
Hacienda El
silvicultura, (agricola y
Hato
Abrevaderos pecuaria)
15 5653 La Abastecimient Captacion, 2.43 0 0 0 2.43
Yeguera San Cayetano o doméstico depuracion y
lote 1a distribucion de
agua
16 7917 Rincón Abastecimient Actividad mixta 0.01 0 0.12 0 0.13
Santo Chamomile o doméstico (agricola y
pecuaria)
17 10186 Altania Abastecimient Otros tipos de 1 0 0 0 1
El Paramo
o doméstico alojamiento n.c.p.
18 10186 Altania Abastecimient Otros tipos de 1 0 0 0 1
El Paramo
o doméstico alojamiento n.c.p.
19 10195 Rincón Abastecimient Actividades de 0.01 0.01 0.1 0 0.12
Santo o doméstico, servicios, agricolas
Riego y y ganaderos
Finca Versalles
silvicultura, exceptos las
Abrevaderos actividades
veterinarias
20 12296 Canica Riego y Explotación mixta 0 0.01 0.19 0 0.2
El Curubo silvicultura, (agrícola y
Abrevaderos pecuaria).
21 12336 Casa Asociación de Abastecimient Captación, 0.39 0 0 0 0.39
Blanca usuarios del o doméstico tratamiento y
acueducto de distribución de
la quebrada agua.
chiqueros
22 13784 Canica Riego y Explotación mixta 0 0.01 0.95 0 0.96
San Pablo silvicultura, (agrícola y
Abrevaderos pecuaria).

111
23 16174 Rincón Abastecimient Explotación mixta 1 4 15 0 20
Santo o doméstico, (agrícola y
Bellavista Riego y pecuaria).
silvicultura,
Abrevaderos
24 16410 Guamal Riego y Actividad mixta 0 0.01 0.11 0 0.12
silvicultura, (agricola y
Las Delicias Abrevaderos pecuaria)

25 16412 Guamal Abastecimient Actividad mixta 0.01 0.01 0.03 0 0.05

Finca San o doméstico, (agricola y
Miguel, san Riego y pecuaria)
miguel silvicultura,
Abrevaderos
26 16418 Guamal Abastecimient Actividad mixta 0.1 0.1 0.6 0 0.8
o doméstico, (agricola y
Bellavista Riego y pecuaria)
silvicultura,
Abrevaderos
27 16419 Guamal Abastecimient Actividad mixta 0.01 0.01 0.01 0 0.03
o doméstico, (agricola y
El Valle Riego y pecuaria)
silvicultura,
Abrevaderos
28 16425 Guamal Abastecimient Explotación mixta 0.01 0.01 0.01 0 0.03
o doméstico, (agrícola y
San Pablo Riego y pecuaria).
silvicultura,
Abrevaderos
29 16426 Guamal Abastecimient Actividad pecuaria 0.01 0.01 0.24 0 0.26
o doméstico, no especializada
La Veguita Riego y
silvicultura,
Abrevaderos
30 18420 Galdame Riego y Cría de ganado 0 0.01 0.01 0 0.02
z El Encanto silvicultura, bovino y bufalino.
Abrevaderos

112
31 19220 Canica Abrevaderos Explotación mixta 0 0.02 0.33 0 0.35
El Cerezo II, La
(agrícola y
Esperanza
pecuaria).
32 19328 Centro Abastecimient Captacion, 1.1 0.07 0 0 1.17
o doméstico, depuracion y
La Ilusión
Abrevaderos distribucion de
agua
33 19457 Guamal Abastecimient Cria especializada 0.01 0.01 0.013 0 0.033
o doméstico, de ganado vacuno
Aconcagua La
Riego y
Cabrera
silvicultura,
Abrevaderos
34 19457 Guamal Abastecimient Cria especializada 0.01 0.01 0.013 0 0.033
o doméstico, de ganado vacuno
Aconcagua La
Riego y
Cabrera
silvicultura,
Abrevaderos
35 20480 Santuario Abastecimient Explotación mixta 0 0 0.03 0 0.03
La Maria 3 o doméstico (agrícola y
pecuaria).
36 21159 La Abrevaderos Actividad pecuaria 0 0.01 0 0 0.01
Tocaima
Cuesta no especializada
37 24755 Galdamez Riego y Actividades de 0 0 0 0 0
Villa Tatiana silvicultura, apoyo a la
Abrevaderos ganadería.
38 25337 Chinga Abrevaderos Actividad pecuaria 0 0.02 0 0 0.02
Los Arrayenes
Frío no especializada
39 25399 Casa Abastecimient Captacion, 1.23 0 0 0 1.23
Blanca o doméstico depuracion y
San Carlos
distribucion de
agua
40 25434 Guamal Abastecimient Captación, 0.59 0.04 0 0 0.63
o doméstico, tratamiento y
Asoguamal
Abrevaderos distribución de
agua.
41 25591 El Tobal Abastecimient Ganadería. 0.01 0.01 0.06 0 0.08
Casas o doméstico,

113
 Riego y
silvicultura
42 25854 Galdamez Abastecimient Actividades no 0.005 0 0 0 0.005
o doméstico diferenciadas de
los hogares
El Cucharo individuales como
productores de
bienes para uso
propio.
43 25854 Galdamez Abastecimient Actividades no 0.005 0 0 0 0.005
o doméstico diferenciadas de
los hogares
El Cucharo individuales como
productores de
bienes para uso
propio.
44 26523 Canica Riego y Explotación mixta 0 0.03 0.05 0 0.08
Baja Finca Deja-vu silvicultura, (agrícola y
Abrevaderos pecuaria).
45 26961 Canica Abastecimient Cría de ovejas y 0.011 0.039 0.69 0 0.74
Baja o doméstico, cabras.
La Danielita Riego y
silvicultura,
Abrevaderos
46 27097 Centro El Bosque de Riego y Otros cultivos 0 0 0.033 0 0.033
kirpal silvicultura transitorios n.c.p.
47 27153 Rincón Abastecimient Cría de ganado 0.02 0.01 0.22 0 0.25
Santo o doméstico, bovino y bufalino.
El Arrayan Riego y
silvicultura,
Abrevaderos
48 27347 Rincón Abastecimient Explotación mixta 0 0 0 0 0
Santo o doméstico, (agrícola y
Lote No. 4
Riego y pecuaria).
silvicultura

114
49 27348 Rincón Abastecimient Actividad mixta 0.01 0 0.01 0 0.02
Santo o doméstico, (agricola y
Los Pinos
Riego y pecuaria)
silvicultura
50 27396 Rincón Abastecimient Actividad mixta 0.021 0.028 0.948 0 0.997
Santo o doméstico, (agricola y
El Cerezo Riego y pecuaria)
silvicultura,
Abrevaderos
51 27451 Santa Abastecimient Actividad mixta 0 0.03 0.6 0 0.63
Rita El Tiber , La o doméstico, (agricola y
Violeta, Lote a Riego y pecuaria)
Santa Rita silvicultura,
Abrevaderos
52 27792 Guamal Abastecimient Cultivos agrícolas 0.01 0 0 0 0.01
o doméstico, transitorios.
El Granadillo Riego y
silvicultura,
Abrevaderos
53 27793 Guamal Abastecimient Hogares privados 0.01 0 0 0 0.01
El Lote o doméstico con servicios
domesticos
54 27891 Canica Riego y Explotación mixta 0 0 0 0 0
Baja Quincha silvicultura, (agrícola y
Abrevaderos pecuaria).
55 27892 Canica Riego y Cultivos agrícolas 0 0.01 0.03 0 0.04
Los Cerezos
Baja silvicultura, permanentes.
San Juanito
Abrevaderos
56 28439 Guamal Abastecimient Explotación mixta 0.01 0.01 0 0 0.02
Los Cerezos o doméstico, (agrícola y
Abrevaderos pecuaria).
57 28753 El Tobal Abastecimient Cría de ganado 0.01 0.01 0 0 0.02
El Paraiso o doméstico, bovino y bufalino.
Abrevaderos
58 28754 Canica Abrevaderos Explotación mixta 0 0.01 0 0 0.01
Baja El Vergel (agrícola y
pecuaria).

115
59 29071 La San Cayetano Riego y Explotación mixta 0 0.03 0.26 0 0.29
Yeguera lote 4 y santa silvicultura, (agrícola y
ines Abrevaderos pecuaria).
60 29212 Galdame Abastecimient Produccion 0.01 0.01 0 0 0.02
z o doméstico, especializada de
San Francisco
Riego y hortalizas y
silvicultura legumbres
61 29628 Galdame Abastecimient Explotación mixta 0.01 0 0 0 0.01
z Altamira o doméstico (agrícola y
pecuaria).
62 29628 Galdame Abastecimient Explotación mixta 0.01 0 0 0 0.01
z Altamira o doméstico (agrícola y
pecuaria).
63 29633 Guamal Abastecimient Captacion, 0.12 0 0 0 0.12
Bella Vista, el
o doméstico depuracion y
ensueño, el
distribucion de
porvenir
agua
64 30122 Canica Abastecimient Cria especializada 0.01 0.01 0.01 0 0.03
El Eden
Baja o doméstico de ganado vacuno
65 30134 Casa Abastecimient Explotación mixta 0.01 0.01 0 0 0.02
Blanca Casablanca y o doméstico, (agrícola y
el retiro Riego y pecuaria).
silvicultura
66 30314 Rincón Abastecimient Explotación mixta 0.01 0 0 0 0.01
Santo El Vergel o doméstico (agrícola y
pecuaria).
67 30343 Rincón Abastecimient Cría de ganado 0.01 0.01 0 0 0.02
Santo Manzanares o doméstico, bovino y bufalino.
Abrevaderos
68 30660 Guamal Riego y Explotación mixta 0.016 0.01 0.001 0 0.027
El Recuerdo silvicultura (agrícola y
pecuaria).
69 30660 Guamal Riego y Explotación mixta 0.016 0.01 0.001 0 0.027
El Recuerdo silvicultura (agrícola y
pecuaria).

116
70 30662 Buenavis Riego y Actividad mixta 0 0 0 0 0.54
ta Barley II silvicultura (agricola y
pecuaria)
71 30681 Cascajal Abastecimient Captacion, 2.36 0 0 0 2.36
o doméstico, depuracion y
Los Curies
Abrevaderos distribucion de
agua
72 31548 Canica Acuicultura y Explotación mixta 0 0 0.14 0 0.14
Pensilvania pesca (agrícola y
pecuaria).
73 31548 Canica Acuicultura y Explotación mixta 0 0 0.14 0 0.14
Pensilvania pesca (agrícola y
pecuaria).
74 31901 El Valle Abastecimient Captacion, 0.94 0.03 0 0 0.97
o doméstico depuracion y
La Ilusion
distribucion de
agua
75 32169 Altania Riego y Actividad mixta 0 0.02 0.58 0 0.6
El Escondite 1 silvicultura, (agricola y
Abrevaderos pecuaria)
76 32298 Guamal Abastecimient Actividad mixta 0.01 0.03 0.15 0 0.19
o doméstico, (agricola y
Buena Vista Riego y pecuaria)
silvicultura,
Abrevaderos
77 32299 Santa Abastecimient Captacion, 2.04 0 0 0 2.04
Rosa o doméstico depuracion y
Altamira
distribucion de
agua
78 32360 La Abastecimient Captación, 0.71 0 0 0 0.71
Pradera o doméstico tratamiento y
Los Cerezos
distribución de
agua.
79 32651 Canica Riego y Cría de otros 0 0.01 0 0 0.01
La Primavera
Baja silvicultura animales n.c.p.
80 32790 Rincón Abastecimient Captacion, 0.54 0.1 0 0 0.64
San Martin
Santo o doméstico depuracion y

117
 distribucion de
agua
81 33084 Galdame Abastecimient Otras actividades 0.01 0.01 0.01 0 0.03
z o doméstico, de servicios NCP
Finca el Rupe Riego y
silvicultura,
Abrevaderos
82 33434 Páramo Abastecimient Captación, 1.42 0 0 0 1.42
o doméstico tratamiento y
Uval
distribución de
agua.
83 33436 Canica Abastecimient Actividad mixta 0.01 0.02 0.76 0 0.79
Baja o doméstico, (agricola y
Hato Rey Riego y pecuaria)
silvicultura,
Abrevaderos
84 33438 Santuario Riego y Explotación mixta 0 0.02 0.53 0 0.55
Finca San
silvicultura, (agrícola y
Francisco
Abrevaderos pecuaria).
85 33440 Canica Riego y Actividad mixta 0 0.02 0.58 0 0.6
La Esperanza silvicultura, (agricola y
Abrevaderos pecuaria)
86 33459 Canica Riego y Produccion agricola 0 0 0.03 0 0.03
Baja Villa Teresa silvicultura en unidades no
especializadas
87 33867 Pantano Abastecimient Actividad pecuaria 0 0.01 0 0 0.01
de Arce El Diamante o doméstico, no especializada
Abrevaderos
88 33958 La Abastecimient Cria especializada 0.02 0.04 1.46 0 1.52
Yeguera o doméstico, de ganado vacuno
El Tribol, las
Riego y
hortensias
silvicultura,
Abrevaderos
89 33986 Rincón Abastecimient Cría de ganado 0.01 0.03 0 0 0.04
Santo Caballerizas o doméstico, bovino y bufalino.
Abrevaderos

118
90 33989 Rincón Riego y Explotación mixta 0 0 0 0 0
Santo Buenavista silvicultura, (agrícola y
Abrevaderos pecuaria).
91 34098 Buenavis Riego y Actividad mixta 0 0.01 0.09 0 0.1
ta Finca los Alisos silvicultura, (agricola y
Abrevaderos pecuaria)
92 34102 Galdamez Riego y Actividad mixta 0 0.01 0.13 0 0.14
Lote 01 las
silvicultura, (agricola y
acacias
Abrevaderos pecuaria)
93 34111 Páramo Abastecimient Ganadería. 0.01 0.02 0 0 0.03
El Portachuelo
o doméstico
94 34165 Galdamez Abastecimient Actividades de los 0 0 0 0 0
o doméstico hogares
individuales como
El Refugio
empleadores de
personal
doméstico.
95 34166 Tibagota Riego y Explotación mixta 0 0.02 0.76 0 0.78
Tannen Bush silvicultura, (agrícola y
Abrevaderos pecuaria).
96 34243 Guamal Abastecimient Explotación mixta 0.01 0.01 0.02 0 0.04
o doméstico, (agrícola y
Finca el
Riego y pecuaria).
Paraiso
silvicultura,
Abrevaderos
97 34576 Galdamez Abastecimient Explotación mixta 0.01 0.01 0.01 0 0.03
o doméstico, (agrícola y
La Loma Riego y pecuaria).
silvicultura,
Abrevaderos
98 34585 Cascajal Abastecimient Actividad pecuaria 0.01 0.02 0 0 0.03
Curies o doméstico, no especializada
Abrevaderos
99 34602 Guamal Abastecimient Explotación mixta 0.01 0 0 0 0.01
Lote Jerico o doméstico, (agrícola y
Abrevaderos pecuaria).

119
100 34906 Canica Riego y Explotación mixta 0 0.03 0.12 0 0.15
San Rufino silvicultura, (agrícola y
Abrevaderos pecuaria).
101 34913 Tibagota Riego y Produccion agricola 0 0 0.1 0 0.1
Tibagota Alta,
silvicultura en unidades no
Tibagota Baja
especializadas
102 35140 La Riego y Actividad mixta 0 0.08 2.14 0 2.22
Nacedero San
Pradera silvicultura, (agricola y
Cayetano
Abrevaderos pecuaria)
103 35216 Santuario Riego y Actividad pecuaria 0 0.02 0.66 0 0.68
tenerife silvicultura, no especializada
Abrevaderos
104 35641 Canica Abastecimient Produccion agricola 0.02 0 0 0 0.02
finca la
o doméstico en unidades no
primavera
especializadas
105 35642 El Abastecimient Actividad mixta 0.01 0.01 0.06 0 0.08
Pantano o doméstico, (agricola y
finca el angel Riego y pecuaria)
silvicultura,
Abrevaderos
106 35984 Canica Riego y Produccion agricola 0 0 0.07 0 0.07
Baja el cerezo silvicultura en unidades no
especializadas
107 35987 Guamal Abastecimient Explotación mixta 0 0 0.01 0 0.01
sinay o doméstico (agrícola y
pecuaria).
108 35990 Guamal Abastecimient Ganadería. 0.01 0 0 0 0.01
las villas o doméstico,
Abrevaderos
109 36070 Guamal Abastecimient Explotación mixta 0 0 0 0 0
el mirador o doméstico, (agrícola y
Abrevaderos pecuaria).
110 36075 La Industrial Captación, 0 0 0 1 1
Pradera tratamiento y
finca la gloria
distribución de
agua.

120
111 36078 Galdame Abastecimient Actividad mixta 0.02 0 0 0 0.02
z bodegas o doméstico (agricola y
pecuaria)
112 36083 Altania Abastecimient Actividad mixta 0.01 0 0 0 0.01
el molino o doméstico (agricola y
pecuaria)
113 36158 Canica Abastecimient Cultivo de 0 0 0 0 0
Baja o doméstico, hortalizas, raíces y
tres esquinas Riego y tubérculos.
silvicultura,
Abrevaderos
114 36388 Canica Abastecimient Explotación mixta 0 0 0 0 0
el picacho
Baja o doméstico, (agrícola y
central
Abrevaderos pecuaria).
115 36541 El Tobal Abastecimient Explotación mixta 0.01 0.01 0.71 0 0.73
o doméstico, (agrícola y
aguas calientes Riego y pecuaria).
silvicultura,
Abrevaderos
116 36542 Altania Abastecimient Actividad mixta 0 0 0 0 0
o doméstico, (agricola y
la lomita de la
Riego y pecuaria)
gaviota
silvicultura,
Abrevaderos
117 36575 La Union Abastecimient Actividad mixta 0.01 0.01 0.03 0 0.05
mirador alto o doméstico (agricola y
pecuaria)
118 36923 Guamal Riego y Explotación mixta 0 0.01 0.1 0 0.11
san jose silvicultura, (agrícola y
Abrevaderos pecuaria).
119 36925 Pueblo Abastecimient Explotación mixta 0.01 0.01 0.2 0 0.22
Viejo o doméstico, (agrícola y
la victoria Riego y pecuaria).
silvicultura,
Abrevaderos

121
120 36926 Guamal Riego y Explotación mixta 0 0 0.42 0 0.43
la esperanza y
silvicultura, (agrícola y
lote 16
Abrevaderos pecuaria).
121 37046 Galdamez Abastecimient Hogares privados 0.01 0 0 0 0.01
el refugio o doméstico, con servicios
Abrevaderos domesticos
122 37047 Galdamez Abastecimient Ganadería. 0.01 0 0 0 0.01
o doméstico,
san isidro Riego y
silvicultura,
Abrevaderos
123 37048 Galdamez la gloria , el Abastecimient Explotación mixta 0.01 0 0 0 0.01
retamal, o doméstico, (agrícola y
altamira Abrevaderos pecuaria).
124 37052 Galdamez Abastecimient Actividad mixta 0.02 0 0 0 0.02
san miguel o doméstico (agricola y
pecuaria)
125 37053 Galdamez rancho Abastecimient Actividad mixta 0.01 0 0 0 0.01
antamaria, o doméstico, (agricola y
rancho santa Abrevaderos pecuaria)
maria
126 37227 Centro Riego y Actividad mixta 0 0.01 0.07 0 0.08
lote zaragoza
silvicultura, (agricola y
caracoli
Abrevaderos pecuaria)
127 37393 Páramo Abastecimient Actividades no 0.01 0 0 0 0.01
o doméstico diferenciadas de
los hogares
ventorillo individuales como
productores de
bienes para uso
propio.
128 37396 Galdame Riego y Explotación mixta 0 0.01 0.04 0 0.05
z covadonga silvicultura, (agrícola y
Abrevaderos pecuaria).
129 37398 La Abastecimient Explotación mixta 0.01 0.01 0.09 0 0.11
el mortiño, el
Pradera o doméstico, (agrícola y
mortiño
Riego y pecuaria).

122
 silvicultura,
Abrevaderos
130 37598 Guamal Abastecimient Otros tipos de 0.01 0 0 0 0.01
finca villaluz
o doméstico alojamiento n.c.p.
131 37599 Tibagota Abastecimient Actividad mixta 0.01 0.01 0.22 0 0.24
o doméstico, (agricola y
chimita dos,
Riego y pecuaria)
chimita dos
silvicultura,
Abrevaderos
132 37898 Galdame Abastecimient Actividades de 0.01 0 0 0 0.01
z san jose o doméstico apoyo a la
ganadería.
133 38300 El Riego y Ganadería. 0 0.01 0 0 0.01
Pantano el carrizal silvicultura,
Abrevaderos
134 38369 Guamal Abastecimient Actividades de 0 0 0 0 0
o doméstico servicios, agricolas
y ganaderos
san carlos
exceptos las
actividades
veterinarias
135 38455 Rincón miraflores , san Riego y Actividades de 0 0 0.1 0 0.1
Santo jorge , san luis silvicultura apoyo a la
, san manuel ganadería.
136 38520 Centro Riego y Cría de ganado 0 0.01 0.16 0 0.17
loma laguna ii silvicultura, bovino y bufalino.
Abrevaderos
137 38566 Galdame Abastecimient Actividades de 0.01 0.01 0.1 0 0.12
z o doméstico, servicios, agricolas
Riego y y ganaderos
el eucalipto
silvicultura, exceptos las
Abrevaderos actividades
veterinarias
138 38805 El Tobal Abastecimient Explotación mixta 0.01 0.01 0 0 0.02
la regadera 1 o doméstico, (agrícola y
Abrevaderos pecuaria).

123
139 38887 Cascajal Abastecimient Explotación mixta 0.01 0.02 0.07 0 0.1
o doméstico, (agrícola y
el mulato Riego y pecuaria).
silvicultura,
Abrevaderos
140 38888 Galdame alarcon - la Abrevaderos Cria especializada 0 0 0 0 0
z esperanza de ganado vacuno
141 39123 Canica Abastecimient Ganadería. 0 0 0 0 0
o doméstico,
finca los
Riego y
copetones
silvicultura,
Abrevaderos
142 39187 Tibagota Abastecimient Actividad mixta 0.01 0.01 0.21 0 0.23
o doméstico, (agricola y
chimita uno Riego y pecuaria)
silvicultura,
Abrevaderos
143 39510 El Tobal Abrevaderos Explotación mixta 0.01 0.01 0.08 0 0.1
el conflicto (agrícola y
pecuaria).
144 39726 Canica Abastecimient 0.01 0.01 0.01 0 0.03
agua fresca o doméstico,
Abrevaderos
145 39729 Guamal finca la Abastecimient Explotación mixta 0.01 0 0.02 0 0.03
primavera (san o doméstico (agrícola y
gabriel) pecuaria).
146 39989 Guamal Abrevaderos Actividad mixta 0.01 0 0 0 0.01
jerico (agricola y
pecuaria)
147 39990 Casa Abastecimient Hogares privados 0.01 0 0 0 0.01
Blanca san roque o doméstico con servicios
domesticos
148 39991 El Tobal Abrevaderos Actividad mixta 0 0.01 0 0 0.01
el volcan (agricola y
pecuaria)

124
149 40107 Guamal Riego y Ganadería. 0.01 0.01 0 0 0.02
las nieves silvicultura,
Abrevaderos
150 40214 Páramo Abastecimient Explotación mixta 0.01 0.01 0 0 0.02
lote 2 o doméstico, (agrícola y
Abrevaderos pecuaria).
151 40480 Pantano Riego y Actividad mixta 0 0.01 0 0 0.01
lote 12 la
de Arce silvicultura, (agricola y
reserva
Abrevaderos pecuaria)
152 40535 Pantano Abrevaderos Cria especializada 0 0.01 0 0 0.01
de Arce el mortiño de ganado vacuno
153 40721 Guamal Abastecimient Cría de ganado 0.01 0.01 0 0 0.02
el recuerdo o doméstico, bovino y bufalino.
Abrevaderos
154 40939 Pantano Abrevaderos Explotación mixta 0 0.01 0 0 0.01
de Arce el tunal (agrícola y
pecuaria).
155 41095 Altania Abastecimient Captación, 11.16 0 0 0 11.16
o doméstico tratamiento y
lote b
distribución de
agua.
156 41335 Guamal Abastecimient Cría de ganado 0.01 0.01 0 0 0.02
tres esquinas
o doméstico, bovino y bufalino.
matepaja
Abrevaderos
157 41336 La Union san agustin Abrevaderos Ganadería. 0 0.01 0 0 0.01
158 41654 Centro Abastecimient Cria especializada 0 0 0 0 0
santa ines
o doméstico de ganado vacuno
159 41682 Canica Abastecimient Captacion, 0.63 0.07 0 0 0.7
o doméstico, depuracion y
finca coscojita
Abrevaderos distribucion de
agua
160 41725 Centro Riego y Ganadería. 0 0.03 0.21 0 0.24
lote no. 5 san
silvicultura,
luis
Abrevaderos
161 41782 Galdame Abastecimient Hogares privados 0 0 0 0 0
z lote 2 el pencil o doméstico con servicios
domesticos

125
162 41783 Galdame Abastecimient Explotación mixta 0 0 0 0 0
z el pino o doméstico (agrícola y
pecuaria).
163 41784 Pantano Abrevaderos Cria especializada 0 0.01 0 0 0.01
de Arce providencia ii de ganado vacuno
164 42195 Galdame Abastecimient Actividades de los 0.01 0 0 0 0.01
z o doméstico hogares
individuales como
el jazmin
empleadores de
personal
doméstico.
165 42240 Guamal Abastecimient Cría de ganado 0.01 0.01 0.07 0 0.09
o doméstico, bovino y bufalino.
el recreo Riego y
silvicultura,
Abrevaderos
166 42242 Galdame Abastecimient Hogares privados 0.02 0 0 0 0.02
z villa luz o doméstico con servicios
domesticos
167 42245 Galdame Abrevaderos Actividad pecuaria 0 0 0 0 0
el peñon
z no especializada
168 42246 Galdame Abastecimient Explotación mixta 0.02 0.01 0 0 0.03
el cucharo lote
z o doméstico, (agrícola y
24
Abrevaderos pecuaria).
169 42266 Canica Riego y Actividad mixta 0 0.12 0.01 0 0.13
finca la
silvicultura, (agricola y
panchita lote 6
Abrevaderos pecuaria)
170 42271 Casa Abrevaderos Ganadería. 0 1 0 0 1
el triunfo
Blanca
171 42318 Galdame Abastecimient Actividades de 0.01 0.01 0 0 0.02
z o doméstico, servicios, agricolas
Abrevaderos y ganaderos
el recuerdo
exceptos las
actividades
veterinarias

126
172 42319 Galdamez Abastecimient Actividad pecuaria 0.01 0.01 0 0 0.02
el susca lote 26 o doméstico, no especializada
Abrevaderos
173 42322 Galdamez Abastecimient Actividades de los 0.01 0 0 0 0.01
o doméstico hogares
individuales como
el laurel
empleadores de
personal
doméstico.
174 42550 Guamal Abastecimient Produccion agricola 0.01 0 0.01 0 0.02
o doméstico, en unidades no
agua bonita
Riego y especializadas
silvicultura
175 42551 Galdamez Abastecimient Cría de ganado 0.01 0 0 0 0.01
o doméstico, bovino y bufalino.
la corona
Riego y
silvicultura
176 42552 Galdamez Abastecimient Actividad pecuaria 0.02 0.01 0 0 0.23
o doméstico, no especializada
la farfalla Riego y
silvicultura,
Abrevaderos
177 42840 Galdamez Abastecimient Cría de ganado 0.01 0.01 0 0 0.02
la dalia, la dalia o doméstico, bovino y bufalino.
Abrevaderos
178 43281 Galdamez Abastecimient Explotación mixta 0.01 0.01 0 0 0.02
el ciruelo o doméstico, (agrícola y
Abrevaderos pecuaria).
179 43656 Los Riego y Explotación mixta 0 0.01 0.01 0 0.02
Llanitos la providencia silvicultura, (agrícola y
Abrevaderos pecuaria).
180 43870 La Abrevaderos Cría de ganado 0 0.01 0 0 0.01
santana
Pradera bovino y bufalino.
181 43874 La Abrevaderos Cria especializada 0 0.01 0 0 0.01
san rafael
Pradera de ganado vacuno

127
182 43887 Centro Abrevaderos Explotación mixta 0.02 0.01 0 0 0.03
vista bella (agrícola y
pecuaria).
183 44464 La Union Riego y Actividad mixta 0.01 0.01 0 0 0.02
el porvenir , la
silvicultura, (agricola y
casona
Abrevaderos pecuaria)
184 44469 La Union Abastecimient Actividades no 0.01 0 0 0 0.01
o doméstico diferenciadas de
los hogares
villa gabriela individuales como
productores de
bienes para uso
propio.
185 44470 La Riego y Explotación mixta 0.01 0.01 0.01 0 0.03
Pradera los laureles silvicultura, (agrícola y
Abrevaderos pecuaria).
186 44598 La Union Riego y Cultivos agrícolas 0.01 0.01 0 0 0.02
casa blanca,
silvicultura, transitorios.
currucuy
Abrevaderos
187 44601 La Union Abrevaderos Cría de ganado 0.01 0.01 0 0 0.02
el recuerdo
bovino y bufalino.
188 44602 La Union Abastecimient Cría de ganado 0 0.01 0 0 0.01
las margaritas
o doméstico bovino y bufalino.
189 44603 La Union Abrevaderos Cría de ganado 0 0.01 0 0 0.01
migajitas
bovino y bufalino.
190 44664 Guamal Abrevaderos Actividades no 0 0.01 0 0 0.01
diferenciadas de
los hogares
lote san jorge individuales como
productores de
servicios para uso
propio.
191 44665 Guamal Abrevaderos Actividades no 0 0.01 0 0 0.01
diferenciadas de
la helida los hogares
individuales como
productores de

128
 bienes para uso
propio.

192 44666 Guamal Riego y Cría de ganado 0 0.01 0.01 0 0.02

el cerezo silvicultura, bovino y bufalino.
Abrevaderos
193 44667 Guamal Riego y Explotación mixta 0 0.01 0.01 0 0.02
la reserva silvicultura, (agrícola y
Abrevaderos pecuaria).
194 44853 Tibagota los sauces Abastecimient Actividad mixta 0 0.01 0 0 0.01
(antes finca el o doméstico, (agricola y
diamante) Abrevaderos pecuaria)
195 45114 Guamal Abastecimient Actividades de 0.01 0 0 0 0.01
o doméstico servicios, agricolas
y ganaderos
lote 16
exceptos las
actividades
veterinarias
196 45636 Pantano Abastecimient Explotación mixta 0.01 0.01 0 0 0.02
de Arce santa ines o doméstico (agrícola y
pecuaria).
197 45675 Canica Riego y Explotación mixta 0 0.01 0.17 0 0.18
Baja el medio silvicultura, (agrícola y
Abrevaderos pecuaria).
198 45792 La Union Abastecimient Explotación mixta 0.01 0.01 0 0 0.02
el triunfo o doméstico (agrícola y
pecuaria).
199 45802 Galdame Abastecimient Actividades de los 0 0 0 0 0
z o doméstico hogares
individuales como
los pinos
empleadores de
personal
doméstico.
200 45828 Guamal Abastecimient Explotación mixta 0.01 0.01 0.11 0 0.13
piedra colorada o doméstico (agrícola y
pecuaria).

129
201 45850 La Union Abastecimient Explotación mixta 0.01 0.01 0 0 0.02
bodegas o doméstico (agrícola y
pecuaria).
202 46285 Guamal Abrevaderos Cría de ganado 0 0.01 0 0 0.01
altamira
bovino y bufalino.
203 46337 Guamal predio Abrevaderos Cría de ganado 0 0.01 0 0 0.01
palestina bovino y bufalino.
204 46338 Puente Riego y Cría de ganado 0 0 0.01 0 0.01
de Piedra la pocholita silvicultura bovino y bufalino.
205 46950 Rincón Abastecimient Explotación mixta 0 0 0 0 0
Santo o doméstico, (agrícola y
santillana del
Riego y pecuaria).
lago
silvicultura,
Abrevaderos
206 47504 Galdame Riego y Cultivo de 0 0.03 0 0 0.03
las acacias lote
z silvicultura hortalizas, raíces y
1
tubérculos.
207 47972 Rincón Abastecimient Actividad mixta 0.01 0.01 0 0 0.02
Santo el triunfo o doméstico (agricola y
pecuaria)
208 48034 Guamal Riego y Explotación mixta 0 0.01 0.433 0 0.443
linaza ii silvicultura (agrícola y
pecuaria).
209 48349 Altania Riego y Cría de ganado 0 0 0.15 0 0.15
loma laguna i silvicultura, bovino y bufalino.
Abrevaderos
210 48357 Guamal Abastecimient Hogares privados 0.01 0 0 0 0.01
el triangulo o doméstico con servicios
domesticos
211 49285 Galdame Abrevaderos Cría de ganado 0 0.01 0 0 0.01
el nogal
z porcino.
212 49763 Centro Abastecimient Cría de ganado 0.01 0.01 0 0 0.02
o doméstico, bovino y bufalino.
san isidro Riego y
silvicultura,
Abrevaderos

130
213 49813 Guamal Abastecimient Explotación mixta 0.05 0.02 0.02 0 0.09
o doméstico, (agrícola y
san kitts Riego y pecuaria).
silvicultura,
Abrevaderos
214 49981 Guamal Abastecimient Explotación mixta 0.0002 0.0001 0.0007 0 0.001
o doméstico, (agrícola y
lote 2 tablanca
Riego y pecuaria).
1
silvicultura,
Abrevaderos
215 50099 El Valle Riego y Cría de ganado 0 0.01 0.15 0 0.16
san fruto silvicultura, bovino y bufalino.
Abrevaderos
216 50151 Guamal el ático , el Abastecimient Explotación mixta 0.0001 0 0 0 0.0001
recuerdo, los o doméstico, (agrícola y
encenillos, Riego y pecuaria).
santa juana silvicultura
217 50193 Guamal Abastecimient Explotación mixta 0.01 0 0 0 0.01
o doméstico, (agrícola y
vista hermosa Riego y pecuaria).
silvicultura,
Abrevaderos
218 51110 Canica Riego y Cultivo de 0 0 0.16 0 0.16
Baja el oasis silvicultura hortalizas, raíces y
tubérculos.
219 51473 Galdame Abastecimient 0.0002 0 0 0 0.0002
el triangulo
z o doméstico
220 51570 Galdame Abastecimient Actividades no 0.0001 0 0 0 0.0001
z o doméstico diferenciadas de
los hogares
veraguita individuales como
productores de
bienes para uso
propio.
221 51571 Galdame Abastecimient Actividades no 0.01 0 0 0 0.01
z el pino o doméstico diferenciadas de
los hogares

131
 individuales como
productores de
servicios para uso
propio.
222 52010 La Riego y Actividades no 0 63 2 0 65
Pradera silvicultura, diferenciadas de
Abrevaderos los hogares
villa nueva de
individuales como
la pradera
productores de
bienes para uso
propio.
223 52811 Canica Abastecimient Actividades no 1 1 3 0 5
Baja o doméstico, diferenciadas de
Riego y los hogares
lote 1 los tunos silvicultura, individuales como
Abrevaderos productores de
bienes para uso
propio.
224 53140 Rincón Abastecimient Produccion agricola 0 0 0.02 0 0.02
lote no. 2
Santo o doméstico, en unidades no
buena vista
Abrevaderos especializadas
225 53568 Rincón Abastecimient Actividades de los 0 0.02 0.51 0 0.53
Santo o doméstico, hogares
Riego y individuales como
los cerezos
silvicultura, empleadores de
Abrevaderos personal
doméstico.
226 54889 Casa Abastecimient Actividades no 0.01 0.11 0 0 0.12
Blanca o doméstico, diferenciadas de
el recuerdo
Abrevaderos los hogares
(actualmente
individuales como
denominado
productores de
limay).
bienes para uso
propio.
227 55201 Canica Riego y Explotación mixta 0 0.01 0.17 0 0.18
Baja palmira silvicultura, (agrícola y
Abrevaderos pecuaria).

132
228 55210 Canica Riego y Actividades de 0 0 0.06 0 0.06
Baja silvicultura servicios, agricolas
lote de terreno y ganaderos
el capricho. exceptos las
actividades
veterinarias
229 55367 Rincón Abastecimient Explotación mixta 0.01 0.01 0.05 0 0.07
Santo o doméstico, (agrícola y
la campiña Riego y pecuaria).
silvicultura,
Abrevaderos
230 55414 Tibagota Riego y Explotación mixta 0 0.09 0.24 0 0.33
la cabanga silvicultura, (agrícola y
Abrevaderos pecuaria).
231 55627 Galdame Abastecimient Explotación mixta 1 1 16 0 18
z o doméstico, (agrícola y
la cima Riego y pecuaria).
silvicultura,
Abrevaderos
232 55750 Canica Abrevaderos Cultivo de frutas 0 0 0.14 0 0.14
Baja finca la ilusion tropicales y
subtropicales.
233 55752 Canica Riego y Cultivo de frutas 0 0 0.04 0 0.04
finca santa
Baja silvicultura tropicales y
lucia
subtropicales.
234 55753 Canica Abastecimient Cría de ganado 0 0.01 0 0 0.01
Baja los alisos o doméstico, bovino y bufalino.
Abrevaderos
235 55754 Canica Riego y Explotación mixta 0 0.02 0.16 0 0.18
Baja la piedra silvicultura, (agrícola y
Abrevaderos pecuaria).
236 55755 Canica Riego y Explotación mixta 0 0.01 0.03 0 0.04
Baja la camelia silvicultura, (agrícola y
Abrevaderos pecuaria).
237 55756 Canica Riego y Explotación mixta 0 0.01 0.02 0 0.03
Baja el mirador silvicultura, (agrícola y
Abrevaderos pecuaria).

133
238 55757 Canica Riego y Cultivo de frutas 0 0 0.04 0 0.04
Baja los laureles silvicultura tropicales y
subtropicales.
239 55758 Canica Riego y Explotación mixta 0 0.01 0.15 0 0.16
Baja san antonio silvicultura, (agrícola y
Abrevaderos pecuaria).
240 55759 Canica Riego y Explotación mixta 0 0.01 0.05 0 0.06
Baja macanas silvicultura, (agrícola y
Abrevaderos pecuaria).
241 55760 Canica Abrevaderos Cría de ganado 0 0.06 0 0 0.06
la cabaña
Baja bovino y bufalino.
242 55761 Canica Riego y Explotación mixta 0 0.02 0.16 0 0.18
Baja san jorge silvicultura, (agrícola y
Abrevaderos pecuaria).
243 56263 Altania Abastecimient 0.0001 0 0 0 0.0001
el triunfo
o doméstico
244 56419 Altania Abastecimient Actividades no 0.012 0 0 0 0.012
o doméstico diferenciadas de
los hogares
la hoya individuales como
productores de
bienes para uso
propio.
245 56947 Canica Abastecimient Ganadería. 0.0001 0.0001 0.0002 0 0.0004
Baja lote 3 la o doméstico,
danielita Riego y
silvicultura
246 57242 Páramo Abastecimient 0.0011 0.001 0.0019 0 0.004
el refugio o doméstico,
Abrevaderos
247 57284 Canica Abastecimient Explotación mixta 0 2 0 0 2
las mercedes o doméstico, (agrícola y
Abrevaderos pecuaria).
248 57384 Canica Abastecimient Cría de ganado 0.02 0.04 0 0 0.06
Baja el cigarral o doméstico, bovino y bufalino.
Abrevaderos

134
249 58084 Galdame Abastecimient Actividades no 0.0001 0.0001 0 0 0.0002
z o doméstico, diferenciadas de
Abrevaderos los hogares
el rosal individuales como
productores de
servicios para uso
propio.
250 58304 Altania Abastecimient Ganadería. 0.01 0.02 0 0 0.03
globo el cajon o doméstico,
Abrevaderos
251 58492 Canica Abastecimient Explotación mixta 0.0001 0.0001 0.0013 0 0.0015
Baja o doméstico, (agrícola y
el tunal Riego y pecuaria).
silvicultura,
Abrevaderos
252 60873 La Abastecimient Actividades de 0.02 0.04 0.72 0 0.78
Pradera o doméstico, apoyo a la
raices Riego y agricultura.
silvicultura,
Abrevaderos
253 61785 Páramo Abastecimient Explotación mixta
manitas la
o doméstico, (agrícola y
reserva
Abrevaderos pecuaria).
254 61802 Páramo Abrevaderos Explotación mixta 0 0.01 0 0 0.01
la reforma (agrícola y
pecuaria).
255 62436 Centro Abastecimient Explotación mixta 0.02 0.02 0.88 0 0.92
o doméstico, (agrícola y
el rodeo Riego y pecuaria).
silvicultura,
Abrevaderos
256 62561 Canica Abrevaderos Explotación mixta 0 0.01 0 0 0.01
Baja el estribo (agrícola y
pecuaria).
257 62659 Galdame lote 1 santa Abrevaderos Cría de ganado 0 0.01 0 0 0.01
z maria bovino y bufalino.

135
258 62660 Galdame Abastecimient Actividades no 0.01 0 0 0 0.01
z o doméstico diferenciadas de
los hogares
agua caliente individuales como
productores de
bienes para uso
propio.
259 62711 Canica Abastecimient Cría de ganado 0.012 0.01 0 0 0.022
Baja aposento o doméstico, bovino y bufalino.
Abrevaderos
260 63533 Galdame Abastecimient Actividades no 0.01 0 0 0 0.01
z o doméstico diferenciadas de
los hogares
Canica Baja individuales como
productores de
servicios para uso
propio.
261 63534 Galdame Abrevaderos Explotación mixta 0 0.01 0 0 0.01
lote 2 la
z (agrícola y
victoria
pecuaria).
262 63565 La Union Riego y Explotación mixta 0 0.01 0.04 0 0.05
la cabaña silvicultura, (agrícola y
Abrevaderos pecuaria).
263 64857 El Valle Abastecimient Actividades de 0.01 0.01 0.18 0 0.2
o doméstico, apoyo a la
bariloche Riego y agricultura.
silvicultura,
Abrevaderos
264 64923 Rincón Abastecimient Explotación mixta 0.02 0.016 0.467 0 0.503
Santo o doméstico, (agrícola y
potosí Riego y pecuaria).
silvicultura,
Abrevaderos
265 64924 Rincón Abastecimient Explotación mixta 0.015 0.012 0.5 0 0.527
Santo o doméstico, (agrícola y
el cerezo Riego y pecuaria).
silvicultura,
Abrevaderos

136
266 65152 Casa Abastecimient Actividades de 0.0005 0 0.0029 0 0.0034
Blanca o doméstico, apoyo a la
potrero relator
Riego y agricultura.
silvicultura
267 65796 Rincón Abastecimient Actividades no 0.014 0 0 0 0.014
Santo o doméstico diferenciadas de
sin dirección los hogares
lote 2 san individuales como
francisco productores de
bienes para uso
propio.
268 69401 Páramo Abrevaderos Cría de ganado 0 0.016 0 0 0.016
el tirol
bovino y bufalino.

137
ANEXO D. Calculo de la demanda por sectores económicos

Demanda Demanda Demanda Demanda

MUNICIPIO Domestica Agricola Pecuario Industrial
(litro/año) (litro/año) (litro/año) (litro/año)
Bojaca 2343755520 13352610456 2557859731 110344464
El Rosal 2350551528 13488338246 2561707123 110344464
Funza 2192871528 13335130051 2554706131 110344464
Madrid 2352128328 13508521286 2563283923 110344464
Mosquera 79921684.8 11530324771 2277126259 0
Subachoque 2352759048 13594929926 2567572819 110344464

138