Петрология, 2020, T. 28, № 1, стр. 23-54

Флюиды гранулитов высоких давлений

С. А. Бушмин a*, Е. А. Вапник b**, М. В. Иванов a, Ю. М. Лебедева a, Е. В. Савва a

a Институт геологии и геохронологии докембрия РАН
199034 Санкт-Петербург, наб. Макарова, 2, Россия

b Department of Geological and Environmental Sciences, Ben-Gurion University of the Negev
84105 Beer-Sheva, Israel

* E-mail: s.a.bushmin@ipgg.ru
** E-mail: vapnik@bgu.ac.il

Поступила в редакцию 26.02.2019
После доработки 24.04.2019
Принята к публикации 27.05.2019

Полный текст (PDF)

Аннотация

Особенности флюидов гранулитов высоких давлений изучены на примере НР гранулитов (~8.7–11 кбар, ~800–900°С) и сингранулитовых инфильтрационных НР метасоматитов (~11–9 кбар, ~920–850°С) Лапландского гранулитового пояса Фенноскандинавского щита. Исследования включали крупномасштабное картирование, микротермометрию флюидных включений, мультиравновесную термобарометрию, расчет активности Н2О по минеральным равновесиям. В основных пироксеновых гранулитах и сингранулитовых метасоматитах (кварцевые бластомилониты с ортопироксеном, силлиманитом, гранатом; жилы и жилоподобные тела ортопироксен-гранатовых и диопсид-скаполитовых пород) обнаружены одинаковые ассоциации сингенетичных типов флюидных включений (преимущественно в кварце, также в гранате, ортопироксене, скаполите) контрастного химического состава: почти чистые включения СО2 (заметно преобладают), рассолы (главные соли CaCl2 и NaCl) и N2 ± H2O. Эти три типа включений сосуществуют в одних и тех же генерациях ранних включений: более редких первичных (p) и преобладающих первично-вторичных (ps). Включения СО2 – высокой и низкой плотности, включения N2 – низкой плотности. Включения рассолов обнаруживают широкий диапазон общего содержания солей (до 30–35 мас. %) и разное соотношение концентраций главных солей: p-включения с соленостью 20 мас. % CaCl2 + + 10 мас. % NaCl; ps-включения с соленостью 5 мас. % CaCl2 + 20 мас. % NaCl; p- и ps-включения с соленостью 5–23 мас. % в экв. NaCl; p-включения с галитом (до 35 мас. % NaCl). В целом CaCl2 является преобладающим компонентом среди солей в группе ранних p- и ps-включений изученных пород. Совокупность имеющихся данных (в том числе, изотопные системы Sr, Nd, кислорода) свидетельствует о том, что при гранулитовом петрогенезисе внешний флюидный поток, по-видимому, был мантийного происхождения. На пике Р-Т параметров захват включений происходил из гетерогенного флюида, в котором одновременно сосуществовали несмесимые водно-солевой и богатый СО2 флюиды, изначально содержащие азот. Данные по химическому составу включений и концентрациям солей, аН2О = 0.40–0.51 сравниваются с теоретическим предсказанием фазового состояния флюида и свойств сосуществующих несмесимых флюидных фаз при установленных Р-Т параметрах гранулитового петрогенезиса на основе численных моделей тройных систем H2O–CO2–NaCl и H2O–CO2–CaCl2. Для цели термодинамического предсказания рассчитаны сольвусы и положение коннод. Обсуждаются сходство и причины отличия теоретических составов возникающих флюидных фаз от состава флюидных включений, геохимические следствия, связанные с фазовой сепарацией гранулитового флюида (возникновение щелочных концентрированных рассолов и потенциально кислотной богатой СО2 флюидной фазы, величины массовых и объемных долей этих фаз в зависимости от вариаций состава исходного гомогенного флюида и др.). Таким образом, в диапазоне Р-Т параметров НР гранулитов существует обширная область составов Н2О флюидов с разным содержанием СО2 и хлоридов Na и Са, в которой гомогенный флюид распадается на контрастные по химическому составу и свойствам флюидные фазы, и эта область значительно расширяется с увеличением содержания CaCl2. Следовательно, нижняя кора уровня НР гранулитовой фации может являться областью зарождения высокотемпературных несмесимых флюидов. С одной стороны, это более плотная флюидная фаза относительно щелочных рассолов, с другой – менее плотная потенциально кислотная фаза Н2О–СО2 флюидов, богатых СО2. Поступая по региональным проницаемым зонам наверх, эти флюидные фазы глубинного происхождения могут играть важнейшую роль в магматическом, метаморфическом, метасоматическом и рудном петрогенезисе средней и верхней коры.

Ключевые слова: HP гранулиты, сингранулитовые НР метасоматиты, флюидные включения, рассолы, азот, активность воды, двухфазный флюид

ВВЕДЕНИЕ

Гранулиты, присутствующие во многих метаморфических поясах, особенно докембрийских щитов, составляют одну из главных частей нижней коры. Поэтому роль и всесторонние характеристики флюидов, участвовавших в глубинном НТ-НР петрогенезисе, – предмет продолжающихся дискуссий (например, Aranovich, Newton, 1996; Shmulovich, Graham, 1996; Newton et al, 1998; Touret, Huizenga, 2012; Aranovich et al., 2016; Safonov et al., 2018). Ни у кого не вызывает сомнений тот факт, что внешний флюидный поток участвует в гранулитовом петрогенезисе. Однако часто, как очевидное, рассматривается СО2 флюид и недооценивается роль водно-солевых флюидов из-за распространенных в гранулитах богатых СО2 флюидных включений. На протяжении длительного периода времени для гранулитового метаморфизма предполагался флюид Н2О–СО2 с низкой активностью воды за счет очень высокого содержания углекислоты (например, Newton et al., 1980). Позднее активно развивалось представление о водно-солевом флюиде с низким содержанием углекислоты, в котором низкая активность Н2О связана с повышенным содержанием солей, причем не только при умеренных, но и при достаточно высоких давлениях (Aranovich, Newton, 1996; Shmulovich, Graham, 1996; Newton et al., 1998; Аранович, 2007; Newton, Manning, 2010; Кориковский, Аранович, 2010; Aranovich et al., 2013; Newton et al., 2014; Manning, Aranovich, 2014). В последние годы также активно развивается представление о достаточно высоком содержании воды при ее низкой активности в гранулитовых флюидах и одновременном сосуществовании одновозрастных рассолов с высокой концентрацией солей и богатых СО2 флюидов, изотопный состав углерода которых указывает на мантийный источник (например, Touret, Huizenga, 2012; Koizumi et al., 2014). Однако возможность инфильтрации богатого Н2О флюида, его фазовый состав, физико-химические характеристики сосуществующих фаз, происхождение Н2О флюидов с высоким содержанием солей, их сосуществование с СО2 флюидами остается предметом активных исследований и полемики (например, Аранович, 2017). Дискуссии касаются причин зарождения расплавов, мигматизации и гранитизации, глубинного метасоматоза или гидротермального рудообразования из этих флюидных потоков в верхних частях коры. Не упоминая все многочисленные исследования флюидов гранулитов умеренных и низких давлений (например, Fonarev et al., 1998; Hisada et al., 2005; Takahashi et al., 2018 и ссылки там), отметим, что данных по флюидам НР гранулитов, формирующихся при нижнекоровом петрогенезисе, значительно меньше (например, Tsunogae et al., 2002, 2008; Santosh et al., 2004; Touret, Huizenga, 2012 и ссылки там; Кориковский, Аранович, 2010, 2015; Higashino et al., 2013; Fu, Touret, 2014).

Предшествующие исследования флюидных включений в породах Лапландского гранулитового пояса (ЛГП) немногочисленны. По данным (Klatt, Schoch, 1974; Barbey, Raith, 1990; Фонарев, Крейлен, 1995), известно, что в российской части северо-западного фрагмента ЛГП НР гранулиты содержат многочисленные включения СО2 и включения N2, в Финской Лапландии отмечены богатые СО2 включения и включения рассолов (5–25 мас. % в экв. NaCl). В основу нашего исследования положены результаты изучения в юго-восточном Кандалакшско-Умбинском фрагменте Лапландского гранулитового пояса Фенноскандинавского щита (юг Кольского полуострова, Порья губа Белого моря, Бушмин, Глебовицкий, 2008) флюидных включений в минералах HP гранулитов и сингранулитовых инфильтрационных метасоматитов, их термобарометрии, расчетов активности воды по минеральным парагенезисам и сравнения этих данных с предсказанием свойств флюидов на основе термодинамических моделей систем H2O–NaCl–CO2 и H2O–CaCl2–CO2. Важно подчеркнуть, что перечисленными методами исследовались одни и те же образцы из названных типов пород. Ранее полученные и опубликованные результаты исследования содержания изотопов кислорода в минералах, использованные при обсуждении результатов настоящего исследования, также были получены при исследовании этих же образцов.

ГЕОЛОГИЧЕСКИЕ СВИДЕТЕЛЬСТВА ВНЕШНИХ ФЛЮИДНЫХ ПОТОКОВ

Палеопротерозойские тектонические покровы преимущественно основных или кислых НР гранулитов Кандалакшско-Умбинского фрагмента ЛГП сложены пакетами пластин, разделенных зонами сдвиговых деформаций c интенсивной бластомилонитизацией (рис. 1). Для этих сдвиговых зон характерна высокая концентрация жил и жилоподобных тел инфильтрационных сингранулитовых метасоматитов разного размера и состава, которые являются свидетельством интенсивных локальных флюидных потоков (Бушмин и др., 2007, 2009) во время гранулитового метаморфизма. Это кварцевые жилы, жилоподобные тела богатых кварцем пород, в том числе кварцевых бластомилонитов, с силлиманитом, ортопироксеном и гранатом (Прияткина, 1977); мафических пород существенно гранатового, ортопироксен-гранатового, ортопироксен-кордиерит-гранатового и кордиерит-ортопироксенового состава (Беляев, 1981); жилы диопсид-скаполитовых, карбонатных и карбонат-диопсидовых пород. На рис. 1 отмечена сдвиговая зона в Порьегубском покрове (преимущественно основные гранулиты) протяженностью более 15 км, в пределах которой были проведены наиболее детальные исследования сингранулитовых метасоматитов. Например, на участке о. Паленый ширина этой зоны около 100 м. На рис. 2 и 3 приведены примеры фотографий жильных тел инфильтрационных сингранулитовых метасоматитов разного минерального состава и примеры детального крупномасштабного картирования сдвиговых зон, в которых они сосредоточены. Эти примеры приведены, прежде всего, потому, что жильные тела отражают пути инфильтрации глубинных флюидных потоков и результаты их взаимодействия с вмещающими породами. Наблюдаемая в сдвиговых зонах отчетливая пространственная дифференциация вещества, выраженная в интенсивном окварцевании и кварцевых жилах, жилоподобных телах мономинеральных, биминеральных и маломинеральных пород, резко обедненных или обогащенных разными породообразующими элементами, свидетельствует о водном составе флюидных потоков, способных растворять, переносить и переотлагать петрогенные элементы.

Рис. 1.

(а) Положение района исследования на схеме тектонического строения Фенноскандинавского щита (Bushmin, Glebovitsky, 2016). ЛГП – Лапландский гранулитовый пояс, КУ – Кандалакшско-Умбинский фрагмент ЛГП. (б) Схема тектонических покровов HP гранулитов Кандалакшско-Умбинского фрагмента ЛГП в районе Порьей губы Белого моря (Бушмин и др., 2007). Колвицкий покров: основные гранулиты (метагаббро-анортозиты); Порьегубский покров: преимущественно основные гранулиты (метаандезито-базальты и метабазальты); Умбинский покров: преимущественно кислые гранулиты (метаосадки). Пунктирные линии: зоны сдвиговых деформаций на границах крупных пакетов тектонических пластин. Пунктирная линия с точками: сдвиговая зона, где были исследованы сингранулитовые НР метасоматиты, в том числе участки с ортопироксен-силлиманитовыми породами (красные звездочки). Точечный пунктир: сдвиговая зона с интенсивной локальной эндербитизацией основных гранулитов – метасоматическая Na–Ca фельдшпатизация и окварцевание, парциальное плавление (мигматиты). Детали см. на рис. 16.

Рис. 2.

Примеры жил разного состава и размера в гранулитах Порьегубского покрова. (а) – Qz жила с Sil, в контактах жилы гранулиты бластомилонитизированы, окварцованы – пропитаны кварцем и содержат сетку кварцевых прожилков, зерна кварца могут быть ленточной формы (о. Высокий); (б) – Qz жила и жилоподобное тело Qz бластомилонитов (о. Наумиха); (в), (г) – Opx-Grt жилы (о. Наумиха); (д) – пример крупномасштабной карты фрагмента сдвиговой зоны (см. рис. 1) с жилоподобными телами и жилами метасоматитов (губа Костариха): 1 – мезократовые Px гранулиты; 2 – меланократовые Px гранулиты; метасоматиты: 3 – богатые Qz породы, Qz бластомилониты с разным содержанием Grt, Pl, Bt ± Spl ± Opx ± Sil; 4 – Opx-Grt жила.

Рис. 3.

Примеры крупномасштабных карт фрагментов сдвиговой зоны, в которой сосредоточены жилы и жилоподобные тела метасоматитов разного состава (см. рис. 1). (а) – о. Наумиха: 1 – Grt-Opx-Cpx и Grt-Opx гранулиты мигматизированные; метасоматиты: 2 – богатые Qz породы с разным содержанием Opx, Sil, Grt, Bt, Pl, Crd ± Spl; 3 – Grt-Opx породы с Bt, Pl, Crd, Qz ± Sil ± Spl. (б) – о. Паленый: 1 – Px гранулиты с редким мелким гранатом, мигматизированные; 2 – Grt-Px гранулиты, мигматизированные; сингранулитовые метасоматиты: 3 – богатые кварцем породы, кварциты, Qz бластомилониты с разным содержанием Sil, Opx, Grt, Crd, Bt, Spl, Spr; 4 – Opx-Grt породы с разным содержанием Sil, Qz, Crd, Spl, Spr; 5 – Di-Scp породы.

ОПИСАНИЕ ПОРОД

Гранулиты

НР гранулиты, в которых развиты сдвиговые зоны с телами изученных инфильтрационных метасоматитов, представляют собой гранат-плагиоклаз-двупироксеновые и гранат-плагиоклаз-ортопироксеновые кристаллические сланцы основного и среднего состава, иногда с тонкими прослоями ортопироксеновых плагиогнейсов. На некоторых участках встречаются прослои Sil-Grt-Bt ± Crd ± Opx22 гнейсов. Гранулиты часто слагают пачки “пластинчатых гранулитов” (рис. 4) и могут быть сильно бластомилонитизированы в контактах пластин. В бластомилонитизированных гранулитах при полной перекристаллизации пироксен может образовывать удлинненные пластинообразные зерна. Мигматизация гранулитов проявлена только в отдельных зонах. Двупироксеновые гранулиты представляют собой среднезернистые породы, состоящие из плагиоклаза (30–60%), ортопироксена (2–10%), клинопироксена (5–25%), граната (5–20%), биотита (5–10%), иногда 5–10% амфибола, иногда кварца (1–5%) и калиевого полевого шпата (до 3%). Кварц слагает редкие зерна или сильно вытянутые линзы. Ортопироксеновые гранулиты представляют собой среднезернистые породы, состоящие из плагиоклаза (30–60%), ортопироксена (10–40%), граната (5–20%), кварца (1–7%), биотита (2–10%), иногда калишпата (до 5%). Типичные минеральные ассоциации этих пород (магнезиальность минералов и содержание анортитового минала в Pl указаны в подстрочных индексах в мольных %): Pl43–81 + Opx52–63 + Cpx44–77 + Grt19–30 + + Bt48–59 ± Hbl44–71 ± Qz ± Kfs и Pl27–91 + Opx54–73 + + Grt20–47 + Bt62–83 ± Qz ± Kfs. Все гранаты содержат кальций: в Cpx-Opx гранулитах ХСа = 0.18–0.20, в Opx гранулитах ХСа = 0.04–0.14.

