DOI: 10.7256/2453-8922.2019.4.31645
Received:
08-12-2019
Published:
17-12-2019
Abstract:
A review of published isotopic data on Yamal Peninsula massive ice near Marre-Sale polar station showed that there is a number of difficulties in indicating the genesis of the ice formation. The published data of δ18О and δ2Н of sea water in the Yamal Peninsula coast are disconnect to idea of participation of this sea water as a water source for ice formation. Comparison with the correct δ18О and δ2Н values of waters of the Barents Sea made it possible to establish the process of mixing seawater with isotopically light precipitation and showed that seawater was involved in the formation of ice. Another difficulty was the erroneous idea of ice formation in a closed system. Formation of ice occurred in an open system with a large amount of free water. Massive ice was most likely formed at the clay-sand boundary from free water, which in large quantities entered the sandy aquifer from lake talik. The water in the lake at the first stage was a mixture of sea and atmospheric waters, and then it was desalinated. The process of desalination and isotope depletion of lake water is responsible for a decrease in the δ18О and δ2Н values of ice with depth.
Keywords:
massive ice, Marresale, oxygen isotopes, hydrogen isotope, initial water, open system, closed system, deuterium excess, lake water, sea water
Введение
Разрез четвертичных отложений в районе полярной станции Марре-Сале, расположенной на западном побережье полуострова Ямал, является достаточно хорошо изученным [1 - 13]. Для этого разреза были описаны многочисленные дислокации вмещающих отложений (рис. 1), сам разрез неплохо датирован по радиоуглероду растительных остатков и торфа.
Рис. 1. Дислокации, складчатые структуры в береговых отложениях близ Марре-Сале. Фото О.Инголфссона
Однако остается ряд ключевых вопросов, по которым возникают разногласия. В первую очередь это касается условий образования мощных залежей пластовых льдов (ПЛ) и дислокаций. Толща четвертичных отложений, вскрываемая в обнажениях близ Марре-Сале, состоит из двух основных частей. Верхняя часть представлена незасоленными континентальными отложениями. Для нижней, более древней толщи, генезис и возраст осадков остается дискуссионным. К нижней толщи относятся в разной степени засоленные отложения прибрежно-морского и морского генезиса [12; 8]. В верхнем континентальном комплексе отложений описаны пластовые льды (рис. 2, а) с лакколитами и штоками [5]. Пластовые льды в нижнем комплексе отложений (рис. 2, б) имеют видимую мощность до 5 м и залегают согласно вмещающим породам [12; 5; 10].
Рис. 2. Два типа пластового льда в обнажении Марре-Сале: а – верхний, из Е.А.Слагоды и др. [5, c. 14], б – нижний, из Г.Е. Облогова [10, с.91]
В данной статье мы рассмотрим опубликованные изотопные данные по пластовым залежам Марре-Сале [5, 7, 10] и проанализируем основные гипотезы формирования льда с позиций изотопной индикации. Целью работы является анализ опубликованных данных изотопного состава кислорода и водорода льда пластовых залежей, а также рассмотрение основных противоречий и трудностей, возникающих при интерпретации изотопных данных. Опубликованные данные изотопного состава кислорода и водорода пластового льда Марре-Сале В обнажениях Марре-Сале было описано два разных типа пластового льда, который назван «пластовый лед 1 (ПЛ1)» и «пластовый лед 2 (ПЛ2)» [7, 10] или «верхняя залежь льда» и «нижняя залежь льда» [5]. Горизонт ПЛ1 представляет собой чередование субгоризонтальных, волнистых, вплоть до субвертикально-изогнутых (дислоцированных) прослоев льда толщиной от нескольких сантиметров до 1,5 м и прослоев грунта, которые представлены песками и супесями с суммарной льдистостью 40–60%. В 4,5 км к югу от п/с Марре-Сале в береговом обрыве высотой 13–15 м под слоем каргинско-сартанского песка и казанцевских глин с глубины около 5 м вскрывается толща ПЛ2 мощностью более 8 м (лед уходит под уровень моря). Лед пластовой залежи 2 представляет собой чередование слоев чистого прозрачного стекловидного льда и льда с примесью минеральных включений (глин, торфяных остатков, песчаных частиц) в виде взвеси. С глубиной лед становится прозрачнее, включает в себя мелкие пузырьки воздуха, которые находятся под давлением. Пластовый лед 2 залегает в засоленных отложениях, датированных поздним плейстоценом. Наиболее интересные его черты проявляются в распределении изотопных характеристик с глубиной и по отношению к более молодому пластовому льду 1. Пластовый лед 1 находится на больших абсолютных высотах, датируется концом позднего плейстоцена [7] и имеет заметно более легкие значения δ18О и δ2Н (табл. 1, рис. 3).
Таблица 1. Изотопные параметры пластового льда Марре-Сале. По данным из Е.А. Слагоды и др. [5] и Г.Е. Облогова [10]
Полевой номер образца
|
Глубина, м
|
δ18О, ‰
|
δ2Н, ‰
|
dexc, ‰
|
ПЛ2 - нижний ярус [10]
|
11MS7
|
6,8
|
–10,02
|
–76,7
|
3,5
|
11MS7
|
7,6
|
–10,85
|
–83,3
|
3,5
|
11MS7
|
8,1
|
–8,82
|
–67,9
|
2,7
|
11MS7
|
8,8
|
–13,61
|
–104,4
|
4,5
|
11MS7
|
9,3
|
–13,36
|
–102,3
|
4,6
|
11MS7
|
9,7
|
–12,76
|
–98,2
|
3,9
|
11MS7
|
10,2
|
–15,56
|
–119,8
|
4,7
|
11MS7
|
14
|
–17,74
|
–137,2
|
4,7
|
12MS05
|
4,9
|
–9,03
|
–71,7
|
0,6
|
12MS05
|
5,2
|
–9,4
|
–74
|
1,3
|
12MS05
|
5,5
|
–9,32
|
–74,9
|
–0,3
|
12MS05
|
5,8
|
–9,49
|
–77
|
–1,1
|
12MS05
|
6,1
|
–14,81
|
–116,2
|
2,3
|
12MS05
|
6,4
|
–13,45
|
–105,5
|
2,1
|
12MS05
|
6,7
|
–13,64
|
–105,4
|
3,7
|
12MS05
|
7
|
–17,94
|
–137,8
|
5,7
|
13MS07
|
4,3
|
–10,1
|
–78
|
2,8
|
13MS07
|
4,6
|
–11
|
–85
|
3
|
13MS07
|
4,9
|
–12,6
|
–96
|
4,8
|
13MS07
|
4,9
|
–12,4
|
–95
|
4,2
|
13MS07
|
4,9
|
–10,7
|
–82
|
3,6
|
ПЛ1 – верхний ярус [10]
|
11MS1
|
4,4
|
–20,58
|
–155,3
|
9,3
|
11MS1
|
4,6
|
–20,45
|
–154,4
|
9,2
|
12MS01
|
4,4
|
–19,67
|
–147
|
10,3
|
12MS01
|
4,6
|
–20,18
|
–152,9
|
8,6
|
12MS01
|
4,8
|
–18,51
|
–138,6
|
9,5
|
12MS03
|
9,2
|
–20,62
|
–156
|
8,9
|
12MS03
|
9
|
–19,88
|
–150,2
|
8,9
|
12MS03
|
9,4
|
–19,7
|
–150,1
|
7,5
|
12MS03
|
9,4
|
–19,16
|
–146,1
|
7,1
|
12MS04
|
9,5
|
–19,18
|
–147
|
6,5
|
12MS04
|
?
|
–18,53
|
–141,2
|
7
|
12MS04
|
?
|
–18,6
|
–142,1
|
6,7
|
13MS06
|
8
|
–21,1
|
–161
|
7,8
|
13MS09
|
7
|
–17,3
|
–129
|
9,4
|
13MS04
|
7
|
–19,4
|
–147
|
8,2
|
13MS05
|
7
|
–19
|
–144
|
8
|
13MS02
|
7
|
–18,9
|
–143
|
8,2
|
Верхняя залежь льда [5]
|
МС-2/2-08
|
21,0
|
–21,1
|
–163,2
|
5,9
|
МС-1/31-08
|
13,0
|
–24,1
|
–185,9
|
7,2
|
МС-2/2-08
|
25,6
|
–21,0
|
–161,8
|
6,1
|
МС-2/2-08
|
25,6
|
–23,7
|
–183,6
|
5,7
|
МС-2/2-08
|
21,0
|
–15,2
|
–112,4
|
8,8
|
МС-2/2-08
|
22,0
|
–14,3
|
–108,6
|
5,5
|
Нижняя залежь льда [5]
|
МС-3-09
|
13,6
|
–9,8
|
–77,9
|
0,2
|
МС-1-09
|
15,7
|
–21,5
|
–162,8
|
9,0
|
Рис. 3. Соотношение δ2Н–δ18О для пластового льда и природных вод региона: 1 – озерная вода по [8], 2 – морская вода по [8], 3 – пластовый лед ПЛ2 по [10], 4 – пластовый лед ПЛ1 по [10], 5 – пластовый лед нижней залежи по [5], 6 – пластовый лед верхней залежи по [5].
Точки изотопных значений пластовых льдов нижнего яруса ПЛ2 и двух образцов нижней залежи из работы [5] формируют на изотопной диаграмме линейный тренд, описывающийся уравнением δ2Н = 7.55 × δ18О − 2.4 (см. рис. 3). Нижний ярус пластовых льдов характеризуется довольно высокими значениями δ2Н и δ18О по сравнению с вышележащим ярусом пластовых льдов. Существует несколько гипотез формирования пластовых льдов Ямала, рассмотрим основные из них и обсудим основные противоречия, возникающие при сопоставлении этих гипотез с полученными изотопными данными, опубликованными в работах [5, 7, 10].
Изотопные противоречия в ледниковой гипотезе
По представлениям Ф.А. Каплянской и В.Д. Тарноградского [2] нижняя толща представлена преимущественно позднеплейстоценовым ледниковым комплексом в виде сильнодислоцированной основной морены. Эти мореноподобные отложения содержат пластовые льды ледникового генезиса, эрратический каменный материал в виде мелкого гравия и гальки, а также мелко раздробленные растительные остатки (см. рис. 1). Эти факты, по мнению ряда исследователей [2; 11; 14] указывают на ледниковое происхождение четвертичных отложений, слагающих рассматриваемые участки западного побережья Ямала.
Обычно ледниковый лед арктических районов (и атмосферных осадков) характеризуется связью значений δ18О и δ2Н с угловым коэффициентом около от 7 до 8 [15]. Поскольку по пластовому льду ПЛ2 получена близкая величина наклона (7.55) в координатах δ2Н–δ18О, можно оценить, могут ли пластовые льды быть остатками плейстоценовых ледников.
На ближайших арктических островах и архипелагах проводились изотопные исследования ледникового льда по материалам кернового бурения, однако в большинстве случаев, анализировался только изотопный состав кислорода [16; 17]. Парные изотопные определения доступны по ледниковому керну с купола Академии наук на Северной Земле [18]. Линия регрессии для этого льда имеет вид δ2Н = 7.57 × δ18О + 2.13, R2 = 0.98. Точки изотопных значений ледникового льда практически лежат на линии метеорных вод, но в области более легких значений δ18О и δ2Н, чем у ПЛ2 (рис. 4). Несмотря на очевидное совпадение наклонов для ледникового льда и льда ПЛ2 самым существенным различием являются намного более высокие величины δ2Н и δ18О пластового льда. В ПЛ2 значения δ18О изменяются от –8.8 до –17.7‰, в ледниковом льду купола Академии наук от –16.8 до –24.6‰ (см. рис. 4). Ледниковый лед имеет значения дейтериевого эксцесса выше, чем ПЛ2 (dexc ледникового льда от 6.8 до 14.1 и dexc пластового льда от 0.3 до 7.2‰). На другом ближайшем арктическом архипелаге Новая Земля ледниковый лед имеет более высокие значения δ18О и δ2Н: от –15.4 до –17.6 ‰ и от –110.1 до –127.1 ‰ соответственно, однако для всех проб льда характерны высокие величины дейтериевого эксцесса, которые не опускаются ниже 10‰ [19]. Это заставляет серьезно сомневаться в возможности ледниковой природы пластового льда, так как если лед ледниковый, т.е. атмосферный, то он должен отражать изотопный сигнал атмосферных осадков, для которых характерно положение на линии метеорных вод либо несколько выше линии метеорных вод в координатах изотопной диаграммы.
Рис. 4. Соотношение δ2Н–δ18О для пластового и ледникового льда: 1 – пластовый лед ПЛ2, 2 – пластовый лед ПЛ1, 3 – ледниковый лед керна Академии наук, Северная земля, по [18], 4 – пластовый лед нижней залежи, 5 – пластовый лед верхней залежи, по [5]
Возможно есть спорное допущение в том, что мы используем изотопные параметры современных ледниковых льдов для сопоставления с пластовым льдом позднеплейстоценового возраста, т.к. пластовый лед залегает в отложениях возрастом более 40-60 тыс. лет [11]. Однако, судя по положению изотерм зимнего периода, реконструированных для периода 40–12 тыс. лет назад для территории России [20, 21], общая картина циркуляции атмосферы в позднем плейстоцене заметно не менялась. Поэтому мы предполагаем, что атмосферные осадки, хотя и были изотопно легче современных, все равно соответствовали глобальной линии метеорных вод. Об этом же говорит и величина дейтериевого эксцесса (dexc) во льду ПЖЛ позднеплейстоценового возраста [10], в которых dexc варьирует от 6 до 11‰. Поскольку ПЖЛ формируются из атмосферных осадков зимнего и весеннего сезонов, то можно полагать, что атмосферные осадки того периода соответствовали глобальной линии метеорных вод. Пластовые льды нижнего яруса, относимые многими авторами к ледниковым [22, 23], в действительности не несут в себе очевидных признаков погребения ледникового льда: «…они залегают согласно вмещающим отложениям и деформированы совместно с ними. Даже если предположить, что карский диамиктон является мореной плейстоценового ледникового щита, то лёд пластов вероятнее всего сформировался внутри грунта при замерзании подледниковых вод» [24]. Пластовые льда верхнего яруса, несмотря на практически полное совпадение значений δ18О и δ2Н с ледниковым льдом купола Академии наук, скорее отражают условия открытой системы и сегрегационного сингенетического льдообразования, чем захоронения глетчерного льда. На это указывает само строение пластового льда, которое в работе называют ледогрунтом [7], поскольку он состоит из чередующихся слоев более чистого льда и льда с большим количеством грунтовых частиц. Нижний контакт ледогрунтовой толщи подчёркивается мощным выдержанным ледяным шлиром толщиной 5 см. Самым слабым местом ледниковой гипотезы образования пластовых льдов нижнего яруса Марре-Сале по мнению авторов являются высокие значения δ18О и δ2Н льда при низких величинах дейтериевого эксцесса, значительно меньших, чем для ледникового льда.
Изотопные противоречия в морской гипотезе
Поскольку лед пластовый лед нижнего яруса (ПЛ2) залегает под засоленными глинами и изотопно более тяжелый, была предложена гипотеза формирования этого льда с участием морской воды, которая осталась в лагунах и изолированных озерах после трансгрессии моря [3; 7, 10]. Предполагается, что морская вода была опреснена, как например, современная вода Баренцева моря у побережья Ямала со значениями δ18О = –6.8 и –5.7‰ [8]. Однако, мы полагаем, что эти воды не могли участвовать в формировании льда ПЛ2, поскольку они достаточно далеко отстоят по своим изотопным параметрам ото льда (см. рис. 3) и отличаются высокими значениями дейтериевого эксцесса: 7.4 и 11.8‰ соответственно. Эти значения dexc слишком велики для морской воды, поскольку даже для опресненных морских вод с соленостью около 11 г/л величина дейтериевого эксцесса должна не превышать 4‰ [25]. Значения дейтериевого эксцесса выше 7‰, полученные в работе [8], заставляют усомниться в верности этих изотопных данных для морской воды. Если формирование льда происходило из морской воды, даже опресненной, то образование льда с минерализацией от 0.032 до 0.218 г/л [10] должно было обеспечить криогенное концентрирование солей при льдообразовании и формирование мощного соленого криопэга ниже или сбоку от пластовой залежи. Еще одним противоречием, содержащимся в работе [7; 10] является предложенная ими модель формирования льда из озера с опресненной морской водой в условиях закрытой системы. Эта модель должна означать, что промерзает обводненный подозерный талик, а слой засоленных глин, перекрывающий пластовый лед – это озерные донные осадки. В этом случае, при промерзании водной линзы в основании талика, должна быть реализована закрытая система. А полученное отношение для льда δ2Н = 7,55 × δ18О – 2,401 (см. рис. 4) указывает не на закрытую систему, а скорее на формирование льда в открытой системе. В дискуссии по вопросу об участии соленых морских вод в образовании мощных пластовых залежей активно участвовал С.М. Фотиев [26, 27], который указывал на формирование пластовых льдов из озерной пресной воды. Основным его доводом против морской гипотезы является гидрохимический состав пластовых льдов, которые помимо низкой минерализации (пресные льды) имеют другой состав преобладающих ионов – возрастает концентрация континентальных гидрокарбонат-иона и кальция [26, 27]. Здесь мы должны отметить, что состав ионов пластовых льдов ПЛ2 неизменно хлоридно-натриевый [7], что в целом отвечает морскому источнику, а в отложениях ПЛ1 он меняется снизу вверх на преобладание континентальных, что свидетельствует об участии атмосферной влаги как основного источника пластовых льдов.
Изотопные противоречия в озерной гипотезе
На озерные воды, как источник мощных пластовых льдов Ямала указывали Ю.К.Васильчук [28] и С.М. Фотиев [26, 27], хотя ими рассматривались два разных сценария поступления воды к месту образования пластового льда. В первом случае, пластовые льды формировались из воды несквозного талика под озером, водяная линза которого пополнялась атмосферными осадками (как и вода самого озера). При промерзании такого водонасыщенного талика, как пишет Ю.К.Васильчук [28]: «дальнейшее промерзание могло приводить образованию мощных пластовых ледяных залежей, перекрытых озёрными осадками» (рис. 5, А). В таком случае реализуется классическая закрытая система, что должно приводить к изотопному фракционированию, отражающемуся в значениях δ2Н и δ18О льда. В случае с пластовыми льдами Марре-Сале как нижнего, так и верхнего яруса, эта модель не реализуется.
Во втором случае, согласно доводам С.М.Фотиева [26, 27], пластовые льды формировались в результате инъекции пресной озерной воды термокарстовых котловин. Проникновение озерной воды в толщу мерзлых отложений могло быть возможным лишь тогда, когда абсолютные отметки дна промерзающего озера опускались ниже абсолютных отметок зоны контакта пачки глин и пачки песка в бортах термокарстовой котловины. Согласно его гипотезе, формирование на озере сплошного ледового покрова приводило к образованию в озерной котловине замкнутой системы. Промерзающая в замкнутой котловине озерная вода приобретала огромное криогенное давление, под воздействием которого через зону контакта по криогенному водоупору она проникала в толщу мерзлых отложений. Водоупором, по его мнению, служила кровля пачки мерзлых песков. «Если толщина озерного льда была меньше глубины озера, то из замкнутой озерной котловины под мощным криогенным давлением к подошве пласта льда поступала прозрачная озерная вода (рис. 5, Б), а если больше, то поступала мутная вода или разжиженные осадки из водоносных таберальных отложений подозерного талика. Наличие в отдельных слоях пласта песка, гравия, гальки и водорослей, а также слоев ледогрунта свидетельствует о том, что вода, поступающая из озера или подозерного талика, не фильтровалась через песчаные слои. Поток воды по кровле криогенного водоупора под напором проникал вдоль подошвы пласта на различные (иногда значительные) расстояния. Только при таком допущении песок, гравий, галька и водоросли могли оказаться внутри отдельных слоев льда» [27]. Нам представляется маловероятным, чтобы незамерзшая вода могла раздвинуть мерзлые отложения и внедриться между мерзлыми плотными глинами и подстилающими ее мерзлыми засоленными песками. И сделать это на большие по площади расстояния. Если вода внедрялась в пески, то она должна была растворить содержащиеся в нем соли и стать соленой, т.к. согласно его схеме (см. рис. 5, Б) в песках были криопэги, но согласно его же доводам, вода, формирующая пластовые льды, была пресная.
Довольно странно, что криогенным водоупором являлись мерзлые пески. Как правило, мерзлые пески имеют малую льдистость, и меньшую плотность по сравнению с перекрывающими их глинами. Глины в данном случае, многие авторы считают важным фактором формирования льдов, поскольку именно на их границе происходит миграция поровой, капиллярной и связанной воды, что важно для сегрегационного льдообразования [29].
Рис. 5. Формирование пластовых льдов на дне несквозного подозерного талика (А) при атмосферном питании озера по Ю.К.Васильчуку [28, с. 62], и (Б) при инъекции озерной воды в стенку котловины на литологическую границу глины/пески по С.М.Фотиеву [26]
Однако в случае озерной гипотезы был предложен инъекционный механизм поступления воды к фронту промерзания. Если рассматривать предложенную С.М.Фотиевым модель, то при льдообразовании должна быть закрытая система. Но она в случае пластовых льдов нижнего яруса Марре-Сале (как, впрочем, и верхнего яруса пластовых льдов) не реализуется. Это является основным изотопным противоречием в озерной гипотезе, если применять ее к пластовым льдам нижнего яруса Марре-Сале.
Изотопные противоречия в гипотезе льдообразования в закрытой систем
Интригующей особенностью ПЛ2 является наклон линии регрессии, описывающей значения δ18О и δ2Н льда. Как видно на рис. 3, все установленные значения по ПЛ2 аппроксимируются линий с наклоном 7.5 (R2 = 0.99). В работе [7] для ПЛ 2 приводится уравнение регрессии вида δ2Н = 4,9 × δ18О − 64,2. Очевидно, в этом опубликованном уравнении содержится ошибка, либо данные приведены не полностью, поскольку в опубликованной работе [10] содержится фактический материал (см. табл. 1), линейный тренд по которому имеет вид δ2Н = 7.55 × δ18О − 2,4 (рис. 6, а). Установленная величина наклона не соответствует замерзанию свободной воды в условиях закрытой системы.
Рис. 6. Расчетные величины δ18О, δ2Н и dexc льда в закрытой системе и несоответствие им реальных значений δ18О, δ2Н и dexc льда ПЛ2: 1 – изотопные характеристики ПЛ2 (из [10]); 2 – расчет модели закрытой системы; 3 – изотопные характеристики нижней залежи льдов (из [5])
Расчет модели закрытой системы показал, что динамика значений изотопного состава кислорода и водорода льда должна определяться соотношением δ2Н = 6.62 × δ18О − 12,2 при использовании коэффициентов, опубликованных в [30] и в рэлеевском процессе при значениях δ18О и δ2Н начальной воды –11.6 и –89‰. Значения начальной воды были рассчитаны по максимальной величине значений δ18О и δ2Н пластового льда (–8.8‰ и –67.9‰, см. табл.1) исходя из δ18Ольда = δ18Оводы + ɛ, аналогично для δ2Н, величина ɛ18О=2.9‰, ɛ2Н =21‰ согласно коэффициентам фракционирования в системе вода-лед [30]. Не только несоответствие расчетного наклона реальным образцам, но и установленный диапазон изменения значений δ18О и δ2Н льда пластовой залежи не соответствует расчетному. Выполненный расчет замерзания воды в закрытой системе имеет диапазон изменения δ18О в 8.5‰ от фракции льда f=0.01 (1% льда по отношению к воде) до 0.95 (95% льда по отношению к замерзающей воде). А во льду ПЛ2 общий разброс значений δ18О составил 9.1‰. Это же касается пластового льда нижней залежи, исследованного Е.А.Слагодой и др. [5] (см. рис. 6, а), где разброс значений δ18О в образцах льда, удаленных друг от друга на 2.1 м по высоте составил 11.7‰ (см. табл. 1). При этом точки изотопных значений этого льда также находятся на линии регрессии с наклоном 7.55. Также на несоответствие изотопных характеристик льда закрытой системе указывает отношение δ2Н-dexc(рис. 6, б), которое используется для подобной диагностики [15]. Возможно, этот лед формировался в условиях открытой системы, а опробованные фрагменты представляют разные генерации льда и разные фракции льда, сформированного из атмосферной влаги, меняющей свои изотопные параметры во времени. Поэтому изотопный состав кислорода и водорода льда ПЛ2 имеет большой диапазон. Следовательно, необходимо понять почему изотопный состав воды, питающей пластовый лед, менялся (поскольку значения δ18О льда сверху вниз по пластовой залежи меняются от –8.8 до –17.7‰, следовательно менялся изотопный состав воды, из которой он сформирован) и каким образом была образована залежь.
Дискуссия
Источник воды. Если принять, что пластовые льды нижнего яруса формировались в условиях открытой системы, а на это указывает наклон линии регрессии равный 7.55, то поскольку коэффициенты изотопного фракционирования вода-лед всегда одинаковые (максимальные α18О = 1.0029 и α2Н = 1.021, а минимальные равны нулю при очень высоких скоростях замерзания воды), то изменение величин δ18О и δ2Н льда связано с изменением изотопных характеристик воды, из которой он формировался. При этом сами значения δ18О и δ2Н льда очень высокие, во всяком случае, намного выше, чем в пластовом льду верхнего яруса. Таким образом, нам следует вернуться к рассмотрению возможности участия опресненных морских вод в качестве источника. Последние исследования опреснения морской воды пресным компонентом в Карском море показали, что опресненная морская вода со значениями δ18О = –10.5 и –7.8‰ имеет значения dexc = 3.4 и 1.8‰ соответственно [25]. Точки этой опресненной морской воды на изотопной диаграмме (рис. 7) лежат на линии, соединяющей неопресненные воды Баренцева моря (δ18О=+0,26‰ и δ2Н=+1,55‰) и пластовые льды верхнего яруса, которые в свою очередь нами принимаются как сформированные из атмосферной влаги (см. рис. 7). Эта линия смешения полностью описывает распределение δ18О и δ2Н пластового льда ПЛ2, точнее точки пластового льда соответствуют линейному тренду смешения морской и пресной воды. Атмосферные осадки на Ямале и Гыданском п-ове в позднем плейстоцене по данным [21] характеризовались изотопным составом кислорода от –21 до –23‰, что также вполне соответствует линии смешения на рис. 7.
Это дает нам основание считать, что опресненные морские воды участвовали в формировании пластового льда ПЛ2, при этом при продолжающемся процессе формирования пластового льда эти воды все более и более опреснялись пресным атмосферным компонентом – атмосферными осадками.
Рис. 7. Соотношение δ2Н-δ18О для морской воды и пластового льда: 1 – пластовый лед ПЛ2 (из [10]), 2 – опресненные воды Карского моря (по [25]), 3 – пластовый лед ПЛ1 (из [10]), 4 – неопресненные воды Баренцева моря (по [25])
Опресненные морские воды на ранних этапах формирования пластового льда, судя по положению точек опресненных вод Карского моря на изотопной диаграмме (см. рис. 7) должны были иметь соленость 11.7 и 17 промилле, т.е. 11-17 грамм на килограмм морской воды (изотопный состав кислорода и водорода этих вод, а также соответствующая им величина солености взята из работы Дубининой и др. [25]). А пластовые льда ПЛ2 имеют общую минерализацию до 0.2 г/л. Это означает, что морские соли были практически полностью отжаты при льдообразовании, что возможно при медленном образовании льда из морской воды. Это сопровождается концентрированием солей в остающейся воде. Именно из такой воды образуются криопэги.
Особенности формирования пластового льда нижнего яруса. Такое криогенное отжатие солей, по нашему мнению, возможно при одном условии: если вода поступала к литологической границе глины/пески по песчаному водоносному горизонту. Существующее озеро с опресненной морской водой имело глубокий талик, подстилавшийся песчаным водоносным горизонтом (рис. 8). Образование озерного льда зимой создавало давление в талике и приводило к увеличению водности горизонта в зимний период. Свободная вода (горизонт был довольно обширный) поступала к месту литологической границы (точка ПЛ2) и здесь на верхней границе формировался конжеляционный лед. Соленость начальной воды в горизонте несколько повышалось из-за отжатия солей при образовании льда. Поскольку воды в горизонте было очень много, и она была засоленной (не пресной) лед образовывался до тех пор, пока концентрация солей не становилась выше температурного предела замерзания. Льда образовывалось значительно меньше, чем оставалось воды в горизонте. Поэтому закрытая система не реализовывалась.
Рис. 8. Схема возможного образования пластового льда нижнего яруса Марре-Сале: 1 – песчаный водоносный горизонт; 2 – глины; 3 – граница мерзлых пород; 4 – пластовый лед; 5 – потоки свободной воды; 6 – криогенное отжатие солей; 7 – криопэги
В последующие сезоны происходило все большее опреснение озерной воды, она же поступала в водоносный песчаный горизонт. Повторялась та же картина. Поскольку лед образовывался на верхней границе, то все отложения выпучивались на высоту образованного льда. Вода, поступающая из озера в горизонт, становилась все более пресной. Но в водоносном горизонте должен был существовать слой соленой воды, отжатой при льдообразовании. Возможно ли предположить, что внутри очень водного горизонта существовали условия меромиксии, когда поступающая все более пресная вода находилась в верхней части, а более соленая и плотная (тяжелая) в нижней части? Тогда это объясняет, что лед (который образовывался на верхней границе горизонта) становился все более пресным (изотопно легким), а в нижней части существовал засоленный слой, который, собственно, и являлся одним большим криопэгом. Сочетание на Ямале в одном разрезе пластовых льдов и криопэгов, образовавшихся в результате криогенного концентрирования солей было рассмотрено в работе [31] и показана их генетическая связь. Эта же связь прослеживается и при нашей интерпретации изотопного состава кислорода и водорода льда (см. рис. 8).
Предложенная нами схема образования пластовых льдов нижнего яруса имеет общие черты с озерной гипотезой, уточняя ее с позиции изотопной идентификации. Эта модель описывает установленные значения δ18О и δ2Н льда Марре-Сале, однако, широкое многообразие условий образования пластовых льдов Ямала предполагает, что такой механизм является далеко не единственным. Очевидно, что многие пластовые залежи могли формироваться при промерзании несквозных таликов в классической закрытой системе. Если в талике был большой объем воды в линзе, то образованные из нее пластовые залежи обнаруживают сочетания разного механизма образования льда – например, при промерзании сверху такой линзы формируется слой льда, оставшаяся вода под большим давлением внедряется в него, формируя инъекционные штоки [5, 32]. Также распространены на Ямале пластовые льды, сформированные в условиях открытой системы при сингенетическом промерзании обводненных пространств с одновременным формированием повторно-жильных льдов. Все образованные в разных условиях типы льда несут свою изотопную метку, зависящую не только от литологических условий и разного источника воды, но и от разного количества воды, пошедшей на формирование льда. Еще одной интересной задачей является оценка того, какая вода образовывала лед – в случае сегрегационного льда это должна быть капиллярная и пленочная вода, мигрирующая к фронту промерзания. Рассмотренные нами в данной работе льды Марре-Сале, вероятнее всего, формировались из свободной воды, присутствовавшей в большом объеме, и лед этот, скорее всего, конжеляционный.
Выводы
Рассмотрение опубликованных изотопных данных по пластовым льдам Ямала в районе полярной станции Марре-Сале показало, что индикация генезиса пластовых льдов по изотопному составу кислорода и водорода льда имеет ряд трудностей. Одной из таких сложностей при интерпретации изотопных вариаций пластового льда нижнего яруса неожиданно стали ошибочные данные величин δ18О и δ2Н опресненной морской воды вблизи побережий Ямала [8]. Использование этих данных при интерпретации величин δ18О и δ2Н льда не позволило уверенно диагностировать участие этой морской воды в образовании льда. Сопоставление с неопресненными морскими водами позволило установить процесс смешения с изотопно легкими атмосферными осадками и показало, что прибрежная морская вода участвовала в образовании льда.
Другой трудностью стала идея об образовании пластового льда нижнего яруса Марре-Сале в закрытой системе. Моделирование распределения значений δ18О и δ2Н при льдообразовании в условиях закрытой системе показало, что наклон линии регрессии для такого льда должен составлять 6.3. Установленная величина наклона (7.55) линии регрессии описывающего соотношение значений δ2Н и δ18О пластового льда нижнего яруса Марре-Сале не соответствует замерзанию воды в условиях закрытой системы. Скорее, образование пластового льда происходило в условиях открытой системы с большим количеством свободной воды, из которой образовывался лед. Наклон ниже 8 в координатах изотопной диаграммы связан не с закрытой системой, а с участием опресненной морской воды на первых стадиях образования залежи.
С позиции ледниковой гипотезы образования пластовых льдов нижнего яруса Марре-Сале труднообъяснимыми являются высокие значения δ18О и δ2Н льда при низких величинах дейтериевого эксцесса, значительно меньших, чем для ледникового льда.
Рассмотренные нами в данной работе пластовые льды Марре-Сале, вероятнее всего, формировались из свободной воды, в большом объеме поступавшей по песчаному водоносному горизонту из озерного талика. Лед образовывался на границе глины/пески при промерзании сверху. Вода в озере на первой стадии являлась смесью морских и атмосферных вод, а затем все больше опреснялась атмосферными осадками. С этим процессом опреснения и изотопного облегчения озерной воды связано уменьшение значений δ18О и δ2Н льда нижнего яруса с глубиной. Слабыми местами такой версии образования льда, является сочетание в одном водоносном горизонте большого количества пресной воды и воды криопэгов. Вторым слабым местом является отсутствие сведение о протяженности и конфигурации песчаного горизонта, являвшемся водоносным, как впрочем и в целом параметров пластового льда нижнего яруса, включая его протяженность по латерали и по глубине. Третьим слабым местом является редкое и недостаточно детальное опробование, и как следствие, небольшое количество изотопных определений.
Недостаток изотопных определений осложняет и генетическую интерпретацию верхней ледяной залежи в Марре-Сале.
Благодарности
Авторы благодарят коллег, опубликовавших данные изотопных определений из пластовывх льдов Марре-Сале, послуживших основным материалом для размышлений над изотопными противоречиями.
References
1. Kuzin I.L., Astaf'ev N.F. Kriogennye dislokatsii na zapadnom poberezh'e poluostrova Yamal // Izv. VGO. 1975. Tom. 107. Vyp. 6. S. 510–515.
2. Kaplyanskaya F.A., Tarnogradskii V.D. Reliktovye gletchernye l'dy na severe Zapadnoi Sibiri i ikh rol' v stroenii raionov pleistotsenovogo oledeneniya kriolitozony // Doklady AN SSSR. 1976. Tom 231. №5. S. 1185–1187.
3. Danilov I.D. Dislokatsii v merzlykh, soderzhashchikh plastovye l'dy pleistotsenovykh otlozheniyakh severa Zapadnoi Sibiri // Formirovanie merzlykh porod i prognoz kriogennykh protsessov. M.: Nauka. 1986. S. 28–41.
4. Vasil'ev A.A., Rogov V.V. Plastovye l'dy v raione Marre-Sale, Zapadnyi Yamal // Materialy Vtoroi konferentsii geokriologov Rossii, tom 1, chast' 2. Litogeneticheskaya geokriologiya. Inzhenernaya geokriologiya. MGU im. M.V. Lomonosova 6-8 iyunya 2001 g. M.: Izd-vo Mosk. un-ta. 2001. S. 188–194.
5. Slagoda E.A., Opokina O.L., Rogov V.V., Kurchatova A.N. Stroenie i genezis podzemnykh l'dov v verkhneneopleistotsen-golotsenovykh otlozheniyakh mysa Marre-Sale (Zapadnyi Yamal) // Kriosfera Zemli. 2012. Tom XVI. №2. S. 9–22.
6. Streletskaya I.D., Kanevskii M.Z., Vasil'ev A.A. Plastovye l'dy v dislotsirovannykh chetvertichnykh otlozheniyakh Zapadnogo Yamala // Kriosfera Zemli. 2006. Tom X. №2. S. 68–78.
7. Streletskaya I.D., Vasil'ev A.A., Oblogov G.E., Matyukhin A.G. Izotopnyi sostav podzemnykh l'dov Zapadnogo Yamala (Marre-Sale) // Led i Sneg. 2013. №2. S. 83–92.
8. Kritsuk L.N. Podzemnye l'dy Zapadnoi Sibiri. M.: Nauchnyi mir. 2010. – 351 s.
9. Opokina O.L., Slagoda E.A., Kurchatova A.N. Stratigrafiya razreza «Marre-Sale» (Zapadnyi Yamal) s uchetom novykh dannykh radiouglerodnogo analiza // Led i Sneg. 2015. T.55. №4. S.87–94.
10. Oblogov G.E. Evolyutsiya kriolitozony poberezh'ya i shel'fa Karskogo morya v pozdnem neopleistotsene – golotsene. Dissertatsiya na soiskanie uchenoi stepeni kandidata geologo-mineralogicheskikh nauk. Tyumen'. 2016. – 197 s.
11. Forman S.L., Ingolfsson O., Gataullin V., Manley W., Lokrantz H. Late Quaternary stratigraphy, glacial limits, and paleoenvironments of the Marresale Area, western Yamal Peninsula, Russia // Quaternary Research. 2002. Vol.57. №3 P. 355–370.
12. Kanevskii M.Z., Streletskaya I.D., Vasil'ev A.A. Zakonomernosti formirovaniya kriogennogo stroeniya chetvertichnykh otlozhenii Zapadnogo Yamala (na primere raiona Marre-Sale) // Kriosfera Zemli. 2005yu Tom IX. № 3. S. 16–27.
13. Forman, S., Ingólfsson, Ó, Gataullin, V., Manley, W., & Lokrantz, H. Late Quaternary Stratigraphy, Glacial Limits, and Paleoenvironments of the Marresale Area, Western Yamal Peninsula, Russia // Quaternary Research. 2002. 57(3). P. 355–370.
14. Solomatin V.I., Belova N.G. Dokazatel'stva pogrebennogo gletchernogo proiskhozhdeniya plastovykh l'dov // Trudy desyatoi mezhdunarodnoi konferentsii po merzlotovedeniyu TICOP. Salekhard. 2012. S.493–497.
15. Lacelle D. On the δ18O, δD and D-excess relation in meteoric precipitation and during equilibrium freezing: theoretical approach and field examples // Permafrost and Periglacial Processes. 2011. № 22. P. 13–25.
16. Kotlyakov V.M., Arkhipov S.M., Henderson K.A., Nagornov O.V. Deep drilling of glaciers in Eurasian Arctic as a source of paleoclimatic records // Quaternary Science Reviews. 2004. Vol. 23. P. 1371–1390.
17. Arkhipov S.M., Kotlyakov V.M., Punning J-M.K., Zogorodnov V.S., Nikolayev V. I., Macheret Yu.Ya., Vaikmaye R., Barkov N.I., Korsun S.A., Korotkevich V., Morev V.A., Evseyev A.V., Vostokova T.A., Andreev A.A., Klementyev O.L., Korotkevitch Y. S., Stiévenard M., Sinkevich S.A., Samoylov O.Yu., Gordienko F. G., Korsun A.V., Tiugu K. R. Deep drilling of glaciers: Russian projects in the Arctic (1975-1995). Institute of Geography, Russian Academy of Sciences, 2008. Moscow, PANGAEA, https://doi.org/10.1594/PANGAEA.707363
18. Opel T., Fritzsche D., Meyer H., Schuëtt R., Weiler K., Ruth U., Wilhelms F., Fisher H. 115 year ice-core data from Akademii Nauk ice cap, Severnaya Zemlya: high-resolution record of Eurasian Arctic climate change // Journal of Glaciology. 2009. Vol. 55(189). P. 21–31.
19. Chizhova Yu.N., Dubinina E.O., Kossova S.A., Avdeenko A.S., Miroshnikov A.Yu. Formirovanie izotopnykh parametrov lednikov arkhipelaga Novaya Zemlya. Materialy XXII Simpoziuma po geokhimii izotopov im. A.P. Vinogradova (29-31 oktyabrya 2019), Rasshirennye tezisy dokladov, M: GEOKhI RAN, 2019. S.482-489
20. Vasil’chuk Yu., Vasil’chuk A. Spatial distribution of mean winter air temperatures in Siberian permafrost at 20–18 ka BP using oxygen isotope data // Boreas. 2014. Vol. 43. Iss. 3. P. 678–687. doi: 10.1111/bor.12033
21. Vasil'chuk Yu.K. Prostranstvenno–vremennoe raspredelenie sredneyanvarskikh paleotemperatur vozdukha na territorii Rossiiskoi Arktiki dlya perioda 30–12 tys. let nazad s vysokim vremennym razresheniem // Arktika i Antarktika. — 2016. – № 1. – S. 86–103. DOI: 10.7256/2453-8922.2016.1.21310 (Vasil'chuk, Y.K. 2017. Spatio-temporal distribution of mean January air temperature over the Russian Arctic during 30-12 ka BP with high temporal resolution. Arctic and Antarctic. No. 1. P. 86–103. DOI: 10.7256/2453-8922.2016.1.21310)
22. Tarnogradskii V.D. O proiskhozhdenii plastovykh zalezhei podzemnykh l'dov na Karskom poberezh'e p‑ova Yamal // Plastovye l'dy kriolitozony. Yakutsk: izd. IMZ SO AN SSSR, 1982. S. 80–89.
23. Gataullin V.N. Plastovye l'dy zapadnogo poberezh'ya p-va Yamal: ikh stroenie, sostav i proiskhozhdenie // Geokriologich. issledovaniya v Arkticheskikh raionakh: Mezhdunar. simpozium. Yamburg, avgust 1989 g. Vyp. 1.Tyumen': izd. IPOS SO ANSSSR, 1990. S. 3–11.
24. Belova N.G. Plastovye l'dy yugo-zapadnogo poberezh'ya Karskogo morya. M.: MAKS Press, 2014. – 180 s.
25. Dubinina E.O., Kossova S.A., Miroshnikov A.Yu., Fyaizullina R.V. Izotopnye (D, O) parametry i istochniki opresnennykh vod Karskogo morya // Okeanologiya. 2017. Tom 57. № 1. S. 38–48.
26. Fotiev S.M. Mekhanizm formirovaniya in''ektsionnykh zhil l'da i gidrolakkolitov // Kriosfera Zemli. 2011. Tom 15. №2. S. 44–55.
27. Fotiev S.M. Genezis i mekhanizm formirovaniya plastov povtorno-in''ektsionnogo l'da // Kriosfera Zemli. 2015. Tom 19. №1. S. 30-40.
28. Vasil'chuk Yu.K. Izotopnye metody v geografii. Chast' 2: Geokhimiya stabil'nykh izotopov plastovykh l'dov. V 2-kh tomakh. M.: Izdatel'stvo Moskovskogo universiteta, 2012 Tom I – 472 s. (Vasil’chuk Y.K. 2012. Isotope Ratios in the Environment. Part 2: Stable isotope geochemistry of massive ice. Moscow: Moscow University Press. 472 p.).
29. Vtyurina E.A., Vtyurin B.I. L'doobrazovanie v gornykh porodakh. M., Nauka, 1970. 279 s.
30. Lehmann M., Siegenthaler U. Equilibrium oxygen– and hydrogen–isotope fractionation between ice and water // Journal of Glaciology. 1991. Vol. 37. N.125. P. 23–26.
31. Streletskaya I.D., Leibman M.O. Kriogeokhimicheskaya vzaimosvyaz' plastovykh l'dov, kriopegov i vmeshchayushchikh ikh otlozhenii Tsentral'nogo Yamala // Kriosfera Zemli. 2002. Tom 6. №3. S. 15–24.
32. Vasil'chuk Yu.K. Gomogennye i geterogennye plastovye ledyanye zalezhi v mnogoletnemerzlykh porodakh // Kriosfera Zemli. 2011. Tom XV. № 1. S. 3–17.
|