RETROCESO GLACIAR E IMPLICANCIAS PARA LA DISPONIBILIDAD HÍDRICA: CASO DE LOS RIOS LECHERIA Y

COLORADO (LIMATAMBO-CUSCO-PERÚ).
Feliciano Escobedo Silva
Escuela académico Profesional de Ingeniería de Minas, Universidad Nacional Micaela Bastidas de Apurímac
Abancay, Apurímac, Perú
Walquer Huacani Calsin
Escuela académico Profesional de Ingeniería de Minas, Universidad Nacional Micaela Bastidas de Apurímac
Abancay, Apurímac, Perú
Darío Dante Sánchez Castillo
Escuela académico Profesional de Ingeniería de Minas, Universidad Nacional Micaela Bastidas de Apurímac
Abancay, Apurímac, Perú
Pablo Rubén Zuloaga Candia
Escuela académico Profesional de Ingeniería de Minas, Universidad Nacional Micaela Bastidas de Apurímac
Abancay, Apurímac, Perú
Elmer Elio Calizaya Llatasi.
Escuela académico Profesional de Ingeniería Topográfica y Agrimensura, Universidad Nacional del Altiplano
Puno Perú

RESUMEN El forzamiento radiativo por el incremento de la temperatura

ambiental, está generando la sublimación directa en los glaciares
La presente investigación es sobre el retroceso glaciar de los [ CITATION Sua15 \l 10250 ],son una antesala a escala
Nevados Oq’obamba y Moyoq (Ubicado al Este de la cordillera fractal de lo que sucederá con los glaciares de mayor altitud, tal
de Vilcabamba) y sus implicancias en la disponibilidad hídrica como está sucediendo con el Nevado Salkantay, cuyo ELA
de las cuencas de los ríos Lechería y Colorado (Limatambo), en (Línea de equilibrio altitudinal) se incrementa en forma
el periodo de 2000 al 2015, ha permitido determinar que los
ascendente cada año [CITATION Mon14 \l 10250 ]
referidos nevados han perdido 24.75 Km2 de su superficie glaciar
A nivel mundial, la disminución de las extensiones de glaciares
que corresponde al 98.21 %. El objetivo principal de la
y nieve de la Criosfera concuerda también con el calentamiento
investigación fue la determinación de la relación de las
global. Datos satelitales obtenidos desde 1978 indican que el
desglaciaciones con la disponibilidad hídrica de las cuencas
promedio anual de la extensión de los hielos marinos árticos ha
directas receptoras que fueron monitoreadas en la estación
disminuido en un 2,7 km2 (entre 2,1 y 3,3) % por decenio; los
hidrométrica del Puente Kunyaq durante el periodo 2000 al
glaciares de montaña y la cubierta de nieve han disminuido en
2015. Ambas variables fueron sometidas al análisis de
ambos hemisferios.
correlación, obteniéndose el coeficiente de correlación de
Pearson r=0.467, deduciéndose que no existe una relación En el continente americano desde Alaska, Parque Nacional de
significativa entre el retroceso glaciar y la disponibilidad hídrica los Glaciares en EE. UU, Volcanes Mexicanos, glaciares
ya que el volumen de agua que no llega a ser medido en la tropicales Andinos desde Venezuela hasta Bolivia y los glaciares
estación hidrométrica, por lo tanto se están sublimando tal como patagónicos comprendidos entre Chile y Argentina en
se pudo observar en las visitas de campo, evidenciándose que los Sudamérica también están en un franco proceso de retroceso.
últimos glaciares se están fragmentando, encapsulando y no [ CITATION Rau05 \l 10250 ]. En el Perú se encuentra el
aportan torrente visible alguno. Los glaciares de los nevados de 71% de los glaciares tropicales del mundo distribuidos en 19
menores altitudes inferiores a 5150 msnm., están ubicados al cordilleras, En el periodo 1962 - 1970 se contabilizó 3044
extremo Este de la cordillera de Vilcabamba son mas glaciares que cubrían 2041.82 km2, en el año 2009 se contabilizó
susceptibles al incremento de la temperatura debido al 2679 glaciares disminuyendo su área a 1298.59 Km 2 (pérdida
calentamiento global. La disminución de la disponibilidad del 42% en 40 años).
hídrica de los ríos Lechería y Colorado en un futuro no muy
lejano generará problemas de estrés hídrico debido a la demanda La cordillera de Vilcabamba, contaba con 355 glaciares
hídrica para cultivos de palta con fines de agroexportación en el comprendidos en una extensión de 169.15 km2, sin embargo, en
valle de Limatambo. el año 2014 se redujo dramáticamente y una extensión de 15.53
Km2 habiendo perdido 22.21 Km2 equivalente al 58.85 %.
Palabras clave: Retroceso glaciar, disponibilidad hídrica y (Autoridad Nacional del Agua (ANA, 2007).
sublimación de glaciares.
1. INTRODUCCIÓN Los Nevados Oq’obamba y Moyoq constituían, la reserva
hidrológica ya que la fusión de nieve alimentaban a los ríos
La motivación del presente estudio es entender el fenómeno del Lechería y Colorado en el distrito de Limatambo (Lado Sur) y
retroceso glaciar y su relación con la disponibilidad hídrica a sus del Riachuelo Silque en Ollantaytambo (Lado Norte), debido a
cuencas receptoras. directas: Llamawasi, Yawarmaki. que estos nevados se encuentran a altitudes muy próximas a la
Q’elwaq’ocha y Ch’akilq’asa. ELA (Altitud de la línea de equilibrio) determinadas en el
glaciar Inkachiriasq’a (Nevado Salkantay) en 5030 msnm., en el
En cuanto al estado de arte, el presente estudio puntualiza el año 2012.[CITATION Mon14 \l 10250 ] El retroceso glaciar
análisis de glaciares de menor altitud ya que la mayoría de las iniciado en los finales de Pequeña Edad del Hielo (PEH) en la
investigaciones tratan sobre los nevados más representativos de cordillera de Vilcabamba (Licciardi et al, 2009) fue gradual
las cadenas montañosas, los glaciares de menor altitud son los intensificándose en la década del 1970 al 2007 [ CITATION
más susceptibles y constituyen los mejores indicadores del
Mar09 \l 10250 ].
cambio climático identificándose 2341 glaciares que
corresponde a un (87.38%), (superficies menores o iguales a 1 Todos creemos por sentido común que a mayor retroceso mayor
km2, y 338 glaciares mayores a 1 km 2.)[ CITATION Mel \l aporte hídrico[ CITATION Bar \l 10250 ]; Sin embargo, los
10250 ] resultados del presente proyecto de investigación indican que en

1
la actualidad los glaciares ya no se derriten ni aportan agua a las relación a los elementos vacíos suelo-roca. [ CITATION Pel10
cuencas receptoras; caudales de los torrentes glaciares han \l 10250 ]
dejado de discurrir, en similitud a lo que está ocurriendo en el
Se realizó el proceso de clasificación no supervisada de
glaciar Julcal Chile[ CITATION Mel \l 10250 ].Como
imágenes satelitales mediante las correcciones: geométrica,
medida de solución a este fenómeno se propone iniciar políticas
radiométrica y la atmosférica, asimismo se ha georeferenciado
de adecuación al inevitable cambio climático y hacer conocer a
cada imagen con coordenadas geográficas Latitud/Longitud con
la comunidad internacional sobre el retroceso acelerado de los
resolución espacial de 30 metros y en el Datum WGS84.
glaciares en la cordillera de Vilcabamba en la región Cusco-
Perú. Se consideró como cubierta nival celda que sean superiores a un
2. ÁREA DE ESTUDIO valor de 0,4 µm, este valor se considera como una orientación,
puesto que el valor del Índice Diferencial Normalizado de Nieve
Específicamente los Nevados Oq’obamba (5050 msnm) y Nevado (NDSI), depende del ángulo de elevación solar y las condiciones
Moyoq (5150 msnm), ocupan el extremo Este de la cordillera de atmosféricas. [ CITATION ting \l 10250 ]. Para el análisis
Vilcabamba, de la cadena oriental Andina comprendidos entre multitemporal se calculó su área en (km 2) para toda la superficie
latitud: 13° 33´ 47.17” - 13° 21’28.4” Sur y longitud 72°34 cubierta por glaciares y nevadas de los cuatro glaciares en el
´30.71”- 72°32´43.88” Oeste. periodo y cumpliendo las recomendaciones de Dozier, se utilizó la
formula conocida:

( B 2−B 5)
NDSI= Ec.
( B 2+ B5)
(1)

Con lo que se obtuvo resultados y cálculo la tasa se retroceso

glaciar y la tendencia del retroceso glaciar.

Tabla 1. Áreas de la cubertura en el periodo 2000-2015.

Año Sensor Área (Km2) Tasa

Variacio Retroceso
    Neta
n (%)
Fig.1 Ubicación de los nevados Oq’obamba y Moyoq
2000-2001 4 TM 1.398 - -
En el área se han verificado la existencia de lagunas de origen
glaciar, algunas con potencialidades de convertirse en Glofs 2001-2002 4 TM 5.148 -3.750 268.240
(Flujos por el desborde violento de una laguna glaciar) con
2002-2003 5 TM 0.390 4.760 -92.420
incrementos de volumen de agua que evidencian el fuerte
retroceso glaciar en las últimas décadas [ CITATION Gua16 \l 2003-2004 5 TM 2.609 -2.220 568.970
10250 ]. 2004-2005 5 TM 0.088 2.520 -96.630
2005-2006 8 TM 5.857 -5.770 6,555.680
3. METODOLOGÍA
2006-2007 5 TM 0.025 5.830 -99.570
Se planteo como objetivo el análisis por correlación de la data de
2007-2008 5 TM 0.124 -0.100 396.000
la desglaciación (Variable independiente) y los caudales
generados (Variable dependiente) en el periodo 2000-2015, de 2008-2009 5 TM 0.901 -0.780 626.610
cuyo resultado se obtendría el valor de estadímetro que 2009-2010 5 TM 0.288 0.610 -68.040
permitiría la interpretación de resultados.
2010-2011 5 TM 3.465 -3.180 1,103.130
3.1 Determinación del Índice normalizado de diferencia de 2011-2012 5 TM 8.465 -5.000 144.300
nieve (NDSI)
2012-2013 8 OLÍ 0.041 8.420 -99.520
Para la obtención del Índice Diferencial Normalizado de Nieve
2013-2014 8 OLÍ 0.915 -0.870 2,131.710
(NDSI), se realizó la descarga de las imágenes satelitales
Landsat de la United States Geological Survey (USGS) mediante 2014-2015 8 OLÍ 0.028 0.890 -96.940
la: www.usgs.gov del periodo 2000 al 2015, se consideró las Fuente: Imágenes de satélite Landsat.
imágenes satelitales de los meses de Mayo a Julio debido a la
mayor ablación en los referidos meses [ CITATION Hal10 \l 3.2 Tendencia del Retroceso y las precipitaciones nivales
10250 ].En el proceso de descarga de las imágenes de satélite
Landsat, obtenidos con sensores MSS, TM y OLI, se ha La tendencia de secuencias de las precipitaciones nivales a los
considerado imágenes libres de nube, y que tienen la capacidad glaciares y cumbres temporales de nieve es negativa, indica que
de distinguir altas reflectancia de los glaciares comparadas con cada año está disminuyendo 0.05 Km 2. Esta retroalimentación es
las zonas circundantes no nevadas, vegetación y agua, siendo el insuficiente para alcanzar un estado de equilibrio que permita la
objetivo principal determinar la reflectancia homogénea acumulación, por lo que se encuentra en permanente estado de
circundantes no nevadas, vegetación y agua, y determinar la ablación [ CITATION Viu1 \l 10250 ]
reflectancia homogénea (Hielo/Nieve) con gran contraste en Los parámetros obtenidos fueron:

2
 ∑ y ∑ x 2−∑ x ∑ xy rio Vilcanota) y actualmente se encuentran a finales de su