Рис. 4.

Пример типичных Px гранулитов Порегубского покрова, вмещающих тела сингранулитовых метасоматитов. (а) – пластины Px гранулитов в масштабе обнажения; (б) – образец Grt-Opx гранулита; (в) – шлиф Grt-Opx гранулита, в Pl матрице присутствует кварц.

Сингранулитовые метасоматиты

Кварцевые бластомилониты

Среди богатых кварцем Mg–Al–Si метасоматических пород с силлиманитом, высокомагнезиальными ортопироксеном и гранатом, часто с биотитом, иногда с кордиеритом, калиевым полевым шпатом для исследования были выбраны Qz бластомилониты (pис. 5а), содержащие парагенезис Qz + Opx + Sil ± Grt, который является хорошим минералогическим индикатором высоких давлений (например, Aranovich, Podlesskii, 1989; Подлесский, 2003; Бушмин, Глебовицкий, 2008; Bushmin, Glebovitsky, 2016). Временная связь Qz бластомилонитов с пластическими деформациями в региональных сдвиговых зонах среди НР гранулитов доказана микроструктурными исследованиями (Козлова и др., 1991). В составе Qz бластомилонитов главными минералами, кроме кварца, являются Sil, Opx и Grt, в меньшем количестве может присутствовать Bt, Crd, Kfs, Ab и в виде редких зерен шпинель и сапфирин. Они отличаются тонкополосчатой текстурой и обладают всеми признаками бластомилонитов со “струйчатым” распределением минералов, иногда пластинчатой формой зерен пироксена и граната и отличаются неоднородным отчетливо зональным распределением минералов от полиминеральных до мономинеральных зон (pис. 5б). Метасоматическая зональность (Коржинский, 1982) в Qz бластомилонитах проявлена парагенезисами, в которых число минералов по мере их последовательного замещения и растворения (выщелачивания) закономерно уменьшается вплоть до образования гранат-ортопироксен-силлиманитовых, (Grt: XMg = = 0.57–0.68; Opx: XMg = 0.76–0.81, Al = 0.36–0.42), гранат-силлиманитовых (Grt: XMg = 0.60–0.66), ортопироксен-силлиманитовых (Opx: XMg = 0.77–0.79, Al = 0.37–0.40) и силлиманитовых кварцитов с участками или зонами кварцевого и силлиманитового состава (рис. 5б). Кальция в гранатах практически нет (XCa = 0–0.03, иногда 0.05). В кварцитах с ортопироксеном и гранатом присутствуют зерна шпинели, контактирующие с кварцем.

Ортопироксен-гранатовые и диопсид-скаполитовые жилы

С богатыми кварцем метасоматическим породами ассоциируют в сдвиговых зонах разнообразные по минеральному составу, богатые в разных количествах Fe, Mg, Ca, гигантозернистые, крупнозернистые и среднезернистые порфиробластические и жильные породы: гранатовые, силлиманит-ортопироксен-гранатовые, кордиерит-силлиманит-ортопироксен-гранатовые, а также диопсид-скаполитовые породы (рис. 2, 3, 6).

Выбранные для исследования ортопироксен-гранатовые жилы, как и Qz бластомилониты, содержали парагенезис Opx + Sil + Qz, который являлся минералогическим индикатором высоких давлений. В составе этих жил главными минералами являются гранат (XMg = 0.57–0.71) и ортопироксен (XMg = 0.77–0.82, Al = 0.38–39), может присутствовать Crd (XMg = 0.87–0.93), Bt (XMg = = 0.82–0.94, XAl = 0.32–0.33). Кордиерит может быть как матричным минералом, так и входить в состав минеральных оторочек и симплектитов вокруг граната. Содержание кальция в гранатах, как и в гранатах Qz бластомилонитов, очень низкое (XCa = 0.02–0.05). В подчиненном количестве присутствуют Sil, Qz, в небольшом количестве может присутствовать калиевый полевой шпат, иногда присутствуют шпинель и сапфирин. Изученные диопсид-скаполитовые жилы на участке о. Паленый ассоциируют с телами Qz бластомилонитов и Opx-Grt пород (pис. 3) и по геологическим взаимоотношениям являются относительно более поздними. В составе этих жил главными минералами являются скаполит и диопсид, присутствуют Ttn и Qz. Преимущественно жилы сложены парагенезисом Scp + Di, в котором скаполит более кальциевый (содержание мейонита 77–81 мол. %), а диопсид более магнезиальный (XMg = 0.7). Однако в жилах присутствуют и локальные небольшие участки и прожилки, сложенные более поздним парагенезисом этих же минералов. В них скаполит более натриевый (содержание мейонита 71–74 мол. %) и диопсид более железистый (XMg = 0.4).

ФЛЮИДНЫЕ ВКЛЮЧЕНИЯ В МИНЕРАЛАХ

Метод исследования

Выделены первичные, первично-вторичные и вторичные генетические типы флюидных включений. Первичные включения (p) представлены отдельными, единичными включениями или группами включений, объемно расположенными в зернах минералов. Они образуются во время кристаллизации минералов на стадиях роста. Первично-вторичные (ps) и вторичные (s) включения образуют плоскости или линейные зоны. Зоны первично-вторичных включений не выходят за пределы минерального зерна и возникают, как и первичные, в зонах микродефектов во время кристаллизации минерала или за счет первичных включений при их перераспределении в другие дефекты или микроразрывы (например, Van den Kerkhof, Hein, 2001). Таким образом, они относятся к категории относительно ранних включений, первоначально образовавшихся при кристаллизации минерала. Вторичные включения расположены в плоскостях и зонах пересекающих границы минеральных зерен. Их химический состав может быть аналогичен составу первичных включений в результате, например, миграции включений внутри кристаллов (Зонова и др., 1996; Baumgartner et al., 2014 и ссылки там), но может и отличаться от состава первичных включений в результате преимущественной потери Н2О или СО2 из включений за счет селективной диффузии и дефектов кристалла (Baumgartner et al., 2014; Bakker, Doppler, 2016; Bakker, 2017). В большинстве случаев вторичные включения появляются в связи с совершенно другими процессами позднего флюидного воздействия, и поэтому их состав отличается от составов первичных и первично-вторичных включений (например, Roedder, 1984).

Исследования включений выполнено методом микротермометрии. По температуре плавления (–56.6°С) замороженного включения выделены включения СО2. Однако определяемые температуры плавления углекислотных включений обычно занижены до –58.0°С, что предполагает присутствие азота либо метана в веществе включений (Van den Kerkhof, 1988). Гомогенизация включений СО2 наблюдалась в жидкую фазу. По температуре гомогенизации СО2 (ThCO2) определялась плотность СО2. Иногда присутствующие во включениях СО2 кристаллы твердой фазы идентифицированы рамановской спектроскопией как карбонат (кальцит). Включения азота были выделены по их гомогенизации в газовую фазу вблизи критической точки азота (–147.0°С). У части таких включений гомогенизация в газовую фазу наблюдалась при повышенных температурах (до –120.0°С), что предполагает примесь углекислоты во включениях. В водно-солевых включениях определялись начальные и конечные температуры плавления льда. Начальные температуры плавления льда (Tm1) позволяют определять тип водно-солевой системы (Борисенко, 1977), конечные температуры плавления льда (Tmf) позволяют судить о концентрации солей во включении. Измерения Tm1, главным образом, от −60 до −30°С при понимании, что возможно присутствие и катионов Mg, Fe, K, позволили предполагать систему CaCl2–NaCl–H2O, и концентрация этих солей определялась по тройной диаграмме CaCl2–NaCl–H2O (Crawford, Hollister, 1986; Oakes et al., 1990; Zwart, Touret, 1994; Steele-MacInnis et al., 2011; Chu et al., 2016 и ссылки там). Для включений, имеющих Тmf выше –21°С, концентрация соли на основании конечных температур плавления льда представлена в эквивалентных значениях NaCl (Bodnar, Vityk, 1994). Ряд включений содержит кристаллы галита. По температуре растворения галита определялась концентрация солей в таких включениях (Bodnar, Vityk, 1994).

Включения в гранулитах

Включения СО2

Первичные и первично-вторичные включения СО2 высокой плотности. Высокоплотные p-включения СО2 редки, многочисленные ps-включения СО2 повышенной плотности в кварце, редко в гранате, иногда содержат кристаллики карбоната. TmCO2 около −57…−58°С, пик ThCO2 от −24 до −20°С. Первичные включения очень мелкие, и поэтому не удалось сделать для них количественные измерения, получены данные только для ps‑включений33, см. ESM_1.pdf (Suppl. 1)a. Так как Tm чистого включения СО2 = −56.6°С, то снижение Tm до −58°С, по-видимому, связано с присутствием СН4 или азота. Высокоплотные включения ассоциируют с ps-включениями азота и с включениями рассола с высокой концентрацией соли.

Вторичные включения СО2 низкой плотности. Вторичные включения с ThCO2 от +2 до +28°С в кварце. В некоторых s-включениях присутствует Н2О с низкой концентрацией солей, см. ESM_1.pdf (Suppl. 1). Вода с низким содержанием солей характерна для вторичных водно-солевых включений и ее обнаружение в низкоплотных включениях СО2, скорее всего, связано с поздним перезаполнением.

Включения рассола и азота

Включения рассола редки. Первично-вторичные включения c Tm1 −60°C и Tmf от −18 до −16°C (вероятно, растворы CaCl2 с 21–20 мас. % в экв. NaCl), а также включения с разным соотношением солей CaCl2 и NaCl (Tm1 от −53 до −42°С и Tmf от −36 до −24°C: соленость 25 мас. % CaCl2 + + 5 мас. % NaCl и 7 мас. % CaCl2 + 18 мас. % NaCl соответственно) наблюдаются в кварце между высокоплотными р- и рs-включениями СО2 и многочисленными ps-включениями азота. Вторичные включения рассола низкой солености (Tmf около −0.5°С: сульфатные и карбонатные растворы) обычно встречаются в ассоциации с s-включениями азота низкой плотности, гомогенизирующими в газовую фазу.

Рис. 5.

Сингранулитовые метасоматиты. (а) – пример обнажения с жилоподобным телом Qz бластомилонита на контакте с вмещающим Px гранулитом (обнажение БЛГ, о. Паленый, см. рис. 3б). (б) – зональный образец Qz бластомилонита с отчетливо выраженной минеральной метасоматической зональностью из этого обнажения.

Рис. 6.

Сингранулитовые метасоматиты: мафические жильные и жилоподобные породы, ассоциирующие с богатыми кварцем породами в сдвиговых зонах. (а) – жила крупнозернистой Opx-Grt породы на контакте с Qz бластомилонитом (о. Паленый), S-форма крупного порфиробласта граната с S-образным распределением минеральных включений свидетельствует о синкинематических условиях кристаллизации в сдвиговой зоне; (б), (в) – гигантские гранаты (диаметр до 12–15 см) в жиле крупнозернистой Opx-Grt породы (о. Паленый); (г) – жильная крупнозернистая Crd-Opx-Grt порода (о. Костариха); (д) – Opx-Grt жила на контакте с Qz бластомилонитом (о. Паленый); (е) – Opx-Grt жилы в Sil кварците (о. Паленый).

Включения в сингранулитовых метасоматических породах

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

Рис. 7.

Сингранулитовые метасоматиты (Qz бластомилонит): ассоциации первичных и первично-вторичных включений СО2, рассола и N2 в матричном кварце из зоны кварцита с гранатом и силлиманитом (см. рис. 5б). (а) – ps-включения СО2 (L) высокой плотности, рядом группа включений рассола (G + L), обр. Л4-1в; (б) – ps-включения СО2 с высокой (L) и низкой (L + G) плотностью, обр. Л4-1в; (в) – включение СО2 (L+G) и включения рассола (L+G), обр. Л4-1в; (г) – ps-включения СО2 (L) высокой плотности и включения рассола со слабо анизотропной фазой, обр. Л4-1б; (д) – ассоциация ps-включений СО2 высокой плотности, рассола и N2, обр. БЛГ-г; (е) – включения рассола и крупное включение рассола с галитом среди включений азота, обр. БЛГ-1а. L – жидкая фаза, G – газовая фаза, Crb – карбонат.

Рис. 8.

Сингранулитовые метасоматиты (Qz бластомилонит): первично-вторичные флюидные включения в зоне кварцита с гранатом, ортопироксеном и силлиманитом (см. рис. 5б), обр. БЛГ-1б. (а), (б) – ассоциация включений СО2 высокой плотности, рассола и N2 в зерне кварца, заключенном в гранат; (в), (г) – ассоциация включений СО2 высокой плотности, рассола и N2 (часто с Н2О) в зерне кварца, заключенном в ортопироксен. L – жидкая фаза.

Рис. 9.

Сингранулитовые метасоматиты (Opx-Grt жила): включения СО2 высокой плотности, рассола и N2. (а) – p-включение СО2 с формой отрицательного кристалла в кварце, обр. Л4-2б; (б) – ps-включения СО2 с формой отрицательного кристалла в кварце, обр. Л4-2б; (в) – ассоциация ps-включений рассола и ps-включений N2 в кварце, многие включения N2 содержат Н2О, обр. Л4-2б; (г) – ps-включения рассола низкой солености в кварце, обр. Л4-2е; (д) – ассоциация ps-включений рассола высокой солености и N2 в кварце, обр. Л4-2е; (е) – червеобразные ps-включения СО2 в крупном порфиробласте граната, обычные для включений такого состава в гранате, обр. Л4-2г. L – жидкая фаза, G – газовая фаза, Crb – карбонат.

Рис. 10.

Сингранулитовые метасоматиты (Opx-Grt жила): (а) – ассоциация ps-включений СО2 высокой плотности и ps-включений N2 в кварце, заключенном в крупное зерно граната, обр. Л4-3; (б) – увеличенный участок с включениями N2. L – жидкая фаза, G – газовая фаза.

Рис. 11.

Сингранулитовые метасоматиты (диопсид-скаполитовая жила): включения СО2 высокой плотности в скаполите, обр. Б880-22б. (а) – p-включения СО2, некоторые содержат карбонат; (б) – ps-включения СО2, с карбонатом и без него, в ассоциации с включениями рассола. L – жидкая фаза, G – газовая фаза, Crb – карбонат.

Включения в кварцевых бластомилонитах

Включения СО2

Первичные и первично-вторичные включения СО2 высокой плотности. Среди p- и ps-включений СО2 (наблюдались только в кварце) высокой плотности (Th до −32°С) преобладают ps-включения, а p-включения (Th –18…–20°C) единичны, см. ESM_2.pdf (а) (Suppl. 2)b. Иногда присутствуют кристаллы карбоната. Зоны таких ps-включений обычно близки к направлениям блаcтомилонитовой полосчатости.

