See discussions, stats, and author profiles for this publication at: https://www.researchgate.

net/publication/306013913

MODELACIÓN HIDROLÓGICA Y ESCENARIOS DE CAMBIO CLIMÁTICO EN

CUENCAS DE SUMINISTRO DE AGUA DE LAS CIUDADES LA PAZ Y EL ALTO,
BOLIVIA HYDROLOGIC MODELING AND SCENARIOS OF CLIMATE C....

Article · September 2013

CITATIONS READS

3 1,918

5 authors, including:

Marisa Escobar Nilo A. Lima

Stockholm Environment Institute Universidad Mayor de San Andres
28 PUBLICATIONS   333 CITATIONS    2 PUBLICATIONS   14 CITATIONS   

SEE PROFILE SEE PROFILE

David Purkey David N. Yates

Stockholm Environment Institute National Center for Atmospheric Research
63 PUBLICATIONS   1,829 CITATIONS    92 PUBLICATIONS   4,024 CITATIONS   

SEE PROFILE SEE PROFILE

Some of the authors of this publication are also working on these related projects:

restrained metal beam columns under biaxial load View project

http://www.aecom.com/projects/partnering-adaptation-resilience-para-agua/ View project

All content following this page was uploaded by Marisa Escobar on 10 August 2016.

The user has requested enhancement of the downloaded file.

 Aqua-LAC - Vol. 5 - Nº 2 - Set. 2013. pp. 2 - 4
MODELACIÓN HIDROLÓGICA Y ESCENARIOS DE CAMBIO CLIMÁTICO EN CUENCAS DE
SUMINISTRO DE AGUA DE LAS CIUDADES LA PAZ Y EL ALTO, BOLIVIA

HYDROLOGIC MODELING AND SCENARIOS OF CLIMATE CHANGE ON WATER SUPPLY

CATCHMENTS OF THE CITIES LA PAZ AND EL ALTO, BOLIVIA

Marisa Escobar1, Nilo Lima1, David Purkey1, David Yates2 y Laura Forni1

Resumen
En este artículo se presenta un enfoque útil para la modelación hidrológica de cuencas de suministro de agua potable
incluyendo incertidumbres como el cambio climático. La modelación se realizó bajo la plataforma Water Evaluation And
Planning System (WEAP) y comprende aporte de glaciares, precipitación-escorrentía, infraestructura hidráulica y de-
manda urbana.
El estudio caso está referido a los sistemas de agua de las ciudades de La Paz y El Alto que están ubicadas en la región
Occidental de Bolivia. Estas ciudades son abastecidas de agua por pequeñas cuencas glaciares de la Cordillera Real. Por
el crecimiento acelerado de la población especialmente en la ciudad de El Alto existen problemas en el abastecimiento
de agua, y esto podría exacerbase por los impactos negativos del cambio climático por ejemplo el retroceso de los glaci-
ares. Existe una necesidad urgente de conocer los posibles impactos del cambio climático, para esto es necesario tener
un mayor conocimiento del potencial hídrico de las cuencas. Sin embargo, existen bastantes limitaciones respecto a la
disponibilidad y calidad de los datos hidro-climáticas. Asimismo, no se cuenta con datos de escenarios de clima futuro a
la escala de las cuencas que permitan modelar los posibles cambios esperados. En ese sentido en esta investigación se
desarrolló una herramienta útil para la planificación de recursos hídricos bajo la incertidumbre de datos históricos y futuros.
Los resultados para un horizonte de tiempo 2010-2050 muestran que los glaciares podrían desaparecer en los próximos
30 años, y esto podría tener implicaciones en la respuesta hidrológicas de las cuencas con la reducción de la oferta de
agua. La ciudad de El Alto sería la más afectada, ya que se advierten serios problemas en el abastecimiento de agua. Se
debe pensar en estudiar estrategias de nuevo almacenamiento de agua como medidas de adaptación al cambio climático.
Palabras clave: Cambio climático, modelación hidrológica, retroceso glaciar, planificación de recursos hídricos, Cordillera
Real

Abstract
In this paper we present a useful approach for modeling watershed hydrological water supply including uncertainties such
as climate change. The modeling was carried out under the platform Water Evaluation And Planning System (WEAP) and
includes contribution of glaciers, rainfall-runoff, hydraulic infrastructure and urban demand.
The study case framework focuses on water systems in La Paz and EL Alto cities in the Westside region of Bolivia. These
cities have as water sources, small basins with glacier influence because of the “Cordillera Real” (The Andean Mountains).
The fast population growth rate in these cities, especially El Alto City creates problems in drinking water availability, and it
could increase because of negative effects of climate change. For these reasons is necessary to develop knowledge about
the potential availability in these basins, however there is a lack of reliable weather data which make hard an ideal analysis.
Information about future climatic sceneries in micro scale is not available also (at level of basins in our study area). In that
sense this study has developed useful tools for water resource planning policies under the uncertainty of projected and
historical data.
The results for threshold in time of 2010 to 2050 shows that glaciers could disappear in the next 30 years, and this could
have implications in water resource availability, with the consequently reduction of water supply volume, El Alto city would
be the most affected. Policy makers must to define new water sources to deal with climate change effects in long term.
Keywords: Climate change, hydrologic modeling, glacier retreat, water resources planning, Cordillera Real

INTRODUCCIÓN
El documento técnico del Panel Intergubernamental al ciclo hidrológico, alterando la intensidad y la dis-
sobre el Cambio Climático (IPCC según sus siglas tribución temporal y espacial de la precipitación, de
en el idioma inglés) (Bates et al., 2008) indica que el la escorrentía de superficie y de la recarga de agua,
cambio climático podría afectar de manera notable produciendo impactos diversos sobre diferentes

1
Stockholm Environment Institute–US Center, Davis, California, Water Group; [email protected]
2
National Center for Atmospheric Research; [email protected]

Recibido: 09/08/2013
Aceptado: 04/11/2013

Aqua-LAC - Vol. 5 - Nº. 2 - Sep. 2013 23

