Dynamics of the Regional Climatic Conditions over the Past 2000 Years on the Basis of Lithogeochemical Analysis of the Bottom Sediments of Lake Karakyol (Western Caucasus)

Full Text

Введение. В эпоху глобальных изменений климата детальные региональные палеоклиматические реконструкции являются принципиально важными для понимания антропогенного вклада и прогнозирования возможных климатических изменений.

Климатическая история Кавказа – огромной горной страны с разнообразными природными условиями – остается изученной лишь в общих чертах. Имеющиеся реконструкции в основном носят качественный характер [6, 10, 18, 23] или охватывают незначительный временной интервал [3, 16].

Важную роль в построении реконструкций занимают озерные донные отложения, играющие роль палеоархивов как биологических, так и литолого-геохимических индикаторов климата. На протяжении последних десятилетий озерные осадки по всему миру стали одним из наиболее полных и универсальных источников палеоклиматической информации. Современные высокопроизводительные аналитические методы (микро-РФА) с высоким пространственным разрешением позволили климатическим реконструкциям выйти на новый уровень. Озерные осадки способны обеспечить детальные непрерывные количественные и качественные реконструкции целого ряда параметров изменчивости среды. На Западном и Центральном Кавказе исследованием современных озерных донных осадков занимаются сотрудники отдела гляциологии Института географии РАН (Москва) в сотрудничестве с коллегами из Института геологии и минералогии СО РАН (Новосибирск). Озеро Каракель, рассматриваемое в настоящей работе, стало пионерным объектом. Другим важным объектом для палеолимнологических исследований стало озеро Донгуз-Орун в Приэльбрусье [3]. В настоящий момент ведется работа по анализу кернов, полученных в других озерах Кавказа в период 2014–2017 гг.

Постановка задачи исследований. Впервые отбор донного осадка озера Каракель был произведен в 2010 г. [11]. Впоследствии были установлены важные взаимосвязи различных параметров осадка – в частности, была выявлена устойчивая связь содержания пыльцы широколиственных видов с концентрацией брома в осадке – двух параметров, связанных с теплообеспеченностью района [12]. Дальнейшие исследования, однако, позволили выявить неточности существующей хронологии за счет получения новых кернов донного осадка в 2014 г. Настоящая работа ставит своей целью создание объединенной мастер-хронологии двух последних тысячелетий и восстановление тенденций климатической изменчивости в регионе на основе исследований состава донного осадка оз. Каракель.

Методика исследований. Отбор кернов. Озеро Каракель располагается в пределах западной высокогорной провинции Большого Кавказа, в долине р. Теберда, берущей свое начало на склонах Главного Кавказского хребта. В этом районе долина реки простирается преимущественно с юго-запада на северо-восток и имеет ширину около 1000–1200 м. Озеро расположено на абсолютной высоте 1335 м, имеет продолговатую форму и вытянуто в направлении простирания долины. Размеры озера составляют около 140 × 280 м, глубина на основной площади озера колеблется незначительно и составляет в среднем 9 м (по состоянию на 2017 г.). Озеро бессточное, однако, часть воды, вероятно, фильтруется сквозь рыхлую толщу подстилающей его морены. Со склона, особенно весной, в озеро стекают ручьи и сходят небольшие сели. На склоне над озером прослеживаются ложбины двух ручьев, отложения которых в нижней части образуют слившийся конус выноса. Основное питание озера – за счет атмосферных осадков и снеготаяния, расход – путем испарения и фильтрации. Вероятно, в питании озера существенна роль подземной составляющей. Зимой озеро покрывается льдом лишь в отдельные годы.