Первичные, первично-вторичные и вторичные включения СО2 низкой плотности. В целом, по сравнению с включениями высокой плотности, включения низкой плотности заметно преобладают. Группа низкоплотных p- и ps-включений CO2 (Th до +32°С с пиком +22°С) объединяет, см. ESM_2.pdf (b) (Suppl. 2) многочисленные ps-включения в матричном кварце и в кварце, включенном в гранат, десятки ps-включений в кварце, включенном в ортопироксен и редкие ps-включения в гранате. Первичные включения низкой плотности единичны и встречены только в кварце, заключенном в ортопироксен. Группа низкоплотных s‑включений (Th от +6 до +32°С с пиком +22°С) представлена включениями в кварце и только редкие из них заключены в кварце внутри граната, см. ESM_2.pdf (с) (Suppl. 2). Большинство низкоплотных включений СО2 распределяется не только в зонах ps-включений как вдоль, так и перпендикулярно бластомилонитовой полосчатости, но часто присутствуют в зонах s-включений, расположенных только ортогонально бластомилонитовой полосчатости.

Включения рассола

Исследовано около 100 водно-солевых включений (CaCl2 и NaCl, Tm1 до −50°С), см. ESM_3.pdf (а, b) (Suppl. 3)c. Первичные и ps-включения присутствуют в матричном кварце и в кварце, заключенном в гранат и ортопироксен. Наряду с включениями повышенной солености (p- и ps-включения с соленостью в диапазоне 10–23 мас. % в экв. NaCl) и включениями высокой солености (p-включения в кварце с 20 мас. % CaCl2 + 10 мас. % NaCl, ps-включения в кварце с 5 мас. % CaCl2 + 20 мас. % NaCl и с 5–23 мас. % в экв. NaCl), встречены редкие включения с галитом. Эти первичные включения рассола с галитом находятся не только в матричном кварце, но и в кварце, заключенном в зерна граната и ортопироксена. Включения с галитом в кварце, заключенном в гранат, находятся среди азотных включений. В азотных включениях довольно часто вокруг газового пузыря азота присутствует водная фаза, но обычен также только газовый пузырь азота без водной фазы. В водно-азотных включениях водная фаза имеет низкую соленость, которая в целом типична для этого типа включений. В этих включениях с галитом растворение галита наблюдалось при Т = 170°С, что соответствует солености 30.5 мас. %. Включения с галитом в кварце, заключенном в ортопироксен, найдены среди ps-включений СО2 и водно-солевых ps-включений с низкой соленостью. Во включениях этого типа растворение галита наблюдалось при Т = 260°С, что соответствует солености 35 мас. %. Включения с галитом в матричном кварце найдены среди азотных и СО2 включений. Во включениях этого типа растворение галита наблюдалось при Т = 175°С, что соответствует солености 30.5 мас. %. Исследованные включения с галитом были мелкими, наблюдать в них удалось только окончание растворения галита и по этой температуре получить информацию о концентрации NaCl (30–35 мас. %). Случай получения начальной температуры плавления льда около –50°С указывает на присутствие CaCl2. Если принять во внимание обычное присутствие двух солей CaCl2 и NaCl во включениях рассола, то общая соленость во включениях с галитом может превышать 30–35 мас. %.

Включения солей кальция и натрия, вероятнее всего, отражают солевой состав флюида, так как хлорид магния достоверно обнаружен только во включениях в магнезиальном гранате, а соли кальция и натрия во включениях в минералах, не содержащих кальций и натрий. Во включениях с Tm1 = −31°С может быть чистый хлорид магния, а может быть и смесь натриево-кальциевых растворов. Включения с Tm ниже −31°С, вероятнее всего, содержат смеси солей. Такая неопределенность в составе солей часто позволяет говорить о величинах солености только в эквиваленте NaCl.

Включения азота

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

Первично-вторичные включения азота. Преобладающая часть ps-включений азота.  см. ESM_3.pdf (с) (Suppl. 3) присутствует в матричном кварце, реже они встречается в кварце, заключенном в гранат и ортопироксен. Азотные включения низкоплотные, с Th в газовую фазу часто выше −147°С (критическая точка азота), что, видимо, связано с присутствием СО2 во включениях. Также в них часто видна Н2О и они находятся в ассоциации с первично-вторичными включениями рассола низкой солености до 3−4 мас. % в экв. NaCl.

Вторичные включения азота. Низкоплотные ps-включения азота, см. ESM_3.pdf (d) (Suppl. 3), встреченные только в матричном кварце, также как и низкоплотные ps- и s-включения СО2 распределены по линейным зонам. Как и в ps-включениях, в s-включениях азота также можно предполагать наличие СО2, часто присутствует Н2О и они находятся в ассоциации с включениями рассола низкой солености до 3–4 мас. % в экв. NaCl.

Включения в ортопироксен-гранатовых жильных породах

Включения СО2

Так же как и в кварцевых бластомилонитах, среди p- и ps-включений СО2 в минералах Opx-Grt жил преобладают ps-включения, многочисленны s-включения, а p-включения единичны. Они представлены группами включений высокой и низкой плотности, но преобладают включения низкой плотности.

Первичные и первично-вторичные включения СО2 высокой плотности. Высокоплотные p- и ps-включения СО2 (Th до −22°С) встречены только в матричном кварце между преобладающими в породе зернами ортопироксена и граната, см. ESM_4.pdf (а) (Suppl. 4)d.

Первично-вторичные и вторичные включения СО2 низкой плотности. Так же как и в кварцевых бластомилонитах, многочисленные низкоплотные ps-включения (Th до +32°С, пик +12°С) найдены как в матричном кварце, так и в главном матричном минерале – гранате или кварце, заключенном в гранат, см. ESM_4.pdf (а) (Suppl. 4). Также многочисленные низкоплотные s-включения (Th до +30°С, пик +26°С) присутствуют в кварце и реже встречаются в кварце, заключенном в гранат, см. ESM_4.pdf (b) (Suppl. 4).

Включения рассола

В Opx-Grt жилах обнаружены редкие ps-включения рассола в матричном кварце со смесью солей CaCl2 и NaCl (Tm1 = −42 и −26°C) и соленостью 21–10 мас. % в экв. NaCl. Как и в кварцевых бластомилонитах, найдены редкие включения с галитом среди ps-азотных включений в кварце, заключенном в гранат. Эти включения декрипитировали при Т ~ 220°С, значит температура растворения галита больше 220°С и соленость рассола более 33 мас. % NaCl. Вероятно, в них присутствует и CaCl2. Также в матричном кварце присутствуют поздние s-включения низкой солености (до 4 мас. % в экв. NaCl), часто встречающиеся в ассоциации с включениями азота.

Включения азота

Во многих образцах Opx-Grt жил немногочисленные p- и ps-включения азота присутствуют как в матричном кварце, так и в кварце, заключенном в ортопироксен и гранат, см. ESM_4.pdf (с) (Suppl. 4). Так же как и в кварцевых бластомилонитах, по-видимому, часть включений содержит примесь СО2, часто вокруг газового пузыря азота видна Н2О и они ассоциируют с водными низкосолеными включениями (до 3–4 мас. % в экв. NaCl). Некоторые включения азота не содержат Н2О.

Включения в диопсид-скаполитовых жильных породах

Включения СО2

По геологическим данным диопсид-скаполитовые жилы более поздние по отношению к кварцевым и ортопироксен-гранатовым породам. Включения найдены только в скаполите. Первичные и первично-вторичные включения СО2 в скаполите небольшой плотности (Th от −4 до +20°С), близкой к низкоплотным вторичным включениям ранних метасоматитов, см. ESM_5.pdf (Suppl. 5)e. Поэтому разуплотнение включений ранних метасоматитов, видимо, связано по времени с деформациями во время образования этих жильных пород.

Включения рассола и азота

Включения рассола и азота в диопсид-скаполитовых жилах изучить не удалось. Были обнаружены только очень мелкие включения, в которых не было видно фазовых изменений, но которые по форме, внешнему виду аналогичны включениям Н2О-соль и азота, широко проявленным в других типах изученных пород.

ТЕРМОБАРОМЕТРИЯ

Мультиравновесная термобарометрия

Термобарометрия пород выполнена методом TWEEQU (Berman, 1991), который позволяет не только рассчитать искомые величины Р и Т, но и проверить равновесность использованных при этом составов минералов. Для расчетов применена взаимосогласованная база данных Р.Дж. Бермана и Л.Я. Арановича (Berman, Aranovich, 1996) версии 2.02b, которая наиболее подходит для гранулитовых ассоциаций. Результаты мультиравновесной термобарометрии гранулитов и сингранулитовых метасоматитов района исследования опубликованы в работах (Бушмин и др., 2007; Лебедева и др., 2012). Для тех образцов, в которых изучались флюидные включения, получены следующие интервалы вариаций Р-Т параметров: для гранулитов ~8.7–11 кбар и ~800–900°С, для сингранулитовых метасоматических кварцевых бластомилонитов ~11–9 кбар и ~920–850°С, для жил мафических ортопироксен-гранатовых пород ~11–8 кбар и ~925–850°С. Примеры результатов расчета Р-Т параметров методом мультиравновесной термобарометрии показаны на рис. 12, 13. Пиковые Р-Т условиях кварцевых и мафических метасоматитов подтверждены анализом парагенезисов минералов – во всех изученных породах сохранялся неизмененный НР парагенезис Opx + Sil + Grt + Qz.

Рис. 12.

Примеры расчета методом TWEEQU Р-Т параметров гранулитов: (а) – Px гранулит (Pl, Grt, Opx, Qz, Bt), обр. Б1017-9-1, о. Костариха; (б), (в) – Px гранулит (Grt, Opx, Pl, Bt, Qz), обр. Б856-4_F2 и Б856-4_F1, о. Паленый. IR – число независимых реакций, СКО – среднеквадратичное отклонение.

Рис. 13.

Примеры расчета методом TWEEQU Р-Т параметров сингранулитовых метасоматитов: (а), (б) – кварцевые бластомилониты (Qz, Opx, Sil, Grt, Bt), обр. БЛГ-2_F2 и БЛГ-2_F3, о. Паленый; (в), (г) – мафические жилы: (в) – Opx-Grt порода с Bt, Sil, Pl, Qz, обр. Б801-56, о. Паленый и (г) – Crd-Opx порода с Pl, Qz, Sil, Bt, Grt, обр. Б1021-13-1, о. Костариха. IR – число независимых реакций, СКО – среднеквадратичное отклонение.

Изохоры флюидных систем

Изохоры построены на основании данных по температурам гомогенизации СО2 в первичных и первично-вторичных включениях, см. ESM_6.pdf (Suppl. 6)f. Величины давления оценены по наиболее плотным включениям. Для вмещающих гранулитов Т = 800–900°С соответствуют Р ~ 7–7.5 кбар, см. ESM_6.pdf (Suppl. 6). Для метасоматических кварцевых бластомилонитов Т = 850–920°С соответствуют Р ~ 8–8.5 кбар, см. ESM_6.pdf (Suppl. 6). Эти более низкие оценки давления с использованием изохор близки к таковым в гранулитах (~8–9 кбар) и к величинам давления начала ретроградного минералообразования в метасоматитах (~9 кбар), полученным методами мультиравновесной термобарометрии. Такая разница в величинах (1–1.5 кбар) флюидного давления (включения) и метаморфического давления (термобарометрия) систематически отмечается во многих гранулитах (например, Touret, 2001). Все остальные первично-вторичные и вторичные включения, обладающие пониженной и низкой плотностью, связываются нами с поздними ретроградными стадиями – с процессами разуплотнения при вскрытии включений во время их пространственного перераспределения в зернах минералов во вторичные цепочки поздних структурных дефектов. А также с поздними гидротермальными процессами в гранулитах. К последним относятся проявленные в районе исследования низкотемпературные процессы образования рудных карбонатных жил на месторождениях серебра и гидротермальные процессы, сопровождавшие внедрение даек лампроитов при эксгумации гранулитов.

РАСЧЕТ АКТИВНОСТИ Н2О (аН2О) ПО МИНЕРАЛЬНЫМ ПАРАГЕНЕЗИСАМ

Расчеты величин активности воды по равновесным минеральным парагенезисам с биотитом и калиевым полевым шпатом в гранулитах и сингранулитовых метасоматитах выполнены методом TWEEQU, описание которого приводилось выше в разделе “Термобарометрия” и который позволяет также оценить активность воды во флюиде. И в гранулитах, и в метасоматитах встречены только единичные парагенезисы, подходящие для расчета аН2О. Примеры результатов расчета приведены на рис. 14.

Рис. 14.

Пример расчета активности Н2О в гранулитах и сингранулитовых метасоматитах на участке о. Костариха: (а) – Px гранулит (Pl, Grt, Opx, Qz, Bt, Kfs), обр. Б1017-9-1; (б) – метасоматическая Crd-Opx жильная порода (Opx, Crd, Qz, Sil, Bt, Kfs, Pl, Grt), обр. Б1021-13-1. IR – число независимых реакций, СКО – среднеквадратичное отклонение.

Рассчитанные величины активности воды при формировании гранулитов (аН2О ~ 0.49), вмещающих зоны сингранулитовых инфильтрационных метасоматитов, и самих метасоматитов (аН2О ~ 0.51) оказались близкими и достаточно большими. Полученные значения аН2О согласуются с устойчивостью магнезиального и высокотитанистого глиноземистого биотита в гранулитах и его широкой устойчивостью в метасоматических породах. Магнезиальные биотиты с небольшим содержанием Cl могут сосуществовать с флюидом, богатым Cl при Р-Т условиях гранулитовой фации (Aranovich, Safonov, 2018). В ряде случаев исследованные биотиты содержат Cl, который наряду с Ti, также стабилизирует биотит при высоких температурах (например, Hammerli, Rubenach, 2018). Устойчивость биотита, содержащего Cl в НР гранулитах, отмечалась в ряде исследований (например, Higashino et al., 2013) и связывалась с инфильтрацией богатого Cl флюида в локальных сдвиговых зонах при Р-Т параметрах, близких к пиковым (850°С и 11 кбар). Обращает на себя внимание, что вычисленные значения активности воды во флюиде близки к оценкам (0.38–0.40), полученным ранее в районе, где проведены наши исследования (Кориковский, Аранович, 2015), но оказались значительно выше оценок (0.10–0.35), полученных для минеральных ассоциаций гранулитов в северо-западном фрагменте ЛГП (Фонарев, Крейлен, 1995). Скорее всего, более высокие оценки аН2О в работах последних лет связаны с ревизией Р-Т условий реакций дегидратации биотита, проведенных в (Aranovich, Newton, 1998; Berman et al., 2007). В то же время для гранулитов других районов мира имеются данные как о низких, так и о высоких значениях активности воды (Newton et al., 2014 и ссылки там). Например, для пиковых Р-Т условий гранулитового комплекса Лимпопо (~900°С, 9 кбар) были получены величины аН2О = = 0.40–0.55 (Hisada et al., 2005; Koizumi et al., 2014). Полученные нами оценки согласуются также с оценками аН2О = 0.4–0.6 при образовании аналогичных ортопироксен-гранат-калишпат-кварцевых и калишпат-ортопироксен-силлиманит-гранат-кварцевых метасоматических пород среди НР гранулитов Кольской гранулито-гнейсовой области в районе пос. Кица (Бушмин, 1996; Доливо-Добровольский, 2002).

ТЕРМОДИНАМИЧЕСКИЕ МОДЕЛИ: ФАЗОВЫЙ СОСТАВ И ХАРАКТЕРИСТИКИ ФЛЮИДА

Мы выполнили термодинамическое предсказание возможного фазового состояния гранулитовой флюидной системы на основе моделей тройной системы Н2О–СО2–NaCl (Аранович и др., 2010б) и Н2О–СО2–CaCl2 (Ivanov, Bushmin, 2018; Иванов, Бушмин, 2019), рассчитав положение сольвуса при установленных Р-Т параметрах изученных НР гранулитов и сингранулитовых метасоматитов (рис. 15).