Marisa Escobar, Nilo Lima, David Purkey, David Yates y Laura Forni
ecosistemas naturales y actividades humanas. En el los glaciares más grandes van a seguir contrayéndo-
quinto informe del IPCC se reporta que la tempera- se (Ramirez et al., 2001; Francou et al., 2003; Bates
tura mundial se ha incrementado en 0.8 °C desde et al., 2008). Los cambios en las contribuciones de
principios del siglo XX y se estima que a mediados agua de los glaciares afectarán al suministro de agua
del presente siglo de continuar las emisiones de CO2 para múltiples usos, incluyendo la agricultura a pe-
como viene aconteciendo en la actualidad habrá un queña y gran escala, los servicios urbanos de agua y
aumento adicional de temperatura de 1.4 a 2.6°C. La la energía hidroeléctrica.
principal conclusión del mencionado informe es que A lo largo de la región Andina existen además otros
el calentamiento del sistema climático es inequívoco agentes de cambio y presiones sobre los sistemas
y desde 1950 son muchos los cambios observados hidrológicos que proveen agua, como la urbaniza-
en todo el sistema climático que no tienen preceden- ción y explosión demográfica por migración del cam-
tes en los últimos decenios a milenios (IPCC, 2013). po a la ciudad, por ejemplo las ciudades de La Paz
El Cuarto Informe de Evaluación del IPCC muestra y El Alto (Figura 1); las mismas están ubicadas en
que durante el siglo XX, los cambios en la precipita- la Provincia Murillo del Departamento de La Paz del
ción en Latinoamérica fueron irregulares, con dismi- Estado Plurinacional de Bolivia (Figura 2).
nuciones en zonas tales como el flanco occidental de
Los Andes y aumentos en otras zonas como el Cono
Sur y en algunas zonas geográficas de la Región An-
dina en Ecuador, Bolivia y Argentina (IPCC, 2007).
Las temperaturas medias por década muestran una
tendencia ascendente en el Cono Sur y en la Región
Andina, y también las temperaturas mínimas, mien-
tras que las temperaturas máximas no muestran una
tendencia uniforme.
La consecuencia más importante del aumento de las
temperaturas en los últimos 30 años ha sido el retro-
ceso crítico de glaciares en Bolivia, Perú, Colombia
y Ecuador. Esto ha aumentado la escorrentía, pero
los glaciares son importantes depósitos de agua, y
se han disminuido e inclusive agotado en algunos lu-
gares. Los estudios indican que en los próximos 15
años, los pequeños glaciares andinos podrían des- Figura 1. Crecimiento de población en las ciudades
aparecer, afectando la disponibilidad de agua y la de La Paz y El Alto (Fuente: Instituto Nacional de
generación de energía hidroeléctrica, mientras que Estadística)

Figura 2. Área de estudio de la modelación hidrológica del sistema de agua de La Paz y El Alto
(Fuente: Elaboración propia)

24 Aqua-LAC - Vol. 5 - Nº. 2 - Sep. 2013

 Modelación hidrológica y escenarios de cambio climático en cuencas de suministro de agua de las ciudades La Paz y El Alto, Bolivia
En la actualidad operan cuatro sistemas de agua po- para contribuir a la planificación del recurso hídrico de
table en La Paz y El Alto: Achachicala, Pampahasi, la región bajo el contexto de incertidumbre de datos
El Alto y Tilata (Figura 2). La fuente de suministro de históricos y futuros? Para responder a esta pregun-
agua de los tres primeros son pequeñas cuencas hi- ta, en esta investigación se realizó una modelación
drográficas ubicadas en la Cordillera Real; el aporte hidrológica bajo la plataforma Water Evaluation And
de estas cuencas se almacena en represas y se con- Planning System (WEAP) en las cuencas de actual
duce a las plantas de tratamiento. El sistema Tilata suministro de las ciudades La Paz y El Alto, y en las
aprovecha aguas subterráneas mediante 32 pozos posibles cuencas alternativas planteadas en el Plan
de una profundidad promedio de 90 metros. Desde el Maestro de Agua Potable del año 1994. Asimismo se
2007, el abastecimiento de agua a ambas ciudades buscó evaluar cuales podrían ser las implicaciones
ha estado a cargo de la Empresa Pública Social de del cambio del clima futuro en el abastecimiento de
Agua y Saneamiento (EPSAS). agua en un horizonte de tiempo 2011-2050 bajo un
Las cuencas de suministro de El Alto cuentan con ambiente de incertidumbres. El enfoque y el tipo de
una importante presencia de glaciares, el derreti- modelo desarrollado podrían ser útiles en el proceso
miento acelerado de glaciares debido a los cambios de toma de decisiones en la planificación de recur-
globales podrían tener implicaciones en la respuesta sos hídricos y adaptación al cambio climático.
hidrológica de las cuencas (Ramírez, 2008). El aná-
lisis multitemporal de los glaciares de la Cordillera METODOLOGÍA
Real sobre restitución de fotografías aéreas muestra
que los glaciares perdieron el 43% de su volumen Modelo WEAP
entre 1963 y 2006 (Soruco et al., 2009). En el caso La plataforma WEAP como herramienta para la pa-
específico de las cuencas Condoriri ha perdido 44% nificación de recursos hídricos incluye un modelo de
de su superficie glaciar y la cuenca Tuni el 55% entre balance hídrico para representar los procesos hidro-
1956 y 2006 y las tendencias indican que desapa- lógicos dentro de un sistema de cuencas (Yates et
recerían en los próximos 30 años (Ramírez 2008). al., 2005a, 2005b). Este algoritmo permite represen-
Un análisis para el periodo 1939-1998 determinaron tar el comportamiento del balance hídrico desde las
que la temperatura en los Andes tropicales ha au- precipitaciones en la cuenca alta con características
mentado en 0.1 a 0.11 °C/década desde 1939 (Vuille fisiográficas de montaña con pendientes escarpadas
& Bradley, 2000). a muy escarpadas y valles en forma de V, y cómo
Existe una necesidad urgente de comprender los im- estas contribuyen al flujo base de las aguas subterrá-
pactos del cambio climático sobre el suministro de neas. Segundo, en la cuenca baja, zona de llanura,
agua, de manera que esto permita tomar acciones contribuyen a los acuíferos aluviales y estos pueden
para garantizar la cobertura de los requerimientos de contribuir a los ríos o pérdidas por filtraciones, es de-
agua futuros. Para evaluar el potencial hidrológico cir existe una interacción.
de las cuencas existen limitaciones referentes a la Se trata de un modelo cuasi físico unidimensional
disponibilidad y calidad de los datos. El monitoreo que puede utilizarse a diferentes pasos de tiempo.
de variables climáticas como temperatura y precipi- En el caso de este estudio, se planteó el uso de un
tación es muy deficiente, ya sea a nivel temporal y paso de tiempo mensual que permite una mirada
espacial. Por iniciativa del programa GREAT-ICE del hacia las proyecciones de cambio climático futuras.
Instituto francés de Investigación para el Desarrollo El algoritmo requiere datos de entrada como precipi-
(IRD) y el Instituto de Hidráulica e Hidrología de la tación, temperatura, humedad relativa, viento, punto
Universidad Mayor de San Andrés (IHH-UMSA) en de derretimiento, punto de congelamiento de nieve,
el año 1991 se inició el monitoreo de algunos glacia- latitud, y cobertura vegetal; además se debe delimi-
res tropicales en la Cordillera Real (Zongo), y pos- tar las unidades básicas de modelación, catchments,
teriormente instalaron sensores para el monitoreo en base a las cuencas y cobertura vegetal (CCG-UC
de niveles de agua en las cuencas Condoriri, Tuni y & SEI, 2009). De esta manera se estiman los compo-
Huayna Potosí (sistema El Alto). Por otro lado, la EP- nentes del balance hidrológico como evapotranspira-
SAS en el año 2000 instaló pluviómetros cerca de los ción, infiltración, escorrentía superficial, escorrentía
embalses. Sin duda estas acciones fueron un gran sub-superficial y flujo base. Los datos de caudales
avance para poder comprender el ciclo hidrológico observados permiten la calibración del modelo. En
de las cuencas, es así que en base a estos esfuer- la zonificación de las unidades básicas de modela-
zos se han realizado una serie de trabajos que per- ción en zonas montañosas se debe incluir la altitud,
mitieron generar datos y conocimiento inicial sobre lo cual permite una mejor representación de los pro-
el comportamiento de la hidrología regional. La otra cesos hidrológicos (CCG-UC & SEI, 2009).
limitación es la falta de datos de escenarios de clima En las cuencas de suministro de La Paz y El Alto
futuro a la escala de cuencas pequeñas, que permi- existen glaciares tropicales, por tanto existe la nece-
tan modelar los aportes de agua futuro y estimar las sidad de evaluar el aporte de los mismos, el WEAP
perdidas probables (Ramírez, 2008). no cuenta con un modelo hidroglaciológico. Sin em-
En este contexto, la pregunta de investigación es bargo se han desarrollado algunos metodologías en
¿qué tipo de herramientas se pueden desarrollar la cuenca del río Santa-Perú modificando el modelo