a= =2.37 existencia.
N ∑ x 2−¿ ¿
En el lado derecho dos cuerpos permanecen aún anclados a las
∑ y ∑ x 2−∑ x ∑ xy concavidades topográficas de la cumbre y por estar ubicados
b= =-0.05 Ec. (2)
N ∑ x 2−¿ ¿ en el flanco sur de la cumbre contrarios a los rayos solares en
horas de la tarde. Estos glaciares en épocas pasadas aportaban
Siendo la recta resultante: Y=2.37- 0.05 X sus torrentes directos la quebrada Yawarmaki (cuenca del río
Apurímac). Continúan ocupando una superficie de 0.108 km2.
6
ärea (mk2)

0
1 3 5 7 9 1 3 5
0 0 0 0 0 1 1 1
0 0 0 0 0 0 0 0
-2 -2 -2 -2 -2 -2 -2 -2
0 2 04 06 8 0 12 4
00 00 0 0 00 01 0 01
2 2 2 2 2 2 2 2
Periodo(Años)

Fig.2. Variaciones del retroceso glaciar periodo 2000-2015.

El análisis de la tendencia del retroceso de los cuatro glaciares

analizados y de las precipitaciones nivales dentro de las áreas de
la cuenca consideradas, indican una tendencia negativa de
menos 0.05 Km2/año, esta tendencia permite pronosticar los Fig.3.Glaciar Yawarmaki en el Nevado Oq’obamba.
nevados Oq’obamba y Moyoq perderán toda su cobertura glaciar
en el año 2030. La variación de las áreas de retroceso y la Sin embargo, todos estos cuerpos no aportan flujos superficiales,
frecuencia de nevadas en las áreas post retroceso en el periodo sus causes directos donde años atrás eran visibles permanecen
2000-2015 se observa que e l índice Oceánico El Niño (ONI) , secos. Sin embargo, no se descarta flujos subsuperficiales,
es uno de los más utilizadas donde se distinguen los fenómenos quedando pendiente mayores investigaciones para comprobar la
del Niño –Oscilación (ENOS y la “Niña”), se ha verificado la referida hipótesis.
existencia de cierta coincidencia del proceso de las
precipitaciones nivales que han retroalimentado a los glaciares
de los nevados Oq’obamba y Moyoq , con los fenómenos del
“Niño” y la “Niña” ,indican que durante los años de predominio
del “niño” disminuye las precipitaciones nivales y se incrementa
cuando ocurre la “Niña” de lo cual se desprende que durante la
“Niña” se recargan los glaciares y cumbres temporales de Nieve
y durante el “Niño” disminuyen.