Впервые керновое бурение донного осадка оз. Каракель было произведено нами в августе 2010 г. Ему предшествовал пробный отбор осадка торфяным буром. Было показано, что на дне озера имеется достаточно мощный слой рыхлых отложений, причем отдельные слои различаются по своим минералогическим и физическим характеристикам. Батиметрическая съемка, предшествовавшая бурению, установила в среднем равномерную глубину на всей площади озера – 9–9.5 м. Керновое бурение производилось модифицированным поршневым буром ударного типа конструкции Атле Несье [21]. В качестве платформы использовались рекреационные металлические катамараны. Бурение производилось в центральной части озера с глубины 9.5 м. Было получено два керна с ненарушенной структурой, перекрывающих друг друга по глубине отбора (первый – с глубины от 0–150 мм до 1000–1100 мм, второй – от ~350–400 мм до 1800–1850 мм). Повторное бурение озера проводилось в августе 2014 г. также сотрудниками полевой экспедиции Отдела гляциологии Института географии. В качестве платформы был использован надувной катамаран для сплава. Получены два керна: Kar-14-1, мощность около 600 мм предположительно от границы вода–осадок (бурение без поршня); Kar-14-2, около 700 мм на глубину максимального проникновения (бурение с поршнем). На рис. 1 представлены места отбора кернов 2010 и 2014 гг.

 

Рис. 1. Русло реки Теберда и озеро Каракель с участками отбора кернов 2010 и 2014 гг.

 

Стратиграфия. Основываясь на стратиграфическом анализе двух кернов оз. Каракель 2010 г., можно проследить три характерных горизонта (рис. 2). Первый (горизонт А) до глубины 520–540 мм представлен желеобразным тонкослоистым сапропелем; второй (горизонт B) до глубины 570 мм – темно-серой, с зеленоватым оттенком, тонкослоистой и сильно обогащенной органическим веществом гиттией; третий (горизонт C) (с глубины 580 мм до конца керна 1810 мм) – голубовато-серой глиной, иногда с коричневыми прослоями, мощностью 1–2 мм. Вниз по разрезу плотность глин возрастает. Керн 2014 г. имеет длину 62 см и сложен темно-коричневым сапропелем. На глубине 20–21 см расположен характерный белесый прослой. Примечательно, что в оригинальном керне 2010 г. на глубине около 3 см присутствует очень похожий белесый прослой. Стратиграфический анализ керна 2010 г., а также имеющиеся радиоуглеродные датировки с уверенностью говорят о двух принципиально разных режимах осадконакопления на протяжении последних двух тысяч лет и в раннем и среднем голоцене в период с 11 до 4 тыс. лет назад. Настоящая работа ставит своей целью детальное рассмотрение верхней части керна 2010 и керна 2014 гг., соответствующих последним 2000 лет.

 

Рис. 2. Фотографии кернов 2010 и 2014 гг. Точками указаны участки отбора материала для датирования и соответствующий радиоуглеродный возраст. Представлены кривые концентраций элементов, обсуждаемых в настоящей работе. Полупрозрачными горизонтальными областями отмечены участки резкого роста концентраций терригенных элементов, соответствующие светлым прослоям в кернах.

 

Микро-РФА. Первые опыты применения сканирующего рентгенофлуоресцентного микроанализа с возбуждением синхротронным излучением для исследования распределения элементов в геологических образцах показали большую перспективность метода [14], в том числе для изучения морских осадочных последовательностей [15]. Впоследствии анализ широко применялся для исследования как морских, так и озерных донных осадков, и в настоящее время все методики проведения измерений аттестованы [1, 2].

Подготовка образцов для измерений проводилась по отработанной методике с изготовлением плоскопараллельных пластинок длиной 170 мм, шириной 15 мм и толщиной 2 мм из твердого препарата, пропитанного эпоксидной смолой [1].

