CC BY
2
Научная статья
Происхождение карбонатно-силикатных пород Порьей губы (Лапландско-Колвицкий гранулитовый пояс) по результатам анализа стабильных изотопов (8180, 813С)
Д.П.КРЫЛОВ Н, Е.В.КЛИМОВА
Институт геологии и геохронологии докембрия РАН, Санкт-Петербург, Россия
Как цитировать эту статью: Крылов Д.П., Климова Е.В. Происхождение карбонатно-силикатных пород Порьей губы (Лапландско-Колвицкий гранулитовый пояс) по результатам анализа стабильных изотопов (5180, 513С) // Записки Горного института. 2024. Т. 265. С. 3-15. EDN GISHQG
Аннотация. В составе гранулитов Порьей губы Лапландско-Колвицкого пояса Фенноскандинавского щита наблюдаются карбонатно-силикатные породы неясного происхождения. Цель настоящей работы - реконструкция возможных источников вещества и условий метаморфических преобразований этих пород на основе изучения изотопного состава кислорода и углерода в сочетании с моделированием условий фазовых равновесий. Результаты изотопного анализа и литохимических реконструкций показывают, что карбонатно-силикатные породы Порьей губы сформировались в результате метаморфизма осадков (возможно, мергелей), изотопный состав которых соответствовал докембрийским диагенетически преобразованным карбонатам (5180 = 17,9 %о, SMOW и 513С = -3,4 %о, PDB). Химический состав варьировал в зависимости от соотношения карбонатного, глинистого и обломочного компонентов. Существенные изменения изотопного состава при метаморфизме обусловлены реакциями разложения первичных карбонатов (доломита, сидерита, анкерита) с образованием СО2 и последующей дегазацией. При этом 5180 и 513С кальцита в изотопном равновесии с СО2 понижаются до 15 %о (SMOW) и -6 %о (PDB) соответственно. Контроль изотопного состава ограничивается локальными реакциями в пределах отдельных разновидностей пород, что исключает значимое влияние магматических и (или) метасоматических процессов.
Ключевые слова: карбонатно-силикатные породы; стабильные изотопы; метаморфизм; неорганический углеродный цикл; изотопное фракционирование; дистилляция Рэлея
Благодарность. Исследование выполнено за счет гранта Российского научного фонда № 22-27-00275, https://rscf.ru/project/22-27-00275/.
Поступила: 07.02.2023 Принята: 20.06.2023 Онлайн: 06.09.2023 Опубликована: 29.02.2024
Введение. Карбонатно-силикатный геохимический цикл (неорганический углеродный цикл [1]) включает процессы выветривания, осадкообразования (с преобразованием силикатов в карбонаты), метаморфизм и вулканическую деятельность (с преобразованием карбонатов в силикаты). В интервале нескольких миллионов лет или больше данный цикл контролирует баланс СО2 при взаимодействии атмосферы, гидросферы, коры и мантии Земли и во многом определяет глобальные изменения климата [2-4]. Карбонатно-силикатные породы в метаморфических комплексах (далее для краткости используется устаревшее наименование - кальцифиры) представляют собой последовательность фазовых трансформаций на этапах цикла, связанную с образованием карбонатного субстрата, метаморфизмом, дегазацией и выделением значительных объемов СО2 в атмосферу [5-7]. Фазовые реакции в карбонатно-силикатном цикле сопровождаются изотопными реакциями, и отношения стабильных изотопов (18О/16О, 13С/12С), а также «нетрадиционных» изотопов магния, кальция, кремния и титана могут служить важнейшими геохимическими трассерами [8-10] и давать возможность оценивать влияние различных резервуаров на состав газов, температуры процессов, баланс СО2/Н2О в атмосфере, гидросфере и газах глубинного происхождения. Тем не менее изучение состава стабильных изотопов в основном ограничивается карбонатными осадочными породами низких степеней метаморфических преобразований, источники вещества
Статья опубликована в открытом доступе по лицензии ^ BY 4.0
и условия преобразования кальцифиров в пределах глубокометаморфизованных (гранулитовых) комплексов изучены меньше [11, 12].
Древнейшие гранулитовые комплексы, широко развитые на щитах и в кристаллическом основании платформ, представляют важные свидетельства о составе и степени преобразования континентальной коры, сформированной на ранних этапах развития Земли [13-15]. Лапландско-Кол-вицкий гранулитовый пояс является элементом строения Фенноскандинавского щита, который на протяжении последних лет является научным полигоном для проверки различных геологических моделей и концепций [16-18]. Палеопротерозойские события, приведшие к образованию комплексов пород Лапландско-Колвицкого пояса, также широко проявились в смежных тектонических структурах, например, в породах Беломорского подвижного пояса [19-21].
В пределах Лапландско-Колвицкого гранулитового пояса в числе прочих образований наблюдаются редкие выходы кальцифиров, происхождение которых до сих пор не установлено. Общее описание, отдельные измерения химического состава [22-24] и изотопных отношений углерода [25] позволили предположить осадочное происхождение и возможную биогенную составляющую некоторых пород пояса, однако чередование кальцифиров с предположительно метаинтрузив-ными породами ставит под сомнение правильность их отнесения к метаосадочным. Цель данной работы - изучение изотопного состава кислорода и углерода кальцифиров Порьей губы в сочетании с моделированием условий фазовых равновесий, реконструкция возможных источников вещества и условий метаморфических преобразований. Результаты изотопных измерений дополняются новыми данными по химическому и фазовому составам пород.
Методы. Определение изотопного состава. Для выделения CO2 из карбонатов (100 % H3PO4) и O2 из силикатов (фторирование с последующей конверсией в CO2) использовались стандартные процедуры [26, 27]. Полнота выделения (не менее 98 %) контролировалась измерением выхода газа, полнота конверсии - измерением остаточного давления кислорода. Изотопные отношения CO2 измерены на масс-спектрометре SIRA-9 (VG-Isogas). Весь цикл выделения и измерения для каждого образца выполнялся как минимум дважды. Кроме того, через каждые 4-5 опытов определялся изотопный состав кислорода внутрилабораторного стандарта CO2. Для контроля результатов измерения карбонатов использованы образцы кальцита, откалиброванные относительно стандарта NBS-18. Принятый кислотный фактор фракционирования кальцита (T = 25 °C) равен 1,01049 [28]. Для контроля результатов измерения силикатов использованы образцы кварца и калиевого полевого шпата, откалиброванные относительно NBS-28. Изотопные данные приведены согласно общепринятым 5-обозначениям относительно стандартов V-SMOW (для кислорода) и PDB (для углерода): 5 = [Кобр^ст - 1]1000, где Лобр/Лст - отношение содержания тяжелого и легкого изотопов в образце и стандарте. Средняя погрешность значений 518О составляет 0,07 %о, 513C - 0,05 %о (1о). Подготовка образцов и измерение изотопного состава кислорода выполнены в Институте петрологии и минералогии Бонна.
Фракционирование изотопов между фазами A и B определяется величиной aAB = RA/RB = = (5a + 1000)/(5B + 1000). Значение aAB можно определить на основе так называемых Р-факто-ров, так что при изотопном равновесии 103ln«AB(T) = 1000lnPA(T) - 1000lnPB(T). Температуры изотопного равновесия определяются по измеренным изотопным отношениям фаз при известных температурных зависимостях Р-факторов. В настоящей работе использованы зависимости P(T) в виде кубических полиномов, определенных на основе изотопных сдвигов частот колебаний фаз [29].
Расчет фазовых равновесий. Фазовые равновесия моделировались для системы SiO2-TiO2--Al2O3-FeO-MgO-CaO-Na2O-K2O-H2O-CO2 (STAFMCNK-H2O-CO2) с использованием комплекса программ Perple_X (версия 6.9, http://www.perplex.ethz.ch/) [30, 31] и GeoPS [32]. Для флюида (H2O-CO2) использовано уравнение состояния и термодинамическая база данных ([33], с изменениями 2004 г.). При вычислениях учитывались присутствующие в породах фазы, их возможные предшественники и продукты преобразования. В качестве исходных данных для расчета фазовых равновесий (псевдосекций) использованы результаты рентгенофлуоресцентной спектрометрии (прибор Philips PW-1480) с аналитической погрешностью не более 3 %. Состав отдельных фаз
определен на приборе JEOL JXA-8230 (ИГГД РАН) при ускоряющем напряжении 15 кВ, токе электронного зонда 20 нА и диаметре пучка 10 мкм. Содержание С02 в кальците определялось после высушивания образца (1 ч при T = 105 °С) по величине выхода газа при взаимодействии с 10 %-ным раствором соляной кислоты. Окончание процесса разрушения карбоната фиксировалось по прекращению выделения газа.
Геологическое положение и состав пород. Метакарбонатные породы выделены в составе нижней толщи Колвицкой зоны гранулитов (Плоскотундровская толща [22]), а также отмечены на юго-восточном окончании Лапландского блока (Сальные и Туадаш-тундры). Наибольшее количество этих пород наблюдается в районе Порьей губы Кандалакшского залива Белого моря на о-ве Медвежий и западном побережье Порьей губы, выходы кальцифиров прослеживаются вдоль простирания (СЗ-ЮВ) на 10-15 км. В работе приводятся данные по породам о-ва Медвежий и о-ва Тамаркина луда (запад Порьей губы, рис.1). В пределах обнажения на о-ве Медвежий кальцифиры образуют прослои (по простиранию прослеживающиеся до нескольких десятков метров) с мощностью от долей сантиметра до 50-60 см, чередующиеся с кристаллическими сланцами (двупироксеновыми, гранат-двупироксеновыми, амфибол-диоп-сидовыми, актинолитовыми, скаполитсодержащими и т.д.), кислыми гранулитами и в различной степени мигматизированными гнейсами. Участками сланцы и гнейсы содержат сульфидную минерализацию. Для пород толщи переслаивающихся кальцифиров и кристаллических сланцев характерны структуры будинажа и разлинзования, обусловленные различной резистентностью карбонатных и силикатных пород.
Оценки возраста метаморфизма пород гранулитового комплекса колеблются в пределах 1925 млн лет (начало метаморфизма) до 1870 млн лет (его завершение) [34]. Образование прото-лита относится к палеопротерозою, хотя имеются свидетельства и его допротерозойского (неоархейского?) образования [35].
Кальцифиры представлены средне-крупнозернистыми серыми полосчатыми породами обычно гранобластовой (гетерогранобластовой) структуры. По минеральному составу в пределах слоев кальцифиров можно выделить две основные группы:
• Кальцифиры, в которых преобладают кальцит, моноклинный пироксен (диопсид), гранат (до 13 об.%), кварц (0,1-0,5 мм, 5-15 %), иногда плагиоклаз (0,01-1,5 мм). Кальцит встречается в виде крупных полисинтетических зерен до 2-3 мм. Крупные выделения диопсида (0,3-2 мм, до 30 % объема породы) бесцветные, замещаются тремолитом и флогопитом. Пойкилобласты граната (до 1 см) преимущественно гроссуляр-андрадитового ряда с небольшим содержанием пиропа и альмандина. Однако в породах с существенным содержанием плагиоклаза (до 10-15 об.%) в составе граната (иногда реликтового, до 5-7 об.%) преобладают альмандин с пиропом. В качестве второстепенных (< 5 %) присутствуют магнезит, доломит, иногда оливин (форстерит), титанит или ильменит, а также апатит. Для краткости данная разновидность далее по тексту будет обозначаться как диопсидовые кальцифиры (^/-кальцифиры).
• Кальцифиры с преобладанием форстерита, кальцита с небольшим содержанием моноклинного пироксена. В некоторых породах встречаются идиобластовые (до нескольких миллиметров) выделения кальцита. Зерна оливина (0,3-0,5 мм) замещаются серпентином (0,2-0,3 мм). Наблюдаются также другие вторичные гидроксилсодержащие минералы - в основном слюды, хлорит, эпи-дот, амфибол (актинолит-тремолит, до 0,5 мм). Несмотря на значительное содержание MgO в оли-винсодержащей разновидности изученных кальцифиров, доломит встречается обычно в виде реликтов. В качестве акцессорных минералов наблюдаются апатит, титанит, ильменит, шпинель (герцинит). Данная разновидность карбонатно-силикатных пород далее обозначается как оливи-новые кальцифиры (^/-кальцифиры).
Помимо минерального состава отличия двух разновидностей кальцифиров проявляются и по содержанию главных окислов (табл.1). Оливиновые кальцифиры имеют существенно более высокое отношение MgO/(MgO + СаО) (около 0,4) по сравнению с кальцифирами первой группы (менее 0,1). Содержание FeO и МпО в о/-кальцифирах меньше, чем в ^/-кальцифирах, более высокие содержания щелочей отмечаются в оливинсодержащих разностях. Отличия содержаний других компонентов (в том числе С02, SiO2, А1203) незначительны.
б
>*-
>-
гхз
8
9
10
12 13
Рис. 1. Геологическая карта района Порьей губы (а) и план детального участка на о-ве Медвежий (б) [22], с изменениями: а - Порья губа: 1 - четвертичные отложения; 2-7 - супракрустальные образования: 2 - силлиманит-гранат-биотитовые
гнейсы, 3 - кислые гранулиты, 4 - плагиоклаз-ортопироксеновые кристаллические сланцы, 5 - двупироксеновые кристаллические сланцы, 6 - плагиоклаз-гранат-пироксеновые кристаллические сланцы, 7 - гранатовые амфиболиты;
8, 9 - гнейсовый комплекс обрамления гранулитов: 8 - биотит-амфиболовые гнейсы, 9 - гнейсы (нерасчлененные); 10, 11 - диафториты по гранулитам: 10 - гранатсодержащие гнейсы с реликтовыми пироксенами, 11 - гранат-пироксеновые амфиболиты; 12-19 - магматические образования: 12 - друзиты (13 - плагиоклазиты, метагаббро-анортозиты,
14 - нориты и пироксениты, 15 - эндербиты, 16 - чарнокиты, 17 - роговообманковые перидотиты, 18 - нерасчлененные ультрабазиты, 19 - щелочные гнейсо-граниты); 20-24 - другие породы: 20 - мигматиты, 21 - кальцифиры, 22 - конгломераты, 23 - эруптивная брекчия, 24 - тектоническая брекчия с сульфидным цементом; 25 - тектонические нарушения; 26 - элементы залегания (а - сланцеватость, б - линейность); 27 - участки отбора проб
(о-в Медвежий, о-в Тамаркина луда);
б - участок о-ва Медвежий: 1 - пегматиты; 2 - кислые гранулиты; 3 - чарнокит-мигматиты; 4 - амфибол-диопсидовые плагиосланцы с гранатом; 5 - гранат-двупироксеновые плагиосланцы; 6 - плагиоклаз-гиперстеновые кристаллические сланцы; 7 - кварциты; 8 - кальцифиры; 9 - биотитовые гнейсы с силлиманитом; 10 - клинопироксениты; 11 - актинолитовые сланцы; 12 - габбро-диабазы; 13 - точки отбора образцов (табл. 1)
Таблица 1
Валовый состав кальцифиров и пространственно ассоциированных кристаллических сланцев Порьей губы
Порода ¿г-кальцифиры о/-кальцифиры Кристаллические сланцы
Номер образца
7 28 13 8 1 3 12 9 10
SiO2 29,97 20,36 19,97 33,07 19,15 22,05 31,36 57,53 52,61
TiO2 0,27 0,07 0,17 0,13 0,20 0,61 0,47 0,76 0,67
A12O3 4,04 2,69 0,10 2,55 1,45 1,55 2,20 13,41 9,66
FeO 5,77 4,40 3,77 2,91 2,56 2,89 1,71 7,66 10,56
MnO 0,39 0,23 0,38 0,18 0,33 0,63 0,76 0,09 0,15
MgO 2,55 1,87 1,69 18,52 19,52 17,61 16,62 6,87 11,32
CaO 36,87 42,33 43,25 25,38 28,87 29,11 25,55 5,94 10,12
Na2O 0,25 0,15 0,27 0,31 0,67 1,11 1,29 3,94 2,29
K2O 0,10 0,62 0,05 0,10 0,59 0,99 0,57 1,62 0,93
P2O5 0,18 0,09 0,20 0,05 0,04 0,05 0,33 0,05 0,07
CO2 19,60 25,77 28,00 16,79 20,63 20,02 17,54 - -
LOI - 0,29 - 0,97 4,70 3,82 0,44 2,12 1,59
Сумма 99,99 98,87 97,85 100,96 98,71 100,44 98,84 99,99 99,97
MgO/(MgO + CaO) 0,06 0,04 0,04 0,42 0,40 0,38 0,39 0,54 0,53
Примечания: обр. 13 - о-в Тамаркина луда, остальные - о-в Медвежий; LOI - потери при прокаливании (105 °С).
Изотопы кислорода и углерода. Измеренные изотопные отношения кальцита (5180, 513С относительно SMOW и PDB соответственно) и силикатов (5180) приведены в табл.2. Величина
5180 кальцита в кальцифирах варьирует от 15,0 до 17,8 %о (среднее значение 16,44±0,92 %о) при вариации 513С от -3,6 до -6,2 %о (среднее -4,90±0,89 %о), что свидетельствует в пользу осадочной природы пород. При этом 5180 кальцита из диопсидовых кальцифиров в среднем составляет 17,05 (±0,61) %о, тогда как в оливинсодержащих разновидностях 5180 кальцита содержание ниже и в среднем составляет 15,68 (±0,59) %о. Средние значения 513С кальцита составляют -4,35 (±0,65) и -5,58 (±0,67) %о соответственно. Наблюдается тенденция увеличения 5180 и 513С кальцита с увеличением его содержания в породе (которая проявляется, например, в виде линейной зависимости 518Оса/ ~ 0,15 С02 + 13,02 мас.%, Я2 = 0,7). В пределе, если порода состоит из чистого кальцита (содержание С02 44 мас.%), то 518Оса/, согласно этой зависимости, оценивается в 19,6 %о, что соответствует изотопному составу докембрийских диафторирован-ных карбонатных осадков [36, 37]. Другая отмеченная закономерность - взаимосвязь 5180 и 513С кальцита - может обуславливаться метаморфической дегазацией пород. Представительные изотопные анализы кальцифиров, отобранные для сравнения на участке Тамаркина луда, следуют той же закономерности.
Температурные зависимости фракционирования стабильных изотопов определены на основе так называемых Р-факторов (приведенных отношений статистических сумм):
ln р=N"1
Nq {q}
- I ln
N h
N ^ Sinh(hVg,i /2kT) ^
vVq, г Slnh(hv*q, г /2kT)
где vq, i - частоты колебаний с волновым вектором q и индексом фононной ветви i от 1 до 3Nat; Nat - количество атомов в единичной ячейке; T - температура, K; h, к - постоянные Планка и Больцмана соответственно; N - число атомов, подвергающихся изотопному замещению; Nq - количество векторов q, учитываемых при суммировании; индекс * относится к более тяжелому изотопу.
Р-факторы (кальцита, кварца, моноклинного пироксена, оливина и граната) в интервале от 0 до 1500 °C интерполированы кубическими полиномами от х = 106/T2:1000lnp180 = ax + bx2 + cX [29]. Частоты колебаний изотопологов изученных фаз определены методом замороженных фононов (программа CRYSTAL: https://www.crystal.unito.it/) в рамках теории функционала плотности (DFT).
Таблица 2
Изотопный состав мономинеральных фракций
Номер Измеренный изотопный состав (5180ш0№ 513СРш) Т, °С F С02(/)
образца 5180са/ 513аса/ 5180„й 5180ф, 5180„, 5180др qtz-ca/ о/-са/ срх-са/ Другие
Л'-кальцифиры
7 16,50 -4,62 17,14 13,10 - 14,47^ 990 - 620 856^?^ 0,70 27,96
28 17,54 -4,20 18,07 - - - 1110 - - - 0,90 28,54
13 17,78 -3,57 19,85 16,45 - 16,32^ 430 - 810 701 0,96 29,25
21 16,41 -5,31 17,85 - - - 570 - - - 0,69 39,69
34 17,03 -4,05 17,65 - - - 1010 - - - 0,80 38,24
о/-кальцифиры
8 15,02 -6,06 - - 12,04 - - 520 - - 0,49 34,35
1 16,31 -5,32 - - 12,51 - - 410 - - 0,67 30,82
3 16,03 -4,76 - - 14,75 - - 990 - - 0,63 32,02
Кристаллические сланцы
12 15,37 -6,18 - - 12,07 - - 470 - - 0,53 32,95
9 - - - 12,11 - 13,42р/ - - - 670р/-срх - -
10 - - - 11,60 - 13,38р/ - - - 480р/-срх - -
Значения коэффициентов температурной зависимости Р-факторов
Р180са/
11,6603 -0,3544 0,00908
Р18Од(г 12,6764 -0,3543 0,00841
Р18Осрт 10,0399 -0,205 0,00401
Р1800/ 9,50807 -0,1654 0,0018
Р180р/ 11,478 -0,457 0,0077
Р18Оаз2 15,1896 -0,7068 0,0206
Примечания: F - доля С02, оставшаяся в кальците после дегазации; С02(/) - содержание СО2 в породе до дегазации (см. табл.1). Для обр. 21, 34 измеренные содержания СО2 составляют 27,22 и 30,50 соответственно.
Значение 5180 силикатной части кальцифиров (кварц, пироксен, оливин, гранат, табл.2) определяется как изотопным составом карбонатной части (чем выше 5180са/, тем в целом выше 5180 каждой конкретной фазы), так и величиной фракционирования 180/160 между различными фазами в зависимости от температуры. Вариация 5180 кварца в ^/-кальцифирах составляет 17,1-19,9 %о, пироксена 13,1-16,5 %. Значения 5180 оливина из о/-кальцифиров меняется от 12,0 до 14,8 %. Уменьшение 5180 происходит в ряду кварц-кальцит-(гранат-пироксен). Изотопное равновесие между минералами не наблюдается, на что указывает большой разброс вычисленных температур по распределению изотопов между разными фазами каждой пробы (табл.2). Например, температуры, вычисленные с применением изотопного геотермометра кварц-кальцита, варьируют от 430 °С (обр. 13) до 990 °С (обр. 7) и 1110 °С (обр. 34), кальцит-клинопироксена - от 620 до 810 °С (по двум образцам).
В породах без карбонатов как у контактов с кальцифирами, так и в некотором удалении (5-10 м) также отмечаются повышенные 5180 (до 11,6-13,5 % пироксена и плагиоклаза).
Реконструкция первичного состава метаморфических пород. Величины изотопных отношений кислорода и углерода указывают на первично-осадочную природу как кальцифиров Порьей губы, так и, возможно, всей изученной толщи пород (северо-западное побережье Порьей губы). Поэтому проведена реконструкция фазового состава исходных осадков на основе химического состава пород. С помощью комплекса програм Регр1е_Х рассчитан возможный (модельный) фазовый состав пород в условиях накопления осадков (25 °С, 1 бар) (табл.3).
Осадки, по химическому составу соответствующие d/-кальцифирам, содержали 70-85 мас.% карбонатной фракции с преобладанием кальцита (до 78 %), с небольшим количеством доломита (до 10 %) и анкерита (до 15 %). Кластическая компонента представлена кварцем (до 22 мас.%), в составе глинистой компоненты (до 9 мас.%) преобладал каолин; содержание детритовой составляющей (рутил, титанит) 0,1-0,2 %. Исходный состав осадков, преобразованных в о/-кальцифиры (преимущественно форстеритовые), отличался существенной долей доломита в карбонатной фракции (до 62 мас.%), а также анкерита (до 8 %), сидерита (до 4 %) и магнезита (менее 1 %). При этом содержание кальцита не превышало 35 мас.%.
Таблица 3
Реконструированный минеральный состав пород в условиях осадконакопления, мас.%
Номер образца
Порода Л-кальцифиры о/-кальцифиры Кристаллические сланцы
7 28 13 8 1 3 12 9 10
Карбонатная фракция
Кальцит 44,3 74,4 78,0 - 35,1 32,4 20,4 - 5,3
Доломит 10,0 - - 61,7 24,9 20,5 33,1 15,4 2,0
Сидерит - - - 3,5 - - - 9,7 -
Магнезит - - - 0,4 - - - 4,3 -
Анкерит 14,9 - 5,5 - 7,2 7,9 4,6 - 3,3
Всего карбонатов 69,3 74,4 83,5 65,6 67,2 60,8 58,1 29,4 10,6
Обломочная фракция
Кварц 21,7 0,3 14,6 22,2 - - 9,2 43,3 5,8
Рутил 0,2 - 0,2 0,1 - 0,6 0,4 0,6 -
Клинопироксен - 4,8 - - - - - - -
Титанит - 0,2 - - - - - - 1,5
Цеолиты - 18,3 - - - - - - -
Ильменит - - - - 0,4 - - - -
Всего обломочной
фракции 21,9 23,6 14,9 22,3 0,4 0,6 9,6 43,9 7,3
Глинистая фракция, амфибол, слюда, другие филлосиликаты
Каолинит 8,8 - 0,8 4,8 - - - 26,7 -
Тальк - - - - - 1,9 5,3 - -
Глаукофан - - - - - 3,3 10,8 - 26,7
Хлорит - 2,0 0,8 - - - - - 2,9
Слюда - - - 7,4 - 2,1 - - 7,3
Серпентин - - - - 20,3 10,3 - - -
Стильпномелан - - - - 12,1 21,0 16,2 - -
Хлоритоид - - - - - - - - 10,2
Актинолит - - - - - - - - 35,0
Всего фракции 8,8 2,0 1,6 12,1 32,4 38,6 32,3 26,7 82,1
Примечания. Состав рассчитан по данным табл.1 при поверхностных условиях T = 298,15 °C, P = 1 бар (модуль Werami программы Perple_X [30]). Содержание CO2 рассчитано по измеренному составу с учетом степени дегазации F и приведено в табл.2; при отсутствии данных вычисления проводились в условиях насыщения летучими (H2O > 5 %).
В составе осадочных пород, в настоящее время представленных кристаллическими сланцами, содержание карбонатов не превышало 30 мас.% (реконструкция проведена при условии избытка CO2), вероятно, преобладали или обломочная фракция (до 43 % кварца, обр. 9), или глинистая фракция (обр. 10). Оценки вариации изотопного состава кислорода осадков при взаимодействии с морской водой при 25 °С составляют 31,8±2,3 %о (вероятный изотопный состав кислорода до ли-тификации осадков).
Фазовые равновесия кальцифиров при метаморфизме. Диаграммы фазовых равновесий с полями устойчивости фазовых ассоциаций в координатах P-T (псевдосекции) построены с учетом химического состава представительных проб диопсидовых кальцифиров о-ва Медвежий и о-ва Тамаркина луда, оливиновых кальцифиров и гранат-пироксеновых плагиосланцев (см. табл.1). Построение псевдосекций дает оценку P-T-условиям фазовых преобразований и, в частности, вероятную температуру выделения компонентов флюида, что существенно при реконструкции изотопного состава пород. Учитывая близкое пространственное расположение исследованных пород, вероятны сходные P-T-условия метаморфизма, которые должны соответствовать пересечению полей стабильности наблюдаемых в разных образцах минеральных ассоциаций. Для наглядности результаты расчета псевдосекций представлены без учета компонент, присутствующих в незначительных количествах (рис.2).
8000
7000 ■
я
ю 6000 а,"
5000 ■
4000
500 600 700 800
Т, °с
900 1000
8000 1
7000 ■
6000 ■
&
ю
сС
5000 -
4000 ■
3000
grs-ttn-cpx(di+hed)--qtz-cal(CO2)
300 400 500 600 700 800 900
т, °с
ю
8000
7000
6000
5000
4000
3000
500
600
700
т оС
800
900
б
а
в
Рис.2. Фазовые P-T-диаграммы (псевдосекции), построенные по результатам анализа представительных кальцифиров
Порьей губы:
а: обр. 7. Безводная проекция. Устойчивые фазы: 1 - alm-grs-ttn-di-qtz-cal-dol; 2 - alm-ttn-di-qtz-hc-cal-dol; 3 - alm-ttn-di-qtz-spr-cal-dol; 4 - alm-ttn-di-qtz-spr-cal-CO2; 5 - alm-di-qtz-ilm-hc-cal-dol; 6 - alm-di-qtz-ilm-hc-cal-CO2; 7 - alm-crd-di-qtz-ilm-cal-dol; 8 - alm-crd-di-qtz-ilm-cal-CO2; 9 - alm-ttn-di-qtz-hc-cal-CO2; 10 - alm-grs-ttn-di-qtz-cal-CO2;
11 - grs-ttn-di-hed-wo-cal-CO2; б: обр. 13. Безводная проекция. Устойчивые фазы: 1 - grs-ttn-di-hed-qtz-arg-ank; 2 - grs-ttn-di-hed-qtz-cal-ank; 3 - alm-ttn-di-hed-qtz-cal-ank; 4 - grs-ilm-di-hed-qtz-cal-ank; 5 - alm-ilm-di-hed-qtz-cal-ank; 6 - alm-ilm-di-fa-qtz-cal-ank; 7 - alm-ilm-di-fa-qtz-cal-COi; 8 - alm-ilm-di-hed-qtz-cal-COi; 9 - alm-ttn-di-hed-qtz-cal-COi; 10 - sil-ttn-di-hed-qtz-cal-COi;
11 - grs-ttn-di-hed-wo-cal-COi; в: обр. 8. Проекция, включая H2O. Устойчивые фазы: 1 - alm-ky-di-tr-qtz-ilm-mgs-dol; 2 - alm-an-di-tr-qtz-ilm-mgs-dol; 3 - alm-en-di-tr-an-qtz-ilm-dol; 4 - alm-en-di-act-an-qtz-ilm-dol; 5 -fa-en-di-an-qtz-ilm-dol-act; 6 - fa-en-di-an-ilm-dol-COi-act; 7 - alm-en-di-an-ilm-dol-CO2-act; 8 -fa-en-di-tr-an-ilm-dol-CO2; 9 - alm-en-di-tr-an-ilm-dol-CO2; 10 -fo-fa-di-tr-an-ilm-dol-CO2; 11 - alm-en-di-tr-ilm-hc-dol-COi; 12 -fa-en-di-tr-ilm-hc-dol-CO2; 13 -fo-fa-di-tr-ilm-hc-dol-CO2; 14 -fo-fa-di-clin-ilm-hc-dol-CO2; 15 -fo-fa-di-clin-ilm-hc-cal-COi; 16 -fo-fa-di-ilm-hc-dol-H2O-CO2; 17 (наблюдаемая ассоциация) -fo-fa-di-ilm-hc-cal-H2O-CO2;
18 -fo-fa-di-hc-usp-cal-H2O-CO2 Обозначения фаз. Карбонаты: ank - анкерит; cal - кальцит; dol - доломит; mgs - магнезит; sid - сидерит; arg - арагонит. Силикаты: qtz - кварц; ol - оливин (fo - форстерит, fa - фаялит); cpx - моноклинный пироксен (di - диопсид, hed - геденбергит); opx - ромбический пироксен (en - энстатит, hyp - гиперстен); grt - гранат (alm - альмандин, grs - гроссуляр); pl - плагиоклаз (an - анортит); spr - сапфирин; crd - кордиерит; sil - силлиманит; ky - кианит; wo - волластонит; amp - амфибол (act - актинолит, gl - глаукофан, hrn - роговая обманка, tr - тремолит); слюда (6t - биотит, mu - мусковит, phl - флогопит); окислы: spl - шпинель (hc - герцинит, usp - ульвошпинель); rt - рутил;
ttn - титанит; ilm - ильменит 1 - линии моновариантных равновесий, ограничивающих появление CO2; 2 - линии, ограничивающие появление кальцита
Реакции декарбонатизации, которые сопровождают фазовые преобразования всех изученных кальцифиров, заслуживают особого внимания для объяснения изменения изотопного состава пород, так как данные реакции вызывают изотопное фракционирование как кислорода, так и углерода; факторы изотопного фракционирования между С02 и твердыми фазами значимы даже при повышенных температурах; большое увеличение объема, связанное с образованием компонентов флюида, вызывает полное или частичное удаление С02.
Для d/-кальцифиров выделение С02 обуславливается реакциями разложения доломита или анкерита, который становится неустойчивым при повышении температуры:
CaMg(C0з)2(do/) + 8102 = СаС0э(са/) + MgSl0з(cpx) + С02; (1)
CaMg(C0з)2(do/) + 28102 = > CaMgSl0з(cpx) + 2С02. (2)
Если при низкой температуре в составе исходной породы преобладают карбонаты железа (анкерит, обр. 13), то выделение углекислого газа происходит, например, по реакции
CaFe(C0з)2 ^¡с1) + 8102 = СаС03(са/) + FeSl0з(cpx) + С02. (3)
В качестве источника кремнезема может выступать как кварц, так и другие силикаты, в том числе следующие:
CaMg(C0з)2(do/) + FeзAl2(Sl04)з(a/ш) = = СаС0з(са/) + 2FeSl0з(cpx) + MgSl0з(cpx) + FeAl204(sp/) + С02; (4)
3CaMg(C0з)2(do/) + CaMg(SiOз)2(d0 = 4СаС0з(са/) + 2Mg2S104(/0) + 2С02. (5)
В d/-кaльцифирaх о-ва Медвежий (обр. 7) при давлении до 6 кбар реакции декарбонатизации начинаются в интервале 580-630 °С (рис.2, а, красная линия). Для наглядности показана упрощенная безводная проекция системы). При более высоком давлении температура начала выделения С02 также увеличивается (например, если Р = 8 кбар, Т начала выделения С02 сдвигается до ~700 °С). В результате последовательности фазовых преобразований из осадочной породы (см. табл.1) возникает наблюдаемая ассоциация grs-ttn-di-hed-qtz-ca/, которая при давлении 8 кбар устойчива от 790 до 970 °С (при другом давлении поле устойчивости сдвигается приблизительно на 30 °С/кбар). При более низкой температуре вместо альмандина устойчив геденбергит, а при более высокой температуре вместо кварца появляется волластонит.
Ассоциация di-кaльцифиров о-ва Тамаркина луда (обр. 13): са/-срх(И + hed)-grt(grs)-qtz-ttn. Начало декарбонатизации при пониженном давлении возможно при температуре 560 °С (Р = 3 кбар). При более высоком давлении температура появления С02 увеличивается до 705 °С (Р = 8 кбар), при этом градиент дТ/дР оценивается приблизительно в 29 °С/кбар (рис.2, б). Как и для диопсидовых кальцифиров о-ва Медвежий (обр. 7), поле стабильности наблюдаемой ассоциации grs-ttn-cpx-qtz-ca/ при повышении температуры ограничивается ассоциацией с волластонитом вместо кварца. Примечательно, что для данной породы равновесные карбонаты до температуры приблизительно 570-650 °С состоят из кальцита и анкерита, а не доломита (как в обр. 7) и, соответственно, декарбонатизация обуславливается разложением анкерита по реакции, подобной (3).
0/-кальцифиры (обр. 8, рис.2, в). Наблюдаемая ассоциация/в-/а^Ы/т^с-са/ находится в равновесии с компонентами флюида (Н20 + С02), начиная от температуры 690 °С при давлении 3 кбар или 860 °С при давлении 8 кбар; верхняя граница устойчивости 750-905 °С соответственно. Более низкотемпературная ассоциация включает доломит, реликты которого наблюдаются в шлифах, в более высокотемпературной ассоциации fo-fa-di-usp-hс-ca/-(H20-а02) за счет ильменита образуется ульвошпинель. Можно отметить, что в отличие от диопсидовых кальци-фиров, начало декарбонатизации (выделено красной линией на рис.2, в) не приводит к полному исчезновению доломита, и кальцит становится устойчивым только при температуре от 680 °С (3 кбар) до 860 °С (8 кбар), что совпадает с нижней температурной границей устойчивости наблюдаемой ассоциации (рис.2, в, синяя линия).
8000
7000
&
ю
сС
6000
5000 -
4000 ■
3000
600
700
800 Т, °с
900
1000
Рис.3. Поля устойчивости минеральных ассоциаций
кальцифиров и кристаллических сланцев. Область пересечения всех наблюдаемых ассоциаций заштрихована 1 - gr?-px-слaнцы; 2 - о/-кальцифиры; 3 - di-кaльцифиры
О
сл
ю
-3
-4
-5
-6
-7 ■
-8
513Сса/ = 0,90815180са/ - 19,828
К2 = о,:
700
14
15
16
5180сг4
17
18
Рис.4. Изотопные отношения 180/160 и 13С/12С кальцита из кальцифиров Порьей губы 1 - фигуративные точки анализов; 2 - тренд, построенный по результатам изотопных анализов; 3 - результаты вычислений по уравнению (6) для изотопного фракционирования (значения кислородного фактора фракционирования а18/16 са/-С02 при Т = 600, 700 и 900 °а составляют 1,00411, 1,00311 и 1,00242 соответственно; значения а18/16 dol-а02 - 1,00303, 1,0025, 1,0018 соответственно [29]; значения углеродного фактора фракционирования а13/12 са/-С02 при заданных Т - 1,0035, 1,0033, 1,0027; а13/12 dol-C02 - 1,00306, 1,0029, 1,00246 [29]);
4 - вычисления при изотопном равновесии С02 с карбонатами и силикатами при 600 °С; 5 - вычисления в условиях изотопного равновесия С02 с доломитом при Т = 500, 700 °С
Р-Т-условия, при которых стабильны все наблюдаемые минеральные ассоциации, определяются пересечением полей равновесий по разным породам (рис.3). При минимальных возможных давлениях (5,5-6 кбар) температура, ограничивающая стабильность всех ассоциаций, составляет 740-870 °С. Если давление метаморфизма выше, оценки температуры возрастают (дТ/дР ~ 30-40 °С/кбар), и при Р = 8 кбар интервал температуры, при котором могут существовать все наблюдаемые ассоциации, составляет 860-910 °С.
Обсуждение. Изменение изотопного состава при метаморфизме. Наблюдается общая (для образцов всех разновидностей изученных кальцифи-ров) значимая корреляция между значениями 513С и 5180 кальцита; Я2 = 0,86 (рис.3). Реакции выделения С02, характерные при метаморфических преобразованиях предшественников кальцифиров, сопровождаются фракционированием изотопов кислорода и углерода с обеднением рестита (остатка твердых фаз после выделения газа) тяжелыми изотопами (10001па между С02 и карбонатами и (или) силикатами всегда больше нуля при равновесном фракционировании изотопов кислорода и температуре выше ~150-200 °С при фракционировании изотопов углерода [29]). При условии постоянного объема вероятно непрерывное удаление образующегося С02 из системы, так что изменение изотопного состава твердых фаз можно оценить по уравнению дистилляции Рэлея, которое для изотопного фракционирования выражается следующим соотношением [38]:
5/ -5гг = 1000(^а - Х) - 1), (6)
где F - мольная доля элемента, которая остается в рестите после выделения газа; 5 - изотопные отношения элемента (кислорода, углерода) в твердой фазе; а - фактор изотопного фракционирования между флюидом и твердой фазой -
(180/160)флюид/(180/160)рестит; / / - величины после выделения газа и до начала выделения газа
соответственно.
Линейный тренд одновременного изменения изотопного состава кислорода и углерода согла-
суется с аналитической зависимостью 513Сса/ и
18
5 0са/ (рис.4) при следующих условиях: при дегазации поддерживалось изотопное равновесие С02 и образующегося карбоната (кальцита и (или) доломита); изотопные реакции С02 с карбонатами, вероятно, не сопровождались изотопными реакциями с силикатами. Второе допущение можно обосновать тем, что тренд, рассчитанный для изотопного равновесия с участием силикатов, не соответствует аналитическим данным (пунктирная линия на рис.4). Кроме того, измеренное
распределение изотопов кислорода между кальцитом и силикатами (см. табл.2) существенно отклоняется от равновесного. Без изотопного равновесия с силикатами величины F для кислорода и углерода совпадают (FC = FO [39]). Величина изотопного фракционирования cal-CO2 соответствует температуре 700 °C, что полностью подтверждается оценками температур декарбонатиза-ции, полученными в результате моделирования фазовых соотношений (с построением P-T-псев-досекций).
Полученная усредненная зависимость объединяет аналитические данные по всем исследованным образцам, включая диопсидовые кальцифиры и оливиновые кальцифиры. Некоторые отклонения аналитических точек могут обуславливаться изменением величин фракционирования (вариациями температур декарбонатизации, изменением состава карбонатов в равновесии с CO2 оливиновых кальцифиров). В первую очередь, однако, измеренные вариации 518O и 513C обуславливаются степенью дистилляции F соотношения (6) при едином изотопном составе карбоната на начальной стадии дистилляции (F = 1, начало выделения CO2) и близкой температуре дегазации (при T = 700 °C значения 518Ocal = 17,96, 513Ccal = -3,4 %о хорошо воспроизводят все аналитические данные).
Пониженные значения 518O кальцита d-кальцифиров относительно 518O кальцита Ji-кальцифиров с общим трендом 518O vs 513C можно интерпретировать как следствие большей степени дегазации оливиновых кальцифиров. Так, максимальные значения 518Ocal = 17,78 и 513Ccal = -3,57 % (обр. 13) отражают начальные стадии дистилляции (F = 0,96), тогда как минимальные значения 518Ocal = 15,02 и 513Ccal = -6,06 % (обр. 8) отражают более продвинутую дистилляцию (F = 0,5). Данные результаты согласуются с процессом декарбонатизации с разложением доломита (или Fe-карбоната, обр. 13) и образованием кальцита (см. рис.2, б) в результате серий последовательных реакций с удалением выделяющегося CO2. Температура дегазации, полученная по результатам изотопного анализа, соответствует результатам фазового моделирования (T = 700 °С при P ~ 8 кбар).
Хорошее совпадение аналитических данных на диаграмме 513C vs 518O с моделью дистилляции Рэлея показывает, что изотопный состав удаляемого CO2 и остающегося карбоната контролировался равновесным обменом между остающимся после выделения углекислого газа карбонатом (кальцитом) и непрерывно удаляемым газом. При этом изотопного обмена с силикатами не было. Изотопный состав образующегося CO2 в равновесии с кальцитом варьировал в зависимости от степени дегазации от 518O = 21,9 и 513C = -1,0 % (F ~ 1) до 18,9 и -3,7 % (F ~ 0,5) соответственно. Интегрированный (batch) изотопный состав образующегося CO2 находится между этими значениями и определяется соотношением материального баланса. При незначительном изменении величины фракционирования в ходе декарбонатизации изотопный состав всего образовавшегося CO2 оценивается величинами 518O = 19,9 и 513C = -2,2 % (F = 0,5).
Равновесие выделяющегося газа с кальцитом без участия силикатов [40] не ограничивает F так называемым кальций-силикатным пределом (F = 0,6 [38]), и минимальные значения измеренных изотопных отношений соответствуют степени дегазации около 0,5 (ol-кальцифиры).
Заключение. Результаты изотопного анализа кислорода и углерода показывают, что кальци-фиры Порьей губы и переслаивающиеся с ними кристаллические сланцы с повышенными значениями 518O и 513C сформировались при метаморфизме осадочных пород (известняков, мергелей, граувакков). По результатам литохимических реконструкций карбонатно-силикатные породы по составу соответствуют мергелистым осадкам, химический состав которых варьировал в зависимости от соотношения глинистой, карбонатной и обломочной фракций. Доля карбонатов составляла 58-84 мас.% (в осадках, из которых образовались кальцифиры) и 10-30 мас.% (в возможных про-толитах кристаллических сланцев). Изотопный состав осадков (518O 17,9 % SMOW и 513C -3,4 % PDB) находится в пределах значений, характерных для докембрийских диагенетически преобразованных осадочных карбонатов [36, 37, 41] без какого-либо влияния как органического углерода [25], так и углерода магматических образований [42].
Наиболее существенные изменения изотопного состава при метаморфизме пород обусловлены последовательностью реакций с образованием CO2 при разложении доломита или других карбонатов (сидерит, анкерит), которые становятся неустойчивыми при повышении температуры приблизительно до 700 °С (P ~ 8 кбар) и замещаются силикатами магния/железа (например,
моноклинным пироксеном, гранатом). В процессе дегазации пород 5180 и 513С кальцита одновременно понижались за счет изотопного фракционирования (и дистилляции Рэлея) до значений 15,0 и -6,2 % соответственно (при удалении до половины С02) с образованием линейной взаимосвязи 5180 vs 513С. Изотопный состав высвобождаемого (и, возможно, добавляемого в атмосферу) углекислого газа оценивается в пределах 20 % (SM0W) и -2 % (PDB).
Таким образом, изменение изотопного состава кальцита отражает разную степень дегазации, которая, в свою очередь, зависит от состава исходных пород. Например, чем выше содержание доломита до метаморфизма было в осадочной породе, тем выше степень дегазации. В оливиновых кальцифирах (Mg0 до 20 мас.%) значения 5180 и 513С ниже, чем в диопсидовых кальцифирах. При этом вариации изотопного состава кальцита переслаивающихся пород согласуются с внутренним (локальным) контролем выделяющегося флюида [40] и противоречат крупномасштабному взаимодействию с каким-либо внешним (глубинным или поверхностным) резервуаром. В том числе наши результаты не подтверждают возможность смешения (ассимиляции) исследованных карбонатных пород с магматическими и (или) метасоматические преобразования при образовании скарнов.
Можно отметить, что при температуре метаморфизма ниже 650-700 °С доломит остается устойчивым, и наблюдаемые доломитовые кальцифиры Колвицкой зоны гранулитов [24], по химическому составу сходные с оливиновыми кальцифирами, рассмотренными в настоящей работе, могут быть предшественниками оливиновых кальцифиров при меньших максимальных температурах метаморфизма.
ЛИТЕРАТУРА
1. Berner R.A., Lasaga A. C., Garrels R.M. The carbonate-silicate geochemical cycle and its effect on atmospheric carbon dioxide over the past 100 million years // American Journal of Science. 1983. Vol. 283. Iss. 7. P. 641-683. DOI: 10.2475/ajs.283.7.641
2. Groppo C., Rolfo F., Castelli D., Connolly J.A.D. Metamorphic CO2 production from calc-silicate rocks via garnet-forming reactions in the CFAS-H2O-CO2 system // Contributions to Mineralogy and Petrology. 2013. Vol. 166. P. 1655-1675. DOI: 10.1007/s00410-013-0947-5
3. Catling D.C., Kasting J.F. Atmospheric Evolution on Inhabited and Lifeless Worlds. Cambridge: Cambridge University Press, 2017. 579 p. DOI: 10.1017/9781139020558
4. Сорохтин Н.О., Никифоров С.Л., Козлов Н.Е. Коромантийная ветвь глобального цикла углерода и происхождение глубинных углеводородов // Вестник МГТУ. 2018. Т. 21. № 1. С. 61-79. DOI: 10.21443/1560-9278-2018-21-1-61-79
5. Brovarone A.V., Tumiati S., Piccoli F. et al. Fluid-mediated selective dissolution of subducting carbonaceous material: Implications for carbon recycling and fluid fluxes at forearc depths // Chemical Geology. 2020. Vol. 549. № 119682. DOI: 10.1016/j. chemgeo.2020.119682
6. Deep Carbon: Past to present / Ed. by Orcutt B.N., Daniel I., Dasgupta R. Cambridge: Cambridge University Press, 2020. 669 p. DOI: 10.1017/9781108677950
7. Stagno V. Carbon, carbides, carbonates and carbonatitic melts in the Earth's interior // Journal of the Geological Society.
2019. Vol. 176. Iss. 2. P. 375-387. DOI: 10.1144/jgs2018-095
8. Skublov S.G., Rumyantseva N.A., Qiuli Li et al. Zircon xenocrysts from the Shaka Ridge record ancient continental crust: New U-Pb geochronological and oxygen isotopic data // Journal of Earth Science. 2022. Vol. 33. № 1. P. 5-16. DOI: 10.1007/s12583-021-1422-2
9. Левашова Е.В., Скублов С.Г., Ойцева Т.А. и др. Изотопно-геохимические особенности циркона из постколлизионных гранитов: на примере рибекитовых гранитов Верхнее Эспе, Восточный Казахстан // Геохимия. 2022. Т. 67. № 1. С. 3-18. DOI: 10.31857/S0016752522010083
10. Румянцева Н.А., Скублов С.Г., Ванштейн Б.Г. и др. Циркон из габброидов хребта Шака (Южная Атлантика): U-Pb возраст, соотношение изотопов кислорода и редкоэлементный состав // Записки Российского минерального общества. 2022. Ч. CLI. № 1. С. 44-73. DOI: 10.31857/S0869605522010099
11. Satish-Kumar M., Miyamoto T., Hermann J. et al. Pre-metamorphic carbon, oxygen and strontium isotope signature of highgrade marbles from the Lutzow-Holm Complex, East Antarctica: apparent age constraints of carbonate deposition // Geological Society, London, Special Publications. 2008. Vol. 308. Iss. 1. P. 147-164. DOI: 10.1144/SP308.7
12. Thaworndumrongsakul P., Booth J., Nantasin P., Kim Y. Metamorphic Evolution of Calc-silicate Body Enclosed in Char-nockitic Gneiss at West Ongul Island, Lutzow-Holm Complex, East Antarctica // IOP Conference Series: Earth and Environmental Science. 2021. Vol. 837. № 012015. DOI: 10.1088/1755-1315/837/1/012015
13. Jacobs J., Mikhalsky E., Henjes-Kunst F. et al. Neoproterozoic geodynamic evolution of easternmost Kalahari: Constraints from U-Pb-Hf-O zircon, Sm-Nd isotope and geochemical data from the Schirmacher Oasis, East Antarctica // Precambrian Research.
2020. Vol. 342. № 105553. DOI: 10.1016/j.precamres.2019.105553
14. Гусев Н.И., Сергеева Л.Ю., Ларионов А.Н., Скублов С.Г. Реликты эоархейской континентальной коры Анабарского щита, Сибирский кратон // Петрология. 2020. Т. 28. № 2. С. 115-138. DOI: 10.31857/S086959032002003X
15. Гусев Н.И., Сергеева Л.Ю., Скублов С.Г. Свидетельства субдукции палеопротерозойской океанической коры в Хапчанском поясе Анабарского щита Cибирского кратона // Петрология. 2021. Т. 29. № 2. С. 115-135. DOI: 10.31857/S0869590321020047
16. ИльченкоВ.Л., АфанасьеваЕ.Н., Каулина Т.В. и др. Лицевское рудопроявление урана (арктическая зона Фенноскан-динавского щита): новые результаты петрофизических и геохимических исследований // Записки Горного института. 2022. Т. 255. С. 393-404. DOI: 10.31897/PMI.2022.44
17. Кудряшов Н.М., Удоратина О.В., Калинин А.А. и др. U-Pb (SHRIMP-RG) возраст циркона из редкометалльных (Li, Cs) пегматитов месторождения Охмыльк зеленокаменного пояса Колмозеро-Воронья (северо-восток Фенноскандинавского щита) // Записки Горного института. 2022. Т. 255. С. 448-454. DOI: 10.31897/PMI.2022.41
18. Смолькин В.Ф., Мокрушин А.В., Баянова Т.Б. и др. Магмаподводящий палеоканал в Мончегорском рудном районе: геохимия, изотопный U-Pb и Sm-Nd анализ (Кольский регион, Россия) // Записки Горного института. 2022. Т. 255. С. 405-418. DOI: 10.31897/PMI.2022.48
19. Салимгараева Л.И., Скублов С.Г., Березин А.В., Галанкина О.Л. Фальбанды Керетского архипелага Белого моря: характеристика состава пород и минералов, рудная минерализация // Записки Горного института. 2020. Т. 245. С. 513-521. DOI: 10.31897/ PMI.2020.5.2
20. Melnik A.E., Skublov S.G., Rubatto D. et al. Garnet and zircon geochronology of the Paleoproterozoic Kuru-Vaara eclogites, northern Belomorian Province, Fennoscandian Shield // Precambrian Research. 2021. Vol. 353. № 106014. DOI: 10.1016/j.precamres.2020.106014
21. Скублов С.Г., Березин А.В., Салимгараева Л.И. Эклогиты Беломорского подвижного пояса: геолого-петрологические и изотопно-геохимические критерии возраста // Геохимия. 2022. Т. 67. № 7. С. 621-638. DOI: 10.31857/S0016752522070044
22. КозловН.Е., Иванов А.А., Нерович Л.И. Лапландский гранулитовый пояс - первичная природа и развитие. Апатиты: КНЦ АН СССР, 1990. 168 с.
23. Виноградов Л.А., Богданова М.Н., Ефимов М.М. Гранулитовый пояс Кольского полуострова. Ленинград: Наука, 1980. 208 с.
24. Сафронов В.Т., Розен О.М. Метакарбонатные породы (кальцифиры) Лапландско-Колвицкого гранулитового пояса Балтийского щита // Литология и полезные ископаемые. 2004. № 5. С. 491-503.
25. Ивлиев А.И. Геология метаморфических комплексов Лапландского гранулитового пояса (Сальные тундры, Кольский полуостров): Автореф. дис. ... канд. геол.-минерал. наук. М.: ИМГРЭ, 1977. 23 с.
26. McCrea J.M. On the Isotopic Chemistry of Carbonates and a Paleotemperature Scale // The Journal of Chemical Physics. 1950. Vol. 18. № 6. P. 849-857. DOI: 10.1063/1.1747785
27. Clayton R.N., Mayeda T.K. The use of bromine pentafluoride in the extraction of oxygen from oxides and silicates for isotopic analysis // Geochimica et Cosmochimica Acta. 1963. Vol. 27. Iss. 1. P. 43-52. DOI: 10.1016/0016-7037(63)90071-1
28. Sang-Tae Kim, O 'Neil J.R. Equilibrium and nonequilihrium oxygen isotope effects in synthetic carbonates // Geochimica et Cosmochimica Acta. 1997. Vol. 61. № 16. P. 3461-3475. DOI: 10.1016/s0016-7037(97)00169-5
29. Крылов Д.П. Фракционирование стабильных изотопов в Сa-Mg карбонатах: расчет Р-факторов методом «замороженных фононов» // Геохимия. 2022. Т. 67. № 10. С. 942-960. DOI: 10.31857/S0016752522100065
30. Connoly J.A.D. Multivariable phase-diagrams - an algorithm based on generalized thermodynamics // American Journal of Science. 1990. Vol. 290. Iss. 6. P. 666-718. DOI: 10.2475/ajs.290.6.666
31. Connolly J.A.D. A Primer in Gibbs Energy Minimization for Geophysicists // Petrology. 2017. Vol. 25. № 5. P. 526-534. DOI: 10.1134/S0869591117050034
32. HuaXiang, Connolly J.A.D. GeoPS: An interactive visual computing tool for thermodynamic modelling of phase equilibria // Journal of Metamorphic Geology. 2022. Vol. 40. Iss. 2. P. 243-255. DOI: 10.1111/jmg.12626
33. Holland T.J.B., Powell R. An internally consistent thermodynamic data set for phases of petrological interest // Journal of Metamorphic Geology. 1998. Vol. 16. № 3. P. 309-343. DOI: 10.1111/j.1525-1314.1998.00140.x
34. Бибикова Е.В., Мельников В.Ф., Авакян К.Х. Лапландские гранулиты: петрология, геохимия и изотопный возраст // Петрология. 1993. Т. 1. № 2. С. 215-234.
35. МинцМ.В., ГлазневВ.Н., Конилов А.Н. и др. Лапландско-Колвицкий гранулито-гнейсовый пояс // Ранний докембрий северо-востока Балтийского щита: палеогеодинамика, строение и эволюция континентальной коры. М.: Научный мир, 1996. С. 112-138.
36. SchidlowskiM, Eichmann R., Junge C.E. Precambrian sedimentary carbonates: carbon and oxygen isotope geochemistry and implications for the terrestrial oxygen budget // Precambrian Research. 1975. Vol. 2. Iss. 1. P. 1-69. DOI: 10.1016/0301-9268(75)90018-2
37. Shields G., Veizer J. Precambrian marine carbonate isotope database: Version 1.1 // Geochemistry, Geophysics, Geosystems. 2002. Vol. 3. № 6. P. 1-12. DOI: 10.1029/2001GC000266
38. Baumgartner L.P., Valley J. W. Stable Isotope Transport and Contact Metamorphic Fluid Flow // Reviews in Mineralogy and Geochemistry. 2001. Vol. 43. № 1. P. 415-467. DOI: 10.2138/gsrmg.43.1.415
39. Valley J. W. Stable isotope geochemistry of metamorphic rocks // Stable Isotopes in High Temperature Geological Processes. 1986. P. 445-490. DOI: 10.1515/9781501508936-018
40. Lattanzi P., Rye D.M., Rice J.M. Behavior of 13C and 18O in carbonates during contact metamorphism at Marysville, Montana: implications for isotope systematics in impure dolomitic limestones // American Journal of Science. 1980. Vol. 280. P. 890-906. DOI: 10.2475/ajs.280.9.890
41. Graf D.L. Geochemistry of Carbonate Sediments and Sedimentary Carbonate Rocks: Part IA-A. Isotopic Composition Chemical Analysis. Urbana: Illinois State Geological Survey, 1960. 42 p.
42. Hoefs J. Stable Isotope Geochemistry. Springer, 2015. 389 p. DOI: 10.1007/978-3-319-19716-6
Авторы: Д.П.Крылов, д-р геол.-минерал. наук, ведущий научный сотрудник, [email protected], https://orcid.org/0000-0001-6654-8659 (Институт геологии и геохронологии докембрия РАН, Санкт-Петербург, Россия), Е.В.Климова, младший научный сотрудник, https://orcid.org/0000-0002-9771-9518 (Институт геологии и геохронологии докембрия РАН, Санкт-Петербург, Россия).
Авторы заявляют об отсутствии конфликта интересов.
