Hostname: page-component-78c5997874-g7gxr Total loading time: 0 Render date: 2024-10-31T23:30:59.326Z Has data issue: false hasContentIssue false

Recent retreat of the Elbrus glacier system

Published online by Cambridge University Press:  14 March 2016

IULIAN-HORIA HOLOBÂCĂ*
Affiliation:
Department of Physical and Technical, Faculty of Geography, Babeş-Bolyai University, 5-7 Clinicilor, Cluj-Napoca, Romania
*
Correspondence: Iulian-Horia Holobâcă <holobaca@geografie.ubbcluj.ro>
Rights & Permissions [Opens in a new window]

Abstract

The glacier system covering Europe's highest mountain, Elbrus, has exhibited an accelerated retreat since 1980. Some studies have related this retreat to a significant summer temperature increasing trend. Relief- and aspect-related parameters for the glacierized area have an important impact on glacier changes. In this paper, the changes in glacier area are identified, quantified and correlated with relief parameters for the period 1985–2007. Spatial analysis was performed using the GLAM-CD (Glacier Mapper – Change Detector) algorithm. The input data for this algorithm were Landsat 5 images, the Aster Global Digital Elevation and the glacier outlines from the GLIMS project (Global Land Ice Measurements from Space). Regression analyses between glacier area losses and relief-related parameters indicate a significant positive relation with the altitude and a significant negative relation with the glacier surface area. In this context, we used a correlated component regression to model these relations. The model explains >50% of the total variation.