Aqua-LAC - Vol. 5 - Nº. 2 - Sep. 2013 25

Marisa Escobar, Nilo Lima, David Purkey, David Yates y Laura Forni
grado-día a grado-mes (Condom et al., 2012), otro cas en estudio. Por otro lado, el SENAMHI cuenta
en los glaciares de Himalaya con modificación al con puntos de medición en la zona aguas abajo de
anterior usando la radiación neta en vez del factor las cuencas Hampaturi y Kaluyo, pero debido a su
grado-mes (Mehta et al., 2010) y en la cuenca Pau- ubicación estos puntos son más representativos de
te-Ecuador simplificando la representación del mo- zonas donde existe bastante intervención humana,
delo grado-mes, denominándolo Ice-KISS (Pouget, como agricultura y zonas peri-urbanas.
2011). Estos algoritmos son acoplados al modelo de Los datos de precipitación fueron procesados para
precipitación-escorrentía en WEAP. En esta investi- el análisis de homogeneidad por el método Vector
gación se utilizó modelo desarrollado en la cuenca Regional (Hiez, 1977; Brunet-Moret, 1979) que está
del rio Santa. incorporado en el software Hydraccess (Vauchel,
El modelo glaciar grado mes modela la dinámica de 2005). Las estaciones Hichucota, Tuni, Milluni, Alto
la extensión de la superficie del glaciar y el aporte de Achachicala, Incachaca y Hampaturi fueron selec-
agua a la salida de la cuenca glaciar (Condom et a, cionadas como estaciones base para la modelación
2012). Cada banda de elevación en cada cuenca se en WEAP, ya que estas se encuentran ubicadas al
divide en zonas glaciares y no glaciares, los glacia- interior de las cuencas. Los vacíos de datos en los
res son conceptualizados por el modelo grado-mes. registros fueron rellenados por correlación múltiple
El área inicial del glaciar que puede ser obtenido de con estaciones con registros extensos y continuos
imágenes satelitales o fotografías aéreas es usada como San Calixto, Laykacota y El Alto-Aeropuerto.
para el cálculo del volumen glaciar inicial mediante Con base en la información procesada se genera-
una ecuación empírica (Bahr et al., 1997). Para cada ron isoyetas mediante interpolación con funciones de
mes del año hidrológico el volumen de escurrimiento base radial (Figura 3). Por otra parte, para el cálculo
superficial en la cuenca es la suma de la fusión de de la temperatura a diferentes altitudes se usó un
nieve, y glaciar y el escurrimiento modelado por Soil gradiente térmico vertical de 0.800°C/100 m, esti-
Moisture (modelo hidrológico de WEAP) en zonas mado después de analizar las estaciones Achumani,
no glaciares. Al final del año hidrológico calcula el Laykacota, El Alto-Aeropuerto, Alto Achachicala, Hi-
balance de masas para evaluar los cambios en el vo- chucota y Mevis (Zongo).
lumen total del glaciar (Condom et al., 2012). El mo-
delo glaciar utiliza solo tres parámetros: T0, que es la
temperatura crítica para la fusión glacial, aice y asonw,
factor grado-mes de hielo y nieve respectivamente.

Disponibilidad y procesamiento de datos de

entrada
En las cuencas de suministro de agua potable (ver-
tiente Oeste de la Cordillera Real) existe una defi-
ciencia en el monitoreo de las variables de base
climática. EPSAS cuenta con un base de datos de
pluviometría desde el año 2000, pero esta red está
ubicada en las represas (zona baja). En general, no
existe información en la zona alta, sin embargo exis-
ten datos en la cuenca glaciar Zongo desde el año
1991, la cual queda en la vertiente Este de la Cordi-
llera Real. En las ciudades La Paz y El Alto existen
datos extensos de precipitación, temperatura, hume-
dad relativa y vientos en las estaciones meteoroló-
gicas San Calixto, Laykacota y El Alto-Aeropuerto
desde el año 1917, 1945 y 1942 respectivamente. Figura 3. Isoyetas interpoladas para las cuencas de
De las cuatro posibles nuevas cuencas de suministro estudio (Fuente: Elaboración propia)
ubicadas al Norte del área de modelación, en tres no
existente ninguna información climática. En la cuen- Para delimitar las unidades básicas de modelación,
ca Khara Khota donde si existen datos, se encuentra se utilizó un modelo numérico de elevación con reso-
información sobre precipitación, humedad relativa, lución espacial de 20 m elaborado con la Carta Na-
temperatura y viento. cional escala 1:50.000 del Instituto Geográfico Militar
Respecto a la información sobre caudales medidos, (IGM-Bolivia); el modelo permitió delimitar las cuen-
en la mayoría de las cuencas de suministro actual cas y realizar una clasificación de bandas de eleva-
especialmente en Tuni, Condoriri y Huayna Potosí ción con intervalos iguales de 300 m. La cobertura
existen base de datos del programa GREAT ICE-IRD. y uso de la tierra se obtuvo mediante interpretación
Estos puntos de medición son de mucha importancia visual y digital sobre una imagen Landsat TM 5 de
ya que están ubicados a una altitud por encima de fecha septiembre/1994, donde se definieron cinco
4500 m donde existe poca intervención humana que categorías: glaciares, vegetación dispersa de gramí-
en general es la característica de todas las cuen- neas y arbustos, bofedales, cuerpos de agua y aflo-