Измерения образцов проводились на станции “Элементный анализ” ЦКП “Сибирский центр синхротронного и терагерцового излучения” [2] с помощью монохроматичного пучка синхротронного излучения с энергией 21 кэВ. Ширина пучка (поперек образца) составляла 10 мм, высота пучка (вдоль образца) – 1 мм. Образец перемещался вертикально с шагом, совпадающим с высотой пучка. Время измерения в точке составляло от 10 до 25 сек, в зависимости от загрузки спектрометрического тракта. Одновременно определялись содержания 25 породообразующих и микроэлементов: Cl, S, K, Ca, Ti, V, Cr, Mn, Fe, Ni, Cu, Zn, Ga, Ge, As, Se, Br, Rb, Sr, Y, Zr, Nb, Mo – по К сериям, U и Th – по L сериям. Кроме этого проводилось измерение отношения интенсивностей упруго- неупругорассеянного на образце излучения, как потенциально климатически обусловленного параметра [4].

Все регистрируемые величины были геометрически привязаны к положению на шкале глубины керна в миллиметрах от верха (граница вода–осадок). Верх керна датировался годом отбора и служил началом отсчета при построении временных рядов по актуальной модели: глубина–возраст осадка.

Датировки. Скорость осадконакопления. Датирование донного осадка озера происходило в несколько этапов. В ходе первого этапа были получены две датировки методом AMS (радиоуглеродная лаборатория г. Познань, Польша) из верхнего керна 2010 г. (30–31 и 52 см) и одна из базального горизонта нижнего керна 2010 г. (глубина 143 см по нижнему керну или 180.5 см от границы вода–осадок). Калиброванный возраст составил соответственно 1528–1377, 2278–2153 и 11348–11063 лет назад (табл. 1). Первые датировки подтвердили голоценовый возраст осадков горизонта С и уточнили минимальный возраст каракельской морены в долине р. Теберда [11, 12].

Впоследствии был сделан еще ряд радиоуглеродных датировок для керна 2010 и 2014 гг. Все полученные датировки представлены в табл. 1.

 

Таблица 1. Результаты датирования кернов 2010 и 2014 гг.

№	Керн	Глубина сводная, см	Возраст 14C, BP (1σ)/pMC, %	Материал
IGAN-5497	Каракель 2014	7–8	100.740±0.292	массовый образец (ТОС)
IGAN-5496	Каракель 2014	13.5–14.5	340±30	массовый образец (ТОС)
IGAN-5493	Каракель 2014	13.5–14.5	725±20	Древесный макроостаток
IGAN-5495	Каракель 2014	18–19	640±20	массовый образец (ТОС)
IGAN-5494	Каракель 2014	21–22	640±20	массовый образец (ТОС)
IGAN-5333	Каракель 2010	28.5–29	880±20	массовый образец (ТОС)
IGAN-5332	Каракель 2010	34.5–35	1030±20	массовый образец (ТОС)
IGAN-5331	Каракель 2010	42.5–43	1235±20	массовый образец (ТОС)
Poz-42587	Каракель 2010	52.5–53.5	1550±30	массовый образец (ТОС)
Poz-42588	Каракель 2010	74–75	2235±35	массовый образец (ТОС)

 

Результаты. Построение мастер-керна. Отбор образцов донных осадков с помощью поршневого озерного бура имеет свою специфику. Для того, чтобы поршень начал продвигаться в пробоотборной трубке, осадок должен быть достаточно плотным. В случае, когда верхняя часть осадка находится в полужидком состоянии, возможны потери части керна. Изначальная работа над керном 2010 г. не предполагала значительных потерь верха, а ограниченное количество радиоуглеродных датировок не позволило точно установить это на первом этапе работы. Таким образом, хронология, описанная в [11] опиралась на предположение о том, что верх керна датируется 2010 г.

Впоследствии, новые датировки, полученные для верхней части керна 2010 г., поставили под сомнение существующую возрастную модель и послужили предпосылками к предположению о том, что верхняя часть керна 2010 г. была потеряна в процессе отбора. Бурение 2014 г. производилось без поршня, что исключало вероятность потери верхней слабо консолидированной части осадка. Сравнение двух кернов показало, что на глубине 20–30 мм в керне 2010 г. и на глубине 200–220 мм в керне 2014 г. наблюдается характерный белесый линзовидный прослой, по всей видимости, аллохтонного происхождения (см. рис. 2).

