Рус Eng Cn Перевести страницу на:  
Please select your language to translate the article


You can just close the window to don't translate
Библиотека
ваш профиль

Вернуться к содержанию

Арктика и Антарктика
Правильная ссылка на статью:

Возможные изменения изотопного состава повторно-жильного льда в результате самодиффузии

Васильчук Юрий Кириллович

ORCID: 0000-0001-5847-5568

доктор геолого-минералогических наук

профессор, Московский государственный университет им. М.В. Ломоносова (МГУ), географический факультет, кафедра геохимии ландшафтов и географии почв

119991, Россия, г. Москва, ул. Ленинские Горы, 1, оф. 2009

Vasil'chuk Yurij Kirillovich

Doctor of Geology and Mineralogy

Professor, Lomonosov Moscow State University, Faculty of Geography, Department of Landscape Geochemistry and Soil Geography

119991, Russia, Moscow, Leninskie Gory str., 1, of. 2009

vasilch_geo@mail.ru
Другие публикации этого автора
 

 

DOI:

10.7256/2453-8922.2021.3.36603

Дата направления статьи в редакцию:

08-10-2021


Дата публикации:

02-11-2021


Аннотация: Предметом исследования является оценка изменения изотопного состава повторно-жильного льда в ледяных жилах в результате самодиффузии. Предложен сравнительно простой метод учета влияния процессов самодиффузии на перераспределение палеоизотопного состава, который позволяет быстро, просто и с приемлемой точностью оценить изменение значений δ18О и δ2H в ледяных жилах на различные временные интервалы четвертичного этапа (1000 лет, 10000 лет, 100000 лет и 1 000000 лет). Закономерности протекания процесса диффузии в физико-химических системах описываются двумя дифференциальными формами законов диффузии Фика, которые преобразованы в интеграл Гаусса, при расчете граничных значений которого использовались шестизначные таблицы Чемберса. При расчете граничных значений полученного интеграла Гаусса использованы экспериментально определенные реальные значения коэффициента самодиффузии во льду – D=(2–10)×10-15 м2/с. Рассчитаны значения интеграла Гаусса, для ряда временных отрезков, на разных расстояниях между образцами с различной концентрацией. Для периода в 1000 лет его значения во всех случаях больше 0,99. Для более же длительных временных интервалов они существенно меняются. Расчетные вероятностные значения изменения концентрации тяжелых изотопов кислорода, в повторно-жильных льдах под действием самодиффузии на различных расстояниях и при разных градиентах концентрации заметны только для периода времени более 100000 лет.


Ключевые слова:

повторно-жильные льды, поздний плейстоцен, голоцен, многолетнемерзлые породы, изотопы кислорода, самодиффузия, законов диффузии Фика, интеграл Гаусса, изменения изотопного состава, градиент концентрации изотопов

Исследование выполнено при поддержке РНФ (проект № 19-17-00126) и Междисциплинарной научно-образовательной школы Московского государственного университета имени М.В.Ломоносова «Будущее планеты и глобальные изменения окружающей среды». This research has been supported by the Russian Scientific Foundation (grant № 19-17-00126) and Interdisciplinary Scientific and Educational School of M.V.Lomonosov Moscow State University «Future Planet and Global Environmental Change»

Abstract:   The subject of this research is the assessment of changes of the isotopic composition of ice wedge as a result of self-diffusion. The author offers a relatively simple method of considering the impact of the process of self-diffusion upon redistribution of paleo-isotopic composition that allows accurately and quickly estimating the changes in the values δ18О and δ2H values in ice wedges at different periods of the quaternary (1000 years, 10,000 years, 100,000 years, and 1,000,000 years). The patterns of diffusion process in physico-chemical systems are described by two differential forms of Fick's diffusion laws, which are transformed into Gaussian integral, using six-digit Chambers charts for calculation of the threshold values δ18О. The article applies the experimentally determined real values of self-diffusion coefficient in the ice – D=(2–10)×10-15 m2/s. The values of the Gaussian integral are calculated for time intervals, at different distances between samples with different concentration; for the period of 1,000 years its values in all cases exceed 0.99. For longer time intervals, these values change significantly. The calculated probabilistic values of changes in the concentration of heavy oxygen isotopes in ice wedge under the impact of self-diffusion at different distances and concentration gradients are noticeable only for the time period of over 100,000 years.  


Keywords:

Ice wedge, Late Pleistocene, Holocene, permafrost, oxygen isotope, self-diffusion, Fick's diffusion laws, Gaussian integral, change of the isotopic abundance ratio, isotope concentration gradient

Введение

Использование данных по вариациям значений δ18О и δ2H в древних плейстоценовых и голоценовых повторно-жильных льдах требует оценки масштабов постгенетических изменений в изотопно-кислородном составе льда, то есть к вопросу о длительности консервации первичной изотопно-кислородной «информации». Отметим, что в настоящий момент эта задача не может быть решена строго, то есть в аналитическом виде. Можно предполагать в дальнейшем достаточно корректное математическое моделирование процессов диффузии с использованием итерационных процедур, однако такой подход потребовал бы большого количества дорогостоящих прецизионных масс-спектрометрических измерений в образцах по всему объему льда, что пока не представляется реальным. Поэтому мы предлагаем сравнительно простой метод учета влияния процессов самодиффузии на перераспределение палеоизотопного состава, который позволяет быстро, просто и с приемлемой точностью оценить изменение значений δ18О и δ2H в ледяных жилах на различные временные интервалы четвертичного этапа (1000 лет, 10000 лет, 100000 лет и 1 000000 лет). Необходимо уточнить, что, хотя вопросы диффузии в различных физико-химических системах рассмотрены достаточно полно, укажем на монографии Р. Бэррера [1], В. Зайта [2], однако в палеогеографической литературе они рассматривались весьма ограниченно [3,4] , некоторые аспекты раскрыты в палеогляциологических работах [5-7], а в палеогеокриологических – не затронуты вообще.

