Известия РАН. Серия географическая, 2021, T. 85, № 1, стр. 97-108

Количественная реконструкция годовых температур воздуха Алтайского региона за последние 1400 лет по данным аналитической микростратиграфии ленточных глин оз. Кучерлинское

А. В. Дарьин a*, Г. Чу b, Ц. Сан c, В. В. Бабич a, И. А. Калугин a, Т. И. Маркович a, В. С. Новиков a, М. А. Максимов d, Ф. А. Дарьин e, Д. С. Сороколетов e, Я. В. Ракшун e, А. А. Гогин f, Р. А. Сенин f

a Институт геологии и минералогии им. В.С. Соболева СО РАН
Новосибирск, Россия

b Институт геологии и геофизики АН Китая
Пекин, Китай

c Национальный исследовательский центр геоанализа
Пекин, Китай

d Институт нефтегазовой геологии и геофизики им. А.А. Трофимука СО РАН
Новосибирск, Россия

e Институт ядерной физики имени Г.И. Будкера СО РАН
Новосибирск, Россия

f Национальный исследовательский центр “Курчатовский институт”
Москва, Россия

* E-mail: avd@igm.nsc.ru

Поступила в редакцию 24.04.2020
После доработки 24.08.2020
Принята к публикации 30.10.2020

Полный текст (PDF)

Аннотация

Приледниковое озеро Кучерлинское (Алтай) содержит ленточные глины, что позволяет строить точные возрастные модели подсчетом слоев и верифицировать их изотопными данными Cs-137, Pb-210 и C-14. Результаты сканирующего микро-РФА на пучках синхротронного излучения были использованы для построения временных рядов более 20 породообразующих и следовых элементов в разрезе донных отложений. Сравнение с метеоданными показало, что коэффициент корреляции годовой температуры с реконструкцией на интервале обучения (1940–2016 гг.) равен +0.59, что является значимой величиной (n = 76, p = 0.99). При использовании 10-летних сглаженных данных коэффициент корреляции достигает величины +0.84, при этом погрешность температурной реконструкции снижается до величины ±0.52°C. Сравнение полученной реконструкции с глобальными для Северного полушария для последних 1400 лет показывает наличие общих трендов и экстремумов.

Ключевые слова: Алтай, озеро Кучерлинское, донные осадки, геохимия, микро-РФА, синхротронное излучение, реконструкция температуры

ПОСТАНОВКА ПРОБЛЕМЫ

Реконструкции температурных аномалий высокого временного разрешения для Северного полушария за последние тысячелетия, опубликованные в XXI в. [17, 21, 23, 24, 27, 28, 30], показывают принципиальное изменение представлений о климате прошлого, ранее обобщенных М. Манном в 1999 г. [26].

Переход от “хоккейной клюшки” Манна с так называемым беспрецедентным глобальным потеплением XX в. до почти синусоидальных (квазициклических) вариаций температуры последних тысячелетий обусловлен появлением количественных региональных температурных реконструкций высокого разрешения длиной более 500 лет, которые были использованы при построении глобальных климатических реконструкций. При этом, по нашему мнению, используемые реконструкции должны отвечать следующим требованиям: воспроизводимость исходных данных; высокое временнóе разрешение (год на тысячелетней шкале); количественная калибровка по инструментальным метеоизмерениям; реалистичная оценка погрешностей как временной шкалы, так и реконструируемых климатических параметров.

Озерные донные осадки во многих случаях являются палеоархивами, позволяющими получить заданные параметры реконструкций. Для исследований с высоким временн${\text{ы}}'$м разрешением необходимо использовать тонкослоистые донные осадки, поскольку временнóе разрешение для неслоистых отложений обычно превышает десятки лет. Особенный интерес исследователей во всем мире вызывают донные отложения с годовыми слоями (варвами), обеспечивающими при должном подходе и правильной интерпретации требуемый уровень детализации и точности реконструкций [33]. Ленточные глины (varved clays) встречаются во многих кайнозойских отложениях Горного Алтая и представляют большой интерес для гляциологических исследований [10].

В приледниковых озерах ленточные глины образуются сезонными поставками ледниковой мути и состоят из крупных алевритовых частиц (весенне-летний слой) и мелкозернистой глины, осаждающейся в зимний период. Использование современных аналитических методов позволяет исследовать такие донные отложения с годовым, а иногда и сезонным разрешением, что открывает новые подходы в палеоклиматических исследованиях [13, 15]. Необходимым требованием к таким работам является использование инструментальных аналитических методик с субмиллиметровым пространственным разрешением, таких как гиперспектральная спектроскопия [16], микро-РФА [18, 29] и микро-РФА с синхротронным возбуждением [5].