Рис. 15.

Фрагменты фазовых диаграмм флюидных систем H2O–CO2–NaCl (а) и H2O–CO2–CaCl2 (б) при Р-Т параметрах НР гранулитовой фации. Толстые линии – линии сольвуса, разделяющие области гомогенного и двухфазного флюида, конноды сосуществующих флюидных фаз – тонкие линии (аН2О = 0.51) и пунктир (аН2О = 0.40).

В поле двухфазного флюида возможно провести конноды, принимая во внимание состав флюидных включений, содержания солей (CaCl2 и NaCl) во флюидных включениях рассолов и величины аН2О, полученные нами и известные по литературным данным. Учитывая содержания солей во включениях рассолов (до ~30–35 мас. % соли) и значение аН2О ~ 0.51, рассчитанное по минеральным равновесиям, проведены конноды (рис. 15), отражающие теоретические составы сосуществующих несмесимых флюидных фаз. В системе NaCl–H2O–CO2 (рис. 15а) нижний конец конноды соответствует СО2–Н2О флюидной фазе с малым количеством соли (8 мас. %) и близкими величинами мольных долей Н2О (~0.5) и СО2 (0.4). Верхний конец конноды соответствует рассолу (30 мас. % соли) с низким содержанием СО2 (x ~ 0.2).

Поле двухфазного флюида в системе CaCl2–CO2–H2O заметно шире и конноды с аН2О = 0.51 могут быть проведены для всех трех представленных на рис. 15б линий сольвуса. Нижние концы коннод соответствует Н2О–СО2 флюидной фазе также с малым количеством соли (5–12 мас. %) и близкими величинами мольных долей Н2О (~0.44–0.55) и СО2 (~0.42–0.55). Верхние концы коннод соответствуют концентрированному рассолу (34–50 мас. % соли) и крайне низкому содержанию СО2 (x ~ 0.05–0.21). Возможные причины отклонения модельного состава флюидных фаз от состава наблюдавшихся чистых включений рассола и чистых включений СО2 будут рассмотрены ниже. По-видимому, именно углекислотно-водный флюид низкой солености (5–12 мас. % соли) участвовал в формировании зон изученных сингранулитовых метасоматитов (см. рис. 3). Литературные данные (см. выше) позволяют также провести конноды с более низкой величиной аН2О = = 0.40 (рис. 15а, 15б). Нижние концы этих коннод, в особенности для системы с CaCl2, отвечают флюидной фазе с массовой долей СО2 до 90% и крайне низкой соленостью, что обращает на себя внимание в связи с типичной картиной преобладания чистых включений СО2 в гранулитах.

На примере системы с CaCl2 рассчитаны плотности возникающих и участвующих в НР гранулитовом петрогенезисе несмесимых флюидных фаз, ограниченных коннодами на рис. 15а: ~1.3–1.7 г/см3 для рассола и ~1.1 г/см3 для богатого СО2 флюида. Величины этих плотностей соответствуют жидкому состоянию флюидных фаз. Водная флюидная фаза, сосуществующая c менее плотным богатым СО2 флюидом, имеет чрезвычайно высокую концентрацию соли и, следовательно, исключительную способность к фильтрации по сравнению с СО2-содержащей флюидной фазой (например, Watson, Brenan, 1987; Holness, 1992; Аранович, 2017) как в объеме межзернового пространства, так и в локальных зонах сдвиговых деформаций. Это позволяет предполагать (см. ниже) участие подобного концентрированного рассола в процессах интенсивной метасоматической Na-Ca фельдшпатизации в локальных зонах гранитизации основных гранулитов (рис. 1, 16). Участие рассолов в процессах гранитизации многих гранулитовых комплексов ранее было убедительно доказано Л.Л. Перчуком и его коллегами (Perchuk, Gerya, 1993; Перчук и др., 1994; Perchuk et al., 2000; Кориковский, Аранович, 2015; Аранович, 2017).

Рис. 16.

Зоны эндербитизации (светло-серый цвет) в основных Px гранулитах с метасоматическими явлениями окварцевания, Na–Ca фельдшпатизации (олигоклаз, кислый андезин) и парциальным плавлением более легкоплавких минеральных ассоциаций, инициированные инфильтрацией рассолов (Столбовые Луды).

ОБСУЖДЕНИЕ РЕЗУЛЬТАТОВ

Геологическая документация флюидных потоков

Мы привели примеры жильных тел инфильтрационных сингранулитовых НТ-НР метасоматитов разного минерального состава и примеры детального картирования фрагментов сдвиговых зон, в которых они сосредоточены. Прежде всего, потому, что они отражают пути инфильтрации глубинных флюидных потоков и результаты их взаимодействия с вмещающими породами, но, по нашему предположению, до сих пор задокументированы недостаточно, хотя чрезвычайно важны для понимания роли флюидов в петрогенезисе нижней коры. Присутствие в гранулитах, по крайней мере, двухфазного несмесимого метасоматизирующего флюида подтверждается этими многочисленными геологическими наблюдениями. Свидетельством инфильтрации несмесимых одновозрастных флюидов являются локализованные в одних и тех же или в соседних зонах сдвиговых деформаций тела сингранулитовых инфильтрационных метасоматитов: кварцевых бластомилонитов, разнообразных по минеральному составу кварцевых жил и богатых кварцем пород с Sil, Grt, Opx, Crd, жил и жилоподобных тел мафических гранат-ортопироксеновых пород с Sil, Crd и Qz, диопсид-скаполитовых пород, карбонат-диопсидовых и карбонатных пород мощностью от первых самтиметров до первых метров.

К настоящему времени известно больше полевых подтверждений флюидной инфильтрации в гранулитах низких и умеренных давлений как во всем объеме пород во время метаморфизма, так и в виде локальных флюидных потоков, вызывавших метасоматоз (например, Коржинский, 1947; Судовиков, 1956). Для НР гранулитов такие примеры единичны (например, Бушмин, 1996; Доливо-Добровольский, 2002, Ходаревская, Кориковский, 2007; Бушмин и др., 2007, 2009).

Три главных типа включений

В этом разделе при обобщении полученных данных и обсуждении состава флюидных включений мы обращаем внимание преимущественно на результаты, полученные по первичным и первично-вторичным включениям, так как предполагаем, что их химический состав, в отличие от состава вторичных включений, соответствует или близок к составу флюида при кристаллизации минералов гранулитов и сингранулитовых метасоматитов при установленных для них пиковых и послепиковых Р-Т параметрах. Первично-вторичные включения сохраняют такую же информацию о составе минералообразующего флюида, что и первичные включения и, прежде всего, это касается включений в кварце (например, Crawford, Hollister, 1986).

В исследованных НР гранулитах (~8.7–11 кбар, ~800–900°С) и сингранулитовых инфильтрационных НР метасоматитах (жилоподобные тела кварцевых бластомилонитов и жилы ортопироксен-гранатовых пород с Р-Т параметрами 11–9 кбар, ~920–850°С, диопсид-скаполитовые жилы) ЛГП присутствуют одинаковые ассоциации сингенетичных типов включений с флюидными фазами контрастного химического состава: СО2 (заметно преобладают), рассолы (главные соли CaCl2 и NaCl) и N2 ± H2O. Эти три типа включений сосуществуют в одних и тех же генерациях ранних включений: более редких первичных (p) и преобладающих первично-вторичных (ps). Включения СО2 – высокой и низкой плотности, включения N2 – низкой плотности. В гранулитах присутствуют ps-включения рассолов богатых CaCl2 (до 21–20 мас. % в экв. NaCl) и включения концентрированных рассолов с разным соотношением солей CaCl2 и NaCl: соленость 25 мас. % CaCl2 + 5 мас. % NaCl и 7 мас. % CaCl2 + 18 мас. %. NaCl. В метасоматитах p- и ps-включения рассолов обнаруживают широкий диапазон общего содержания солей (до 30–35 мас. %) и разное соотношение концентраций главных солей: p-включения с соленостью 20 мас. % CaCl2 + 10 мас. % NaCl; ps-включения с соленостью 5 мас. % CaCl2 + 20 мас. % NaCl; p- и ps-включения с соленостью 5–23 мас. % в экв. NaCl; p-включения с галитом (до 35 мас. % NaCl). В целом CaCl2 является преобладающим компонентом среди солей в группе ранних p- и ps-включений изученных пород, что не является редкостью для гранулитов Фенноскандинавского щита (Touret, 1985; Fonarev et al., 1998). Все эти особенности включений позволяют предполагать внешний флюидный поток и свидетельствуют о том, что на пике Р-Т параметров захват включений происходил из гетерогенного флюида, в котором одновременно сосуществовали несмесимые водно-солевой и богатый СО2 флюиды, изначально содержащие азот. При одинаковом химическом типе включений и в гранулитах, и в сингранулитовых метасоматитах общее количество присутствующих флюидных включений в метасоматитах значительно выше, что указывает на более высокое отношение флюид/порода в локальных флюидопроводящих сдвиговых зонах.

Другие районы ЛГП с НР гранулитами, где были изучены флюидные включения, немногочисленны. Такие же типы включений присутствуют в НР гранулитах российской части северо-западного фрагмента ЛГП (Фонарев, Крейлен, 1995). В этом районе ранние типы включений в гранулитах представлены высокоплотными включениями, богатыми СО2, с примесью или заметным количеством CH4 или N2 и включениями азота с примесью СН4. Причем, азотные включения присутствуют как одиночные, так и в виде скоплений. В Финнской Лапландии (Klatt, Schoch, 1974; Barbey, Raith, 1990) отмечены богатые СО2 включения и включения рассолов (5–25 мас. % в экв. NaCl). Другие примеры относятся к архейской Центрально-Кольской гранулитовой области Фенноскандинавского щита, где исследованы зоны НР гранулитового метаморфизма и метасоматоза с пиковыми параметрами 9–10 кбар, 800–900°С (Доливо-Добровольский, 2002). Эти гранулиты содержат такие же главные типы сингенетичных включений (или как их называют авторы – “группы синхронных включений”): включения богатые СО2 с 5–35 мол. % N2, включения рассолов CaCl2 (до 28 мас. % в экв. NaCl) или NaCl (до 30 мас. %), чистые включения N2, иногда с 8 мол. % CH4 (Fonarev et al., 1998.). Позднее в этом районе были детально исследованы подобные флюидные включения в сингранулитовых сапфириновых породах метасоматического происхождения в зонах НР гранулитового метаморфизма (неопубликованные данные Е.А. Вапника). Многочисленные находки ассоциации сингенетичных включений, богатых СО2 (больше 50 мол. %), и концентрированных рассолов (до 50% и более мас. % в экв. NaCl) известны давно и в других регионах, где распространены метаморфические породы с высокими Р-Т параметрами (например, Sisson et al., 1981; Trommsdorff et al., 1985). Все три типа этих включений обычны для гранулитов, а устойчивая ассоциация СО2 включений и NaCl–CaCl2 рассолов с N2 включениями типична особенно для эклогитов и НР гранулитов (например, Touret, 1985, 2001 и ссылки там).

Захват включений из гетерогенного флюида

Анализ фазового состояния гранулитовой флюидной системы на основе моделей тройных систем Н2О–СО2–NaCl (Аранович и др., 2010б) и Н2О–СО2–CaCl2 (Ivanov, Bushmin, 2018; Иванов, Бушмин, 2019) при установленных Р-Т параметрах изученных НР гранулитов показал, что флюид мог находиться в гетерогенном состоянии. Сопоставление результатов термодинамического предсказания с данными по флюидным включениям убедило нас в том, что захват включений происходил из гетерогенного флюида, в котором одновременно сосуществовали несмесимые водно-солевой и богатый СО2 флюиды, а, обращая внимание на данные по включениям азота, предполагаем, что богатая СО2 флюидная фаза изначально содержала и N2.

При захвате из гетерогенного флюида вероятно присутствие включений разных типов в одном и том же кластере или микротрещине кристаллизующегося минерала. Это вероятно в том случае, если петрогенезис протекал сразу после фазовой сепарации гомогенного флюида и несмесимые флюиды не успели разделиться в геологическом пространстве в силу многих причин (деформация, стресс, смачиваемость и повышенная способность к фильтрации, гравитационная сепарация и др.). Часто включения несмешивающихся флюидов встречаются в соседних микротрещинах. В отношении СО2 включений и рассолов это нашло объяснение в экспериментальных исследованиях синтетических включений, в результате которых было показано, что возможен и преимущественный захват только одного из сосуществующих флюидов (Sterner, Bodnar, 1984). Присутствие в гранулитах среди преобладающих включений СО2 высокой плотности одновозрастных включений рассолов, в том числе с CaCl2 и NaCl, известно давно (например, Newton, 1986; Touret, 1985, 1995), хотя в ряде других исследований они не обнаруживались или не принимались во внимание (например, Фонарев, Крейлен, 1995). Более того, только редкие включения рассолов и множество газовых включений отмечалось даже в зонах метасоматической переработки, где очевидно участие водных флюидов (например, Damman et al., 1996). Также известны и примеры формирования гранулитов в результате притока концентрированных водных рассолов (Nijland et al., 1998). В связи с этим при всестороннем исследовании роли рассолов высокой концентрации в нижнекоровом гранулитовом метаморфизме была принята во внимание очевидная недооценка данных по присутствию немногочисленных, по сравнению с СО2, включений рассолов (Newton et al., 1998).

То, что концентрированные водные рассолы с небольшим количеством СО2 могут устойчиво сосуществовать с флюидами, богатыми СО2, при LP гранулитовом метаморфизме (до 5–7 кбар) давно показано на примере системы H2O–NaCl–CO2 (Bowers, Helgeson, 1983; Duan et al., 1995). Хотя в основу моделей и расчетов были положены фрагментарные экспериментальные данные, эти представления послужили основанием предполагать одновозрастность ассоциаций включений рассолов и СО2 и в более широком диапазоне Р-Т параметров глубинного петрогенезиса (Crawford, Hollister, 1986; Touret, 1985), что и было позднее подтверждено также на примере системы H2O–NaCl–CO2 (Аранович и др., 2010б).

Что же касается включений с азотом, то здесь нам многое неясно. Петрологическая проблема, связанная с источником азота во флюидных включениях, известна давно. Существует много сообщений о присутствии в гранулитах включений азота (например, Touret, Dietvorst, 1983). При этом обращалось внимание на совместное присутствие чистых включений N2 и чистых включений СО2 в одной и той же зоне (микротрещине). На основании того, что эти летучие неполярные газы должны хорошо смешиваться между собой в широком диапазоне температур, предполагалось, что они могли поступать в гранулиты в виде отдельных струй флюидов на поздних стадиях. Отсутствие поля их несмесимости в системе CO2–N2 при рассматриваемых Р-Т условиях (Kooi et al., 1998) позволяло интерпретировать сосуществование включений СО2 и N2 в соседних зонах (микротрещинах) за счет сближенной по времени инфильтрации таких разных по составу флюидов (например, Fonarev et al., 1998). Однако в наших исследованиях установлено присутствие N2 совместно с ранними включениями СО2 и рассола в одних и тех же зонах (микротрещинах), что вроде бы необъяснимо поздней инфильтрацией азота. Конечно, можно предположить, что включения азота связаны с остатками поровых флюидов осадочного протолита. Однако осадочного протолита для исследованных гранулитов и метасоматитов нет. Протолитом гранулитов являлись андезибазальты и базальты (Козлов, 1990), протолитом сингранулитовых метасоматитов – вышеназванные Pl-Рx гранулиты (Бушмин и др., 2007). Несмесимость в системе H2O–N2 проявляется при температурах ниже 400°С при низких давлениях. Расслоение смеси Н2О–N2 возможно при более высоких температурах только в сочетании с высоким давлением, например, в мантийных условиях зон субдукции (Аранович, 2013 и ссылки там). При Р ~ 10 кбар расслоение можно ожидать при Т не выше ~500°С. Однако солевая нагрузка приводит к расширению этой области расслоения. Также можно принимать во внимание, что в азотсодержащих хлоридных растворах возможно присутствие соединения азота с хлором NH4Cl, для которого с уменьшением давления возможна реакция: NH4Cl = NH3 + HCl. При фазовой сепарации азотсодержащего флюида продукты этой реакции могут фракционировать в малоплотную газовую фазу (Петренко и др., 1989).