Предлагаемый подход к расчету самодиффузии

При решении задачи диффузии каких-либо веществ (в том числе тяжелых изотопов кислорода и водорода в повторно-жильных льдах) возникает три главных аспекта: 1) поиск уравнения, дающего сравнительно простое решение, но все же с достаточной точностью аппроксимирующего реальный процесс диффузии во льду; 2) поиск значений необходимых коэффициентов, входящих в уравнения диффузии (в данном случае коэффициента самодиффузии во льду); 3) решение выбранных уравнений. Рассмотрим названные три аспекта.

Основные закономерности протекания процесса диффузии в физико-химических системах описываются двумя дифференциальными формами законов диффузии Фика:

m = –q×D×∂R/∂x (1)

∂R/∂t = D×∂2R/∂x2 (2),

где m – количество вещества, диффундирующее в единицу времени через площадь q, при коэффициенте диффузии D, расстоянии х и градиенте концентрации ∂R/∂x (R – концентрация вещества). Знак минус в правой части уравнения обозначает направление диффузии, поскольку она идет в сторону убывания концентрации.

Дифференциальные уравнения диффузии Фика, были преобразованы нами [8] в интеграл Гаусса (консультировали нас в этом д.х.н. Г.Р.Энгельгардт и д.х.н. Ю.Ф.Опруненко), при расчете граничных значений которого использовались шестизначные таблицы Чемберса [9].

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

Чтобы использовать второе уравнение для решения нашей задачи, его необходимо представить в интегральной форме с заданными определенными граничными условиями. Простейшие граничные условия можно рассмотреть в одномерной задаче с длинным цилиндром (или стержнем) с постоянным поперечным сечением, разделенным в середине плоскостью, перпендикулярной к оси цилиндра на два объема. Граничные условия и математический аппарат преобразования второго уравнения Фика из дифференциальной формы в интегральную можно отыскать, например, в работе В.Зайта [2].

Путем сравнительно несложных алгебраических действий, используя интегральное преобразование и известное выражение расчета величины δ18Oy в образце используя концентрации в образце () и в стандарте (:

δ18Oy = ( (3),

второе уравнение Фика мы представили в интегральной форме [8].

δ18Oy = (4),

где δ18Oy – относительная концентрация тяжелых изотопов кислорода в точке y в промилле, – концентрации в стандарте, относительная концентрация в объемах 1 и 2.

Нам представилось более целесообразным привести расчет не значений относительной концентрации в какой-либо точке стержня, а ее изменений через время t, т.е. Δδ18Ot.

Для точки с большими значениями концентрации (в нашей задаче положительная область цилиндра – x > 0)

Δ1δ18O = Δδ18Ot – Δδ18O1 (5),

Подставив это выражение в уравнение 4, получаем:

Δ1δ18O = ) (6),

Для точки с меньшими значениями концентрации (x < 0)

Δ2δ18O = Δδ18Ot – Δδ18O2 (7),

Подставив это выражение в уравнение 4, получаем:

Δ2δ18O = ) (8),

Нетрудно видеть, что выражения (6) и (8) являются симметричными, т.е. для нашей задачи приращение концентрации в какой-либо точке х объема 2 в момент времени tn будет равно убыванию ее в объеме 1.

Выведенный нами интеграл в правой части уравнений является известным интегралом вероятностей Гаусса, промежуточные значения которого приведены в соответствующих математических таблицах [9].