Высокое временнóе разрешение позволяет проводить прямое сравнение набора полученных прокси с данными региональных метеонаблюдений на интервале последнего столетия с нахождением климатических индикаторов и построением трансферных функций для реконструируемых параметров [2]. Использование фокусирующей рентгеновской оптики в РФА позволяет перейти на уровень субмикронного пространственного разрешения и исследовать внутреннее строение отдельных годовых слоев с получением информации о сезонных вариациях климатических прокси и построением литолого-геохимических моделей годового цикла осадкообразования [6].

Исследование современных озер с ледниковыми глинами позволяет изучать изменения окружающей природной среды на временны́х интервалах последних тысячелетий с годовым временны́м разрешением. Одним из ключевых преимуществ таких объектов является возможность построения временны́х шкал, по точности и правильности сопоставимых с дендрохронологическими. Находки новых озер и быстрое развитие современных методов анализа открывают высокий потенциал ленточных глин для палеоклиматических реконструкций [31].

В данной работе представлена 1400-летняя количественная реконструкция среднегодовых температур воздуха Алтайского региона, полученная при исследовании ежегодно ламинированных донных осадков приледникового оз. Кучерлинское с использованием современных методик аналитической микростратиграфии – сканирующего рентгенофлуоресцентного анализа с возбуждением синхротронным излучением (РФА СИ) с субмиллиметровым пространственным разрешением.

Район исследований, геология, геоморфология

Озepo Кучерлинское располагается в южной части Горного Алтая на северном макросклоне Катунского хребта. Катунский хребет (2800–4500 м) – один из крупнейших центров оледенения Алтая. Высота снеговой линии на склонах Катунского хребта составляет порядка 2600 м [8]. С вершины и склонов высочайшей точки Алтая г. Белухи (4506 м) спускаются несколько крупных долинных ледников, один из которых – ледник Кучерлинский.

Кучерлинское озеро расположено в верховьях долины р. Кучерлы в гребневой части Катунского хребта на территории природного парка “Белуха” и имеет ледниковое питание. Основной объем воды поступает с Кучерлинского ледника. Озеро расположено на высоте 1790 м над уровнем моря. Длина озера составляет 5 км, ширина – 1 км, максимальная глубина доходит до 45 м. Общая площадь 321 га. Озеро является проточным; в него впадает р. Кучерла, а также три небольших ручья, которые стекают с западного и восточного склонов котловины. Вытекает из озера одна река – Кучерла (рис. 1). На расстоянии 20 км от озера находится метеостанция Кара-Тюрек с непрерывными рядами метеонаблюдений с 1940 г.

Рис. 1.

Космический снимок района оз. Кучерлинское. Звездочкой отмечено место отбора керна в 2018 г.

Породы территории водосборного бассейна Кучерлинского озера в основном представлены аккемской и кучерлинской свитами, сложенными слоистыми зеленовато-серыми полевошпатово-кварцевыми известковистыми песчаниками, алевропесчаниками и алевролитами [12].

Климат региона

Горный Алтай имеет умеренно-континентальный климат с холодной зимой и теплым летом. Существенное влияние на климат региона оказывает рельеф, который образует вертикальную климатическую зональность – зону низкогорного климата (до 500–600 м), зону среднегорного климата (от 500 до 1500 м и более) и зону высокогорного климата (свыше 2000–2500 м). Необходимо отметить, что в пределах Горного Алтая с северо-запада на юго-восток происходит аридизация климата, наблюдается уменьшение количества выпадающих осадков [9].