Если рассмотреть тройные системы H2O-хлориды с разными неполярными газами (CO2, CH4), то для них отмечалось одинаковое поведение в отношении флюидной несмесимости при высоких Р и Т (до 5 кбар), с расширением границы сольвуса при высокой концентрации соли, особенно в системах с CaCl2. При отсутствии экспериментальных данных допускалось, что принципиальные фазовые соотношения этих систем можно транслировать в систему и с неполярным газом N2 (Heinrich, 2007; Liebscher, 2010 и ссылки там). Принимая во внимание, что более поздними исследованиями показано существование поля несмесимости флюидов в системах H2O–NaCl–CO2 и H2O–CaCl2–CO2 при значительно больших Р и Т в диапазоне НР гранулитов (Аранович и др., 2010б; Иванов, Бушмин, 2019), следует ожидать возможность распада гомогенного флюида при близких высоких Р-Т условиях в тройной системе с азотом. Такие геологические примеры известны. Например, ассоциация одновозрастных включений концентрированных рассолов, включений СО2, включений СО2–N2 и СО2–CH4–N2 обнаружена в Px гранулитах на пике НР гранулитового метаморфизма в Танзании (Herms, Schenk, 1998).

Так как в нашем случае в изученных породах присутствуют не только поздние s-включения азота, но и водно-азотные ps-включения, которые в одних и тех же зонах ассоциируют с высокоплотными включениями СО2 и рассолами, мы считаем, что наблюдаемые азотные включения являются свидетельством присутствия азота в несмесимых флюидных фазах во время их захвата во включения и не связаны с поздней инфильтрацией азотсодержащего флюида в гранулиты. В связи с этим обращает на себя внимание тот факт, что CO2 и N2 флюиды обычно присутствуют и несмесимы с CaCl2-рассолами в эклогитах (например, Xiao et al., 2000). Включения Н2О–N2, весьма вероятно, могут быть свидетельством глубинного внешнего источника флюидов, так как в алмазах не только устанавливается примесь азота, но и присутствуют Н2О–N2 флюидные включения (Aranovich et al., 2016 и ссылки там). По нашему мнению, скорее, азот был растворен в захваченной во включения богатой СО2 флюидной фазе и после захвата флюидных включений был отделен при высоких давлениях и более низких температурах от СО2 (Tsunogae et at., 2008 и ссылки там) и от Н2О (Аранович, 2013).

Происхождение и состав внешнего материнского флюида

Наряду с полевыми геологическими свидетельствами на внешний флюидный поток указывает и одинаковый состав флюида во включениях в гранулитах (просачивание во всем объеме породы) и сингранулитовых метасоматитах (фокусированный поток). Если бы в гранулитах состав флюида отличался, можно было предполагать сценарий гранулитового метаморфизма с участием только “местного” флюида, выделившегося при реакциях дегидратации. И наряду с этими свидетельствами на внешний источник флюида указывают изотопные данные. На мантийный генезис флюидного потока указывают параметры изотопных систем в изученных метасоматитах (неопубликованные данные С.А. Бушмина): (87Sr/86Sr)0 = 0.7031–0.7060 и εNd(Т) = +1…+4.3. Учитывая, что изотопные параметры деплетированного мантийного резервуара 2 млрд лет назад (возраст процессов: Бушмин и др., 2009) составляли εNd(Т) = +5…+6 и (87Sr/86Sr)0 = 0.701–0.7015 (Фор, 1989), можно предполагать, что в качестве источников метасоматических флюидов выступало вещество с изотопными характеристиками верхней мантии.

На существование внешнего фокусированного потока НТ флюидов в исследованных сингранулитовых НР метасоматитах указывает и анализ распределения изотопов кислорода (δ18Ο) в сосуществующих минералах сингранулитовых метасоматитов (Аранович и др., 2009, 2010а) тех же образцов, в которых изучены флюидные включения. Были не только подтверждены петрологические данные о больших температурах (около 900°С) при НР гранулитовом петрогенезисе, но и сделан вывод, что метасоматические процессы происходили на пике метаморфизма под воздействием внешнего флюида. При неодинаковом интегральном объемном отношении флюид/порода в различных зонах весьма интенсивный поток быстро остывал, был весьма кратковременным и дискретно распределенным в сдвиговых зонах в гранулитах.

На внешний мантийный источник флюидов косвенно могут указывать и ассоциация ранних включений рассолов, СО2 и N2, так как они известны во многих породах, принадлежавших верхней мантии. В большинстве мантийных ксенолитов присутствуют почти чистые включения СО2, в ряде случаев содержащие Н2О и азот (например, Frezzotti, Touret, 2014 и ссылки там). Возможным источником внешнего флюида мантийного происхождения могли быть карбонатитовые расплавы верхней мантии, содержащие СО2, Н2О, кремнезем, хлориды кальция и щелочей (например, Frezzotti, Touret, 2014 и ссылки там). Эти исследования показали, что с понижением давления карбонатные расплавы превращаются во флюиды, в которых преобладают СО2 и рассолы с растворенными хлоридами, кремнеземом, щелочами. Данные по флюидным включениям предполагают промежуточные свойства между карбонатно-силикатно-водным расплавом и флюидом. Приблизительно такого состава мог быть и “материнский” гомогенный флюид для исследованных гранулитов.

Сравнение термодинамического предсказания и данных по включениям

Состав включений, их количество и результаты моделирования

Для объяснения преобладающего присутствия в гранулитах богатых СО2 включений, наряду с общепризнанным механизмом обогащения флюида газами (СО2, СН4), в результате преимущественного удаления Н2О в расплав, также обосновывается возможность поступления СО2 флюида из очагов карбонатных расплавов верхней мантии (например, Frezzotti, Touret, 2014 и ссылки там; Duan, 2014 и ссылки там).

Другой возможный механизм, объясняющий присутствие многочисленных СО2 включений в гранулитах, основывается на сценарии фазовой сепарации гомогенного флюида (см. выше), при которой состав возникающих несмесимых фаз (Н2О-соль и богатая СО2), участвующих в минералообразовании, определяется положением сольвуса. Экспериментально было показано, что эти возникающие флюидные фазы обладают свойствами, определяющими их разную транспортную способность и разную предрасположенность к захвату в минералы. Были выполнены экспериментальные исследования (Watson, Brenan, 1987; Holness, 1992 и ссылки там) углов смачивания между флюидами разного состава и породообразующими минералами (Qz, Opx, Cpx, Cal и др.). Установлено, что в зависимости от величины угла смачивания флюид образует или изолированные в межзерновом пространстве “карманы” (при углах больше 60°), или взаимосвязанные каналы вдоль границ зерен (при углах меньше 60°), что определяет способность флюида к фильтрации. Например, в системе кварц–флюид в условиях гранулитовой фации при 10 кбар для Н2О–СО2 флюида угол смачивания больше 60° и возрастает до ~100° в чистом СО2 флюиде. Напротив, для рассола Н2О–NaCl угол смачивания меньше 60° и значительно снижается при увеличении содержания соли. Таким образом, из-за чрезвычайно низкой способности к смачиванию богатый СО2 флюид, образуя изолированные “карманы”, малоподвижен и преимущественно захватывается во флюидные включения, а Н2О-соль флюид, имея высокую способность смачивания, образует сетку взаимосвязанных каналов, определяющих повышенную способность к фильтрации и удалению из породы. При этом концентрированные рассолы при более высокой способности к фильтрации будут иметь очень низкую предрасположенность к захвату в виде включений при высоких Р-Т параметрах гранулитов и это свойство может объяснять редкость включений рассолов по сравнению с богатыми СО2 включениями. А если существовал геологически значимый разрыв во времени между возникновением несмесимых флюидов и началом формирования породы, то в результате вероятного их разделения в пространстве породообразующим флюидом может оказаться одна из несмесимых флюидных фаз (богатый СО2 флюид или рассол), химический состав которой и будет преимущественно отражен во флюидных включениях.

Существуют также экспериментальные исследования, установившие возможность исчезновения уже захваченных включений Н2О-соль за счет их миграции внутри кристаллов (Зонова и др., 1996). Было показано, что миграция, движущей силой которой является градиент Т и векторного Р, могла не ограничиваться только отдельными кристаллами, но и приводить к выходу включений в межзерновое пространство за реальное время метаморфических процессов. В результате таких процессов р-включения рассола могли переходить в состояние ps- и s-включений или вообще удаляться из минеральных зерен.

Чистые включения

Включения, найденные и в гранулитах и в метасоматитах – это СО2 включения без Н2О и включения рассолов без СО2. Получение на сольвусе относительно чистого CO2 флюида и столь же свободного от CO2 водного флюида, близких по составу к бинарным системам, теоретически возможно в поле двухфазного флюида тройной системы, но при определенных составах исходного материнского гомогенного флюида (рис. 15а, 15б). Такой же результат может быть получен и при снижении Р-Т параметров. Детально механизм образования флюида, например, с высоким содержанием CO2 и низким содержанием соли при снижении Р-Т параметров можно проиллюстрировать на примере тройной системы с CaCl2 (рис. 17а). Предположим, точка 1 (аН2О = 0.50) на сольвусе 11 кбар 900°С соответствует менее плотной малосоленой фазе, образовавшейся из гомогенного флюида при пиковых условиях. Флюид в точке 1 содержит 10 мас. % CaCl2. Дальнейшее снижение давления до 10 кбар при 900°С вызывает новый распад получившегося флюида. Состав новой образовавшейся фазы в точке 2 (аН2О = 0.52) определится коннодой сольвуса с новыми параметрами – 10 кбар и 900°С. Флюидная фаза, отвечающая точке 2, будет содержать только 7 мас. % CaCl2. Дальнейшее снижение Р и Т до 9 кбар и 850°С приведет к флюиду с составом, отвечающим точке 3 (аН2О = 0.54) и содержащим 5 мас. % CaCl2. Таким образом, снижение Р и Т от 11 кбар и 900°С до 9 кбар и 850°С вызывает двукратное сокращение содержания соли в малоплотной СО2 содержащей флюидной фазе. Содержание воды в этой фазе также сокращается. Изменение Р-Т параметров в более широких пределах приведет к еще большему уменьшению солености и концентрации воды в СО2 содержащей флюидной фазе. Также следует принимать во внимание, что в обычно встречающихся чистых включениях СО2 на их стенках присутствует тонкая невидимая пленка воды (Н2О до 20 мол. %), которую трудно обнаружить под микроскопом (Touret, 1985; Crawford, Hollister, 1986 и ссылки там). Недавними исследованиями методом рамановской спектроскопии Н2О была установлена в тех флюидных включениях гранулитов, которые ранее считались чистыми включениями СО2 (например, Lamadrid et al., 2014).

Рис. 17.

(а) Фрагмент фазовой диаграммы системы Н2О–СО2–CaCl2: сценарий выделения из исходного гомогенного флюида фазы с высоким содержанием CO2 при снижении Р-Т параметров. Толстые линии – линии сольвуса, разделяющие области гомогенного и двухфазного флюида, тонкие линии – конноды сосуществующих флюидных фаз. (б) Фрагмент фазовой диаграммы системы Н2О–СО2–NaCl: толстая линия – линия сольвуса, разделяющая гомогенный и двухфазный флюид, пунктир – коннода двухфазного флюида с аН2О = 0.40, точки – произвольно выбранные составы исходного гомогенного флюида с разным содержанием соли, соответствующие разным количественным соотношениям возникающих при сепарации флюидных фаз. (в) Сценарий расслоения тройной системы H2O–NaCl–CO2. Тенденция изменения псевдобинарных составов флюида на сольвусе при снижении давления и температуры: водно-солевой флюид сдвигается в более “соленую” и менее “углекислотную” область, а богатый СО2 флюид – в более “углекислотную” область. x – мольные доли, m – массовые доли.

Многие годы существует представление, что СО2 и СН4, имея очень низкую растворимость в гранитных расплавах, остаются во флюиде, а Н2О, напротив, перемещается в мигматитовый расплав и этот сценарий давно используется для объяснения богатых СО2 включений в гранулитах (например, Touret, Huizenga, 2012 и ссылки там). Поскольку исследованные чистые включения СО2 обнаруживают меньшую плотность, соответствующую давлениям, которые на ~1–1.5 кбар меньше пиковых значений по данным термобарометрии (рис. 13), они, скорее всего, образовались в результате “перераспределения” флюида, захваченного в более раннюю пиковую стадию. Предполагалось, что причиной изменения плотности, отсутствия Н2О и примеси СН4 в СО2 включениях может быть диффузионная потеря водорода, который расходуется согласно реакции: СН4 + Н2О = СО2 + Н2. В зависимости от степени диффузии Н2 плотность флюида во включении может меняться, а состав приближаться к чистому СО2 флюиду (Roedder, 1984). Различные причины потери воды из включений неоднократно обсуждались (напримеер, Cordier et al., 1994; Carlson et al., 2015). Предполагалось и подтверждено экспериментально (например, Diamond, Tarantola, 2015), что перераспределение флюида с выборочной утечкой Н2О из включения, богатого СО2, может происходить даже при небольшом изменении Р и Т. Принято считать, что включения обычно могут выдерживать Р = 1–1.5 кбар, максимально до 2 кбар (Crawford, Hollister, 1986 и ссылки там). Это избыточное давление вызывает образование микротрещин, по которым преимущественно Н2О уходит из включений в силу большой способности к смачиванию и инфильтрации по сравнению с СО2. Одновременно к увеличению плотности включений могло приводить и уменьшение их объема при стрессе на фоне сдвиговых деформаций, сопровождавших образование кварцевых бластомилонитов. Возможность преимущественной потери Н2О или СО2 из включений за счет селективной диффузии и деффектов кристалла экспериментально подтверждена на примере синтетических включений H2O–NaCl, H2O–CO2 и природных включений H2O–CO2–NaCl в кварце (например, Sterner, Bodnar, 1989; Vityk et al., 2000; Baumgartner et al., 2014; Bakker, Doppler, 2016; Bakker, 2017).

При рассмотрении сценариев потери Н2О в СО2 включениях при селективной диффузии возникает вопрос – куда исчезает солевая составляющая этих включений? Ниже будет рассмотрено важное геохимическое следствие: в поле несмесимости сосуществующих флюидных фаз и в концентрированном рассоле, и в Н2О–СО2 флюидной фазе с крайне низким содержанием соли активность соли одинакова, что не зависит от ее количества. Поэтому при взаимодействии с окружающими породами Н2О–СО2 флюида низкой солености, но с такой же величиной активности, как в концентрированном расcоле, соль может быстро расходоваться в минеральных реакциях. В связи с проблемой солевой составляющей также можно обратить внимание на присутствие кальцита во включениях СО2. Включения карбонатов в СО2 включениях разных пород не являются редкостью и известны даже в ассоциации с алмазами в СО2 включениях в мантийных ксенолитах (например, Murakami et al., 2008). Что касается гранулитов, то включения карбонатов во включениях СО2 известны, например, в гранулитах Шри Ланки (Bolder-Schrijver et al., 2000). Эти авторы предполагали образование карбоната (в данном случае магнезита) из гомогенной флюидной фазы, состоящей из СО2 и MgO. В нашем случае мог быть близкий механизм образования карбоната (кальцита). В результате гидролиза CaCl2 возникающий гидроксид кальция реагировал с СО2 с образованием кальцита, согласно реакции СО2 + Са(ОН)2 = = СаСО3 + Н2О. Растущие кристаллы карбоната не использовали вещество кристаллизующегося минерала, так как это был кварц.

