Верификация палеореконструкций климатических условий позднеплейстоценовых летних сезонов, на основе палинологического анализа едомных отложений и ледяных жил

Васильчук Алла Константиновна доктор географических наук ведущий научный сотрудник, Московский государственный университет им. М.В. Ломоносова (МГУ) 119991, Россия, г. Москва, ул. Ленинские Горы, ГСП-1, 1, географический факультет, НИЛ геоэкологии Севера Vasil'chuk Alla Constantinovna Doctor of Geography Leading Research Fellow, Laboratory of Geoecology of the Northern Territories, Faculty of Geography, Lomonosov Moscow State University 119991, Russia, g. Moscow, Leninskie Gory, GSP-1, 1,, geograficheskii fakul'tet, NIL geoekologii Severa alla-vasilch@yandex.ru Другие публикации этого автора

Васильчук Юрий Кириллович ORCID: 0000-0001-5847-5568 доктор геолого-минералогических наук профессор, кафедра геохимии ландшафтов и географии почв, Московский государственный университет им. М.В.Ломоносова 119991, Россия, Г Москва, область, г. Г Москва,, ул. Ленинские Горы, 1, оф. 2009 Vasil'chuk Yurij Kirillovich Doctor of Geology and Mineralogy Professor, the department of Geochemistry of Landscapes and Geography of Soils, M. V. Lomonosov Moscow State University 119991, Russia, G Moscow, oblast', g. G Moscow,, ul. Leninskie Gory, 1, of. 2009 alla-vasilch@yandex.ru

Работа выполнена при поддержке РНФ (проект № 19-17-00126 – полевые исследования) и в рамках фундаментальной госбюджетной темы "Эволюция, современное состояние и прогноз развития береговой зоны Российской Арктики (ГЗ), номер ЦИТИС: 121051100167-1.

Введение

В последнее время палеотемпературные реконструкции на основании палинологических данных становятся все более продвинутыми и требуют все более сложного математического аппарата. Поскольку ископаемая пыльца использовалась в качестве палеоэкологического и палеоклиматологического показателя в течение многих десятилетий ^[1], интенсивные исследования в этом направлении привели к появлению большого количества работ с оценкой сильных и слабых сторон интерпретации ископаемых палиноспектров в качестве палеоклиматического индикатора ^[1-8]. Развивающееся моделирование палеоклиматических условий демонстрирует существенные расхождения для позднего плейстоцена ^[9]. Некоторые из этих расхождений могут быть вызваны тем фактом, что многие палеоклиматические данные отражают определенный сезон, в то время как сигнал часто интерпретируется как среднегодовое значение ^[10], или потому, что даже в регионах, расположенных близко друг к другу, регистрируются существенные климатические градиенты. Часто при моделировании к данным по палеорастительности подходят со значительным осреднением, что не способствует точности результатов моделирования, в частности это касается летних температур.

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

Для арктических территорий важно определить непосредственно границу Арктики. Для интерпретации данных палинологического анализа как представляется автору, больше всего подходит определение этой границы у А.И. Толмачева ^[11], который считает границей Арктики северную границу непрерывного распространения редколесий (рис. 1).

Рис. 1. Граница растительности Арктики по А.И.Толмачеву ^[11]: I – подзона высокоарктических тундр, IА – приледниковый гумидный вариант, IБ – приледниковый криоаридный вариант, I – внеледниковый вариант, II – подзона арктических тундр, III – подзона северных гипоарктических тундр, IV – подзона южных гипоарктических тундр, V – подзона крупных стлаников, VI – океанические аналоги подзон IV и V за пределами собственно тундровой зоны; 1 – границы подзон, 2 – границы абиссальных глубин Арктического бассейна, 3 – границы вариантов подзоны высокоарктических тундр, 4 – граница Гренландского ледникового щита

