CC BY-NC
15Арктика и Антарктика
Правильная ссылка на статью:
Васильчук А.К., Буданцева Н.А., Суркова Г.В., Чижова Ю.Н. — О надежности палеотемпературно-изотопных уравнений Васильчука и становлении изотопной палеогеокриологии // Арктика и Антарктика. - 2021. - № 2. DOI: 10.7256/2453-8922.2021.2.36145 URL: https://nbpublish.com'llbrary_read_article.php?id=36145
О надежности палеотемпературно-изотопных уравнений Васильчука и становлении изотопной палеогеокриологии
Васильчук Алла Константиновна
доктор географических наук
ведущий научный сотрудник, Московский государственный университет им. МВ. Ломоносова (МГУ)
119991, Россия, г. Москва, ул. Ленинские Горы, ГСП-1, 1, географический факультет, НИЛ геоэкологии
Севера
И alla-vasilch@yandex.ru Буданцева Надежда Аркадьевна
кандидат географических наук
старший научный сотрудник, Московский государственный университет имени МВЛомоносова, географический факультет, кафедра геохимии ландшафтов и географии почв
119991, Россия, г. Moscow, ул. Ленинские Горы, 1, оф 2007
И nadin.budanceva@mail.ru Суркова Галина Вячеславовна
доктор географических наук
профессор, кафедра метеорологии и климатологии, Московский государственный университет имени
М.В. Ломоносова, географический факультет
119991, Россия, г. Moscow, ул. Ленинские Горы, 1, оф 2001 И galina_surkova@mail.ru
Чижова Юлия Николаевна
кандидат географических наук
старший научный сотрудник, Институт геологии рудных месторождений, петрографии, минералогии и геохимии Российской академии наук (ИГЕМ РАН), Лаборатория изотопной геохимии и геохронологии
119017, Россия, г. Москва, ул. Старомонетный Пер., 35 И eacentr@yandex.ru
Статья из рубрики "Многолетнемерзлые породы и подземные льды Арктики, Антарктики и горных регионов" DOI:
10.7256/2453-8922.2021.2.36145
Дата направления статьи в редакцию:
22-07-2021
Аннотация: Статья посвящена оценке вклада проф. Васильчука в развитие изотопной геокриологии, а также тестированию надежности палеогеокриологических реконструкций, основанных на изучении изотопного состава полигонально-жильных льдов. Его открытие прямой зависимости между изотопным составом ростков современных ледяных жил и температурными характеристиками зимнего сезона в 1989 г положило начало перспективному изучению голоценовых и позднеплейстоценовых сингенетических повторно-жильных льдов, как надежного палеоклиматического архива. Впервые им были получены характеристики зимнего периода для позднего плейстоцена и голоцена, и по результатам исследования созданы карты распределения палеотемператур для ключевых периодов позднего плейстоцена. Полученные позже разными исследователями данные по изотопному составу повторно-жильных льдов все без исключения вписываются в распределение палеотемпературы в пределах криолитозоны, реконструированной проф. Васильчуком. Для установления степени надежности реконструкций палеотемпературы протестированы соотношения, предложенные Ю .К.Васильчуком и уравнения регрессии из работ Коняхина, Майераи Облогова. По результатам тестирования установлено, что значения среднезимней и среднеянварской температуры при реконструкции по соотношению Васильчука всегда попадают в доверительный интервал во всем современном температурном диапазоне. Тестирование современной температуры с использованием уравнений регрессии часто обнаруживают существенное отклонение (с ошибкой нередко более 3-4°С) от реальных значений среднезимних и среднеянварских температур.
Ключевые слова: повторно-жильные льды, поздний плейстоцен, голоцен, изотопы кислорода, па леотемпература , тестирование уравнений регрессии, изотопная палеогеокриология, модель гетероциклического развития, стратегия датирования едомы, позднеплейстоценовый криохрон
Работа выполнена в рамках фундаментальной госбюджетной темы "Эволюция, современное состояние и прогноз развития береговой зоны Российской Арктики (ГЗ), номер ЦИТИС: 121051100167-1.
Введение
30 лет назад Ю.К. Васильчуком в его докторской диссертации и вышедшей следом монографии была опубликована концепция палеоклиматических реконструкций по сингенетическим повторно-жильным льдам благодаря чему сформировалось целое научное направление - изотопная геокриология. Впервые он получил датированные по радиоуглероду детальные изотопные диаграммы полигонально-жильных комплексов 1982
г.-Ш Всего c 1975 по 2021 гг. он лично опробовал и изучил более 50 опорных разрезов от Воркуты до Уэлена, что обеспечило возможность прямого хронологического сопоставления изотопных записей подземных и наземных льдов и их надёжную возрастную привязку для периода более чем в 40 тыс. лет.
В этой статье авторы хотели бы оценить вклад проф. Васильчука в развитие изотопной геокриологии, а также подтвердить надежность палеогеокриологических реконструкций, основанных на изучении изотопного состава полигонально-жильных (повторно-жильных) льдов (ПЖЛ).
Ооответствие среднезимней и среднеянварской температуры воздуха изотопному
составу кислорода жильных ростков
Впервые понятие о сингенетическом формировании клиновидных льдов в геокриологическую терминологию ввел Г. Гальвиц в 1949 г.^2^. После того, как А.И. Попов[3,4] сформулировал гипотезу о сингенетическом накоплении мощных ледяных жил в аллювиальных отложениях, а Б.Н. Достовалов^ рассчитал возраст ледяных жил, интерес к ним сильно вырос. Отложения с повторно-жильными льдами активно изучались, однако, по мере накопления данных из-за разочарования в возможностях прочтения палеогеографической информации, заключенной в ПЖЛ в конце 70-х и в 80-е годы ХХ века интерес к изучению сингенетических повторно-жильных льдов заметно
угас. И даже очень яркая монография Н.Н. Романовскогс^, изданная по материалам его докторской диссертации, не вернула интерес к сингенетическим ПЖЛ. Определенным ренессансом в этой области можно считать открытие прямой зависимости между изотопным составом ростков современных ледяных жил и среднеянварской температурой воздуха, которое сделал Ю.К. Васильчук в период с 1981 по 1989 гг. Это открыло широкую перспективу для изучения голоценовых и позднеплейстоценовых сингенетических ПЖЛ, как надежного палеоклиматического архива. В отличие от ледниковых кернов сингенетические повторно-жильные льды содержат исключительно зимний палеотемпературный сигнал. Изучение стабильных изотопов в повторно-жильных
льдах Ю.К. Васильчук начал в 1981 г. когда получил первые данные по содержанию 18О в повторно-жильных льдах едомы в устье р. Сеяха с помощью В.А. Полякова и Л.Н.
Крицук на масс-спектрометре ВСЕГИНГЕО. Первая изотопная диаграмм^ из 9 образцов была получена в сентябре 1982 г. при содействии Л.Д. Сулержицкого и А.Д. Есикова. Измерения велись на масс-спектрометре Института Водных проблем АН СССР.
Одновременно по этому разрезу был получен ряд 14С датировок в радиоуглеродной лаборатории ГИН АН СССР. Эта первая датированная изотопная диаграмма по
сингенетическим повторно-жильным льдам была опубликована спустя 2 года в ДАНШ. В дальнейшем Ю.К. Васильчуком была развернута работа по изотопному исследованию сингенетических повторно-жильных льдов в других районах криолитозоны. Все образцы льда из ледяных жил для изотопного анализа и органического материала для датирования отбирались им собственноручно. При отборе образцов из мощных жил всегда использовалось правило: детальный отбор по вертикали не реже чем через 1 м (позже интервал был сокращен до 0,5 м и даже до 0,1 м), и, одновременно, детальный отбор по горизонтали, там, где имеются признаки сингенеза, например, выраженные «плечики». В период 1982-1988 гг. было опробовано несколько десятков разрезов в Западной Сибири, Якутии, Забайкалье и на Чукотке: Гыда (1983 г.), Матюй-Сале (1983 г.) Зеленый Мыс (1983-1988 гг.), Плахинский Яр, Омолон и Дуванный Яр (1984 г.), Мамонтова Гора и Быковский п-ов (1985 г.), Анадырь, о.Айон и Кулар (1986 г.), оз.Коолень, Ледовый обрыв и У сть-Алганский разрезы на р.Майн (1987 г.), Чара (1988 г.) и др. Выполнение анализа стабильных изотопов кислорода производилось на современных масс-спектрометрах Института водных проблем РАН, Института геологии (Таллин), Хельсинского университета, геологической службы Нижней Саксонии (Ганновер), лаборатории "Арсенал" (Вена). Помимо детального отбора позднеплейстоценовых и голоценовых ПЖЛ Ю.К. Васильчуком одновременно решалось несколько принципиальных задач. Во-первых, была теоретически оценена скорость диффузии и обоснована сохранность изотопного сигнала в течение сотен тысяч лет в
повторно-жильных льдах [7,8]. Во-вторых, велся поиск основы для интерпретации изотопных данных в палеоклиматическом аспекте. Для сравнения во всех точках опробования и, особенно вблизи полярных метеостанций, исследовались ростки современных жильных ростков на севере Западной Сибири, в Якутии, Забайкалье,
Магаданской области и на Чукотке. Уже в середине 1988 г. основная часть программы по полевому опробованию современных ростков, голоценовых и плейстоценовых ледяных жил на огромной территории от Ямала до Восточной Чукотки была закончена. До настоящего времени эта коллекция является самым представительным массивом данных по изотопному составу жильных ростков для криолитозоны России. Впоследствии эти данные были использованы в качестве основного ядра изотопно-температурных
расчетов М.А.Коняхиным^, В.А.Николаевым и Д.В.Михалевы Х. Майером^11^,
И.Д.Стрелецкой-^12!, Г.В. Облоговым^13!. Надо отметить, что в уравнениях Ю.К. Васильчука также учтены результаты изотопных исследований, выполненных в конце 80-х годов А.А. Архангеловым, М.А. Коняхиным, Д.В. Михалевым.
Первые 2-3 года интерпретация изотопных данных производилась Ю.К. Васильчуком по уравнениям В. Дансгора^14^ однако уже к 1985 г. ему стало ясно, что необходимо получение физически более обоснованных уравнений применительно к изучаемому объекту, поскольку у В. Дансгора изотопные данные преобразуются в среднегодовые температуры, а для ПЖЛ температуры и осадки летнего сезона незначимы. Поэтому для сравнения с изотопным составом кислорода в ПЖЛ были выбраны температурные
характеристики зимнего сезона: сумма среднесуточных температур ниже 0оС - типичная метеорологическая характеристика, приводимая во всех климатических справочниках, (количество суток в справочниках также приводится); среднезимняя температура была
рассчитана как отношение суммы среднесуточных температур ниже 0оС, к количеству дней с отрицательной температурой- т.е. расчет среднезимней температуры относится не к календарной зиме, а к температурной; в качестве еще одной важной характеристики была выбрана среднеянварская температура, как средняя за период наблюдений температура самого холодного месяца (в некоторых точках, где февральская температура ниже, то не января, а февраля). Для этого использовались данные по 250
станциям из Справочника по климату СССГ^1^ . Для анализа были выбраны температурные ряды тех метеостанций, которые находились ближе всего к точкам отбора современных жильных ростков. Ю .К.Васильчуком было установлено эмпирическое соответствие среднезимней и среднеянварской температуры изотопному составу кислорода жильных ростков для всех регионов, в которых был исследован изотопный состав жильных ростков. Поскольку жильные ростки формируются примерно 50-100 лет, в 80-е годы это совпадало с периодом функционирования полярных станций, основанных в начале 30-х годов ХХ века и метеонаблюдений в наиболее крупных городах, расположенных в криолитозоне - Якутск, Салехард, и др., т.е. анализировались температурные ряды метеостанций длительностью 100-50 лет, реже анализировались температурные ряды длительностью 30 лет. Ю .К. Васильчуком были проанализированы эти 3 параметра для 250 метеостанций, расположенных на территории криолитозоны
России. К 1989 г им была построена карта распределения величин б18О в современных жильных ростках, формировавшихся приблизительно 50-100 лет, выполнено сравнение
значений б18О с зимними температурами приземного слоя воздуха, рассчитаны эмпирические уравнения, связывающие среднеянварскую, среднезимнюю, и сумму среднезимних температур с изотопным составом современных жильных ростков. Эти результаты вместе с полученными к этому времени палеотемпературными реконструкциями для ключевых этапов голоцена и позднего плейстоцена были доложены Ю.К. Васильчуком на III конференции «Изотопы в гидросфере» в Каунасе в 1989 г^16 и опубликованы в журнале «Водные ресурсы» в 1990 г. (версия на русском)!-17!, и в 1991 (версия на английском)^18!. В 1991 г. эти материалы были обобщены в докторской
диссертации «Позднечетвертичные синкриогенные толщи севера Евразии: Строение,
изотопно-кислородный состав и условия формирования»-19!, которая была блестяще защищена в Якутске в Институте мерзлотоведения 29 октября 1991 г. Оппонентами были проф. Н.Н. Романовский, проф. Б.И. Втюрин и проф. С.М.Фотиев. По результатам выполненной работы было сформировано новое научное направление - изотопная палеогеокриология. Установлена длительная сохранность изотопно-кислородного состава в подземных льдах на протяжении десятков тысяч лет, обусловленная относительно незначительным обменом изотопов кислорода под действием самодиффузии-7!.
Палеотемпературные карты для ключевых периодов позднеплейстоценового
криохрона
Заключительный этап позднего плейстоцена - период 40-10 тыс. лет назад - был вперые выделен Ю .К.Васильчуком как единый криохрон с очень суровыми зимами, внутри которого осцилляции температур воздуха и мёрзлых толщ были незначительны. Был определен тренд распределения изотопа кислорода-18 в ледяных жилах различного возраста на севере Евразии. Было установлено, что в позднеплейстоценовых повторно-жильных льдах тренд изменений значений б18О был подобен современному: величина
б18О уменьшается при продвижении с запада на восток на 8-10%о, от -19 до -25%о в западно-сибирских полигонально-жильных системах до -30, -35% в северо-якутских, а затем вновь возрастает на 2-3% : до -28, -33% на севере Чукотки и на 6-8% , до -23, -29% на востоке Чукотки-19!; в голоценовых повторно-жильных льдах сохранилась та же тенденция распределения значений б18О - уменьшение на 6-8% на восток: от -14, -20% в западно-сибирских полигонально-жильных системах до -23, -28% в североякутских, однако существенно положительнее стали значения б18О в ледяных жилах на
Чукотке: до -15, -21%[19]. В результате было определено, что характер воздушного переноса на большей части евроазиатской Субарктики в конце позднего плейстоцена
был подобен современному, преобладал западный перенос воздушных масс-19!. Влияние Атлантики (на всем пространстве от Ямала до северной Якутии) было существенным, хотя и несколько слабее современного. В самых восточных районах (Чукотка) влияние тихоокеанских воздушных масс было заметно меньше современного. Было установлено, что сопоставление изотопного состава повторно-жильных льдов со среднегодовыми температурами не соответствует процессу формирования повторно-жильных льдов т.к. весь летний сезон морозобойные трещины закрыты или залиты попавшей в них замерзшей водой. Разработан новый сценарий эволюции температуры воздуха и мерзлых толщ в верхних горизонтах криолитосферы в позднечетвертичное время на основании данных, полученных при использовании сопряжённого изотопно-палеогеокриологического анализа: - массовые данные палинологических (более 2000), радиоуглеродных (около 400) и изотопно-кислородных определений (более 1500) послужили основой для составления новой шкалы палеоклиматических изменений на
севере Евразии для последних 40 тысяч лет-19!. Ю.К. Васильчук довел исследования до очень важного результата - до стадии реконструкции палеотемператур на ключевые периоды позднего плейстоцена (рис. 1). Отметим, что полученные позже данные разных авторов[20-30 и др. ] по изотопному составу повторно-жильных льдов все без исключения вписываются в распределение палеотемпературы в пределах криолитозоны, впервые полученное Ю.К.Васильчуком-19!.
а
40 60 80 100 120 140 160 180
б
40 60 80 100 120 140 160 180
Рис. 1. Карта распределения реконструированной среднеянварской температуры для ключевых периодов позднеплейстоценового криохрона: а - 30-28 тыс. лет назад; б - 2422 тыс. лет назад. Из учебника Ю.К. Васильчука и В.М. КотляковЦифры на карте: в знаменателе - величина б18О, в числителе - реконструированная температура
Результаты этих исследований были столь новы, что потребовалась длительная экспертиза на кафедре мерзлотоведения с участием приглашенных специалистов. Работа Ю.К. Васильчука была поддержана ведущими учеными криолитологами и геокриологами: И.Д. Даниловым, Т.Н. Каплиной, К.А. Кондратьевой и другими выступавшими. На защите диссертации также возникла дискуссия, был отмечен интерес к работе не только крупнейших советских ученых-изотопистов и геокриологов: Я-М Пуннинга, А.В. Раукаса, В.И. Соломатина, В.В. Рогова, С.М. Фотиева, В.И. Втюрина, Н.Н. Романовского, Т.Н Каплиной, но и мэтров зарубежной геокриологии Дж. Росс Маккая, А.Л. Уошборна, Э. Костера, которые прислали положительные отзывы, где отмечалось, что, к сожалению, для криолитозоны в Канаде и США на тот период подобных работ не проводилось. К сожалению, их нет и до сих пор, хотя работы в этом направлении ведутся.
Обсуждение изотопно-температурных соотношений с Росс Маккаем, Уошборном и
Дансгором
Выдающийся канадский геокриолог, профессор Росс Маккай из Университета Британской Колумбии в отзыве на автореферат докторской диссертации Ю.К. Васильчука (рис. 2, а) написал: "Я читал английский вариант автореферата Ю.К. Васильчука. По моему мнению, применение изотопного состава для расчета геотемператур для последних 40 тысяч лет очень важно и в работе продемонстрировано очень убедительно. Я бы хотел подчеркнуть два факта. Первое: хотя я и хорошо знаком с литературой по изотопии, я не являюсь авторитетом в этой области. Второе, географически область моих исследований - Канадская Арктика, а не Северная Евразия. По-моему, рисунок 1, который суммирует огромное количество данных (рис. 3) показывает очень близкую взаимосвязь между
б18О в ледяных жилах и зимними температурами очень интересен. Так же можно прокомментировать другие рисунки и таблицы. Расчет прошлых зимних температур по данным изотопного анализа является наиболее значительным вкладом автора. Сопоставимых данных или исследований такого направления по Северной Америке не существует. Вкратце резюмируя: на меня произвели впечатление и вклад автора и выполненные им исследования. Я надеюсь вскоре увидеть эту диссертацию опубликованной".
Выдающийся американский геокриолог, профессор А.Л. Уошборн из Центра четвертичных исследований Вашингтонского Университета в отзыве на автореферат докторской диссертации Ю.К. Васильчука (рис. 2, б) написал следующее: "Я считаю, что исследования региональных колебаний изотопов кислорода в ледяных жилах потенциально очень перспективны в геокриологии. Хотя я не геохимик, но я в состоянии оценить методологические аспекты обширной работы Ю.К.Васильчука по Северной Евразии. Я уверен, что результаты этого основательного исследования охватывающего обширнейший регион будут важны не только для геокриологии, но также внесут значительный вклад в изучение четвертичного периода в целом. Они также имеют особо важное значение для международных программ связанных с глобальными климатическими изменениями. Я желаю доктору Васильчуку успехов".
TfWUMOWJAI. 011 Tu.K. VtiBlI'clHih preprint, fflf 1
deposit« dî North ot Bnvnlni ^tjructur«, and formation conditions", pr«*"wiit-®d s« dagree of Üoctor Science in gavlc/yr.
I.ntn Qnat^niRtt ejficryeeniUe
pifgiMt 1яо(»р* eompoeltiitrt оопПгмлИо» of Sclunc*
I have read the English abstract of the thesis by Yu. K. Vastlr clink. 10 Illy opinion, his use of isotopes in estimating geotemperatures for the last forty thousand scars is met important ami viry convincing, However» I would Like to stress two facts, Plrat, although 1 am reasonably familiar with tli* Isotope Literature, i am not an authority oil the subject. Second» iny geographic area of familiarity is the Canadian Apetii, not northern Eurasia- To as, Figure 1, tihlflh RUiiiMariaeR a great deal of data, showa a very close relationship» and an interesting one, between 0" in ice wedges and winter temperabures, A aimLiar type of eomm-Ht could be im-le for Che other figures and tables, "Hie estimation of past winter temperatures fnn lavtopic itnaljfft-i» is a most Blgnirlcani con tri bub ion. 'il nii'e are HO comparable data or studies in Worth Aiiwrlca, ill summary, I am Lmprenaeil with the research ami contributions uf the author, i hope to see, at flow* future 'Irtle, the published them In•
g Testimonial
on Yu.K.Vasil'chuk's preprint of the Thesis "Late Quaternary syncryogenlc deposits ol North of Eurasia: structure, oxygen Isotope composition and formation conditions ", presented on confirmation ol Science degree 01 Doctor science in geology
I regard research on regional variations in oxygen-isotope ratios 1n ice wedges of comparable age as constituting a potentially very rewarding contribution to geocryology, Although 1 am not a geochemlst and in a position to evaluate the methodoliglcal details of Dr. Vasil'chuk's extensive work on this subject in Northern Eurasia, I m confident that the results of sound research encompassing this extensive region will lie not only highly important for geocryology but would also contribute significantly to both Quaternary studies and the present international emphasis on programs involving climate and global change.
I wish Dr. Vasil 'chuk well.
s тл. u P of til«
"К^ЪоюМккзу Dtpaurncnt ol GestfMpny TI» University Ы ENilsu OolumbU Vaniouvf, ЗЛ, V6T 122
jj.ßtti/nlfufe^-
S I G N A I II J S Professor, Dr. J-Itosa Паска/ 1Л (Kirim*tit or &*agriipj)j\ Uni orally of British CelunbiR Canada
Ояотагу Besseret, Сет«., akjo
Stamp oi the olilce
sàifiïîra
Frolessor A. L.Washburn, 5
Quaternary Research Center University oi Washington Seattle, UEi
Рис. 2. Отзывы проф. Росс Маккая (а) и проф. А. Уошборна (б) на автореферат диссертации Ю. К. Васильчука "Позднечетвертичные синкриогенные толщи севера
Евразии: строение, изотопно-кислородный состав и условия формирования"^32!
Рис» I. Сопоставление характера вариаций содержания кислорода-18 в современных ростках ледяных жж [формировавшихся на территории Евроазиатской криолитозоны в последние 1СЮ лет) и суш зимних температур воздуха, осреднённых за период последних 60-100 лет: изолинии: I - уверенно проводимые, 2 - предполагаемые. 3 - южная граница современного активного роста ледяных жил
Именно сопоставление характера изотопных диаграмм в мощных повторно-жильных льдах и в текстурных льдах из вмещающих их син-криогенных пород послужило нам основанием для гипотезы о циклическом характере формирования мощных ледяных жил. Сущность этой гипотезы отражена на рис. 2.
16
Рис. 3. Рисунок 1 из автореферата докторской диссертации Ю.К. Васильчука-32! (1991), на который ссылается проф. Дж. Росс Маккай
С профессором Дж. Россом Маккаем Ю.К. Васильчук встречался лично, во время его визита в Москву. С ним удалось обсудить проблему отбора образцов из повторно-жильных льдов в связи с процессами растрескивания и заполнения трещин. Профессор Дж. Росс Маккай изучал процессы растрескивания в течение многих полевых сезонов и поэтому обсуждение данного вопроса его заинтересовало. В результате дискуссии вертикальный отбор повторно-жильных льдов был признан гораздо более
информативным по сравнению с горизонтальным в целях изотопной индикации палеотемпературных условий, принимая во внимание, что растрескивание не всегда происходит в осевой части жилы. Ю.К. Васильчуком была разработана процедура отбора образцов из крупных ледяных жил, в которой сочетался и вертикальный и горизонтальный отбор. Вертикальный отбор производится вдоль осевой части ледяной жилы, а горизонтальный на уровнях выраженных поясковых структур.
Палеотемпературная интерпретация изотопного состава кислорода повторно-жильных льдов обсуждалась с наиболее компетентным ученым в области изотопных исследований - с проф. Дансгором. В 1995 г. проф. Васильчук специально поехал на конференцию в г. Рос Килле (Дания), посвященную изотопным исследованиям, на которую почетным гостем был приглашен профессор Вилли Дансгор (рис. 4). Именно тогда у Ю.К. Васильчука появилась возможность представить палеотемпературную интерпретацию изотопного состава повторно-жильных льдов наиболее компетентному исследователю в этой области.
Рис. 4. Конференция в Роскилле (дат. Roskilde), 1995 (Дания ), слева направо М. Дансгор, В. Дансгор, Ю. Васильчук, А. Васильчук
Рис. 5. В. Дансгор и Ю. Васильчук в низкотемпературной камере института Нильса Бора с керном GRIP
Проф. Дансгор пригласил Ю.К. Васильчука для дальнейшей дискуссии и демонстрации керна GRIP (рис. 5) и уникального масс-спектрометра (рис. 6) в институт им. Нильса Бора.
Рис. 6. В. Дансгор и Ю. Васильчук рядом с масс-спектрометром в лаборатории Отдела геофизики Института им. Нильса Бора в Университете Копенгагена. Этот прибор позволяет измерять более 180 образцов в сутки, на нем выполнены более сотни тысяч определений стабильных изотопов из ледников Гренландии (например, только из керна Саммит - более 50 тыс. образцов, а из керна Дай-3 - 60 тыс. образцов). В центре снимка система пробоподготовки, позволяющая также измерять одновременно 8 эталонов, что делает измерения исключительно точными.
В ходе дискуссии проф. Дансгор задал вопрос, почему бы не использовать его уравнение-14! для интерпретации изотопного состава ПЖЛ. Но потом, он полностью согласился с тем, что ледяные жилы являются результатом осадков зимнего сезона, а его уравнения посвящены сопоставлению среднегодовых изотопных вариаций в осадках и среднегодовых климатических характеристик.
При сопоставлении изотопных характеристик с зимними температурами возник принципиальный вопрос - уравнениями какого вида описывается имеющееся соотношение. Уже с самого начала Ю .К. Васильчуку, имеющему неплохое базовое математическое образование (матшкола при мехмате МГУ) стало ясно, что традиционные уравнения регрессии мало применимы для решения данной задачи по двум причинам. Во-первых, температура воздуха, даже осредненная за год или за несколько лет, является весьма неустойчивой характеристикой, в течение 20-30 лет среднезимняя и среднеянварская температура могут варьировать в пределах 4-6 и более градусов. Во-вторых, уравнения регрессии работают только в пределах выбранных числовых множеств и как инструмент для палеореконструкций позднего плейстоцена малоприменимы, поскольку речь идет о широком диапазоне более низкой температуры и, соответственно, о других коэффициентах уравнения регрессии и существенно меняют свою конфигурацию при добавлении каждой новой точки измерения. Одно из основных свойств уравнений регрессии заключается в том, что они достоверно описывают только
использованные для их составления пары чисел, в нашем случае значения б180 жильных ростков и среднеянварской (£ л) и среднезимней (£ температуры. Для
значений, выходящих за пределы массива данных как в положительную, так и отрицательную стороны, эти уравнения перестают работать и дают искаженные
результаты. При реконструкциях по значениям б180 из позднеплейстоценовых ПЖЛ, которые отличаются более легким составом (на 4-8%) по сравнению с современными, это особенно актуально. Поэтому Ю.К. Васильчуком, были предложены на первый взгляд очень простые соотношения без свободного члена, но с допустимой погрешностью, величина погрешности рассчитана на основании имеющихся в климатических справочниках-15! температурных рядов:
для средней температуры января t л = 1.5 б180 (±3); для средней зимней температуры t щ = б180(±2).
В этих соотношениях важны все члены уравнений: обязательно использование цифр в скобках, как меры точности реконструкций. Действительно, даже если в течение двух-трех десятков лет вариации современной среднеянварской и среднезимней температуры достигают 4-6°С, то для голоцена и, тем более для позднего плейстоцена, при палеотемпературных реконструкциях необходимо указывать доверительный интервал (т.е. пределы в которых находится точное значение определяемого показателя). Зависимость с учетом доверительного интервала на сегодняшний день является более
корректной, чем любые из пока предложенных уравнений регрессии (так как уравнения регрессии выводятся для конкретной, достаточно узкой области соотношений температуры и изотопного состава, в которую температура и изотопный состав для позднего плейстоцена естественно не включаются). Как один из путей совершенствования уравнений регрессии для целей палеореконструкций в дальнейшем можно предложить в качестве исходных условий комбинирование наблюдаемых
сочетаний значений б180 жильных ростков с зимней температурой и моделируемых предполагаемых значений этих параметров для позднего плейстоцена.
Тестирование уравнений регрессии связь изотопного состава кислорода жил с
температурами воздуха
В обзоре Т. Опеля с соавторами-30! формула Ю .К.Васильчука (б18О = ^±2) ошибочно трактуется, как зависимость, характеризующая связь изотопного состава кислорода жил с температурами воздуха периода с декабря по февраль. Однако, данная формула получена для климатического зимнего (а не календарного) периода - примерно с октября по май (т.е. того периода, когда среднесуточные температуры воздуха ниже 0°С). Там же неверно указано количество образцов современных ростков жил, на результатах изучения которых получены формулы. В действительности для жильных ростков было изучено 69 ключевых участков, на каждом участке отбиралось жильные ростки в разных фациальных обстановках т.е. 3-8 образцов на участок. Попытки уточнить
эти уравнения-9'11'13! часто подходят для современного температурного интервала, что касается позднеплейстоценовых палеореконструкций, то реконструкция до десятых долей градуса нам представляется излишней.
Для тестирования степени надежности реконструкций мы использовали уравнения, которые скопированы из работ М.А. Коняхина [9!, Х. Майера -Ш!, Г.Е. Облогова [13!:
Уравнение М.А.Коняхина (копия из автореферата, стр.15[9!):
Ь180 = 1,29 ТЗШНф- 4,23 (Г = 0,864).
Уравнения Х.Майера (копия из диссертации, стр. 93-11!):
mean winter temperature (TJ 51sO = 1.02 Tw - 0.72 (R2=0.81) (1)
mean January temperature (TJan) 51B0 = 0.59 TJan - 3.46 (R2=0.72) (2)
The mean winter temperature is related to the period, when the precipitation fails as snow. Since most of the samples were published by Vasil'chuk (1992), it is evident that the new equations are similar to his calculations.
Уравнения Г.Е.Облогова (копия из автореферата, стр. 16[13!):
Тср.янв = 1,12 5"0 - 6,43, Я2 = 0,745, а = 2,6, Тер.зим = 1,15 5180 - 4,6, Я2 = 0,754, а = 2,7, Тср.хол - 0,885 - 2,55, Я2 - 0,674, а = 2,7. Разница оценок по сравнению с данными Ю.К. Ваеильчука (1992, 2006) составляет менее 2 °С. С вероятностью 0,85 пределы варьирования частных
Таблица 1. Тестирование изотопно-температурных уравнений разных авторов при изменении значений б18О от -15 до -27%о и реконструированная палеотемпература (°С) позднего плейстоцена при изменении значений б180 от -30 до -35%
автор формулы, год формулы для реконструкции среднеянварской и средне зимней температуры б180 = -15 °/оо (Амдер-ма, совр.) б180 = -18%о (Ямал, совр. Р.Щу-чья) б180 = -20%о (п-ов Таймыр, совр. Мыс Саблера) б180 = -25%о (долина р.Колымы -совр. Анюйск) б180 = -27%о, ( Мамонтова Гора, совр.) б180 = -30%о (долин р.Колы мы палео)
Ю.К. Васильчук[16-18]* t ] = 1.5 б180 (±3) -22.5 (±3) 27(±3) -30(±3) -37.5 (±3) -43.5 (±3) -45(±3)
t „ = б180(±2) -15(±2) 18(±2) -20(±2) -25(±2) -29(±2) -30(±2)
М.А. Коняхин^9-** б180 = 1.29£зимн -4.23 -8.35 -10.67 -12.22 -16.10 -19.2 -19.97
Н.МеуегШ!*** б180= 0,59£ г3,46 -19.55 -24.64 -28.03 -36.5 -43.29 -44.98
б180 = 1,02£ м -0,72 -14 -16.94 -18.90 -23.8 -27.73 -28.70
Г.Е.Облогов[13!**** 1ср.янв= 1,12 б180 -6,43 -23.23 -26.59 -28.83 -34.43 -38.91 -40.03
1ср.хол= 0,885 б180 -2,55 -15.83 -18.48 -20.25 -24.68 -28.22 -29.1
Климатический справочник^15! Среднеянварская темпера тура, °С -19 -24 -33 -35 -45
Среднезимняя темпера тура (период когда среднесут. темп. ниже 0°) °С -13 -16 -21 -23 -26
^Трактовка понятия «среднезимняя температура» Ю.К.Васильчуком в соответствии с климатическим справочником - климатический период отрицательных среднесуточных температур.
**Трактовка понятия «среднезимняя температура» М.А.Коняхиным в соответствии с
климатическим справочником.
***Трактовка среднезимней температуры у Х. Майера - период, когда осадки выпадают в твердом виде.
****Трактовка г.Е. Облоговасредней температуры холодного периода совпадает с определением среднезимней температуры, которое используется в климатических справочниках, под среднезимней температурой Г.Е.Облогов понимает среднюю температуру трех месяцев: декабря, января и февраля.
Результаты реконструкций палеотемператур по формулам разных авторов демонстрируют особенности реконструкций с использованием уравнений регрессии, а именно их использование дает надежные результаты, если для их построения применены те натурные данные, которые подчиняются тем же закономерностям, что и реконструируемые параметры. В двух нижних строках табл. 1 приведены реальные температуры для оценки степени надежности реконструкции.
Как видно из тестирования (табл. 1) значения среднезимних и среднеянварских температур при реконструкции по формулам Ю.К. Васильчука всегда попадают в доверительный интервал во всем современном температурном диапазоне. При
реконструкции палеотемпературы со значениями б18О от -30 до -35% реконструируемые значения примерно совпадают при использовании формул Ю.К. Васильчука и Х. Майера, значения среднеянварских температур, полученных по уравнениям Г.Е.Облогова сильно отличаются более, чем на 5-7%. Как отмечал сам Х. Майер поскольку уравнения основаны на данных Ю.К.Васильчука результаты реконструкций часто идентичны. Однако даже уравнения Х. Майера все же
обнаруживают существенные отклонения при значении б18О равном -20%: абсолютная величина отклонения от реальных значений среднеянварской температуры 4.70С и 2.1 0С от реальных значений среднезимней температуры. Возможно, это связано с тем, что, работая в сравнительно более узком регионе на севере Таймыра и Арктических островах, он исключил точки по более восточным и западным регионам из рассматриваемого множества. В формуле Г.Е.Облогова, связывающей среднезимнюю температуру и изотопный состав ростков присутствует ошибка, заключающаяся в том, что коэффициент, связывающий в его же уравнении среднеянварскую температуру
(коэффициент 1,12) со значением б18О меньше, чем коэффициент для среднезимней температуры (коэффициент 1,15) при близких значениях свободного члена. При низких
(позднеплейстоценовых) значениях б18О (-35 - -37 %) среднезимняя температура (напомним, что под среднезимней температурой Г.Е. Облогов понимает среднюю температуру трех месяцев: декабря, января и февраля) будет отрицательнее, чем значение среднеянварской температуры. При этом чем ниже значения б18О, тем более неправдоподобны результаты. Так при значениях -37 % значение среднеянварской температуры -47.9 ос, а средняя температура периода с декабря по февраль по уравнениям Г.Е. Облогова -47.2 ос. Это совершенно нереалистично. Так, например вблизи современного полюса холода в современном Батагае значение среднеянварской температуры равно -47.9 ос, а средняя температура периода с декабря по февраль равна -40.6 ос. По результатам тестирования уравнение Г.Е.Облогова для реконструкции среднеянварской температуры дает результаты, совпадающие с натурными
наблюдениями при значении б18О -25 %, при других значениях б18О разница с данными справочников более заметная, ее абсолютные значения составляют 4.230С при
величине б18О равной -15%, 2.590С при значении б18О равном -18%, 6.090С при
значении б18О равном -27%.
Поэтому значения позднеплейстоценовых среднеянварских температур, по уравнению Г.Е. Облогова
получаются для центральных районов Якутии выше, чем даже современные, например для позднеплейстоценовых жил Батагайского едомного комплекса характерны значения
б18О -35 %, реконструированная средняя температура января равна -45.630С, а период с 1978 по 2014 г она 6 раз была ниже этих значений. Заметим, что для этого
района значения б18О в ростках ледяных жил составляет -29 %[19,31!. Такая же существенная ошибка отмечается и для расчетов по Куларской едоме. Таким образом реконструкции позднеплейстоценовых среднеянварских температур, рассчитываемые по уравнению Г.Е. Облогова неадекватны реальности.
Рассчитываемые по уравнению Г.Е. Облоговатестируемые температуры холодного периода (в понимании Ю.К.Васильчука и Х.Майера - это среднезимняя температура)
достоверны в интервале значений б18О от -20 до -27%, полученные значения близки к реальным, но для значений -15% и -18% реконструируемая температура не совпадает с данными справочника в абсолютном выражении на 2.83 и 2.48 соответственно (табл.
1). В диапазоне низких значений б18О (от -30 до -35%), согласно уравнению Г.Е. Облогова получены значения палеотемпературы холодного периода близкие к результатам по формулам Ю .К.Васильчука и Х.Майера. Согласно расчетам Г.Е. Облогова
разница оценок по сравнению с данными Ю.К. ВасильчукаГ16-18! составляет менее 2 °С, а пределы варьирования частных значений определяемой по уравнениям регрессии температуры составляют ±3,8 °С с вероятностью 0,85Г13!. Уравнение М.А. Коняхина содержит ошибку, тестирование его уравнения дает неудовлетворительный результат во всем диапазоне, абсолютная величина отклонения от реальных значений составляет 4.65-8.78°С. На наш взгляд построение уравнений регрессии на данном уровне знаний скорее является шагом назад, поскольку имеет существенные и непредсказуемые ограничения при интерпретации, которые невозможно определить, интерпретируя данные по плейстоценовым повторно-жильным льдам.
Результаты обратного тестирования, приведенные в табл. 2, для формул реконструкции среднезимней температуры, демонстрируют, что для уравнений регрессии имеются интервалы, в которых реконструкция вполне совпадает с натурными данными, но есть интервалы, где расхождение с натурными данными очень существенно.
Таблица 2. Значения б180 в современных ростках ледяных жил, полученные в результате тестирования уравнений разных авторов и их реальные измеренные значения (%).
а вто р формул ы, год формулы для тестирования среднезимней температуры Современные среднезимние температуры (период когда среднесут. темп. ниже 0°) по метеостанциям, °С Г15, 33!
-13 Амдерма -16 Устье р.Сеяха -21 Мыс Саблера -22 П.Черский -24 П.Ко-лымское -27 Мамонтова Гора
Значения б180 в современных ростках ледяных жил рассчитанные по тестированию формул среднезимних те мпе ра тур, %
Ю.К. t „ = б180(±2) -13(±2) - -21(±2) -22(±2) -24(±2) -27(±2)
Васильчук^16-18! 16(±2)
М.А. Коня-хин!9- б180=1.29£ зимн -4.23 -21 -25.07 -31.32 -32.61 -35.19 -39.06
Н.МеуегШ! б180=1.02£ м -0.72 -13.98 -17.04 -22.14 -23.16 -25.2 -28.26
Г.Е. Облогов^19 ^с р.хол = °.885 б180 -2.55 -11.81 -15.20 -20.85 -21.98 -24.24 -27.63
Значения б180 измеренные в современных ростках ледяных жил,°/оо -15.2 -18 -20.4 -24 -26 -26.3
Для уравнения Х. Майера результаты, отличающиеся от реальных значений б18О в жильных ростках отмечены при значениях среднезимней температуры -21 и -27°С (1,74 и 1,96 % соответственно), а для уравнения Г.Е. Облогова - неудовлетворительные результаты тестирования характерны для интервала относительно высокой среднезимней
температуры от -13 до -18°С (разница с реальными значениями б18О составляет 3.42.8%) Отметим также, что соотношение Ю .К.Васильчука работает во всем диапазоне. Уравнение М.А. Коняхина демонстрирует неправдоподобные результаты противоречащие реальным данным для всего диапазона температур характерных для Арктических районов (из-за содержащейся в нем алгебраической ошибки).
Соотношения, предложенные Ю.К. Васильчуком, были дополнительно протестированы на современных данных. Для повышения точности палеотемпературных реконструкций на современном этапе потребовалось внести некоторые корректировки в интерпретацию изотопных данных по повторно-жильным льдам, что выполнено проф. Ю .К.Васильчуком с
климатологом проф.Г.В. Сурковой^331. Поскольку для интерпретации результатов изотопных исследований в основном были использованы данные по росткам жилок, сформировавшихся до 1990-х годов и соответствующие им температуры из справочников 1960-х годов, то потребовалось обновление фактологической базы и уточнение соотношения изотопного состава современных ростков повторно-жильных льдов и температуры холодного периода с учетом климатических изменений последних 50 лет. Для верификации корреляционной зависимости современных зимних температур воздуха (среднеянварской и среднезимней) были использованы данные метеорологических
наблюдений с 1930 г. по 2017 г. и значения б18О в ростках ледяных жил по каждому региону криолитозоны с целью выявления наиболее явной и устойчивой зависимости изотопного состава и зимней температуры (табл. 3). Рассчитанные коэффициенты корреляции (достоверность 95%) между данными наблюдений на станциях Росгидромета
РФ и содержанием 18О в ростках повторно-жильных льдов показывают очень высокий уровень связи между этими показателями и свидетельствуют о высоком вкладе температурного режима холодного сезона в процесс формирования ПЖЛ.
Таблица 3. Распределение современной среднеянварской и среднезимней температуры, вблизи точек определения изотопного состава ростков ПЖЛ для трех периодов: а) с
1930 по 1966 гг., б) с 1967 по 2000 гг., в) с 2001 по 2019 гг. (по 1341 Васильчук, 1992 для периода 1930-1966 гг.). Фрагмент таблицы из статьи Ю.К. Васильчука, Г.В. Сурковой
[33!(2020).
Местоположение современных жилок(Recent veinlets) 518О , veinlets' %о ^р.з. tя ^р.з. tя ^р.з. tя
Период метеонаблюдений 1930-1966 1967-2000 2001-2017
г. Амдерма -15.2 -13 -19 -12 -20 -11 -17
р. Щучья -18.2 -16 -24 -15 -25 -14 -24
Мыс Саблера -20.4 -21 -33 -21 -33 -21 -30
Мамонтова Гора -26.3 -27 -45 -27 -44 -26 -42
Для реконструкции температурных условий были привлечены все имеющиеся данные наблюдений на метеостанциях, ближайших к точкам определения изотопного состава ростков, от мыса Шприндлера и Воркуты на западе криолитозоны до Уэлена и Анадыря на востоке Чукотки и от арктических островов на севере, до Чарской и Тоджинской
котловин на юге криолитозоны с созданием электронной базы данных значений б18О в ростках современных повторно-жильных льдов. По результатам сопоставления
изотопного состава ростков повторно-жильных льдов (518О^П|6^) и показателей зимней
температуры (£ ср.з., t я) получены значения соотношений для каждого периода-33-, которые очевидно следует учитывать при палеореконструкциях тех периодов, когда зимы были существенно теплее современных.
Говоря о возможном искажении первичного изотопного сигнала снега перед заполнением морозобойной трещины в результате изотопного фракционирования в снежном покрове в течение зимы, можно отметить, что подобные эффекты ожидаемы в районах с резко-континентальным климатом (например, Центральная Якутия). Но даже там далеко не во всех случаях исследования изотопного состава как жильных ростков, так и жил более древнего возраста, была зафиксирована изотопная трансформация, обусловленная испарением или сублимацией снежного покрова. Так, например, в северном Забайкалье, в долине р. Чара, в районе с резко континентальным климатом, где снег может практически полностью испариться в течение зимы, во льду
современного жильного ростка было получено среднее значение б18О -21.5 % Г34,35! что совпадает со среднезимней температурой воздуха в этом районе (-21 °С), а со среднеянварской температурой (составляющей в этом районе -33°С) соотносится с коэффициентом 1.6. Выполненное исследование современных жилок, сформировавшихся в условиях холодных и сухих зим Центральной Якутии и Забайкалья, и жилок, образовавшихся в более мягких зимних условиях на Ямале и Чукотке показывает, что изотопный состав кислорода жилок напрямую соответствует зимним температурам. Таким образом, можно уверенно говорить о том, что изотопный "термометр", созданный для современных ледяных жилок, может адекватно работать и при палеотемпературных реконструкциях прошлых эпох-31!. Об изменении начальных изотопных характеристик снега в результате испарения можно судить по соотношению б180-б2Н. В случае
активных процессов испарения снега, значения б180 и б2Н заметно отклоняются от глобальной линии метеорных вод. Обобщение большей части имеющихся в настоящее время данных об изотопном составе кислорода и водорода ростков ледяных жил полигональных массивов на территории севера Евразии от устья реки Нгарка-Тамбъяха на северо-востоке Европы до пос. Лорино на восточном побережье Чукотки, показало, что для большинства образцов наклон линий соотношения б180-б2Н близок к наклону
для глобальной линии метеорных вод. Также было установлено, что в меридиональном распределении значений как 618O так и б2Н в жильных ростках выявлен тренд снижения значений от побережья Байдарацкой губы на севере Европейской части России до побережья моря Лаптевых (от -141,9% до -193,5%) и далее на восток - повышение значений (до -99^-122%) на восточном побережье Чукотки, что отражает изотопное облегчение осадков при перемещении атлантических воздушных масс над большей частью Российской криолитозоны и преобладающее влияние воздушных масс Тихого
океана на Чукотке-^36!. С 2009 г все данные по содержанию стабильных изотопов кислорода и водорода в повторно-жильных льдах и осадках получены в лаборатории стабильных изотопов кафедры геохимии ландшафтов и географии почв Географического факультета МГУ под руководством проф. Васильчука. Измерения проводятся на масс-
спектрометре Delta V Plus. Измеренные величины 6*8O и б2Н калибруются относительно международных стандартов. В редких случаях отмечаются пробы со значениями, которые могут указывать на изотопную трансформацию снега перед заполнением морозобойных трещин. Кроме того, анализ изотопного состава снежных осадков и снежного покрова показывает, что вариации значений dexc варьируют в широком диапазоне (как правило,
от 3 до 20%, дополнительную информацию можно получить, проанализировав изотопные данные по зимним осадкам на станциях GNIP в российской Арктике), что может объясняться особенностями испарения влаги в очаге формирования воздушных масс и особенностями их передвижения над континентом. При этом значения dexc в
жильном льду, как в современных жилках, так и в более древних, варьируют, как правило, в диапазоне от 4 до 15%, и в основном находятся внутри диапазона для снега, что еще раз подтверждает отсутствие заметного влияния изотопного фракционирования на изотопный состав снега.
Модель гетероциклического развития сингенетических повторно-жильных льдов
Важным моментом в развитии изотопной геокриологии явилась разработанная Ю.К. Васильчуком модель гетероциклического развития сингенетических повторно-жильных
льдов^37!. На основе анализа собственных полевых наблюдений и опубликованных материалов в виде описаний, фотографий и рисунков многих исследователей им была предложена новая структурная модель формирования мощных сингенетических повторно-жильных льдов. П роцесс сингенетического формирования мощных ледяных жил рассматривается не в рамках принятой парадигмы - как процесс непрерывного формирования льда, а как пульсирующий - циклический процесс. Главное отличие этой модели от имеющихся моделей циклического развития жил состоит в не климатическом запускающем механизме выявленной цикличности (т.е. цикличность проявляется независимо от климатических ритмов потеплений или похолоданий). Основным детерминирующим механизмом является неоднократная повторяющаяся смена характера осадконакопления на поверхности полигонального массива - субэрального и субаквального. Модель гетероциклического развития сингенетических повторно-жильных льдов обосновала вертикальный отбор ПЖЛ вдоль оси и горизонтальный отбор на уровнях, где заметны признаки субаэрального роста жил: прослои торфа, «плечики» и др., а также позволила более адекватно оценить скорости роста ПЖЛ, надежность радиоуглеродных датировок, поскольку вероятность накопления переотложенной органики в субаквальных отложениях намного выше.
Прямые радиоуглеродные датировки из повторно-жильных льдов
Ю.К. Васильчук получил прямые радиоуглеродные датировки из повторно-жильных
льдов с применением ускорительной масс-спектрометрии впервые в мире (в 1998 г.), что для развития изотопной геокриологии имело большое значение. Были датированы микроорганические включения 200-400 мкм, включения растительных макроостатков, и общий органический углерод во льдах Сеяхинской и Щучьинской толщ на Ямале, разрезов Бизон, Плахинский Яр, Зелёный Мыс, Дуванный Яр, Батагай, Мамонтова Гора в
Якутии, Анадырь на Чукотке I"38-421. Это позволило не только определить возраст жил, но и доказать их сингенетичность вмещающим отложениям, а также продемонстрировало омоложение возраста жил снизу-вверх. Прямое датирование льда позволило привязать имеющиеся данные по изотопному составу повторно-жильных льдов к радиоуглеродной шкале и сопоставить изотопный состав повторно-жильных льдов и гренландских ледниковых кернов. Прямое датирование повторно-жильного льда снижает остроту дискуссии по поводу того как производить отбор ПЖЛ вдоль оси или по горизонтали. Иллюстрацией могут служить результаты датирования ледяных жил по горизонтали в центральном Юконе-431. Поскольку растрескивание происходит хаотично, как было
показано натурными наблюдениями Дж.Р.Маккая [44-481 14С даты с обеих сторон от центральной части жилы оказались асимметричными (31 500 лет cal BP слева, 13 685 лет cal BP справа). Кроме того, датировка в центре ледяной жилы была старше (6360 лет cal BP), чем на правой стороне (2450 cal BP). Три из четырех ледяных жил, которые были отобраны горизонтально, продемонстрировали ассиметричные и значительно
различающиеся значения 618O, б2Н и d-эксцесса в левой и правой сторонах ледяного клина. Эти исследования продемонстрировали, что растрескивание и заполнение талой водой не обязательно происходит вблизи центральной части ледяной жилы.
Типы цикличность формирования сингенетических повторно-жильных льдов
За прошедшие 30 лет после того как Ю.К. Васильчук создал новое научное направление, он развил его, получив важные результаты. Во-первых, он разработал модель циклически-пульсирующего формирования сингенетических толщ с мощными
повторно-жильными льдами-371, во-вторых разработал новую стратегию радиоуглеродного датирования позднеплейстоценовых сингенетических повторно-жильных льдов и точной привязки во времени диаграмм распределения стабильных изотопов включающая использование ускорительной масс-спектрометрии (AMS) для прямого датирования возраста ледяных жил по микровключениям разнообразного
органического материала —501—511, датировал по углероду позднеплейстоценовые повторно-жильные льды впервые в мире-41,421. Им установлена вертикальная и латеральная гетерохронность и гетерогенность едомных массивов, что имеет важное значение для оценки их возраста. Для оценки процессов формирования синкриогенных толщ с повторно-жильными льдами им выделено три типа цикличности: микро-, мезо и макроцикличность. Микроцикличность, связана с сезонной периодичностью изменения глубины деятельного слоя. Мезоцикличность связана с пульсирующим изменением уровня водоема, по берегам которого, или на отмелях которого идет формирование жил. Макроцикличность связана с коренной перестройкой режима седиментации или реже (в основном на юге ареала повторно-жильных льдов) с крупными климатическими
осцилляциями [49,50]. Разработал стратегию датирования сингенетических многолетнемерзлых отложений с повторно-жильными льдами, предложив принцип
выбора наиболее молодой датировки ISÜ. Для изотопной записи в сингенетических повторно-жильных льдах криолитозоны Северной Евразии были установлены единичные спектральные пики, длительностью от 1.5 тыс. лет, идентифицируемые как события Дансгора-Эшгера.За прошедшие 30 лет появилось много работ посвященных изотопному
составу повторно-жильных льдов. Изучены арктические острова, дельта р. Лены,
разрезы Якутии [11,18, 19,20~30,52,53]. Отметим, что практически все полученные данные по изотопному составу подземных льдов, имеющие возрастную привязку, вписываются в распределения зимних температур на территории криолитозоны Евразии, которые были
опубликованы в 1992 г. (см. рис.1). Например, Портер и Опель^53 обобщили данные по изотопному составу повторно-жильных льдов из 82 точек в Арктике. Во-первых, было обнаружено, что прибрежные повторно-жильные льды изотопически тяжелее, по сравнению с более континентальными; это объясняется их удаленностью от источника
влаги, что было установлено в работе Ю.К. Васильчука в 1991 г.^19^. Во-вторых Опель и
Портер также пришли к выводу, полученному ранее Ю.К. ВасильчукоД34!, о том, что преобладание западного переноса влаги из Атлантики и закономерности ведет к
последовательному облегчению изотопного состава (значений ö18O и б2Н) с запада на восток в Северной Евразии.
Развитие Ю.К. Васильчуком направления изотопной геокриологии еще раз д е м о н с тр и р у е т приоритетность российских геокриологических исследований сингенетических отложений с повторно-жильными льдами - этого важнейшего и одного из самых сложных в изучении типов подземных льдов. Разработанные новые позиции соз да ют более объективный подход как к анализу получаемых данных, так и к самому процессу полевого изучения сложно построенных массивов с полигонально-жильными льдами и родственных им геологических образований.
Библиография
1. Васильчук Ю.К., Трофимов В.Т. Изотопно-кислородная диаграмма повторно-жильных льдов Западной Сибири, ее радиологический возраст и палеогеокриологическая интерпретация // Доклады АН СССР. 1984. Том 275. N 2, с.425-428.
2. Gallwitz H. Eiskeile and glaziale sedimentation // Geoiogica. 1949. Bd 2. 24 s.
3. Попов А. И. Морозобойные трещины и проблема ископаемых льдов // Труды Ин-та мерзлотоведения им. В.А.Обручева АН СССР, том IX. Вечная мерзлота различных районов СССР. М.: Изд-во АН СССР. 1952, с. 5-18.
4. Попов А.И. Особенности литогенеза аллювиальных равнин в условиях сурового климата // Известия АН СССР. Серия географическая. 1953. №2. С. 29-41.
5. Достовалов Б.Н. О физических условиях образования морозобойных трещин и развития трещинных льдов в рыхлых породах // В сб.: Исследование вечной мерзлоты в Якутской республике, вып. 3. М.: Изд-во АН СССР, 1952. С.162-194.
6. Романовский Н.Н. 1977 Формирование полигонально-жильных структур Новосибирск: Наука, 1977, 212 с.
7. Васильчук Ю.К. Изотопно-кислородный состав подземных льдов // Материалы гляциологических исследований, вып.66. 1989, с. 196-210.
8. Васильчук Ю.К. Условия формирования позднеплейстоценовых и голоценовых повторно-жильных льдов Чукотки (изотопно-криолитохронологический анализ)// Доклады АН СССР. 1989. Том 309. №4, с. 920-924.
9. Коняхин М.А Изотопно-кислородный состав полигонально-жильных льдов как показатель условий их формирования и генезиса / Автореферат дисс. на соискание ученой степени кандидата географических наук. МГУ. Москва. 1988. 24 с.
10. Nikolaev V.I., Mikhalev D.V. An oxygen isotope paleothermometer from ice in Siberian permafrost // Quat. Res. 1995 Vol. 43. Iss. 1 p. 14-21. doi: 10.1006/qres.1995.1002.
11 Meyer H. Late Quaternary climate history of Northern Siberia - evidence from ground ice PhD thesis. 2003 Ber. Alfred-Wegener-Institut Polar- und Meeresforschung Forschungsstelle Potsdam Polarforsch. Meeresforsch. 461, 110 p.
12. Стрелецкая И.Д., Васильев А.А., Облогов Г.Е., Токарев И.В. Реконструкция палеоклимата Российской Арктики в позднем неоплейстоцене-голоцене на основе данных по изотопному составу полигонально-жильных льдов // Криосфера Земли. 2015. Том XIX, № 2, с. 98-106.
13. Облогов Г.Е. Эволюция криолитозоны побережья и шельфа Карского моря в позднем неоплейстоцене - голоцене / Автореферат дисс. на соискание ученой степени кандидата геолого-минералогических наук. ИКЗ РАН. Тюмень. 2016. 24 с.
14. Dansgaard W. Stable isotopes in precipitation // Tellus. 1964. Vol. 16. N4. P.436 -468.
15. Справочник по климату СССР. Л. Гидрометеоиздат. Во всех выпусках - Часть 2. Температура воздуха и почвы. Вып. 1. Архангельская и Вологодская области, Карельская и Коми АССР. 1965. 359 с. Вып. 17. Тюменская и Омская области. 1965. 400 с. Вып. 20. Томская, Новосибирская, Кемеровская области и Алтайский край. 1966. 396 с. Вып. 21. Красноярский край и Тувинская АССР. 1966. 504 с. Вып. 23. Бурятская АССР и Читинская область. 1966. 319 с. Вып. 24. Якутская АССР. 1966. 403 с. Вып. 33. Чукотский национальный округ и Магаданская область. 1967. 288 с.
16. Васильчук Ю.К. Корреляция изотопно-кислородного состава повторно-жильных льдов со среднезимними и среднеянварскими температурами воздуха // Изотопы в гидросфере: Тезисы докладов 3-го Всесоюзного симпозиума. г.Каунас. 29 мая - 1 июня 1989 г. М.: Изд-во ИВП АН СССР. 1989. С. 82-83.
17. Васильчук Ю.К. Реконструкции палеоклимата позднего плейстоцена и голоцена на основе изотопных исследований подземных льдов и вод криолитозоны // Водные ресурсы. 1990. №6. С. 162 - 170.
18. Vasil'chuk Yu. K. Reconstruction of the palaeoclimate of the Late Pleistocene and Holocene of the basis of isotope studies of subsurface ice and waters of the permafrost zone // Water Resources. 1991. Vol. 17. №. 60, p. 640-647.
19. Васильчук Ю.К. Позднечетвертичные синкриогенные толщи севера Евразии: Строение, изотопно-кислородный состав и условия формирования / Диссертация на соиск. уч. степени доктора геол.-минерал. наук. М. 1991. В 2-х томах. Том 1 - 622 с., том 2. - 244 с.
20. Schirrmeister, L., Grosse, G., Schwamborn, G., Andreev, A. A., Meyer, H., Kunitsky, V. V., Kuznetsova, T. V., Dorozhkina, M. V., Pavlova, E. Y., Bobrov, A. A., Oezen, D. Late Quaternary history of the accumulation plain north of the Chekanovsky Ridge (Lena Delta, Russia): a multidisciplinary approach // Polar Geography. 2003. Vol. 27, p.2 77319. doi: 10.1080/789610225.
21. Schirrmeister, L., Grosse, G., Kunitsky, V., Magens, D., Meyer, H., Dereviagin, A., Kuznetsova, T., Andreev, A., Babiy, O., Kienast, F., Grigoriev, M., Overduin, P. P., Preussner, F. Periglacial land scape evolution and environmental changes of Arctic lowland areas for the last 60 000 years (Western Laptev Sea coast, Cape Mamontov Klyk) // Polar Research. 2008. Vol. 27, p. 249-272. doi: 10.1111/j.1751-8369.2008.00067.
22. Popp S., Diekmann B., Meyer H., Siegert C., Syromyatnikov I., Hubberten H. W. Palaeoclimate signals as inferred from stable-isotope composition of ground ice in the Verkhoyansk foreland, Central Yakutia // Permafrost Periglacial Processes. 2006. Vol. 17, p. 119-132. doi. 10.1002/ppp.556.
23. Murton Ju.B., Goslar T., Edwards M.E., Bateman M.D., Danilov P.P., Savvinov G.N., Gubin S.V., Ghaleb B., Haile J., Kanevskiy M., Lozhkin A.V., Lupachev A.V., Murton D.K., Shur Yu., Tikhonov A., Vasil'chuk A.C., Vasil'chuk Yu.K., Wolfe S.A. Palaeoenvironmental Interpretation of Yedoma Silt (Ice Complex) Deposition as Cold-Climate Loess, Duvanny Yar, Northeast Siberia // Permafrost and Periglacial Processes. 2015. Vol. 26. Iss. 3, p. 208-288. doi: 10.1002/ppp.1843.
24. Wetterich S., Tumskoy V., Rudaya N., et al. Ice Complex formation in arctic East Siberia during the MIS3 Interstadial // Quaternary Science Reviews, 2014. Vol. 84, p. 39-55. doi: 10.1016/j.quascirev.2013.11.009.
25. Wetterich, S., Meyer, H., Fritz, M., Mollenhauer, G., Rethemeyer, J., Kizyakov, A., et al. Northeast Siberian permafrost ice-wedge stable isotopes depict pronounced last Glacial maximum winter cooling // Geophysical Research Letters. 2021. Vol. 48, e2020GL092087. doi: 10.1029/2020GL092087.
26. Wetterich, S., Kizyakov, A., Fritz, M., Wolter, J., Mollenhauer, G., Meyer, H., et al. (2020). The Cryostratigraphy of the Yedoma Cliff of Sobo-Sise Island (Lena delta) Reveals Permafrost Dynamics in the central Laptev Sea Coastal Region during the Last 52 Kyr // The Cryosphere. 2020. Vol. 14, p. 4525-4551. doi:10.5194/tc-14-4525-2020.
27. Wetterich S., Rudaya N., Nazarova L., Syrykh L., Pavlova M., Palagushkina O., Kizyakov A., Wolter J., Kuznetsova T., Aksenov A., Stoof-Leichsenring K.R., Schirrmeister L., Fritz M. Paleo-Ecology of the Yedoma Ice Complex on Sobo-Sise Island (EasternLena Delta, Siberian Arctic) // Front. Earth Sci. 2021. Vol. 9: 681511. doi: 10.3389/feart.2021.681511.
28. Opel T., Wetterich S., Meyer H., Dereviagin A. Y., Fuchs M. C., Schirrmeister L. Ground-ice stable isotopes and cryostratigraphy reflect late Quaternary palaeoclimate in the Northeast Siberian Arctic (Oyogos Yar coast, Dmitry Laptev Strait) // Climate of the Past. 2017. Vol.13, p. 587-611. doi: 10.5194/cp-13-587-2017.
29. Opel T, Meyer H, Wetterich S, Laepple T, Dereviagin A, Murton J. Ice wedges as archives ofwinter paleoclimate: A review // Permafrost and Periglacial Processes. 2018. Vol.29, p. 199-209. doi: 10.1002/ppp.1980.
30. Murton J.B., Opel T., Toms P., Blinov A., Fuchs M., Wood J., Gärtner A., Merchel S., Rugel G., Savvinov G., Wetterich S. A multimethod dating study of ancient permafrost, Batagay megaslump, east Siberia // Quaternary Research. 2021. p. 1-22. doi: 10.1017/qua.2021.27.
31. Васильчук Ю.К., Котляков В.М. Основы изотопной геокриологии и гляциологии. М., Изд-во Моск. ун-та. 2000. 616 С.
32. Васильчук Ю.К. Позднечетвертичные синкриогенные толщи севера Евразии: Строение, изотопно-кислородный состав и условия формирования / Автореф. диссертации докт. геолого-минералогических наук. М. 1991. 48 с.
33. Васильчук Ю.К., Суркова Г.В. Верификация соотношения изотопного состава повторно-жильных льдов и температуры холодного периода за последние 80 лет на севере криолитозоны России // Метеорология и гидрология. 2020. №11. С. 84-91 (Vasil'chuk, Y.K., Surkova, G.V. 2020. Verification of the Relationship between the Isotopic Composition of Ice Wedges and Cold-season Temperature over the Recent 80 Years in the Northern Permafrost Zone of Russia. Russ. Meteorol. Hydrol. 45, 791-796. doi: 10.3103/S1068373920110060).
34. Васильчук Ю.К. Изотопно-кислоpодный состав повтоpно-жильных льдов (опыт палеогео^иологических pеконстpукций): В 2 т. - М.: изд. Отдел теоpетических пpоблем PАН, МГУ, ПНИИИС, 1992, Т. 1. 420 https://e.mail.ru/messages/inbox/с.; Т.
2, 264 с.
35. Vasil'chuk Yu K., Vasil'chuk A.C., Stanilovskaya Ju V. Early Holocene climate signals from stable isotope composition of ice wedge in the Chara Basin, Northern Transbaikalia, Russia // Geoscience Frontiers, 2018. Vol. 9, № 2, p. 471-483. doi: 10.1016/j.gsf.2017.04.008.
36. Буданцева Н.А., Васильчук Ю.К. Вариации изотопов кислорода в ростках современных сингенетических повторно-жильных льдов в низовьях реки Колымы // Арктика и Антарктика, 2019. № 3, с. 39-53. doi: 10.7256/2453-8922.2019.3.30744.
37. Васильчук Ю.К. Модель циклически-пульсирующего формирования сингенетических толщ с мощными повторно-жильными льдами // Криосфера Земли. 1999. Том 3. N2. С. 50-61.
38. Vasil'chuk Yu. K., van der Plicht J., Jungner H., Vasil'chuk A.C. AMS-dating of Late Pleistocene and Holocene syngenetic ice-wedges // 8th International Conference on Accelerator Mass Spectrometry. Abstracts. Palais Auersperg, Vienna, Austria 6-10 September. 1999, p. 141.
39. Vasil'chuk Yu.K., van der Plicht J., Jungner H., Sonninen E., Vasil'chuk A.C. First direct dating of Late Pleistocene ice-wedges by AMS // Earth and Planetary Svience Letters. 2000. Vol.179. N2. P. 237-242. doi: 10.1016/S0012-821X(00)00122-9.
40. Vasil'chuk Yu.K., van der Plicht J., Jungner H., Vasil'chuk A.C. AMS-dating of Late Pleistocene and Holocene syngenetic ice-wedges // Nuclear Instruments and Methods in Physics Research. Section B: Beam Interactions with Materials and Atoms. 2000. Vol.172. P. 637-641. doi: 10.1016/S0168-583X(00)00212-3.
41. Васильчук Ю.К., ван дер Плихт Й., Васильчук А.К., Юнгнер Х., Соннинен Э. Первые радиоуглеродные датировки сингенетических позднеплейстоценовых повторно-жильных льдов // Доклады Российской Академии Наук. 2000. Том 371. № 1. С. 114117.
42. Васильчук Ю.К., Ким Ч.Ч., Васильчук А.К. Радиоуглеродное AMS-датирование изотопных диаграмм позднеплейстоценовых повторно-жильных льдов // Доклады Российской Академии Наук. 2002. Том 383. №3. С. 390-396.
43. Grinter M., Lacelle D., Baranova N., Murseli S., Clark I. Late Pleistocene and Holocene ice-wedge activity on the Blackstone Plateau, central Yukon, Canada // Quaternary Research. 2019. Vol. 91(1), p. 179-193. doi:10.1017/qua.2018.65.
44. Mackay J. R. Ice-wedge cracks, Garry Island, Northwest Territories // Canadian Journal of Earth Sciences. 1974. Vol. 11, p. 1366-1383. doi: 10.1139/e74-133.
45. Mackay J. R. Some observations on the growth and deformation of epigenetic, syngenetic and anti-syngenetic ice-wedges // Permafrost and Periglacial Processes, 1990. Vol. 1, p. 15-29. doi: 10.1002/ppp.3430010104.
46. Mackay J. R. The frequency of ice-wedge cracking (1967-1987) at Garry Island, western Arctic coast, Canada // Canadian Journal of Earth Sciences. 1992. Vol. 29, p. 236-248. doi: 10.1139/e92-022.
47. Mackay J. R. Thermally induced movements in ice-wedge polygons, western arctic coast: a long-term study // Geographi e physique et Quaternaire. 2000. Vol. 54 (1), p. 41-68. doi: 10.720 2/004846ar.
48. Mackay J.R., Burn C. R. The first 20 years (1978-1979 to 1998-1999) of ice-wedge growth at the Illisarvik experimental drained lake site, western Arctic coast, Canada // Canadian Journal of Earth Sciences. 2002. Vol. 39(1), p. 95-111. https://doi.org/10.1139/e01-048
49. Васильчук Ю.К. Повторно-жильные льды: гетероцикличность, гетерохронность,
гетерогенность. М.: Изд-во Моск. ун-та. 2006. 404 с.
50. Vasil'chuk Y. K. Syngenetic ice wedges: cyclical formation, radiocarbon age and stable-isotope records // Permafrost and Periglacial Processes. 2013. Vol. 24(1), p. 82-93. doi:10.1002/ppp.1764.
51. Vasil'chuk Y.K., Vasil'chuk, A.C. Validity of radiocarbon ages of Siberian yedoma // GeoResJ. 2017. Vol. 13, p. 83-95. doi: 10.1016/j.grj.2017.02.004.
52. Коняхин М.А., Михалев Д.В., Соломатин В.И. Изотопно-кислородный состав подземных льдов. М.: Изд-во МГУ. 1996. 156 с.
53. Porter T. J., Opel T. Recent advances in paleoclimatological studies of Arctic wedge-and pore-ice stable-water isotope records // Permafrost and Periglacial Processes. 2020. Vol. 31(3), p. 429-441. doi: 10.1002/ppp.2052.
Результаты процедуры рецензирования статьи
В связи с политикой двойного слепого рецензирования личность рецензента не раскрывается.
Со списком рецензентов издательства можно ознакомиться здесь.
Рецензия на статью «О надежности палеотемпературно-изотопных уравнений Васильчука и становлении изотопной палеогеокриологии» Предмет исследования
Предметом исследования статьи «О надежности палеотемпературно-изотопных уравнений Васильчука и становлении изотопной палеогеокриологии» является анализ истории изотопной геокриологии, надежности палеогеокриологических реконструкций, основанных на изучении изотопного состава полигонально-жильных (повторно-жильных) льдов (ПЖЛ), а также оценка вклада работ Ю. К. Васильчука в развитии данных научных направлений. Методология исследования
Основным методологическим приёмом, использованным в рецензируемой статье, является ретроспективный анализ истории развития положений изотопной геокриологии. Такой подход позволил выделить основные этапы развития данного научного направления, и роль существенную роль в его становлении Ю.К. Васильчука. Актуальность
Рассматриваемая в статье проблема достоверности и надёжности методик определения климатических характеристик по изотопным данным чрезвычайно актуально в настоящее время. Неопределённость в понимании развития природной среды в будущем, особенно важна для территорий занятых многолетнемёрзлыми породами. Открытие прямой зависимости между изотопным составом ростков современных ледяных жил и среднеянварской температурой воздуха, которое сделал Ю .К. Васильчук открыло широкую перспективу для изучения голоценовых и позднеплейстоценовых сингенетических ПЖЛ, как надежного палеоклиматического архива. Предложенные им палеотемпературно-изотопные уравнения в наибольшей мере позволяют достоверно оценить динамику изменений климата в прошлом и прогнозировать тенденции его изменения в будущем Научная новизна
Ю.К. Васильчуком впервые предложена модель циклически-пульсирующего формирования сингенетических толщ с мощными повторно-жильными льдами. Разработана стратегия радиоуглеродного датирования позднеплейстоценовых сингенетических повторно-жильных льдов и точной привязки во времени диаграмм распределения стабильных изотопов. Установлена вертикальная и латеральная
гетерохронность и гетерогенность едомных массивов, что имеет важное значение для оценки их возраста. Им выделено три типа цикличности процессов формирования синкриогенных толщ с повторно-жильными льдами: микро-, мезо и макроцикличность. Микроцикличность, связана с сезонной периодичностью изменения глубины деятельного слоя. Мезоцикличность связана с пульсирующим изменением уровня водоема, по берегам которого, или на отмелях которого идет формирование жил. Макроцикличность связана с коренной перестройкой режима седиментации или с крупными климатическими осцилляциями. Разработал стратегию датирования сингенетических многолетнемерзлых отложений с повторно-жильными льдами, предложив принцип выбора наиболее молодой датировки.
Материал изложен хорошим научным языком, структура и содержание работы позволили полностью раскрыть заявленную тему. К положительным аспектам работы можно отнести подробное изложение научных дискуссий в ведущими специалистами в которых участвовал Ю.К. Васильчук.
В рассматриваемой статье анализ результатов ведущих исследователей является одним из основных методов исследования, дан критический анализ основных работ в области изотопной геокриологии. Список использованных работ насчитывает 53 источника, это позволило не только оценить общее состояние данной отрасли науки, но и значительный вклад в её развитие работ Ю.К. Васильчука.
Апелляция к оппонентам уважительная, подробно изложены позиции различных авторов, корректно приведены доводы участников научной дискуссии по проблемам изотопной геокриологии.
Представленные материалы существенно дополняют, научные знания в области изотопной геокриологии. Приведённые в статье данные убедительно демонстрирует приоритетность российских геокриологических и исследований сингенетических отложений с повторно-жильными льдами. Подчёркивается важная роль в создании и совершенствовании этого направления Ю.К. Васильчука. Предложенные им методики создают более объективный подход как к анализу получаемых данных, так и к самому процессу полевого изучения сложно построенных массивов с полигонально-жильными льдами. Публикация будет интересна не только широкому кругу специалистов, профессионально занимающихся проблемами изотопной геокриологии, но и всем интересующимся историей развития научных взглядов в области криологических исследований.
Рассмотрение материалов изложенных в статье позволяют рекомендовать её приинять к публикации.