Debido a que los glaciares analizados están a altitudes

inferiores al ELA[ CITATION Viu1 \l 10250 ] presentan
un comportamiento caótico y su retroalimentación se da Fig.4.Quebrada Yawarmaki sin flujo alguno.
únicamente durante las nevadas esporádicas y con mayor
incidencia durante la presencia del fenómeno de la “Niña”. El termino quechua “Yawarmaki” (mano de sangre) hace
referencia, que antiguamente los torrentes discurrían por la
quebrada del mismo nombre sedimentos ferrosos de color
3.3 Cálculo de la superficies glaciares existentes al año 2015. rojo(óxidos) y en la actualidad permanecen secos.
A medida que avanza el retroceso glaciar, se observa el progreso
Como se puede evidenciar mediante imágenes satelitales del
de gramíneas, musgos y líquenes, en las áreas desglaciadas.
Google Earth los glaciares de nevados Oq’obamba y Moyoq
se han reducido dramáticamente y continúan fragmentándose.

Durante las visitas de campo se observó que el nevado

Oq’obamba (5050 msnm), se ha fragmentado en cuatro
cuerpos encapsulados dos ubicados en el lado izquierdo, y los
restantes en el lado derecho. Se están subdividiendo en
fragmentos cada vez menores en definitiva se ha roto el
equilibrio entre la acumulación (Ganancia de la masa glaciar)
y la ablación (pérdida de la masa glaciar); Las escasas
retroalimentaciones de las precipitaciones nivales no son
suficientes para recuperar su equilibrio. [CITATION
Mar08 \n \l 10250 ]
En el lado izquierdo durante las épocas pasadas descargaban Fig. 5. Los nevados Oq’obamba (Color Verde fosforescente) y
sus torrentes a los orígenes del rio Pampaqahuana (cuenca del Moyoq (Color rojo) fines del año 2015. Fuente: Google Earth

3
 histórico.
CURVA DE CALIBRACION DE CAUDALES EN EL PUENTE
El nevado Moyoq, como se puede observar en la fig.5 (color rojo), para el KUNYAQ
año 2015, el glaciar Ch’akilq’asa mantiene únicamente 0.003 km2 (0.33 600
Has) ha perdido gran parte de su cubertura glaciar en las áreas
pertenecientes a la cuenca del río Apurímac, los restantes (0.493 k m2) se 500 f(x) = 41.65 x^1.77

CAUDL (m3/seg)
ubican en la cuenca del rio Vilcanota. El análisis del Google Earth histórico 400 R² = 0.99
indica que el glaciar Q’elwaq’ocha perdió toda su superficie glaciar en 300
2006.
200
100
0
0.5 1 1.5 2 2.5 3 3.5 4 4.5
Tabla 2. Área y perímetros de los fragmentos de glaciares de los nevados
Oq’obamba y Moyoq al año 2015. Altura de mira (m)

Nevado Moyoq Nevado Oq'obamba  

Fracción Perímetro Área Fracción Perímetro Área Total Fig.6. Curva de calibración obtenida de las campañas de aforo en la
estación hidrométrica del puente Kunyaq.
  (Km) (Km2)   (Km) (Km2) (Km2)
Una vez determinado el caudal en el puente Kunyaq se calcularon,
1 3.40 0.390 1 1.16 0.064  
los caudales en el rio Alto Apurímac, rio Bermejo también previa
2 1.35 0.100 2 0.96 0.044   discriminación del rio Blanco que colecta las aguas provenientes del
3 0.26 0.003 3 0.33 0.008  
nevado Salkantay y alrededores; Con la misma lógica se pudo
establecer los caudales de los ríos Colorado y Lechería aplicándose
  4.14   4 0.68 0.021   la Ley de extrapolación de áreas y caudales. Una vez precisado los
Sub total   0.493   0.137 0.630 caudales en los ríos Colorado y Lechería se calculó los caudales de las
microcuencas directamente colectoras de los glaciares Oq’obamba y Moyoq:
Llamawasi, Yawarmaki, Q’elwaq’ocha y Ch’akil. (Ver tabla 5).
Durante el periodo de estudio, los nevados Oq’obamba y Moyoq se
han reducido a 0.630 Km2, de los cuales el glaciar Llamawasi, ha
perdido su superficie glaciar mucho antes del inicio de la
investigación.

3.4 Determinación de la ecuación de calibración de Curva

Siendo la ecuación de calibración de Curva, la representación

gráfica de las variaciones de las descargas(Q) en la sección
transversal de un río, en función de las lecturas de mira (H) que se
produce teóricamente en un régimen permanente [ CITATION
Fig. 7. Estación hidrométrica del puente Kunyaq en el Rio Apurímac (Km 100
16Un \l 10250 ], se ha considerado el registro de las descargas Carretera Cusco-Abancay)
diarias durante los 15 años en el periodo 2000-2015, anotadas en la
estación hidrométrica del puente Kunyaq, mediante un Limnímetro,
Las campañas de aforos son realizadas por personal técnico
dado las dificultades de una medición directa y continua del caudal
especializado del servicio nacional de meteorología e Hidrología
de un río.
de Perú, (Senamhi) región Cusco.
La toma de información del limnímetro se realizan tres veces por día
Tabla 3. Cálculo de caudales de los ríos Colorado y Lechería.
y se contrasta con campañas de aforos periódicos para lo cual se
utilizan equipos de última tecnología como Q-Liner y procesamiento
Cálculo de caudales de las cuencas relacionadas a partir de la
en campo mediante software QTT. (Quantum technology tollbox).
calibración de curva en la estación hidrométrica del puente
La Ecuación de calibración de curva obtenida en la estación Kunyaq
Periodo
hidrométrica del puente Kunyaq, para el rio Apurímac es: Caudales (m3/seg) Ríos
(Años)
Y=41.647*H1.7652 Ec. (3) Alto Bermej Colorad Lecherí
 
Apurímac o o a
Donde: Y=Caudal del rio Apurímac (m3/s) y H=Altura del 2000-2001 203.07 4.2053 2.5482 0.5415
limnímetro, en una sección determinada (m). 2001-2002 149.92 3.1046 1.8812 0.3997
2002-2003 203.07 4.2053 2.5482 0.5415
2003-2004 152.95 3.1674 1.9193 0.4078
2004-2005 113.5 2.3504 1.4242 0.3026
2005-2006 197.63 4.0927 2.4799 0.5269
2006-2007 181.45 3.7576 2.2769 0.4838
2007-2008 141.79 2.9363 1.7792 0.3781
2008-2009 139 2.8785 1.7442 0.3706

4
 2009-2010 191.8 3.9719 2.4068 0.5114 2005-2006 5.857 2.479
2010-2011 190.26 3.94 2.3874 0.5073 2006-2007 0.025 2.276
2011-2012 256.31 5.3078 3.2162 0.6834
2007-2008 0.124 1.779
2012-2013 206.47 4.2757 2.5908 0.5505
2008-2009 0.901 1.744
2013-2014 149.4 3.0939 1.8747 0.3983
2009-2010 0.288 2.406
2014-2015 179.4 3.7151 2.2512 0.4783
2010-2011 3.465 2.387
Con la misma lógica de la extrapolación de áreas y caudales se 2011-2012 8.465 3.216
calculó los caudales generados en las cuencas directas de los
2012-2013 0.041 2.59
glaciares analizados.
La validación de los resultados se ha realizado comparativamente 2013-2014 0.915 1.874
con los resultados obtenidos en el glaciar Inkachiriasq’a en el 2014-2015 0.028 2.251
nevado Salkantay, realizados por la Autoridad nacional del agua
(ANA) ubicado a 16 Km del nevado Oq’obamba y a 22 km del 3.6 Resultados del análisis de correlación de las variables independiente y
nevado Moyoq.[CITATION Mon14 \l 10250 ] dependiente.

El análisis de correlación del coeficiente de Pearson mediante el

Estos resultados fortalecen al estado de arte propuesto
inicialmente, ya que el análisis de los glaciares de menor altitud ha software statgraphics[CITATION Bat \n \l 10250 ] entre las
permito obtener los resultados mostrados debido a que: Son de variables arrojó un valor r= 0.467, lo que indica que no existe una
fácil acceso, economizan los costos de investigación, evitan las relación significativa entre el retroceso de los glaciares y la
costosas campañas en equipos y personal especializado, los disponibilidad hídrica [CITATION Lah03 \l 10250 ]sin embargo,
resultados son visibles y comparables con medios disponibles en el se puede analizar que dicho valor es bajo debido a la Sublimación de
medio. los glaciares, los caudales de los ríos están relacionado
estrechamente con el aporte hídrico de los deshielos de los glaciares.
Tabla 4. Cálculo de caudales de los riachuelos Llamawasi, Yawarmaki, Para fines de comprobación del resultado se utilizó el SPSS
Q’elwaq’ocha y Ch’akil. (Statistical Package for the Social Sciences ) con el que obtuvo igual
Periodo Caudales en las Cuencas evaluadas (m3/seg) valor r=0.467.
Llamawas
(Años)
i
Yawarmaki Q’elwaq’ocha Ch'akil 3.7 Discusión de resultados
El valor de r=0.467 indica que no existe una relación significativa en
2000-2001 0.1246 0.0698 0.1389 0.0434
entre el retroceso glaciar y la disponibilidad hídrica. Este hecho
2001-2002 0.092 0.0515 0.1025 0.032 confirma las evidencias observadas en el campo, la falta torrente
2002-2003 0.1246 0.0698 0.1389 0.0434 directo a las cuencas receptoras y se plantea la hipótesis de la
sublimación glaciar presente en esta parte de la cordillera de
2003-2004 0.0938 0.0526 0.1046 0.0327
Vilcabamba. La sublimación es la transformación directa del hielo
2004-2005 0.0696 0.039 0.0776 0.0243 del estado sólido al gaseoso, proceso en la requiere una cantidad de
2005-2006 0.1212 0.0679 0.1352 0.0422 energía equivalente a 2834 Joules/kg, más de 8.50 veces más que la
fusión, entonces indica que las condiciones de temperatura
2006-2007 0.1113 0.0538 0.1241 0.0388
ambiental se han incrementado considerablemente en los nevados
2007-2008 0.087 0.0487 0.097 0.0303 analizados, la perdida de la superficie glaciar ha implicado que no se
2008-2009 0.0853 0.0478 0.0951 0.0297 está reflejando la luz solar y este hecho hace que se acumule la
2009-2010 0.1177 0.0659 0.1312 0.041
energía del sol en las rocas generando un incremento de calor a
nivel local y esta energía es la causa de la sublimación asociado a
2010-2011 0.1167 0.0654 0.1301 0.0407 otras causas materia de investigación. Este fenómeno no es
2011-2012 0.1572 0.0881 0.1753 0.0548 exclusivo en la cordillera de Vilcabamba, la sublimación de
2012-2013 0.1267 0.071 0.1412 0.0441 glaciares has sido reportado en el glaciar Juncal (Chile) donde se ha
reportado que hasta el 50 % de los glaciares se están evaporando
2013-2014 0.0916 0.0514 0.1022 0.0319
[ CITATION Mel \l 10250 ], indican que el concepto preliminar
2014-2015 0.11 0.0617 0.1227 0.0383 de la sublimación glaciar debido a los efectos de la baja humedad, en
ambientes montañosos áridos, ha quedado desfasado ya que las
Tabla 5. Resultados finales de las variables independientes y últimas investigaciones realizadas indican que la elevación del
dependientes glaciar juega un factor preponderante, debido a que las condiciones
frías, secas y de exposición al viento; los glaciares se incrementan
Años Área(km2) Caudal(m3/seg)
con la altitud (msnm) y están asociadas a elevadas tasas de
2000-2001 1.398 2.548 sublimación.
2001-2002 5.148 1.881
2002-2003 0.39 2.548 4. CONCLUSIONES
2003-2004 2.609 1.919 La disminución progresiva y significativa de los caudales de los ríos
2004-2005 0.088 1.424 Lechería y Colorado se debe al retroceso glaciar que han sufrido los

5
glaciares:Llamawasi,Yawarmaki(NevadoOq’obamba),Q’elwaq’ocha
, Ch’akilq’asa (Nevado Moyoq) cuya área 24.75 Km2 corresponde
al 98.21 % de su superficie glaciar.

La correlación entre variables de áreas y los caudales cuyo resultado

es r=0.467, indica que no existe relación significativa entre el
retroceso glaciar y la disponibilidad hídrica. La inexistencia de una
correlación directa entre el retroceso glaciar y la disponibilidad
hídrica se debe a la sublimación glaciar, hecho que se ha podido
verificar en los trabajos de campo.

En el valle de Limatambo existe una política agresiva de

plantaciones de paltas para fines de exportación, cuya demanda
hídrica se incrementa sin embargo la oferta hídrica ha disminuido;
Por lo que urgente tomar medidas de manejo responsable del uso del
agua.
Agradecimientos

Expresamos nuestra gratitud a las siguientes instituciones y

personas: Escuela pos grado Universal Nacional San Antonio Abad
del Cusco, Universidad Nacional Micaela Bastidas de Apurímac,
Servicio Nacional de Meteorología e Hidrología, Servicio Nacional
de Áreas Naturales Protegidas por el Estado, Autoridad Nacional del
Agua Zona XI. Pampas Apurímac, Comunidad Campesina de
Sondorf, Ing. Juan Jiménez, Biól. Américo Cabrera, Franco
Ubalde,Rosmel Catalán, Bernabé Valenzuela, Nelson
Onzueta,Egoavil Gallardo, Edu Arias, Vladimir . Escobedo,
Elizabeth Tevés y William Ccanahuire. y Demetrio Gallegos.
Referencias bibliográficas:

Ayala, A., Pellicciotti, F., Peleg, N., & Burlando, P. (2017). Melt and surface sublimation across a glacier in a dry environment:
distributed energy-balance modelling of Juncal Norte Glacier, Chile. Cambridge University UK: Jornal of Glaciology.

Baraer, M., Mark, B. G., Mckenzie, J. M., Condom, T., Huh, K.-I., Portocarrero , C., . . . Gómez , J. (2012). Glacier recession and water
resources in Peru’s Cordillera Blanca. Lima-Perú.: Journal of Glaciology, Vol. 58, No. 207. doi:doi: 10.3189/2012JoG11J186

Batanero B., C., & Diaz B., C. (2008). Análisis de datos con Statgráphics. Granada España: Universidad de Granada.

Campos, N. O. (2012). Glacier evolutiónin the South West slope of Nevado Coropuna (Cordillera de Ampato, Perú. Madrid-España:
Universdad Complutense de Madrid.

Colindres Estrada, W. (2016). Curvas de calibración de caudales de las estaciones hidrológicas de la cuenca Coatán a la cuenca madre
vieja, correspondiente a la vertiente del pacífico para el periodo 2010-2013. Tegucigalpa-Honduras: Universidad de San Carlos.
Obtenido de http://biblioteca.ingenieria.usac.edu.gt/

Dozier, J. (1989). Spectral Signature of Alpine Snow Cover from Lansat Thematic mapper”.Remote sensing of environment. Remote
sensing of enviroment ,volumen 28 Elsevier. doi:doi.org/10.1016/0034-4257(89)90101-6

Guardamino, L., & Drenkhan, F. (2016). ).” Evolución y potencial amenaza de lagunas glaciares en la cordillera de Vilcabamba (Cusco y
Apurímac, Perú) entre 1991 y 2014. Zurich-Suiza: University of Zurich.

Hall k, D., George A, R., James L, F., & Sujay V, K. (2010). Development and evaluation of a cloud-gap-filled MODIS daily snow-cover
product. Remote Sensing of Environment 114-Elsivier.

Kozhikkodan Veettil, B., Pereira S, F., Wang, S., Valente, P., Grondona , A., Rondon, A., . . . Simooes, J. (2016). Un análisis comparativo
del retroceso glaciar en los andes tropicales usando Teledetección. SantiagoChile: Dialnet Nº 51, 2016, págs. 3-36. Obtenido de
SSN 0716-0364,

Lahura, E. (2003). El coefeciente de correlacón y las correlaciones espúreas. Lima-Perú: PUCP. Obtenido de
https://core.ac.uk/download/pdf/6445817.pdf

Licciardi M, J., Schaefer M, J., Taggart R, J., & Lund C, D. (2009). Holocene Glacier Fluctuations in the Peruvian Andes Indicate
Northern Climate. Sciencie Vol. 325 N° 5948 págs.1677-1679. doi:10.1126 / science.1175010

6
Marangunic D., C. (2008). Manual de Nieve y Nivometría Volumen 1. Santiago Chile: Gobierno de Chile.

Martin de la Calle, A. (2009). Evolución reciente de los Glaciares de la vertiente sur del Nevado Salcantay (Perú). Madrid-España:
Universidad Complutense de Madrid.

Pellikka, P., & Gareth W, R. (2010). Remote sensing of glaciers Techinques for topographic spatial and thermatic mappiing of glaciae.
Londres UK.

Rabatel, A., Francou, B., Soruco, A., Gémez, J., Cáceres, B., Cevallos L, J., . . . Vuille, M. (2013). Current state of glaciaers in the tropical
Andes:a multi-century perspective on glaciaer evolution and climate change. The Cryosphere, 7, 81–102, 2013.
doi:https://doi.org/10.5194/tc-7-81-2013

Rau, F., Mauz, F., Vogt, S., Siri Jodha, S., & Raup, B. (2005). Illustrated GLIMS Glacier Classification Manual -Glacier Classification
Guidance for the GLIMS Glacier Inventory. Freiburg (Germany): Institut für Physische Geographie.

Suarez, W., Macedo, N., Montoya, N., Arias, S., Simin, Schauwecker, S., . . . Condom, T. (2015). Balance energético neto (2012-2014) y
evolución temporal del nevado Quisoquipina en la región Cusco (1990-2010). Lima-Perú: Senamhi (Servicio Nacional de
meterología e hidrología del Perú).

Ubeda Palemque, J. (2011). El impacto del cambio climático en los glaciares del complejo volcánico de Coropuna (Cordillera occidental
de los andes centrales). Madrid España: Universidad Complutense de Madrid.

Unidad de Glaciología y Recursos Hídricos -Servicio nacional de áreas naturales protegidas. (2014). Monitoreo Glaciar Inkachiriasq´a.
Cusco” Santuario Histórico de Machupicchu. Lima -Perú: Servicio nacional de áreas naturales protegidas (SERNANP).

Viulle, M., Bradley, R., Francou, B., Kaser, G., & Mark, B. (2012). Climate Change in the tropical Andes - Impacts and consequences for
glaciation and water resources Part I: The scientific basis”. Massachusetts-USA: PRAA.

Vuille, M. (2013). Climate Change and Water Resources in the Tropical Andes. IDB Interamerican Development Bank.