Севернее этой границы арктические и аркто-альпийские виды составляют значительную часть флоры, а южнее они играют незначительную роль ^[11]. В область Арктики также им включены: области крупных стлаников в бассейнах Анадыря и Пенжины, северная окраина Скандинавии, Исландия, Гренландия. Но здесь в сложении флоры аномально высокая роль принадлежит бореальным, океаническим гипоарктическим альпийским и аркто-альпийским видам (сформированные здесь субфоссильные палиноспектры практически неотличимы от спектров северной тайги). Широкое признание такой границы обусловлено тем, что она тесно связана с важнейшим геоботаническим рубежом – южной границей зоны тундр, или северной границей древесной растительности. На положение северного предела леса воздействует вторжение холодных арктических вод и льдов. Например, длительное стояние льдов в Гудзоновом заливе обусловливает расположение северной границы леса около 50° с.ш. Из-за этого различие между самым южным и самым северным расположением границы леса в северном полушарии составляет 20°. В Евразии северная граница леса, проходя параллельно береговой линии, колеблется от 72°30΄с.ш. в низовьях Лены до 58°55΄с.ш. на побережье Охотского моря. Именно эта граница может быть определена палинологически, как граница существенного падения концентрации пыльцы древесных пород на фоне высокой концентрации пыльцы типичных тундровых таксонов (Artemisia, Ericaceae, Poaceae, Betulasect. Nanae и др.), за исключением области крупных стлаников. На наш взгляд районирование растительного покрова Арктики по А.И.Толмачеву ^[11] лучше всего подходит для определения аналогов при палеореконструкциях, хотя и требует дополнительных данных, в частности, характеристик температуры и влажности, как основных факторов.

Для количественной оценки летних температур и фенологических связей удобно использовать показатель суммы среднесуточных положительных температур воздуха. Для растений средней полосы в качестве показателя используется сумма эффективных температур (сумма температур выше или равных +5⁰С). Для начала продуцирования пыльцы лиственных пород (Populus, Betula, Salix, Acer, Ulmus,) эта сумма составляет 70–80 ⁰С, исключение – Tilia (645– 690 ⁰С), хвойных (Picea, Pinus) – 180–200 ⁰С, – Poaceae – 400 ⁰С, а, например, для Tussilago farfara L. (мать-и-мачехи) всего +5.3⁰С, т.е. всего один день со средней температурой выше 5⁰С ^[12]. Продолжительность вегетационного периода зависит главным образом от географической широты и климатических условий. Вегетационный период - это время, необходимое для развития растения: для однолетнего - от прорастания семян до созревания семени; для многолетнего - от всхода семян или набухания почек до созревания семени; для деревьев и кустарников от начала сокодвижения и распускания почек до опадания листьев. Тундровые растения в основном многолетники, поскольку это позволяет в кратчайшие сроки переходить к активному росту. Для растений тундровой зоны таким показателем, объединяющим длительность сезона вегетации и температурный режим, является сумма положительных температур, поскольку, например, в пределах арктических тундр среднесуточные температуры могут не превышать +5⁰. Эта характеристика является нормативной, она приводится во всех климатических справочниках. Как оказалось, связь между составом палиноспектров и суммами среднесуточных положительных температур воздуха имеет детерминированный характер, хотя выразить это строго математически довольно сложно (рис. 2). На основании позднеплейстоценовых палинореконструкций именно этот показатель можно оценить достаточно адекватно.

Рис. 2. Изотермы суммы среднесуточных положительных температур для криолитозоны России ^[13].

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

В настоящее время реконструкции по палинологическим данным производятся в основном с использованием понятия биом, который является ключевой экологической и биогеографической единицей. Введение понятия биома привело к смещению акцентов при моделировании растительности прошлых эпох акцента с описания чистых моделей на понимание функциональности растительного покрова. Это также сформировало современное понимание концепции биома как динамической биологической сущности, которая имеет глубокие корни в эволюционном прошлом и претерпевающая изменения в настоящем ^{[14, 15, 2]}.

Разделение спорово-пыльцевых диаграмм из Арктики на информационные блоки лучше всего проработано в проекте САРЕ (Circum Polar Arctic Environment - циркумполярная арктическая среда обитания), участники которого подготовили специальную процедуру биомизации для интерпретации палинологических данных в терминах палеотемператур ^[16]. Под биомизацией подразумевается объединение палиноспектров в группы согласно сочетаниям функциональных типов растений ^[17]. Собственно биомизация палинологических данных - это процесс, посредством которого предполагаемая растительность может быть представлена относительно небольшим числом биомов, описывающих основные структурные, функциональные и флористические различия. Для процедуры биомизации в модели BIOME4 отобрано ограниченное количество функциональных типов растений. Для интерпретации палинологических данных севера Евразии мы использовали классификацию биомов ^[18] и сопоставление состава палиноспектров с функциональными типами растений по Н. Биджелоу с соавторами ^[19] (табл. 1). Очевидно, что поскольку выделенные биомы не отражают всех компонентов палиноспектра, то коррекция значения суммы температур производилась с учетом состава конкретного палиноспектра выделенной палинозоны, с учетом классификации растительного покрова Арктики, предложенной А.И.Толмачевым ^[11].

