Реконструкция гидрометеорологических условий на северном кавказе по дендрохронологическим данным за период с
На правах рукописи
Долгова Екатерина Антоновна Реконструкция гидрометеорологических условий на Северном Кавказе по дендрохронологическим данным за период с 1800-2005 гг.
Специальность 25.00.25 – Геоморфология и эволюционная география
Автореферат диссертации на соискание ученой степени кандидата географических наук
Москва – 2011 г.
Работа выполнена в Учреждении Российской академии наук Институте географии РАН
Научный консультант:
чл.-корр. РАН, доктор географических наук Ольга Николаевна Соломина
Официальные оппоненты:
доктор географических наук, ведущий научный сотрудник Борисова Ольга Кимовна кандидат географических наук, доцент Виктор Владимирович Поповнин
Ведущая организация:
Учреждение Российской академии наук Институт проблем экологии и эволюции им. А.Н. Северцова РАН
Защита состоится «28 » октября 2011 г. в на заседании Диссертационного совета Д 002.046. в Учреждении Российской академии наук Институте географии РАН по адресу: 119017 Москва, Старомонетный переулок,
С диссертацией можно ознакомиться в библиотеке Учреждения Российской академии наук Института географии РАН Текст объявления и автореферат размещены на сайте Института географии РАН www.igras.ru Автореферат разослан « 28 » сентября 2011 г.
Отзывы на автореферат (в двух экземплярах), заверенные печатью, просим отправлять по адресу: 119017 Москва, Старомонетный переулок, 29.
Факс: (495) 959-00-33;
E-mail: [email protected], [email protected]
Ученый секретарь диссертационного совета, кандидат географических наук И.С. Зайцева Актуальность темы. Годичные кольца деревьев являются важным источником косвенной информации о климате прошлого.
Дендрохронологический метод позволяет реконструировать температуру, осадки, речной сток с годовым или сезонным разрешением, и таким образом, удлинить ряд наблюдений и использовать эти данные для прогноза, оценить частоту и суровость экстремальных явлений, тренды и ритмику изменений климатических параметров и связанных с ними природных процессов.
Высокогорные районы дают уникальную возможность для определения и анализа изменчивости климата в прошлом по годичным кольцам деревьев. Это связано с тем, что деревья, растущие на пределе своего существования, обладают повышенной чувствительностью к изменению тех климатических параметров, которые ограничивают их рост (лимитирующие факторы). Именно в горах встречаются высоковозрастные деревья, которые обеспечивают построение реконструкций продолжительностью в несколько столетий, а иногда даже и тысячелетий. Несмотря на относительно высокую изученность Кавказского региона в отношении современных изменений климата, климатические реконструкции высокого разрешения, основанные на результатах дендрохронологического метода, здесь до сих пор отстутствуют.
Абсолютно сдатированных древесно-кольцевых хронологий для этой территории почти нет. Поэтому создание таких хронологий, исследование зависимостей прироста деревьев от климатических колебаний и использование этих данных для реконструкции климатической изменчивости на Кавказе в прошлом, представляют интерес, как с научной, так и с практической точки зрения. Этим проблемам посвящена настоящая работа.
Цель работы – создание абсолютно сдатированных хронологий по ширине и максимальной плотности колец сосны и пихты и выявление их возможностей для реконструкции температуры тёплого периода, баланса массы ледников и речного стока для территории Северного Кавказа.
Задачи исследования. Для достижения поставленной цели в диссертации решались следующие взаимосвязанные задачи:
1. Создание сети абсолютно сдатированных древесно-кольцевых хронологий по ширине и плотности годичных колец сосны и пихты в долинах рек Теберда и Баксан;
пространственно-временной анализ этих дендрохронологических рядов.
2. Построение сводных древесно-кольцевых хронологий, отражающих региональные особенности прироста.
3. Анализ связи метеорологических параметров с древесно-кольцевыми хронологиями.
4. Создание количественной реконструкции температуры воздуха теплого периода для территории центральных районов северного макросклона Кавказа и оценка ее достоверности.
5. Анализ связей индексов древесно-кольцевых хронологий со значениями баланса массы ледников Приэльбрусья. Оценка возможности применения хронологий для реконструкции баланса массы ледников.
6. Исследование возможности создания реконструкции стока р. Теберды по дендрохронологическим данным.
Научная новизна работы заключается в том, что в ней впервые представлены абсолютно сдатированные хронологии сосны и пихты по ширине годичных колец и максимальной плотности древесины в долинах рек Баксан и Теберда продолжительностью 205-456 лет. В результате этого исследования была создана первая количественная реконструкция летних температур для Кавказа с годичным разрешением за период 1800-2005 гг. Впервые показана возможность применения дендрохронологического метода для рекнструкции баланса массы ледника Гарабаши и стока р. Теберды.
Защищаемые положения • На Кавказе на верхней границе леса на радиальный прирост сосны текущего года оказывают влияние температуры мая, июня и ноября предыдущего года.
На формирование ширины годичного кольца пихты отрицательно влияют температуры августа и сентября предыдущего года. Устойчивой статистически значимой связи (пригодной для реконструкции) между индексами ширины годичных колец сосны и пихты со среднемесячными температурами и количеством осадков текущего года не обнаружено.
• Максимальная плотность древесины сосны вблизи верхней границы леса на Кавказе более тесно связана с климатической изменчивостью, чем ширина годичных колец. Самая высокая статистически значимая связь, которая объясняет более 50% изменчивости, обнаружена между индексами хронологии по максимальной плотности колец и среднемесячной температурой апреля-сентября высокогорной метеостанции Северный Клухор (2047 м). Такая связь статистически значима и для ряда наблюдений на гмс Пятигорск продолжительностью 110 лет, что свидетельствует о её устойчивости во времени и подтверждает достоверность полученной реконструкции.
• Согласно выполненной реконструкции, за последние два столетия самое холодное лето за 200 лет отмечалось в 1817 году, возможно, в связи с крупным извержением вулкана Тамбора в 1815 году, самое теплое – 1872 г.
Частота экстремумов за 200 лет не обнаруживает статистически значимого группирования во времени. Периоды летних похолоданий на Кавказе отмечались в 1832-1838, 1895-1906, 1926-1936 гг. Линейных трендов в ходе реконструированных температур за последние 50, 100, 150, 200 лет не обнаружено. Тренды температуры апреля-сентября отсутствуют и в ряду инструментальных наблюдений г. Пятигорска за последние 110 лет.
• Ширина и максимальная плотность годичных колец сосны коррелируют с речным стоком р.Теберды и балансом массы ледника Гарабаши.
Работа основана на обширном Личный вклад автора.
дендрохронологическом материале, полученным автором в ходе полевых работ в 2006-2007 гг., а также сотрудниками ИГРАН с 2002 по 2009 гг. Автор проводил полевые работы по сбору дендрохронологического материала, лабораторную обработку образцов, измерение ширины годичных колец, химическую обработку образцов и измерение максимальной плотности колец, перекрёстную датировку полученных древесно-кольцевых серий и построение хронологий, статистическую обработку метеорологической и дендрохронологической информации, построение реконструкций гидрометеорологических параметров и оценку их точности и надежности, сопоставление построенной реконструкции температуры с реконструкциями других авторов по Северному полушарию.
Практическая значимость работы. Полученные данные о температурном режиме тёплого периода на Кавказе позволяют уточнить характер и особенности климатической изменчивости этого района в прошлом, что обеспечивает более надежный прогноз, причем как длительных изменений, так и короткопериодных флуктуаций. Этот прогноз, в свою очередь, позволит обеспечить более обоснованное долгосрочное планирование хозяйственной и рекреационной деятельности на территории Северного Кавказа. Практический интерес представляют и полученные данные об изменчивости в прошлом стока реки Теберды, притока Кубани, питающей район с сельскохозяйственной специализацией. Построенные хронологии могут быть в дальнейшем использованы для датировки погребённой древесины и исторических памятников. Информация о климатическом отклике деревьев может использоваться в лесоводстве при прогнозе реакции радиального прироста деревьев на изменения климата.
Апробация работы. Основные положения диссертации были доложены на российских и международных конференциях: симпозиум «Гляциология в канун Международного Полярного Года», Пушкинские Горы (октябрь 2006), международная конференция «Вклад России в МПГ 2007/2008» в г. Сочи (октябрь 2009), конференция «Лёд и снег в климатической системе» в г. Казань (май 2010), на международной конференции «WorldDendro 2010» в г.Рованиеми (Финляндия, июнь 2010), тезисы представлены на конференции РусДендро, Екатеринбург (октябрь 2011).
Структура и объем работы. Диссертационная работа изложена на страницах машинописного текста и состоит из введения, пяти глав, заключения и двух приложений. Список литературы включает 120 наименования, в том числе 51 - на иностранных языках. Текст иллюстрирован 70 рисунками и содержит 25 таблиц.
Глава 1. Район работ и состояние его изученности в отношении дендроклиматологии.
Работы проводились на Северном Кавказе на территории Тебердинского Государственного Биосферного заповедника и в национальном парке Приэльбрусье (рис. 1).
Рис. 1. Расположение пробных площадей отбора дендрохронологического материала (жёлтые кружки), используемых метеостанций (зелёные кружки), некоторых ледников (голубые квадраты). Площадки отбора в Тебердинском заповеднике: КYZ – долина реки Кызгыч, KHTP и KHAT – долина р. Тебрда, ALI – предполье ледника Алибек. В Приэльбрусье:
BAZ – конечная морена ледника Большой Азау, GAR – предполье ледника Гарабаши, CHE – верхняя граница леса склона г. Чегет восточной экспозиции г. Чегет, CHS - верхняя граница леса склона г. Чегет южной экспозиции, KV – морена 14 в. ледника Большой Азау, TERS – верхняя граница леса близ ледника Терскол.
Тебердинский Государственный Биосферный заповедник, представленный двумя участками Тебердинским и Архызским, располагается на северном макросклоне Главного Кавказского хребта. На юге заповедник ограничен отрезком главного хребта между вершиной Джаловчат (3870 м) и Клухор-баши (3450 м), на западе граница проходит по водораздельному гребню между реками Теберда и Малая Марка, восточная – охватывает водораздельную цепь между реками Теберда и Даут. Река Теберда образуется при слиянии рек Аманаус и Гоначхир и является левым притоком реки Кубани. Приэльбрусье вместе с Эльбрусом (5642 м) – самая западная часть Центрального Кавказа (Природные процессы..., 2004). Большой Кавказ входит в состав Альпийского (Средиземноморского) геосинклинального пояса (Гвоздецкий и др., 1987).
Современный рельеф изучаемой территории создавался под влиянием длительного воздействия разнообразных факторов, главнейшими из которых являются геологическое строение, влияние древнего оледенения, водная и лавинная эрозия (Тушинский, 1957). Существенной особенностью орографии является большая разница относительных высот, что, оказывает влияние на формирование природных особенностей района, а также обуславливает вертикальную поясность ландшафтов.
Климат района определяется его положением в умеренном широтном климатическом поясе, с чем связаны особенности радиационного режима и циркуляции атмосферы. Зимой на процессы атмосферной циркуляции на Северном Кавказе влияют Исландский минимум и Азиатский антициклон.
Средние температуры воздуха и суммы выпавших осадков за холодный период (ноябрь-март) составляют:
-2,9 0С и 733 мм, -0,4 0С и 271 мм, -4,7 0С и 294 мм для гмс Северный Клухор, Теберда и Терскол, соответсвенно. В тёплый период года на циркуляционные процессы на Северном Кавказе большое влияние оказывает Азорский максимум. На эти макроциркуляционные процессы накладывается влияние Чёрного и Каспийского морей. Средние значения температуры воздуха и суммы выпавших осадков за тёплый период (апрель октябрь) для гмс Северный Клухор, Теберда и Терскол составляют: +8,6 0С и 900,5 мм, +11,5 0С и 504,7 мм, +7,9 0С и 648,9 соответственно.
На высоте 1300 м (гмс Теберда) переход температур через ноль наблюдается в конце февраля - в марте. На гмс Северный Клухор и Терскол, расположенных выше 2000 м, температуры становятся положительными в середине марта.
Ход средних летних температур на гмс Северный Клухор не обнаруживает статистически значимого тренда за период 1957-2005 гг. Для гмс Теберда характерно увеличение среднегодовых температур (r=0,33, p0,05) за счёт их роста в летние месяцы (r=0,35;
0,48;
0,44;
p0,05) и в сентябре (r=0,38;
p0,05).
Статистически значимый отрицательный тренд обнаружен для ряда среднегодовых температур на гмс Терскол (r=-0,35;
p0,05) за период 1951 2005 гг. Из трёх метеостанций лишь для гмс Терскол обнаружен статистически значимый положительный тренд в ходе осадков за год, а также в июне (r=0,33;
p0,05) и сентябре (r=0,35;
p0,05).
В исследуемом районе ледники относятся к бассейнам двух морей:
Азовского и Каспийского. В бассейне реки Кубань, в истоках которой находится Тебердинский заповедник, площадь оледенения достигает 23,9% от общей площади оледенения Северного Кавказа, в бассейне реки Терек - 69,9% (Панов, 2008). Ледникам и снежникам принадлежит большая роль в питании рек, берущих начало из ледников. Так, в Тебердинском заповеднике, где проводились работы по реконструкции речного стока, вклад ледникового питания составляет 26% от общего объема речного стока.
Основные дендрохронологические работы по отбору образцов велись на пределе распространения сомкнутого высокоствольного леса. Условия произрастания растений в пределах экотона верхней границы леса характеризуются сочетанием низких температур и значительного увлажнения, короткого вегетационного периода при интенсивной солнечной радиации;
почвы имеют малую мощность и низкую плодородность. На Северном Кавказе с увеличением степени континентальности климата с запада на восток возрастает количество солнечных дней и, следовательно, интенсивность радиации, что способствует продвижению вверх лесной растительности, поэтому высота верхней границы леса на Эльбрусе выше, чем в Тебердинском заповеднике. На Западном Кавказе верхнюю границу леса составляют разнотравно-злаковые буковые криволесья, сменяющиеся берёзовыми криволесьями. В более континентальных условиях Центрального Кавказа верхнюю границу леса на высоте 2400 м образуют сосновые леса с травяным ярусом с доминированием осоки низкой, овсяницы овечьей и других видов сухих лугов (Биота экосистем Большого Кавказа, 1990).
Дендрохронологические работы в высокогорье Центрального Кавказа начались еще в 1970-х гг. В.И.Турманина (1971, 1972, 1979) установила, что прирост сосны в Приэльбрусье снижается в годы существенного понижения температур вегетационного периода и/или в периоды снижения нормы годовых осадков при повышении температур вегетационного периода. Обширные дендрохронологические работы на Кавказе в 1970-е-1980-е гг. велись Институтом Ботаники АН Литовской ССР. Результаты этих работ (ширина колец и ее индексированные значения) публиковались в сборнике «Дендрошкалы Советского Союза». Позже П.В. Ковалёвым и др. (1984) была предпринята попытка использования одного спила пихты для реконструкции среднегодовых температур для гмс Теберда на Западном Кавказе. В последние годы на кафедре ландшафтоведения МГУ дендрохронологический метод используется для индикации динамики изменений ландшафтов в высокогорьях Кавказа (Бочкарёв, Дьяконов, 2009). Одна хронология для тиса ягодного (Taxus baccata L.) за период с 1556-1980 гг. представлена для территории Грузии (N=42, E=45,17;
H=2000 м) в Международном банке дендрохронологических данных http://www.ncdc.noaa.gov/paleo/treering.html (Kuniholm, Steele, Tovar).
Несмотря на интересные результаты, полученные исследователями в конце 20-го века, детальные работы по выявлению сопряженности динамики приростов в пространстве, по идентификации климатического сигнала в ширине годичных колец у разных древесных пород, по реконструкции колебаний климата и других параметров среды все еще составляют задачи будущего. Некоторые аспекты этих проблем будут рассмотрены ниже в этой работе.
Глава 2. Материалы и методы Дендрохронологический материал.
С 2002 года вблизи верхней границы леса и на предпольях ледников на высоте 1900-2500 м сотрудниками отдела гляциологии ИГРАН было отобрано около 200 кернов сосны обыкновенной (Pinus sylvestris L.) и пихты кавказской (Abies nordmaniana (Steven) Spach.);
они послужили основой для этой работы (Долгова и др., 2007) (рис. 1). Все дендрохронологические площадки, использованные в диссертации, расположены в природоохранной зоне, что обеспечивает минимальное антропогенное влияние на местообинания деревьев.
В работе использовались данные метеостанций Теберда, Северный Клухор, Терскол и Пятигорск (табл. 1). Отбор образцов был организован вблизи высокогорных метеостанций для более успешного поиска климатического сигнала. Данные гмс Пятигорск были использованы при анализе устойчивости модели.
Таблица 1. Использованные в работе метеорологические и гидрологические данные Координаты Метеостанция Высота, м Длина ряда Метеопараметры (с.ш., в.д.) Теберда 4327, 4144. 1313 1956- Северный Клухор 4315., 4150 2047 1956-2005 Среднемесячные температуры и осадки Пятигорск 4405, 43 05 538 1891- Терскол 4315, 4230 2214 1951- Гидрологический Среднемесячный сток 4326, 4144 2210 1927- пост Теберда р. Теберда Методика исследования. Отбор и анализ образцов древесины (керны и спилы) проводился в соответствии с общепринятыми методическими требованиями древесно-кольцевого анализа (Шиятов и др., 2000;
Cook and Kairiukstis, 1990). Измерения ширины годичных колец произведены на полуавтоматической установке LINTAB с точностью 0,01 мм. Далее данные по ширине годичных колец подвергались перекрёстной датировке в программе TSAP-Win с последующей проверкой в программе COFECHA (Holmes, 1983).
Оценка возрастного тренда производилась для каждого индивидуального дерева, возрастная кривая аппроксимировалась негативной экспонентой или прямой с отрицательным наклоном, относительно которой рассчитывались индексы прироста. Для целей стандартизации и удаления возрастного тренда использовалась программа ARSTAN (Cook, 1991). Для сохранения долгопериодной изменчивости в рядах в последние годы многими исследователями (Esper et al., 2002;
Esper, 2003;
Frank et al., 2005) возрастной тренд удаляется при помощи так называемой стандартизации региональной кривой (Regional Curve Standardization). В данной работе использованы стандартные хронологии, а метод RCS не был применён в виду отсутствия субфоссильной древесины с датами гибели, отличными от современных живых деревьев.
Известно, что максимальная плотность поздней древесины является более интегральным показателем температуры периода вегетации, чем ширина колец (Schweingruber, 1988;
D’Arrigo, 1992). С одной стороны, максимальная плотность поздней древесины отражает температурные условия того времени, когда она формируется (конец лета), а с другой – находится под гормональным контролем условий начала вегетационного периода (апрель-май) (Hughes, 2001;
Vaganov et al., 2006). На основе максимальной плотности колец восстановлены температуры для Северной Америки (Davi et al., 2003;
D’Arrigo et al., 2004;
Luckman and Wilson, 2005), Сибири (Briffa et al. 2001;
Kirdyanov et al., 2008), Европы (Briffa et al., 1988, 2004;
Bntgen et al., 2008).
Анализ плотности древесины, методические основы которых впервые предложены Польжем (Polge, 1966), является сложным, трудоёмким и дорогостоящим процессом (Ваганов, 2000). В связи с этим возникло новое альтернативное направление – измерение оптической плотности древесины с использованием сканеров высокого разрешения и специального программного обеспечения для анализа (Sheppard et al., 1996;
Campbell et al., 2007). Профили максимальной плотности сосны для территории северной Финляндии, полученные традиционным и альтернативным способами, согласуются между собой и практически одинаково отражают изменчивость температуры (Campbell et al., 2007). Эти теоретические и технологические разработки послужили основанием для применения альтернативного метода измерения плотности древесины («оптическая плотность») в данной работе.
Древесина сосны содержит значительное количество смол (Shweingruber, 1988), поэтому для получения корректных результатов керны были предварительно обработаны в аппарате Сокслета в растворах этанола и толуола (Jensen, 2007), так как наличие смолы в древесине может серьёзно повлиять на результаты измерений плотности. Затем поверхность кернов была зачищена с помощью микротома. Отсканированные керны в высоком разрешении (1000 dpi) обрабатывались в програме Lignovision (Rinntech), которая использует изменчивость яркости колец для построения профиля плотности. Удаление смолы из кернов и измерения плотности производились в лаборатории древесно-кольцевых исследований, Аризона, США в 2009 году.
Для оценки качества древесно-кольцевых рядов использовались такие показатели как коэффициент чувствительности и стандартное отклонение рядов. Качество хронологии оценивается с помощью разных показателей, в том числе EPS (Expressed Population Signal (Briffa et al., 1995). Значения EPS в нашем случае были посчитаны для каждой хронологии с помощью 30-летнего «плавающего» окна с 29-летним перекрытием. Значение индекса EPS 0, гарантирует, что ширина кольца за данный год адекватно отражает прирост всей совокупности (Wigely et al., 1984). В данной работе пороговое значение индекса было принято 0,8, что позволило удлинить ряд с незначительной потерей качества.
Для определения статистических связей между различными климатическими параметрами и индексами древесно-кольцевых хронологий использовался корреляционный анализ. Рассчитывались коэффициенты корреляции Пирсона дендрохронологических данных с ежемесячными или ежегодными параметрами за период январь–октябрь текущего года прироста и май-декабрь предыдущего года (Fritts, 1976). После того, как была установлена связь метеорологическх показателей с древесно-кольцевыми хронологиями, строилась регрессионная модель, позволяющая перейти от показателей ширины и плотности кольца к климатическим параметрам. В нашем случае это была модель линейной регрессии.
При создании статистической модели, описывающей связь дендрохронологических и метеорологических показателей, калибровочный ряд обычно делят на две части, собственно калибровочную и верификационную, для того, чтобы проверить устойчивость модели. В нашем случае, учитывая короткий калибровочный ряд (25 лет), для оценки устойчивости и надежности нашей модели мы дополнительно использовали ряд гмс Пятигорск 1891-1997 гг.
Короткий ряд инструментальных наблюдений за балансом массы ледника Гарабаши не позволил использовать стандартную методику деления на калибровочную и верификационную части. В данном случае для оценки качества модели был использован метод скользящего контроля: модель строилась по (n-1) точкам, а в оставшейся точке рассчитывалось отклонение от прогноза, и описанная процедура повторялась n раз (т.е. для каждой точки) (Вапник, 1974).
Глава 3. Влияние климатических факторов на ширину и плотность годичных колец и реконструкции летних температур Хронологии ширины и плотности годичных колец. Всего в ходе этого исследования построено 7 хронологий по ширине годичных колец сосны, 1 по ширине годичных колец пихты и 1 по максимальной плотности колец сосны.
Длительность хронологий составляет от 264 до 456 лет для сосны и 205 лет для пихты (рис. 2, табл.2). Статистические характеристики хронологий показывают, что ограничений для их использования в целях климатических реконструкций нет. Хронология по максимальной плотности колец получена для деревьев с площадок урочища Хатипара (KHTP и KHAT), которые впоследствии были объеденены в хронологию MaxD. Корреляционный анализ, а также применение метода главных компонент показали, что первый фактор 6 хронологий по ширине годичных колец сосны объясняет почти 50% дисперсии (рис. 3).
Второй фактор объясняет 13% общей изменчивости и может быть выделен только для хронологии KYZ. Очевидно, что условия произрастания сосен в долине реки Кызгыч отличаются от остальных и являются уникальными.
Сводная хронология сосны, таким образом, не включает хронологию KYZ и построена путем осреднения 6 локальных хронологий. Региональная хронология покрывает период с 1708 года (при наличии 5 образцов в ее начале) по 2009 год.
Влияние климатических факторов на прирост сосны и пихты. Для хронологий по ширине годичных колец сосны и пихты Тебердинского района, а также для Эльбруских хронологий, расположенных на верхней границе леса, обнаружена положительная статистически значимая связь индексов прироста и температуры ноября предыдущего года (r = 0,36-0,42;
p0,05). На прирост пихты отрицательно сказываются высокие температуры августа и сентября предыдущего года (r = -0,33;
-0,43;
p0,05). Величина прироста сосны вблизи верхней границы леса не обнаруживает статистически значимых связей с температурой лета и количеством осадков текущего года. Для некоторых хронологий по ширине колец (BAZ, KV и KHTP) обнаружена слабая положительная связь с осадками июня и июля. Очевидно, что в самые жаркие месяцы лета отсутствие атмосферных осадков негативно сказывается на приросте деревьев. На формирование максимальной плотности древесины сосны положительно влияют высокие среднемесячные температуры тёплого периода года марта – сентября (рис. 4).
1550 1600 1650 1700 1750 1800 1850 1900 1950 образцов Индексы Кол-во KYZ 0 Кол-во образцов KHTP Индексы 1 0 1, Кол-во образцов CHS Индексы 1, 0,5 Кол-во образцов Индексы GAR 0 2, Кол-во образцов 1, Индексы KV 1, 0, 0,0 Кол-во образцов Индексы BAZ 0 Кол-во образцов Индексы TERS 0 Кол-во образцов Индексы ALI 0 образцов Индексы Кол-во MaxD 1550 1600 1650 1700 1750 1800 1850 1900 1950 Рис.2. Индексы хронологий по ширине (KYZ, KHTP, CHS, GAR, KV, BAZ, TERS) и плотности (MaxD) годичных колец сосны, по ширине годичных колец пихты (ALI) и количество используемых образцов.
Таблица 2. Статистические характеристики хронологий Название хронологии KYZ KV KHTP GAR CHS BAZ ALI TERS MaxD Длина хронологии 1550- 1640- 1678- 1693- 1738- 1660- 1800- 1714- 1759 2006 2002 2005 2002 2002 2004 2005 2009 Кол-во образцов 20 8 57 13 15 18 24 15 (деревьев), вошедших (14) (4) (29) (7) (8) (11) (13) (9) (43) в хронологию Межсериальный 0,4 0,55 0,51 0,61 0,53 0,44 0,28 0,58 0. коэффициент корреляции Стандартное 0,23 0,18 0,22 0,22 0,15 0,14 0,13 0,16 0. отклонение Средняя 0,16 0,13 0,16 0,13 0,11 0,12 0,12 0,12 0. чувствительность Средняя 0,66 0,55 0,51 0,61 0,52 0,44 0,28 0,58 0. автокорреляция 1-ого порядка Средняя длина 203 283 175 160 160 233 114 233 сегмента EPS0,8 (год) 1783 1805 1753 1892 1961 1752 1883 1776 Доля первой главной 47,1 46,0 42,0 39,7 44,2 42,6 38 32,0 компоненты, % Рис. 3. Поля распределения корреляционных связей между локальными хронологиями по ширине годичных колец в проекции двух главных компонент за общий период 1738-2002.
Рис. 4. Корреляция индексов максимальной плотности сосны и среднемесячной температуры воздуха гмс Северный Клухор.
Реконструкция температур вегетационного периода. Для реконструкции была выбрана наиболее устойчивая, статистически значимая и экологически обусловленная зависимость между индексами хронологии сосны по максимальной плотности древесины и среднемесячными температурами апреля-сентября гмс Северный Клухор (R2 =0,63;
p0.05). При сравнении восстановленных и инструментальных данных можно отметить, что практически все пики обоих параметров совпадают во времени и по амплитуде (рис. 5). Для независимой проверки модели инструментальный период наблюдений был разделён на 2 части: верификации (1956-1980 гг.) и калибровки (1981-2005 гг.). Все полученные коэффициенты корреляции для этих отдельных периодов достигают статистически значимых значений (табл.3).
Таблица 3. Результаты верификации модели Период Rsquare RE CE 1956-2005 0,5196 - 1956-1980 0,7044 0,2481 0, 1981-2005 0,3702 0,5747 0, Из реконструкции температуры апреля-сентября, представленной в виде аномалий относительно климатической нормы 1961-1990 гг., следует, что самые холодные вегетационные периоды за последние 200 лет наблюдались в 1817 (-2,6 С), 1810 (-2,1С), 1884 (-2,1С) гг. (в скобках - аномалии температур апреля-сентября, рассчитанные относительно 1961-1990 гг.). Периоды положительных аномалий отмечались в 1820-1831, 1839-1894, 1907-1925, 1937 1944, 1948-1955 гг., отрицательных - в 1832-1838, 1895-1906, 1926-1936, 1956 1962 гг. (рис. 5). На рис. 5 также заметно отсутствие долгопериодных трендов за последние 200 лет и за более короткие вековые и полувековые интервалы времени.
Рис 5. Аномалии температур апреля-сентября, рассчитанные относительно климатической нормы 1961-1990 гг., полученные по максимальной плотности годичных колец сосны (синяя линия), по инструментальным данным гмс Северный Клухор (красная линия), и по инструментальным данным гмс Пятигорск (зелёная линия). Ряды восстановленных и инструментальных температур гмс Пятигорск сглажены 7-летним скользящим средним. Корреляция посчитана между несглаженными рядами среднемесячных температур апреля-сентября по инструментальным данным гмс Пятигорск и восстановленных по дендрохронологии.
Гмс Пятигорск (1891-1997 гг.) расположена в предгорьях Кавказа на расстоянии около 150 км от площадок отбора дендрохронологических образцов. Возможность использования температурных данных этой станции обусловлена их высокой согласованностью с температурами высокогорных метеостанций. Сравнение температуры апреля-сентября гмс Пятигорск с температурой, реконструированной по дендрохронологическим данным (рис.
5), показывает высокую согласованность рядов (r = 0,63;
p0,05), в том числе для периода 1891-1960 гг. (r = 0,59;
p0,05), который не был использован для построения модели. Самая низкая корреляция между рядами, которая обнаружена для самых ранних периодов наблюдений в начале 20 века (1905 1934 гг., r = 0,42;
p0,05), все же остаётся на статистически значимом уровне.
Это сравнение с более длинным рядом наблюдений подтверждает адекватность построенной модели и устойчивость сигнала во времени, а также означает, что хронология по максимальной плотности колец сосны, полученная для верховьев р. Теберды, отражает климатический сигнал не только на локальном уровне, но значима и для всего центрального региона Северного Кавказа. Это также подтверждается сравнением реконструкции с данными реанализа (рис. 6), которое показывает положительную корреляцию этой реконструкции с температурой всего причерноморского региона. Особенно высока эта корреляция для западной и центральной части Кавказа. Интересно, что корреляция отмечается и для южного макросклона Кавказа, и даже для Закавказья. Область положительных корреляций выделяется также в Малой Азии, в восточном Средиземноморье и в Северной Африке, откуда на Кавказ часто приходят воздушные массы летом.
Рис. 6. Пространственная корреляция реконструированной температуры апреля сентября с данными реанализа NCEP/NCAR доступными с 1948-2005 гг.
Хронологии были проанализированы на предмет цикличности динамики прироста с помощью вэйвлет-анализа и спектрального (Фурье) анализа. Во всех хронологиях (по ширине и плотности) выделяется ритм в 24-27 лет. В сводной хронологии по ширине годичных колец сосны для Кавказа также обнаружены циклы в 60 и 10 лет. Циклы более высокого порядка достоверно не идентифицируются в виду относительно небольшой длины рядов.
Глава 4. Подходы к реконструкции речного стока и баланса массы ледников на Кавказе. Прирост древесины может быть использован для реконструкции речного стока, так как в засушливых районах оба показателя находятся под влиянием сходного набора климатических факторов (например, осадки и суммарное испарение) (Fritts, 1976;
Meko et al., 1995). В связи с этим древесно-кольцевые хронологии успешно использовались для реконструкции стока рек в аридных районах Северной Америки (например, Woodhouse et al., 2004;
Diaz et al., 2001), в горах Тянь-Шаня (Yuan et al., 2007), в Монголии (Davi et al., 2006). По-видимому, некоторый дефицит увлажнения наблюдается и в районе нашего исследования, так как ширина колец сосны коррелирует здесь с такими интегральными показателями засушливости, как объем речного стока в летние месяцы.
Ширина годичных колец сосны положительно коррелирует с величиной речного стока р. Теберды за отдельные месяцы, хотя теснота этой связи и невелика. До 60 % стока р. Теберды приходится на лето, 17 % на осень, 5 % на зиму и 18 % на весну. 55,8 % стока, измеренного на гидрологической станции Теберда, приходится на таяние снега и льда. Среднегодовой расход воды здесь равен 26,8 м3/с. Ежемесячный сток реки Теберды коррелирует с температурой весеннего сезона (r до 0,62 в апреле) и с осадками осеннего сезона (r до 0,72 в октябре). Значимые отрицательные тренды присутствуют и в среднем стоке за летние месяцы и в среднегодовом стоке за период с 1927 по 2005 г.
Для создания моделей реконструкции стока р. Теберды мы использовали метод регрессии на главные компоненты. Модели были сгенерированы для каждого из 12 месяцев и для всего гидрологического года. Семь стандартных древесно-кольцевых хронологий использовались в качестве предикторов. Для определения наилучшего набора предикторов для каждой из моделей, были испробованы все возможные их комбинации (полный перебор) и выбран набор, показавший максимальный коэффициент детерминации R2 на скользящем контроле. Реконструкции годичного разрешения показали недостаточно хорошие результаты на скользящем контроле. Поэтому мы сгладили реконструированные ряды 10-летним скользящим средним. В результате обнаружены связи: майский сток (лучший набор предикторов – хронологии ALI, BAZ, CHS и KHAT, R2 = 0,6487), июльский сток (ALI, BAZ и KHAT, R2 = 0,6567) и августовский сток (ALI, CHS и KHAT, R2 = 0,6615) (рис. 7, a).
Согласно инструментальным данным, на эти месяцы приходится 12 %, 22 % и 17 % годового стока, соответственно. Был проведен анализ вейвлет-когеренции между инструментальными и реконструированными (не сглаженными) данными для оценки возможности реконструкции воспроизводить низкочастотную изменчивость. Майская и июльская реконструкции показали сильную и стабильную когеренцию с инструментальными данными на периодах колебаний более 25 лет, в отличие от августовской реконструкции (рис. 7, б).
Нарушение декадной изменчивости майского и июльского стока заметно при визуальном анализе сглаженных реконструкций (рис. 7, a). Вейвлет-анализ может определить точный год этих изменений. Непрерывное вейвлет преобразование реконструкций стока (не сглаженных) (см. рис. 7, г) показывает ослабление долгопериодных циклов (периоды от 16 до 32 лет) после 1929 г. для июля и 1938 г. для мая. Мы не проводили вейвлет-анализ августовской реконструкции из-за слабой вейвлет-когеренции реконструкции и инструментальных данных. Сток всех трех реконструированных месяцев (май, июль и август), также как и годовой сток, имеет значимый отрицательный тренд в период с 1927 по 2005 г. Значения трендов равны –0,063, –0,140, –0, и –0,429 м3/с/год, соответственно.
Рис. 7. Реконструкции стока реки Tеберды и их спектральные характеристики. а значения стока по (1) инструментальным данным, (2) реконструкция по данным дендрохронологии;
б – вейвлет-когеренция реконструированного и инструментального стока за 1927–2002 гг., черные границы обозначают области, значимые с вероятностью более 95 % на фоне красного шума;
скользящие тренды реконструированного стока (движущееся окно в 76 лет) за 1850–2002 гг., год соответствует первому году в окне, пунктирными линиями обозначены границы значимости на уровне 95 %.
Возможности реконстуркции баланса массы ледников (на примере ледника Гарабаши).
Данные баланса массы и колебаний концов ледников могут быть использованы в качестве источников информации об изменчивости климата в прошлом. Прирост деревьев, растущих вблизи верхней границы леса на Кавказе, определяется набором параметров сходным с тем, что и для баланса массы ледников. В первую очередь это - летние температуры и осадки, в том числе и твердые. Это сходство было замечено ранее и во многих других горных странах (Альпы, Скалистые Горы, Скандинавия), где дендрохронологический метод успешно применялся для реконструкции баланса массы ледников (Linderholm et al., 2007;
Luckman et al., 1993;
Watson et al., 2003;
Nicolussi et al., 2000;
Wood et al., 2011).
В данной работе предпринята попытка использования древесно-кольцевых хронологий сосны по ширине и максимальной плотности колец для восстановления баланса масс ледников Джанкуат и Гарабаши. Для территории Северного Кавказа подобная работа выполнена впервые.
Наиболее детальные балансовые наблюдения на Кавказе проводятся на леднике Джанкуат с 1967 года по настоящее время. Джанкуат – типичный долинный ледник северо-западной экспозиции площадью около 3 км2, расположенный в верховьях долины Адылсу в Приэльбрусье (см. рис.1).
Хронология сосны по максимальной плотности коррелирует с балансом массы ледника Джанкуат (r=-0,44;
p0,05). Однако, использование одной хронологии по максимальной плотности для реконструкции баланса массы ледника невозможно не только из-за слабости корреляционной связи, но и потому что в этом случае в балансе окажется выраженным лишь сигнал, связанный с летней температурой, от которой зависит плотность древесины. Применив модель множественной регрессии, мы обнаружили, что самый высокий коэффициент корреляции (r=0,6;
p0,05) достигается при наборе всех 8 хронологий. Однако на скользящем контроле значения корреляции становятся незначимыми, что делает эту модель недостоверной.
Более удачными были эксперименты с балансом массы другого ледника – Гарабаши. Ледник Гарабаши расположен на южном склоне Эльбруса в непосредственной близости от Эльбрусских дендрохронологических площадок.
Это – часть ледниковой шапки Эльбруса южной экспозиции. Его площадь составляет 4,47 км ( Рототаева и др., 2000). За весь период инструментальных наблюдений (1983-2005 гг.) (www.geo.uzh.ch/microsite/wgms/), включающий года, значения аккумуляции и абляции ледника оказались почти равными ( и 1330 мм в.э. соответственно). Среднеквадратическое отклонение величин аккумуляции 174,7 мм в.э., относительная изменчивость Cv = 14,3%, для абляции, соответственно, 426,6 мм в.э. и 32,1%. Таким образом, летнее таяние на леднике характеризуется большей изменчивостью, чем аккумуляция. Летние температуры вносят больший вклад в изменчивость баланса массы ледника Гарабаши, по сравнению зимними осадками. Анализ изменчивости аккумуляции, абляции и баланса массы не выявил статистически значимых линейных трендов.
Корреляционный анализ показывает более тесную, чем для ледника Джанкуат, связь баланса массы ледника Гарабаши и хронологии по максимальной плотности (r=-0,62;
p0,05). В этом случае также наиболее адекватной для целей реконструкции баланса оказывается модель множественной регрессии, которая включает несколько предикторов, ответственных за разные климатические параметры. Поиск наилучшей модели (с наивысшими коэффициентами корреляциями) производился путём перебора всех возможных наборов предикторов (индексы хронологий) для каждого отклика (аккумуляция, абляция, баланс массы). В модель вошли 2 хронологии по ширине годичных колец сосны, а именно KYZ и KHTP, и хронология по максимальной плотности. Отметим, что параметры годичных колец сосны (плотность и ширина) с высокой точностью отражают амплитуду и пики хода баланса массы ледника. Особенно это заметно при самом низком значении баланса массы за инструментальный период в 1998 году. В то же время отмечаются некоторые различия. Так, например, в 2000 году значение баланса массы недооценено по дендрохронологическим данным, т.е. реальные значения несколько ниже, чем реконструированные. Отметим, однако, что даже в случае 2000 года, для которого заметно существенное различие модели и данных измерений, реконструированный баланс массы остается отрицательным (рис.
8).
Рис. 8. Сравнение хода инструментальных наблюдений за балансом массы ледника Гарабаши (красная линия) с балансом массы, реконструированным по дендрохронологии до (голубая линия) и после (голубая пунктирная линия) скользящего контроля.
Обнаруженная статистическая связь теоретически позволяет реконструировать баланс массы ледника Гарабаши до 1800 г. (рис. 9). Для того чтобы оценить полученную реконструкцию, мы сравнили её с реконструкцией О.В. Рототаевой и др. (2000) за 1905-2000 гг. Для продления ряда абляции авторы использовали метеостанцию Пятигорск с 1891 г. Реконструкция аккумуляции на леднике опиралась на косвенные расчёты осадков по данным метеостанции Сухими-Хаиши для 1904/05-1937/38 гг. и Местия - 1938/39 1950/51 гг. Данные прямых наблюдений, использованные для калибровки модели, охватывают 1983/84-1994/95 гг.
Реконстуркции баланса массы ледника Гарабаши, основанные на дендрохронологических и метеорологических данных, в целом имеют высокую взаимную погодичную согласованность (r = 0,53;
p0,05), особенно для периода 1975-2000 гг. Однако при сравнении двух реконструкций обнаружены и существенные различия. Так, в 1920-е гг., согласно реконструкции О.В.Рототаевой с коллегами, значения баланса массы были близки к нулю, в то время как в реконструкции по дендрохронологическим данным в этот период преобладали отрицательные значения баланса массы. С 1929 г. обе балансовые кривые идут вниз вплоть до 1950 гг. Это согласуется с ходом летней температуры по данным гмс Терскол, которая в это время была в целом высокой. В 1960-е гг. произошла резкая смена циркуляционных условий, которая благоприятно сказалась на состоянии ледников Кавказа. Однако кривые реконструкций существенно различаются в это время. Так ряд, полученный по дендрохронологии, имеет отрицательные значения, а реконструкция Рототаевой и др. - положительные. Интересно, что дендрохронологическая реконструкция в этот период времени согласуется с реконструкцией баланса массы Джанкуат (Дюргеров и др. 1981). После 1975г.
ряды по Гарабаши имеют высокую корреляцию и значимых различий не обнаружено. Несогласованность хода двух реконструкций в отдельные периоды может быть связана с рядом причин. Как одну из наиболее вероятных можно назвать влияние на баланс лет с аномально высокой аккумуляцией, которая не всегда находит отражение в дендрохронологических рядах. К недостаткам реконструкции можно также отнести короткий инструментальный ряд наблюдений, на котором строится модель, и возможная нестационарность сигнала во времени (Leonelli et al., 2011).
На рис. 9 видно, что у двух рассмотренных здесь реконструкций баланса массы ледников Грабаши и Джанкуат, основанных на метеорологических данных, также есть существенные расхождения, и коэффициент их взаимной корреляции составляет всего r = 0,53;
p0,05. Это связано с различной методикой их построения. Интересно, что дендрохронологическая реконструкция в период 1880-1920 гг. больше похожа на баланс массы ледника Джанкуат и достаточно точно воспроизводит его внтуривековую изменчивость (рис. 9).
Несмотря на некоторые различия при сравнении с гляцио метеорологическими реконструкциями, можно утверждать, что в целом полученная дендрохронологическая реконструкция адекватно отражает тенденции изменений ледников на Кавказе. Так, она показывает накопление массы ледников в 1830-х- 1860-х гг., когда граница питания на ледниках Кавказа была примерно на 150 м ниже, чем в конце 20-го века, а ледники были существенно длиннее и имели большую толщину (Серебрянный и др., 1988;
Золотарёв, 2009).
Рис. 9. Сравнение реконструкций баланса массы ледника Гарабаши, полученные автором по дендрохронологии (зелёная линия), Рототаевой и др. (2000) построенные путем продления рядов метеорологических наблюдений (голубая линия), Дюргеровым и др. (1981) для ледника Джанкуат. Ряды нормированы и сглажены 7-летним скользящим средним.
Глава 5. Сравнение дендрохронологических данных на Кавказе с другими данными об изменениях климата.
Известно, что эксплозивные вулканические извержения вызвают кратковременные, но сильные похолодания летом, что часто отражается в структуре, ширине и плотности колец либо в тот же год, либо с запаздыванием на 1-3 года (LaMarche et al., 1984;
Jones et al., 1995;
Hughes et al., 1999;
Robock, 2000).
Для оценки влияния вулканизма на климат Северного Кавказа древесно кольцевые хронологии и реконструкция температуры сравнивалась с рядом лет климатически эффективных вулканических извержений, восстановленных по концентрации сульфатов в слоях ледовых колонок Гренландии (GISP2) (http://www.ncdc.noaa.gov/paleo/pubs/ipcc2007/fig614.html).
Отметим, что при визуальном обследовании кернов сосны и пихты, полученной на Кавказе, морозобойных колец обнаружено не было. Однако, было установлено, что в сводной Кавказской хронологии по ширине годичных колец 1818 год отличался самой большой депрессией прироста за весь период (индекс 0,67). Лишь к 1821 прирост стал положительным. Возможно, это понижение прироста связано с понижением летних температур в Северном полушарии после извержения вулкана Тамбора (Индонезия) в 1815 году. Такой же эффект трёхлетнего запаздывания отмечен в хронологии по ширине колец в Монголии (D’Arrigo al., 2003) и в Восточных Карпатах (Popa et al., 2008).
В хронологии по максимальной плотности колец сосны на Кавказе самая низкая температурная аномалия отмечена на год раньше, чем в хронологии по ширине годичных колец – т.е. в 1817 году. Запаздывание отклика на один год отмечено в реконструкции летних температур по максимальной плотности для Европейских Альп (Bntgen et al., 2006): самое холодное за 1250 лет лето зафиксировано в 1816 году (аномалия -4,5С). Интересно, что реконструкция температур тёплого периода для Северного Кавказа в целом согласуется с кривой вулканической активности Северного полушария (рис. 10).
Рис.10. Сравнение восстановленной температуры апреля-сентября на Северном Кавказе и вулканической активности Северного полушария.
Сравнение с региональными реконструкциями.
Грузия, Кахетия. Ближайшая к району нашего исследования дендрохронологическая площадка расположена в Грузии, Кахетия. Сравнение этой хронологии со сводной по сосне, построенной нами, показывает некоторое их сходство, как на уровне межгодовой изменчивости (r=0,22;
p0,05, N=341), так и при сравнении сглаженных средних (рис. 11). Самая высокая корреляция между рядами отмечена в 1640-1670 гг. (r=0,78;
p0,05), в один из периодов похолодания и наступания ледников на Кавказе в малом ледниковом периоде (Соломина, 1999). В то время как в 20 веке корреляция значительно снижается и достигает отрицательных значений.
Рис. 11. Сравнение хронологии «Кахетия» и сводной Кавказской хронологии, сглаженые 11-летним средним.
Статистически значимого сходства между Альпы и Татры.
реконструкциями летней температуры в ближайших районах - в Альпах (Frank et al, 2008) и в Татрах (Bntgen et al., 2007) - с кавказской реконструкцией не обнаружено. Корреляция становится значимой после 10-летнего сглаживания рядов (рис. 12). После 1976 года для рядов Альп и Татр заметно значительное повышение летних температур. Такая тенденция для Кавказской реконструкции не обнаружена. Альпийская реконструкция имеет долгопериодный тренд, который отсутствует у двух остальных кривых. Это может отражать рельные климатические тенденции, но может быть связано и с особенностями стандатризации рядов Кавказа и Татр (применение стандартизации региональной кривой для Альпийской хронологии).
Рис. 12. Сравнение реконструкций средних температур апреля-сентября для Татр (голубая линия), июня-сентября для Альп (зелёная линия), апреля-сентября для Кавказа (красная линия). Все ряды сглажены 10-летним скользящим средним. Значение коэффициентов корреляции и их статистическая значимость указаны на графике.
Реконструкции температуры Северного полушария. До сих пор в реконструкциях температуры Северного полушария не были представлены кавказские хронологии. Между тем - это важный в климатическом отношении район, расположенный на границе умеренной и субтропической зоны.
Сравнение аномалий годовых температур для Северного полушария, полученных по базе данных HadCRU3v (Brohan et al., 2006) и реконструированных для периода апрель-сентябрь по максимальной плотности колец на Кавказе за общий период 1850-2005 гг. показывет их существенное различие и, прежде всего, отсутствие долгоперидных флуктуаций у кавказской реконструкции (рис. 13).
Рис. 13. Аномалии температур апреля-сентября Северного полушария из базы данных HadCRU3v (синия линия) и реконструированные температуры апреля-сентября на Северном Кавказе (красная линия), рассчитанные относительно климатической нормы 1961-1990 гг.
Реконструкции температур по Cеверному полушарию (Wilson et al., 2007;
Moberg et al., 2005;
D’Arrigo et al., 2006) существенно различны на уровне внутривековых и межгодовых колебаний, однако их объединяет наличие положительного линейного тренда, особенно заметного на реконструкции Вильсона (Wilson et al., 2007). Кавказская реконструкция – единственная, которая не показывает этого положительного тренда. Скорее всего, это связано с ее положением в низких широтах, где и по инструментальным данным положительный тренд выражен гораздо слабее, а в некоторых районах и вовсе отсутствует (IPCC, 2007).
Заключение В климатических реконструкциях, построенных на основе денрохронологических данных, реконструкции из низких широт представлены в очень ограниченном объеме. В частности, Кавказ до сих пор не был включен в эти модели. Настоящая работа дает возможность в какой-то мере заполнить этот пробел. В ней представлена первая количественная реконструкция температуры тёплого периода, построенная по дендрохронологическим данным за период 1800-2005 гг., которая характеризует условия на севером макросклоне Центрального и Западного Кавказа. В ходе этого исследования было построено семь абсолютно сдатированных хронологий сосны и хронология пихты по ширине годичных колец, а также хронология по максимальной плотности колец сосны. Обнаружена положительная статистически значимая связь индексов прироста сосны и температуры ноября предыдущего года. Установлено, что на формирование максимальной плотности древесины влияют среднемесячные температуры тёплого периода года (апрель-сентябрь). С помощью модели линейной регрессии построена реконструкция температур апреля-сентября. Верификация модели показала, что реконструированный сигнал является устойчивым. Обнаружена также корреляционная связь между индексами хронологии по максимальной плотности колец сосны с балансом массы ледника Гарабаши. Корреляция ширины колец с речным стоком Теберды за лето позволила реконструировать и этот параметр за период с 1850-2005 гг.
В результате этой работы мы пришли к выводу, что открываются большие перспективы для дендрохронологических реконструкций на Кавказе. Мы планируем в ближайшем будущем представить наши хронологии в международную базу данных, продлить реконструкцию по плотности, по крайней мере, еще на два века, удлинить хронологии за счет исторической древесины, сохранившейся в старых домах Осетии и Кабардино-Балкарии, использовать стабильные изотопы для реконструкции других климатических параметров. Все эти данные будут использованы для лучшего понимания и более надежного моделирования изменений климата на Кавказе в прошлом и будущем.
Благодарности Автор выражает благодарность научному руководителю О.Н. Соломиной, сотрудникам отдела гляциологии ИГ РАН, и другим людям, помогавшим в сборе и обработке дендрохронологического материала и консультировавшим автора: Г.А.Носенко И.Ф.Хмелевскому, Ю.О.Юриной, В.Н.Михаленко, И.И.Лаврентьеву, И.С.Бушуевой, В.В.
Мацковскому, С.С.Кутузову, О.Е.Максимовой;
М.Хьюзу за помощь проведении лабораторных работ, за предоставление метеорологических данных А.Д. Салпагарову и А.Боку за помощь в проведении полевых работ в Тебердинском заповеднике.
Список опубликованных работ по теме диссертации:
В изданиях, рекомендованных ВАК 1. Е.Долгова, О.Соломина (2010) Первая количественная реконструкция температуры тёплого периода на Кавказе по дендрохронологическим данным. // Доклады Академии Наук, 2010, том 431, № 2.
2. Мацковский В.В., Долгова Е.А., Соломина О.Н. (2011) Применение дендрохронологических данных для реконструкции стока реки Теберды за 1850-2005 гг. // Лед и снег, 2011-1 (119) – С. 119-123.
В других изданиях 1. Долгова Е.А., Соломина О.Н, Жомелли В., Юрина Ю.О., Олейников А.Д., Володичева Н.А.(2007) Дендрохронологическое датирование морен ледников Большой Азау, Шхельда, Терскол, долина р. Баксан, Приэльбрусье. // Новые методы в дендроэкологии: Материалы Всероссийской научной конференции с международным участием (Иркутск, 10-13 сентября 2007 г.) – Иркутск:
Издательство Института географии им. В.Б. Сочавы СО РАН, 2007, 87-89с.
2. Долгова Е.А., Соломина О.Н., Салпагаров А.Д., Бок А.А.(2007) Датировка засух в период 18-20 вв. в древесно-кольцевой хронологии сосны (Pinus Hamata) в долине реки Теберды, Западный Кавказ. // VI Пущинская международная школа-семинар по экологии «Экология 2007: эстафета поколений» 13-15 декабря 2007 г.
3. Dolgova E.A., Solomina O.N., Bok A.C., Salpagarov A.D. (2007) Glacier retreat and climate change in Teberda valley, West Caucasus, Russian Federation. // Geophysical Research Abstracts, Vol. 9, 4. Solomina, O.N., Jomelli, V., Dolgova, E.A. (2008) Last millennium climate and glacier variations in the Northern Caucasus. // Geophysical Research Abstracts Vol.
10, EGU2008-A-12114.
5. Matskovskiy V.V., Dolgova E.A., Solomina O.N. (2009) Teberda valley runoff variability (AD 1797-2003) based on tree-ring reconstruction (Northern Caucasus, Russia). // Geophysical Research Abstracts, Vol. 11, EGU2009-1047-4.
6. Matskovskiy V.V., Dolgova E.A., Solomina O.N. (2009) Exploring climatic signal in pine ring width chronologies at the high elevation sites in the Northern Caucasus, Russia. // PAGES 3rd Open Science Meeting. Program and Abstracts, OC2-5.
7. Мацковский В.В., Кононова Н.К., Соломина О.Н., Долгова Е.А. (2009) Связь ширины годичных колец на Северном Кавказе с типами атмосферной циркуляции. // Современные проблемы климатологии. Материалы Всероссийской конференции посвященной 100-летию профессора О.А.
Дроздова (1909-2001) 20-22 октября 2009 г. СПб: ВВМ. 2009. С. 94-96.
8. Мацковский В.В., Долгова Е.А. (2009) Реконструкция стока реки Теберды по ширине годичных колец сосны на Северном Кавказе 1797-1999 гг. // Третья международная научная конференция молодых ученых и талантливых студентов «Водные ресурсы, экология и гидрологическая безопасность».
Сборник трудов. М.: ИВП РАН. 2009. С. 81-84.
9. Е.Долгова, О.Соломина. (2009) Реконструкция температуры тёплого периода на Кавказе по дендрохронологическим данным. Тезисы конференции «Вклад России в МПГ 2007/2008», Сочи, октябрь 2009, с.18.
10. Мацковский В.В., Е.А.Долгова. (2010) Перспективы использования годичных колец для реконструкции баланса массы ледников и речного стока на Северном Кавказе // Международный гляциологический симпозиум, 31 мая – 4 июня, 2010 г., Казань.
11. Matskovsky V.V., E.A. Dolgova and O.N. Solomina (2010) Teberda River runoff variability (AD 1850–2005) based on tree ring reconstruction (Northern Caucasus, Russia). // IOP Conference Series: Earth and Environmental Science Volume 9, 2010.
http://iopscience.iop.org/1755-1315/9/1/012017.
12. Dolgova E.A. (2010) The first quantitative warm period temperature reconstruction in the Caucasus mountains derived from tree-rings // Abstracts of the 8-th International Conference on Dendrochronology WorldDendro 2010.