Одной из сложностей при расчете граничных значений полученного интеграла Гаусса было определение реальных значений коэффициента самодиффузии. Число экспериментальных определений коэффициента самодиффузии во льду весьма ограничено. Одна из наиболее обстоятельных публикаций принадлежит В. Куну и М. Тюркауфу [10], которые измеряли диффузию между двумя ледяными блоками, один из которых состоял из природной воды, а второй был получен из воды, обогащенной дейтерием и тяжелым кислородом. Соединенные (смерзшиеся) блоки были помещены на 4–5 недель в термостат, где поддерживалась температура от –1,5 до –2,0°С. После этого была измерена концентрация содержания тяжелых изотопов кислорода и дейтерия в обоих блоках и в разделяющей их тонкой пластинке из натурального льда. Были получены значения коэффициента самодиффузии D (равный 10×10–15 м2/с, который согласно [10] в 105 раз меньше, чем в воде. Важен и вывод о том, что диффузия во льду происходит посредством перемещения целых молекул, а не ионов и атомов его составляющих, т. о. скорость перемещения для 2Н, 1Н, 18О и 16О одинакова. Впоследствии О. Денгель и Н. Риль [11] измеряли диффузию трития в синтетических (искусственных) монокристаллах льда в диапазоне температур –10 °С – –35°С и получили величину D = 2×10-15 м2/с, а значения энергии активации процесса самодиффузии – Еа = 13,5 ± 1 ккал/моль. К. Итагаки [12, 13] провел весьма интересное исследование самодиффузии во льду с помощью паров, насыщенных тритием (получаемых при испарении тритированной воды с удельной активностью 30 милликюри на мл.) в образцах льда, привезенных с ледника Менденхалле (Аляска). Продолжительность опытов составляла от 1 до нескольких сотен дней. Замеры велись последовательно в слоях, снимаемых с помощью микротома, при этом, толщина их составляла 10 микрон. Полученный таким способом коэффициент самодиффузии при температуре –10°С составил (1,7–7,8) ×10-15 м2/с (столь широкий диапазон приведенных значений D объясняется тем, что это диффузия и вдоль длинной оси кристаллов и вкрест ее простирания, средняя же величина коэффициента самодиффузии, если его рассчитать по приведенной в [12] энергии активации Еа = 15×7 ± 2 ккал/моль, составил D=2,51×10-15 м2/с). Для температуры –30 °С ее оказалось невозможным измерять, ввиду слишком малой скорости диффундирования трития при этой температуре. Можно констатировать, что температурный режим в котором проводились замеры диффузии во льду, близки к природным, свойственным многим районам Субарктики и поэтому указанные значения можно считать достаточно приемлемым для решения задач изотопной палеогеокриологии. Предполагая, что недооценка самодиффизии может существенно исказить результаты наших палеопостроений, мы сочли необходимым рассчитать изменение концентрации для всего диапазона значений коэффициента самодиффузии, вернее, для его крайних значений D=(2–10)×10-15 м2/с. Эти значения коэффициента самодиффузии было использовано нами в расчетах граничных значений интеграла Гаусса.

Результаты

Мы рассчитывали величину самодиффузии с применением интеграла вероятностей Гаусса, промежуточные значения которого определяются с помощью математических таблиц Чемберса, значения которых взяты из работы Комри [9]. Использованы значения коэффициента самодиффузии, определенные лабораторным путем.

В табл. 1 приведены значения интеграла Гаусса (определенные с помощью таблиц Чемберса [9], входящего в выражение (8) при коэффициенте D=(2–10)×10-15 м2/с, для ряда временных отрезков, на разных расстояниях между образцами с различной концентрацией. Для периода в 1000 лет его значения во всех случаях больше 0,99. Для более же длительных временных интервалов они существенно меняются.

Таблица 1. Пределы изменения интеграла Гаусса при коэффициентах самодиффузии, изменяющихся от 2×10-15 (первая цифра в таблице) до 10×10-15 м2/сек (вторая цифра)

Расстояние, м

Период годы

10000

100000

1000000

0,1

0,79–1,00

0,31–0,63

0,10–0,22

0,2

0,99–1,00

0,57–0,93

0,20–0,43

0,3

1,00

0,77–0,99

0,30–0,60

0,4

1,00

0,80–1,00

0,39–0,74

0,5

1,00

0,95–1,00

0,48–0,84

0,6

1,00

0,98–1,00

0,55–0,91

0,7

1,00

0,99–1,00

0,62–0,95

0,8

1,00

1,00

0,68–0,98

0,9

1,00

1,00

0,74–0,99

1,0

1,00

1,00

0,79–1,00

В табл. 2 (с использованием данных табл. 1 и значений полученных при подстановке в интегральные уравнения Гаусса) приведены расчетные данные изменения значений δ18О в жильных льдах под действием самодиффузии для различных расстояний между точками с реально встреченными значениями δ18О в повторно-жильных льдах с запасом, укладывающимся в диапазоне 0–45‰, что соответствует диапазону абсолютной концентрации (т. е. отношению 18О/16О= (1,97–1,90)×10-3.

Таблица 2. Расчетные вероятностные значения изменения концентрации тяжелых изотопов кислорода (δ18О, ‰), в повторно-жильных льдах под действием самодиффузии на различных расстояниях и при разных градиентах концентрации

Гра-диент

Период, лет

10000

100000

1000000

Расстояние, м

1,0

0,5

0,1

1,0

0,5

0,1

1,0

0,5

0,1

–5

0

0

0–0,53

0

0–0,13

0,93–1,73

0–0,58

0,4–1,3

1,95–2,25

–10

0

0

0–1,05

0

0–0,25

1,85–3,45

0–1,06

0,8–2,60

3,90–4,50

–15

0

0

0–1,58

0

0–0,38

2,78–5,18

0–1,58

1,2–3,9

5,85–6,75

–20

0

0

0–2,10

0

0–0,50

3,70–6,90

0–2,12

1,6–5,20

7,80–9,00

–25

0

0

0–2,63

0

0–0,63

4,63–8,63

0–2,63

2,0–6,5

9,75–11,25

–30

0

0

0–3,15

0

0–0,75

5,55–10,35

0–3,16

2,40–7,80

11,90–13,50

–35

0

0

0–3,68

0

0–0,88

6,48–12,08

0–3,68

2,8–9,1

13,65–15,75

–40

0

0

0–4,20

0

0–1,00

7,40–13,80

0,4,24

3,2–10,40

15,60–18,00

–45

0

0

0–4,73

0

0–1,13

8,33–15,53

0–4,73

3,6–11,7

17,55–20,25

Примечания: значение градиента дано в виде разницы между точками с меньшим значением δ18О и точкой с большим ее значением; в тех случаях, когда приведено значение 0, изменения составляют менее 2‰; первая цифра в таблице - при коэффициенте самодиффузии D=2×10-15 м2/с, вторая при коэффициенте самодиффузии D=10×10-15 м2/с.

Проведенные расчеты показали, что заметное изменение концентрации стабильных тяжелых изотопов кислорода (которые удобнее выражать через Δδ18О = δ118О— δ218О) в повторно-жильных льдах происходят весьма медленно.

Для периода в 10000 лет изменения будут сколь-нибудь существенными лишь на расстоянии 0,1 м, во всех остальных случаях они не составят и 0,02‰ (табл. 1). Даже за период 100 000 лет на расстоянии 1,0 м при максимально возможных (точнее, мыслимых) градиентах концентрации между двумя точками в теле ледяных жил происходит перераспределение значений δ18О не более, чем на 0,02‰ (см. табл. 2).

Несколько большие изменения первичной концентрации за 100 000 лет на расстоянии 0,5 м. Они уже заметно превосходят ошибку масс-спектрометрических измерений, составляя при градиенте 20‰ около 0,5‰ (см. табл. 2). На более близких расстояниях даже при градиенте в 5‰ перераспределение тяжелых изотопов кислорода за этот период настолько ощутимо, что может очень сильно исказить характер первичного их распределения и повлиять на достоверность палеореконструкций.

Это подтверждают и данные натурных наблюдений. В трех жилах разного возраста на севере Якутии – более 40 тысяч лет (обнажение Дуванный Яр), 6–3 тысяч лет и менее 3 тысяч лет – нами были отобраны образцы по горизонтали через 3–5 сантиметров. Отбор был произведен так, что каждый образец включал не более 4–7 элементарных жилок. Проведенные измерения дали следующие значения δ18О: в позднеплейстоценовой жиле Дуванного Яра подряд слева направо –31,3; –32,4; –32,1; –32,2; –31,9; –31,1 ‰; то есть здесь сохранились различия в изотопном составе соседних частей жилы, разделенных 3–5 см при градиентах, превышающих 1‰. В голоценовой жиле (в толще аласа) возрастом 6–8 тысяч лет аналогичное опробование дало следующие результаты: –25,1; –26,1; –27,0; –26,1; –26,4; –26,2‰, то есть здесь тоже соседствуют образцы, различающиеся на 1‰ по значениям δ18О. В жиле (в толще поймы р. Колымы) возрастом менее 3 тысяч лет замеры таковы: –27,9; –27,6; –27,0; –26,7‰ [8].

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

Дискуссия

На вопрос о возможном обмене реликтовых жил с современной средой воздушным путем по трещинам достаточно определенный ответ получен при исследованиях содержания трития в мерзлой породе. Ф. Майкл и П. Фритц [14] исследовали тритий в породах из шести глубоких 6–22 м скважин, причем в одной из них отбор произведен с интервалом 2–3 см. Во всех скважинах тритий был обнаружен только в верхней 3-метровой толще; ниже залегающие породы трития не содержали и следовательно залегают вне зоны активного обмена. Резкое снижение содержания трития было обнаружено Ф. Майклом [15] и на подошве сезонно-талого слоя в голоценовых толщах Юкона, сходные данные получены Д. В. Михалевым, исследовавшим концентрацию трития в мерзлых толщах Нижнеколымского района [16]. Исследования трития в различных генетических типах подземных льдов Южной Якутии показали, что в отличие от формирующихся в настоящее время повторно-жильных льдов, концентрации трития в которых превышает 230 ТЕ, реликтовые сингенетические льды (даже относительно молодые - голоценовые) и вмещающие их мерзлые породы содержат тритий в незначительных количествах (менее 15 ТЕ), свидетельствуя о том, что в настоящее время древние жилы находятся в стадии консервации, и современные морозобойные процессы их не затрагивают [17].

Заметим, что на отсутствие взаимодействия большинства жил с современной средой указывает прежде всего однородность их текстуры по вертикали, в подавляющем числе разрезов отсутствуют признаки более поздних внедрений. Различие в изотопном составе голоценовых и плейстоценовых жил – также признак весьма надежный, если бы взаимодействие происходило, то за столь длительное время эти сравнительно небольшие (в геологическом и не в палеотемпературном плане) различия несомненно бы выровнялись. Все эти факты указывают на длительную консервацию реликтовых сингенетических жил и на сохранность первичного распределения в них стабильных изотопов [18].

Если перераспределение изотопно-кислородного состава в повторно-жильных льдах происходит в результате самодиффузии (при отсутствии иных причин, побуждающих к миграции изотопов), то оно ведет только к уравниванию значений δ18О и δ2H по всему объему жил, поэтому если к моменту определения сохранились существенные различия значений δ18О и δ2H по вертикали или по горизонтали жил, то можно безусловно утверждать, что лед формировался в различных температурных условиях.

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

Для голоценовых жил отбор образцов можно проводить практически через любые интервалы, и процессами диффузии при палеогеокриологической оценке колебаний значений δ18О и δ2H можно во всех случаях пренебрегать. Для позднеплейстоценовых жил оптимальным шагом опробования является 1,0 м, при этом диффузию можно также не учитывать, при шаге менее 0,5 м уже необходимо допускать, что некоторое выравнивание значений между соседними образцами произошло.

В наземных и подземных льдах по-разному проявляются процессы, ведущие к преобразованию первичного изотопного сигнала. Изотопная “летопись” в подземных льдах сохраняется почти неизменной на протяжении десятков тысяч лет, а в наземных льдах первичная сезонная изотопная дифференциация исчезает через 10-15 тыс. лет, хотя полного выравнивания изотопного состава не происходит на протяжении сотен тысяч лет.

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

Если перераспределение изотопно-кислородного состава в повторно-жильных льдах происходит в результате самодиффузии (при отсутствии других причин, побуждающих к миграции изотопов), оно ведет к уравниванию содержания 18O и 2H по всему объему жил. Если к моменту определения сохранились существенные различия значений δ2H и δ18О по вертикали или горизонтали жил, можно утверждать, что лед формировался в разных температурных условиях.

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

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

Отметим недавнюю работу В.В. Малаховой и А.В. Елисеева [19], изучивших влияние диффузии солей на состояние и распространение многолетнемёрзлых пород и зоны стабильности метан-гидратов шельфа моря Лаптевых. Учёт диффузии соли в задачах исследования динамики субаквальной мерзлоты необходим при определении положения верхней границы субаквальной мерзлоты, а также при расчёте скорости её деградации. Перенос солей может изменить положение и скорость смещения верхней границы многолетнемёрзлых пород в несколько раз по сравнению со случаем неизменной во времени солёности и, следовательно, постоянной во времени температурой замерзания. Вместе с тем В.В. Малахова и А.В. Елисеев [19] отметили, что перенос солей незначительно влияет на положение нижней границы многолетнемёрзлых пород.

А.А. Екайкин и В.Я. Липенков с соавторами подчеркивают [20], что возможность реконструкции климатических событий плейстоцена по образцам ледяного керна определяется не только возрастом слагающего его льда, но и степенью сохранности климатического сигнала, который подвергается существенному диффузионному сглаживанию в придонных, наиболее древних слоях ледника. А.А. Екайкин и др. [20] уделили особое внимание моделированию процесса диффузии молекул воды во льду и оценке степени ослабления климатического сигнала в изотопном профиле ледяного керна, который может быть получен при осуществлении проекта глубокого бурения Антарктического ледника.

Модельные расчёты для Купола В выполненные в работе [20] позволили рассчитать распределение длины диффузии для различных сценариев распределения возраста и температуры в леднике. В работе [20] показано, что значение начальной длины диффузии практически не влияет на итоговую длину диффузии в древнем льду. Для всех сценариев длина диффузии быстро убывает при увеличении возраста льда до 100 тыс. лет за счёт утончения слоёв. Длина диффузии тесно связана с температурой, на величину этого параметра первостепенное влияние оказывает геотермальный поток тепла (Geothermal Heat Flux – GHF). Для основного сценария длина диффузии достигает 3 см при возрасте льда 1,4 млн лет. Максимальная длина диффузии, которая может быть достигнута при сценарии 4, когда температура льда на ложе достигает точки плавления, составляет 4 см [20].

В работе [21] показано, что для ледяного керна, из скважины пробуренной на Куполе С, характерна «дополнительная диффузия», при которой наблюдаемое диффузионное сглаживание за счёт неустановленных пока причин сильнее того, которое предполагается моделью диффузии. Для керна скважины 5Г на станции Восток было показано, что в условиях станции Восток это дополнительное сглаживание невелико и составляет не более 30% [20] расчётной величины длины диффузии. С учётом этого фактора максимально возможная оценка длины диффузии в древнем льду Купола В будет равна 5,2 см. В плейстоцене длительность основного климатического цикла составляла около 40 тыс. лет [22]. Во льду возрастом 1,4 млн лет на глубине 2440 м в районе Купола В временнóй интервал длиной 40 тыс. лет заключён в слое льда толщиной около 290 см. При длине волны 290 см и длине диффузии 5,2 см ослабление амплитуды сигнала составляет около 0,6%. Таким образом, при отсутствии складчатости и перемешивании слоёв в ледяных кернах должен сохраняться ненарушенный климатический сигнал возрастом до 1,4 млн лет [20].

Существенно влияют на изотопный состав снежной толщи и верхних частей ледников процессы фирнизации. По расчетам С.Йонсена [23] Johnsen, 1977], суммарное среднее расстояние распространения диффузии составляет около 80 мм и почти не зависит от температуры и скорости аккумуляции. Это уравнение показывает, что амплитуда в годовых слоях льда толщиной 0,6; 0,3 и 0,2 м понижается соответственно на 30, 75 и 96%. Таким образом годовые циклы δ2H и δ18О мало сглаживаются в процессе фирнизации и длительное время сохраняются во льду. Вводя модель течения, температурный профиль и диффузионную константу для условий ст. Кемп Сенчури С. Йонсен [23] получил, что для уменьшения амплитуды А на 10% потребуется 5 тыс. лет, на 50% – 8 тыс. лет, на 90% – 10 тыс. лет.

Интерпретация сезонных вариаций δ18О в верхней части ледяного керна со ст. Бэрд в Западной Антарктиде гораздо сложнее, чем для керна Кемп Сенчури, поскольку они постепенно исчезают в процессе фирнизации и не отслеживаются в плотном льду.

Выводы

Масштаб возможных изменений изотопного состава в повторно-жильных льдах в результате самодиффузии зависит прежде всего от возраста жил и от градиента концентрации стабильных изотопов кислорода.

1. В повторно-жильных льдах возраст которых составляет 10000 лет изменения изотопного состава крайне незначительны и при реально возможных градиентах существ в 5‰ на расстоянии 0,1 м изменения концентрации тяжелых изотопов кислорода δ18О не превысят 0,5‰.

2. В повторно-жильных льдах возраст которых составляет 100000 лет изменения изотопного состава незначительны и при реально возможных градиентах в 5‰ на расстоянии 0,1 м изменения концентрации тяжелых изотопов кислорода δ18О могут составить 0,93-1,73‰.

3. В повторно-жильных льдах возраст которых составляет 100000 лет изменения изотопного могут быть достаточно заметны и при реально возможных градиентах в 5‰ на расстоянии 0,1 м изменения концентрации тяжелых изотопов кислорода δ18О превысят 1,95-2,25‰, т.е приведут к заметному выравниванию изотопного состава жил.

Библиография
1. Бэррер Р.М. Диффузия в твердых телах. Пер. с англ. З. С. Суворовой; Под ред. Б. Д. Тазулахова. Москва : Гос. изд-во иностр. лит., 1948 504 с.
2. Зайт В. Диффузия в металлах. Процессы обмена мест. Пер. со второго немецк. издания Г.С. Куликова и Р. Малковича, под ред. Б.И. Болтакса. М.: Гос. изд-во иностр. лит., 1958 382 с.
3. Боуэн Д. Четвертичная геология: Стратиграфическая основа междисциплинарных исследований. Пер. с англ. И.И. Спасской; предисл. Л. Серебрянного. М.: Мир. 1981, 272 с.
4. Юри Г., Лоуэнштам Г., Эпстайн С., Мак-Кинни К. Определение палеотемператур, в частности температур верхнего мела Англии, Дании и Юго-Восточных штатов США // Изотопы в геологии. М.: Гос. изд-во иностр. лит., 1954, с.543-572.
5. Dansgaard W,. Johnsen S.J., Møller J., Langway C.C., Jr. One thousand centuries of climatic record from Camp Century on the Greenland Ice Sheet // Science. 1969. Vol. 166, N3903. P.377–381.
6. Johnsen, S.J., Dansgaard W., Clausen H.B., Langway C.C., Jr. Oxygen isotope profiles through the Antarctic and Greenland ice sheets // Nature, 1972. Vol. 235 (5339), p. 429-434,
7. Johnsen S., Clausen, H. B., Cuffey K.M., Hoffmann G., Schwander J., Creyts T. Diffusion of stable isotopes in polar tirn and ice: the isotope effect in tirn diffusion // Physics of Ice Core Records. Edited by T. Hondoh Hokkaido University Press, Sapporo, 2000, p. 121-140.
8. Васильчук Ю.К. Изотопно-кислородный состав подземных льдов (опыт палеогеокриологических реконструкций). М.: Изд. Отдел. Теоретических проблем РАН. МГУ, ПНИИИС. 1992. В 2–х томах. Т.1. – 420 с. Т.2 – 264 с. (Vasil'chuk Yu. K. 1992: Oxygen isotope composition of ground ice (application to paleogeocryological reconstructions). Volume 1, 420 pp., Volume 2, 264 pp. Theoretical Problems Department, Russian Academy of Sciences and Lomonosov Moscow University Publications, Moscow (in Russian with English contents section).
9. Комри Л.Дж.Шестизначные математические таблицы Чемберса. М.: Наука. 1964. 574 с
10. Kuhn W., Thürkauf M. Isotopentrennung beim Gefrieren von Wasser und Diffusionskonstanten von D and im Eis. Mit Diskussion der Möglichkeit einer Multiplikation der beim Gefrieren auftretenden Isotopentrennung in eine Haarnadelgegenstromvorrichtung // Helvetica Chimica Acta. 1958. Vol. 41. N4. S.938–971.
11. Dengel O., Riehl N., Diffusionen von Protonen (Tritonen) in Eiskristallen. Phys. Kondens. Materie, 1963, vol. 1, p. 191-196
12. Itagaki K. Self-Diffusion in Single Crystals of ice // J. of the Physical Society of Japan. 1964. Vol. 19. N6. P.1081.
13. Itagaki K. Self-Diffusion in Single Crystal Ice // Journal of the Physical Society of Japan, 1967, vol. 22, no. 2 : p. 427-431. doi: 10.1143/JPSJ.22.427
14. Michel F.A., Fritz P. Environmental isotopes in permafrost related waters along the Mackenzie valley corridor // Permafrost. Third International Conference, Proceedings. Edmonton. Alberta, Canada. Ottawa: National Research Council of Canada, 1978. P.207–211
15. Michel F.A. Isotope geochemistry of frost-blister ice, North Fork Pass, Yukon, Canada // Canadian J. of Earth Sciences. 1986. Vol. 23. N4. P.543–549.
16. Михалев Д.В. Изотопно-кислородный анализ текстурообразующих льдов (на примере Колымской низменности и Енисейского Севера) / Автореферат дисс. на соиск. уч. степ. канд. геогр. наук. М. 1990. 24 с.
17. Морковкина И.К., Романов Б.В., Чижов А.Б., Чижова Н.И. Исследование содержания трития в мерзлых породах и подземных льдах // Вод. ресурсы. 1982. № 5. С. 92–97.
18. Vasil'chuk, Yu.K., Yesikov, A.D., Oprunenko, Yu.F., Petrova, Ye.A., Vasil'chuk, A.C. and Sulerzhitskiy, L.D. 1985. New data of stable oxygen isotopes composition in syngenetic Late Pleistocene ice wedge of the lower Kolyma River. Transactions (Doclady) of the USSR Academy of Sciences. Earth Science Sections. Published by Scripta Technica, Inc. A Wiley Company. New York. Vol. 281. N2. P. 91–94.
19. Malakhova V.V., Eliseev A.V. Salt diffusion effect on the submarine permafrost state and distribution as well as on the stability zone of methane hydrates on the Laptev Sea shelf. Led i Sneg. Ice and Snow. 2020. 60 (4): 533–546. [In Russian]. doi: 10.31857/S2076673420040058.
20. Ekaykin A.A, Lipenkov V.Ya., Tchikhatchev K.B. Preservation of the climatic signal in the old ice layers at Dome B area (Antarctica). Led i Sneg. Ice and Snow. 2021. 61 (1): 5–13. [In Russian]. doi: 10.31857/S2076673421010067.
21. Pol K., Masson-Delmotte V., Johnsen S., Bigler M., Cattani O., Durand G., Falourd S., Jouzel J., Minster B.,Parrenin F., Ritz C., Steen-Larsen H.C., Stenni B. New MIS 19 EPICA Dome C high resolution deuterium data: Hints for a problematic preservation of climate variability at sub-millennial scale in the 'oldest ice'// Earth and Planetary Science Letters. 2010, 298: 95–103.
22. Berends C.J., de Boer B., van de Wal R.S.W. Reconstructing the Evolution of Ice Sheets, Sea Level and Atmospheric CO2 During the Past 3.6 Million Years. Clim. Past Disc. 2020: 1–22. https://doi.org/10.5194/cp-2020-52.
23. Johnsen S.J. Stable isotope homogenization of polar firn and ice // Isotopes and Impurities in Snow and Ice.1977, vol. 118, p. 210–219.
References
1. Berrer R.M. Diffuziya v tverdykh telakh. Per. s angl. Z. S. Suvorovoi; Pod red. B. D. Tazulakhova. Moskva : Gos. izd-vo inostr. lit., 1948 504 s.
2. Zait V. Diffuziya v metallakh. Protsessy obmena mest. Per. so vtorogo nemetsk. izdaniya G.S. Kulikova i R. Malkovicha, pod red. B.I. Boltaksa. M.: Gos. izd-vo inostr. lit., 1958 382 s.
3. Bouen D. Chetvertichnaya geologiya: Stratigraficheskaya osnova mezhdistsiplinarnykh issledovanii. Per. s angl. I.I. Spasskoi; predisl. L. Serebryannogo. M.: Mir. 1981, 272 s.
4. Yuri G., Louenshtam G., Epstain S., Mak-Kinni K. Opredelenie paleotemperatur, v chastnosti temperatur verkhnego mela Anglii, Danii i Yugo-Vostochnykh shtatov SShA // Izotopy v geologii. M.: Gos. izd-vo inostr. lit., 1954, s.543-572.
5. Dansgaard W,. Johnsen S.J., Møller J., Langway C.C., Jr. One thousand centuries of climatic record from Camp Century on the Greenland Ice Sheet // Science. 1969. Vol. 166, N3903. P.377–381.
6. Johnsen, S.J., Dansgaard W., Clausen H.B., Langway C.C., Jr. Oxygen isotope profiles through the Antarctic and Greenland ice sheets // Nature, 1972. Vol. 235 (5339), p. 429-434,
7. Johnsen S., Clausen, H. B., Cuffey K.M., Hoffmann G., Schwander J., Creyts T. Diffusion of stable isotopes in polar tirn and ice: the isotope effect in tirn diffusion // Physics of Ice Core Records. Edited by T. Hondoh Hokkaido University Press, Sapporo, 2000, p. 121-140.
8. Vasil'chuk Yu.K. Izotopno-kislorodnyi sostav podzemnykh l'dov (opyt paleogeokriologicheskikh rekonstruktsii). M.: Izd. Otdel. Teoreticheskikh problem RAN. MGU, PNIIIS. 1992. V 2–kh tomakh. T.1. – 420 s. T.2 – 264 s. (Vasil'chuk Yu. K. 1992: Oxygen isotope composition of ground ice (application to paleogeocryological reconstructions). Volume 1, 420 pp., Volume 2, 264 pp. Theoretical Problems Department, Russian Academy of Sciences and Lomonosov Moscow University Publications, Moscow (in Russian with English contents section).
9. Komri L.Dzh.Shestiznachnye matematicheskie tablitsy Chembersa. M.: Nauka. 1964. 574 s
10. Kuhn W., Thürkauf M. Isotopentrennung beim Gefrieren von Wasser und Diffusionskonstanten von D and im Eis. Mit Diskussion der Möglichkeit einer Multiplikation der beim Gefrieren auftretenden Isotopentrennung in eine Haarnadelgegenstromvorrichtung // Helvetica Chimica Acta. 1958. Vol. 41. N4. S.938–971.
11. Dengel O., Riehl N., Diffusionen von Protonen (Tritonen) in Eiskristallen. Phys. Kondens. Materie, 1963, vol. 1, p. 191-196
12. Itagaki K. Self-Diffusion in Single Crystals of ice // J. of the Physical Society of Japan. 1964. Vol. 19. N6. P.1081.
13. Itagaki K. Self-Diffusion in Single Crystal Ice // Journal of the Physical Society of Japan, 1967, vol. 22, no. 2 : p. 427-431. doi: 10.1143/JPSJ.22.427
14. Michel F.A., Fritz P. Environmental isotopes in permafrost related waters along the Mackenzie valley corridor // Permafrost. Third International Conference, Proceedings. Edmonton. Alberta, Canada. Ottawa: National Research Council of Canada, 1978. P.207–211
15. Michel F.A. Isotope geochemistry of frost-blister ice, North Fork Pass, Yukon, Canada // Canadian J. of Earth Sciences. 1986. Vol. 23. N4. P.543–549.
16. Mikhalev D.V. Izotopno-kislorodnyi analiz teksturoobrazuyushchikh l'dov (na primere Kolymskoi nizmennosti i Eniseiskogo Severa) / Avtoreferat diss. na soisk. uch. step. kand. geogr. nauk. M. 1990. 24 s.
17. Morkovkina I.K., Romanov B.V., Chizhov A.B., Chizhova N.I. Issledovanie soderzhaniya tritiya v merzlykh porodakh i podzemnykh l'dakh // Vod. resursy. 1982. № 5. S. 92–97.
18. Vasil'chuk, Yu.K., Yesikov, A.D., Oprunenko, Yu.F., Petrova, Ye.A., Vasil'chuk, A.C. and Sulerzhitskiy, L.D. 1985. New data of stable oxygen isotopes composition in syngenetic Late Pleistocene ice wedge of the lower Kolyma River. Transactions (Doclady) of the USSR Academy of Sciences. Earth Science Sections. Published by Scripta Technica, Inc. A Wiley Company. New York. Vol. 281. N2. P. 91–94.
19. Malakhova V.V., Eliseev A.V. Salt diffusion effect on the submarine permafrost state and distribution as well as on the stability zone of methane hydrates on the Laptev Sea shelf. Led i Sneg. Ice and Snow. 2020. 60 (4): 533–546. [In Russian]. doi: 10.31857/S2076673420040058.
20. Ekaykin A.A, Lipenkov V.Ya., Tchikhatchev K.B. Preservation of the climatic signal in the old ice layers at Dome B area (Antarctica). Led i Sneg. Ice and Snow. 2021. 61 (1): 5–13. [In Russian]. doi: 10.31857/S2076673421010067.
21. Pol K., Masson-Delmotte V., Johnsen S., Bigler M., Cattani O., Durand G., Falourd S., Jouzel J., Minster B.,Parrenin F., Ritz C., Steen-Larsen H.C., Stenni B. New MIS 19 EPICA Dome C high resolution deuterium data: Hints for a problematic preservation of climate variability at sub-millennial scale in the 'oldest ice'// Earth and Planetary Science Letters. 2010, 298: 95–103.
22. Berends C.J., de Boer B., van de Wal R.S.W. Reconstructing the Evolution of Ice Sheets, Sea Level and Atmospheric CO2 During the Past 3.6 Million Years. Clim. Past Disc. 2020: 1–22. https://doi.org/10.5194/cp-2020-52.
23. Johnsen S.J. Stable isotope homogenization of polar firn and ice // Isotopes and Impurities in Snow and Ice.1977, vol. 118, p. 210–219.

Результаты процедуры рецензирования статьи

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

Рецензия
на статью «Возможные изменения изотопного состава повторно-жильного льда в результате самодиффузии»

Предмет исследования статьи являются оценка перераспределения значений δ18О и δ2H в древних плейстоценовых и голоценовых повторно-жильных льдах за счёт диффузии.
В рассматриваемой статье применён комплексный анализ процессов перераспределение изотопного состава δ18О и δ2H повторно-жильных льдах за счёт диффузии. В этот комплекс включены преобразования дифференциальных уравнений диффузии Фика, в интеграл Гаусса, при расчете граничных значений которого использовались шестизначные таблицы Чемберса. Затем были использованы результаты опытов проведённых В. Куном и М. Тюркауфом, по измерению диффузии между двумя ледяными блоками, один из которых состоял из природной воды, а второй был получен из воды, обогащенной дейтерием и тяжелым кислородом. Эти значения коэффициента диффузии были использовано в расчетах граничных значений интеграла Гаусса. Затем был произведён расчёт величины изменения значений δ18О в жильных льдах под действием самодиффузии для различных расстояний между точками с реально встреченными в природе значениями δ18О. Проведенные расчеты показали, что заметное изменение концентрации стабильных тяжелых изотопов кислорода в повторно-жильных льдах происходят весьма медленно. И в заключении расчётные значения диффузии сравнивались с данными натурных наблюдений.
Проблема оценки масштабов постгенетических изменений в изотопно-кислородном составе льда весьма актуальна. Данный метод широко применяется в геокриологии для оценки времени и палеогеогафических условий формирования повторно-жильных льдов. Поэтому проблема оценки интенсивности перераспределения и выравнивания изотопного состава за счёт диффузии в повтороно-жильных льдах, является ключевой. Недооценка самодиффизии может существенно исказить результаты палеогеографических построений, Исследователю необходимо знать, на сколько может изменяться первичная изотопно-кислородной «информация» во времени, в условиях неравномерного соотношения δ18О и δ2H в жиле. В статье справедливо отмечается, что в настоящий момент эта задача не может быть решена строго, то есть в аналитическом виде. Математическое моделирование процессов диффузии с использованием итерационных процедур, требует большого количества дорогостоящих прецизионных масс-спектрометрических измерений в образцах по всему объему льда, что в полевых условиях осуществить затруднительно.
Исследований диффузного выравнивания изотопного состава в подземных льдах, помимо работ автора практически не проводилось. Известны лишь работы по диффузным изменениям химического состава в многолетнемерзлых толщах. В статье предлагаем сравнительно простой метод учета влияния процессов самодиффузии на перераспределение палеоизотопного состава, который позволяет быстро, просто и с приемлемой точностью оценить изменение значений δ18О и δ2H в ледяных жилах на различные временные интервалы четвертичного этапа (1000 лет, 10000 лет, 100000 лет и 1000000 лет). Чрезвычайно важными и новыми выводами статьи являются утверждения о том, что в повторно-жильных льдах, возраст которых составляет 100000 лет изменения изотопного состава и диффузионное сглаживание незначительны. При реально возможных градиентах в 5‰ на расстоянии 0,1 м изменения концентрации тяжелых изотопов кислорода δ18О могут составить 0,93-1,73‰. Только при достижения возраста повтороно-жильных льдов около 1000000 лет изменения изотопного могут быть достаточно заметны и при реально возможных градиентах в 5‰ на расстоянии 0,1 м изменения концентрации тяжелых изотопов кислорода δ18О превысят 1,95-2,25‰, т.е приведут к заметному выравниванию изотопного состава жил. Приведённые данные позволяют найти интервал возраста повтороно-жильных льдов с допустимыми значениями диффузионного сглаживания.
Результаты, приводимые в статье, показывает её высокий научный уровень. Работа написана понятным для специалиста языком. Стиль изложения соответствует принятому для научных статей и свидетельствует о высоком уровне и широком научном кругозоре авторов. Структура работы логичная, изложение последовательное. Содержание статьи отражает проведённые исследования и позволяет сравнить полученные результаты с данными других исследователей.
Список источников составляет 23 наименования. Автором проделана большая работа по анализу материалов предыдущих специалистов от первых классических работ по данной тематике до современных последних публикаций. Широко представлены работы зарубежных исследователей.
Расчеты и полевые наблюдения за диффузией стабильных изотопов указывают на длительную консервацию реликтовых сингенетических жил и сохранность первичного распределения в них стабильных изотопов, что позволяет предполагать корректность палеотемпературных реконструкций, выполняемых по повторно-жильным льдам, возраст которых не превышает 100 тыс. лет. Когда же дело касается более древних льдов, с процессами самодиффузии уже необходимо считаться, по крайней мере вводя новые пределы точности.
Приведенные в статье материалы привлекут внимание специалистов геокриологов, палеогеографов, экологов, почвоведов, климатологов, физико-географов, а также всех интересующихся проблемами формирования криолитозоны.