CC BY
31
DOI: https://doi.org/10.23670/IRJ.2019.90.12.027
ОЖИДАЕМОЕ ИЗМЕНЕНИЕ ПОВЕРХНОСТНОГО БАЛАНСА МАССЫ ЛЕДНИКОВОГО КОМПЛЕКСА ЭЛЬБРУСА В УСЛОВИЯХ ГЛОБАЛЬНОГО ПОТЕПЛЕНИЯ
Научная статья
Рыбак О.О.1' *, Рыбак Е.А.2, Корнева И.А.3
1 ORCID: 0000-0003-3923-7163;
2 ORCID: 0000-0003-3920-4176;
3 ORCID: 0000-0002-6453-8315;
1 2 Сочинский научно-исследовательский центр РАН, Сочи, Россия;
1 2 3 Филиал Института природно-технических систем, Сочи, Россия;
1 Институт водных проблем РАН, Москва, Россия;
3 Институт глобального климата и экологии им. академика Ю.А. Израэля, Москва, Россия
* Корреспондирующий автор (o.o.rybak[at]gmail.com)
Аннотация
Площадь ледникового комплекса Эльбруса, который состоит из 24 ледников, составляет в настоящее время около 10% площади оледенения всего Кавказа. Изменения климата оказывали значительное влияние на параметры оледенения Эльбруса в течение последнего столетия. Учитывая ощутимый вклад талой воды ледников Эльбруса в общий речной сток в горных и предгорных регионах Северного Кавказа и его стабилизирующую роль, становится понятным интерес исследователей к прогностическим расчетам поверхностного баланса массы и ледникового стока. В настоящей статье мы исследуем указанные параметры оледенения при реализации экстремального климатического сценария RCP8.5 к концу текущего столетия. Показано, в частности, что одновременное увеличение приземной температуры воздуха, суммарной солнечной радиации и сокращение осадков, рассчитанные на модели глобального климата и приведенные к региону Центрального Кавказа, будут причиной существенного отступания снеговой линии и сокращения площади аккумуляции, что будет причиной усиливающейся деградации оледенения Эльбруса в грядущие годы.
Ключевые слова: горный ледник, Кавказ, Эльбрус, поверхностный баланс массы, ледниковый сток, математическая модель, численный эксперимент, прогноз.
EXPECTED CHANGE IN SURFACE BALANCE OF MASS IN GLACIUM COMPLEX OF ELBRUS UNDER
CONDITIONS OF GLOBAL WARMING
Research article
Rybak O.O.1, *, Rybak E.A.2, Korneva I.A.3
1 ORCID: 0000-0003-3923-7163;
2 ORCID: 0000-0003-3920-4176;
3 ORCID: 0000-0002-6453-8315;
1 2 Sochi Research Center, Russian Academy of Sciences, Sochi, Russia;
1 2 3 Branch of the Institute of Natural Technical Systems, Sochi, Russia; institute for Water Problems of the Russian Academy of Sciences, Moscow, Russia 3.Institute for Global Climate and Ecology named after Yu. A. Israel, Moscow, Russia
* Corresponding author (o.o.rybak[at]gmail.com)
Abstract
The area of the Elbrus glacier complex, which consists of 24 glaciers, currently makes up about 10% of the glaciation area of the entire Caucasus. Climate change has a significant impact on the parameters of the glaciers of Elbrus over the past century. Given the tangible contribution of meltwater from Elbrus glaciers to the total river flow in the mountain and foothill regions of the North Caucasus and its stabilizing role, the interest of researchers in the prognostic calculations of the surface mass balance and glacial runoff becomes clear. In this paper, the authors examine the indicated glaciation parameters during the implementation of the extreme climatic scenario RCP8.5 by the end of this century. It was shown, in particular, that a simultaneous increase in surface air temperature, total solar radiation and a decrease in precipitation calculated on the global climate model and reduced to the Central Caucasus region will cause a significant retreat of the snow line and a reduction in the accumulation area, which will cause the increasing degradation of the Elbrus glaciation in coming years.
Keywords: mountain glacier, Caucasus, Elbrus, the surface balance of mass, glacial runoff, mathematical model, numerical experiment, forecast.
Введение
Фактическая сторона масштабов сокращения оледенения Эльбруса достаточно хорошо известна и отражена как в работах отечественных [1], [2], [3], так и зарубежных исследователей [4]. Эволюция оледенения Эльбруса в целом соответствует тенденции уменьшения площади оледенения на Кавказе в целом в последние два десятилетия 20-го века [5] и в начале 21-го века [6]. За период 1960-2014 площадь оледенения Эльбруса сократилась на 14,7±2,4%, темпы сокращения составили в среднем 0,27% в год [7], а за два десятилетия, 1997-2017 гг., - на 10,8% [8]. По состоянию на 2017 г. площадь оледенения Эльбруса составляла 112,2±0,58 км2 [8] (по другой оценке 115,4 км2, а за вычетом скальных выступов и морен - 108,8 км2 в 2015 г. [9]), его объем оценивается в 5,05±0,85 км3 [8]. Реки, берущие начало на склонах Эльбруса, относятся к бассейнам Черного и Каспийского морей. Доля ледникового стока в
общем объеме стока рек северного склона Большого Кавказа оценивается в 1-9%, однако в теплый период года на отдельных реках его доля может достигать 40-45% [10]. Очевидно, что оледенение Эльбруса вносит значительный вклад в формирование режима стока, поскольку служит своего рода стабилизатором, поставляя воду в речную сеть в теплый период года, когда потребность в ней наиболее велика. Последнее обстоятельство обуславливает необходимость аккуратного расчета ледникового стока, который формируется на ледниковом комплексе Эльбруса, и, следовательно, расчет поверхностного баланса массы. Наиболее удобным инструментом для этого будет математическое моделирование, тем более, если стоит задача прогнозирования изменений стока с большой заблаговременностью.
Для настройки математической модели и ее валидации требуется определенный набор метеорологической и климатической информации в регионе исследования.
Начиная с 30-х годов XX века и до недавнего времени на склонах Эльбруса выполнялись фрагментарные метеонаблюдения различной продолжительности [11]. Инструментальные наблюдения за компонентами баланса массы систематически проводились, начиная с 1983 г., только на Гарабаши - одном из ледников Эльбруса на его южном склоне [12] (рис. 1). Для доинструментального периода ряд баланса массы был реконструирован с помощью дендрохронологических методов [13]. Таким образом, прямую калибровку и валидацию модели можно провести только на ограниченном участке ледникового комплекса, на котором проводились инструментальные наблюдения. В качестве альтернативы, для косвенной валидации модели возможно использовать космические снимки ледникового комплекса в разные моменты времени, однако это довольно трудоемкая задача, являющаяся, вообще говоря, самостоятельным научным исследованием.
Рис. 1 - Цифровая модель рельефа ледникового комплекса Эльбруса (показан серым цветом) и окружающей области
Примечание: светло-серым выделен ледник Гарабаши. Зеленые квадратики указывают расположение автоматических метеостанций на леднике Гарабаши (3850 м над у.м.) и Западном плато (5150 м над у.м.). Здесь и на рис. 3 и 4 боковая оцифоровка показывает расстояние в метрах от условной угловой точки области; изолинии
высоты поверхности проведены через 200 м
Математическая модель поверхностного баланса массы горного ледника подробно разобрана ранее, поэтому в настоящей статье мы не будем останавливаться на ее описании. Упомянем лишь, что модель в качестве составной части (суб-модели) комплексной динамической модели применялась для исследования горных ледников Марух (Западный Кавказ) [14], Джанкуат (Центральный Кавказ) [15], [16], Сары-Тор (Внутренний Тянь-Шань) [17], Туюксу (Заилийский Алатау, не опубл.). В приведенных работах содержится подробное описание модели.
В настоящей статье рассматривается вероятные изменения поверхностного баланса массы ледникового комплекса Эльбруса в результате реализации экстремального по версии Межправительственной группы экспертов по изменениям климата (МГЭИК) сценария ЯСР8.5 [18].
Особенности постановки численных экспериментов
Климатический форсинг
Как и в работах [15], [16] в качестве исходной климатической информации использованы ряды 1999-2008 гг. приземной температуры воздуха на метеостанции Терскол, расположенной в полутора десятках километров от южных отрогов Эльбруса, и сумм осадков на метеостанции Местиа, расположенной несколько дальше, на территории Грузии на южном склоне Главного Кавказского хребта. Выбор станции Местиа обусловлен тем, что ряды осадков на ней имеют гораздо меньше пропусков, чем на станции Терскол, а общее количество и распределение осадков в течение
года очень схожи с таковыми на Терсколе. Исходные ряды среднесуточной температуры воздуха и сумм осадков осреднялись за балансовое десятилетие, в результате были получены два ряда длиной из 365 среднесуточных (температура) и суточных сумм (осадки) значений.
Значения температуры воздуха в каждом узле регулярно сетки 25x25 м на леднике были получены масштабированием значений на Западном плато Эльбруса (5150 м над у.м.) с использованием высотного градиента -0,6 °С/км, что соответствует среднему температурному градиенту для этого региона [11]. В свою очередь, значения температуры воздуха на Западном плато были получены из ряда среднесуточных значений на станции Терскол с помощью простого линейного соотношения ТЗападное плато =а+ЬТТерскол, рассчитанному в ходе анализа параллельных наблюдений на двух метеостанциях - Терскол и автоматической метеостанции (АМС), которая функционировала в августе и сентябре 2007 г. на Западном плато Эльбруса на абсолютной высоте 5150 м над у.м. (рис. 1 и 2а). Характерно, что значение коэффициента b практически совпадает с аналогичным значением в линейном уравнении связи рядов температуры воздуха на станции Терскол и АМС на леднике Джанкуат, расположенном в 20 км от Терскола в противоположном (по отношению к Эльбрусу) направлении [16]. В то же время, значение коэффициента b в выражении ТГарабаши =а+ЬТТерскол, связывающем температуру воздуха на АМС на леднике Гарабаши (3850 м над у.м., функционировала в июле и августе 2013 г., см. рис. 1), и на метеостанции Терскол, несколько больше предыдущего (рис. 2б).
-2
-4
о
I- пЗ -6
с;
с:
<и -8
О
X
rci -10
с
го -12
р
н -14
-16
• •
f ••• • . • :
f • • W • у* • . / • • i ••• ■i
f • -• / 'г •
• • • • I
• .......... • ......... • Т(Зап.плато) = а+Ь~1 а = -16.49 ±0.77 b = 0.81 ±0.07 ............................. (Терскол): : .............г
^ 4
3
га ю
fü о О. ^
m
О
8 10 12 14 16 Т, °С (Терскол)
-2
-4
: ~'"(Гарабаши)-а+'ЭТ(Терскол) ra = -10.45 ±1.40 : b = 1.05 ±0.12 • • • •• 1 t
г • • • t/4 • ____ / • ■ / • -
« •• • • > •
Г • • / / • • •
/* ;.....• • ......... ......... ......... 1 's
8 10 12 Т, °С (Терскол)
14
Рис. 2 - Линейные соотношения между среднесуточной приземной температурой воздуха на сетевой метеостанции Терскол и АМС на Западном плато Эльбруса (а) и АМС на леднике Гарабаши (б)
Высотный градиент осадков в диапазоне высот 3300-4000 м составил 1,3х 10-4 м/м в холодное полугодие (октябрь-март) и 0,5х 10-4 м/м в теплое полугодие (апрель-сентябрь). Такой выбор значений высотного градиента обеспечил максимальное значение осадков и аккумуляции в поясе 3900-4000 м (1327 и 1275 мм водн. экв. соответственно), что в целом соответствует наблюдениям на леднике Гарабаши [9,] [12]. Опять же согласно наблюдениям [12] выше 4000 м аккумуляция (но не количество осадков) быстро снижается из-за перемещения снега вниз по склону лавинами и выдувания ветром. Учитывая значительную межгодовую изменчивость аккумуляции, в модели были заданы такие высотные градиенты осадков выше 4000 м, какие бы обеспечили среднегодовую аккумуляцию выше 5500 м на уровне 445 мм водн. экв., что качественно соответствует обрывочным данным наблюдений. Набор остальных параметров модели обеспечил итоговый модельный баланс массы всего комплекса в контрольном эксперименте очень близкий к нулевому (-3,1 мм водн. экв.). Фактически, этот баланс массы соответствует климатическим условиям десятилетия 1999-2008 гг.
2.2 Формализация климатических изменений
Климатические изменения в регионе исследования получены в результате расчетов на климатической модели INMCM Института вычислительной математики РАН и последующей регионализации полученных данных к исследуемой области. Методика регионализации и полученные результаты для Черноморско-Каспийского региона целом и для отдельных регионов на Кавказе детально разобраны в работе [19]. Прогноз внутригодового роста температуры и суммарной радиации в среднем за период 2071-2100 относительно средних 1971-2000 согласно [19] приведен в табл. 1. Данные расчетов были использованы нами для оценки вероятных изменений поверхностного баланса массы (SMB) и ледникового стока (RO), а также отношения площади аккумуляции к общей площади ледникового комплекса (AAR) вследствие роста приземной температуры воздуха и суммарной радиации (эксп. 2 и 3 в табл. 2) Изменение осадков к концу текущего столетия в Черноморско-Каспийском регионе согласно [19] в целом статистически незначимо. Однако в отдельных областях в пределах региона прогнозируется их снижение, причем как в зимние, так и в летние месяцы (например, в области ледников Марух и Джанкуат). Поэтому нами проведены также и эксперименты, в которых задано сокращение осадков на 10% в холодное и теплое полугодие и в течение всего года (эксп. 4-6 в табл. 2). В эксп. 7 задействованы изменения всех трех факторов.
Таблица 1 - Модельные оценки изменения приземной температуры воздуха для Центрального Кавказа (ДГ, °С) и _суммарной радиации (Д^, Вт м-1) в 2071-2100 гг. по сравнению с периодом 1971-2000 гг. [19]_
Мес. I II III IV V VI VII VIII IX X XI XII
AT 3,1 2,9 2,8 2,8 3,5 4,6 5,2 4,6 3,8 3,0 3,5 3,5
AS -0,3 0 0 1,2 18,3 21,0 15,8 2,2 10,7 9,9 2,0 -0,5
Таблица 2 - Характеристика численных экспериментов и их результаты.
Эксп. 1 2 3 4 5 6 7
AT - + + - - - +
AS - - + - - - +
APW - - - + - - -
APs - - - - + - -
APws - - - - - + +
AAR 0,73 0,22 0,19 0,70 0,72 0,68 0,14
SMB -3 -854 -960 -96 -56 -148 -1103
RO 1072 1757 1868 1095 1051 1073 1863
Обозначения: знак «+» означает, что данный вид форсинга присутствует, знак «-» - отсутствует, AT -приземная температура, AS" - суммарная солнечная радиация, aAPw - зимние осадки, APS - летние осадки, aAPws - годовые осадки, AAR - соотношение между площадью аккумуляции и общей площадью ледника (accumulation area ratio), SMB - поверхностный баланс массы (мм водн. экв.), RO - сток (мм водн. экв.).
Эксперимент 1 - контрольный
Обсуждение результатов
В контрольном эксперименте (эксп. 1) абляция фактически прекращается выше отметки 5100 м над у.м. (рис. 3а), ниже 4300 м сезон абляции продолжается более трех месяцев. Заметим, что в контексте настоящего исследования за день с абляцией принимается день, в который таяние наблюдалось хотя бы в течение одного часа. z Годовой баланс массы, как было упомянуто выше, близок к нулю, AAR=0,73 (то есть в зоне абляции, где расход массы больше, чем ее приход, находится немногим более четверти площади оледенения). Высота снеговой линии (границы между зонами аккумуляции и абляции) находится на высоте около 3600 м над у.м., иногда опускаясь на высоту до 3500 м над у.м. (рис. 4а). В зоне с абляцией более 1000 мм водн. экв. находятся только концевые части ледников.
О 2000 4000 6000 8000 10000 12000 14000 16000 18000 0 2000 4000 6000 8000 10000 12000 14000 16000 18000
О 2000 4000 6000 8000 10000 12000 14000 16000 18000 0 2000 4000 6000 8000 10000 12000 14000 16000 18000
Рис. 3 - Количество дней с абляцией в контрольном эксперименте (Эксп. 1 в табл. 2) (а)
и в Эксп. 7 (б)
Подключение температурного форсинга в эксп. 2 (табл. 1) приводит к резкому снижению баланса массы до -854 мм водн. экв. и подъему снеговой линии до высоты 4400-4600 м над у.м. (рис. 4б). Соответственно более, чем в три раза сокращается площадь зоны аккумуляции.
Рис. 4 - Поверхностный баланс массы (мм водн. экв.) в контрольном эксперименте (Эксп. 1) (а), эксп. 2 (б), эксп. 3 (в), эксп. 4 (г), эксп. 5 (д), эксп. 6 (е), эксп. 7 (ж) - см. табл. 2. Коричневым цветом показана снеговая линия
При подключении радиационного форсинга (эксп. 3) происходит дальнейшее снижение баланса массы (до -960 мм водн. экв.). Снеговая линия поднимается еще на 200 м.
Использование в качестве форсинга снижение сумм осадков на ожидаемую характерную величину 10% - зимних (эксп. 4), летних (эксп. 5) и годовых (эксп. 6) влечет за собой не столь разрушительные для ледника последствия, как ожидаемый рост приземной температуры воздуха. AAR остается без существенных изменений, баланс массы, разумеется, несколько понижается за счет сокращения его приходной части баланса массы. Характерно, что в эксп. 5 снижение летних осадков вызывает снижение ледникового стока, который состоит не только из талой ледниковой воды, но и из суммы выпавших жидких осадков, которые не подверглись вторичному замерзанию при просачивании в толщу снега/фирна.
При воздействии всех трех факторов (эксп. 7) изменения наиболее внушительны. Снеговая линия поднимается до высоты 4800-4900 м над у.м., а AAR сокращается в пять раз по сравнению с контрольным экспериментом, так что
зона аккумуляции составляет лишь 14% общей площади комплекса. Величина ледникового стока увеличивается на 80% по сравнению с контрольным экспериментом.
Здесь, однако, необходимо сделать очень важное примечание. В настоящей работе все оценки сделаны для конфигурации ледника по состоянию на 2014 г. Нами не рассматривалось течение льда, которое приводит к перераспределению массы льда вниз по склону. Соответственно, все оценки сделаны, исходя из современной площади оледенения. Разумеется, по мере отступания ледников Эльбруса AAR, баланс массы и сток необходимо будет рассчитывать, исходя из меньшей площади. Тем не менее, прогностическое положение снеговой линии, а также ожидаемое количество дней с абляцией (за исключением показателя для участков, которые должны исчезнуть) изменится несущественно.
Благодарности
Авторы благодарят П.А. Торопова (географический факультет МГУ им. М.В. Ломоносова) за предоставленные метеорологические данные, И.И. Лаврентьева и С.С. Кутузова (Институт географии РАН) за данные по рельефу исследуемого региона. Конфликт интересов
Не указан
Acknowledgements
Authors appreciate P.A. Toropov (Department of Geography, M.V. Lomonosov Moscow State University) for meteorological data, I.I. Lavrentiev and S.S. Kutuzov (Institute of Geography of RAS) for data on the relief in the region of interest.
Conflict of interest
Not declared
Список литературы / References
1. Золотарев Е.А. Сокращение оледенения Эльбруса в ХХ столетии / Е.А. Золотарев, А.А. Алейников, Е.Г. Харьковец // Материалы гляциологических исследований. - 2005. - Вып. 98. - C. 162-166.
2. Золотарев Е.А. Оценка деградации оледенения Эльбруса методами цифрового картографирования / Е.А. Золотарев, Е.Г. Харьковец // Вестник МГУ, Серия географическая. - 2007. - №5. - С. 45-51.
3. Золотарев Е.А. Картографо-аэрокосмический мониторинг баланса массы оледенения Эльбруса после МГГ / Е.А. Золотарев, Е.Г. Харьковец // Материалы гляциологических исследований. - 2007. - Вып. 103. - С. 123-128.
4. Holobaca I.-H. Recent retreat of the Elbrus glacier system / I.-H. Holobaca // Journal of Glaciology. - 2016. - V. 62(231). - P. 94-102. doi: 10.1017/jog.2016.15
5. Stokes C.R. Late-20th-century changes in glacier extent in the Caucasus Mountains, Russia/Georgia / C.R. Stokes, S.D. Gurney, M. Shahgedanova, V. Popovnin // Journal of Glaciology. - 2006. - V. 52(176). - P. 99-109.
6. Shahgedanova M. Deglaciation of the Caucasus Mountains, Russia/Georgia, in the 21st century observed with ASTER satellite imagery and aerial photography / M. Shahgedanova, G. Nosenko, S. Kutuzov et. al. // The Cryosphere. - 2014. - V. 8.
- P. 2367-2379. doi:10.5194/tc-8-2367-2014
7. Tielidze L.G. The Greater Caucasus Glacier Inventory (Russia, Georgia and Azerbaijan) / L.G. Tielidze, R.D. Wheate // The Cryosphere. - 2018. - V. 12. -P. 81-94. https://doi.org/10.5194/tc-12-81-2018
8. Kutusov S. Volume Changes of Elbrus Glaciers From 1997 to 2017 / S. Kutusov, I. Lavrentiev, A. Smirnov et. al. // Frontiers in Earth Sciences [Электронный ресурс]. URL: https://doi.org/10.5194/tc-12-81-2018. 2018. V. 7. Article 153.
9. Калов Х.М. Комплексный мониторинг динамики ледников Эльбруса / Х.М. Калов, М.Ю. Беккиев, М.Д. Докукин // Материалы Всероссийской конференции «Устойчивое развитие: проблемы, концепции, модели». -Нальчик: 2017. - С. 156-159.
10. Лурье П.М. Влияние изменения климата на современное оледенение и сток рек северного склона Большого Кавказа / П.М. Лурье, В.Д. Панов // Устойчивое развитие горных территорий. - 2013. - Т. 5. - №2. - С. 70-77.
11. Торопов П.А. Температурный и радиационный режим ледников на склонах Эльбруса в период абляции за последние 65 лет / П.А. Торопов В.Н. Михаленко, С.С. Кутузов и др. // Лед и Снег. - 2016. - Т. 56. - №1. - С. 5-19. doi:10/15356/2076-6734-2016-1-5-19
12. Рототаева О.В. Изменения баланса массы ледника Гарабаши (Эльбрус) на рубеже XX-XXI вв. / О.В. Рототаева, Г.А. Носенко, А.М. Керимов и др. // Лед и Снег. - 2019. - Т. 59. - №1. - С. 5-22. doi: 10.15356/2076-67342019-1-5-22
13. Долгова Е.А. Реконструкция баланса массы ледника Гарабаши (1800-2005 гг.) по дендрохронологическим данным / Е.А. Долгова В.В. Мацковский, О.Н. Соломина и др. // Лед и Снег. - 2013. - Т. 53. - №1. - С. 34-42.
14. Рыбак О.О. Калибровка математической модели динамики ледника Марух, Западный Кавказ / О.О. Рыбак, Е.А. Рыбак, С.С. Кутузов и др. // Лед и ^ег. - 2015. - Т. 55. - №2. - С. 9-20.
15. Rybak O.O. Model-based calculations of surface mass balance of mountain glaciers for the purpose of water consumption planning: focus on Djankuat Glacier (Central Caucasus) / O.O. Rybak., E.A Rybak. // IOP Conference Series: Earth and Environmental Science. - 2018. - V. 107. - doi :10.1088/1755- 1315/107/1/012041.
16. Рыбак О.О. Применение данных сетевых метеорологических станций для расчета баланса массы ледников (на примере ледника Джанкуат, Центральный Кавказ) / О.О. Рыбак, Е.А. Рыбак // Системы контроля окружающей среды.
- 2017. - Вып. 9 (29). - C. 100-108.
17. Рыбак О.О. Модельные исследования эволюции горных ледников на примере ледника Сары-Тор (Внутренний Тянь-Шань) / О.О. Рыбак, Е.А. Рыбак., Н.А. Яицкая и др. // Криосфера Земли. - 2019. - Т. 23. - №3. - С. 33-51. doi:10.21782/KZ1560-7496-2019-3(33-51)
18. Climate Change 2013: The Physical Science Basis. Contribution of Working Group I to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change (AR5) / Еds .: T.F. Stocker, D. Qin, G. K. Plattner et al . Cambridge, United Kingdom and New York, NY, USA: Cambridge University Press. - 2013. 1535 p .
19. Морозова П.А. Регионализация данных глобального климатического моделирования для расчета баланса массы горных ледников / П.А. Морозова, О.О. Рыбак // Лед и Снег. 2017. Т. 57. №4. С. 437-452. doi: 10.15356/20766734-2017-4-437-452
Список литературы на английском языке / References in English
1. Zolotarev E.A. Sokrashchenie oledeneniya El'brusa v XX stoletii [Shrinking of Elbrus glaciation in the XX century] / E.A. Zolotarev, A.A. Alejnikov, E.G. Har'kovec // Materialy glyaciologicheskih issledovanij [Data of Glaciological Studies]. -2005. - No. 98. - P. 162-166 [in Russian]
2. Zolotarev E.A. Ocenka degradacii oledeneniya El'brusa metodami cifrovogo kartografirovaniya [Estimate of degradation of Elbrus glaciation using methods of numerical mapping] / E.A. Zolotarev, E.G. Har'kovec // Vestnik MGU, Seriya geograficheskaya [Herald of the Moscow State University, Geographical Series]. - 2007. - №5. - P. 45-51. [in Russian]
3. Zolotarev E.A. Kartografo-aerokosmicheskij monitoring balansa massy oledeneniya El'brusa posle MGG [Mapping and aerial-and-space monitoring of Elbrus glaciation mass balance after IGY] / E.A. Zolotarev, E.G. Har'kovec // Materialy glyaciologicheskih issledovanij [Data of Glaciological Studies]. - 2007. - No. 103. - P. 123-128. [in Russian]
4. Holobaca I.-H. Recent retreat of the Elbrus glacier system / I.-H. Holobaca // Journal of Glaciology. - 2016. - V. 62(231). - P. 94-102. doi: 10.1017/jog.2016.15
5. Stokes C.R. Late-20th-century changes in glacier extent in the Caucasus Mountains, Russia/Georgia / C.R. Stokes, S.D. Gurney, M. Shahgedanova, V. Popovnin // Journal of Glaciology. - 2006. - V. 52(176). - P. 99-109.
6. Shahgedanova M. Deglaciation of the Caucasus Mountains, Russia/Georgia, in the 21st century observed with ASTER satellite imagery and aerial photography / M. Shahgedanova, G. Nosenko, S. Kutuzov et. al. // The Cryosphere. - 2014. - V. 8. - P. 2367-2379. doi:10.5194/tc-8-2367-2014
7. Tielidze L.G. The Greater Caucasus Glacier Inventory (Russia, Georgia and Azerbaijan) / L.G. Tielidze, R.D. Wheate // The Cryosphere. - 2018. - V. 12. -P. 81-94. https://doi.org/10.5194/tc-12-81-2018
8. Kutusov S. Volume Changes of Elbrus Glaciers From 1997 to 2017 / S. Kutusov, I. Lavrentiev, A. Smirnov et. al. // Frontiers in Earth Sciences [Электронный ресурс]. URL: https://doi.org/10.5194/tc-12-81-2018. 2018. V. 7. Article 153.
9. Kalov H.M. Kompleksnyj monitoring dinamiki lednikov El'brusa [Complex monitoring of Elbrus glaciers dynamics] / H.M. Kalov, M.Yu. Bekkiev, M.D. Dokukin // Materialy Vserossijskoj konferencii «Ustojchivoe razvitie: problemy, koncepcii, modeli» [Materials of the All-Russian Conference "Sustainable development: problems, concepts, models"]. -Nal'chik: 2017. - P. 156-159. [in Russian]
10. Lur'e P.M. Vliyanie izmeneniya klimata na sovremennoe oledenenie i stok rek severnogo sklona Bol'shogo Kavkaza [Impact of climate change on the modern glaciation and river runoff on the northern slope of the Greater Caucasus] / P.M. Lur'e, V.D. Panov // Ustojchivoe razvitie gornyh territory [Sustainable development of mountain territories]. - 2013. - V. 5. -№2. - P. 70-77. [in Russian]
11. Toropov P.A. Temperaturnyj i radiacionnyj rezhim lednikov na sklonah El'brusa v period ablyacii za poslednie 65 let [Temperature and radiation regime of glaciers on slopes of the Mount Elbrus in the ablation period over last 65 years] / P.A. Toropov V.N. Mihalenko, S.S. Kutuzov et al. // Led i Sneg [Ice and Snow]. - 2016. - V. 56. - №1. - P. 5-19.doi:10/15356/2076-6734-2016-1-5-19 [in Russian]
12. Rototaeva O.V. Izmeneniya balansa massy lednika Garabashi (El'brus) na rubezhe XX-XXI vv. [Changes of the mass balance of the Garabashy Glacier, Mount Elbrus, at the turn of 20th and 21st centuries] / O.V. Rototaeva, G.A. Nosenko, A.M. Kerimov et al. // Led i Sneg [Ice and Snow]. - 2019. - V. 59. - №1. - P. 5-22. doi: 10.15356/2076-6734-2019-1-5-22 [in Russian]
13. Dolgova E.A. Rekonstrukciya balansa massy lednika Garabashi (1800-2005 gg.) po dendrohronologicheskim dannym [Reconstruction of mass balance of Garabashi Glacier (1800-2005) using dendrochronology data] / E.A. Dolgova V.V. Mackovskij, O.N. Solomina et al. // Led i Sneg [Ice and Snow]. - 2013. - V. 53. - No. 1. - P. 34-42. [in Russian]
14. Rybak O.O. Kalibrovka matematicheskoj modeli dinamiki lednika Maruh, Zapadnyj Kavkaz [Calibration of a mathematical model of Marukh Glacier, Western Caucasus ] / O.O. Rybak, E.A. Rybak, S.S. Kutuzov et al. // Led i sneg [Ice and Sow]. - 2015. - V. 55. - №2. - P. 9-20. [in Russian]
15. Rybak O.O. Model-based calculations of surface mass balance of mountain glaciers for the purpose of water consumption planning: focus on Djankuat Glacier (Central Caucasus) / O.O. Rybak., E.A Rybak. // IOP Conference Series: Earth and Environmental Science. - 2018. - V. 107. - doi :10.1088/1755- 1315/107/1/012041.
16. Rybak O.O. Primenenie dannyh setevyh meteorologicheskih stancij dlya rascheta balansa massy lednikov (na primere lednika Dzhankuat, Central'nyj Kavkaz) [Application of observations on the regular meteorological stations for calculation of mass balance of mountain glaciers (case study of Djankuat Glacier, Central Caucasus)] / O.O. Rybak, E.A. Rybak // Sistemy kontrolya okruzhayushchej sredy [Systems of Environmental Control]. - 2017. - No. 9 (29). - P. 100-108. [in Russian]
17. Rybak O.O. Model'nye issledovaniya evolyucii gornyh lednikov na primere lednika Sary-Tor (Vnutrennij Tyan'-SHan') [Modeling the evolution of mountain glaciers: a case study of Sary-Tor glacier, Inner Tien Shan] / O.O. Rybak, E.A. Rybak., N.A. YAickaya et al. // Kriosfera Zemli [Earth's Cryosphere]. - 2019. - V. 23. - No 3. - P. 33-51. doi:10.21782/KZ1560-7496-2019-3(33-51) [in Russian]
18. Climate Change 2013: The Physical Science Basis. Contribution of Working Group I to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change (AR5) / Eds .: T.F. Stocker, D. Qin, G. K. Plattner et al . Cambridge, United Kingdom and New York, NY, USA: Cambridge University Press. - 2013. 1535 p.
19. Morozova P.A. Regionalizaciya dannyh global'nogo klimaticheskogo modelirovaniya dlya rascheta balansa massy gornyh lednikov [Downscaling of the global climate model data for the mass balance calculation of mountain glaciers] / P.A. Morozova, O.O. Rybak // Led i Sneg [Ice and Snow]. 2017. V. 57. No. 4. P. 437-452. doi: 10.15356/2076-6734-2017-4-437452 [in Russian]