26 Aqua-LAC - Vol. 5 - Nº. 2 - Sep. 2013

 Modelación hidrológica y escenarios de cambio climático en cuencas de suministro de agua de las ciudades La Paz y El Alto, Bolivia
ramientos rocosos. La intersección de cada banda modelos climáticos, utilizando un modelo estadísti-
de elevación con los límites de la cuenca y la cober- co Bayesiano para sintetizar la información de varios
tura de la tierra, permitió definir las unidades básicas GCMs en una función de densidad de probabilidad
de modelación. de cambio en la temperatura y precipitación (Tebaldi
et al., 2005).
Escenarios de cambio climático Primero se realizó un análisis estadístico de la es-
El dominio de la modelación hidrológica es muy pe- tructura de datos históricos de la estación meteoro-
queño, ya que el área de las cuencas modeladas lógica con registro más extenso, que en este caso
es menor a 100 km2. El enfoque utilizado para de- es la estación San Calixto con datos desde el año
sarrollar los escenarios de cambio climático fue el 1917 (Figura 4), la cual está ubicada a una altitud de
uso directo de GCMs, en este caso la probabilidad 3.658 msnm. A partir de esto se generó una matriz
de condiciones futuras con Regional Climate-Chan- de transición de dos estados de probabilidad (Tabla
ge Projection from Multi-Model Ensembles (RCPM, 1). Con base en este análisis se realizó la extensión
http://rcpm.ucar.edu/), que son resultados de un de series de tiempo para las estaciones base usada
análisis estadístico de proyecciones de diferentes en WEAP.

Figura 4. Precipitación histórica de la estación San Calixto (Fuente: Elaboración propia).

Tabla 1. Matriz de transición de estados de probabilidad en la estación meteorológica San Calixto

AÑO ACTUAL PROBABILIDAD SIGUIENTE AÑO PROBABILIDAD SIGUIENTE
DEBAJO DE PROMEDIO (1) AÑO ENCIMA DE PROMEDIO (2)
DEBAJO DE PROMEDIO 55% 45%
ENCIMA DE PROMEDIO 43% 56%

Fuente: Elaboración propia

Segundo, se usó el RCPM para la región de interés

(Figura 5). Esta es una técnica estadística para cal-
cular los cambios esperados del clima en el futuro
para una región determinada sobre la base de datos
de observación y los resultados de un conjunto de
Modelos Climáticos Globales (GCM) y suposiciones
sobre emisiones futuras. Finalmente se aplicaron
los cambios esperados del análisis probabilístico del
RCPM a las series sintéticas mensuales de precipita-
ción. Se realizó una combinación de gamas de opcio-
nes, desde escenarios secos y muy calientes (Figura
6), que en adelante se denominan como pesimistas,
hasta escenarios húmedos y calientes (Figura 7), de-
nominados como optimistas de acuerdo con la forma
en que estos escenarios podrían beneficiar al siste-
ma al generar mayor disponibilidad de agua.

Figura 5. Dominio del análisis en RCPM

(Fuente: Elaboración propia)

Aqua-LAC - Vol. 5 - Nº. 2 - Sep. 2013 27

Marisa Escobar, Nilo Lima, David Purkey, David Yates y Laura Forni

Figura 6. Combinación de escenario seco y muy caliente (Fuente: Elaboración propia)

Figura 7. Combinación de escenario húmedo y caliente (Fuente: Elaboración propia)

Se generaron 18 escenarios de temperatura y 18 gundo la calibración de la precipitación-escorrentía

escenarios de precipitación, de los cuales se selec- para las cuencas Condoriri, Tuni y Huayna Potosí.
cionaron seis combinaciones de temperatura y pre- La cuenca glaciar Zongo tiene un área de 3.5 km2
cipitación, excluyendo los más extremos. Con este con cobertura glaciar del 65% y el resto conformada
procedimiento se generaron las proyecciones climá- por afloramientos rocosos constituidos de granodio-
ticas para el futuro. rito impermeable, por tanto casi la totalidad de agua
de fusión transita por la estación limnimétrica ubica-
RESULTADOS da sobre una altitud de 4830 m (Leonardini, 2010).
Para el periodo de tiempo de calibración (2004-2010)
Calibración y validación del modelo los parámetros de calibración obtenidos fueron: T0
La estrategia de calibración para el área de estudio en 0.5 °C, aice en 350 mm mes-1 °C-1 y asonw en 550
consistió primero en la calibración del modulo gla- mm mes-1 °C-1. Para evaluar la bondad de ajuste
ciar grado-mes para la cuenca glaciar Zongo y se- del modelo se utilizó el índice de eficiencia de Nash-

28 Aqua-LAC - Vol. 5 - Nº. 2 - Sep. 2013

 Modelación hidrológica y escenarios de cambio climático en cuencas de suministro de agua de las ciudades La Paz y El Alto, Bolivia
Sutcliffe, que es una medida estadística normalizada datos. Para verificar su representación adecuada el
que determina la magnitud relativa de la varianza re- área modelada se comparó con el área digitalizada
sidual en comparación con la varianza de datos de sobre una imagen del satélite Landsat TM 5. En la
medición (Nash & Sutcliffe, 1970). Cuando el índice Figura 9 se puede observar dicha comparación en
de eficiencia es igual a 1 el ajuste es perfecto entre un diagrama de dispersión para el año 2005 y 2010,
los datos modelados y observados. El valor de índice donde los resultados tienden a alinearse a una lí-
de calibración obtenido es 0.6 que corresponde a un nea 1:1 indicando una buena correspondencia. Sin
nivel adecuado de calibración considerando las limi- embargo, considerando la resolución espacial de la
taciones en datos. Por otro lado, la desviación relati- imagen usada que es 30 m, pueden existir errores en
va de los caudales o el sesgo fue estimado en 3.7%, el cálculo del área digitalizada u observada. La evo-
demostrando una adecuada calibración. En la Figura lución del área glaciar indica que la cuenca Condoriri
8 se muestra la distribución mensual de caudales ca- ha perdido 36% de su superficie glaciar y la cuenca
librados, donde se evidencia que existe una adecua- Tuni el 47% entre 1995 y 2010.
da representación de los patrones temporales.

Tabla 2. Parámetro de calibración en las cuencas

Condoriri, Tuni y Huayna Potosí.
Vegetación
Parámetro Roca Bofedales
dispersa
Coeficiente de
0.95 0.9 0.9
cultivo, Kc
Capacidad de
almacenamiento
500 800 1000
en la zona de
raíces, Sw [mm]
Capacidad de
almacenamiento
500 500 500
de agua en la zona
profunda, Dw [mm]
Figura 8. Caudal observado y simulado medio men- Factor de
sual en la cuenca glaciar Zongo durante el periodo resistencia a la 5, 1 5, 1 5, 1
2004-2010 (Fuente: Elaboración propia) escorrentía, RRF
Conductividad de
zona de raíces, Ks 400, 50 400, 50 400, 50
[mm/mes]
Conductividad de
zona profunda, Kd 50 50 50
[mm/mes]
Dirección
preferencial de 0.9 0.7 0.3
flujo, f
Z1 [%] 10 10 10
Z2 [%] 30 30 30
Fuente: Elaboración propia

En el proceso de calibración de precipitación-esco-

rrentía se recomienda primero entender todos los pa-
rámetros que intervienen en Soil Moisture Model de
WEAP y segundo conceptualizar las caracterísitcas
hidrológicas particulares de las cuencas asociadas a
la cobertura y uso del suelo, así tambien a otras ca-
racterísticas fisiográficas. En el caso de las cuencas
Figura 9. Diagrama de dispersión de áreas glaciares en estudio se realizaron algunos ajustes en la Con-
observados y simulados para el año 2005 y 2010 ductividad en la Zona de Raíces (Ks) y Factor de Re-
(Fuente: Elaboración propia). sistencia a la Escorrentía (RRF). En el primer caso,
Ks, es un parámetro no linear que puede definirse
Los parámetros de calibración de Zongo fueron ex- en función de la precipitación ya que en estas cuen-
trapolados a los glaciares donde no se cuenta con cas el suelo al inicio del año hidrológico (septiembre)

Aqua-LAC - Vol. 5 - Nº. 2 - Sep. 2013 29

Marisa Escobar, Nilo Lima, David Purkey, David Yates y Laura Forni
se encuentra muy seco. Cuando existe precipitación modelo está basado en el modelo hidrológico FIES-
esta no escurre directamente sino que el escurri- TA (Mulligan & Burke, 2005; Bruijnzeel et al., 2011)
miento empieza una vez el suelo está lo suficiente- y procesos físicos; además no requiere calibración.
mente saturado para liberarlo. Entonces para la ca- Los componentes que utiliza son la distribución de la
libración de este parámetro se planteó la siguiente precipitación a través de la interacción con el viento,
hipótesis: menos conductividad en la época seca (50 radiación solar, evapotranspiración potencial, evapo-
mm/mes) y mayor conductividad en la húmeda (400 transpiración real y cubierta vegetal (http://www.po-
mm/mes). Asimismo para el RRF, cuando el suelo licysupport.org/aguaandes). El modelo opera sobre
está seco existe baja escorrentía (RRF=5), ya que una plataforma WEB y a una resolución espacial de
gran parte de la precipitación se infiltra y en la época 1 km2 o 1 ha. Cuando se procesa para alguna región
húmeda debido a que el suelo se encuentra satura- de la Cordillera de Los Andes el sistema se denomina
do existe mayor escorrentía (RRF=1). El umbral de AguAAndes, ahora si se aplica en otras partes del
precipitación calibrado para realizar estos ajustes en mundo se conoce como WaterWorld (Mulligan, 2012).
las cuencas fue 100 mm. Los restantes parámetros Los caudales modelados por AguAAndes correspon-
se pueden observar en la Tabla 2. den al periodo 1950-2000, mientras que la modela-
Los resultados de la calibración son aceptables. El ción en WEAP se realizó para el periodo 1995-2010.
parámetro de eficiencia Nash-Sutcliffe (Nash & Sut- Los resultados para algunas cuencas estudiadas
cliffe, 1970) en las tres cuencas oscila entre 0.6 y podemos ver en la Tabla 4. En el caso de cuencas
0.73, pero en Tuni hay una subestimación del 17.57% con datos de caudales observados (El Alto) los re-
según la medida del sesgo (Tabla 3). Los parámetros sultados de comparación son relativamente acepta-
de calibración del fueron extrapolados para las cuen- bles, al igual que en las cuencas potenciales; pero
cas sin caudales observados bajo los siguientes cri- en las cuencas de suministro de La Paz AguAAndes
terios: las cuencas son de régimen natural, es decir sobreestima los caudales. Sin embargo no es co-
no existe mucha intervención humana y por tanto la rrecto dicha comparación, ya que la oferta calculada
cobertura y uso de la tierra es prácticamente similar corresponde a distintos periodos de tiempo. A objeto
en todos; están sobre una altitud mayor a 4300 m; de validar se efectuó la prueba de U Mann Whitney,
tienen características geológicas similares, donde se que es un test no paramétrico que permite comparar
observan rocas sedimentarias del Silúrico y Ordoví- la distribución de dos poblaciones.
cico, batolitos del Triásico y depósitos de morrena El test considera que la hipótesis alterna (Ha) esta-
del Cuaternario (GEOBOL, 1995); contienen forma- blece que las distribuciones son distintas, mientras
ción de suelos in situ de textura franco arenoso y que la hipótesis nula (H0) que las poblaciones tienen
franco arcillo arenoso. Según los criterios anteriores, la misma distribución. En base a los resultados mos-
se consideró que era posible extrapolar la calibración trados en la Tabla 5, el valor U adoptado es 44.5 el
de las cuencas con información a aquellas cuencas cual es más grande que el nivel de significancia por
exploradas como nuevas fuentes. Aunque esto pue- tanto se rechaza la Ha y se acepta la H0, es decir que
de generar fuentes de incertidumbre en las predic- las poblaciones provienen de distribuciones idén-
ciones, esta incertidumbre no se ignora, sino que se ticas. Esta técnica de comparación puede ser una
asume de forma que los datos obtenidos constituyen forma de validación de la extrapolación de los pará-
estimados y no datos precisos. metros de calibración.

Tabla 3. Índices/métricas de la calibración del mo- Tabla 4: Comparación entre caudales modelados
delo hidrológico por WEAP y AguAAndes
Estación hidrométrica Nash-Sutcliffe Bias [%] AguAAndes WEAP
Fuente Ciudad Cuenca
Tuni 0.73 -17.57 [m3/s] [m3/s]
Condoriri 0.60 -0.44 El Alto Tuni 0.26 0.29

Huayna Potosí 0.65 6.42 Condoriri 0.33 0.36

Fuente: Elaboración propia Actual Huayna Potosí 0.73 0.70

La Paz Milluni 0.93 0.74
Comparación de caudales modelados por WEAP Chacaltaya 0.21 0.12
y AguAAndes
Kaluyo 0.95 0.43
El AguAAndes es un modelo sofisticado de balance
Incachaca 0.72 0.32
hídrico distribuido espacialmente que ha sido desa-
rrollo para cuencas con escases de datos y carac- Khara Khota 1.01 1.10
terísticas muy heterogéneas. Además provee herra- Potencial
El Alto Taypichaca 1.37 1.33
mientas para el análisis de escenarios de cambio
climático y uso de la tierra; asimismo es posible el Jachawaquiwiña 0.30 0.33
análisis de opciones de políticas o intervención. El Fuente: Elaboración propia

30 Aqua-LAC - Vol. 5 - Nº. 2 - Sep. 2013

 Modelación hidrológica y escenarios de cambio climático en cuencas de suministro de agua de las ciudades La Paz y El Alto, Bolivia
Tabla 5: Resultados del Test U Mann Whitney existen pérdidas por evaporación en los embalses.
Además la capacidad máxima de la infraestructura
Modelo N R U α=0.05 de conducción del embalse hasta la planta de tra-
AguAAndes 10 110.5 44.5 Se acepta H0 tamiento es de 0.968 m3/s. Es así que la cobertura
modelada es solo del 87%. El sistema está al límite
WEAP 10 99.5 55.5 Se rechaza Ha
de su capacidad de producción, por esta razón se
N: Tamaño de la muestra extrae agua del embalse Milluni (sistema Achachi-
R: Suma de rangos cala) en un total de 0.055 m3/s y ésta se inyecta al
U: Valor de Mann Whitney sistema El Alto.
Considerando la demanda actual, es necesario y ur-
Oferta y demanda de agua gente la inversión en nueva infraestructura hidráuli-
Para la modelación de la demanda se utilizó la base ca en otras cuencas planteadas en el Plan Maestro
de datos de EPSAS que incluye: población, consu- de Agua Potable del año 1994, como Khara Khota y
mo per cápita y volúmenes de embalses observados. Taypichaca, donde el caudal estimado fue de 1.10
Los datos finales de entrada al modelo se pueden y 1.33 m3/s respectivamente. En estas cuencas la
observar en la Tabla 6. La oferta de agua determina- presencia de glaciares es mucho más importante en
da con el modelo WEAP para las cuencas en estudio términos de área, y por tanto el aporte glaciar esta en
para el periodo 1995-2010 se muestra en la Tabla 7. el orden del 49 y 43% respectivamente.
En el caso de las cuencas de suministro El Alto la En las cuencas de suministro del sistema Pampa-
oferta oscila de 0.29 a 0.70 m3/s donde el aporte gla- hasi, Hampaturi e Incachaca, el caudal es de 0.49
ciar representa entre 20 a 28%. En la época de lluvia y 0.32 m3/s respectivamente, haciendo un total de
el caudal máximo medio está entre 0.77 y 1.83 m3/s; 0.81 m3/s (Tabla 7). El aporte glaciar modelado en
mientras en la época seca 0.04 y 0.11 m3/s (Figura la cuenca Hampaturi es del orden del 10%. La de-
10). El caudal de ingreso al embalse Tuni en el año manda modelada para el año 2010 fue de 0.74 m3/s.
2010 fue de 1.3 m3/s. La cobertura modelada considerando todas las pér-
didas en el sistema es del 100%.
Tabla 6: Datos para modelación de la demanda de En el sistema Achachicala la oferta de la cuenca
los sistemas de agua de la Paz y El Alto Milluni es de 0.74 m3/s. La demanda calculada para
el 2010 es 0.75 m3/s (Tabla 7). Además existe una
Sistema Ciudad Población Consumo Perdidas aducción en la cuenca baja del rio Kaluyo (Choque-
[hab] per cápita [%] yapu) con una capacidad máxima de 0.2 m3/s cerca
[m3/año]
a la planta de tratamiento, es así que se logra una
Achachicala 195,000 81 33 cobertura del 100%.
La Paz
Pampahasi 300,000 52 33
El Alto El Alto 1’150,000 19 38 Tabla 7: Oferta de agua en las cuencas de sumi-
nistro de agua actual y potencial de las ciudades La
Fuente: en base a base de datos EPSAS
Paz y El Alto para el periodo 1995-2010.
Abastece Caudal Aporte
Sistema la ciudad Cuenca total Glaciar
de [m3/s] [m3/s]
0.07
Tuni 0.29
(24%)
0.10
El Alto El Alto Condoriri 0.36
(28%)
0.14
Huayna Potosí 0.70
(20%)
0.06
Achachicala Milluni 0.74
(8%)
Incachaca 0.32
La Paz
Pampahasi 0.05
Hampaturi 0.49
(10%)
Figura 10. Caudal de aporte de las cuencas de 0.54
suministro de agua de las ciudades de La Paz y El Khara Khota 1.10
(49%)
Alto (Fuente: Elaboración propia) 0.57
Nuevas Taypichaca 1.33
(43%)
cuencas
La demanda de agua modelada para el sistema El El Alto Janchallani 0.14
Alto en el año 2010 fue del 1.11 m3/s, del cual el 38% 0.14
Jachawaquiwiña 0.33
corresponde a pérdidas en el sistema. Si bien la ofer- (42 %)
ta en las cuencas es mayor que la demanda, pero Fuente: Elaboración propia

Aqua-LAC - Vol. 5 - Nº. 2 - Sep. 2013 31

Marisa Escobar, Nilo Lima, David Purkey, David Yates y Laura Forni
Impacto del cambio climático en el En la Figura 11 se muestra la evolución del retro-
abastecimiento de agua ceso glaciar de la cuenca Tuni bajo escenarios de
cambio climático, donde los resultados indican que
El análisis del impacto del cambio climático en el
los glaciares en dicha cuenca desaparecerían en los
abastecimiento de la demanda se realizó para un pe-
próximos 30 años. Esto definitivamente podría tener
riodo de tiempo a mediano plazo 2011-2026 y a largo
implicaciones en el abastecimiento de agua en el sis-
plazo 2035-2050. Los cambios de precipitación y tem-
tema El Alto, donde la demanda modelada hasta el
peratura son variables para cada estación base usa-
año 2026 fue en 1.29 m3/s y hasta el año 2050 en
da en WEAP. En la Tabla 8 se muestran los cambios
1.55 m3/s. Con la actual infraestructura la cobertura
esperados en el clima respecto al periodo histórico
de la demanda para el escenario de clima pesimista
1995-2010, en el caso de la temperatura a mediano
a mediano y largo plazo seria del 71% y 51% respec-
plazo se espera un incremento hasta 1.1 °C y a largo
tivamente; y para el caso del clima futuro optimista
plazo 1.8 °C. Respecto a la precipitación, en la esta-
75% y 62% para ambos casos.
ción base Tuni se esperan reducciones hasta 11.4% y
8.7% para los respecticos periodos de análisis.
Tabla 9. Cambios probables en la oferta de agua
para el periodo 2011-2026 respecto al periodo de
Tabla 8: Cambios de precipitación (estación base referencia 1995-2010.
Tuni) y temperatura a mediano y largo plazo.
Escenario de Cambio Climático
Periodo Cuenca
1 2 3 4 5 6
Clima 2011-2026 2035-2050 Tuni -21% -17% -24% -24% -21% -28%
 P [%]  T[°C] P [%] T [°C] Condoriri -19% -14% -19% -19% -17% -22%
1 +0.8 +0.0 +6.3 +0.8 Huayna
2 +1.1 +0.4 +8.7 +1.1 Potosi -14% -9% -14% -16% -13% -19%

3 -4.9 +0.5 +3.1 +1.2 Milluni +5% 0% -3% -14% -12% -14%
Incachaca +3% -6% -9% -16% -16% -28%
4 -5.5 +0.6 -0.4 +1.4
Hampaturi +6% 0% -6% -12% -16% -22%
5 -4.8 +0.8 -8.1 +1.7
Khara Khota -28% -15% -15% -13% -7% -7%
6 -11.4 +1.1 -8.7 +1.8
Taypichaca -27% -17% -17% -17% -13% -14%
Fuente: Elaboración propia
Fuente: Elaboración propia
Nota: la descripción de los seis escenarios de clima se
Respecto a la demanda de agua se realizaron pro- puede ver en la Tabla 8.
yecciones de población en base a los datos históri-
cos y tendencias del Instituto Nacional de Estadística
(INE). Donde para El Alto se estimó una población de Tabla 10. Cambios probables en la oferta de agua
1’600.000, para Pampahasi 360.000 y para Achachi- para el periodo 2035-2050 respecto al periodo de
cala 197.000 hasta el año 2050. Respecto al consu- referencia 1995-2010.
mo per cápita y pérdidas se consideraron los mismos
Escenario de Cambio Climático
datos del periodo 1995-2010. Los cambios probables Cuenca
en la oferta de agua por el cambio climático en las 1 2 3 4 5 6
cuencas de suministro de agua se muestran en la Tuni -14% -14% -21% -24% -38% -31%
Tabla 9 y Tabla 10.
Condoriri -14% -17% -22% -25% -36% -33%
Los resultados indican que las reducciones de cau-
dales son más críticos en las cuencas de suministro Huayna Potosi -6% -6% -14% -17% -30% -27%
del sistema El Alto. Estos cambios marcados están Milluni +9% +24% 11% -14% -15% -15%
directamente relacionados con el aporte de los gla-
ciares, ya que las cuencas son pequeñas en com- Incachaca +9% +3% +0% -13% -25% -28%
paración de otras cuencas de suministro como Inca- Hampaturi +8% +8% +6% -10% -20% -27%
chaca, Hampaturi y Milluni. Para el periodo histórico
según la modelación hubo un retroceso acelerado Khara Khota -11% -18% -22% -28% -35% -33%
especialmente en las bandas de altitud 4700-5000 Taypichaca -14% -21% -25% -32% -38% -35%
m y 5000-5300 m, según la reclasificación de la alti- Fuente: Elaboración propia
tud para la modelación. Asimismo el año hidrológico Nota: la descripción de los seis escenarios de clima se pu-
1997-1998 fueron los más cálidos de los últimos seis ede ver en la Tabla 8.
décadas (Vuille & Bradley, 2000). Considerando que
el periodo histórico de modelación es simplemente
dieciséis años, el aporte glaciar es mayor en compa- Con base en la demanda modelada tanto para los
ración a los aportes modelados hacia el futuro, y por sistemas de Pampahasi y Achachicala y los cambios
ende se dan cambios notables en los caudales. esperados en los caudales de las cuencas, la cober-

32 Aqua-LAC - Vol. 5 - Nº. 2 - Sep. 2013

 Modelación hidrológica y escenarios de cambio climático en cuencas de suministro de agua de las ciudades La Paz y El Alto, Bolivia
tura de la demanda para el escenario de clima futuro REFERENCIAS BIBLIOGRÁFICAS
pesimista a mediano plazo sería de 97% y 78% para
ambos casos. A largo plazo las coberturas se reduci- AguAAndes. http://www.policysupport.org/aguaandes
rían a 94% y 68% respectivamente. Bahr, D.B., M.F. Meier and S.D. Peckham. 1997.
The physical basis for glacier volume-area scaling.
Journal of Geophysical Research 102 (B9): 20355–
20362.
Bates, B.C., Z.W. Kundzewicz, S. Wu & J.P. Palu-
tikof, Eds. 2008. Climate Change and Water. Techni-
cal Paper of the Intergovernmental Panel on Climate
Change, IPCC Secretariat, Geneva, 210 pp.
Bruijnzeel, L. A., M. Mulligan & F.S.Scatena. 2005.
Hydrometeorology of tropical montane cloud forests:
emerging patterns. Hydrological Processes 25: 465–
498.
Figura 11. Retroceso glaciar en la cuenca
Tuni para diferentes escenarios de clima futuro Brunet-Moret, Y.1979. Homogénéisation des precipi-
(Fuente: Elaboración propia) tations. Cahier ORSTOM, Série Hydrologie 16: 3–4
Nota: la descripción de los seis escenarios de clima se Centro de Cambio Global, Universidad de Chile &
puede ver en la Tabla 8 Stockholm Environment Institute. 2009. Guía Meto-
dológica – Modelación Hidrológica y de Recursos Hí-
CONCLUSIONES dricos Con El Modelo WEAP. Santiago, Chile, y So-
merville, MA, EEUU. http://weap21.org/downloads/
En esta investigación se desarrolló una herramienta Guia_modelacion_WEAP_Espanol.pdf
que permite la modelación de los recursos hídricos
incluyendo incertidumbres como cambio climático Condom, T., M. Escobar, D. Purkey, J.C. Pouget, W.
que puede ser útil en la planeación. Sin embargo, Suarez, C. Ramos, J. Apaestegui, A. Tacsi & J. Go-
se deben hacer algunas consideraciones. El modulo mez. (2012): Simulating the implications of glaciers’
glaciar usado es simple debido a las limitaciones en retreat for water management: a case study in the
datos, sería adecuado usar un modelo basado en el Rio Santa basin, Peru. Water International 37 (4):
balance energético. 442-459
La modelación hidrológica en zonas de montaña Francou, B., M. Vuille, P. Wagnon, J. Mendoza & J.E.
como el área estudiada es bastante compleja, de- Sicart. 2003. Tropical climate change recorded by a
bido a la escasez de datos y a la alta variabilidad glacier in the central Andes during the last decades
climática. En las cuencas de suministro de La Paz of the twentieth century: Chacaltaya, Bolivia, 16°S.
y El Alto, existen escasos datos climáticos, y en el Journal of Geophysical Research: Atmospheres
caso de los glaciares fue necesario extrapolar datos 108(D5). DOI:10.1029/2002JD002959.
del glaciar Zongo. Con estas consideraciones los re-
GEOBOL. 1995. Carta Geologica Hoja 5945 - Milluni.
sultados presentados de los aportes glaciares deben
La Paz, Bolivia.
ser usados como referencia y no como valores exac-
tos. En el caso de la extrapolación de los parámetros Hiez, G., 1977. L’homogénéité des données pluvio-
de calibración del modelo hidrológico a las posibles métriques. Cahier ORTSTOM, Série Hydrologie 14:
cuencas de suministro, fue validado con el balance 129–172.
hídrico en AguAAndes. Los resultados modelados
IPCC. 2007. Climate Change 2007: Impacts, Adapta-
tanto en WEAP como AguAAndes son aceptables.
tion and Vulnerability. Contribution of Working Group
De esta manera esta herramienta puede ser útil en
II to the Fourth Assessment Report of the Intergov-
cuencas donde no se cuenta con datos suficientes.
ernmental Panel on Climate Change. M. L. Parry, O.
El sistema El Alto en el futuro podría sufrir serios pro- F. Canziani, J. P. Palutikof, P. J. van der Linden, y C.
blemas en el abastecimiento de debido al incremento E. Hanson (eds.). Cambridge University Press, Cam-
de la demanda y la reducción de la oferta de agua bridge, UK.
por los cambios esperados en el clima futuro. Asimis-
IPCC. 2013: Summary for Policymakers. In: Climate Change
mo, el sistema Pampahasi también podría presentar
2013: The Physical Science Basis. Contribution of Working Group
problemas pero menos preocupantes que en el caso
I to the Fifth Assessment Report of the Intergovernmental Panel on
de El Alto. Finalmente en el caso del sistema Acha-
Climate Change [Stocker, T.F., D. Qin, G.-K. Plattner, M. Tignor, S.
chicala, si bien las cuencas podrían sufrir cambios en
K. Allen, J. Boschung, A. Nauels, Y. Xia, V. Bex and P.M. Midgley
la disponibilidad de agua ésta no presentaría proble-
(eds.)]. Cambridge University Press, Cambridge, United Kingdom
mas en la cobertura de la demanda.
and New York, NY, USA.

Aqua-LAC - Vol. 5 - Nº. 2 - Sep. 2013 33

 Marisa Escobar, Nilo Lima, David Purkey, David Yates y Laura Forni
Leonardini G. 2010. Análisis mensual de los balan- Soruco, A., C. Vincent, B. Francou & J.F. Gonzalez.
ces de masa glaciológico e hidrológico del glaciar 2009. Glacier decline between 1963 and 2006 in the
tropical Zongo (Bolivia). Tesis de licenciatura. Uni- Cordillera Real, Bolivia. Geophysical Research Let-
versidad Mayor de San Andrés. Bolivia. ters 36: L03502. DOI:10.1029/2008GL036238
Mehta, V., M. Escobar & D. Decondapa. 2010. UN- Tebaldi, C., R. Smith, D. Nychka & L. Mearns. 2005.
DERSTANDING GLACIATED HIMALAYAN HEAD- Quantifying Uncertainty in Projections of Regional
WATERS - An Integrated Climate-Energy-Water As- Climate Change: A Bayesian Approach to the Analy-
sessment Framework. White Paper. sis of Multimodel Ensembles. J. Climate 18: 1524–
1540. DOI: http://dx.doi.org/10.1175/JCLI3363.1
Molina, F. & R. Pereira. 2005. Crecimiento demográ-
fico y el problema del agua en El Alto. La Prensa. 310 Vauchel, P., 2005. HYDRACCESS: Software de pro-
ed. Bolivia. cesamiento de datos hidro-meteorologico. Disponi-
ble en: http://www.mpl.ird.fr/hybam/outils/logiciels.
Mulligan, M. 2012. WaterWorld: a self-parameteri-
htm.
sing, physically-based model for application in data-
poor but problem-rich environments globally. Hydro- Vuille, M. & R. Bradley. 2000. Mean annual tempera-
logy Research, in press. DOI: 10.2166/nh.2012.217 ture trends and their vertical structure in the tropical
Andes. Geophysical Research Letters 27 (23): 3885-
Mulligan, M. & S. M. Burke. 2005. FIESTA: Cloud wa-
3888.
ter interception for the Enhancement of Streamflow in
Tropical Areas. Report to UK DfID, 174 pp. Yates, D., J. Sieber, D. Purkey & A. Huber-Lee. 2005a.
WEAP21—A Demand-, Priority-, and Preference-
Nash, J. E. & J.V. Sutcliffe. 1970. River flow forecast-
Driven Water Planning Model. Water International 30
ing through conceptual models part I — A discussion
(4): 487–500. DOI:10.1080/02508060508691893.
of principles. Journal of Hydrology 10(3): 282–90.
DOI: 10.1016/0022-1694(70)90255-6. Yates, D., D. Purkey, J. Sieber, A. Huber-Lee & H.
Galbraith. 2005b. WEAP 21--A Demand, Priority,
Pouget, J. C. 2011. Propuesta del modelo Ice KISS
and Preference-Driven Water Planning Model: Part
como un componente de un sistema de apoyo a la
2, Aiding Freshwater Ecosystem Service Evaluation.
planificación de los recursos hídricos. Quito-Ecuador.
Water International 30 (4): 501–2.
Ramírez, R., E. 2008. Impactos del cambio climá-
tico y gestión del agua sobre la disponibilidad de
recursos hídricos para las ciudades de La Paz y El AGRADECIMIENTOS
Alto. Revista Virtual REDESMA 2 (3)Ramirez, E., Este trabajo de investigación forma parte del apoyo
B. Francou, P. Ribstein, M. Descloitres, R. Guérin, técnico al “Programa de Resiliencia Climática para
J. Mendoza, R. Gallaire, B. Pouyaud & E. Jordan. los Sistemas de Agua y Saneamiento de las Áreas
2001. Small glaciers disappearing in the tropical Metropolitanas de La Paz y El Alto” (PPCR/SPCR).
Andes: a case-study in Bolivia: Glaciar Chacaltaya Agradecemos el apoyo del personal del Banco In-
(16° S). Journal of Glaciology 47(157): 187–94. teramericano de Desarrollo (BID) en Washington y
DOI:10.3189/172756501781832214. Bolivia, especialmente a Alfred Grunwaldt, Alejandro
Deeb, Ana Iju, Martin Kerres, y Prem Jai Vidaurre de
la Riva.

34 Aqua-LAC - Vol. 5 - Nº. 2 - Sep. 2013

View publication stats