Таблица 1. Классификация циркумполярных биомов (по Kaplan et al. ^[18])

Эта простая классификация позволяет выделить реакцию растительности на изменения летней температуры и влажности.

Таблица 2. Сопоставление классификации арктических биомов ^{[16, 18,19]}, палинологическая характеристика выделенных биомов, соответствующие суммы летних температур (по данным авторов)

Наши исследования по изучению соотношения длительности вегетационного сезона, а именно изучение массива субфоссильных спектров из ряда регионов Арктики ^[20,21], позволили включить в характеристику соотношения биомов и палиноспектров сумму положительных температур в качестве синтетического показателя, отражающего и средние летние температуры и длительность летнего сезона (табл. 2). Сумма эффективных температур менее удобна, поскольку для некоторых районов температура летнего сезона не поднималась выше +5°С.

Применение данной методики реконструкции суммы положительных температур позволяет определить относительные изменения палиноспектров в отдельных образцах и выразить их в единых терминах суммы положительных температур. Следовательно, предложенная интерпретация позволяет сопоставлять между собой датированные спорово-пыльцевые диаграммы разрезов, расположенных на значительном удалении друг от друга. Разумеется, реконструкции климатических показателей с помощью переходных функций и статистических показателей ^[21-24], учитывающих соотношение каждого палинотаксона с температурами июля и среднегодовыми температурами позволяет реконструировать, палеотемпературы более точно, особенно для голоцена. Преимущество предложенной здесь методики реконструкций сумм положительных температур позволяет точнее интерпретировать опубликованные ранее спорово-пыльцевые диаграммы по разрезам едомных толщ, по которым информация представлена исключительно в графическом виде, и отследить резкие изменения в составе палиноспектров интерпретируя их в терминах суммы положительных температур (рис. 3).

Рис. 3. Корреляция палеотемпературных кривых сумм положительных температур, полученных в результате авторской интерпретации датированных по ¹⁴С спорово-пыльцевых диаграмм по разрезам: Мыс Саблера - А.А.Андреева и Г.Н.Шиловой (Andreev et al. ^[25,26]); Мус-Хая Н.О.Рыбаковой и Л.Л.Козачихиной ^[27], Томирдиаро, Чёрненький ^[28]); Мамонтова Хайота - А.А.Андреева ^[29]; Воронцовский Яр - Н.С.Болиховской ^[30]; по Ойгосский Яр и Мыс Чукочий - Л.Л. Козачихиной ^[28]; по Алёшкинской террасе - И.Н. Федоровой и А.В. Ложкина ^[28]; по Эликчан-4 - А.В.Ложкина, П.Андерсон и Л.Брубейкер ^[31]; по разрезам Гыда, Сеяха-II и III, Монгаталянг, Зелёный Мыс, Плахинский Яр и Дуванный Яр – А.К. Васильчук ^[31-36]. Колебания содержания изотопов кислорода в керне NorthGRIP и привязка событий Хайнриха и Дансгора – Эшгера, согласно данным В.Дансгора, Х.Хайнриха, Д.Даль-Енсен, Дж. Бонда, Дж. Доудсвелла, Л.Хиннов ^{[37, 38, 39, 40, 41, 42]}