В завершение вышесказанного отметим, что сложность и многовариантность в объяснении причин появления почти чистых включений СО2 и включений рассолов без СО2, возможно, связана и с недостатком знаний физико-химических свойств природных флюидов, содержащих, кроме NaCl, такую важную соль, как CaCl2, которая значительно влияет на положение тройного сольвуса по отношению к бинарным системам. В тройной системе с увеличением содержаний CaCl2 и СО2 линии сольвуса быстро приближаются к краевым (простым) бинарным системам (pис. 15а, 15б).

Состав исходного гомогенного флюида и соотношение возникающих флюидных фаз

При фазовой сепарации вариации состава исходного гомогенного флюида существенно влияют на величины массовых и объемных долей возникающих флюидных фаз. В качестве примера рассмотрим вариант, представленный на рис. 17б. Изображенная на рисунке коннода соответствует двухфазному флюиду – рассолу и Н2О–СО2 флюиду с содержанием соли 72 и 2 мас. % соответственно. Образование этих фаз является результатом распада исходного гомогенного флюида при заданных P-T параметрах. Состав этого исходного флюида, в принципе, мог отвечать любой точке составов на этой конноде. Однако соотношение количеств и, соответственно, объемов, занимаемых рождающимися при распаде флюидными фазами, резко зависит от положения этих точек составов. Приведенным для примера составам с mNaCl = 0.05, mNaCl = = 0.08 и mNaCl = 0.16 соответствует соотношение массовых долей более плотного (рассол) и менее плотного (Н2О–СО2) флюидов 1/50, 1/12 и 1/3 и, соответственно, объемных долей 1/30, 1/8 и 1/2. Эта зависимость также может учитываться при анализе преобладающего присутствия в минералах флюидных включений определенного состава, например включений СО2.

Причины появления первичных включений с галитом

Включения рассола с галитом повышенной солености в изученных метасоматитах встречаются редко. Обычно это одиночные включения и, поэтому мы предполагаем, что они первичные. Еще в работе (Touret, 1985) отмечены аналогичные включения с галитом в гранулитах юга Норвегии. Петрографическое изучение минералов, содержащих включения с галитом, позволяют считать некоторые первичные включения с галитом относительно более поздними, так как они встречаются в кварце крупнозернистых участков перекристаллизации матричного кварца кварцевых бластомилонитов. Или встречаются в кварце, заключенном в гранат или ортопироксен, которые являются более поздними порфиробластами, по сравнению с пластинчатыми формами этих же минералов кварцевых бластомилонитов. Также они встречены в более крупных порфиробластах граната ортопироксен-гранатовых жил. Такие включения с рассолами (до 35 мас. % NaCl), по нашему мнению, могли формироваться при продолжавшемся во время кристаллизации пород относительно быстром снижении P и T от пиковых величин. Мы объясняем их появление следующим образом. Мог продолжаться процесс распада ранее возникших при расслоении флюидных фаз. Так могут появляться все более и более концентрированные рассолы и все более богатые СО2 и бедные солью флюидные фазы, соответствующие концам коннод на рис. 17в. Если предположить, что имела место диффузия Na между p- и ps-включениями флюида в кварце и внешним флюидом поздних зон s-включений в кварцевых бластомилонитах, то более вероятна потеря Na, чем его приобретение, так как этот поздний флюид имеет очень низкие концентрации NaCl. Таким образом, в однажды расслоившемся гомогенном флюиде при снижении P и T будет продолжаться расслоение его флюидных фаз, близких по составу к бинарным системам. То есть состав флюида и несмешивающихся фаз после распада будет контролироваться положением сольвуса в системе, и это дает нам представление об основной тенденции изменения содержания солей и СО2 во флюидных фазах, как еще не захваченных, так и захваченных во включения. Если снижение быстрое и сопоставимо с процессом кристаллизации породы, то в составе первичных флюидных включений новых генераций минералов будет фиксироваться эволюция состава этих краевых флюидных систем.

Возникает вопрос – почему в таких включениях рассолов большая концентрация NaCl? Одно из возможных объяснений – процессы гидролиза CaCl2 и NaCl. Степень гидролиза первой соли значительно больше (Bischoff et al., 1996), и после фракционирования продуктов гидролиза солевая система возможно смещается в сторону увеличения содержания NaCl. Такая тенденция наблюдается, например, в исследованных диопсид-скаполитовых жилах. Скаполитовые жилы, ассоциирующие в одних и тех же сдвиговых зонах с не содержащими кальций кварцевыми бластомилонитами и Opx-Grt жилами, мы считаем зонами комплементарного переотложения кальция, вынесенного из зон кварцевых бластомилонитов и Opx-Grt пород. Это согласуется с преобладанием CaCl2 во флюиде при образовании метасоматитов. Ранний метасоматический парагенезис с более кальциевым скаполитом (мейонита 77–81 мол. %) частично замещается более поздним парагенезисом с менее кальциевым скаполитом (содержание мейонита 71–74 мол. %). Это согласуется как с присутствием первичных Н2О включений с галитом в относительно более поздних генерациях метасоматических граната и ортопироксена, так и с возможным повышением содержания NaCl в смеси солей при снижении Р-Т параметров, исходя из термодинамического предсказания.

Сравнение расчетов активности воды по минеральным равновесиям с термодинамическим предсказанием

Термодинамическое предсказание фазового состояния, состава сосуществующих флюидов и активности Н2О в системе Н2О–СО2-соль сопоставимо с выводами, основанными на данных по флюидным включениям, и с расчетами аН2О по минеральным парагенезисам. Пределы вариаций аН2О, рассчитанные по равновесиям (гранулиты 0.49 и метасоматиты 0.51) вполне соответствуют диапазону максимальных активностей (система с NaCl: 0.50–0.53, система с CaCl2: 0.52–0.57) на сольвусах двухфазовых полей (рис. 15а, 15б). Также хорошо согласуются расчетные величины аН2О по минеральным равновесиям в конкретных образцах, полученные методом мультиравновесной термобарометрии с термодинамическим предсказанием по модельным системам. В поле двухфазного флюида расчетные величины активностей не могут быть больше модельных величин, что показано на приводимых примерах:

Обр. Б1017-9-1, Р = 8659 бар, Т = 786°С (гранулит):

Мультиравновесная термобарометрия: aH2O = = 0.49.

Максимальная активность на сольвусе двухфазного флюида с CaCl2: аН2О = 0.58.

Обр. Б1021-13-1, Р = 9349 бар, Т = 864°С (сингранулитовый метасоматит):

Мультиравновесная термобарометрия: aH2O = = 0.51.

Максимальная активность на сольвусе двухфазного флюида с CaCl2: аН2О = 0.56.

Некоторые геохимические следствия несмесимости в гранулитовом флюиде

Постоянство активностей компонентов несмесимых флюидных фаз

Это важное геохимическое следствие распада гомогенного флюида рассмотрим на примере сольвуса 900°C и 11 кбар тройной системы с CaCl2 (рис. 15б). В поле несмесимости флюидных фаз активность соли (как и других компонентов) и в концентрированном рассоле (75 мас. % соли), и в Н2О–СО2 флюидной фазе, богатой СО2 (2 мас. % соли), одинакова, постоянна (0.31) и не зависит от количества соли. Условие постоянства активностей (химических потенциалов) компонентов флюида, независимо от изменения их количества, близко к концепции Д.С. Коржинского (Коржинский, 1949) о природных физико-химических системах с вполне подвижными компонентами, так как эти параметры определяются фазовым равновесием несмесимых флюидов, а не реакциями флюид–порода.

Тенденции изменения кислотности флюидных фаз при меняющихся Р-Т параметрах

Преобладающие компоненты глубинных водных флюидов СО2, Cl, Si, Ca, K, Na, ±N2, ±CH4 (Frezzotti, Touret, 2014 и ссылки там), поступающих из верхней мантии в нижнюю кору на уровень НР гранулитов, присутствуют как в виде нейтральных молекул, так и в виде ионов, что зависит от степени диссоциации электролитов (соли, кислоты, щелочи) и определяет кислотность-щелочность флюидной фазы. При снижении Р и Т изменение степени диссоциации различно для разных электролитов. Диссоциация HCl и H2CO3 (главных кислот природных флюидов) возрастает с понижением температуры, также возрастает и диссоциация щелочей, солей. Следовательно, важно, какой электролит преобладает в гомогенном флюиде или несмесимых флюидных фазах. И важна относительная скорость изменения диссоциации, которая выше у кислот по сравнению с солями и щелочами, что может приводить к повышению кислотности (снижению рН) флюидов (например, Говоров, 1977; Наумов В., Наумов Г., 1980; Рафальский, 1987; Bushmin, 2000; Кигай, Тагиров, 2010; Кольцов, 2015; Manning, 2018 и ссылки там; Mei et al., 2018). В поле двухфазного флюида может проявляться и эффект гидролиза солей CaCl2 и NaCl – фракционирование продуктов гидролиза гидроксида и кислоты, соответственно, между более плотной и менее плотной флюидными фазами. Степень гидролиза усиливается с повышением температуры, но увеличение давления уменьшает его степень, поэтому заметное проявление гидролиза при высоких давлениях должно смещаться в более высокотемпературную область. Степень гидролиза у CaCl2, относительно NaCl, значительно больше (Bischoff et al., 1996) и в системах с CaCl2 гидролиз, по-видимому, более вероятен, также из-за образования, наряду с HCl, малорастворимого гидроксида Ca(OH)2. Это позволяет предполагать геохимический эффект гидролиза с возможностью обогащения HCl менее плотной флюидной фазы и при значительных давлениях.

При высоких температурах и давлениях HP гранулитов (~700–900°С, 10 кбар) диссоциация NaOH значительно больше, чем у NaCl и HCl, которая в этих условиях является слабой кислотой (Manning, 2018; Mei et al., 2018). Поэтому возникающая при распаде гомогенного флюида фаза рассола NaCl, содержащая NaOH, будет обладать щелочными свойствами, которые будут усиливаться из-за повышения диссоциации с понижением температуры (ниже 700–800°С). А возникающая богатая летучими флюидная фаза (Н2О–СО2 ± рассол), содержащая H2CO3 и HCl, – будет только потенциально кислотной со свойствами, близкими к нейтральным, и кислотные свойства проявятся только при снижении Р-Т параметров.

При этом кислоты HCl и H2CO3 будут вести себя по-разному при изменении давления и температуры. Хотя диссоциация HCl и H2CO3 возрастает с понижением температуры, при 800–600°С (температура гранулитовой-амфиболитовой фаций) степень диссоциации H2CO3 ничтожна, кислотность и ее увеличение при снижении температуры обеспечивается диссоциацией, главным образом, HCl (Mei et al., 2018). При температурах меньше 450–400°С (температура зеленосланцевой фации), когда HCl практически полностью диссоциирована, снижение температуры влияет практически только на диссоциацию H2CO3, в связи с чем будет усиливаться роль углекислоты как кислотного компонента флюида.

В завершении подчеркнем, что сказанное выше будет реализовываться в природных флюидах, придавая им кислотные или щелочные свойства, если несмесимые флюиды, раздельно обогащенные щелочами или кислотами, независимо эволюционируют при снижении Р-Т параметров в процессе удаления от мест своего зарождения. Кислотные или щелочные свойства этих флюидов, прежде всего, определяют растворимость, а значит транспорт и отложение рудных и породообразующих элементов. Примером может служить предполагаемая связь орогенных месторождений золота в низкотемпературных частях региональных сдвиговых зон верхней коры с гранулитовыми флюидами, поступившими из нижней коры (например, Fu, Touret, 2014).

О присутствии мигматитов и плавлении

Высокая активность Н2О (~0.5) как будто бы противоречит наблюдаемому слабому плавлению гранулитов около зон флюидной метасоматической переработки. Однако противоречие может быть снято, если принять во внимание, что в гранулитовой фации благодаря присутствию концентрированных рассолов температура плавления пород повышается на 80–100°С. В этом “температурном окне” при инфильтрации водных флюидов процессы метасоматической дегидратации и гидролиза могут протекать без плавления пород (Newton et al., 1998; Аранович, 2017). Наличие эффекта повышения температуры плавления в присутствии рассола было продемонстрированно экспериментально на примере альбита (Shmulovich, Graham, 1996) и модельного гранита (Aranovich et al., 2013).

Обращает на себя внимание, что свидетельства частичного плавления (мигматиты) присутствуют только в гранулитах, как непосредственно вмещающих зоны сингранулитовых метасоматитов, так и расположенных на удалении от метасоматических зон. Тогда как в жилоподобных телах кварцевых бластомилонитов и Opx-Grt пород, не содержащих полевых шпатов, мигматитов и других признаков плавления пород нет. Исходя из того, что в гранулитовом петрогенезисе участвовали две несмесимых флюидных фазы, мы объясняем это следующим образом. Локализованный в проницаемых зонах флюид повышенной кислотности (менее плотная Н2О–СО2 флюидная фаза) при взаимодействии с вмещающими гранулитами удалял из них Na, Ca, K, Fe и Mg и в результате формировались кварцевые бластомилониты. По мере насыщения флюида этими элементами происходило их комплементарное переотложение в сопряженных зонах “базификации” (Бушмин и др., 2007) в виде Opx-Grt и Di-Scp жил (например, pис. 3б). При рассматриваемых высоких Р-Т параметрах гранулитовой фации такой концентрированный флюид мог приобретать свойства, стирающие грань между флюидом и расплавом (например, Higashino et al., 2013). Признаками небольшого количества расплавной фазы в кварцевых бластомилонитах могут быть присутствующие в изученных образцах редкие каемки альбита или калиевого полевого шпата вокруг граната среди матричного кварца.

Второй (более плотной) несмесимой флюидной фазой, в которой могла быть стерта грань между флюидом и расплавом, были щелочные Na–Ca концентрированные солевые Cl-флюиды (верхние концы коннод с aH2O = 0.40 и в меньшей степени с aH2O = 0.51 на рис. 15а, 15б). Также локализованные в проницаемых зонах эти флюиды-рассолы вызывали метасоматическую Na–Ca фельдшпатизацию основных гранулитов с образованием Qz-Pl метасоматических пород. Образование этой более легкоплавкой минеральной ассоциации инициировало парциальное плавление и мигматизацию, интенсивную эндербитизацию Px гранулитов в протяженных сдвиговых зонах района нашего исследования (pис. 1, 16). Эффекты влияния щелочного Na–Ca метасоматоза под воздействием Cl флюидов на процессы плавления при гранулитовом метаморфизме убедительно продемонстрированы, например, в (Aranovich, Safonov, 2018). Высокая вероятность появления таких несмесимых флюидов в тектонических зонах высокой проницаемости (сдвиговые зоны) позволяет объяснить пространственную ассоциацию дискретных зон с интенсивным проявлением гранитизации и сингранулитового метасоматоза в гранулитовых комплексах.