Type
Papers
Creative Commons
Creative Common License - CCCreative Common License - BYCreative Common License - NCCreative Common License - SA
This is an Open Access article, distributed under the terms of the Creative Commons Attribution-NonCommercial-ShareAlike licence (http://creativecommons.org/licenses/by-nc-sa/4.0/), which permits non-commercial re-use, distribution, and reproduction in any medium, provided the same Creative Commons licence is included and the original work is properly cited. The written permission of Cambridge University Press must be obtained for commercial re-use.
Copyright
Copyright © The Author(s) 2016

1. INTRODUCTION

The decline of mountain glaciers since the Little Ice Age (LIA) is an important topic in glaciology and related sciences. Mountain glaciers are vital sources of water for many regions in the world. Also, the risks associated with rapid ice melting can affect local communities.

Because glacier mass balance is very sensitive to climate variability, the dynamics of glaciers is a key indicator of low-frequency climate change (Oerlemans, Reference Oerlemans2005; Chinn and others, Reference Chinn, Winkler, Salinger and Haakensen2005). Worldwide, glaciers have retreated during the 20th century in response to global warming. Many regional analyses indicate a strong specific negative mass balance since 1990, with rapid losses especially in Patagonia, Alaska (Dyurgerov and Meier, Reference Dyurgerov and Meier2005) and the European Alps (Zemp and others, Reference Zemp, Frauenfelder, Haeberli and Hoelzle2005).

Caucasus glaciers have retreated since the mid-19th century (Solomina, Reference Solomina2000). Zolotarev and Kharkovets (Reference Zolotarev and Kharkovets2010) analysed the evolution of Elbrus glacierization since the LIA using cartographic and proxy methods, providing valuable information about past retreat rates.

The beginning of the 1980s was a landmark in glacial evolution in the Caucasus. Regional studies reported an accelerated retreat of Greater Caucasus glaciers since 1980 (Popovnin and Rozova, Reference Popovnin and Rozova2002; Stokes and others, Reference Stokes, Gurney, Shahgedanova and Popovnin2006). Some studies related this retreat to an increase in summer temperature (Shahgedanova and others, Reference Shahgedanova, Stokes, Gurney and Popovnin2005) or to a significant summer warming trend and a decline in winter precipitation (Holobâcă, Reference Holobâcă2013). Other studies highlighted potential consequences of the rapid glacial melt, such as glacial surges (Kotlyakov and others, Reference Kotlyakov, Rototaeva, Desinov, Zotikov and Osokin2002), multi-phase mass movements (Petrakov and others, Reference Petrakov, Chernomorets, Evans and Tutubalina2008), and catastrophic detachment and high-velocity long-runout flow (Evans and others, Reference Evans2009).

In a previous paper (Holobâcă, Reference Holobâcă2013), we found that climatic variability imposes the general trend for the Elbrus recent glaciation but relief parameters such as altitude, slope and aspect are very important for the glacial dynamics at the local scale. In this paper, we make a further contribution by relating the identified changes to relief parameters. Another important aim of this paper is to model these relations for future applications.

1.2. Study area

Situated in the northwestern part of the Greater Caucasus Mountain, Mount Elbrus is as an excellent candidate for local scale analyses. A well-developed ice cap covers the top of a dormant stratovolcano with two cones (Fig. 1). Mount Elbrus (5642 m), the highest peak in the Caucasus and Russian Federation, is covered by ~110 km2 (in 2007) of ice and perennial snow (Table 1). The ice thickness surpasses 240 m on the Western Plateau (Mikhalenko and others, Reference Mikhalenko, Kutuzov, Lavrentiev and Kunakhovich2008).

Fig. 1. Mount Elbrus glacier system outlines, superposed over the 2007 colour composite of Landsat bands 5, 4, 3.

Table 1. Change of the glaciers surface area, relative losses and rate of change (1985–2007)

Notes: The average values for small glaciers (surface area <1.5 km2), medium glaciers (1.5 km2< surface area>10 km2) and large glaciers (surface area >10 km2) are printed in italic bold; the average values for Elbrus area are printed in bold.

According to the World Glacier Monitoring Service (WGMS) inventory, there are 24 valley glaciers radiating from this ice cap. The surface area of valley glaciers ranges from <0.1 to >22 km2 (Table 1). The mean altitude of the Elbrus glacier system is about 4000 m, but the valley glaciers descend down to 2700 m in the southern and southeastern part (Table 2).

Table 2. Change of glacier altitude (1985–2007)

Notes: Median values for small glaciers (surface area <1.5 km2), medium glaciers (1.5 km2< surface area>10 km2) and large glaciers (surface area >10 km2) are printed in italic bold; the average values for Elbrus area are printed in bold.

The Elbrus valley glaciers, like other Caucasus glaciers, are affected by surge episodes (Kotlyakov and others, Reference Kotlyakov, Rototaeva, Desinov, Zotikov and Osokin2002). The associated rapid increases in ice flow rates can generate catastrophic detachments (Evans and others, Reference Evans2009), and can be a major risk for downstream human structures and habitat. In addition, the Elbrus glacier system is an important source of freshwater for the North Caucasus region. Ice melt feeds the headwater of two major rivers, Kuban and Terek.

2. METHOD AND DATA

2.1. Spatial data

One of the major issues for our multitemporal surface analysis was to ensure the comparability of data. Only representative spatial and temporal data have been chosen. The length of the study period (1985–2007) was also an important concern. The recession trend observed for all Elbrus glaciers during the analysed period is consistent with climatic warming trends. The study period of 22 years appears to be long enough to avoid complications caused by non-climate related glacier perturbations (e.g. surge episodes). Also, Zolotarev and Kharkovets (Reference Zolotarev and Kharkovets2010) used a time period of 22 years (1957–79) to study Elbrus glaciers. This offers a good opportunity for comparison of results.

The Landsat 5 satellite images used in this study were acquired by the Thematic Mapper (TM) sensor, with 30 m equivalent pixel size. Both images were obtained from the end of the ablation season (late-July 1985 and early August 2007), when the maximum area of snow-free ice is exposed (Fountain and others, Reference Fountain, KrimmeI and Trabant1997). The scenes have minimal snow and cloud cover and were acquired before the first snowfall of winter. The source of satellite data is the United States Geological Survey Global Visualization Viewer (GLOVIS) Landsat archive (http://glovis.usgs.gov/). The orthorectified satellite scenes were co-registered using the ERDAS 9.1 software (Intergraph Corporation, Madison, AL, USA). The horizontal RMS error was <3.0 m.

Ice divides were located using Global Land Ice Measurement from Space (GLIMS) project data. This database is accessible at http://glims.colorado.edu/glacierdata/, and stores multitemporal glacier outlines and metadata (processing technics, scalar data and related literature) (Raup and others, Reference Raup2007).

The relief parameters (altitude, slope and aspect) used were derived from an Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER) global DEM. The main advantage of this DEM, available at http://www.ersdac.or.jp/GDEM/E/ is that it has the same 30 m spatial resolution as the Landsat 5 TM images. The DEM was registered using heights from 1:50 000 topographic maps, accessible on the World Wide Web at http://mapstor.com/mapsets/russia-maps.html and 22 ground control points (GCPs) obtained with a GPS receiver in the field. The registration of the DEM was made using ERDAS 9.1 software (Intergraph Corporation, Madison, AL, USA). The method uses a second-order polynomial transformation, which allows a shift across the image that follows the curve of the second-order polynomial. The horizontal RMS error was in the range of 15 m (0.5 of DEM resolution).

2.2. Spatial analysis

Spatial analysis was carried out in multiple steps (Fig. 2): (1) derivation of the Elbrus glacier outlines from multitemporal satellite scenes using the Glacier Mapper – Change Detector (GLAM-CD) algorithm; (2) definition of the ice divides for individual glaciers from GLIMS data; (3) identification of the glacier outlines for individual glaciers combining two categories of spatial data (glacier–non-glacier outlines and ice divides); and (4) identification of the spatial patterns of the glacier changes.

Fig. 2. Automatic extraction and surface analysis algorithm.

2.2.1. Derivation of the Elbrus glacier outlines

The GLAM-CD is a logical structured algorithm that allows evaluation of relief-related changes in the glacial area from two time-comparable satellite images (Holobâcă, Reference Holobâcă2013). This ARC GIS toolbox uses an image ratio method in order to enhance the variation in the slope of the spectral reflectance curves between two different spectral ranges. The GLAM-CD has two tools, one for proper ratio images and the other for the Normalized Difference Snow Index (NDSI) calculation. The NDSI, a widely used index for automated glacier delineation (Paul and others, Reference Paul, Huggel, Kääb, Kellenberger and Maisch2002), was employed in this analysis. As snow strongly reflects visible light (e.g. green) but absorbs middle infrared light, NDSI is high over snow and ice. Acceptable snow cover results can be found with NDSI thresholds in the range, 0.25–0.45 (Hall and others, Reference Hall, Riggs and Salomonson1995). In this study, glacier outlines are derived using a threshold value of 0.4, which provides better results for glacier delineation.

We use a semiautomatic method to extract glacier outlines from the Landsat 5 images. The initial glacier outlines obtained with the GLAM-CD tool were checked visually using the colour composite of Landsat bands 5, 4, 3 and 3, 2, 1. This visual correction for classification errors, debris covered ice and shadow analysis is important due to the limitations of full automated methods (Paul and others, Reference Paul2013).

2.2.2. Definition of ice divides for individual glaciers

Following algorithms developed by Manley, GLIMS uses watershed commands in ArcInfo to separate individual glaciers from large contiguous ice masses (Raup and others, Reference Raup2007). In this multitemporal study, the same ice divides have been used for individual glacier definition. IDs or names from the GLIMS database were adopted for glaciers rather than the local name from topographic maps. In order to avoid misidentifications we also employed the WGMS IDs. The ice divides for individual glaciers were checked visually, and manually improved using colour composite Landsat bands (5, 4, 3 and 3, 2, 1). Special attention was given to internal rocks in the ice divide area.

2.2.3. Identification of glacier outlines for individual glaciers

The individual glacier outlines have been obtained combining two categories of vector spatial data. These vector data were the Elbrus glacier outlines and the individual glacier ice divides.

2.2.4. Identification of the spatial patterns of glacier changes

We analysed the vector files from the two-time positions of the glacier outline to identify changes in the glacierized area.

The relative losses of individual glaciers were correlated with the characteristics of the relief parameters (Table 1). We then tested the correlation coefficient at the 0.05 significance level. In addition, a Correlated Component Regression Linear Model (CCR.LM) was employed to model the relation between glacier changes and relief parameters. This regression method uses fast cross-validation to determine the amount of regularization in order to produce reliable predictions from data with correlated explanatory variables. We use XLSTAT software (Addinsoft, USA) to build the model.

3. RESULTS

Recession of glaciers was the main characteristic during the analysed period. The size of the change varied dramatically from glacier to glacier (Tables 1 and 2). Climate is thought to be the main driver for the changes, but local conditions can alter the response of individual glaciers to the warming trend. In order to better understand these complex relations, a quantitative approach is presented below.

3.1. Glacier areal change

Table 1 indicates that glacier size is important for the intensity of the recession. The statistical relation is negative and significant at 0.05 p-level (Table 3). Small glaciers, with an area <1.5 km2, dramatically diminished when considering relative changes. The relative losses for these glaciers with an average value of 32.2% and ranging from 12 to 64%, surpassed the average value for the Elbrus glacier system as a whole (8%).

Table 3. Correlation coefficient values (R) between glacier losses and relief parameters

Notes: Correlation values statistically significant at 0.05 level of signification are printed in bold, min, minimum; max, maximum; S25%, surface slope averaged over the lowermost 25% of the glacier; E, easting; N, northing.

For the most part, for the medium-sized glaciers, with an area ranging from 1.5 to 10 km2, glacier losses were <8%. The average value for the medium glaciers was 7.8%. The most frequent range of values was 2.5–3.6%. We noticed two major exceptions: Garabashi Glacier, with the smallest relative loss value within the Elbrus glacier system (0.8%) and Glacier 576 (SU4H08004316), which lost 39% of its surface area in the analysed period. The Garabashi Glacier is situated in the southeast part of Mount Elbrus. The redistribution of snow from the northwest to the southeast side of the mountain range by the dominant north-westerly winds can explain this small glacier relative loss.

The values of relative losses for glaciers with an area >10 km2 were greater than expected with an average value of 7.0%. Dzhikiugankez Glacier, the most extensive glacier in the study area (22 km2), lost 19% of its surface area. The relative loss values for Ulluchiran (5.7%) and Bol'shoy Azau (9.3%) glaciers were close to the value for the overall glacierized area. Only Malyi Azau Glacier had a relative loss in the expected range (1.8%). The important relative losses for the Dzhikiugankez Glacier can be explained by the role of the intense post-volcanic activity, especially the influence of thermal and fluid flows in the north-eastern part of the Elbrus volcano (Masurenkov and Sobisevich, Reference Masurenkov and Sobisevich2012).

3.2. Glacierized altitude change

Many glaciological studies suggest that altitude is the most important independent variable and is a fundamental control of both temperature and precipitation (Evans, Reference Evans and Chorley1969; Sugden and John, Reference Sugden and John1969; Flint, Reference Flint1971; Price, Reference Price1981). Changes in altitude also indicated a generalized glacial retreat for all the Elbrus glaciers. The results had some spatial patterns, exhibited in Table 2.

Change in minimum altitude is an objective indicator of frontal glacier retreat. The largest changes, ranging from −93 to −131 m, characterized the glaciers that descended below 3000 m. The loss in minimum altitude for the glaciers situated above 3000 m was complex. The loss values greater than the Elbrus average change (−47 m) were observed for medium and large glaciers that descended at low altitude (below 3200 m). The average values for these glaciers were −46.1 and −71.4 m, respectively. For small glaciers, altitudinal losses >47 m have been observed over a large altitudinal range (3532–3743 m). The medium range of values (−43 to −49 m) characterized the small glaciers with an average value of −34.5 m. The exception was the medium Karachaul glacier, which lost 45 m in minimum altitude. Minor losses in minimum altitude, sometimes below the DEM resolution, characterized the smallest glaciers, but were also observed for some medium and large glaciers (Table 2).

The average altitude of the Elbrus glaciers rose during the analysis period. The average altitudinal change for the Elbrus glacier system was −29 m. Changes >−30 m characterized the glaciers with large area losses (Table 1) having west, northwest, southwest and northeast-facing slopes (Fig. 1). Bityuktyube Glacier had the biggest change in average altitude (−81.9 m).

For the most part, the medium range of average altitude change (from −20 to −30 m) characterized the medium and large glaciers with minor glacial area losses. The only exception was the 310 glacier with an area of 0.1 km2.

The class of minor average altitudinal change values ranged between −0.4 and −20 m. Only three glaciers (532v, Irik and Garabashi) with an area >1 km2 had minor average altitudinal losses. All these glaciers are situated in the south-eastern part of the Elbrus massive. On the other hand, the glaciers with area <1 km2 had the smallest average altitude losses.

Fewer changes were observed at maximum altitude. Most glaciers had no change in maximum altitude during the analysis period (Table 2). Some of them experienced a decrease in altitude. The maximum decrease (−69 m) was observed for the 565a glacier. Small increases in maximum altitude, ranging from 2 to 17 m, were identified for three glaciers. The changes in the maximum altitude seem to be random and uncorrelated with other relief parameters.

3.3. Modelled relations between relief parameters and glacier losses

Table 3 presents correlation coefficients between glacier loss and relief parameters. The relative glacier losses were correlated with initial relief parameters for individual glaciers, e.g. surface altitude, slope and aspect (easting and northing). Correlation coefficients were tested at the 0.05 significance level.

As expected, glacier surface area has a negative correlation with glacier loss. Significant negative correlations between these two parameters were found during the studied period (Table 3). Generally, the small glaciers have been more affected by glacial retreat. Also, we found significant correlations between glacier loss and altitude parameter. Our analysis revealed a significant positive correlation between glacier loss and minimum altitude. This relation can be explained by the position of the small glaciers (with high loss) situated at higher altitude than the medium or larger glaciers. The medium or larger glaciers, particularly the glaciers situated in the south or south-eastern part of the Elbrus massive had lower minimum altitudes and minor glacier losses (Tables 1 and 2).

Change in glacial surface area and the average and maximum altitude are negatively correlated. The glacier losses were higher at low altitudes and the glaciers with high average or maximum altitude have been less affected by deglaciation. Correlation coefficients, ranging from −0.5 to −0.6, were statistically significant.

The difference between maximum and minimum altitude has a significant negative relation with relative glacier loss. This variable had the highest correlation coefficient value (−0.7) indicating a higher climatic sensitivity of the glaciers with smallest altitudinal range.

Of the parameters tested against glacier loss, the lowest (magnitude) correlation coefficients, ranging from −0.3 to 0.2, were observed for maximum slope and average slope. Regression of the surface slope averaged over the lowermost 25% of the glaciers indicated a correlation coefficient, 0.2 which was not statistically significant. On the contrary, the minimum slope had a significant correlation with glacier loss. In our opinion, this statistically significant correlation can be explained by the size of the glaciers. The medium and large glaciers with minor glacier losses had minimum slope values of zero (or very close to zero) and the small glaciers with significant glacial surface area losses had larger minimum slope values.

The aspect parameters included in the model were easting and northing. The easting had a direct relation with relative glacier losses (0.3), but the relation between changes in glacial surface area and northing becomes negative (−0.3). The results seem to indicate a preferred sheltered orientation of the glaciers from direct solar radiation.

A correlated component multiple linear regression model was built in order to relate the relief parameters to glacier losses. The model uses 11 variables, listed in Table 3. The dependent variable is relative glacial loss for individual glaciers.

The model equation is:

(1) $$\eqalign{{\rm PGL}(\% ) & = 113.34 - 10^{ - 6 \ast} S1985 + 0.005^{\ast} {\rm altmin} \cr &\quad - 0.004^{\ast}{\rm altmax} - 0.027^{\ast}{\rm altaver} \cr &\quad + 0.0006^{\ast}({\rm altmax} - {\rm altmin}) + 5.59^{\ast}{\rm slopemin} \cr &\quad + 0.11^{\ast}{\rm slopemax} + 0.49^{\ast}{\rm slopemed} \cr &\quad - 0.50^{\ast}S25\% + 6.34^{\ast}N - 0.42^{\ast}E} $$

where PGL is relative glacier loss, S1985 is the surface area (km2), altmin is the minimum altitude (m), altmax is the maximum altitude (m), altaver is the average altitude (m), slopemin is the minimum slope (°), slopemax is the maximum slope (°), slopemed is the average slope (°), S25% is the surface slope (°) averaged over the lowermost 25% of the glacier, N is northing and E is easting.

The standardized (β) and unstandardized (B) coefficients of the multiple regression are presented in Table 4. The standardized coefficients indicate which of the independent variables has greater effect on the dependent variable in a multiple regression (Schroeder and others, Reference Schroeder, Sjoquist and Stephan1986). The results indicate that average altitude (β = −0.39) and minimum slope (β = 0.43) have the major impact on the relative glacier losses.

Table 4. Unstandardized (B) and standardized (β) coefficients of the multiple regression

Notes: min, minimum; max, maximum; S25%, surface slope averaged over the lowermost 25% of the glacier; E, easting; N, northing.

The regression does a good job of modelling glacier losses. The adjusted r 2 has a value of 0.67. Therefore, more than half the variation in relative loss is explained by the model. Also, the cross-validation value of the model is 0.14. Two usual tests for heteroscedasticity (Breusch and Pagan, Reference Breusch and Pagan1979; White, Reference White1980) were performed on model residuals. The significance level was 0.05. The results of both tests indicate that the null hypothesis (H 0: Residuals are homoscedastic) cannot be rejected at the selected significance level.

4. DISCUSSION

Our technique, combining semiautomatic extraction of glacier outlines with surface analysis has some major advantages. Automatic delineation of glacier outlines on a TM scene is up to 300 times faster than manual techniques (Paul and others, Reference Paul2013) and even our semiautomatic approach is much faster. In addition, there is apparently little loss of accuracy associated with automatic glacier outline delineation. Relative area differences between automatic and manual delineation techniques range from 2 to 5% (Paul and Kääb, Reference Paul and Kääb2005; Bolch and Kamp, Reference Bolch, Kamp, Kaufmann and Sulzer2006; Andreassen and others, Reference Andreassen, Paul, Kääb and Hausberg2008). In our case, differences fall within this range. Larger values have been observed for large and complex glaciers (e.g. Dzhikiugankez glacier – 3.4%, Irik glacier – 2.1%, Kyukyurtlyu glacier – 1.7%, etc.), while the difference was null or insignificant (<0.5%) for the small glaciers.

The integration of individual glacier ice divides with the GLIMS database considerably reduces processing time. Delineation of ice divides requires multiple spatial data sets (detailed DEM, high definition remote sensing data, aerial photography, field survey, etc.) and uses a complex algorithm. Use of the GLIMS ice divides also assures compatibility of multitemporal data. In fact, the reduction, or even elimination, of consistency problems is one of the most important GLIMS aims (Raup and others, Reference Raup2007). Also, automatic standard algorithms can facilitate meaningful comparisons between multitemporal analyses within the same region (Bishop and others, Reference Bishop2004).

Debris-covered areas of glaciers remain the most important source of systematic error in automatic glacier outline delineation. This issue can be solved using semi-automated methods proposed by Paul and others (Reference Paul, Huggel and Kääb2004), utilizing complex processing techniques and an accurate DEM. The obtained data are manually edited. More recent, accurate debris-covered glacier data have been acquired using summer coherence images from microwave sensors (Atwood and others, Reference Atwood, Meyer and Arendt2010; Frey and others, Reference Frey, Paul and Strozzi2012).

Spatial analyses have been designed to detect changes in glaciated areas and to correlate these changes with relief parameters. Our findings support previous studies (e.g. Popovnin and Rozova, Reference Popovnin and Rozova2002; Shahgedanova and others, Reference Shahgedanova, Stokes, Gurney and Popovnin2005; Stokes and others, Reference Stokes, Gurney, Shahgedanova and Popovnin2006 etc.) that reported accelerated retreat of glaciers in the Caucasus Mountains. This is more evident when the relative rate of change of glacier surface area (% a−1) in the analysis period is compared with the Zolotarev and Kharkovets (Reference Zolotarev and Kharkovets2010) data (Fig. 3). The relative rate of change was calculated for three previous periods: 1850–87, 1887–1957 and 1957–79 and computations were made employing the absolute values of the glacier surface area change for the total Elbrus glaciation (source of data: Zolotarev and Kharkovets, Reference Zolotarev and Kharkovets2010, p. 24, Table 3).

Fig. 3. Retreat trends in Elbrus Mountain since 1850 (% a−1). Data source: Zolotarev and Kharkovets (Reference Zolotarev and Kharkovets2010); p. 24, Table 3

The relative rate of change of glacier surface area in our analysis period exhibits the largest decrease since the LIA in Elbrus Mountain. The rate of change of glacier's surface area from 1985 to 2007 for the Elbrus glacier system (−0.37% a−1) is more than double that during the period 1957−97 (−0.16% a−1). The actual rate is even greater when compared with the change rate of −0.15% a−1 during the period 1887–1957 and almost double the rate of −0.20% a−1 observed in the 1850–87 period (Fig. 3).

Surface analysis of glacial changes confirms the findings of previous studies (e.g. Shangguan and others, Reference Shangguan2006; Lemke and others, Reference Lemke, Solomon, Qin, Manning, Chen, Marquis, Averyt, Tignor and Miller2007; Frauenfelder and Kääb, Reference Frauenfelder and Kääb2009; Bolch and others, Reference Bolch, Menounos and Wheate2010; Le Bris and others, Reference Le Bris, Paul, Frey and Bolch2011 etc.), that small glaciers (area <1.5 km2) and low glaciated areas (below 3200 m) are most affected by deglaciation. The major exception was Dzhikiugankez glacier, the larger Elbrus glacier, which exhibits high loss. This exception may, in part, be due to factors other than climate or relief, e.g. endogenic factors. Masurenkov and Sobisevich (Reference Masurenkov and Sobisevich2012) found that intense post volcanic activity, especially the influence of thermal and fluid flows, was responsible for glacial loss in the north-eastern part of Mount Elbrus.

The results of the regression analysis indicated statistically significant relationships between relative glacier losses and glacier surface area and altitude. Correlations with the average and maximum slope and with easting and northing of the aspect were not statistically significant. Only the minimum slope had a statistically significant relation with the melted glacial surface area. This situation may indicate the state of the dynamic equilibrium of the glaciers. The glaciers loosing area least have near zero minimum slope values, while the glaciers affected more by deglaciation have non-zero slope values.

Previous studies indicate some issues when using relative area changes as a parameter for inventory change (Nuth and others, Reference Nuth2013). The most important reported problems are the dependence of this parameter on the terrain geometry and heteroscedasticity of the variability for the glaciers having different sizes. In our approach, the correlated component regression linear model related the parameters of the relief with the relative area change. Also, the heteroscedasticity test revealed that the residuals of the model are homoscedastic.

5. CONCLUSION

The results of the change analysis indicate a general glacier recession trend of 0.4% a−1 during the studied period. This recession is very uneven among the Elbrus glaciers. In our opinion, this situation is caused by the impact of local geographical factors superposed on the general climatic driver of change.

A technique that combines a semiautomatic extraction of glacier outlines with a surface analysis is a fast and consistent method of detecting glacial change. These qualities can be very useful if the analysis is extended to the regional scale.

We developed a regression model between relative glacier loss and the relief parameters, which explains more than 67% of the total variation. These results can help in water supply planning and management of rapid glacial melt-associated risks in the Mount Elbrus area. It can also lead towards further models that can predict future glacial area evolution at both local and regional scale.

Acknowledgements

I gratefully acknowledge the financial support and generosity of The International Glaciological Society (IGS), without which the present study could not have been published.

References

REFERENCES

Andreassen, LM, Paul, F, Kääb, A and Hausberg, JE (2008) Landsat derived glacier inventory for Jotunheimen, Norway, and deduced glacier changes since the 1930s. Cryosphere, 2(2), 131145 (doi: 10.5194/tc-2-131-2008)Google Scholar
Atwood, DK, Meyer, F and Arendt, A (2010) Using L-band SAR coherence to delineate glacier extent. Can. J. Remote Sens., 36(1), 186195 (doi: 10.5589/m10-014)Google Scholar
Bishop, MP and 16 others (2004) Global land ice measurements from space (GLIMS): remote sensing and GIS investigations of the Earth's cryosphere. Geocarto Int., 19(12), 5784 Google Scholar
Bolch, T and Kamp, U (2006) Glacier mapping in high mountains using DEMs, Landsat and ASTER data. In Kaufmann, V and Sulzer, W eds. Proceedings of the 8th International Symposium on High mountain Remote Sensing Cartography, 21–27 March 2005, Karl Franzens University, La Paz, Bolivia, Graz, 1324 (Grazer Schriften der Geographie und Raumforschung 41)Google Scholar
Bolch, T, Menounos, B and Wheate, R (2010) Landsat-based inventory of glaciers in western Canada, 1985–2005. Remote Sens. Environ., 114(1), 127137 (doi: 10.1016/j.rse.2009.08.015)Google Scholar
Breusch, T and Pagan, A (1979) A simple test for heteroscedasticity and random coefficient variation. Econometrica, 47, 12871294 (doi: 10.2307/1911963)Google Scholar
Chinn, TJH, Winkler, S, Salinger, MJ, and Haakensen, N (2005) Recent glacier advances in Norway and New Zealand: a comparison of their glaciological and meteorological causes. Geografiska Annaler, 87A(1), 141157 (doi: 10.1111/j.0435-3676.2005.00249.x)Google Scholar
Dyurgerov, M and Meier, MF (2005) Glaciers and the Changing Earth System: A 2004 Snapshot. In: Occasional Paper 58, Institute of Arctic and Alpine Research, University of Colorado, Boulder, COGoogle Scholar
Evans, IS (1969) The geomorphology and morphometry of glacial and nival areas. In Chorley, RJ ed. Water, Earth and Man, London, Methuen, 369380 Google Scholar
Evans, SG and 6 others (2009) Catastrophic detachment and high-velocity long-runout flow of Kolka Glacier, Caucasus Mountains, Russia in 2002. Geomorphology, 105, 314321 (doi: 10.1016/j.geomorph.2008.10.008)Google Scholar
Flint, RF (1971) Glacial and quaternary geology. John Wiley and Sons, New York Google Scholar
Fountain, AG, KrimmeI, RM and Trabant, DC (1997) A Strategy for Monitoring Glaciers. U.S. Geological Survey Circular 1132, 16 pGoogle Scholar
Frey, H, Paul, F and Strozzi, T (2012) Compilation of a glacier inventory for the western Himalayas from satellite data: methods, challenges and results. Remote Sens. Environ., 124, 832843 (doi: 10.1016/j.rse.2012.06.020)Google Scholar
Frauenfelder, R and Kääb, A (2009) Glacier mapping from multi-temporal optical remote sensing data within the Brahmaputra river basin. In: Proceedings of the 33rd International Symposium on Remote Sensing of Environment, Stresa, Italy, 4–8 May 2009, International Center of Remote Sensing of Environment, Tucson, Arizona, Paper 299, 4Google Scholar
Hall, DK, Riggs, GA and Salomonson, VV (1995) Development of methods for mapping global snow cover using moderate resolution imaging spectroradiometer data. Remote Sens. Environ., 54(2), 127140 (doi: 10.1016/0034-4257(95)00137-P)CrossRefGoogle Scholar
Holobâcă, IH (2013) Glacier Mapper – a new method designed to assess change in mountain glaciers. Int. J. Remote Sens. 34(23), 84758490 (doi: 10.1080/01431161.2013.843804)CrossRefGoogle Scholar
Kotlyakov, VM, Rototaeva, OV, Desinov, LV, Zotikov, IA and Osokin, NI (2002) Causes and effects of a catastrophic surge of Kolka Glacier in the Central Caucasus. Z. Gletscherkd. Glazialgeol., 38(2), 117128 (doi: 10.3189/172756504781829710)Google Scholar
Le Bris, R, Paul, F, Frey, H and Bolch, T (2011) A new satellite derived glacier inventory for Western Alaska. Ann. Glaciol., 52(59), 135143 (doi: 10.3189/172756411799096303)Google Scholar
Lemke, P and 10 others (2007) Observations: changes in snow, ice and frozen ground. In Solomon, S, Qin, D, Manning, M, Chen, Z, Marquis, M, Averyt, KB, Tignor, M and Miller, HL eds. Climate change 2007: the physical science basis. Contribution of Working Group I to the Fourth Assessment Report of the Intergovernmental Panel on Climate Change. Cambridge University Press, Cambridge, United Kingdom and New York, NY, USA, 339378 Google Scholar
Masurenkov, YuP and Sobisevich, AL (2012) Influence of thermal and fluid flows on the dynamics of Elbrus glaciers. Dokl. Earth Sci., 442(1), 144147 (doi: 10.1134/S1028334X12010175)CrossRefGoogle Scholar
Mikhalenko, VN, Kutuzov, S, Lavrentiev, II and Kunakhovich, MG (2008) Paleoclimatic ice core studies in the Caucasus Mountains as a base for a new paleoclimatic reconstruction for the low latitudes of the Northern Hemisphere: prospects and first results. Geophys. Res. Abstr., 10, 10531054 Google Scholar
Nuth, C and 7 others (2013) Decadal changes from a multi-temporal glacier inventory of Svalbard, Cryosphere, 7, 16031621 (doi: 10.5194/tc-7-1603-2013)Google Scholar
Oerlemans, J (2005) Extracting a climate signal from 169 glacier records. Science, 308, 675677 (doi: 10.1126/science.1107046)CrossRefGoogle ScholarPubMed
Paul, F and Kääb, A (2005) Perspectives on the production of a glacier inventory from multispectral satellite data in Arctic Canada: Cumberland Peninsula, Baffin Island. Ann. Glaciol., 42(1), 5966 (doi: 10.3189/172756405781813087)Google Scholar
Paul, F, Huggel, C and Kääb, A (2004) Combining satellite multispectral image data and a digital elevation model for mapping debris-covered glaciers. Remote Sens. Environ., 89(4), 510518 (doi: 10.1016/j.rse.2003.11.007)Google Scholar
Paul, F, Huggel, C, Kääb, A, Kellenberger, T and Maisch, M (2002) Comparison of TM-derived glacier areas with higher resolution data sets. In: Proceedings of EARSeL–LISSIG-Workshop – Observing our Cryosphere from Space. European Association of Remote-Sensing Laboratories, Special Interest Group Land Ice and Snow, Bern, 1521 Google Scholar
Paul, F and 19 others (2013) On the accuracy of glacier outlines derived from remote sensing data. Ann. Glaciol., 54(63), 171182 (doi: 10.3189/2013AoG63A296)Google Scholar
Petrakov, DA, Chernomorets, SS, Evans, SG, and Tutubalina, OV (2008) Catastrophic glacial multi-phase mass movements: a special type of glacial hazard. Adv. Geosci. Eur. Geosci. Union (EGU), 14, 211218 (doi: 10.5194/adgeo-14-211-2008)CrossRefGoogle Scholar
Popovnin, VV and Rozova, AV (2002) Influence of sub-debris thawing on ablation and runoff of the Djankuat glacier in the Caucasus. Nordic Hydrol., 33(1), 7594 (doi: 10.2166/nh.2002.005)Google Scholar
Price, LW (1981) Mountains and man: a study of process and environment. University of California Press, Berkeley Google Scholar
Raup, BH and 5 others (2007) The GLIMS geospatial glacier database: a new tool for studying glacier change. Glob. Planet. Change, 56, 101110 (doi: 10.1016/j.gloplacha.2006.07.018)Google Scholar
Schroeder, LD, Sjoquist, DL and Stephan, PE (1986) Understanding regression analysis. Sage Publications, New York Google Scholar
Shangguan, D and 9 others (2006) Monitoring the glacier changes in the Muztag Ata and Konggur mountains, east Pamirs, based on Chinese Glacier Inventory and recent satellite imagery. Ann. Glaciol., 43(1), 7985 (doi: 10.3189/172756406781812393)Google Scholar
Shahgedanova, M, Stokes, CR, Gurney, SD and Popovnin, V (2005) Interactions between mass balance, atmospheric circulation, and recent climate change on the Djankuat Glacier, Caucasus Mountains, Russia. J. Geophys. Res. – Atmos., 110, 4964 (doi: 10.1029/2004JD005213)Google Scholar
Solomina, ON (2000) Retreat of mountain glaciers of northern Eurasia since the little ice age maximum. Ann. Glaciol., 31, 2630 Google Scholar
Stokes, CR, Gurney, SD, Shahgedanova, M and Popovnin, V (2006) Late-20th-century changes in glacier extent in the Caucasus Mountains, Russia/Georgia. J. Glaciol., 52(176), 99109 (doi: 10.3189/172756506781828827)Google Scholar
Sugden, DE and John, BS (1969) Glaciers and landscape. Edward Arnold, London Google Scholar
White, H (1980) A heteroscedasticity-consistent covariance matrix estimator and a direct test for heteroscedasticity. Econometrica, 48, 817838 (doi: 10.2307/1912934)Google Scholar
Zemp, M, Frauenfelder, R, Haeberli, W and Hoelzle, M (2005) Worldwide glacier mass balance measurements: general trends and first results of the extraordinary year 2003 in Central Europe. In: XIII Glaciological Symposium, Shrinkage of the Glaciosphere: Facts and Analyses, St. Petersburg, Russia [Science, R.A.O. (ed.)]. Data of Glaciological Studies [Materialy glyatsiologicheskikh issledovaniy], Moscow, Russia, 312 Google Scholar
Zolotarev, YA and Kharkovets, YG (2010) Evolution of the Elbrus glaciation since the mid XIX century under changing climate. Key findings of the glacio-cartographical monitoring. Geogr. Environ. Sustainability, 2, 1332 Google Scholar
Figure 0

Fig. 1. Mount Elbrus glacier system outlines, superposed over the 2007 colour composite of Landsat bands 5, 4, 3.

Figure 1

Table 1. Change of the glaciers surface area, relative losses and rate of change (1985–2007)

Figure 2

Table 2. Change of glacier altitude (1985–2007)

Figure 3

Fig. 2. Automatic extraction and surface analysis algorithm.

Figure 4

Table 3. Correlation coefficient values (R) between glacier losses and relief parameters

Figure 5

Table 4. Unstandardized (B) and standardized (β) coefficients of the multiple regression

Figure 6

Fig. 3. Retreat trends in Elbrus Mountain since 1850 (% a−1). Data source: Zolotarev and Kharkovets (2010); p. 24, Table 3