Для точного совмещения кернов и построения опорного разреза были использованы результаты сканирующего микро-РФА, полученные для разрезов 2010 и 2014 гг. Совмещение профилей сканирования показало совпадение маркирующих пиков для всех аналитических данных. На рис. 3 представлен профиль Ti в совмещенном мастер-керне осадков озера.

 

Рис. 3. Совмещение профилей содержания Ti в кернах Каракель 2010 и Каракель 2014 использовано для учета недостающего верзха в керне Каракель 2010 и построения итогового мастер-керна донных осадков озера Каракель.

 

Итоговое удлинение верхнего керна с коэффициентом 1.17 и сдвиг относительно нижнего на 225 мм позволили объединить оба керна в один опорный разрез. Изменение размеров верхнего керна связано с тем, что отбор проводился в разных точках озера, где линейные скорости осадконакопления немного различаются. Но при этом текстура и состав осадков в кернах совпадают, а их изменения отражают реакцию системы седиментации озера в целом на изменения параметров внешней среды.

Датировка мастер-керна. Создание возрастной модели было произведено с помощью алгоритма Bacon программной среды R (http://chrono.qub.ac.uk/blaauw/bacon.html). Bacon использует Байесовскую вероятность для того, чтобы разбить толщу осадка на множество мелких секций, для каждой из которых дискретно определяется скорость осадконакопления. Для построения возрастной модели с помощью Bacon используется калибровочная кривая IntCal13 [22]. Одна из датировок (датировался макроостаток, тогда как во всех остальных случаях – массовые образцы осадка – total organic carbon (TOC)) была исключена из построения модели. Хотя, как правило, датировки по макроостаткам считаются более точными, в нашем случае оказалось, что именно эта датировка выбивается из общего тренда. Это может быть связано с разными причинами, в том числе и с переотложением материала, взятого для датирования. Отклонение этой датировки от результирующей взвешенной небольшое и радикально не меняет оценку скорости осадконакопления, что позволяет использовать отношение глубины и возраста осадка любого участка керна (рис. 4).

 

Рис. 4. Временные модели мастер-керна донных осадков озера Каракель: (а): модель 1 – без учета турбидитных слоев; (б): модель 2 – с учетом турбидитных слоев.

 

Данные С-14 подтверждают правильность построения мастер-керна и дают возможность построить простую линейную модель: глубина керна – возраст осадка с оценкой средней линейной скорости осадконакопления величиной равной примерно 0.33 мм/год.

Поиск климатических индикаторов. Принятая возрастная модель позволяет датировать верхние слои осадка (в предположении, что верх соответствует году отбора – 2014 г.) и провести сравнение данных аналитической микростратиграфии (профили содержаний элементов) с данными региональных метеонаблюдений – метеостанция Теберда (временной ряд 1927–2010 гг.) (www.meteo.ru).

Шаг опробования при проведении сканирующего анализа составлял величину 1.1 мм, что при средней скорости осадконакопления 0.33 мм/год соответствует 3.6 года. Это определяет минимальную величину временного разрешения полученных результатов. Однако, поскольку осадок не имеет следов ежегодной стратификации, то имеется и усреднение состава осадка в процессе осаждения. Поэтому можно предполагать, что реальное временное разрешение аналитических данных составляет величину более 5 лет. Поэтому для сравнения метеоданных с геохимическими индикаторами использовали средние 5-ти и 11-тилетние величины температуры и количества осадков.

Двигая временную шкалу по глубине керна, находим, что корреляция метеоданных с содержанием элементов достигает максимума при положении начального года сравнения (2010) на глубине 2–3 мм. Можно предположить, что линейная скорость осадконакопления для верхних слоев выше средней и оценивается величиной порядка 1 мм/год, что и обусловливает положение 2010 г. относительно верха керна. На рис. 5 представлен график изменения коэффициентов корреляции для 5-ти и 11-тилетних температур с рядом геохимических индикаторов.

 

Рис. 5. Изменение коэффициентов корреляции для 5-ти и 11-тилетних температур с геохимическими индикаторами в зависимости от положения начала года сравнения (2010 г.) на шкале глубины керна.

 

Важным фактором является наличие высоких (до +0.95, значимых при p < 0.01) коэффициентов корреляции для сопоставляемых временных рядов длиной более 80 точек (1927–2010 гг.), что свидетельствует о наличии устойчивых и значимых связей между наблюдаемыми величинами (метеопараметрами и составом донного осадка на этом же временном интервале). В табл. 2 представлены данные о наборе породообразующих и микроэлементов, имеющих высокую корреляцию с температурой и атмосферными осадками. Установленная связь дает основания для расчета трансферных функций и построения палеореконструкций по данным аналитической микростратиграфии.

Определение трансферных функций. Набор элементов, имеющих максимальную корреляцию с 5-ти и 11-тилетними температурами, был использован для расчета уравнений регрессии количественной связи состава донного осадка и реконструируемого метеопараметра. В данном случае использован набор простых трансферных функций для каждого элемента в отдельности. Пример расчета для Rb и Ti представлен на рис. 6.

 

Рис. 6. Примеры получения трансферных функций для Rb и Ti, устанавливающих количественную связь между содержаниями элементов в датированных слоях донного осадка оз. Каракель и 11-летними температурами региона.

 

Подготовка литолого-геохимических данных для построения реконструкций. Перед аппроксимацией полученные трансферные функции на всю глубину опробования керна необходимо провести дополнительный анализ временных рядов литолого-геохимических данных, используемых для реконструкций. Это связано с тем, что в случае значительных отличий в содержаниях элементов на интервале обучения от остальных участков керна, построенные реконструкции метеопараметров не будут отвечать не только количественным, но и качественным оценкам. Нами не выработаны строгие оценки допустимых отклонений, но предыдущий опыт показывает, что при вариации содержаний элементов по глубине керна более чем на 50% от верхнего интервала обучения, необходимо провести дополнительное исследование для оценки изменения модели осадконакопления в прошлом.

 

Таблица 2. Набор породообразующих и микроэлементов, имеющих максимальную корреляцию с метеоданными. Т – среднегодовая температура, Pr – годовое количество осадков

Параметр	K	Ti	Zn	Rb	Y	Nb
5-ти лет Т	0.69	0.78	0.74	0.82	0.80	0.85
11-ти лет Т	0.72	0.84	0.76	0.93	0.89	0.95
5-ти лет Pr	0.51	0.63	0.72	0.75	0.77	0.72
11-ти лет Pr	0.48	0.73	0.74	0.83	0.83	0.74

 

В нашем случае наблюдаются явные аномальные повышения содержаний большинства элементов в двух интервалах мастер-керна: 164–230 мм и 544–566 мм (см. рис. 2). По-видимому, эти интервалы связаны с разовыми вбросами большого количества терригенного материала – аналогами турбидитных отложений в морях и океанах. Также фактором в пользу событийного происхождения рассматриваемых интервалов служит совпадение оценок возраста для образцов, взятых с глубин 180 и 210 мм. Радиоизотопный возраст для обеих глубин дает одинаковую оценку 640 лет назад, что логично, если это одномоментный вброс вещества одного состава.

Таким образом, есть основания удалить два интервала, нарушающие непрерывную хронологию, из мастер-керна донных осадков оз. Каракель, после чего построить новую временную модель – модель 2 на рис. 4. Все дальнейшие реконструкции сделаны на временных рядах, построенных на этой модели.

Построение реконструкций изменения 11-тилетней температуры региона на интервале опробования. С использованием рассчитанных трансферных функций был построен набор реконструкций изменения 11-тилетней температуры региона на временном интервале последних тысячелетий. На рис. 7а, в качестве примера, представлена часть реконструкций, рассчитанных по временным рядам отдельных элементов: двух породообразующих – K и Ti и одному микроэлементу – Rb. Также представлена суммарная реконструкция, объединяющая все использованные элементы.

 

Рис. 7. (а) – реконструкции средней 11-тилетней температуры по литолого-геохимическим данным донных осадков оз. Каракель; (б) – сравнение одной из реконструкций, полученной по содержанию RB, с метеоданными на интервале обучения (1927–2010 гг.); (в) – температурные реконструкции для Северного полушария по литературным данным.

 

Обсуждение. Вариации терригенной, биогенной и аутигенной составляющих в озерных осадках характеризуют изменения окружающей среды и климата на протяжении времени их формирования, однако интерпретация этой информации очень сложна и, зачастую, неоднозначна. Статистический анализ показывает чрезвычайно высокую корреляцию между содержанием химических элементов в керне оз. Каракель – K, Ti, Zn, Rb, Y, Nb, U и метеорологическими параметрами – среднегодовым количеством осадков и среднегодовой температурой. Эта связь положительна для всех представленных терригенных элементов. Наличие подобных корреляционных связей объяснимо, так как химический состав компонентов глинистой фракции осадка, например коэффициент химического выветривания, позволяет судить о соотношении в осадке первичных и вторичных минералов, и косвенно связан с изменениями климата [7]. В работе [5] показано, что среди геохимических индикаторов особенно информативными являются такие породообразующие элементы, как K, Ca, Ti, Fe, Mn и микроэлементы, характеризующие терригенную (Cu, Ga, Rb, Sr, Y, Zr, Nb, Sn, Ba) и органогенную (Br, I, Mo, U) компоненты осадка, а также соотношение Sr/K, Ca/Sr, Br/Rb и др. На практике эти зависимости удалось использовать для реконструкции температуры на Алтае – по соотношению стронция к рубидию и содержанию брома и титана в осадках Телецкого озера. Было показано, что соотношение Sr/Rb является индикатором доли невыветрелой обломочной фракции и отражает температуру сезона снеготаяния, содержание брома связано с биопродуктивностью на водосборе [5]. Ранее [12] было проведено сопоставление отношения Sr/Rb в приледниковых водах (оз. Турье, р. Алибек) и в оз. Каракель и его притоке. Оно показало, что разница составляет несколько порядков (0.4 и 4.1 в высокогорье и 132 и 215 в Каракеле и впадающем в него ручье). Очевидно, что интенсивность выветривания, о которой косвенно свидетельствует соотношение Sr/Rb, зависит от температуры, а интенсивность поступления продуктов выветривания в озеро – от количества осадков [17]. Сопоставление концентрации брома в сапропеле оз. Каракель с данными пыльцевого анализа показывает явное сходство хода этого показателя с колебаниями общего содержания пыльцы как древесных, так и травянистых растений. Лучше всего общее содержание брома согласуется с ходом кривой концентрации пыльцы широколиственных пород. Обе характеристики могут быть связаны с теплообеспеченностью, которая положительно влияет и на пыльцевую продуктивность, и на содержание брома в осадке [12].

Для каждого водоема набор геохимических индикаторов, которые наиболее тесно связаны с изменениями окружающей среды, индивидуален. В этой работе обнаружены наиболее тесные корреляционные связи температуры и осадков с содержанием в керне терригенных элементов K и Ti и микроэлемента Rb. Показательно, что кривые содержания этих элементов в осадке имеют очень высокое сходство между собой, что дает основание использовать их осреднение для большей надежности результатов.

В нашем случае отмечаются положительные корреляции геохимических параметров и с температурой, и с количеством атмосферных осадков. Эта связь не случайна: высокие температуры и повышенное увлажнение способствуют более высокой интенсивности выветривания пород и отложения продуктов этого выветривания на дне озер. Сравнительный анализ метеопараметров показывает, что среднегодовые температуры и осадки положительно коррелируют между собой за период 1927–2009 гг. с коэффициентом корреляции +0.42. Таким образом, отмечается тенденция к сочетанию: тепло–влажно, сухо–холодно. На Кавказе это характерно для зимних условий: холодная сухая антициклональная погода связана в основном со вторжениями воздуха с севера, а циклоны приносят с собой тепло и осадки.

Одна из существенных проблем при создании высокоразрешающих реконструкций на основе озерных не ламинированных осадков – надежный хронологический контроль. В этой работе мы использовали данные сканирующего РФА СИ для объединения двух кернов в одну мастер-хронологию, которая, в свою очередь, основана на 11 радиоуглеродных датировках. В нашей предыдущей работе мы предполагали, что осадконакопление в оз. Каракель происходит с постоянной скоростью, примерно равной 0.22 мм в год, которая увеличивается только в самой верхней части отложений. Более детальный анализ результатов РФА СИ показал, что этот монотонный процесс периодически нарушался, и в озеро единовременно поступало большое количество обломочного материала. Для того, чтобы хронологическая шкала не нарушалась, эти катастрофические события необходимо идентифицировать и исключить из мастер-хронологии. На основе геохимического анализа удалось выявить и удалить из разреза два таких события – 640 и более 1600 лет назад.

Следует отметить, что ранее в отложениях оз. Каракель [8] выделено несколько переотложенных диатомовых комплексов, которые определялись по характерному изменению таксономических пропорций. Таким образом, мы не можем исключить вероятности и других нарушений в процессе отложения осадков в озере, но, вероятно, если они и случались, то были меньшего масштаба, чем событие, случившееся около 640 л.н. Об этом свидетельствуют и геохимический анализ, и результаты AMS датирования, поскольку вся совокупность датировок в целом хорошо описывается линейной функцией.

Катастрофические события, связанные с дополнительным поступлением вещества в озеро, могут происходить по нескольким причинам. Стабильное формирование современных отложений могло быть нарушено антропогенной деятельностью, но в масштабах, несоизмеримых с природными событиями. Селевые потоки могут быть основным источником нарушения непрерывных осадочных последовательностей. В настоящее время в озеро поступает небольшое количество воды и наносов из р. Теберда. Но периодически, при высоком уровне воды в реке, объем воды и взвеси, в которых преобладает ледниковая составляющая, также мог поступать в озеро и откладываться в озерных осадках, резко увеличивая скорость осадконакопления. Учитывая тектоническую активность района, единовременное поступление терригенного материала в озеро может быть связано с сейсмообвалами, аналогично событиям в долине Баксана, датированным величиной 1440 ± 40 и 1400 ± 60 лет [9].

Сравнение полученной реконструкции температуры воздуха за последние 2 тыс. лет показывает, что в целом она хорошо согласуется с литературными данными. Амплитуда изменчивости температур не превосходит 2 °С, что характерно и для других районов умеренной зоны Евразии [24]. Выделяется два продолжительных периода потепления – современный и на рубеже первого и второго тысячелетия новой эры. Два периода похолодания в V–VII и XIII–XV вв. примерно совпадают с двумя холодными интервалами около 1300–1200 и 400–350 л.н., выделенными Е.В. Квавадзе и Ю.В. Ефремовым [18], которые исследовали озерные и озерно-болотные отложения в долине Архыза. Относительное похолодание отмечается на нашей реконструкции и в XIX в., когда по многочисленным данным было зафиксировано несколько наступлений ледников. XVIII в. по многим косвенным данным был теплым [23]. Несколько неожиданно выглядит на нашей кривой XVII в., который во многих районах северного полушария был холодным. По современным представлениям, в XVII в. на Кавказе было мощное наступание ледников, и оледенение было существенно больше по размерам, чем в период наступания XIX в. Возможно, эти противоречия можно объяснить некоторым сбоем в хронологии осадконакопления в оз. Каракель. Период с VIII по XII вв. соответствует “архызскому перерыву”, т.е. периоду средневекового оптимума, который выделяется на Кавказе по спорово-пыльцевым и археологическим данным [13, 23]. О климате первого тысячелетия в этом регионе известно очень мало и наша реконструкция, видимо, является первой, которая дает некоторое представление о количественных изменениях температуры на северном Кавказе в это время.

По сравнению с глобальными реконструкциями температуры Северного полушария [19, 20] полученная реконструкция температуры на Кавказе выглядит более контрастной, с ярко выраженной низкочастотной составляющей и периодами потеплений и похолоданий (рис. 7в). На ней также гораздо более явственно выражен средневековый оптимум. Также как и на полушарной кривой, на Кавказе примерно с XVII в. начинается тренд к потеплению. Современное потепление по своим масштабам выглядит примерно равным потеплению в период пика средневекового оптимума и в теплое время начала 1 тыс. н.э.

Заключение. Использование высокоразрешающих современных методик исследования (микро-РФА СИ) донных осадков оз. Каракель позволило успешно решить следующие задачи:

• построить мастер-керн путем объединения кернов 2010 и 2014 гг.;
• выявить период массового поступления аллохтонного материала в озеро и внести корректировку в кривую скорости осадконакопления;
• установить наличие корреляционных связей между среднегодовыми температурами воздуха и осадками (гмс Теберда) с одной стороны, и содержанием в озерных отложениях геохимических индикаторов с другой;
• построить реконструкцию среднегодовых температур за последние 2 тыс. лет с 11-летним осреднением и шагом опробования, соответствующим 5–10 годам.

Реконструкция температуры воздуха за 2 тыс. лет с таким разрешением на Кавказе выполнена впервые. В будущем планируется провести высокоразрешающие реконструкции изменений природной среды в этом районе на основе спорово-пыльцевого и диатомового анализов и сопоставить полученные результаты. Это даст возможность уменьшить неопределенности, связанные с новым методом палеореконструкций, использованным в этой работе.

Финансирование. Работа выполнена при поддержке грантов РФФИ № 17-35-50134 и частично № 16-05-00647 (построение временной модели), ГЗ № 0148-2018-0005 (построение реконструкции) текущего Государственного задания (ИГ РАН). Использовалось оборудование ЦКП “СЦСТИ” на базе ВЭПП-3/ВЭПП-4М/НЛСЭ ИЯФ СО РАН, поддержанного Минобрнауки России (уникальный идентификатор проекта RFMEFI62117X0012)”.

Funding. The research was supported by RFBR grant no. 17-35-50134, partially no. 16-05-00647 (creating age-depth model), the Government Contract no. 0148-2018-0005 (creating reconstruction), the current Government Contract (Institute of Geography RAS). The infrastructure of the Shared-Use Center “Siberian Synchrotron and Terahertz Radiation Center (SSTRC)” based on VEPP-3/VEPP-4M/NovoFEL of BINP SB RAS was employed funded by Ministry of Education and Science of the Russian Federation (project RFMEFI62117X0012).

About the authors

M. Yu. Alexandrin

Institute of Geography, Russian Academy of Sciences

Author for correspondence.
Email: alexandrin@igras.ru
Russian Federation, 29, Staromonetny, Moscow, 119017

A. V. Darin

Sobolev Institute of Geology and Mineralogy, Siberian Branch, Russian Academy of Sciences

Email: alexandrin@igras.ru
Russian Federation, 3, Koptyug prospect, Novosibirsk, 630090

A. M. Grachev

Institute of Geography, Russian Academy of Sciences

Email: alexandrin@igras.ru
Russian Federation, 29, Staromonetny, Moscow, 119017

O. N. Solomina

Institute of Geography, Russian Academy of Sciences

Email: alexandrin@igras.ru
Russian Federation, 29, Staromonetny, Moscow, 119017

References

Supplementary files