В качестве стратотипа мы использовали изменение изотопного состава ледникового керна NGRIP Гренландского ледникового щита, интерпретируемое как в большой мере обусловленное изменением среднегодовых температур. Вполне естественно предположить, что глобальные события отражались и на составе растительного покрова Северной Евразии и, следовательно, могут быть запечатлены на спорово-пыльцевых диаграммах. Изучение разрезов криолитозоны с учетом особенностей формирования многолетнемерзлых отложений позволило привязать имеющиеся данные к ритмам глобальных изменений среднего масштаба – событиям Хайнриха (см. рис. 2). Сопоставление кривых сумм положительных температур между собой выявило ряд общих трендов. Во-первых отмечается близкое сходство полученных распределений на региональном уровне. Так, кривые сумм положительных температур Сеяха 3 и Монгаталянг совпадают почти по всем ключевым точкам, как минимумам, так и максимумам. Совпадения по максимумам и минимумам отмечается также для более удаленных друг от друга разрезов: группы западно-сибирских разрезов - Сеяха 3 и Сеяха 2 Монгаталянг- и таймырского Мыс Саблера. Следовательно, можно говорить, что условия летних сезонов в данном временном интервале изменялись практически синхронно на севере Западной Сибири и Таймыре. Палеотемпературные кривые Зеленого Мыса, Дуванного Яра и Плахинского Яра также имеют ряд общих черт. Сходство между кривыми, полученными по данным из разрезов в долине р.Колымы выше, чем с кривыми полученными в едомных разрезах на севере Западной Сибири и Таймыра. Отметим очень близкую палинологическую характеристику разрезов Плахинского Яра и Алешкинской террасы и практически полное совпадение на одинаково датированных фрагментах колебаний сумм положительных температур.

Имеются также общие моменты, зафиксированные почти на всех кривых сумм положительных температур из проинтерпретированных спорово-пыльцевых диаграмм датированных по радиоуглероду разрезов синкриогенных отложений с повторно-жильными льдами. На полученных кривых заметны минимумы, соответствующие похолоданиям, связанным с событиями Хайнриха, это означает, что в эти периоды летние сезоны отличались низкими температурами, а их длительность также снижалась (рис.4). Сумма положительных температур для этих периодов не превышала 300 град. х сут. Распределение сумм положительных температур для стадиалов событий Хайнриха практически не отличалось. В то время как интерстадиальный растительный покров различался по доминирующим компонентам и распределение сумм положительных температур для каждого периода летнего потепления было более индивидуальным.

Рис. 4. Распределение значений сумм положительных температур в минимальную (стадиальную) фазу третьего события Хайнриха 28-27 т.л.н.

Возможно, снижение температур вегетационного периода было связано с увеличением площади ледового покрова в летние сезоны. Например, выделяется минимум летних температур, соответствующий третьему событию Хайнриха примерно 28,5-27 тыс. лет назад. Интересно, что, наиболее заметное увеличение сумм положительных температур, за которым следует резкое их снижение отмечается в начале каждого события Хайнриха.

Исследованиями М Стайвера и П.Грутса ^[43] установлено, что серджи, связанные с событиями Хайнриха, происходили в течение холодного пика отдельных событий Дансгора-Эшгера, выделяемых на изотопных диаграммах гренландских ледниковых кернов, средняя продолжительность цикла Дансгора–Эшгера составляла около 1460 лет, при этом холодная фаза длилась около 600 лет и заканчивалась резким переходом к теплой фазе. То есть, судя по реакции растительного покрова в Северной Евразии и данным изотопного состава ледяных жил ^[44], наиболее выраженным было летнее потепление, которое прерывалось резким похолоданием, связанным с поступлением айсбергового материала в моря Северного Ледовитого океана.