В заключение заметим, что полученные характеристики флюидов НР гранулитов Лапланского пояса Фенноскандинавского щита не могут отражать условия образования всех известных НР гранулитов. Однако полученные результаты убеждают нас в том, что в глубинном гранулитовом петрогенезисе могли играть существенную роль дискретные потоки глубинных высокотемпературных несмесимых флюидов контрастного химического состава. После пространственной сепарации, в том числе из-за разной способности к фильтрации, поступая по региональным проницаемым зонам наверх, эти флюидные фазы глубинного происхождения могут играть важную роль в магматическом, метаморфическом, метасоматическом и рудном петрогенезисе средней и верхней коры. Прежде всего, потому, что они являются важным агентом транспорта породообразующих элементов, рудных элементов и летучих. При снижении Р и Т различия в физико-химических свойствах этих флюидных фаз, включая кислотно-щелочные свойства, должны усиливаться (например, Manning, 2018).

ВЫВОДЫ

1. В диапазоне Р-Т параметров НР гранулитов существует обширная область составов Н2О флюидов с разным содержанием СО2, хлоридов натрия и кальция, в которой гомогенный флюид распадается на контрастные по химическому составу и свойствам несмесимые флюидные фазы и эта область значительно расширяется с увеличением содержания CaCl2.

2. Следовательно, нижняя кора уровня НР гранулитовой фации может являться областью зарождения высокотемпературных несмесимых флюидов. С одной стороны, это более плотная флюидная фаза относительно щелочных рассолов с высоким содержанием воды и хлоридов, гидроксидами щелочных и щелочноземельных металлов, с другой – менее плотная потенциально кислотная фаза Н2О–СО2 флюидов, содержащая кислоты HCl и H2CO3. Поступая по региональным проницаемым зонам наверх, эти флюиды глубинного происхождения могут играть важную роль в магматическом, метаморфическом, метасоматическом и рудном петрогенезисе средней и верхней коры.

Благодарности. Авторы благодарны Л.Я. Арановичу за внимательное прочтение рукописи статьи и конструктивные замечания, А.Б. Кольцову, с которым обсуждались закономерности изменения кислотно-щелочных свойств флюидов, а также Д.П. Крылову, предоставившему один из образцов гранулита для наших исследований.

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

Список литературы

  1. Аранович Л.Я. Флюиды гранулитовой фации: физико-химический аспект // Гранулитовые комплексы в геологическом развитии докембрия и фанерозоя. СПб.: ИГГД РАН, 2007. С. 35–39.

  2. Аранович Л.Я. Флюидно-минеральные равновесия и термодинамические свойства смешения флюидных систем // Петрология. 2013. Т. 21. № 6. С. 588–599.

  3. Аранович Л.Я. Роль рассолов в высокотемпературном метаморфизме и гранитизации // Петрология. 2017. Т. 25. № 5. С. 491–503.

  4. Аранович Л.Я., Бортников Н.С., Бушмин С.А. и др. Флюидные потоки в региональных зонах деформаций // Петрология. 2009. Т. 17. № 4. С. 415–436.

  5. Аранович Л.Я., Дубинина Е.О., Авдеенко А.C. и др. Изотопный состав кислорода сосуществующих минералов силлиманит-гиперстеновых пород Порьей губы: свидетельства участия флюида в метаморфизме гранулитовой фации // Геохимия. 2010а. № 8. С. 787–800.

  6. Аранович Л.Я., Закиров И.В., Сретенская Н.Г., Геря Т.В. Тройная система H2O–CO2–NaCl при высоких Р-Т параметрах: эмпирическая модель смещения // Геохимия. 2010б. № 5. С. 1–10.

  7. Беляев О.А. Кислотное выщелачивание и сопряженный Fe–Mg метасоматоз в условиях гранулитовой фации // Метасоматоз и метасоматиты в метаморфических комплексах докембрия. Апатиты: Изд. Кол. ФАН СССР, 1981. С. 10–19.

  8. Борисенко А.С. Изучение солевого состава газово-жидких включений в минералах методом криометрии // Геология и геофизика. 1977. № 8. С. 16–27.

  9. Бушмин С.А. Метасоматические образования в зонах регионального метаморфизма // Геологическая съемка метаморфических и метасоматических комплексов (ред. В.А. Глебовицкий). 1996. СПб.: Изд-во ВСЕГЕИ, 1996. С. 84–125.

  10. Бушмин С.А., Доливо–Добровольский Д.В., Лебедева Ю.М. Инфильтрационный метасоматоз в условиях гранулитовой фации высоких давлений (на примере ортопироксен-силлиманитовых пород сдвиговых зон Лапландского гранулитового пояса) // Докл. АН. 2007. Т. 412. № 3. С. 383–387.

  11. Бушмин С.А., Глебовицкий В.А. Схема минеральных фаций метаморфических пород // Записки Российского минералогического общества. 2008. Ч. 137. № 2. С. 1–13.

  12. Бушмин С.А., Глебовицкий В.А., Савва Е.В. и др. Возраст высокобарического метасоматоза в зонах сдвиговых деформаций при коллизионном метаморфизме в Лапландском гранулитовом поясе: U-Pb-SHRIMP-II-датирование цирконов из силлиманит-гиперстеновых пород Порьегубского покрова // Докл. АН. 2009. Т. 428. № 6. С. 792–796.

  13. Говоров И.Н. Термодинамика ионно-минеральных равновесий и минералогения гидротермальных месторождений. М.: Наука, 1977. 239 с.

  14. Доливо-Добровольский Д.В. Происхождение и условия образования сапфиринсодержащих пород Центрально-Кольской гранулито-гнейсовой области // Автореф. дисс. … канд. геол.-мин. наук. СПб.: ИГГД РАН, 2002. 24 с.

  15. Иванов М.В., Бушмин С.А. Уравнение состояния флюидной системы H2O–CO2–CaCl2 и свойства флюидных фаз при Р-T параметрах средней и нижней коры // Петрология. 2019. Т. 27. № 4. С. 431–445.

  16. Зонова И.А., Синогейкин С.В., Шмулович К.И. Подвижность флюидных включений в кварце при температурном градиенте // Докл. АН. 1996. Т. 346. № 3. С. 380–382.

  17. Кигай И.Н., Тагиров Б.Р. Эволюция кислотности рудообразующих флюидов, обусловленная гидролизом хлоридов // Петрология. 2010. Т. 18. № 3. С. 270–281.

  18. Козлов Н.Е. Лапландский гранулитовый пояс – первичная природа и развитие. Апатиты: Изд-во КНЦ АН СССР, 1990. 170 с.

  19. Козлова Н.Е., Балаганский В.В., Богова М.Н., Реженова С.А. Структурно-петрологическое изучение ортопироксен-силлиманитовой ассоциации лапландских гранулитов // Изв. АН СССР. Сер. геол. 1991. № 4. С. 66–76.

  20. Кольцов А.Б. Влияние источников и путей эволюции растворов на состав метасоматитов // Геохимия. 2015. № 2. С. 144–161.

  21. Кориковский С.П., Аранович Л.Я. Чарнокитизация и эндербитизация основных гранулитов Лапландского гранулитового пояса (южная часть Кольского полуострова, район Порьей губы): I. Петрология и геотермобарометрия // Петрология. 2010. Т. 18. № 4. С. 340–368.

  22. Кориковский С.П., Аранович Л.Я. Чарнокитизация ортопироксен-клинопироксен-флогопитовых бесполевошпатовых метаультрамафитов в Лапландском гранулитовом поясе (юг Кольского полуострова): изменение состава пород и минералов, Р-Т параметры, флюидный режим // Петрология. 2015. Т. 23. № 3. С. 211–250.

  23. Коржинсккий Д.С. Биметасоматические флогопитовые и лазуритовые месторождения архея Прибайкалья. М.: Изд-во АН СССР, Тр. ИГН АН СССР. 1947. Вып. 29. № 10. 164 с.

  24. Коржинский Д.С. Открытые системы с вполне подвижными компонентами и правило фаз // Изв. АНСССР. Сер. геол. 1949. № 2. С. 3–14.

  25. Коржинский Д.С. Теория метасоматической зональности. М.: Наука, 1982. 104 с.

  26. Лебедева Ю.М., Бушмин С.А., Глебовицкий В.А. Термодинамические условия метасоматоза в высокотемпературных и высокобарических зонах сдвиговых деформаций (Кандалакшско-Умбинская зона, Кольский полуостров) // Докл. АН. 2012. Т. 445. № 2. С. 191–195.

  27. Наумов В.Б., Наумов Г.Б. Минералообразующие флюиды и физико-химические закономерности их эволюции // Геохимия. 1980. № 10. С. 1450–1460.

  28. Перчук Л.Л., Геря Т.В., Корсман К. Модель чарнокитизации гнейсовых комплексов // Петрология. 1994. Т. 2. № 5. С. 451–479.

  29. Петренко Г.В., Малинин С.Д., Арутюнян Л.А. Экспериментальное исследоание растворимости сульфидов тяжелых металлов в гетерогенных системах хлоридный раствор – паровая фаза // Геохимия. 1989. № 6. С. 882–887.

  30. Подлесский К.К. Гиперстен в ассоциации с силлимантом и кварцем как индикатор условий метаморфизма // Докл. АН. 2003. Т. 389. № 1. С. 84–87.

  31. Прияткина Л.А. Кислотное выщелачивание в условиях гранулитовой фации метаморфизма // Докл. АН СССР. 1977. Т. 234. № 5. С . 1179–1182.

  32. Рафальский Р.П. К проблеме кислотности гидротермальных растворов // Геохимия. 1987. № 3. С. 402–415.

  33. Судовиков Н.Г. Железо-магнезиально-кальциевый метасоматоз в архее Алданского щита и некоторые вопросы “основного фронта” // Изв. АН СССР. Сер. геол. 1956. № 1. С. 29–49.

  34. Фонарев В.И., Крейлен Р. Доказательство полистадийного метаморфизма на основе изучения флюидных включений в породах Лапландского гранулитового пояса // Петрология. 1995. Т. 3. № 4. С. 379–396.

  35. Фор Г. Основы изотопной геологии. Пер. с англ. М.: Мир, 1989. 590 с. (Gunter Faure. Principles of isotope geology. Second edeition. Published by John Wiley & Sones, Inc. N.Y. 1986. 589 p.).

  36. Ходоревская Л.И., Кориковский С.П. Метасоматические гранат-клинопироксен-ортопироксен-роговообманковые жилы в метаанортозитах Колвицкого массива, Кольский полуостров: минеральный состав и связь с сингранулитовой гранитизацией // Докл. АН. 2007. Т. 415. С. 539–543.

  37. Aranovich L.Y., Newton R.C. H2O activity in concentrated NaCl solutions at high pressures and temperatures measured by the brucite-periclase equilibrium // Contrib. Mineral. Petrol. 1996. V. 125. P. 200–212.

  38. Aranovich L.Y., Newton R.C. Reversed determination of the reaction: Phlogopite + quartz = enstatite + potassium feldspar + H2O in the ranges 750–875°C and 2–12 kbar at low H2O activity with concentrated KCl solutions // Amer. Mineral. 1998. V. 83. P. 193–204.

  39. Aranovich L.Y, Podlesskii K.K. Geothermobarometry of high-Grade metapelites: simultaneously operating reactions // Geol. Soc. London. Special Publ. 1989. V. 43. № 1. P. 45–61.

  40. Aranovich L., Safonov O. Halogens in High-Grade Metamorphism // The Role of Halogens in Terrestrial and Extraterrestrial Geochemical Processes. Eds. D.E. Harlov, L. Aranovich. Springer Geochemistry. Springer International Publ. 2018. P. 713–757.

  41. Aranovich L.Y., Newton R.C. Manning C.E. Brine-assisted anatexis: Experimental melting in the system haplogranite-H2O-NaCl-KCl at deep-crustal conditions // Earth Planet. Sci. Lett. 2013. V. 374. P. 111–120.

  42. Aranovich L.Y., Makhluf A.R., Manning C.E. et al. Fluids, Melting, Granulites and Granites: A Controversy – Reply to the Commentary of  J.D. Clemens, I.S. Buick and G. Stevens // Precambrian Res. 2016. V. 278. P. 400–404.

  43. Bakker R.J. Re-equilibration processes in fluid inclusion assemblages // Minerals. 2017. V. 7. № 117. P. 1–19.

  44. Bakker R.J., Doppler G. Salinity and density modifications of synthetic H2O and H2O–NaCl fluid inclusions in re-equilibration experiments at constant temperature and confining pressure // Chemical Geol. 2016. V. 424. P. 73–85.

  45. Barbey P., Raith M. The Granulite belt of Lapland // Eds. D. Vielzeuf, Ph. Vidal. Granulites and Crustal Evolution. Kluwer Acad. Publ. Netherlands. 1990. P. 111–132.

  46. Baumgartner M., Ronald J., Bakker R.J., Doppler G. Re-equilibration of natural H2O–CO2-salt-rich fluid inclusions in quartz. Part 1. Experiments in pure water at constant pressures and differential pressures at 600°C // Contrib. Mineral. Petrol. 2014. V. 167. P. 1017.

  47. Berman R.G. Thermobarometry using multiequilibrium calculations: a new technique with petrologic applications // Canadian Mineral. 1991. V. 29. № 4. P. 833–855.

  48. Berman R.G., Aranovich L.Ya. Optimized standard state and solution properties of minerals: I. Model calibration for oli-vine, orthopyroxene, cordierite, garnet, and ilmenite in the system FeO–MgO–CaO–Al2O3–TiO2–SiO2 // Contrib. Mineral. Petrol. 1996. V. 126. № 1–2. P. 1–24.

  49. Berman R.G., Aranovich L.Y., Rancourt P., Mercier P.H.J. Reversed phase equilibrium constraints on the stability of Mg–Fe–Al biotite // Amer. Mineral. 2007. V. 92. P. 139–150.

  50. Bischoff J.L., Rosenbauer R.J., Fournier R.O. The generation of HCl in the system CaCl2–H2O: Vapor-liquid relations from 380–500°C // Geochim. Cosmochim. Acta. 1996. V. 60. P. 7–16.

  51. Bodnar R.J., Vityk O. Interpretation of microthermometric data for H2O–NaCl inclusions // Fluid inclusions in minerals: Methods and applications. Eds. B. De Vivo, M.L. Frezzotti. Virginia: Tech, 1994. P. 117–130.

  52. Bolder-Schrijver L.J.A., Kriegsman L.M., Touret J.L.R. Primary carbonate-CO2 inclusions in sapphirine-bearing granulites from central Sri Lanka // J. Metam. Geol. 2000. V. 18. P. 259–269.

  53. Bowers T.S., Helgeson H.C. Calculation of the thermodynamic and geochemical consequences of nonideal mixing in the system H2O–CO2–NaC1 on phase relations in geologic systems. Equation of state for H2O–CO2–NaCl fluids at high pressures and temperatures // Geochim. Cosmochim. Acta. 1983. V. 47. P. 1247–1275.

  54. Bushmin S.A. Evolutional model of metasomatism in metamorphic cycle // Theophrastus contributions to advanced studies in Geology. V. 3. Models and Modelling of geologic processes and objects. Eds. V.A. Glebovitsky, V.N. Dech. Theophrastus Publication. St.-Petersburg: Athens, 2000. P. 137–140.

  55. Bushmin S.A., Glebovitsky V.A. Scheme of mineral facies of metamorphic rocks and its application to Fennoscandian shield with representative sites of orogenic gold mineralization // Transactions of Karelian Research Centre RAS. Precambrian Geology Series. 2016. № 2. P. 3–27.

  56. Carlson W.D., Hixon D., Garber M., Bodnar J. Controls on metamorphic equilibration: the importance of intergranular solubilities mediated by fluid composition // J. Metamorph. Geol. 2015. V 33. P. 123–146.

  57. Chu H., Chi G., Chou I-M. Freezing and melting behaviors of H2O–NaCl–CaCl2 solutions in fused silica capillaries and glass-sandwiched films: implications for fluid inclusion studies // Geofluids. 2016. V. 16. P. 518–532.

  58. Cordier P, Doukhan J.C., Ramboz C. Influence of dislocations on water leakage from fluid inclusions in quartz: A quantitative reappraisa // Eur. J. Mineral. 1994. V. 6. P. 745–752.

  59. Crawford M.L., Hollister L.S. Metamorphic fluids, the evidence from fluid inclusions // Fluid-rock interactions during metamorphism. Eds. J.V. Walther, B.J. Wood. N.Y.: Springer, 1986. P. 1–35.

  60. Damman A.H., Kars S.M., Touret J.L.R. et al. PIXE and SEM analyses of fluid inclusions in quartz crystals from the K-alteration zone of the Rosia Poeini porphyry-Cu deposit, Apuseni Mountains, Rumania // Eur. J. Mineral. 1996. V. 8. P. 1081–1096.

  61. Diamond L.W., Tarantola A. Interpretation of fluid inclusions in quartz deformed by weak ductile shearing: Reconstruction of differential stress magnitudes and pre-deformation fluid properties // Earth Planet. Sci. Lett. 2015. V. 417. P. 107–119.

  62. Duan Z., Moller N., Weare J.H. Equation of state for the NaCl–H2O–CO2 system. Prediction of phase equilibria and volumetric properties // Geochim. Cosmochim. Acta. 1995. V. 59. P. 2869–2882.

  63. Duan X. A general model for predicting the solubility beha-vior of H2O–CO2 fluids in silicate melts over a wide range of pressure, temperature and compositions // Geochim. Cosmochim. Acta. 2014. V. 125. P. 582–609.

  64. Fonarev V.I., Touret J.L.R., Kotelnikova Z.A. Fluid inclusions in rocks from the Central Kola Granulite area (Baltic Shield) // Eur. J. Mineral. 1998. V. 10. P. 1181–1200.

  65. Frezzotti M.-L., Touret J.L.R. CO2, carbonate-rich melts, and brines in the mantle // Geoscience Frontiers. 2014. V. 5. P. 697–710.

  66. Fu B., Touret J.L.R. From granulite fluids to quartz-carbonate megashear zones: The gold rush // Geoscience Frontiers. 2014. V. 5. P. 747–758.

  67. Hammerli J., Rubenach M. The Role of Halogens During Regional and Contact Metamorphism // The Role of Halo-gens in Terrestrial and Extraterrestrial Geochemical Processes: Surface, Crust, and Mantle. 2018. Ch. 10. P. 10.

  68. Herms P., Schenk V. Fluid inclusions in high pressure granulites of the Pan-African belt in Tanzania (Uluguru Mus): A record of prograde to retrograde fluid evolution // Contrib. Mineral. Petrol. 1998. V. 130. P. 199–212.

  69. Heinrich W.H. Fluid Immiscibility in Metamorphic Rocks // Rev. Mineral. Geochem. 2007. V. 65. P. 389–430.

  70. Hisada K., Perchuk L.L., Gerya T.V. et al. P-T-fluid evolution in the Mahalapye Complex. Limpopo high-grade terrane, eastern Botswana // J. Metamorph. Geol. 2005. V. 23. № 5. P. 313–334.

  71. Higashino F., Kawakami T., Satish-Kumar M. et al. Chlorine-rich fluid or melt activity during granulite facies metamorphism in the Late Proterozoic to Cambrian continental collision zone-An example from the Sor Rondane Mountains, East Antarctica // Precambrian Res. 2013. V. 234. P. 229–246.

  72. Holness M.B. Equilibrium dihedral angles in the system quartz–CO2–H2O–NaCl at 800°C and 1–15 kbar: the effect of pressure and fluid composition on the permeability of quartzites // Earth Planet. Sci. Lett. 1992. V. 114. P. 171–184.

  73. Ivanov M.V., Bushmin S.A. Thermodynamic model of the system H2O–CO2–CaCl2 at high P-T parameters // Experi-ment in Geosciences. 2018. V. 24. №. 1. P. 110–113.

  74. Klatt E., Schoch A.E. Comparison of fluid inclusion characteristics for high-grade metamorphic rocks from Finnish Lapland with medium-grade metamorphic rocks from northwestern Cape Province, South Africa // 27th Intern. Geol. Congr. M., 1974. V. 10. P. 73.

  75. Koizumi T., Tsunogae T., van Reenen D.D. Fluid evolution of partially retrogressed pelitic granulite from the Southern Marginal Zone of the Neoarchean Limpopo Complex, South Africa: Evidence from phase equilibrium modeling // Precambrian Res. 2014. V. 253. P. 146–156.

  76. Kooi M.E., Schouten J.A., Van den Kerkhof A.V. et al. The system CO2–N2 at high pressure and applications to fluid inclusions // Geochim. Cosmochim. Acta. 1998. V. 62. P. 2837–2843.

  77. Lamadrid H.M., Lamb W.M., Santosh M., Bodnar R.J. Raman spectroscopic characterization of H2O in CO2-rich fluid inclusions in granulite facies metamorphic rocks // Gondwana Res. V. 26. 2014. P. 301–310.

  78. Liebscher A. Aqueous fluids at elevated pressure and temperature // Geofluids. 2010. V. 10. P. 3–19.

  79. Manning C.E. Fluids of the lower Crust: Deep is different // Annu. Rev. Earth Planet. Sci. 2018. V. 46. P. 67–97.

  80. Manning C.E., Aranovich L.Y. Brines at high pressure and temperature: Thermodynamic, petrologic and geochemical effects // Precambrian Res. 2014. V. 253. P. 6–16.

  81. Mei Y., Liu W., Brugger J. et al. The dissociation mechanism and thermodynamic properties of HCl(aq) in hydrothermal fluids (to 700°C, 60 kbar) by ab initio molecular dynamics simulations // Geochim. Cosmochim. Acta. 2018. V. 226. P. 84–106.

  82. Murakami T., Wallis S., Enami M., Kagi H. Forearc diamond from Japan // Geology. 2008. V. 36. P. 219–222.

  83. Newton R.C. Fluids of granulite facies metamorphism // Advances in Physical Geochemistry, Fluid-rock interactions during metamorphism. Eds. J.V. Walther, B.J. Wood. N.Y.: Springer-Verlag, 1986. V. 5. P. 36–59.

  84. Newton R.C., Aranovich L.Ya., Hansen E.C., Vandenheuvel B.A. Hypersaline fluids in Precambrian deep-crustal metamorphism // Precambrian Res. 1998. V. 91. P. 41–63.

  85. Newton R.C., Smith J., Windley V. Carbonic metamorphism, granulites and crustal growth // Nature. 1980. V. 288. P. 45–50.

  86. Newton R.C, Manning C.E. Role of saline fluids in deep-crustal and upper-mantle metasomatism: insights from experimental studies // Geofluids. 2010. V. 10. P. 58–72.

  87. Newton R.C., Touret J.L., Aranovich L.Y. Fluids and H2O activity at the onset of granulite facies metamorphism // Precambrian Res. 2014. V. 253. P. 17–25.

  88. Nijland T.G., Touret J.L.R., Visser D. Anomalously low temperature orthopyroxene, spinel, and sapphirine occurrences in metasediments from the Bamble amphibolite-to-granulite facies transition zone (South Norway): Possible evidence for localized action of saline fluids // J. Geol. 1998. V. 106. P. 575–590.

  89. Oakes C.S., Bodnar R.I., Simonson J.M. The system NaCl–CaCl2–H2O: I. The ice liquidus at 1 atm total pressure // Geochim. Cosmochim. Acta. 1990. V. 54. P. 603–610.

  90. Perchuk L.L., Gerya T.V. Fluid control of charnockitization // Chem. Geol. 1993. V. 108. P. 175–186.

  91. Perchuk L.L., Safonov O.G., Gerya T.V. et al. Mobility of components in metasomatic transformation and partial melting of gneisses: an example from Sri Lanka // Contrib. Mineral. Petrol. 2000. V. 140. P. 212–232.

  92. Roedder E. Fluid inclusions // Reviews in Mineralogy. Mineral. Soc. America. 1984. V. 12. 644 p.

  93. Santosh M., Tsunogae T., Yoshikura S. “Ultrahigh density” carbonic fluids in ultrahigh-temperature crustal metamorphism // J. Mineral. Petrol. Sci. 2004. V. 99. P. 164–179.

  94. Safonov O.G., Reutsky V.N., Varlamov D.A. et al. Composition and source of fluids in high-temperature graphite-bearing granitoids associated with granulites: Examples from the Southern Marginal Zone, Limpopo Complex, South Africa // Gondwana Res. 2018. V. 60. P. 129–152.

  95. Shmulovich K.I., Graham C.M. Melting of albite and dehydration of brucite in H2O–NaCl fluids to 9 kbars and 700–900°C: implications for partial melting and water activities during high pressure metamorphism // Contrib. Mineral. Petrol. 1996. V. 124. P. 370–382.

  96. Sisson V.B., Crawford M.L., Thompson P.H. CO2-brine immiscibility at high temperatures, evidence from calcareous metasedimentary rocks // Contrib. Mineral. Petrol. 1981. V. 78. P. 371–378.

  97. Sterner S.M., Bodnar R.J. Synthetic fluid inclusions in natural quartz. I. Compositional types synthesized and applications to experimental geochemistry // Geochim. Cosmochim. Acta. 1984. V. 48. P. 2659–2668.

  98. Sterner S.M., Bodnar R.J. Synthetic fluid inclusions. VII. Re-equilibration of fluid inclusions in quartz during laboratory-simulated metamorphic burial and uplift // J. Metamorph. Geol. 1989. № 7. P. 243–260.

  99. Steele-MacInnis M., Bodnar R.J., Naden J. Numerical model to determine the composition of H2O–NaCl–CaCl2 fluid inclusions based on microthermometric and microanalytical data // Geochim. Cosmochim. Acta. 2011. V. 75. P. 21–40.

  100. Takahashi K., Tsunogae T., Ugwuonah E.N. Fluid-induced high-temperature metasomatism at Rundvågshetta in the Lützow-Holm Complex, East Antarctica: Implications for the role of brine during granulite formation // Geoscience Frontiers. 2018. V. 9. P. 1309–1323.

  101. Touret J.L.R. Fluid regime in southern Norway, the record of fluid inclusionsm // The deep Proterozoic crust in the North Atlantic Provinces. Eds. A.C. Tobi, J.L.R. Touret. Dordrecht: Reidel, 1985. P. 517–549.

  102. Touret J.L.R. The role and nature of fluids in the continental lower crust // Geol. Soc. India. 1995. V. 34. P. 143–160.

  103. Touret J.L.R. Fluids in metamorphic rocks // Lithos. 2001. V. 55. P. 1–25.

  104. Touret J., Dietvorst P. Fluid inclusions in high-grade anatectic metamorphites // J. Geol. Soc. London. 1983. V. 140. P. 635–649.

  105. Touret J.L.R., Huizenga J.M. Fluid-assisted granulite metamorphism: A continental journey // Gondwana Res. 2012. V. 21. P. 224–235.

  106. Trommsdorff V., Skippen G., Ulmer P. Halite and sylvite as solid inclusions in high-grade metamorphic rocks // Contrib. Mineral. Petrol. 1985. V. 89. P. 24–29.

  107. Tsunogae T., Santosh M., Osanai Y. et al. Very high-density carbonic fluid inclusions in sapphirine-bearing granulites from Tonagh Island in the Archean Napier Complex, East Antarctica: implications for CO2 infiltration during ultrahigh-temperature (T > 1100°C) metamorphism // Contrib. Mineral. Petrol. 2002. V. 143. P. 279–299.

  108. Tsunogae T., Santosh M., Dubessy J. Fluid characteristics of high- to ultrahigh-temperature metamorphism in southern India: A quantitative Raman spectroscopic study // Precambrian Res. 2008. V. 162. P. 198–211.

  109. Van den Kerkhof A.M. The system CO2–CH4–N2 in fluid inclusions: theoretical modelling and geological applications. PhD Dissertation. Amsterdam: Free University Press, 1988. 206 p.

  110. Van den Kerkhof A.M., Hein U.F. Fluid inclusion petrography // Lithos. 2001. V. 55. P. 27–47.

  111. Vityk M.O., Bodnar R.J., Doukhan J.C. Synthetic fluid inclusions. XV. TEM investigation of plastic flow associated with reequilibration of fluid inclusions in natural quartz // Contrib. Mineral. Petrol. 2000. V. 139. P. 285–297.

  112. Watson E.D, Brenan J.M. Fluids in the lithosphere: 1. Experimentally-determined wetting characteristics of CO2–H2O fluids and their implications for fluid transport, host-rock physical properties, and fluid inclusion formation // Earth Planet. Sci. Lett. 1987. V. 85. P. 497–515.

  113. Whitney D.L., Evans B.W. Abbreviations for names of rock-forming minerals // Amer. Mineral. 2010. V. 95. P. 185–187.

  114. Xiao Y., Hoefs J., Van den Kerkhof A.M. et al. Fluid history of UHP metamorphism in Dabie Shan, China: a fluid inclusion and oxygen isotope study on the coesite-bearing eclogite from Bixiling // Contrib. Mineral. Petrol. 2000. V. 139. P. 1–16.

  115. Zwart E.W., Touret J.L.R. Melting behavior and composition of aqueous fluid inclusions in fluorite and calcite: applications within the system H2O–CaCl2–NaCl // Eur. J. Mineral. 1994. V. 6. P. 773–786.

