UNIVERSIDAD NACIONAL FEDERICO VILLARREAL INSTITUTO DE INVESTIGACIÓN
FACULTAD DE INGENIERÍA GEOGRÁFICA, AMBIENTAL Y ECOTURISMO EVALUACIÓN DE LA COBERTURA GLACIAR UTILIZANDO TÉCNICAS DE TELEDETECCIÓN ESPACIAL EN LA SUBCUENCA DEL RIO QUELLCAYHUANCA RIO SANTA
–
Dr. Walter Gómez Lo ra Dr. No éZam o ra Talaver an o L ic . Cé s ar Cab r el L a Ro s a
CAMBIO CLIMATICO
En los ultimos 30 años, el Peru ha perdido 1/3 de la superficie de sus glaciares y se ha proyectado en los glaciares situados a 5500 msnm desapareceran. En los ultimos 35 años, se ha producido en el Peru una reduccion del 12% del agua dulce en la zona costera, donde vive el 60% de la poblacion del pais. Se estima una perdida de alrededor de 700 m3 de agua. La costa peruana, donde se desarrolla gran parte de la agricultura, se vera seriamente afectada por su alta dependencia del agua dulce proveniente de los glaciares andinos. Ciudades como Lima estarán bajo seria amenaza debido a la desaparición de los glaciares que ocasionarían un estrés
UBICACIÓN DE LAS 18 CORDILLERAS NEVADAS DEL PAÍS LEYENDA Cordilleras nevadas Ramales de la Cordillera Andina Andina DISTRIBUCION DE LAS CORDILLERAS NEVADAS 1.- Blanca * 2.- Huallanca 3.- Huayhuash 4.- Raura 5.- Huagoruncho 6.- La Viuda 7.- Central 8.- Huaytapallana 9.- Chonta
10.- Ampato 11.- Vilcabamba 12.- Urubamba 13.- Huanzo 14.- Chila 15.- La Raya 16.- Vilcanota 17.- Carabaya 18.- Apolobamba
Cordillera Blanca
* La Cordillera Blanca es la cordillera glaciar tropical más alta y extensa del PERU mundo.
Área Territorial = 1`285,216 1`285,216 Km2 Población = 26`392,000 Hab.
MONITOREO DE GLACIARES EN EL PERÚ Unidad de Glaciología y Recursos Hídricos ARTESONRAJU (INRENA-IRD-INNSBRUCK)
SHALLAP (INRENA - INNSBRUCK)
CORDILLERA BLANCA
URUASHRAJU (INRENA)
YANAMAREY YANAMAREY (INRENA - IRD)
PASTORURI (INRENA)
CORDILLERA CENTRAL
SHULLCON
CORDILLERA VILCABAMBA
SALKANTAY
(INRENA-IRD)
(INRENA – 2007)
• La
capacidad actual de nuestros glaciares es de 43 mil millones de m3 (2007)
• Y
representa el 60% de nuestro ‘stock’ de agua • El 40% de nuestro ‘stock’ de agua -25 mil millones de m3- participa
en el flujo hídrico •
Se proyecta que en 10 años el 40% de nuestro ‘stock’ de agua (glaciares)
se irá en el flujo hídrico
Escenario Altamente Posible Fuente: MINAM
RETROCESO GLA CIAR YANA MAREY (CORDILL ERA BL ANCA – ALTITUD 4786 m sn m )
1987
1981
2006
1997
1981 1987 1997 2006
FORMULACIÓN DEL PROBLEMA La problemática del retroceso del glaciar se inició con el estudio de la tasa de retroceso del glaciar, utilizando las herramientas de la Teledetección el cual permitirá realizar la clasificación de los glaciares por el método del NDSI (Índice Normalizado y Diferenciado de Nieve) y conocer el área de la cobertura glaciar, para posteriormente analizar hidrológicamente los aportes por deshielo de la masa del área glaciar de la Subcuenca Quellcayhuanca. Siendo el problema principal la pérdida de área glaciar del nevado de la subcuenca rio Quellcayhuanca que se constituye hidrológicamente en un aporte de escorrentía a las lagunas y quebradas cuyo rendimiento hídrico debe evaluarse. ¿En qué medida la pérdida del área glaciar influye en el rendimiento hídrico en la subcuenca del rio Quellcayhuanca?
FORMULACIÓN DEL PROBLEMA Siendo los problemas segundarios:
¿Cómo se caracteriza el comportamiento hidrológico y pluviométrico medido en las estaciones hidrometeorológicas?, ¿Cuál es la perdida de área glaciar y la tasa de retroceso de la subcuenca del rio Quellcayhuanca? Y ¿Cuál es el aporte hídrico por perdida de área glaciar en la subcuenca del rio Quellcayhuanca?
Se plantea la hipótesis: La disminución del área glaciar influirá en el rendimiento hídrico de la subcuenca del rio Quellcayhuanca.
OBJETIVO DE LA INVESTIGACIÓN OBJETIVO GENERAL
conocer las pérdidas del área glaciar mediante técnicas hidrológicas y de teledetección con la finalidad de determinar su la influencia en el rendimiento hídrico de la subcuenca del rio Quellcayhuanca utilizando técnicas de teledetección espacial OBJETIVOS ESPECÍFICOS
1) Evaluar los años pluviométricos e hidrológicos característicos medio, húmedo y seco de la subcuenca del rio Quellcayhuanca. 2) Evaluar las pérdidas de área glaciar mediante las técnicas de SIG y Teledetección con la finalidad de determinar la evolución y su tasa de retroceso glaciar 3) Determinar el aporte hídrico por perdida de área glaciar mediante técnicas hidrológicas.
JUSTIFICACIÓN Se justifica en este contexto, el presente estudio que pretende establecer si el retroceso glaciar ha influido en las últimas décadas en los caudales, mediante el comportamiento de los cauces en la subcuenca del rio Quellcayhuanca, ubicado en la cuenca del rio Santa, departamento de Ancash. Utilizando las técnicas del Sistema de información Geográfica y Teledetección del glaciar y métodos hidrológicos. El presente documento aportara al entendimiento del comportamiento hidrológico del nevado, y evaluar la disponibilidad hídrica y de reservas en la subcuenca del rio Quellcayhuanca como fuente de abastecimiento al consumo humano y agrícola. Esta ultima la principal actividad generadora de trabajo en la población rural.
IMPORTANCIA El estudio es importante porque los glaciares son formaciones frágiles de gran importancia como reguladores de la disponibilidad de agua dulce. Hoy grandes sectores de la población mundial no tienen acceso al agua dulce, ya sea porque falta en su región o porque está mal distribuida. Esa escasez se agrava por la contaminación del agua dulce con agroquímicos, efluentes industriales y residuos que han aumentado de forma alarmante, mientras el mundo requiere más agua dulce pura y confiable, y el calentamiento global está modificando el ciclo hidrológico.
UBICACIÓN Y LOCALIZACION
Cuenca SANTA
Ubicación Política Departamento Provincia ANCASH Huaraz Dist. Huaraz
Cuenca Cañete
Latitud Mín. Max 09º15’36”
09º23’20”
Longitud Min
Max
77º19’4.8”
76º36’36”
MATERIALES Y EQUIPOS MATERIAL Cartas
Escala 1:100000, Hojas:
Mapas
19h, 19i, 20h, 20i .
Imágenes LANDSAT ETM ortorectificadas
Temáticos
Mapa Topográfico Mapa Ecológico Mapa de Cobertura Vegetal Mapa Geológico Mapa de Temperatura Mapa de Isoyetas
EQUIPOS
• Imágenes
Nacional
Computadora Core Duo Scaner Hp A3
SOFTWARE
ArcGis v. 10 Autocad 2011 Envi 4.8 Erdas Imagine 2010 Excel - Estadística gráfica
Año
Tipo de Imagen
Path / Row
Fecha
Tipo
1987
Imagen satelital
008/067
31/05/1987
GeoTIFF
1989
Imagen satelital
008/067
30/12/1989
GeoTIFF
2007
Imagen satelital
008/067
22/05/2007
GeoTIFF
2010
Imagen satelital
008/067
14/05/2010
GeoTIFF USGS: US Geological Survey
Las imágenes satelitales LANDSAT se adquirieron orto rectificadas y corregidas geométricamente por el proveedor U.S Geological Survey (USGS) las cuales son las siguientes:
PROGRAMA ArcGis 10.0 (integrado Arc2Earth)
Generación e Integración de las bases de datos, junto a los cálculos asociados a las características de los glaciares.
Envi 4.8
Procesamiento de imágenes satelitales para la generación de polígonos para los glaciares descubiertos
Google Earth Pro
Corroboración de aspectos morfológicos de glaciares rocosos y descubiertos
(NDSI)
NDSI=(TM2-TM5)/(TM2+TM5) Este índice es apropiado para diferenciar nieve
DESCRIPCIÓN
TIPO DE INVESTIGACION El tipo de investigación es Transversal consiste en describir variables y analizar su incidencia e interrelación en un momento dado.
NIVELES DE INVESTIGACIÓN Descriptiva. Se describen las variables en estudio para identificar su relación que existe entre ellas.
Correlacional. Se analiza la relación entre los diferentes variables que intervienen en la en la influencia glaciar en el rendimiento hídrico de la subcuenca Mala.
MÉTODO •
Determinación De La Evolución De La Cobertura Glaciar
Método del NDSI (Índice de Nieve de Diferencia Normalizada).- Este método es utilizado para el cartografiado de superficies glaciares (hielo); y su utilidad consiste en la eliminación de cuerpos de agua, nubes y nieve que podrían cartografiarse cuando se utilizan imágenes satelitales, el álgebra de bandas es el siguiente: …………… Ecuaci ón
Nº 01
METODOLOGIA Análisis Multitemporal del Glaciar Selección y Adquisición de Imágenes Selección de las imágenes Satélite correspondiente al área de estudio.
Adecuación de Cartografía
Delimitación del Área de estudio
Selección de Cartografía Adquisición de Cartografía Actualización de Cartografía
Reconocimiento general de área de estudio Selección del área de estudio Ubicación de coordenadas del área de estudio
Reprocesamiento y correcciones de imágenes
Clasificación no supervisada
Cortes específicos del área de estudio
Modelo de Elevación • Modelo Elev. Digital DEM •Modelo TIN
•Variación
de altura de la capa de hielo nivel del glaciar
•Curvas de
Método NDSI
Modelamiento Espacial Obtención de cartografía temática combinada. Definición de criterios de ponderación por cada temática
INTEGRACIÓN CON CARTOGRAFÍA TEMÁTICA
CARACTERIZACIÓN DE ÁREAS GLACIARES
RESULTADOS DE ESTUDIO ÄREA GLACIAR
VOLUMEN GLACIAR
IMÁGENES LANDSAT 7 (ETM) BANDAS
AMPLITUD DE ONDA
RESOLUCIÓN ESPACIAL
1= Azul
0.45 - 0.52 µm
30 m
2= Verde
0.52 - 0.60 µm
30 m
3= Rojo
0.63 - 0.90 µm
30 m
4= Infrarrojo cercano
0.76 – 0.90 µm
30 m
5= Infrarrojo medio
1.55 – 1.75 µm
30 m
6= Infrarrojo termal
10.40 – 12.50 µm
60 m
7= Infrarrojo medio
2.08 – 2.35 µm
30 m
8= Pancromático
0.52 – 0.90 µm
15 m
IMÁGENES LANDSAT 7 (ETM) BANDAS
APLICACIONES
1= Azul
Aguas costeras , diferenciación entre suelo y vegetación, vegetación conífera y caducifolia.
2= Verde
Mapeo de vegetación, calidad de agua .
3= Rojo
Diferenciación de especies vegetales y áreas urbanas, uso del suelo, agricultura, calidad de agua .
4= Infrarrojo cercano
Cuerpos de agua, mapeo geomorfológico mapeo geológico áreas de incendios, áreas húmedas (agricultura vegetación ).
5= Infrarrojo medio
Uso del suelo, humedad de la vegetación, diferenciación entre nubes y nieve.
6= Infrarrojo termal
Estrés térmico en plantas, corrientes marinas, propiedades termales del suelo.
7= Infrarrojo medio
Identificar minerales, mapeo hidrotermal.
8= Pancromático
Catastro rural, infraestructura, uso de suelo.
Indice NDSI Normalized Difference Snow Index
R: 7 G:4 B:2
2 1 3
31/05/1987
N O
4
E S
La visualización en combinación R7G4B2 para el año de 1987 (izq.) y el NDSI obtenido de la imagen Landsat TM de 1987 (derecha)
Igual se realizó para 1989, 2007 y 2010 31/05/1987
MÉTODO Cálculo del aporte hídrico del nevado en años hidrológicos secos.Método análisis de la curva de frecuencia.- La curva de duración es un procedimiento gráfico para el análisis de la frecuencia de los datos de caudales y representa la frecuencia acumulada de ocurrencia de un caudal determinado. Esta curva puede ser definida para caudales diarios, mensuales, anuales, etc.
Método de análisis adimensional.- El método para hallar esta curva regional, es comparar gráficamente las diferentes curvas de frecuencia, existentes en la zona, adimensionalizadas por el caudal promedio diario correspondiente. La adimensionalización se hace mediante la siguiente expresión: ………………………………………………… .. Ecuaciónón Nº 02
Donde:
Z: Caudal adimensional Q: Caudal registrado Qmedio: Caudal promedio diario multianual
RESULTADOS
ESTACIONES PLUVIOMÉTRICAS Se recopilo la información de las estaciones pluviométricas cercanas a la Subcuenca Honda, el cual nos permitirá generar isoyetas de precipitación para diferentes años característicos, el cual se detalla a continuación (según el Ministerio de Energía y Minas. 1981): Tabla Nº 03.- Estaciones Pluviométricas ESTACIÓN
ALTITUD
PRECIPITACIÓN
LATITUD
LONGITUD
TIPO DE INFORMACIÓN
Santiago de Chuco
3129
580
8°08'
78°10'
Pp. total media mes
Chancos
2840
520
9º 04'
77º39'
Pp. total media mes
Cahuich
4550
950
9°41'
77°17'
Pp. total media mes
Recuay
3394
760
9°43'
77°27'
Pp. total media mes
Estación Huaraz
Altitud 3052
Precipitación 620
Latitud 9°31'
Longitud 77°32'
Tipo de Información Pp. total mensual*
Fuente: Ministerio Energía y Minas, Evaluación Potencial Hidroeléctrica Nacional, Volumen 5, (1980)
CARTOGRAFÍA Tabla Nº 2: Archivos digitales en formato Shape de la
Quellcayhuanca
ARCHIVO Rios.shp Curvas.shp Cotas.shp Lagos.shp Islas.shp Limite.shp Polurbano.shp Txt_rios.shp Txt_cerros,shp Txt_lagos.shp Malla100
subcuenca del rio
DESCRIPCION Hidrografía Curvas de nivel Cotas de terreno Lagunas Islas Divisoria de aguas de la cuenca Casco Urbana Nombre de ríos Nombre de cerros Nombres de lagos Índice hojas al escala 1: 100000
Cartas digitales a escala 1:100000, Hojas 19h, 19i, 20h, 20i las cuales se homogenizaron para el diseño de la subcuenca, siendo su sistema de referencia UTM, datum WGS 84 zona18.
Parámetros fisiográficos Tabla N°02. Parámetros fisiográficos de la subcuenca del rio Quellcayhuanca Parámetros Área Perímetro Longitud del Río Principal Altura del punto más alto del Río Altura de salida de la Cuenca Longitud total de Ríos Pendiente Media del Río Principal Altura Media de la Cuenca Densidad de Drenaje Índice de Compacidad Lado Mayor del Rectángulo Equivalente Lado Menor del Rectángulo Equivalente Índice de Pendiente Factor de Forma
Magnitud
Unidad
248.5
km²
81.28
Km
28.45 6250
Km msnm
3050
msnm
140.49 5.25
Km %
4487.26
M
0.57 1.44
km/km² Adimensional
33.14
Km
7.5
Km
28.01 0.31
% Adimensional
Pluviometría Con las estaciones pluviométricas que están en el ámbito del área de estudio y en la que solo Chancos pertenece a la zona baja de la subcuenca. Se determinó la relación precipitación vs altura, tal como se muestra en la ecuación regional (N °04) que representa el comportamiento pluviométrico de la subcuenca y que es muy similar a la subcuenca del rio Quellcayhuanca, determinándose la ecuación regional del comportamiento pluviométrico donde:
P = 0.000128H2 + 1.205311H……………………………….… Ecuación N ° 02 R = 0.98 Donde: Pp = precipitaci ón H = Altura
Análisis Regional En el análisis regional de la precipitación se tiene que para la cuenca de estudio la precipitación total anual es de 877.2 valor muy similar al obtenido por Loroña (2011) de 899mm., con un error de -1.3% lo que explica que su ubicación le corresponde características similares en el comportamiento pluviométrico regional y por lo tanto su distribución mensual generada debe ser muy parecida; dado que la media, desviación estándar, valores aleatorios son similares. Por lo tanto la información pluviométrica obtenida por Loroña (2011) son valores similares a los encontrados en este estudio como representativa cuyos valores se presenta en la tabla N°04
Tabla N° 04. Precipitaciones media mensual generada para la Subcuenca de Quillcay (1965-66/2009-10) Año Hidrológico 1965-1966 1966-1967 1967-1968 1968-1969 1969-1970 1970-1971 1971-1972 1972-1973 1973-1974 1974-1975 1975-1976 1976-1977 1977-1978 1978-1979 1979-1980 1980-1981 1981-1982 1982-1983 1983-1984 1984-1985 1985-1986 1986-1987 1987-1988 1988-1989 1989-1990 1990-1991 1991-1992 1992-1993 1993-1994 1994-1995 1995-1996 1996-1997 1997-1998 1998-1999 1999-2000 2000-2001 2001-2002 2002-2003 2003-2004 2004-2005 2005-2006 2006-2007 2007-2008 2008-2009 2009-2010
Sep
Oct.
Nov.
Dic.
Ene
Febr.
Mar
Abr.
May
Jun.
Jul.
Ago.
PP. total
69,2 18,9 38,7 40,9 3,7 101,2 13,9 16,1 70,7 31,5 43,9 5,6 39,8 12,8 43,1 58,9 32,8 21,7 30,7 57,9 17,5 47,0 65,4 51,9 43,2 13,3 0,0 55,5 68,9 41,7 20,8 17,2 26,5 10,7 66,3 21,6 15,8 16,5 7,9 51,5 4,0 15,1 8,4 24,9 2,1
113,8 199,1 175,2 93,0 86,6 80,0 85,2 41,6 128,9 29,8 61,8 10,9 26,7 45,6 34,7 26,8 166,5 99,4 65,6 26,1 14,8 56,6 63,5 94,2 155,6 186,1 69,4 86,9 220,9 58,7 84,9 77,2 44,0 142,2 20,8 83,0 51,6 105,8 31,0 108,9 47,2 52,1 107,5 140,0 158,3
54,5 95,2 44,8 36,3 141,3 79,3 59,4 76,6 150,7 37,7 67,2 20,5 134,6 81,7 69,2 58,1 126,0 42,3 63,2 21,1 84,0 84,0 126,7 108,7 45,8 165,2 85,8 50,2 186,5 63,7 90,1 45,9 83,6 74,9 124,9 42,2 111,0 144,1 33,2 95,6 21,1 58,9 60,8 46,2 144,5
161,7 67,0 72,6 40,6 143,9 142,7 154,5 155,7 101,8 66,1 56,2 62,7 111,0 66,7 57,3 68,9 125,5 131,4 91,5 98,9 194,0 144,9 206,9 79,6 10,6 86,2 37,1 52,7 264,2 93,1 140,9 49,9 213,7 83,2 145,6 228,0 52,2 130,8 178,4 110,1 117,0 167,2 65,3 99,8 196,6
154,5 140,2 111,0 78,7 209,0 116,7 218,6 141,6 145,1 156,4 179,2 124,3 68,2 43,3 75,8 105,9 180,2 103,4 109,5 259,3 226,5 197,9 174,5 200,3 124,3 80,9 2,6 170,0 196,9 52,4 152,6 85,5 187,9 120,1 62,7 123,0 122,0 112,7 64,7 158,9 112,6 101,2 123,0 187,7 107,6
93,4 180,4 61,7 111,1 95,3 142,1 154,2 120,9 155,9 179,4 84,4 94,6 167,7 137,8 114,8 30,3 341,0 130,0 134,7 91,9 153,7 141,2 150,7 158,4 48,2 91,6 80,2 149,0 210,2 111,1 181,8 187,5 213,8 290,5 165,0 138,8 179,2 59,5 116,7 93,0 106,0 91,5 155,8 148,9 51,9
88,0 201,8 131,6 169,5 105,6 236,7 280,9 154,6 209,6 250,1 157,3 149,6 64,6 124,3 311,5 169,9 255,7 123,8 186,9 192,8 117,9 85,9 100,6 106,5 85,0 202,3 127,2 278,9 369,0 162,9 235,7 63,7 186,2 146,0 192,2 170,4 178,8 144,6 79,8 220,9 249,1 178,9 113,6 237,2 130,9
78,4 66,3 20,7 95,8 173,1 105,9 153,7 145,6 155,7 86,3 68,8 31,7 75,7 58,1 111,5 40,9 25,9 106,6 101,0 185,4 126,5 84,2 129,5 92,9 55,7 139,6 56,2 175,7 164,1 145,5 70,1 32,9 122,2 104,7 116,2 153,5 109,4 96,4 51,4 36,4 153,6 196,5 90,5 124,5 75,6
26,0 50,9 14,2 11,6 72,3 3,3 27,3 40,8 0,3 86,1 45,5 22,1 22,2 72,7 37,9 8,4 22,0 34,4 33,7 47,8 18,5 7,9 37,4 37,5 37,3 59,1 33,8 32,2 58,3 44,0 31,8 14,5 0,0 63,3 89,3 31,5 17,5 30,1 24,7 3,0 5,0 33,4 2,3 31,8 48,6
0,3 3,3 0,0 1,5 4,4 3,0 0,7 9,1 20,4 7,6 11,1 0,0 0,0 5,6 0,0 5,4 40,7 3,7 10,1 2,2 5,7 0,0 2,8 8,0 3,7 0,0 20,5 0,0 5,5 0,7 0,0 0,0 12,7 7,0 2,0 5,3 0,6 3,8 4,4 0,0 18,1 0,7 7,8 9,6 6,4
0,0 14,6 0,3 2,3 0,2 0,6 1,0 9,6 5,0 0,0 1,2 0,6 2,9 6,2 7,4 0,0 0,7 9,4 10,1 14,1 28,0 0,0 0,0 0,0 0,0 6,6 0,0 6,1 0,9 0,0 0,0 0,0 0,0 3,9 3,7 4,2 0,0 1,5 16,6 0,0 1,7 7,3 10,3 0,3 0,0
0,9 9,4 30,7 0,0 7,2 38,3 33,7 15,6 24,1 19,9 0,8 1,4 0,0 0,0 6,2 17,9 11,7 11,9 23,4 21,9 14,3 7,0 24,5 25,9 5,4 0,0 14,8 0,0 7,0 5,1 7,7 0,0 0,0 0,0 9,0 3,2 0,0 1,7 0,0 9,1 5,5 0,9 3,3 2,3 3,6
840,7 1047,2 701,5 681,4 1042,5 1049,8 1182,9 927,7 1168,2 951,0 777,5 524,0 713,4 654,7 869,5 591,4 1328,7 817,9 860,4 1019,5 1001,5 856,6 1082,5 963,8 614,9 1031,0 527,7 1057,3 1752,5 778,9 1016,4 574,3 1090,6 1046,4 997,8 1004,6 838,0 847,5 608,9 887,6 841,1 903,7 748,6 1053,2 926,1
ECUACIÓN REGIONAL DE DÉFICIT Habiendo obtenido la precipitación anual y teniendo la ecuación de déficit: D = 1.195P – 618.6 ………………………………………………… ..Ecuación N° 04 Se determina el caudal superficial por lluvia (Ver Tabla N° 5), ejemplo: Precipitación: 1965-1966 = 840.66 mm Déficit = 386 mm Es = P – D …………………………………………………………… Ecuación N° 05 Es = 840.7 – 386 = 454.7 Q = Es x A / 31567
…………………………………………………
Q = (454.7 * 248.5) / 31567 = 3.58 m 3/seg.
Ecuación N° 06
Tabla N° 05. Determinación de Déficit, Escorrentía y caudal de nevado y escorrentía por precipitación. Año 1965-1966 1966-1967 1967-1968 1968-1969 1969-1970 1970-1971 1971-1972 1972-1973 1973-1974 1974-1975 1975-1976 1976-1977 1977-1978 1978-1979 1979-1980 1980-1981 1981-1982 1982-1983 1983-1984 1984-1985 1985-1986 1986-1987 1987-1988 1988-1989 1989-1990 1990-1991 1991-1992 1992-1993 1993-1994 1994-1995 1995-1996 1996-1997 1997-1998 1998-1999 1999-2000 2000-2001 2001-2002 2002-2003 2003-2004 2004-2005 2005-2006 2006-2007 2007-2008 2008-2009 2009-2010
Pp. Total (mm) 840,66 1047,17 701,50 681,35 1042,55 1049,81 1182,94 927,69 1168,19 951,01 777,51 524,03 713,37 654,74 869,45 591,42 1328,68 817,89 860,42 1019,50 1001,47 856,60 1082,50 963,80 614,90 1030,96 527,70 1057,30 1752,48 778,90 1016,40 574,30 1090,60 1046,45 997,76 1004,59 838,02 847,48 608,88 887,55 841,07 903,75 748,63 1053,18 926,14
Escorrentía (mm) 454,67 414,40 481,81 485,74 415,30 413,89 387,93 437,70 390,80 433,15 466,99 516,42 479,49 490,93 449,06 503,27 359,51 459,11 450,82 419,80 423,31 451,56 407,51 430,66 498,69 417,56 515,70 412,43 276,87 466,71 420,40 506,61 405,93 414,54 424,04 422,70 455,19 453,34 499,87 445,53 454,59 442,37 472,62 413,23 438,00
Q escor (m3/seg) 3,59 3,27 3,80 3,83 3,28 3,27 3,06 3,45 3,08 3,42 3,68 4,07 3,78 3,87 3,54 3,97 2,84 3,62 3,56 3,31 3,34 3,56 3,22 3,40 3,93 3,29 4,07 3,25 2,18 3,68 3,32 4,00 3,20 3,27 3,35 3,34 3,59 3,58 3,94 3,52 3,59 3,49 3,73 3,26 3,46
Q real (m3/seg) 6,72 7,36 6,28 6,22 7,34 8,97 7,47 6,73 8,19 7,08 6,50 6,92 6,47 7,64 6,14 8,08 7,45 8,65 7,53 5,11 5,90 7,49 7,64 7,43 6,77 7,03 7,04 8,28 7,93 7,24 7,14 4,75 6,19 7,36 7,21 7,23 6,71 6,74 6,00 6,86 6,72 6,91 6,43 7,38 6,98
Q nev (m3/seg) 3,13 4,09 2,48 2,39 4,07 5,70 4,41 3,27 5,10 3,66 2,82 2,85 2,68 3,76 2,59 4,10 4,61 5,03 3,97 1,80 2,56 3,92 4,42 4,03 2,84 3,74 2,98 5,02 5,75 3,56 3,82 0,75 2,99 4,09 3,86 3,89 3,12 3,16 2,05 3,35 3,13 3,42 2,70 4,12 3,53
Aceptando en una forma conservadores el aporte de nevado de la tabla N° 05 con un promedio de 2.2m3/seg, que con relación a lo obtenido por Loroña (2011) es aproximadamente 13% mayor, valor muy aceptable en los trabajos de regionalización hídrica y se dan rangos mayores al 40%.
RETROCESO GLACIAR Se realizó el mapeo de la cobertura glaciar mediante el índice normalizado de nieve (NDSI) con la ayuda de las imágenes satelitales LANDSAT se determinó la superficie del glaciar por un periodo de 4 años, que se presenta en el Tabla N ° 07, a partir del cual se determinó una regresión lineal definiendo la Ecuación Nº 07 que nos permite definir la pérdida del glaciar en un periodo de 26 años. Tabla 07. Áreas de la cobertura glaciar del periodo de 1987 - 2007 Año 1987 1989 2007 2010
Tipo de Imagen Imagen satelital Imagen satelital Imagen satelital Imagen satelital
Path / Row 008/067 008/067 008/067 008/067
Área (km2) 34,91 34,71 24,24 23.11
Variación de Área (km2) ----0,19 10,47 1.13
Tasa de Retroceso (km2/Año) ------0,095 -0,582 -0,377
Aglaciar = -0.5393(T) + 1107 R² = 0.9958 ……………………. Ecuación
N° 07. Correlación Area glaciar con tiempo.
RETROCESO GLACIAR Se realizo el mapeo de la cobertura glaciar mediante el índice normalizado de nieve (NDSI) con la ayuda de las imágenes satelitales LANDSAT se determinó la superficie del glaciar por un periodo de 04 años, que se presenta en el Tabla N° 07 y en los mapas N° 06, 07, 08 y N° 09. Los glaciares de la Subcuenca de Quillcayhuanca muestran una pérdida de superficie de un 33.8 % durante el periodo en estudio. En este glaciar se pueden distinguir que en el primer periodo entre 1987 y 1989 hay un retroceso anual muy marcado y acelerado de 0.19 km2 /año, en el periodo de 1989 y 2007 se muestra mínimo aumento de 10.47 km2 /año de nevado, y el periodo entre 2007 y 2010 se obtuvo como resultado una pérdida de 1.13 km2 /año. A pesar de que existen diferencias entre los tres periodos, esas diferencias no son significativas desde un punto de vista estadístico.
RETROCESO GLACIAR A partir del volumen de glaciar determinado, se determinó una regresión lineal definiendo la Ecuación Nº 07 que nos permite definir la pérdida del glaciar en un periodo de 26 años.
Grafica 2.a.Correlacion del área de la cobertura glaciar en el Tiempo
RETROCESO GLACIAR Año 1987 1988 1989 1990 1991 1992 1993 1994 1995 1996 1997 1998 1999 2000
Área Glaciar (km2) 35.4 34.9 34.3 33.8 33.3 32.7 32.2 31.6 31.1 30.6 30.0 29.5 28.9 28.4
Año 2001 2002 2003 2004 2005 2006 2007 2008 2009 2010 2011 2012 2013
Área Glaciar (km2) 27.9 27.3 26.8 26.2 25.7 25.2 24.6 24.1 23.5 23.0 22.5 21.9 21.4
DISCUSIÓN La cartografía constituye la base principal para el estudio de la Evolución de la cobertura glaciar y la influencia en su Rendimiento Hídrico en la subcuenca del rio Quellcayhuanca, en este sentido, se consideró el mayor cuidado la automatización de la cartografía base habiendo utilizado la cartografía del Instituto geográfico Nacional (IGN) a una escala 1: 100000 la que fue automatizada utilizando las herramientas de sistemas de información geográfica con el programa ArcGis 10. Calculando un área integral de 248.8 Km2 con una perímetro 81.28 Km.
DISCUSIÓN Subcuenca del rio Quellcayhuanca, es uno de los tributarios de la cuenca del rio Santa y aporta recursos hídricos a la cuenca en mención, del mismo modo la subcuenca de Quillcay también es un tributario del rio santa y es un subcuenca vecina de la quebrada honda y por ende tiene características similares. F.Loroña (2011) determina la precipitación promedio anual histórica de 907 mm para la subcuenca de Quillcay, valor muy cercano obtenido en la Subcuenca del rio Quellcayhuanca (877.2 mm). Electro-Perú Volumen V (1980) determina para esta zona una precipitación de 820mm.
DISCUSIÓN Se determinó para la Subcuenca del rio Quellcayhuanca a partir de esta ecuación la evolución del nevado en el periodo comprendido de 1987 al 2012, encontrando que la perdida de nevado en este periodo es de 11,80km2, lo que representa una pérdida de 0.50 km2 /año; valores muy similares a la tasa de retroceso que determina F. Loroña (2011) que es de 0.57 km2 /año.
Se generó cartografía digital e información alfanumérica de Mapas topográfico, retroceso glaciar de la Subcuenca del rio Quellcayhuanca , que serán insumos de mucha utilidad para otras investigaciones.
Se determinó la evolución del nevado en el periodo comprendido de 1987 al 2012, encontrando que la pérdida de nevado en este periodo es de 11,80km 2 , lo 2 /añ o . que representa una pérdida de 0.50km
Instalar estaciones climatológicas en la subcuenca de la quebrada Honda. Instalar estaciónes hidrométrica en la subcuenca quebrada Honda. Instalación de una estación nivometrica
El Cambio climático acelera el retroceso glaciar. Glaciar Yanamarey (Cordillera Blanca - Altitud 4.786 msnm). Periodo 1982 - 2005. Fuente: INRENA De 1989 a 1997 hemos perdido 21,85% de superficie glaciar. (De 2 042 Km2 a 1 595 Km2)
Antes
Hace un año
