MODELLING FUTURE LAKES IN GLACIER BEDS: FIRST EXPERIENCES IN PERU Climate change has caused a dramatic reduction in glacier extent and volume. In Peru, the second National Glacier Inventory documents an overall loss in glacier-covered area of 42.64% in comparison to the total of 204.85 km 2 found in the first inventory of 970 (ANA, 202). Glacier retreat can lead to the formation of new lakes in exposed glacier-bed overdeepenings. These lakes can amplify natural hazards to downstream population but can also become tourist attractions and new possibilities for water management at local, regional or basin level (Haeberli et al., 206a; Haeberli et al., 207; Salzmann et al., 2004). As new and larger reservoirs of fresh water, these lakes could supply benefits to, for example, populations, hydropower and agriculture. Identifying sites of possible future lake formation is thus an essential step in early planning of disaster risk management measures, adaptation to climate change and integrated management of water resources (Haeberli et al., 206a). This planning should therefore already take into account the formation of new high-mountain landscapes with changing ecosystems in the long term. A pioneering work of compiling and locating possible future lakes in Peru was begun by the Unidad de Glaciología y Recursos Hídricos (Glacier and Water Resources Unit, UGRH) 2 of the Autoridad Nacional de Agua (National Water Authority, ANA), and continued by the Instituto Nacional de Investigación en Glaciares y Ecosistemas de Montaña (National Institute for Research in Glaciers and Mountain Ecosystems, INAIGEM) 3. The team of the Glaciares+ Project, initiative of the Global Climate Change Program of the Swiss Cooperation (COSUDE), has accompanied and assisted in this process. LOCATING NEW LAKES: METHODOLOGY The methodology shows the procedures applied for an approximate classification of the most probable time period in which future glacial lakes may form. A Digital Elevation Model (DEM) from the year 2000 (SRTM) with spatial resolution of 90 m was used in combination with glacier outlines identified in the second ANA National Glacier Inventory (ANA, 202). Estimates of possible future lakes with areas of over ha (exceeding the limits of SRTM resolution) were made in three steps using geographic information systems (GIS), as shown in figure (Frey et al., 200). Step I. Slope Classification Step II. Morphological criteria of glacier Step III. Numerical modelling Figure. Steps and criteria taken into account in delimiting future glaciar lakes in Peru (modified according to Colonia et al., 207; Frey et al., 200). This executive brief has been developed based on the scientific article of Colonia, D., Torres, J., Haeberli, W., Schauweckers, S., Braendle, E., Giraldez, C. & Cochachin, A. (207). Compiling an Inventory of Glacier-Bed Overdeepenings and Potential New Lakes in De-Glaciating Areas of the Peruvian Andes: Approach, First Results and Perspectives for Adaptation to Climate Change. Water, 9 (5), 336. (DOI: 0.3390/w9050336) 2 www.ana.gob.pe/gestion-agua/cambio-climatico/page/2 3 www.inaigem.gob.pe/
Step I involved mapping (using the DEM) glacier areas with surface slopes <0, as a pre-selection of sites with potential bed overdeepenings. This step was the basis for focusing the visual analysis of the three morphological indicators (step II, see figure 2). The number of fulfilled morphological criteria (MC) for each site allowed for the attribution of probability or confidence levels concerning the existence of inventoried glacier-bed overdeepenings. No crevasses Flat Slope Lateral narrowing Crevasses present. Flat section More inclined section 2. Wide section Narrow section 3. Area without crevasses Area with crevasses Area CM CM2 CM3 Probability Area A Low Area B Medium Area C High Figure 2. Step II of the methodology: Morphological criteria that indicate the existence of glacier-bed overdeepenings (modified according to Colonia et al., 207; Frey et al., 200). The results of the first two steps were compared to the results of the GlabTop model, which was used to estimate the distribution of ice thicknesses (Linsbauer et al., 202; Linsbauer et al., 2009) based on the analysis and combination of the DEM, the glacier outlines and branch lines (see figure 3). The estimate of ice thickness using GlabTop, with a average range of uncertainty of about ±30%, can only provide an approximation of the order of magnitude of depths and potential volumes of the lakes (Haeberli et al., 206b; Linsbauer et al., 202). GLABTOP MODEL Estimación de espesor de hielo. INPUTS. DEM 2. Glacier limits 3. Branch lines Contour line Glacier outline Flow lines Surface slopes Base points Interpolation of glacier bed Figure 3. Step III of the methodology: Numerical modeling (modified according to Linsbauer et al., 202).
An important aspect concerning lake safety is whether the potential new lakes will become dammed by bedrock or till/moraine would be blocking their water (Zemp et al., 2005). As a first approximation and estimation, clean glaciers with small or absent debris input from surrounding rock walls can be assumed to predominantly have rock beds, while heavily debriscovered glaciers tend to have very thick/elevatedmoraine beds (Haeberli & Fisch, 984). In order to get a rough impression concerning the possible time when the formation of such lakes could initiate, simple extrapolation schemes had to be used, which could be rapidly applied to the many unmeasured glaciers on the basis of inventory information. In case of Peru, with data gathered previously in the first and second inventory, the average annual horizontal and vertical change (along an assumed central flow line) of the lowest point (hmin) of each glacier was calculated. Both corresponding rates of change of the lowest point were extrapolated to the future for each individual glacier. Comparison with the calculated position and elevation of the modeled glacier-bed overdeepenings provided an indicative time for the onset of the of lake formation in case of continued change at constant rates. Values reflecting in a simple way possible future acceleration trends in global warming and glacier retreat were obtained by doubling the average rates of change and by correspondingly dividing the time to the onset of possible lake formation by 2. This acceleration scenario mirrors an increase in area losses from about % per year during past decades to about 2% per year in the near future. Based on estimates for constant and accelerated rates of change, a rough distinction was made between three possible periods for future lake formation: Lake formation already underway or imminent (within the coming about 0 years). Lake formation likely during the first half of the century (within about 0-40 years) Lake formation probable around mid-century or later (after 40 years or more) PERU: FIRST RESULTS Using the aforementioned methodology, 20 sites with glacier-bed overdeepenings exceeding 0 000 m2 were Andes North Centre and South South Basin Santa Marañón Pativilca Huallaga Mantaro Rímac Mala Centre CaÑete Ocoña Camaná Alto Apurímac 23 24 5 2 2 4 3 3 8 6 Urubamba 67 Alto Madre de Dios Inambari Azángaro Suches Total 29 5 7 20 found in the cordilleras of Peru (see table, table 2 and figure 4) with an estimated total future lake volume of about 260 million m3. Compared to the total glacier volume (38 km3) calculated for the years 2003-200 in the second national glacier inventory, the potential lake volume corresponds to about 0.5%-%. This can be explained by the limited extent of remaining flat glacier parts in the Peruvian cordilleras. Flat/clean glacier tongues - in contrast to more slowly downwasting debris-covered ice - have already mostly disappeared in the recent past. The higher number of possible future lakes in the Cordillera Vilcanota on a flat high-altitude plateau than in the Cordillera Blanca with its deeply cut valleys mirrors this effect. Classification Underway or imminent 33 First half of the century 43 Mid-century or later 25 Total 20 Table. Distribution of possible future lakes (>0,000 m2) by Table 2. of future lakes by possible initial time of river basin formation.
Area (km2) Area Cordillera Figure 4. Area and number of possible future lakes in the cordilleras of Peru. A number of mountain ranges exist where no further lake formation appears likely because the glaciers are located on steep slopes and cover small areas. Most of the new lakes will form in bedrock depressions and, hence, have stable dams. However, most of them will also form in the immediate neighbourhood of debuttressed lateral moraine/rock slopes and extremely steep bedrock peaks above 5000 m a.s.l. with warming hanging glaciers and degrading permafrost. This rapidly changing and destabilizing high-mountain environment causes the long-term probability of large rock or ice avalanches into such lakes to increase. Risks from impact and flood waves even for humans and their infrastructure at considerable distances downstream are, therefore, also systematically increasing. (Haeberli et al., 207). CONCLUSIONS The combination of terrain analysis, visual inspection of glacier morphology and geometry with numerical modelling were optimal techniques used for reaching realistic assessments including definitions of probability / confidence levels. With further glacier shrinking and even vanishing, the rate of lake formation can be expected to be on the decline. Most of the anticipated future lakes are likely to come into existence within the next few decades. The most robust predictions refer to the location and approximate area of glacier-bed overdeepenings; estimated morphometries relating to their exact shape, depth or volume are less certain, providing orders of magnitude rather than clearly defined values. The volume of 260 million m3 in 20 modeled bed overdeepenings represents only 0.5% to % of the presently still glacier volume. This small ratio results from the fact that most flat glacier parts have already disappeared in the investigated mountains, leaving mostly small and steeply inclined glaciers. The first results set out for glaciers in the Peruvian cordilleras are part of the initial process of data collection and inventory of possible future lakes nationwide. The data gathered in the second national glacier inventory (2003-200) enabled this first approximation. The progress made by Peru in current data gathering and glacier analysis is worth noting; without it, the methodology used would not have been possible.
RECOMMENDATIONS The input for the methodological process was current glacier data, which can be used in other and future related research leading to decision-making and the development of policies that are sustainable in the long term. It is therefore vital that this information be provided to the relevant authorities in order that they make practical use and application of the data. The results presented for Peru should be understood as a rough order-of-magnitude estimate. It is therefore necessary to apply more sophisticated procedures (e.g. radio-echo soundings) in case of locally-or regionally-enhanced interest. This information is crucial for achieving better estimates of future lake geometries. Given the probability that most new lakes of significant volume are forming now and within a few coming decades, long term planning should strongly emphasize water storage and use. With regard to projects to supply hydroelectric power or provide water, local research will complement plans for dam building. With the formation of new lakes there is the possibility of an increase in the vulnerability and exposure of the population and infrastructure located downstream to extreme events with negative impact. It is recommended that the relevant authorities carry out analyses and assessments of potential related risks, enabling them to establish a list of priorities that integrates the need for more research with actions towards disaster hazard management, primarily in high-mountain regions close to glaciers. The potential of new lakes for water storage and use should be suitably approached through a participatory process in an integrated water resources management, ensuring sustainability and recognizing the many users of water in order to avoid competition and conflict over the scarce resource in some regions. REFERENCES Autoridad Nacional del Agua (ANA). (202). Inventario Nacional de Glaciares y Lagunas. Autoridad Nacional del Agua: Lima. Colonia, D., Torres, J., Haeberli, W., Schauweckers, S., Braendle, E., Giraldez, C. & Cochachin, A. (207). Compiling an Inventory of Glacier-Bed Overdeepenings and Potential New Lakes in De-Glaciating Areas of the Peruvian Andes: Approach, First Results and Perspectives for Adaptation to Climate Change. Water, 9 (5), 336. https://doi.org/0.3390/w9050336 Drenkhan, F., Guardamino, L., Huggel, C. & Frey, H. (in review). Current and future glacier and lake assessment in the deglaciating Vilcanota-Urubamba basin, Peruvian Andes. Frey, H., Haeberli, W., Linsbauer, A., Huggel, C., & Paul, F. (200). A multi - level strategy for anticipating future glacier lake formation and associated hazard potentials. Natural Hazards and Earth System Science, 0 (2), 339-352. https://doi.org/0.594/nhess-0-339-200. Haeberli, W., & Fisch, W. (984). Electrical resistivity soundings of glacier beds: A test study on Grubengletscher Wallis, Swiss Alps. Journal of Glaciology, 30, 373-376. https://doi.org/0.07/s002243000006250 Haeberli, W., Bütler, M., Huggel, C., Lehmann, T., Schaub. Y., & Schleiss, A. (206 a). New lakes in deglaciating high - mountain regiones: Opportunities and risks. Climate Change, 39 (2), 20-24. https://doi.org/0.007/s0584-06-77-5 Haeberli, W., Linsbauer, A., Cochachin, A., Salazar, C., & Fischer, U. H. (206b). On the morphological characteristicas of overdeepenings in high - mountain glacier beds. Earth Surface Processes and Landforms, 4 (3), 980-990. https://doi.org/0.002/esp.3966 Haeberli, W., Schaub, Y., & Huggel, C. (207). Increasing risks related to landslides from degrading permafrost into new lakes in de - glaciating mountain ranges. Geomorphology. https://doi.org/0.06/j.geomorph.206.02.009 Linsbauer, A., Paul, F., & Haeberli, W. (202). Modeling glacier thickness distribution and bed topography over entire mountain ranges with GlabTop: Application of a fast and robust approach. Journal of Geophysical Research: Earth Surface, 7(3), - 7. https://doi.org/0.029/20jf00233
Linsbauer, A., Paul, F., Hoelzle, M., Frey, H., & Haeberli, W. (2009). The Swiss Alps without glaciers: A GIS-based modelling for reconstruction of glacier beds. Proceedings of the Geomorphometry, 243-247. Salzmann, N., Kääb, A., Huggel, C., Allgöwer, B., & Haeberli, W. (2004). Assessment of the hazard potential of ice avalanches using remote sensing and GIS-modelling. Norwegian Journal of Geography, 58(2), 74-84. http://dx.doi.org/0.080/002995040006805 Zemp, M., Kääb, A., Hoelzle, M., & Haeberli, W. (2005). GIS-based modelling of the glacial sediment balance. Zeitschrift für Geomorphologie, 38, 3-29. The development of this document was possible thanks to the Glaciares+ Project, an initiative of the Swiss-Peruvian cooperation within the framework of the Global Climate Change and Environment Program of the Swiss Agency for Development and Cooperation (SDC), executed by CARE Peru and the Swiss consortium led by the University of Zurich and including CREALP, METEODAT GmbH and EPFL. The project is carried out in close coordination with the National Water Authority (ANA), the Ministry of the Environment (MINAM) and the Center for Estimation, Prevention and Reduction of Disaster Risk (CENEPRED); and is implemented by the Water Resources and Glaciology Unit (UGRH) of ANA, regional governments of Ancash, Cusco and Lima, local governments and universities. Partners PG+ Download and view the article Compiling an Inventory of Glacier- Visualice y descargue el Bed Overdeepenings and Potential New Lakes in De-Glaciating Resumen Ejecutivo en Areas of the Peruvian Andes: Approach, First Results and español. Perspectives for Adaptation to Climate Change.