В исследуемом регионе, по данным метеостанции Кара-Тюрек, среднегодовые температуры за период 1940–2017 гг. колебались от –7.6 до –3.4°С. Максимальные среднемесячные летние температуры достигали отметки 9.9°С, а зимние доходили до –25.2°С . При этом количество атмосферных выпадений колебалось от 835 до 368 мм (http://meteo.ru/data). Метеостанция Кара-Тюрек находится на 800 м выше оз. Кучерлинского (2601 м над ур. м.), поэтому абсолютные величины инструментальных метеоданных могут отличаться от территории озера, но относительные изменения и общий тренд погодно-климатических параметров совпадает.

МЕТОДИКИ ИССЛЕДОВАНИЙ

Отбор керна и пробоподготовка

Керн донных осадков оз. Кучерлинское был получен со льда с использованием ударной трубки в марте 2018 г. Длина керна, отобранного в самой глубокой части озера с глубины 45 м, составила 1120 мм. В процессе получения керна особое внимание обращалось на сохранность границы вода–осадок, что позволяло в дальнейшем вести отсчет годовых слоев от года отбора.

После отбора трубка с керном была плотно закрыта снизу и оставлена в вертикальном положении на несколько часов в теплом помещении. Затем вода из верхней части трубки была удалена и верх плотно закрыт без деформации керна. Керн был доставлен в лабораторию Института геологии и минералогии СО РАН (Новосибирск) в вертикальном положении. Проведенные манипуляции позволили максимально сохранить верхнюю часть керна от повреждений.

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

Из твердых препаратов для всего интервала керна были изготовлены плоскопараллельные образцы для микро-РФА толщиной 2 мм и оптические шлифы, перекрывающие три интервала керна с визуально наблюдаемыми слоями: верхний (0–170 мм), средний (480–650 мм) и нижний (940–1110 мм).

Варвохронология

Подсчет слоев проводился по фотографиям исходного (влажного) керна на интервалах с визуально выделяемыми слоями и по фотографиям оптических шлифов на трех выше указанных интервалах. Фотографии оптических шлифов были сделаны в двух вариантах – в отраженном свете без увеличения (оптический сканер) и в проходящем свете в оптическом микроскопе с 2.5-кратным увеличением. Подсчет проводился по визуально выделяемым парам – светлый–темный слой – двумя исследователями, независимо друг от друга. Полученные данные усреднялись и использовались для построения возрастной модели. Детально процедуры описаны в [19].

В ряде спорных случаев в дополнение к визуальным данным при подсчете слоев были использованы данные сканирующего микро-РФА на пучках СИ с субмикронным пространственным разрешением [3]. Исследования отдельных интервалов керна проводились с использованием рентгеновской концентрирующей оптики, позволяющей фокусировать пучок возбуждающего излучения до размеров 15–20 микрон. Используемый шаг сканирования в 20–100 мкм позволял получить профили изменения содержания элементов внутри отдельных годовых слоев и находить границу их раздела.

Наиболее эффективным геохимическим индикатором границ годовых слоев в приледниковых пресноводных озерах является отношение Rb/Sr, отражающее размер частиц в осадке. Обоснование выбора Rb/Sr отношения как индикатора годовых слоев в осадках приледниковых озер рассмотрено ранее в [6]. Подсчет слоев в осадках оз. Кучерлинское с использованием отношения Rb/Sr показан в работе [4].

Изотопные исследования

Распределения активности изотопов Cs-137 и Pb-210, используемые для верификации датировки верхнего интервала керна (0–200 мм), сделаны в Аналитическом центре Института геологии и минералогии СО РАН (Новосибирск) с использованием полупроводниковой низкофоновой гамма-спектрометрии на коаксиальном Ge-детекторе с низкофоновым криостатом по стандартной методике (аналитик М.С. Мельгунов). Результаты представлены на рис. 2.

Рис. 2.

Распределение активности изотопов Cs-137 и Pb-210 в верхней части (интервал 0–200 мм) керна донных отложений оз. Кучерлинское.

Для верификации возрастной модели на всей глубине керна были сделаны анализы двух образцов керна на радиоуглеродный возраст (С-14). Анализ образца, взятого в интервале глубин 1050–1060 мм от верха керна, выполнен в ЦКП “Геохронология кайнозоя”, Новосибирск (https://archaeology.nsc.ru/proekty/ckp_gk/). Образец 1000–1010 мм исследован в лаборатории Beta analytics, Майами, США (https://www.radiocarbon.com/).

Аналитическая микростратиграфия

Методики аналитической микростратиграфии донных осадков с использованием сканирующего РФА на пучках синхротронного излучения (РФА СИ) были развиты авторским коллективом в последние десятилетия и подробно представлены в ряде публикаций [5, 7, 20].

Для анализа были использованы 2-мм плоскопараллельные образцы (пластинки) с пришлифованной поверхностью, приготовленные из твердых препаратов донных осадков. При этом пластинки для микро-РФА и для оптических шлифов были вырезаны из одного блока и имели одинаковые последовательности годовых слоев. Это давало возможность дополнительного контроля за точностью привязки данных микро-РФА к датированным слоям донных осадков.

Всего были получены данные о распределении 24 породообразующих и следовых элементов: K, Ca, Ti, V, Cr, Mn, Fe, Ni, Cu, Zn, Ga, Ge, As, Se, Br, Rb, Sr, Y, Zr, Nb, Mo, Pb, Th, U. Также регистрировалось отношение упруго-/неупруго-рассеянного на образце излучения (Coh/Incoh) как показателя рентгеновской плотности исходного керна. Таким образом, в каждой точке сканирования регистрировался набор из 25 литолого-геохимических индикаторов (мульти-прокси), отражающих изменение состава донных осадков озера во времени.

В данной работе представлены результаты аналитических измерений, полученные в Сибирском центре синхротронных и терагерцовых исследований (ИЯФ СО РАН, Новосибирск, https://ssrc.biouml.org/#!), Курчатовском комплексе синхротронно-нейтронных исследований (КИСИ-Курчатов, Москва, http://kcsni.nrcki.ru/) и Шанхайском синхротронном центре (SSRF, http://www.ssrf.ac.cn/) в 2018–2020 гг.

РЕЗУЛЬТАТЫ ИССЛЕДОВАНИЙ

Построение возрастной модели и временн$ы'$х рядов геохимических данных

Полученные изотопные данные позволяют оценить положение максимальной активности Cs-137 в 160–170 мм от границы вода–осадок (см. рис. 2). Это соответствует времени глобальных атмосферных выпадений после ядерных испытаний на Новоземельском полигоне 1961–1962 гг. [14]. Данные варвохронологии (подсчет слоев) в нескольких вариантах дают для глубины 165 мм разброс оценок от 1957 до 1965 г. При этом среднее значение 1961 г. ±4 полностью совпадает с изотопными данными, что подтверждает предположение о годовой природе наблюдаемых слоев.

Возрастная модель для интервала керна 0–200 мм (рис. 3б) была построена по данным 3-кратного независимого подсчета визуально выделяемых годовых слоев на этом интервале. Правильность полученной модели подтверждается хорошим совпадением с результатами изотопных анализов – распределением активности Cs-137 и Pb-210. Подробно построение возрастной модели на этом интервале рассмотрено в [19].

Рис. 3.

(а) – возрастная модель (глубина керна – возраст слоя осадка) для керна Куч_2018 на полную глубину опробования. Модель построена по данным варвохронологии и верифицирована датировками по С-14. “ЦКП–1055” – анализ образца с глубины 1050–1060 мм выполнен в ЦКП “Геохронология кайнозоя”. “Beta_An-1005” – анализ образца с глубины 1000–1010 мм выполнен в Beta Analytic.inc. (б) – возрастная модель для верхнего интервала керна 0–200 мм по данным варвохронологии и изотопных анализов Cs-137.

Возрастная модель на интервале керна 200–1120 мм (рис. 3а) построена по данным подсчета визуально выделяемых слоев на фото исходного (влажного) керна и исследования шлифов на интервалах 480–650 и 940–1110 мм. Полученные данные были усреднены и аппроксимированы экспоненциальной функцией. Правильность модели подтверждается совпадением с радиоуглеродными датировками.

С использованием полученной возрастной модели был проведен пересчет линейной глубины керна в возрастную (лет назад). После этого дробные значения возраста были переведены в целочисленные значения. Все данные об элементном составе внутри одного года усреднялись и приписывались этому году в виде одного значения. Поскольку отбор керна осуществлялся в марте 2018 г. и не было отмечено потерь верхней части, то начало временной шкалы (верх керна) датируется 2017 г. Таким образом были построены временн${\text{ы}}'$е ряды для всех анализированных элементов.

Создание трансферной функции, реконструкция температуры

Величина шага сканирования определяется средней линейной скоростью (мм/год) осадконакопления на данном интервале керна и подбирается такой, чтобы на интервал одного года приходилось 3–5 точек. В каждой точке получается набор концентраций более 25 элементов. Из них сначала отбраковываются данные по критерию – статистическая значимость аналитического сигнала (отношение пик/фон в исходной спектрограмме). Оставшийся массив аналитических данных, привязанных к линейной шкале (мм глубины керна от границы вода–осадок), пересчитывается во временну́ю шкалу по возрастной модели с годовым временн${\text{ы}}'$м разрешением. Для каждого линейного интервала керна, соответствующего одному году, проводится усреднение аналитических данных, попавших в этот интервал. Таким образом, исходные данные сканирования керна превращаются во временные ряды, где каждому году приписан ряд аналитических данных.

Для поиска геохимических индикаторов, коррелированных со среднегодовыми региональными температурами, были использованы инструментальные данные ближайшей метеостанции Кара-Тюрек. С использованием методических приемов, подробно описанных в [1], методом множественной регрессии построена функция, связывающая среднегодовую температуру с литолого-геохимическими данными. Сравнение проводилось на годовой временнóй шкале. Полученная функция имеет следующий вид: T = 0.602*Ti + + 0.034*Ni – 0.742*Sr + 0.754*Y + 1.46*Nb – ‒ 0.795*Mo – 6.762. Символы элементов обозначают концентрацию этого элемента в г/т (ppm), приписанную интервалу керна, соответствующему годовому слою при пересчете миллиметровой шкалы в годовую по используемой возрастной модели.

В рамках данной статьи авторы не пытаются рассматривать полученную функцию как цифровую литолого-геохимическую модель годового цикла осадкообразования в озере Кучерлинское. Общая модель связи элементного состава донных осадков с внешними (климатическими) воздействиями понятна и сводится к простым объяснениям. Изменение температуры и количества атмосферных осадков приводит к вариациям соотношения терригенной, органогенной, хемогенной (для соленых озер) и аэрозольной компонент донных отложений. Каждая компонента осадка имеет специфический элементный состав, и вариации соотношений компонент приводят к вариациям состава донных осадков. Для данного объекта ничтожна доля органики, озеро пресное – поэтому нет хемогенной компоненты, вклад аэрозоля невелик. Основные вариации состава донных отложений связаны с интенсивностью поступления терригенного материала – ледникового происхождения и сносимого с бортов озера. Основными внешними факторами являются температура и атмосферные осадки. Какой из факторов будет иметь наибольший отклик определяется статистическими методами. В данном случае уравнение с более высокой корреляцией построено для температуры.

Использование полученной зависимости (трансферной функции) дает значимый коэффициент корреляции (+0.54 для n = 77) между инструментальной среднегодовой температурой воздуха на интервале обучения (1940–2016 гг.) и реконструкцией по составу донных осадков для этого же временно́го интервала с ошибкой в 1.3°C для 95% вероятности. Формально можно рассматривать найденную закономерность как свидетельство наличия устойчивой линейной зависимости элементного состава донных осадков оз. Кучерлинское от изменяющихся внешних условий. При этом цифровым параметром внешних условий служит среднегодовая температура воздуха.

Использование сглаженных 10-летних величин позволяет увеличить коэффициент корреляции до +0.84, при этом уменьшив величину погрешности до 0.52°C. На рис. 4 представлен сглаженный метеоряд на интервале обучения, реконструкция по трансферной функции и интервал погрешности.

Рис. 4.

Сглаженные 10-летние температуры на интервале 1940–2016 гг. (метеостанция Кара-Тюрек), реконструкция по геохимическим данным, 95% интервал погрешности реконструкции.

Применение трансферной функции на всю глубину опробования керна позволило получить температурную реконструкцию для временнóго интервала 1940–576 г. (рис. 5а).

Рис. 5.

Сравнение температурной реконструкции с литературными данными. (а) – реконструкция по литолого-геохимическим данным осадков оз. Кучерлинское (данная работа); (б) – Северное полушарие, A. Moberg, 2005 [28]; (в) – Китай, B. Yang, 2002 [32]; (г) – С.‑В. Европа, V. Klimenko, 2014 [25]; (д) – Северное полушарие, B. Christiansen and F.C. Ljungqvist, 2012 [17]; (е) – Северное полушарие, IPCC, 2013 [24].

ОБСУЖДЕНИЕ РЕЗУЛЬТАТОВ

На рис. 5 представлено сравнение полученной нами региональной температурной реконструкции с рядом обобщенных глобальных реконструкций [17, 24, 25, 28, 32]. Для удобства сравнения все литературные реконструкции были нормированы методом минимакс от 0 до 1. Несмотря на значительную разницу в исходных данных, усреднениях и сглаживаниях, все реконструкции имеют ярко выраженные общие тренды и экстремумы.

Можно выделить средневековое потепление, по амплитуде сопоставимое с современным глобальным. Также отчетливо во всех реконструкциях проявлен малый ледниковый период, хотя и имеющий немного отличий на временнóй шкале. Выделяемый температурный тренд имеет примерно одинаковый период. Тонкая структура температурных вариаций совпадает для многих временн${\text{ы}}'$х интервалов.

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

Для более детальных оценок точности полученной реконструкции, как возрастной шкалы, так и абсолютной величины реконструированной температуры, проведем сравнение с имеющимися региональными гляциологическими [22] и дендроклиматическими [11] данными годового временнóго разрешения. На рис. 6а представлены данные о реконструкции температуры воздуха за последние 750 лет, полученные на основе исследования изотопного состава льда ледника Белуха [22], расположенного рядом с оз. Кучерлинское. Рисунок 6б – наша реконструкция, для удобства сравнения абсолютных величин с данными изотопии льда, пересчитанная в аномалии от среднего за 1250–2001 гг. Рисунок 6в – отклонения средней летней температуры воздуха Алтайского региона по дендро данным, калиброванные по инструментальным данным метеостанции Усть-Кокса [11].

Рис. 6.

Сравнение региональных температурных реконструкций высокого временного разрешения. (а) – реконструкция температуры на основе d18O ледника Белуха. Данные представляют собой аномалии, от среднего за 1250–2001 гг. Показаны отдельные значения (серые), средние за 10 лет и интервал погрешности (черные, пунктирные) [22]; (б) – реконструкция по литолого-геохимическим данным осадков оз. Кучерлинское, пересчитанная как аномалии от среднего за 1250–2001 гг.; (в) – отклонения средней летней температуры воздуха Алтайского региона по дендрохронологическим данным, калиброванные по инструментальным данным метеостанции Усть-Кокса. Представлены 7-летнее скользящее среднее и отклонения от периода 1940–1994 гг. [11].

Дендрохронологическая реконструкция дает летние температуры, поэтому не будем сравнивать абсолютные величины. Больший интерес представляет сравнение ярко выраженных температурных экстремумов на возрастной шкале. Возрастные шкалы, построенные подсчетом годовых слоев (дендрохронология), благодаря отработанным методикам и отсутствию значимых систематических погрешностей являются на сегодняшний день наиболее точными и могут рассматриваться в качестве эталонных при проведении региональных исследований. Графическое сравнение наиболее выраженных экстремумов позволяет оценить отклонения нашей возрастной шкалы от дендро на интервале 0–400 лет назад абсолютными величинами 5–10 лет, но не более 5 отн. %. Приняв дендрохронологическую шкалу за абсолютно точную, считая летние температурные экстремумы совпадающими со среднегодовыми, можем утверждать, что точность подсчета годовых слоев в донных осадках оз. Кучерлинское (варвохронология) позволяет строить возрастную шкалу с относительной погрешностью менее 5 отн. %.

Сравнение с ледниковыми записями по временнóй шкале дает больший разброс, что может быть связано с суммированием погрешностей варвной и ледниковой возрастных моделей. При этом временнáя разница сопоставляемых экстремумов на температурных реконструкциях составляет величины 5–35 лет, или 4–12 отн. %.

Интервал разброса реконструируемой величины в случае ледниковых данных составляет 3°C (от –1 до +2) относительно среднего 1250–2001 гг. Для наших данных разброс величин чуть больше: от –1.5 до +1.8°C, что очень близко с учетом диапазона ошибок обеих реконструкций. Существенное отличие нашей реконструкции от ледовой проявляется в слабо выраженном современном потеплении, хотя, с учетом интервала погрешности, различие может быть не столь значимо.

ЗАКЛЮЧЕНИЕ (ВЫВОДЫ)

1. Озepo Кучерлинское имеет тонкослоистый осадок, содержащий отдельные годовые слои (варвы), что позволяет строить возрастную модель глубина керна–возраст слоя осадка с точностью не хуже 5 отн. %. Подсчет слоев может вестись по фотографиям исходного влажного керна, но необходимо дополнительно контролировать результаты подсчета исследованием оптических шлифов и сканирующим РФА СИ с субмикронным пространственным разрешением.

2. Индикатором границ отдельных годовых слоев в донных осадках оз. Кучерлинское служит отношение Rb/Sr по данным сканирующего РФА СИ.

3. Сравнение временн${\text{ы}}'$х рядов литолого-геохимических данных на годовой шкале с региональными метеонаблюдениями позволяет определить набор геохимических индикаторов климата и построить трансферные функции, связывающие состав датированных слоев донных отложений с метеопараметрами.

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

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

Список литературы

  1. Бабич В.В., Рудая Н.А., Калугин И.А., Дарьин А.В. Опыт комплексного использования геохимических особенностей донных отложений и палинологических записей для палеоклиматических реконструкций (на примере оз. Телецкое, Российский Алтай) // Сиб. экол. журн. 2015. Т. 22. № 4. С. 497–506. https://doi.org/10.15372/SEJ20150401

  2. Дарьин А.В., Александрин М.Ю., Калугин И.А., Соломина О.Н. Связь метеорологических данных с геохимическими характеристиками современных донных осадков оз. Донгуз-Орун, Кавказ // ДАН. 2015. Т. 463. № 5. С. 602. https://doi.org/10.7868/S0869565215230176

  3. Дарьин А.В., Бабич В.В., Калугин И.А., Маркович Т.И., Рогозин Д.Ю., Мейдус А.В., Дарьин Ф.А., Ракшун Я.В., Сороколетов Д.С. Исследование геохимических особенностей годовых слоев в донных осадках пресноводных озер методом рентгенофлуоресцентного микроанализа с возбуждением синхротронным излучением // Изв. РАН. Сер. физ. 2019. Т. 83. № 11. С. 1572–1575. https://doi.org/10.1134/S0367676519110085

  4. Дарьин А.В., Калугин И.А., Бабич В.В., Маркович Т.И., Грачев А.М., Дарьин Ф.А., Ракшун Я.В., Сороколетов Д.С. Поиск годично стратифицированных донных осадков в озерах Горного Алтая методом рентгенофлуоресцентного микроанализа с использованием синхротронного излучения // Изв. РАН. Сер. физ. 2019. Т. 83. № 2. С. 243–246. https://doi.org/10.1134/S0367676519020108

  5. Дарьин А.В., Калугин И.А., Ракшун Я.В. Сканирующий рентгеноспектральный микроанализ образцов донных осадков с использованием синхротронного излучения из накопителя ВЭПП-3 ИЯФ СО РАН // Изв. РАН. Сер. физ. 2013. Т. 77. № 2. С. 204. https://doi.org/10.7868/S0367676513020105

  6. Дарьин А.В., Ракшун Я.В., Сороколетов Д.С., Да рьин Ф.А., Калугин И.А., Максимова Н.В., Маркович Т.И. Исследование сезонного геохимического сигнала в годовых слоях донных осадков оз. Донгуз-Орун методом сканирующего РФА с использованием микрокапиллярной рентгеновской оптики // Изв. РАН. Сер. физ. 2015. Т. 79. № 1. С. 137. https://doi.org/10.7868/S036767651501010X

  7. Дарьин Ф.А., Ракшун Я.В., Сороколетов Д.С., Дарьин А.В., Калугин В.М. Разработка методик микро-РФА на пучках синхротронного излучения из накопителя ВЭПП-3 и их применение для исследования распределения элементов в природных образцах // Ядерная физика и инжиниринг. 2017. Т. 8. № 1. С. 86–90. https://doi.org/10.1134/S2079562917010067

  8. Нарожный Ю.К., Осипов А.В. Ороклиматические условия оледенения Центрального Алтая // Изв. РГО. 1999. Т. 131. Вып. 3. С. 49–57.

  9. Ненашева Г.И. Растительность и климат голоцена межгорных котловин Центрального Алтая: монография. Барнаул: Изд-во Алт. ун-та, 2013. 164 с.

  10. Рихванов Л.П., Окишев П.А., Соболева Н.П., Матаев Е.И. Геохимическая характеристика ленточных глин Горного Алтая и возможности их использования при гляциологических исследованиях // Изв. Томск. политех. ун-та. Инжиниринг георесурсов. 2015. Т. 326. № 2. С. 23–36.

  11. Сыромятина М.В., Москаленко И.Г., Чистяков К.В. Тенденции изменения климата на Алтае на фоне глобальных климатических изменений (по инструментальным и дендрохронологическим данным) // Вестн. Санкт-Петербургского ун-та. Сер. 7. Геология. География. 2010. № 3. С. 82–91. https://doi.org/10.18551/rjoas.2015-07.01

  12. Федак С.И., Туркин Ю.А., Гусев А.И., Шокальский С.П. и др. Государственная геологическая карта Российской Федерации. Масштаб 1 : 1 000 000 (третье поколение). Серия Алтае-Саянская. Л. М-45. Горно-Алтайск. Объяснительная записка. СПб.: Картографическая фабрика ВСЕГЕИ, 2011. 567 с.

  13. Alexandrin M., Dolgova E., Grachev A., Solomina O., Darin A., Kalugin I. Annual sedimentary record from lake Donguz-orun (central Caucasus) constrained by high resolution SR-XRF analysis and its potential for climate reconstructions // Frontiers in Earth Sci. 2018. V. 6. P. 158. https://doi.org/10.3389/feart.2018.00158

  14. Appleby P. The use of 210Pb and 137Cs as tracers in modelling transport processes in lake catchment systems // Stud. in Env. Sci. 1997. V. 68. P. 441–448. https://doi.org/10.1016/S0166-1116(09)70124-4

  15. Brauer A. Annually laminated lake sediments and their paleoclimatic relevance // Fischer H. et al. (Eds.). The Climate in Historical Times. Berlin: Springer, 2004. P. 109–127. https://doi.org/10.1007/978-3-662-10313-5_7

  16. Butz C., Grosjean M., Fischer D., Wunderle S., Tylmann W., Rein B. Hyperspectral imaging spectroscopy: A pro-mising method for the biogeochemical analysis of lake sediments // J. App. Remote Sens. 2015. V. 9. № 1. 096 031. https://doi.org/10.1117/1.JRS.9.096031

  17. Christiansen B., Ljungqvist F. The extra-tropical Nor-thern Hemisphere temperature in the last two millennia: Reconstructions of low-frequency variability // Climate Past. 2012. V. 8. № 2. P. 765–786. https://doi.org/10.5194/cp-8-765-2012

  18. Croudace I., Löwemark L., Tjallingii R., Zolitschka B. Current perspectives on the capabilities of high resolution XRF core scanners // Quat. Int. 2019. V. 514. P. 5–15. https://doi.org/10.1016/j.quaint.2019.04.002

  19. Darin A., Chu G., Maksimov M., Novikov V. Layer counting and isotopic analysis of the recent bottom se-diments of the glacial lake Kucherla (Russia, Gorny Altai) // 19th Int. Multidisciplinary Scientific Geo Conference SGEM 2019. Sofia, 2019. P. 257–264. https://doi.org/10.5593/sgem2019V/4.2/S06.035

  20. Darin F.A., Kalugin I.A., Darin A.V., Rakshun Ya.V. The study internal structure of the annual layers in lake sediments using synchrotron radiation with x-ray focusing optics // Acta Geol. Sin. (Engl. Ed.). 2014. T. 88. № S1. P. 5–6. https://doi.org/10.1111/1755-6724.12265_1

  21. Cook E.R., Krusic P.J., Anchukaitis K.J., Buckley B.M., Nakatsuka T., Sano M. Tree-ring reconstructed summer temperature anomalies for temperate East Asia since 800 C.E. // Clim. Dyn. 2013. V. 41. P. 2957–2972. https://doi.org/10.1007/s00382-012-1611-x

  22. Eichler A., Olivier S., Henderson K., Laube A., Beer J., Papina T., Heinz W., Schwikowski M. Temperature response in the Altai region lags solar forcing // Geophys. Res. Lett. 2009. V. 36. L01808. https://doi.org/10.1029/2008GL035930

  23. Feng S., Yang B., Mairesse A., Gunten L., Li J., Bräuning A., Yang F., Xiao X. Northern Hemisphere temperature reconstruction during the last millennium using multiple annual proxies // Clim. Res. 2013. V. 56. P. 231–244. https://doi.org/10.3354/cr01156

  24. IPCC, 2013: Climate Change 2013: The Physical Science Basis. Contribution of Working Group I to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change / T.F. Stocker, D. Qin, G.-K. Plattner, M. Tignor, S.K. Allen, J. Boschung, A. Nauels, Y. Xia, V. Bex, P.M. Midgley (Eds.). Cambridge, UK; New York, NY, USA: Cambridge Univ. Press, 2013. 1535 p. https://doi.org/10.1017/CBO9781107415324

  25. Klimenko V., Matskovsky V., Dahlmann D. Multi-archive temperature reconstruction of the Russian Arctic for the past two millennia // Geogr. Environ. Sustain. 2014. V. 7. № 1. P. 16–29. https://doi.org/10.24057/2071-9388-2014-7-1-16-29

  26. Mann M., Bradley R., Hughes M. Northern hemisphere temperatures during the past millennium: Inferences, uncertainties, and limitations // Geophys. Res. Lett. 1999. V. 26. P. 759–762. https://doi.org/10.1029/1999GL900070

  27. Mann M., Zhang Z., Rutherford S., Bradley R., Hughes M., Shindell D., Ammann C., Faluvegi G., Ni F. Global signatures and dynamical origins of the Little Ice Age and Medieval Climate Anomaly // Science. 2009. V. 326. № 5957. P. 1256–1260.https://doi.org/10.1126/science.1177303

  28. Moberg A., Sonechkin D., Holmgren K., Datsenko N., Karlén W. Highly variable Northern Hemisphere temperatures reconstructed from low-and high-resolution proxy data // Nature. 2005. V. 433. № 7026. P. 613–617. https://doi.org/10.1038/nature03265

  29. Rothwell R., Croudace I. Micro-XRF Studies of Sediment Cores: A Perspective on Capability and Application in the Environmental Sciences / I. Croudace, R. Rothwell (Eds.). Micro-XRF Studies of Sediment Cores. Developments in Paleoenvironmental Res. 2015. V. 17. https://doi.org/10.1007/978-94-017-9849-5_1

  30. Trachsel M., Kamenik C., Grosjean M., McCarroll D., Moberg A., Brázdil R., Büntgen U., Dobrovolný P., Esper J., Frank D., Friedrich M., Glaser R., Larocque-Tobler I., Nicolussi K., Riemann D. Multi-archive summer temperature reconstruction for the European Alps, AD 1053–1996 // Quat. Sci. Rev. 2012. V. 46. P. 66–79. https://doi.org/10.1016/j.quascirev.2012.04.021

  31. Tylmann W., Zolitschka B. Annually Laminated Lake Sediments−Recent Progress // Quaternary. 2020. V. 3. № 1. P. 5. https://doi.org/10.3390/quat3010005

  32. Yang B., Braeuning A., Johnson K.R., Yafeng S. General characteristics of temperature variation in China during the last two millennia // Geophys. Res. Lett. 2002. V. 29. № 9. P. 1324. https://doi.org/10.1029/2001gl014485

  33. Zolitschka B., Francus P., Ojala A.E., Schimmelmann A. Varves in lake sediments – A review // Quat. Sci. Rev. 2015. V. 117. P. 1–41. https://doi.org/10.1016/j.quascirev.2015.03.019

Дополнительные материалы отсутствуют.