Общий позитивный пик отмечен около 29 тыс. лет практически на всех кривых положительных температур. Он соответствует интерстадиалу третьего события Хайнриха. На кривых Сеяха-2 и Зеленый Мыс можно выделить существенное снижение температур вегетационного периода, которое может соответствовать первому событию Хайнриха, а затем последующий резкий рост, который может быть сопоставлен с аллередом. Это потепление на севере Западной Сибири проявилось несколько ярче с большей амплитудой. На кривой Гыда удалось выделить период снижения сумм положительных температур, сопоставляемый с поздним дриасом, и последующий переход к гренландскому периоду голоцена. Можно утверждать, что глобальные климатические колебания, обусловленные изменениями ледовой обстановки в Северной Атлантике – событиями Хайнриха особенно заметно отражались на растительном покрове современной криолитозоны, т.е климатические характеристики летнего сезона изменялись весьма существенно. Разумеется, каждая наземная система по-разному реагирует на климатические изменения, что приводит к различному времени отклика и усложняет хронологию. По результатам проведения конференции «События Хайнриха в наземных архивах” на конгрессе INQUA в Дублине в 2019 году, события Хайнриха отразились на составе ряда наземных объектов: эоловых отложений ^[45], в Ауэль-Маара в Германии, на составе растительности и остракод в оз. Петен-Ица в Гватемале, а также на изотопном составе спелеотем. С. Вететрих и Н. Рудая с ^[46], на основании изучения пыльцы и спор в разрезе Собо-Сисе и сопоставления их с синхронно формировавшимися едомами пришли к заключению, что региональные пыльцевые записи в дельте Лены и на Быковском полуострове подтверждают существование позднеплейстоценовой тундро-степи в Западной Берингии во время MIS 3 с частично более теплым, чем сегодня, летом и более влажными условиями. Наиболее благоприятные условия отмечаются для интерстадиала после похолодания пятого события Хайнриха, что совпадает с максимумом суммы положительных температур для кривой Зеленый Мыс (700-750 град.х дней.). Согласно их выводам С Веттериха ^[46], во время MIS 2 лето, вероятно, было более холодным, чем в течение MIS 3, и более сухим с нестабильными условиями увлажнения в полигональных ландшафтах. Согласно реконструкциям суммы положительных температур (рис. 2) в течение интерстадиалов 2, 3 и 4 событий Дансгора-Эшгера (стадии MIS2) сумма положительных температур была на 100-150 град.хдней ниже чем в течение интерстадиалов с 5 по 15 событий Дансгора –Эшгера.

По результатам реконструкций состава палиноспектров из отложений оз. Лирней ^[47] минимальное значение среднеиюльской температуры 7.7°C установлено 27.90 тыс. кал. лет назад, а максимальное 9.3°C, 15.72–15.47 тыс. кал. лет назад. На основании результатов моделирования с использованием данных палинологического анализа и изменений термогалинной циркуляции поверхностных и глубинных вод Северной Атлантики ^[48] получен вывод, что интерстадиалы событий Дансгора-Эшгера в Западной Европе характеризовались вдвое большей суммой осадков по сравнению со стадиалами, температура самого холодного месяца была на 17°C выше, а самого теплого месяца на 8°C выше, реконструированная сумма положительных температур составила около 1300 °C х дней. Cовокупный эффект низких температур в зимний период, недостаточное увлажнение летом и низкое значение суммы эффективных температур, полученное по результатам моделирования, объясняет по мнению авторов отсутствие лесов к северу от 48°с.ш. во время межстадиала 14 события Дансгора-Эшгера ^[48] несмотря на относительно благоприятные условия вегетационного сезона. На наш взгляд, причиной этому не недостаточная сумма эффективных температур, а прежде всего короткий сезон вегетации. Поэтому в разрезах позднеплейстоценовой палеокриолитозоны в Европе встречается пыльца травянистых растений, требовательных к тепловому режиму и в настоящее время произрастающих в лесной зоне ^[49], а пыльца древесных пород отмечается в небольших количествах, не соответствующих лесным формациям. По мнению авторов данной статьи, при построении моделей для реконструкции растительного покрова и температуры сезона вегетации и длительность сезона вегетации или не учитываются, или учитываются недостаточно.

В обзорной статье М.Шевалье с соавторами ^[2], дана оценка 75 летнему опыту разработки и применения количественных методов реконструкции климата на основании палинологических данных, начиная с пионерных реконструкции Й. Иверсена ^[50]. Несмотря на успешную разработку альтернативных климатических индикаторов в последние годы, в том числе основанных на изотопах и биомаркерах, ископаемая пыльца, по мнению авторов обзора ^[2], остается наиболее важным источником количественных климатических данных для наземных неледниковых регионов, например ^[51]. С развитием новейших статистических способов интерпретации, количественные реконструкции стали намного доступнее.

Выводы

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

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

Сопоставление кривых сумм летних температур воздуха между собой выявило ряд общих трендов: а) близкое сходство полученных распределений на региональном уровне; б) синхронность изменений условий летних сезонов на севере Западной Сибири и Таймыре; в) сходство диаграмм сумм летних температур воздуха, полученных по данным опорных позднеплейстоценовых разрезов в долине р. Колымы.

Установлено хронологическое соответствие минимумов сумм летних температур воздуха на полученных диаграммах и событий Хайнриха.