Дополнительные материалы

скачать ESM_1.zip
Suppl. 1, ESM_1. Гранулиты: гистограмма включений СО2 высокой плотности (ps первично вторичные включения) и низкой плотности (s вторичные включения). Th температура гомогенизации СО2.
скачать ESM_2.zip
Suppl. 2, ESM_2. Сингранулитовые метасоматиты (MS): гистограммы первичных (p), первично вторичных (ps) и вторичных (s) включений CO2 в Qz бластомилонитах (Qz Blm), Th - температура гомогенизации CO2. (а) p и ps включения высокой плотности; (б) p и ps включения низкой плотности; (в) s включения низкой плотности.
скачать ESM_3.zip
Suppl. 3, ESM_3. Сингранулитовые метасоматиты (MS): гистограммы первичных (p), первично вторичных (ps) и вторичных (s) включений рассолов и азота в Qz бластомилонитах (Qz Blm). (а) p и ps включения Н2О соль, Tml начальная температура плавления льда; (б) p и ps включения Н2О соль, Tmf - конечная температура плавления льда, в скобках указан состав включений, с которыми включения рассолов находятся в ассоциации; (в) ps включения азота, Th - температура гомогенизации N2; (г) s включения азота.
скачать ESM_4.zip
Suppl. 4, ESM_4. Сингранулитовые метасоматиты (MS): гистограммы первичных (p), первично вторичных (ps) и вторичных (s) включений СО2 и азота в ортопироксен гранатовых жилах (Opx-Grt), Th температура гомогенизации. (а) p и ps включения СО2 высокой и низкой плотности; (б) s включения СО2 низкой плотности; (в) p и ps включения азота.
скачать ESM_5.zip
Suppl. 5, ESM_5. Сингранулитовые метасоматиты (MS): гистограмма первичных (p) и первично вторичных (ps) включений СО2 в диопсид скаполитовых жилах (Di-Scp). Th температура гомогенизации CO2.
скачать ESM_6.zip
Suppl. 6, ESM_6. Изохоры включений СО2.